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.

Realistic synthetic observations of theoretical source models are essential for our understanding of real observational data. In using synthetic data, one can verify the extent to which source parameters can be recovered and evaluate how various data corruption effects can be calibrated. These studies are important when proposing observations of new sources, in the characterization of the capabilities of new or upgraded instruments, and when verifying model-based theoretical predictions in a comparison with observational data. We present the SYnthetic Measurement creator for long Baseline Arrays (SYMBA), a novel synthetic data generation pipeline for Very Long Baseline Interferometry (VLBI) observations. SYMBA takes into account several realistic atmospheric, instrumental, and calibration effects. We used SYMBA to create synthetic observations for the Event Horizon Telescope (EHT), a mm VLBI array, which has recently captured the first image of a black hole shadow. After testing SYMBA with simple source and corruption models, we study the importance of including all corruption and calibration effects. Based on two example general relativistic magnetohydrodynamics (GRMHD) model images of M87, we performed case studies to assess the attainable image quality with the current and future EHT array for different weather conditions. The results show that the effects of atmospheric and instrumental corruptions on the measured visibilities are significant. Despite these effects, we demonstrate how the overall structure of the input models can be recovered robustly after performing calibration steps. With the planned addition of new stations to the EHT array, images could be reconstructed with higher angular resolution and dynamic range. In our case study, these improvements allowed for a distinction between a thermal and a non-thermal GRMHD model based on salient features in reconstructed images.

We report on the first detection of pulsed radio emission from a radio pulsar with the ALMA telescope. The detection was made in the Band-3 frequency range (85-101 GHz) using ALMA in the phased-array mode developed for VLBI observations. A software pipeline has been implemented to enable a regular pulsar observing mode in the future. We describe the pipeline and demonstrate the capability of ALMA to perform pulsar timing and searching. We also measure the flux density and polarization properties of the Vela pulsar (PSR J0835$-$4510) at mm-wavelengths, providing the first polarimetric study of any ordinary pulsar at frequencies above 32 GHz. Finally, we discuss the lessons learned from the Vela observations for future pulsar studies with ALMA, particularly for searches near the supermassive black hole in the Galactic Center, and the potential of using pulsars for polarization calibration of ALMA.

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

In April 2019, the Event Horizon Telescope (EHT) collaboration revealed the first image of the candidate super-massive black hole (SMBH) at the centre of the giant elliptical galaxy Messier 87 (M87). This event-horizon-scale image shows a ring of glowing plasma with a dark patch at the centre, which is interpreted as the shadow of the black hole. This breakthrough result, which represents a powerful confirmation of Einstein's theory of gravity, or general relativity, was made possible by assembling a global network of radio telescopes operating at millimetre wavelengths that for the first time included the Atacama Large Millimeter/ submillimeter Array (ALMA). The addition of ALMA as an anchor station has enabled a giant leap forward by increasing the sensitivity limits of the EHT by an order of magnitude, effectively turning it into an imaging array. The published image demonstrates that it is now possible to directly study the event horizon shadows of SMBHs via electromagnetic radiation, thereby transforming this elusive frontier from a mathematical concept into an astrophysical reality. The expansion of the array over the next few years will include new stations on different continents - and eventually satellites in space. This will provide progressively sharper and higher-fidelity images of SMBH candidates, and potentially even movies of the hot plasma orbiting around SMBHs. These improvements will shed light on the processes of black hole accretion and jet formation on event-horizon scales, thereby enabling more precise tests of general relativity in the truly strong field regime.

We report the discovery of maser emission in the two lowest rotational transitions of CS toward the high-mass protostar W51 e2e with ALMA and the JVLA. The masers from CS J=1-0 and J=2-1 are neither spatially nor spectrally coincident (they are separated by ~150 AU and ~30 km/s), but both appear to come from the base of the blueshifted outflow from this source. These CS masers join a growing list of rarely-detected maser transitions that may trace a unique phase in the formation of high-mass protostars.

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 astrophysical agenda of the 21st century requires a very sharp view of celestial objects. High angular resolution studies are essential for fundamental studies of a broad variety of astrophysical phenomena ranging from relativistic physics of black holes, their gravitational and electromagnetic imprints, violent transient processes, including those producing detectable gravitational waves, birth and evolution of planetary systems. Over the past decades, radio astronomy made huge leap in achieving ground-breaking angular resolution measured in tens of microarcseconds (one tenth of nanoradian and better). Recently a global Event Horizon Telescope (EHT) collaboration obtained first direct images of the shadow of a super-massive black hole in the nucleus of the active galaxy M87. These observations were conducted at 230 GHz. The two first generation Space Very Long Baseline Interferometry (VLBI) missions, VSOP/HALCA led by the Japan Aerospace Exploration Agency (JAXA) and RadioAstron led by the Russia Roscosmos State Corporation and Russia Academy of Sciences, achieved the highest angular resolution at frequencies from 0.3 to 22 GHz in observations conducted in the period 1997 - 2019. The next step in advancing high angular resolution radio astronomy is in combining high frequency (millimeter and sub-millimeter wavelengths) and interferometric baselines exceeding the Earth diameter. The present THEZA White Paper describes a combination which would unify technology developments in giga-/tera-hertz instrumentation and space-borne radio astronomy. The current preprint version of the THEZA White Paper is slightly re-formatted and edited comparing to the official submitted version.

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.

(Abridged) The CASA software suite, can now reduce very long baseline interferometry (VLBI) data with the recent addition of a fringe fitter. Here, we present the Radboud PIpeline for the Calibration of high Angular Resolution Data (rPICARD), which is an open-source VLBI calibration and imaging pipeline built on top of the CASA framework. The pipeline is capable of reducing data from different VLBI arrays. It can be run non-interactively after only a few non-default input parameters are set and delivers high-quality calibrated data. CPU scalability based on a message-passing interface (MPI) implementation ensures that large bandwidth data from future arrays can be processed within reasonable computing times. Phase calibration is done with a Schwab-Cotton fringe fit algorithm. For the calibration of residual atmospheric effects, optimal solution intervals are determined based on the signal-to-noise ratio (S/N) of the data for each scan. Different solution intervals can be set for different antennas in the same scan to increase the number of detections in the low S/N regime. These novel techniques allow rPICARD to calibrate data from different arrays, including high-frequency and low-sensitivity arrays. The amplitude calibration is based on standard telescope metadata, and a robust algorithm can solve for atmospheric opacity attenuation in the high-frequency regime. Standard CASA tasks are used for CLEAN imaging and self-calibration. In this work we demonstrate the capabilities of rPICARD by calibrating and imaging 7 mm VLBA data of the central radio source in the M87 galaxy. The reconstructed jet image reveals a complex collimation profile and edge-brightened structure. A potential counter-jet is detected that has 10 % of the brightness of the approaching jet. This constrains jet speeds close to the radio core to about half the speed of light for small inclination angles.

Currently, HOPS and AIPS are the primary choices for the time-consuming process of (millimeter) Very Long Baseline Interferometry (VLBI) data calibration. However, for a full end-to-end pipeline, they either lack the ability to perform easily scriptable incremental calibration or do not provide full control over the workflow with the ability to manipulate and edit calibration solutions directly. The Common Astronomy Software Application (CASA) offers all these abilities, together with a secure development future and an intuitive Python interface, which is very attractive for young radio astronomers. Inspired by the recent addition of a global fringe-fitter, the capability to convert FITS-IDI files to measurement sets, and amplitude calibration routines based on ANTAB metadata, we have developed the the CASA-based Radboud PIpeline for the Calibration of high Angular Resolution Data (rPICARD). The pipeline will be able to handle data from multiple arrays: EHT, GMVA, VLBA and the EVN in the first release. Polarization and phase-referencing calibration are supported and a spectral line mode will be added in the future. The large bandwidths of future radio observatories ask for a scalable reduction software. Within CASA, a message passing interface (MPI) implementation is used for parallelization, reducing the total time needed for processing. The most significant gain is obtained for the time-consuming fringe-fitting task where each scan be processed in parallel.

We present a detailed description of the special procedures for calibration and quality assurance of Atacama Large Millimeter/submillimeter Array (ALMA) observations in Very Long Baseline Interferometry (VLBI) mode. These procedures are required to turn the phased ALMA array into a fully calibrated VLBI station. As an illustration of these methodologies, we present full-polarization observations carried out with ALMA as a phased array at 3mm (Band 3) and 1.3mm (Band 6) as part of Cycle-4. These are the first VLBI science observations conducted with ALMA and were obtained during a 2017 VLBI campaign in concert with other telescopes worldwide as part of the Global mm-VLBI Array (GMVA, April 1-3) and the Event Horizon Telescope (EHT, April 5-11) in ALMA Bands 3 and 6, respectively.

We report the discovery of the first 14NH3 (2,2) maser, seen in the Sgr B2 Main star forming region near the center of the Milky Way, using data from the Very Large Array radio telescope. The maser is seen in both lower resolution (3" or ~0.1 pc) data from 2012 and higher resolution (0''.1 or ~1000 AU) data from 2018. In the higher resolution data ammonia (2,2) maser emission is detected toward 5 independent spots. The maser spots are not spatially or kinematically coincident with any other masers in this region, or with the peaks of the radio continuum emission from the numerous ultracompact and hypercompact \hii\, regions in this area. While the (2,2) maser spots are spatially unresolved in our highest resolution observations, they have unusually broad linewidths of several kilometers per second, which suggests that each of these spots consists of multiple masers tracing unresolved velocity structure. No other ammonia lines observed in Sgr B2 Main are seen to be masers, which continues to challenge theories of ammonia, maser emission that predict simultaneous maser emission in multiple ammonia transitions.

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.

The Galactic Center supermassive black hole Sagittarius A* (Sgr A*) is one of the most promising targets to study the dynamics of black hole accretion and outflow via direct imaging with very long baseline interferometry (VLBI). At 3.5 mm (86 GHz), the emission from Sgr A* is resolvable with the Global Millimeter VLBI Array (GMVA). We present the first observations of Sgr A* with the phased Atacama Large Millimeter/submillimeter Array (ALMA) joining the GMVA. Our observations achieve an angular resolution of ~87{\mu}as, improving upon previous experiments by a factor of two. We reconstruct a first image of the unscattered source structure of Sgr A* at 3.5 mm, mitigating effects of interstellar scattering. The unscattered source has a major axis size of 120 $\pm$ 34{\mu}as (12 $\pm$ 3.4 Schwarzschild radii), and a symmetrical morphology (axial ratio of 1.2$^{+0.3}_{-0.2}$), which is further supported by closure phases consistent with zero within 3{\sigma}. We show that multiple disk-dominated models of Sgr A* match our observational constraints, while the two jet-dominated models considered are constrained to small viewing angles. Our long-baseline detections to ALMA also provide new constraints on the scattering of Sgr A*, and we show that refractive scattering effects are likely to be weak for images of Sgr A* at 1.3 mm with the Event Horizon Telescope. Our results provide the most stringent constraints to date for the intrinsic morphology and refractive scattering of Sgr A*, demonstrating the exceptional contribution of ALMA to millimeter VLBI.

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 report the detection of NaCl, KCl, and their $^{37}$Cl and $^{41}$K isotopologues toward the disk around Orion SrcI. About 60 transitions of these molecules were identified. This is the first detection of these molecules in the interstellar medium not associated with the ejecta of evolved stars. It is also the first ever detection of the vibrationally excited states of these lines in the ISM above v = 1, with firm detections up to v = 6. The salt emission traces the region just above the continuum disk, possibly forming the base of the outflow. The emission from the vibrationally excited transitions is inconsistent with a single temperature, implying the lines are not in LTE. We examine several possible explanations of the observed high excitation lines, concluding that the vibrational states are most likely to be radiatively excited via rovibrational transitions in the 25-35 {\mu}m (NaCl) and 35-45 {\mu}m (KCl) range. We suggest that the molecules are produced by destruction of dust particles. Because these molecules are so rare, they are potentially unique tools for identifying high-mass protostellar disks and measuring the radiation environment around accreting young 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.

The compact radio source Sagittarius~A$^*$ (Sgr~A$^*$)in the Galactic Center is the primary supermassive black hole candidate. General relativistic magnetohydrodynamical (GRMHD) simulations of the accretion flow around Sgr\,A$^*$ predict the presence of sub-structure at observing wavelengths of $\sim 3$\,mm and below (frequencies of 86\,GHz and above). For very long baseline interferometry (VLBI) observations of Sgr\,A$^*$ at this frequency the blurring effect of interstellar scattering becomes subdominant, and arrays such as the High Sensitivity Array (HSA) and the global mm-VLBI Array (GMVA) are now capable of resolving potential sub-structure in the source. Such investigations improve our understanding of the emission geometry of the mm-wave emission of Sgr\,A$^*$, which is crucial for constraining theoretical models and for providing a background to interpret 1\,mm VLBI data from the Event Horizon Telescope (EHT). We performed high-sensitivity very long baseline interferometry (VLBI) observations of Sgr\,A$^*$ at 3\,mm using the Very Long Baseline Array (VLBA) and the Large Millimeter Telescope (LMT) in Mexico on two consecutive days in May 2015, with the second epoch including the Green Bank Telescope (GBT). We find an overall source geometry that matches previous findings very closely, showing a deviation in fitted model parameters less than 3\% over a time scale of weeks and suggesting a highly stable global source geometry over time. The reported sub-structure in the 3\,mm emission of Sgr\,A$^*$ is consistent with theoretical expectations of refractive noise on long baselines. However, comparing our findings with recent results from 1\,mm and 7\,mm VLBI observations, which also show evidence for east-west asymmetry, an intrinsic origin cannot be excluded. Confirmation of persistent intrinsic substructure will require further VLBI observations spread out over multiple epochs.

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.

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.

Understanding the formation of massive stars is one of the unsolved problems in modern astronomy. The main difficulty is that the intense radiation from the high-luminosity stars and the thermal pressure from the resulting ionized gas (both insignificant for low-mass stars) may be able to reverse the accretion flow and prevent the star from accreting fresh material. Such feedback effects can naturally be mitigated if accretion proceeds through discs, which is the established mechanism to form sun-like stars. However, recent 3D MHD simulations have shown that accretion on 1000 au scales is through filaments rather than a large disc. This theoretical prediction has never been confirmed via observations owing to the poor linear resolution of previous studies (>1000 au). Here we present the first observational evidence that mass assembly in young high-mass stars forming in protoclusters is predominantly asymmetric and disordered. In particular, we observed the innermost regions around three deeply embedded high-mass protostars with very high spatial resolution (~100 au). We identified multiple massive (several solar masses), warm (50-150 Kelvin) filamentary streamers pointing onto the central sources, which we interpret as multi-directional accretion channels. These structures inhibit the formation of a large, steady disc. Nevertheless, the identification of fast collimated outflows in the three observed systems indicates that (non-steady) compact discs may be present (we measure upper limits on their radii of <80 for one object and <350 astronomical units for the remaining two objects). Our finding contrasts with the simplified classic paradigm of an ordered (and stable) disc/jet system and provides an experimental confirmation of a multi-directional and unsteady accretion model for massive star formation supported by recent 3D (magneto)hydrodynamic simulations.

We report ALMA long-baseline observations of Orion Source I (SrcI) with resolution 0.03-0.06" (12-24 AU) at 1.3 and 3.2 mm. We detect both continuum and spectral line emission from SrcI's disk. We also detect a central weakly resolved source that we interpret as a hot spot in the inner disk, which may indicate the presence of a binary system. The high angular resolution and sensitivity of these observations allow us to measure the outer envelope of the rotation curve of the H$_2$O $5_{5,0}-6_{4,3}$ line, which gives a mass $M_I\approx15\pm2$ Msun. We detected several other lines that more closely trace the disk, but were unable to identify their parent species. Using centroid-of-channel methods on these other lines, we infer a similar mass. These measurements solidify SrcI as a genuine high-mass protostar system and support the theory that SrcI and the Becklin Neugebauer Object were ejected from the dynamical decay of a multiple star system $\sim$500 years ago, an event that also launched the explosive molecular outflow in Orion.

High-mass stars form in much richer environments than those associated with isolated low-mass stars, and once they reach a certain mass, produce ionised (HII) regions. The formation of these pockets of ionised gas are unique to the formation of high-mass stars (M $>8$ M$_\odot$), and present an excellent opportunity to study the final stages of accretion, which could include accretion through the HII region itself. This study of the dynamics of the gas on both sides of these ionisation boundaries in very young HII regions aims to quantify the relationship between the HII regions and their immediate environments.We present high-resolution ($\sim$ 0.5$"$) ALMA observations of nine HII regions selected from the Red MSX Source (RMS) survey with compact radio emission and bolometric luminosities greater than 10$^4$ L$_\odot$. We focus on the initial presentation of the data, including initial results from the radio recombination line H29$\alpha$, some complementary molecules, and the 256 GHz continuum emission. Of the six (out of nine) regions with H29$\alpha$ detections, two appear to have cometary morphologies with velocity gradients across them, and two appear more spherical with velocity gradients suggestive of infalling ionised gas. The remaining two were either observed at low resolution or had signals that were too weak to draw robust conclusions. We also present a description of the interactions between the ionised and molecular gas (as traced by CS (J=5-4)), often (but not always) finding theHII region had cleared its immediate vicinity of molecules. Of our sample of nine, the observations of the two clusters expected to have the youngest HII regions (from previous radio observations) are suggestive of having infalling motions in the H29$\alpha$ emission, which could be indicative of late stage accretion onto the stars despite the presence of an HII region.

We present the first images of the 28SiO v=1, J=2-1 maser emission around the closest known massive young stellar object Orion Source I observed at 86 GHz (3mm) with the VLBA. These images have high spatial (~0.3 mas) and spectral (~0.054 km/s) resolutions. We find that the 3mm masers lie in an X-shaped locus consisting of four arms, with blue-shifted emission in the south and east arms and red-shifted emission in the north and west arms. Comparisons with previous images of the 28SiO v=1,2, J=1-0 transitions at 7mm (observed in 2001-2002) show that the bulk of the J=2-1 transition emission follows the streamlines of the J=1-0 emission and exhibits an overall velocity gradient consistent with the gradient at 7mm. While there is spatial overlap between the 3mm and 7mm transitions, the 3mm emission, on average, lies at larger projected distances from Source I (~44 AU compared with ~35 AU for 7mm). The spatial overlap between the v=1, J=1-0 and J=2-1 transitions is suggestive of a range of temperatures and densities where physical conditions are favorable for both transitions of a same vibrational state. However, the observed spatial offset between the bulk of emission at 3mm and 7mm possibly indicates different ranges of temperatures and densities for optimal excitation of the masers. We discuss different maser pumping models that may explain the observed offset. We interpret the 3mm and 7mm masers as being part of a single wide-angle outflow arising from the surface of an edge-on disk rotating about a northeast-southwest axis, with a continuous velocity gradient indicative of differential rotation consistent with a Keplerian profile in a high-mass proto-binary.

High-mass stars have generally been assumed to accrete most of their mass while already contracted onto the main sequence, but this hypothesis has not been observationally tested. We present ALMA observations of a 3 x 1.5 pc area in the W51 high-mass star-forming complex. We identify dust continuum sources and measure the gas and dust temperature through both rotational diagram modeling of CH3OH and brightness-temperature-based limits. The observed region contains three high-mass YSOs that appear to be at the earliest stages of their formation, with no signs of ionizing radiation from their central sources. The data reveal high gas and dust temperatures (T > 100 K) extending out to about 5000 AU from each of these sources. There are no clear signs of disks or rotating structures down to our 1000 AU resolution. The extended warm gas provides evidence that, during the process of forming, these high-mass stars heat a large volume and correspondingly large mass of gas in their surroundings, inhibiting fragmentation and therefore keeping a large reservoir available to feed from. By contrast, the more mature massive stars that illuminate compact H ii regions have little effect on their surrounding dense gas, suggesting that these main sequence stars have completed most or all of their accretion. The high luminosity of the massive protostars (L > 10^4 solar luminosities), combined with a lack of centimeter continuum emission from these sources, implies that they are not on the main sequence while they accrete the majority of their mass; instead, they may be bloated and cool.

Einstein's General Theory of Relativity (GR) successfully describes gravity. The most fundamental predictions of GR are black holes (BHs), but in spite of many convincing BH candidates in the Universe, there is no conclusive experimental proof of their existence using astronomical observations in the electromagnetic spectrum. Are BHs real astrophysical objects? Does GR hold in its most extreme limit or are alternatives needed? The prime target to address these fundamental questions is in the center of our own Galaxy, which hosts the closest and best-constrained supermassive BH candidate in the Universe, Sagittarius A* (Sgr A*). Three different types of experiments hold the promise to test GR in a strong-field regime using observations of Sgr A* with new-generation instruments. The first experiment aims to image the relativistic plasma emission which surrounds the event horizon and forms a "shadow" cast against the background, whose predicted size (~50 microarcseconds) can now be resolved by upcoming VLBI experiments at mm-waves such as the Event Horizon Telescope (EHT). The second experiment aims to monitor stars orbiting Sgr A* with the upcoming near-infrared interferometer GRAVITY at the Very Large Telescope (VLT). The third experiment aims to time a radio pulsar in tight orbit about Sgr A* using radio telescopes (including the Atacama Large Millimeter Array or ALMA). The BlackHoleCam project exploits the synergy between these three different techniques and aims to measure the main BH parameters with sufficient precision to provide fundamental tests of GR and probe the spacetime around a BH in any metric theory of gravity. Here, we review our current knowledge of the physical properties of Sgr A* as well as the current status of such experimental efforts towards imaging the event horizon, measuring stellar orbits, and timing pulsars around Sgr A*.

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.

We present the results of a closure phase analysis of 3 mm very long baseline interferometry (VLBI) measurements performed on Sagittarius A* (Sgr A*). We have analyzed observations made in May 2015 using the Very Long Baseline Array, the Robert C. Byrd Green Bank Telescope and the Large Millimeter Telescope Alfonso Serrano and obtained non-zero closure phase measurements on several station triangles - indicative of a non-point-symmetric source structure. The data are fitted with an asymmetric source structure model in Sgr A*, represented by a simple two-component model, which favours a fainter component due East of the main source. This result is discussed in light of a scattering screen with substructure or an intrinsically asymmetric source.

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 present new JVLA observations of the high-mass cluster-forming region W51A from 2 to 16 GHz with resolution ${\theta}_{fwhm} \approx$ 0.3 - 0.5". The data reveal a wealth of observational results: (1) Currently-forming, very massive (proto-O) stars are traced by o-H2CO $2_{1,1}-2_{1,2}$ emission, suggesting that this line can be used efficiently as a massive protostar tracer. (2) There is a spatially distributed population of $\sim$mJy continuum sources, including hypercompact H ii regions and candidate colliding wind binaries, in and around the W51 proto-clusters. (3) There are two clearly detected protoclusters, W51e and W51 IRS2, that are gas-rich but may have most of their mass in stars within their inner $\sim$ 0.05 pc. The majority of the bolometric luminosity in W51 most likely comes from a third population of OB stars between these clusters. The presence of a substantial population of exposed O-stars coincident with a population of still-forming massive stars, along with a direct measurement of the low mass loss rate via ionized gas outflow from W51 IRS2, together imply that feedback is ineffective at halting star formation in massive protoclusters. Instead, feedback may shut off the large-scale accretion of diffuse gas onto the W51 protoclusters, implying that they are evolving towards a state of gas exhaustion rather than gas expulsion. Recent theoretical models predict gas exhaustion to be a necessary step in the formation of gravitationally bound stellar clusters, and our results provide an observational validation of this process.

This paper is the third in a series of ammonia multilevel imaging studies in well-known high-mass star forming regions. Using the JVLA, we have mapped the hot and dense molecular gas in W51 Main, with about 0.2 arcsec angular resolution, in five highly-excited metastable inversion transitions of ammonia (ammonia): (J,K)=(6,6), (7,7), (9,9), (10,10), and (13,13). We have identified and characterised two main centers of high-mass star formation in W51-Main: the W51e2 complex and the W51e8 core (6 arcsec southward of W51e2). The former breaks down into three further sub-cores: W51e2-W, which surrounds the well known HC HII region, where hot ammonia is observed in absorption, and two additional dusty cores, W51e2-E (~0.8" to the East) and W51e2-NW (~1" to the North), where hot ammonia is observed in emission. The velocity maps towards the HC HII region show a clear velocity gradient that may indicate rotation, though we do not directly observe a Keplerian velocity profile. The absence of outflow and/or maser activity and the low amount of molecular gas available for accretion (~5 solar masses) with respect to the mass of the central YSO (>20 solar masses), both indicate that the central YSO has already accreted most of its final mass. On the other hand, in the nearby W51e2-E object, the relatively large amount of hot molecular gas available for accretion (~20 solar masses, within about half an arcsecond or 2500 AU), along with strong outflow and maser activity, indicates that the main accretion center in the W51e2 complex is W51e2-E rather than W51e2-W. Finally, W51e2-NW and W51e8, although less dense, are also hot cores and contain a significant amount of molecular gas (~30 and 70 solar masses, respectively). We speculate that they may host high-mass YSOs either at a previous evolutionary stage or with lower mass than W51e2-E and W51e2-W.

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.

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.

We have used the JVLA at the 1 cm band to map five highly-excited metastable inversion transitions of ammonia, (J,K)=(6,6), (7,7), (9,9), (10,10), and (13,13), in W51 IRS2 with ~0.2 angular resolution. We present detections of both thermal (extended) ammonia emission in the five inversion lines, with rotational states ranging in energy from about 400 to 1700 K, and point-like ammonia maser emission in the (6,6), (7,7), and (9,9) lines. The thermal ammonia emits around a velocity of 60 km/s, near the cloud's systemic velocity, is elongated in the east-west direction across 4" and is confined by the HII regions W51d, W51d1, and W51d2. The ammonia masers are observed in the eastern tip of the dense clump traced by thermal ammonia, offset by 0.65" to the East from its emission peak, and have a peak velocity at ~47.5 km/s. No maser components are detected near the systemic velocity. The ammonia masers are separated by 0.65" (3500 AU) from the (rare) vibrationally-excited SiO masers, excited by the deeply-embedded YSO W51-North. This excludes that the two maser species are excited by the same object. Interestingly, the ammonia masers originate at the same sky position as a peak in a submm line of SO2 imaged with the SMA, tracing a face-on circumstellar disk/ring around W51-North. In addition, the thermal emission from the most highly excited ammonia lines, (10,10) and (13,13), shows two main condensations, the dominant one towards W51-North with the SiO/H2O masers, and a weaker peak at the ammonia maser position. We propose a scenario where the ring seen in SO2 emission is a circumbinary disk surrounding (at least) two high-mass YSOs, W51-North (exciting the SiO masers) and a nearby companion (exciting the ammonia masers), separated by 3500 AU. This finding indicates a physical connection (in a binary) between the two rare SiO and ammonia maser species.

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.

Photoevaporation due to high-energy stellar photons is thought to be one of the main drivers of protoplanetary disk dispersal. The fully or partially ionized disk surface is expected to produce free-free continuum emission at centimeter (cm) wavelengths that can be routinely detected with interferometers such as the upgraded Very Large Array (VLA). We use deep (rms noise down to 8 $\mu$Jy beam$^{-1}$ in the field of view center) 3.5 cm maps of the nearby (130 pc) Corona Australis (CrA) star formation (SF) region to constrain disk photoevaporation models. We find that the radio emission from disk sources in CrA is surprisingly faint. Only 3 out of 10 sources within the field of view are detected, with flux densities of order $10^2$ $\mu$Jy. However, a significant fraction of their emission is non-thermal. Typical upper limits for non-detections are $3\sigma\sim 60~\mu$Jy beam$^{-1}$. Assuming analytic expressions for the free-free emission from extreme-UV (EUV) irradiation, we derive stringent upper limits to the ionizing photon luminosity impinging on the disk surface $\Phi_\mathrm{EUV}<1-4\times10^{41}$ s$^{-1}$. These limits constrain $\Phi_\mathrm{EUV}$ to the low end of the values needed by EUV-driven photoevaporation models to clear protoplanetary disks in the observed few Myr timescale. Therefore, at least in CrA, EUV-driven photoevaporation is not likely to be the main agent of disk dispersal. We also compare the observed X-ray luminosities $L_X$ of disk sources with models in which photoevaporation is driven by such photons. Although predictions are less specific than for the EUV case, most of the observed fluxes (upper limits) are roughly consistent with the (scaled) predictions. Deeper observations, as well as predictions spanning a wider parameter space, are needed to properly test X-ray driven photoevaporation.

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.

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.

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.

An international consortium is presently constructing a beamformer for the Atacama Large Millimeter/submillimeter Array (ALMA) in Chile that will be available as a facility instrument. The beamformer will aggregate the entire collecting area of the array into a single, very large aperture. The extraordinary sensitivity of phased ALMA, combined with the extremely fine angular resolution available on baselines to the Northern Hemisphere, will enable transformational new very long baseline interferometry (VLBI) observations in Bands 6 and 7 (1.3 and 0.8 mm) and provide substantial improvements to existing VLBI arrays in Bands 1 and 3 (7 and 3 mm). The ALMA beamformer will have impact on a variety of scientific topics, including accretion and outflow processes around black holes in active galactic nuclei (AGN), tests of general relativity near black holes, jet launch and collimation from AGN and microquasars, pulsar and magnetar emission processes, the chemical history of the universe and the evolution of fundamental constants across cosmic time, maser science, and astrometry.

In Orion BN/KL, proper motions of 7 mm vibrationally-excited SiO masers trace rotation of a nearly edge-on disk and a bipolar wide-angle outflow 10-100 AU from radio Source I, a binary young stellar object (YSO) of ~20 Msun. Here we map ground-state 7 mm SiO emission with the Very Large Array and track proper motions over 9 years. The innermost and strongest emission lies in two extended arcs bracketing Source I. The proper motions trace a northeast-southwest bipolar outflow 100-1000 AU from Source I with a median 3D motion of ~18 km/s. An overlying distribution of 1.3 cm H2O masers betrays similar flow characteristics. Gas dynamics and emission morphology traced by the masers suggest the presence of a magnetocentrifugal disk-wind. Reinforcing evidence lies in the colinearity of the flow, apparent rotation across the flow parallel to the disk rotation, and recollimation that narrows the flow opening angle ~120 AU downstream. The arcs of ground-state SiO emission may mark the transition point to a shocked super-Alfvenic outflow.

Very long baseline interferometry (VLBI) observations of SiO masers in Orion Source I has enabled for the first time to resolve the outflow from a high-mass protostar in the launch and collimation region. Therefore, Source I provides a unique laboratory to study mass-loss and mass-accretion in a high-mass protostar. We numerically simulate the dynamics of the disk-wind inside 100 AU from Source I. This enables us to investigate the balance of different forces (gravitational, magnetic, thermal) regulating gas dynamics in massive star formation. In this work, we adopt magnetohydrodynamic (MHD) disk-wind models to explain the observed properties of the disk-wind from Orion Source I. The central source is assumed to be a binary composed of two 10\,$\msun$ stars in a circular orbit with an orbital separation of 7 AU. High resolution ideal MHD wind launching simulations (which prescribe disk as a boundary) are performed using the PLUTO code. The simulations are allowed to run until a steady state is obtained. MHD driven disk-wind provides a consistent model for the wide-angle flow from Source I probed by SiO masers, reproducing the bipolar morphology, the velocity amplitude and rotational profile, the physical conditions, and the magnetic field strength.

Compared to their centimeter-wavelength counterparts, millimeter recombination lines (RLs) are intrinsically brighter and are free of pressure broadening. We report observations of RLs (H30alpha at 231.9 GHz, H53alpha at 42.9 GHz) and the millimeter and centimeter continuum toward the Becklin-Neugebauer (BN) object in Orion, obtained from the Atacama Large Millimeter/submillimeter Array (ALMA) Science Verification archive and the Very Large Array (VLA). The RL emission appears to be arising from the slowly-moving, dense (Ne=8.4x10^6 cm^-3) base of the ionized envelope around BN. This ionized gas has a relatively low electron temperature (Te<4900 K) and small (<<10 km s^-1) bulk motions. Comparing our continuum measurements with previous (non)detections, it is possible that BN has large flux variations in the millimeter. However, dedicated observations with a uniform setup are needed to confirm this. From the H30alpha line, the central line-of-sight LSR velocity of BN is 26.3 km s^-1.

Using Science Verification data from the Atacama Large Millimeter/Submillimeter Array (ALMA), we have identified and imaged five rotational transitions (J=5-4 and J=6-5) of the three silicon monoxide isotopologues 28SiO v=0, 1, 2 and 29SiO v=0 and 28Si18O v=0 in the frequency range from 214 to 246 GHz towards the Orion BN/KL region. The emission of the ground-state 28SiO, 29SiO and 28Si18O shows an extended bipolar shape in the northeast-southwest direction at the position of Radio Source I, indicating that these isotopologues trace an outflow (~18 km/s, P.A. ~50deg, ~5000 AU in diameter) that is driven by this embedded high-mass young stellar object (YSO). Whereas on small scales (10-1000 AU) the outflow from Source I has a well-ordered spatial and velocity structure, as probed by Very Long Baseline Interferometry (VLBI) imaging of SiO masers, the large scales (500-5000 AU) probed by thermal SiO with ALMA reveal a complex structure and velocity field, most likely related to the effects of the environment of the BN/KL region on the outflow emanating from Source I. The emission of the vibrationally-excited species peaks at the position of Source I. This emission is compact and not resolved at an angular resolution of ~1.5" (~600 AU at a distance of 420 pc). 2-D Gaussian fitting to individual velocity channels locates emission peaks within radii of 100 AU, i.e. they trace the innermost part of the outflow. A narrow spectral profile and spatial distribution of the v=1 J=5-4 line similar to the masing v=1 J=1-0 transition, provide evidence for the most highly rotationally excited (frequency > 200 GHz) SiO maser emission associated with Source I known to date. The maser emission will enable studies of the Source I disk-outflow interface with future ALMA longest baselines.

We present new HST/WFPC2 observations of 20 fields centered around T Tauri stars in the Chamaeleon I star forming region. Images have been obtained in the F631N ([OI]6300A), F656N (Ha) and F673N ([SII]6716A+6731A) narrow-band filters, plus the Johnson V-band equivalent F547M filter. We detect 31 T Tauri stars falling within our fields. We discuss the optical morphology of 10 sources showing evidence of either binarity, circumstellar material, or mass loss. We supplement our photometry with a compilation of optical, infrared and sub-millimeter data from the literature, together with new sub-mm data for three objects, to build the Spectral Energy Distributions (SED) of 19 single sources. Using an SED model fitting tool, we self-consistently estimate a number of stellar and disk parameters, while mass accretion rates are directly derived from our Ha photometry. We find that bolometric luminosities derived from dereddened optical data tend to be underestimated in systems with high alpha(2-24} IR spectral index, suggesting that disks seen nearly edge-on may occasionally be interpreted as low luminosity (and therefore more evolved) sources. On the other hand, the same alpha(2-24) spectral index, a tracer of the amount of dust in the warmer layers of the circumstellar disks, and the mass accretion rate appear to decay with the isocronal stellar age, suggesting that the observed age spread (~0.5-5 Myr) within the cluster is real. Our sample contains a few outliers that may have dissipated their circumstellar disks on shorter time-scale.

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.

The Orion BN/KL complex is the nearest site of ongoing high-mass star formation. Recent proper motion observations provide convincing evidence of a recent (about 500 years ago) dynamical interaction between two massive young stellar objects in the region resulting in high velocities: the BN object and radio Source I. At the same time, Source I is surrounded by a nearly edge-on disc with radius ~50 au. These two observations taken together are puzzling: a dynamical encounter between multiple stars naturally yields the proper motions, but the survival of a disc is challenging to explain. In this paper we take the first steps to numerically explore the preferred dynamical scenario of Goddi et al., in which Source I is a binary that underwent a scattering encounter with BN, in order to determine if a pre-existing disc can survive this encounter in some form. Treating only gravitational forces, we are able to thoroughly and efficiently cover a large range of encounter parameters. We find that disc material can indeed survive a three-body scattering event if 1) the encounter is close, i.e. BN's closest approach to Source I is comparable to Source I's semi-major axis; and 2) the interplay of the three stars is of a short duration. Furthermore, we are able to constrain the initial conditions that can broadly produce the orientation of the present-day system's disc relative to its velocity vector. To first order we can thus confirm the plausibility of the scattering scenario of Goddi et al., and we have significantly constrained the parameters and narrowed the focus of future, more complex and expensive attempts to computationally model the complicated BN/KL region.

Using the Expanded Very Large Array, we have mapped the vicinity of the Orion BN/KL Hot Core with sub-arcsecond angular resolution in seven metastable inversion transitions of ammonia: (J,K)=(6,6) to (12,12). This emission comes from levels up to 1500 K above the ground state, enabling identification of source(s) responsible for heating the region. We used this multi-transition dataset to produce images of the rotational/kinetic temperature and the column density of ammonia for ortho and para species separately and on a position-by-position basis. We find rotational temperature and column density in the range 160-490 K and (1-4)x10^17 cm^-2, respectively. Our spatially-resolved images show that the highest (column) density and hottest gas is found in a northeast-southwest elongated ridge to the southeast of Source I. We have also measured the ortho-para ratio of ammonia, estimated to vary in the range 0.9-1.6. Enhancement of ortho with respect to para and the offset of hot ammonia emission peaks from known (proto)stellar sources provide evidence that the ammonia molecules have been released from dust grains into the gas-phase through the passage of shocks and not by stellar radiation. We propose that the combined effect of Source I's proper motion and its low-velocity outflow impinging on a pre-existing dense medium is responsible for the excitation of ammonia and the Orion Hot Core. Finally, we found for the first time evidence of a slow (5 km/s) and compact (1000 AU) outflow towards IRc7.

We present new 7mm continuum observations of Orion BN/KL with the VLA. We resolve the emission from the protostar radio Source I and BN at several epochs. Source I is highly elongated NW-SE, and remarkably stable in flux density, position angle, and overall morphology over nearly a decade. This favors the extended emission component arising from an ionized disk rather than a jet. We have measured the proper motions of Source I and BN for the first time at 43 GHz. We confirm that both sources are moving at high speed (12 and 26 km/s, respectively) approximately in opposite directions, as previously inferred from measurements at lower frequencies. We discuss dynamical scenarios that can explain the large motions of both BN and Source I and the presence of disks around both. Our new measurements support the hypothesis that a close (~50 AU) dynamical interaction occurred around 500 years ago between Source I and BN as proposed by Gomez et al. From the dynamics of encounter we argue that Source I today is likely to be a binary with a total mass on the order of 20 Msun, and that it probably existed as a softer binary before the close encounter. This enables preservation of the original accretion disk, though truncated to its present radius of ~50 AU. N-body numerical simulations show that the dynamical interaction between a binary of 20 Msun total mass (I) and a single star of 10 Msun mass (BN) may lead to the ejection of both and binary hardening. The gravitational energy released in the process would be large enough to power the wide-angle flow traced by H2 and CO emission in the BN/KL nebula. Assuming the proposed dynamical history is correct, the smaller mass for Source I recently estimated from SiO maser dynamics (>7 Msun) by Matthews et al., suggests that non-gravitational forces (e.g. magnetic) must play an important role in the circumstellar gas dynamics.

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)

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 present multi-epoch VLBA imaging of the 28SiO v=1 & v=2, J=1-0 maser emission toward the massive YSO Orion Source I. Both SiO transitions were observed simultaneously with an angular resolution of ~0.5 mas (~0.2 AU for d=414 pc). Here we explore the global properties and kinematics of the emission through two 19-epoch animated movies spanning 21 months (2001 March 19 to 2002 December 10). These movies provide the most detailed view to date of the dynamics and temporal evolution of molecular material within ~20-100 AU of a massive (~>8M_sun) YSO. The bulk of the SiO masers surrounding Source I lie in an X-shaped locus; emission in the South/East arms is predominantly blueshifted and in the North and West is predominantly redshifted. In addition, bridges of intermediate-velocity emission connect the red and blue sides of the emission distribution. We have measured proper motions of >1000 maser features and find a combination of radially outward migrations along the four arms and motions tangent to the bridges. We interpret the SiO masers as arising from a wide-angle bipolar wind emanating from a rotating, edge-on disk. The detection of maser features along extended, curved filaments suggests that magnetic fields may play a role in launching and/or shaping the wind. Our observations appear to support a picture in which stars with M ~>8 M_sun form via disk-mediated accretion. However, we cannot rule out that the Source I disk may have been formed or altered following a close encounter. (Abridged).

We have used the Very Large Array (VLA) at 7 mm wavelength to image five rotational transitions (J=1-0) from three SiO isotopologues towards Orion BN/KL: 28SiO v=0,1,2; and 29SiO and 30SiO v=0. For the first time, we have mapped the 29SiO and 30SiO J=1-0 emission, established the maser nature of the emission, and confirmed association with the deeply embedded high-mass young stellar object commonly denoted radio Source I. The 28SiO v=0 maser emission shows a bipolar structure that extends over ~700 AU along a northeast-southwest axis, and we propose that it traces a bipolar outflow driven by Source I. The high-brightness isotopic SiO maser emission imaged with a ~0.2 arcsec resolution has a more compact distribution, generally similar to that of the 28SiO v=1,2 emission, and it probably traces bulk gas flows in a region of diameter <100 AU centered on Source I. On small scales of <10 AU, however, compact 29SiO/30SiO v=0 and 28SiO v=1,2 emission features may be offset from one another in position and line-of-sight velocity. From a radiative transfer analysis based on a large velocity gradient (LVG) pumping model, we derive similar temperatures and densities for the optimum excitation of both 29SiO/30SiO v=0 and 28SiO v=1,2 masers, significantly higher than required for 28SiO v=0 maser excitation. In order to account for the small-scale differences among the isotopologues (v=0) and the main species (v=1,2), follow-up radiative transfer modeling that incorporates non-local line overlap among transitions of all SiO isotopic species may be required.

We have surveyed molecular line emission from Orion BN/KL from 42.3 to 43.6 GHz with the Green Bank Telescope. Sixty-seven lines were identified and ascribed to 13 different molecular species. The spectrum at 7 mm is dominated by SiO, SO2, CH3OCH3, and C2H5CN. Five transitions have been detected from the SiO isotopologues 28SiO, 29SiO, and 30SiO. We report here for the first time the spectra of the 29SiO and 30SiO v=0 J=1-0 emission in Orion BN/KL, and we show that they have double-peaked profiles with velocity extents similar to the main isotopologue. The main motivation for the survey was the search of high-velocity (100-1000 km/s) outflows in the BN/KL region as traced by SiO Doppler components. Some of the unidentified lines in principle could be high-velocity SiO features, but without imaging data their location cannot be established. Wings of emission are present in the v=0 28SiO, 29SiO and 30SiO profiles, and we suggest that the v=0 emission from the three isotopologues might trace a moderately high-velocity (~30-50 km/s) component of the flows around the high-mass protostar Source I in the Orion BN/KL region. We also confirm the 7 mm detection of a complex oxygen-bearing species, acetone (CH3COCH3), which has been recently observed towards the hot core at 3 mm, and we have found further indications of the presence of long cyanopolyynes (HC5N and HC7N) in the quiescent cold gas of the extended ridge.

A comprehensive picture of high-mass star formation has remained elusive, in part because examples of high-mass YSOs tend to be relatively distant, deeply embedded, and confused with other emission sources. These factors have impeded dynamical investigations within tens of AU of high-mass YSOs--scales that are critical for probing the interfaces where outflows from accretion disks are launched and collimated. Using observations of SiO masers obtained with the VLA and the VLBA, the KaLYPSO project is overcoming these limitations by mapping the structure and dynamical/temporal evolution of the material 10-1000 AU from the nearest high-mass YSO: Radio Source I in the Orion BN/KL region. Our data include ~40 epochs of VLBA observations over a several-year period, allowing us to track the proper motions of individual SiO maser spots and to monitor changes in the physical conditions of the emitting material with time. Ultimately these data will provide 3-D maps of the outflow structure over approximately 30% of the outflow crossing time. Here we summarize recent results from the KaLYPSO project, including evidence that high-mass star formation is occurring via disk-mediated accretion.

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.

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