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.

We model the ALMA and VLA millimeter radial profiles of the disk around HL Tau to constrain the properties of the dust grains. We adopt the disk evolutionary models of Lynden-Bell \& Pringle and calculate their temperature and density structure and emission. These disks are heated by the internal viscosity and irradiated by the central star and a warm envelope. We consider a dust size distribution $n(a) da \propto a^{-3.5} da $, and vary the maximum grain size in the atmosphere and the midplane, $a_{\rm max}=100\ \mu$m, 1 mm, and 1cm. We also include dust settling and vary the dust-to-gas mass ratio from 1 to 9 times the ISM value. We find that the models that can fit the observed level of emission along the profiles at all wavelengths have an atmosphere with a maximum grain size $a_{\rm max} = 100 \ \mu$m, and a midplane with $a_{\rm max}=1$ cm. The disk substructure, with a deficit of emission in the gaps, can be due to dust properties in these regions that are different from those in the rings. We test an opacity effect (different $a_{\rm max}$) and a dust mass deficit (smaller dust-to-gas mass ratio) in the gaps. We find that the emission profiles are better reproduced by models with a dust deficit in the gaps, although a combined effect is also possible. These models have a global dust-to-gas mass ratio twice the ISM value, needed to reach the level of emission of the 7.8 mm VLA profile.

Understanding planet formation requires to discern how dust grows in protoplanetary disks. An important parameter to measure in disks is the maximum dust grain size present. This is usually estimated through measurements of the dust opacity at different millimeter wavelengths assuming optically thin emission and dust opacity dominated by absorption. However, ALMA observations have shown that these assumptions might not be correct in the case of protoplanetary disks, leading to overestimation of particle sizes and to underestimation of the disk's mass. Here, we present an analysis of high quality ALMA and VLA images of the HL Tau protoplanetary disk, covering a wide range of wavelengths, from 0.8 mm to 1 cm, and with a physical resolution of $\sim$7.35 au. We describe a procedure to analyze a set of millimeter images without any assumption about the optical depth of the emission, and including the effects of absorption and scattering in the dust opacity. This procedure allows us to obtain the dust temperature, the dust surface density and the maximum particle size at each radius. In the HL Tau disk, we found that particles have already grown up to a few millimeters in size. We detect differences in the dust properties between dark and bright rings, with dark rings containing low dust density and small dust particles. Different features in the HL Tau disk seem to have different origins. Planet-disk interactions can explain substructure at the external half of the disk, but the internal rings seem to be associated to the presence of snow lines of several molecules.

In this study we analyze 70 radio continuum sources associated with dust clumps and considered to be candidates for the earliest stages of high-mass star formation. The detection of these sources was reported by Rosero et al. (2016), who found most of them to show weak (${\scriptstyle <}$1 mJy) and compact (${\scriptstyle <}\,$0.6$^{\prime \prime}$) radio emission. Herein, we used the observed parameters of these sources to investigate the origin of the radio continuum emission. We found that at least $\sim 30\%$ of these radio detections are most likely ionized jets associated with high-mass protostars, but for the most compact sources we cannot discard the scenario that they represent pressure-confined HII regions. This result is highly relevant for recent theoretical models based on core accretion that predict the first stages of ionization from high-mass stars to be in the form of jets. Additionally, we found that properties such as the radio luminosity as a function of the bolometric luminosity of ionized jets from low and high-mass stars are extremely well-correlated. Our data improve upon previous studies by providing further evidence of a common origin for jets independently of luminosity.

In recent years, the disk populations in a number of young star-forming regions have been surveyed with ALMA. Understanding the disk properties and their correlation with those of the central star is critical to understand planet formation. In particular, a decrease of the average measured disk dust mass with the age of the region has been observed. We conducted high-sensitivity continuum ALMA observations of 43 Class II young stellar objects in CrA at 1.3 mm (230 GHz). The typical spatial resolution is 0.3". The continuum fluxes are used to estimate the dust masses of the disks, and a survival analysis is performed to estimate the average dust mass. We also obtained new VLT/X-Shooter spectra for 12 of the objects in our sample. 24 disks are detected, and stringent limits have been put on the average dust mass of the non-detections. Accounting for the upper limits, the average disk mass in CrA is $6\pm3\,\rm M_\oplus$, significantly lower than that of disks in other young (1-3 Myr) star forming regions (e.g. Lupus) and appears consistent with the 5-10 Myr old Upper Sco. The position of the stars in our sample on the HR diagram, however, seems to confirm that that CrA has age similar to Lupus. Neither external photoevaporation nor a lower than usual stellar mass distribution can explain the low disk masses. On the other hand, a low-mass disk population could be explained if the disks are small, which could happen if the parent cloud has a low temperature or intrinsic angular momentum, or if the the angular momentum of the cloud is removed by some physical mechanism such as magnetic braking. In order to fully explain and understand the dust mass distribution of protoplanetary disks and their evolution, it may also be necessary to take into consideration the initial conditions of star and disk formation process, which may vary from region to region, and affect planet formation.

We study dust concentration in axisymmetric gas rings in protoplanetary disks. Given the gas surface density, we derived an analytical total dust surface density by taking into account the differential concentration of all the grain sizes. This model allows us to predict the local dust-to-gas mass ratio and the slope of the particle size distribution, as a function of radius. We test this analytical model comparing it with a 3D magneto-hydrodynamical simulation of dust evolution in an accretion disk. The model is also applied to the disk around HD 169142. By fitting the disk continuum observations simultaneously at $\lambda = 0.87$, 1.3, 3.0 mm, we obtain a global dust-to-gas mass ratio $\epsilon_{\rm global} = 1.05 \times 10^{-2}$ and a viscosity coefficient $\alpha = 1.35 \times 10^{-2}$. This model can be easily implemented in numerical simulations of accretion disks.

Polarized dust emission outside of disks reveal the magnetic field morphology of molecular clouds. Within disks, however, polarized dust emission can arise from very different mechanisms (e.g., self-scattering), and each of them are useful for constraining physical properties in the disk. For example, these mechanisms allow us to constrain the disk grain size distributions and grain/disk geometries, independent from current methods of measuring these parameters. To accurately model these features and disentangle the various polarization mechanisms, multiwavelength observations at very high resolution and sensitivity are required. With significant upgrades to current interferometric facilities, we can understand how grains evolve in disks during the planet formation process.

We used ALMA to observe the star-forming region GGD27 at 1.14 mm with an unprecedented angular resolution, 40 mas (56 au) and sensitivity (0.002 Msun). We detected a cluster of 25 continuum sources, most of which are likely tracing disks around Class 0/I protostars. Excluding the two most massive objects, disks masses are in the range 0.003-0.05 Msun. The analysis of the cluster properties indicates that GGD27 displays moderate subclustering. This result combined with the dynamical timescale of the radio jet (10000 years) suggests the youthfulness of the cluster. The lack of disk mass segregation signatures may support this too. We found a clear paucity of disks with Rdisk >100 au. The median value of the radius is 34 au, smaller than the median of 92 au for Taurus but comparable to the value found in Ophiuchus and in the Orion Nebula Cluster. In GGD27 there is no evidence of a distance-dependent disk mass distribution (i. e., disk mass depletion due to external photoevaporation), most likely due to the cluster youth. There is a clear deficit of disks for distances <0.02 pc. Only for distances >0.04 pc stars can form larger and more massive disks, suggesting that dynamical interactions far from the cluster center are weaker, although the small disks found could be the result of disk truncation. This work demonstrates the potential to characterize disks from low-mass YSOs in distant and massive (still deeply embedded) clustered environments.

We present an analysis of radio (Karl G. Jansky Very Large Array (VLA)), optical (HST), and X-ray (Chandra and XMM-Newton) observations and archival data of the Herbig-Haro objects HH 80 and HH 81 in the context of jet-cloud interactions. Our radio images are the highest angular resolution to date of these objects, allowing to spatially resolve the knots and compare the regions emitting in the different spectral ranges. We found that soft X-ray thermal emission is located ahead of the non-thermal radio peak. This result is consistent with a radiative forward shock that heats the shocked gas up to 10 6 K, and an adiabatic reverse shock able to accelerate particles and produce synchrotron radiation detectable at radio frequencies. These high angular resolution radio images also reveal a bow shock structure in the case of HH 80N, being the first time this morphology is detected in a Herbig-Haro object at these frequencies.

In this paper we explore the effects of self-obscuration in protostellar disks with a radially decreasing temperature gradient and a colder midplane. We are motivated by recent reports of resolved dark lanes (`hamburgers') and (sub)mm spectral indices systematically below the ISM value for optically thin dust $\alpha_{\rm ISM} =3.7$. We explore several model grids, scaling disk mass and varying inclination angle $i$ and observing frequency $\nu$ from the VLA Ka band ($\sim 37$ GHz) to ALMA Band 8 ($\sim 405$ GHz). We also consider the effects of decreasing the index of the (sub-)mm dust opacity power law $\beta$ from 1.7 to 1. We find that a distribution of disk masses in the range $M_{\rm disk} = 0.01-2~M_\odot$ is needed to reproduce the observed distribution of spectral indices, and that assuming a fixed $\beta =1.7$ gives better results than $\beta=1$. A wide distribution of disk masses is also needed to produce some cases with $\alpha <2$, as reported for some sources in the literature. Such extremely low spectral indices arise naturally when the selected observing frequencies sample the appropriate change in the temperature structure of the optically thick model disk. Our results show that protostellar disk masses could often be underestimated by $> \times10$, and are consistent with recent hydrodynamical simulations. Although we do not rule out the possibility of some grain growth occurring within the short protostellar timescales, we conclude that self-obscuration needs to be taken into account.

In this contribution we explore the new science that a Next Generation Very Large Array (ngVLA) would be able to perform on the topics of the youngest HII regions and (proto)stellar jets. Free-free continuum and radio recombination line (RRL) emission are often the only way of peering into the dense envelopes surrounding (proto)stars of all masses, and trace their initial feedback in the form of `radio jets', `hypercompact HII regions', or photoevaporating, partially-ionized flows. Properly disentangling free-free from dust emission is also mandatory in studies of protoplanetary and accretion disks. Current VLA research has reached an impasse in which a population of faint ionized radio sources, probably corresponding to the above mentioned objects, is detected, but their nature is mostly unknown. The ngVLA would allow us to resolve the density structure and kinematics of such sources, revolutionizing our knowledge of star formation across the entire stellar-mass spectrum.

The magnetic field plays an important role in every stage of the star-formation process from the collapse of the initial protostellar core to the star's arrival on the main sequence. Consequently, the goal of this science case is to explore a wide range of magnetic phenomena that can be investigated using the polarization capabilities of the Next Generation Very Large Array (ngVLA). These include (1) magnetic fields in protostellar cores via polarized emission from aligned dust grains, including in regions optically thick at wavelengths observable by the Atacama Large Millimeter/submillimeter Array (ALMA); (2) magnetic fields in both protostellar cores and molecular outflows via spectral-line polarization from the Zeeman and Goldreich-Kylafis effects; (3) magnetic fields in protostellar jets via polarized synchrotron emission; and (4) gyrosynchrotron emission from magnetospheres around low-mass stars.

We present broadband polarimetric observations of a sample of high Faraday rotation measure (RM) AGN using the VLA telescope from 1 to 2 GHz, and 4 to 12 GHz. The sample (14 sources) consists of compact sources (linear resolution smaller than ~ 5 kpc) that are unpolarized at 1.4 GHz in the NVSS. Total intensity data have been modelled using combination of synchrotron components, revealing complex structure in their radio spectra. Depolarization modelling, through the so called qu-fitting, have been performed on the polarized data using an equation that attempts to simplify the process of fitting many different depolarization models that we can divide into two major categories: External Depolarization and Internal Depolarization models. Understanding which of the two mechanisms are the most representative, would help the qualitative understanding of the AGN jet environment, whether it is embedded in a dense external magnetoionic medium or if it is the jet-wind that causes the high RM and strong depolarization. This could help to probe the jet magnetic field geometry (e.g. helical or otherwise). This new high-sensitivity data, shows a complicated behaviour in the total intensity and polarization radio spectrum of individual sources. We observed the presence of several synchrotron components and Faraday components in their total intensity and polarized spectra. For the majority of our targets, (12 sources) the depolarization seems to be caused by a turbulent magnetic field. Thus, our main selection criteria (lack of polarization at 1.4 GHz in the NVSS), results in a sample of sources with very large RMs and depolarization due to turbulent magnetic fields local to the source. We show how the new qu-fitting technique can be used to probe the magnetised radio source environment and to spectrally resolve the polarized components of unresolved radio sources.

Radio emission from protostellar jets is usually dominated by free-free emission from thermal electrons. However, in some cases, it has been proposed that non-thermal emission could also be present. This additional contribution from non-thermal emission has been inferred through negative spectral indices at centimeter wavelengths in some regions of the radio jets. In the case of HH 80-81, one of the most powerful protostellar jets known, linearly polarized emission has also been detected, revealing that the non-thermal emission is of synchrotron nature from a population of relativistic particles in the jet. This result implies that an acceleration mechanism should be taking place in some parts of the jet. Here, we present new high sensitivity and high angular resolution radio observations at several wavelengths (in the 3-20 cm range) of the HH80-81 radio jet. These new observations represent an improvement in sensitivity and angular resolution by a factor of $\sim$10 with respect to previous observations. This allows us to resolve the morphology of the radio jet, and to study the different emission mechanisms involved through spectral index maps. We conclude that synchrotron emission in this jet arises from an extended component detected at low frequencies and from the termination points of the jet, where strong shocks against the ambient medium can produce efficient particle acceleration.

Young stellar objects (YSOs) may undergo periods of active accretion (outbursts), during which the protostellar accretion rate is temporarily enhanced by a few orders of magnitude. Whether or not these accretion outburst YSOs possess similar dust/gas reservoirs to each other, and whether or not their dust/gas reservoirs are similar as quiescent YSOs, are issues not yet clarified. The aim of this work is to characterize the millimeter thermal dust emission properties of a statistically significant sample of long and short duration accretion outburst YSOs (i.e., FUors and EXors) and the spectroscopically identified candidates of accretion outbursting YSOs (i.e., FUor-like objects). We have carried out extensive Submillimeter Array (SMA) observations mostly at $\sim$225 GHz (1.33 mm) and $\sim$272 GHz (1.10 mm), from 2008 to 2017. We covered accretion outburst YSOs located at $<$1 kpc distances from the solar system. We analyze all the existing SMA data of such objects, both published and unpublished, in a coherent way to present a millimeter interferometric database of 29 objects. We obtained 21 detections at $>$3-$\sigma$ significance. Detected sources except for the two cases of V883 Ori and NGC 2071 MM3 were observed with $\sim$1$"$ angular resolution. Overall our observed targets show a systematically higher millimeter luminosity distribution than those of the $M_{*}>$0.3 $M_{\odot}$ Class II YSOs in the nearby ($\lesssim$400 pc) low-mass star-forming molecular clouds (e.g., Taurus, Lupus, Upp Scorpio, and Chameleon I). In addition, at 1 mm our observed confirmed binaries or triple-system sources are systematically fainter than the rest of the sources even though their 1 mm fluxes are broadly distributed. We may have detected $\sim$30-60\% millimeter flux variability from V2494 Cyg and V2495 Cyg, from the observations separated by $\sim$1 year.

We present preliminary results of JVLA wideband full polarization observations of a sample of Active Galactic Nuclei (AGN) with very high Rotation Measure (RM) values, a sign of extreme environment. Polarization properties show a complex behaviour such that the polarization angle (PA) and fractional polarization (fp) change dramatically within the wide band. The measured RM is not constant within the wide band. Its complex behaviour reflects the complexity of the medium with the presence of several Faraday components. The depolarization has been studied by modelling the variations of the Stokes parameters Q and U together with the polarization parameters (PA and fp) with wavelength using combinations of the simplest existing depolarization models. With this JVLA study we could spectrally resolve multiple polarized components of unresolved AGN. These preliminary results reveal the complexity of these objects, but improvements to the depolarization modelling are needed to better understand the polarization structure of these sources.

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 a high sensitivity radio continuum survey at 6 and 1.3$\,$cm using the Karl G. Jansky Very Large Array towards a sample of 58 high-mass star forming regions. Our sample was chosen from dust clumps within infrared dark clouds with and without IR sources (CMC-IRs, CMCs, respectively), and hot molecular cores (HMCs), with no previous, or relatively weak radio continuum detection at the $1\,$mJy level. Due to the improvement in the continuum sensitivity of the VLA, this survey achieved map rms levels of $\sim$ 3-10 $\mu$Jy beam$^{-1}$ at sub-arcsecond angular resolution. We extracted 70 centimeter continuum sources associated with 1.2$\,$mm dust clumps. Most sources are weak, compact, and are prime candidates for high-mass protostars. Detection rates of radio sources associated with the mm dust clumps for CMCs, CMC-IRs and HMCs are 6$\%$, 53$\%$ and 100$\%$, respectively. This result is consistent with increasing high-mass star formation activity from CMCs to HMCs. The radio sources located within HMCs and CMC-IRs occur close to the dust clump centers with a median offset from it of 12,000$\,$AU and 4,000$\,$AU, respectively. We calculated 5 - 25$\,$GHz spectral indices using power law fits and obtain a median value of 0.5 (i.e., flux increasing with frequency), suggestive of thermal emission from ionized jets. In this paper we describe the sample, observations, and detections. The analysis and discussion will be presented in Paper II.

We present Submillimeter Array (SMA) 1.35 mm subarcsecond angular resolution observations toward the LkH{\alpha} 234 intermediate-mass star-forming region. The dust emission arises from a filamentary structure of $\sim$5 arcsec ($\sim$4500 au) enclosing VLA 1-3 and MM 1, perpendicular to the different outflows detected in the region. The most evolved objects are located at the southeastern edge of the dust filamentary structure and the youngest ones at the northeastern edge. The circumstellar structures around VLA 1, VLA 3, and MM 1 have radii between $\sim$200 and $\sim$375 au and masses in the $\sim$0.08-0.3 M$_{\odot}$ range. The 1.35 mm emission of VLA 2 arises from an unresolved (r$< 135$ au) circumstellar disk with a mass of $\sim$0.02 M$_{\odot}$. This source is powering a compact ($\sim$4000 au), low radial velocity ($\sim$7 km s$^{-1}$) SiO bipolar outflow, close to the plane of the sky. We conclude that this outflow is the "large-scale" counterpart of the short-lived, episodic, bipolar outflow observed through H$_2$O masers at much smaller scales ($\sim $180 au), and that has been created by the accumulation of the ejection of several episodic collimated events of material. The circumstellar gas around VLA 2 and VLA 3 is hot ($\sim$130 K) and exhibits velocity gradients that could trace rotation. There is a bridge of warm and dense molecular gas connecting VLA 2 and VLA 3. We discuss the possibility that this bridge could trace a stream of gas between VLA 3 and VLA 2, increasing the accretion rate onto VLA 2 to explain why this source has an important outflow activity.

We report ALMA observations of a one-sided, high-velocity ($\sim$80 km s$^{-1}$) CO($J = 2 \rightarrow 1$) jet powered by the intermediate-mass protostellar source Serpens SMM1-a. The highly collimated molecular jet is flanked at the base by a wide-angle cavity; the walls of the cavity can be seen in both 4 cm free-free emission detected by the VLA and 1.3 mm thermal dust emission detected by ALMA. This is the first time that ionization of an outflow cavity has been directly detected via free-free emission in a very young, embedded Class 0 protostellar source that is still powering a molecular jet. The cavity walls are ionized either by UV photons escaping from the accreting protostellar source, or by the precessing molecular jet impacting the walls. These observations suggest that ionized outflow cavities may be common in Class 0 protostellar sources, shedding further light on the radiation, outflow, and jet environments in the youngest, most embedded forming stars.

We report new Karl G. Jansky Very Large Array (JVLA), 0$"$.5 angular resolution observations of linearly polarized continuum emission at 6.9 mm, towards the Class 0 young stellar object (YSO) NGC1333 IRAS4A. This target source is a collapsing dense molecular core, which was resolved at short wavelengths to have hourglass shaped B-field configuration. We compare these 6.9 mm observations with previous polarization Submillimeter Array (SMA) observations at 0.88 mm, which have comparable angular resolution ($\sim$0$"$7). We found that at the same resolution, the observed polarization position angles at 6.9 mm are slightly deviated from those observed at 0.88 mm. Due to the lower optical depth of the emission at 6.9 mm, and the potential effect of dust grain growth, the new JVLA observations are likely probing B-field alignments in regions interior to those sampled by the previous polarization observations at higher frequencies. Our understanding can be improved by more sensitive observations, and observations for the more extended spatial scales.

While most protostellar jets present free-free emission at radio wavelengths, synchrotron emission has been also proposed to be present in a handful of these objects. The presence of non-thermal emission has been inferred by negative spectral indices at centimeter wavelengths. In one case (the HH 80-81 jet arising from a massive protostar), its synchrotron nature was confirmed by the detection of linearly polarized radio emission. One of the main consequences of these results is that synchrotron emission implies the presence of relativistic particles among the non-relativistic material of these jets. Therefore, an acceleration mechanism should be taking place. The most probable scenario is that particles are accelerated when the jets strongly impact against the dense envelope surrounding the protostar. Here, we present an analysis of radio observations obtained with the Very Large Array of the Triple Radio Source in the Serpens star-forming region. This object is known to be a radio jet arising from an intermediate-mass protostar. It is also one of the first protostellar jets where the presence of non-thermal emission was proposed. We analysed the dynamics of the jet as well as the nature of the emission and discuss these issues in the context of the physical parameters of the jet and the particle acceleration phenomenon.

Synchrotron radiation is commonly observed in connection with shocks of different velocities, ranging from relativistic shocks associated with active galactic nuclei, gamma-ray bursts or microquasars to weakly- or non-relativistic flows as those observed in supernova remnants. Recent observations of synchrotron emission in protostellar jets are important not only because they extend the range over which the acceleration process works, but also because they allow to determine the jet and/or interstellar magnetic field structure, thus giving insights on the jet ejection and collimation mechanisms. In this paper, we compute for the first time polarized (synchrotron) and non polarized (thermal X-ray) synthetic emission maps from axisymmetrical simulations of magnetized protostellar jets. We consider models with different jet velocities and variability, as well as toroidal or helical magnetic field. Our simulations show that variable, low-density jets with velocities $\sim$ 1000 km s$^{-1}$ and $\sim$ 10 times lighter than the environment can produce internal knots with significant synchrotron emission, and thermal X-rays in the shocked region of the leading bow shock moving in a dense medium. While models with a purely toroidal magnetic field show a very large degree of polarization, models with helical magnetic field show lower values and a decrease of the degree of polarization, in agreement with observations of protostellar jets.

(abridged) The HH 80/81/80N jet extends from the HH 80 object to the recently discovered Source 34 and has a total projected jet size of 10.3 pc, constituting the largest collimated radio-jet system known so far. It is powered by IRAS 18162-2048 associated with a massive young stellar object. We report 6 cm JVLA observations that, compared with previous 6 cm VLA observations carried out in 1989, allow us to derive proper motions of the HH 80, HH 81 and HH 80N radio knots located about 2.5 pc away in projection from the powering source. For the first time, we measure proper motions of the optically obscured HH 80N object providing evidence that HH 81, 80 and 80N are associated with the same radio-jet. We derived tangential velocities of these HH objects between 260 and 350 km/s, significantly lower than those for the radio knots of the jet close to the powering source (600-1400 km/s) derived in a previous work, suggesting that the jet material is slowing down due to a strong interaction with the ambient medium. The HH 80 and HH 80N emission at 6 cm is, at least in part, probably synchrotron radiation produced by relativistic electrons in a magnetic field of 1 mG. If these electrons are accelerated in a reverse adiabatic shock, we estimate a jet total density of $\lesssim1000$ cm$^{-3}$. All these features are consistent with a jet emanating from a high mass protostar and make evident its capability of accelerating particles up to relativistic velocities.

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.

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

We report 1.3 cm and 6 cm continuum observations toward the massive proto-stellar candidate G11.11-0.12P1 using the Karl G. Jansky Very Large Array (VLA). We detect a string of four unresolved radio continuum sources coincident with the mid-IR source in G11P1. The continuum sources have positive spectral indices consistent with a thermal (free-free) ionized jet. The most likely origin of the ionized gas are shocks due to the interaction of a stellar wind with the surrounding high-density material. We also present NIR United Kingdom Infrared Telescope (UKIRT) archival data which show an extended structure detected only at K-band (2.2 $\mu$m), which is oriented perpendicular to the jet, and that may be scattered light from a circumstellar disk around the massive protostar. Our observations plus the UKIRT archival data thus provide new evidence that a disk/jet system is present in the massive protostellar candidate located in the G11.11-0.12P1 core.

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.

We report multi-epoch VLBI H$_2$O maser observations towards the compact cluster of YSOs close to the Herbig Be star LkH$\alpha$ 234. This cluster includes LkH$\alpha$ 234 and at least nine more YSOs that are formed within projected distances of $\sim$10 arcsec ($\sim$9,000 au). We detect H$_2$O maser emission towards four of these YSOs. In particular, our VLBI observations (including proper motion measurements) reveal a remarkable very compact ($\sim$0.2 arcsec = $\sim$180 au), bipolar H$_2$O maser outflow emerging from the embedded YSO VLA 2. We estimate a kinematic age of $\sim$40 yr for this bipolar outflow, with expanding velocities of $\sim$20 km s$^{-1}$ and momentum rate $\dot M_w V_w$ $\simeq$ $10^{-4}-10^{-3}$ M$_{\odot}$ yr$^{-1}$ km s$^{-1}$$\times (\Omega$/$4\pi)$, powered by a YSO of a few solar masses. We propose that the outflow is produced by recurrent episodic jet ejections associated with the formation of this YSO. Short-lived episodic ejection events have previously been found towards high-mass YSOs. We show now that this behaviour is also present in intermediate-mass YSOs. These short-lived episodic ejections are probably related to episodic increases in the accretion rate, as observed in low-mass YSOs. We predict the presence of an accretion disk associated with VLA 2. If detected, this would represent one of the few known examples of intermediate-mass stars with a disk-YSO-jet system at scales of a few hundred au.

We present multi-epoch Very Long Baseline Array (VLBA) H$_2$O maser observations toward the massive young stellar objects (YSOs) VLA 2 and VLA 3 in the star-forming region AFGL 2591. Through these observations, we have extended the study of the evolution of the masers towards these objects up to a time span of $\sim$ 10 yrs, measuring their radial velocities and proper motions. The H$_2$O masers in VLA 3, the most massive YSO in AFGL 2591 ($\sim$ 30--40~M$_{\odot}$), are grouped within projected distances of $\lesssim$ 40 mas ($\lesssim$ 130 AU) from VLA 3. In contrast to other H$_2$O masers in AFGL 2591, the masers associated with VLA 3 are significantly blueshifted (up to $\sim$ 30 km s$^{-1}$) with respect to the velocity of the ambient molecular cloud. We find that the H$_2$O maser cluster as a whole, has moved westwards of VLA~3 between the 2001 and 2009 observations, with a proper motion of $\sim$ 1.2 mas yr$^{-1}$ ($\sim$ 20 km s$^{-1}$). We conclude that these masers are tracing blueshifted outflowing material, shock excited at the inner parts of a cavity seen previously in ammonia molecular lines and infrared images, and proposed to be evacuated by the outflow associated with the massive VLA 3 source. The masers in the region of VLA 2 are located at projected distances of $\sim$ 0.7$''$ ($\sim$ 2300 AU) north from this source, with their kinematics suggesting that they are excited by a YSO other than VLA 2. This driving source has not yet been identified.

In this paper we analyze multi-epoch Very Long Baseline Interferometry (VLBI) water maser observations carried out with the Very Long Baseline Array (VLBA) toward the high-mass star-forming region AFGL 2591. We detected maser emission associated with the radio continuum sources VLA 2 and VLA 3. In addition, a water maser cluster, VLA 3-N, was detected ~ 0.5" north of VLA 3. We concentrate the discussion of this paper on the spatio-kinematical distribution of the water masers towards VLA 3-N. The water maser emission toward the region VLA 3-N shows two bow shock-like structures, Northern and Southern, separated from each other by ~ 100 mas (~ 330 AU). The spatial distribution and kinematics of the water masers in this cluster have persisted over a time span of seven years. The Northern bow shock has a somewhat irregular morphology, while the Southern one has a remarkably smooth morphology. We measured the proper motions of 33 water maser features, which have an average proper motion velocity of ~ 1.3 mas/yr (~ 20 km/s). The morphology and the proper motions of this cluster of water masers show systematic expanding motions that could imply one or two different centers of star formation activity. We made a detailed model for the Southern structure, proposing two different kinematic models to explain the 3-dimensional spatio-kinematical distribution of the water masers: (1) a static central source driving the two bow-shock structures; (2) two independent driving sources, one of them exciting the Northern bow-shock structure, and the other one, a young runaway star moving in the local molecular medium exciting and molding the remarkably smoother Southern bow-shock structure. Future observations will be necessary to discriminate between the two scenarios, in particular by identifying the still unseen driving source(s).

We present new [SII] images of the HH 30 jet and counterjet observed in 2006, 2007, and 2010 that allowed us to measure with improved accuracy the positions and proper motions of the jet and counterjet knots. Our results show that the motion of the knots is essentially ballistic, with the exception of the farthest knots, which trace the large scale C-shape bending of the jet. The observed bending of the jet can be produced by a relative motion of the HH 30 star with respect to its surrounding environment, caused either by a possible proper motion of the HH 30 star, or by the entrainment of environment gas by the red lobe of the nearby L1551-IRS 5 outflow. Alternatively, the bending can be produced by the stellar wind from a nearby CTTS, identified in the 2MASS catalog. The proper motion velocities of the knots of the counterjet show more variations than those of the jet. In particular, we identify two knots of the counterjet that have the same kinematic age but whose velocities differ by almost a factor of two. Thus, it appears that counterjet knots launched simultaneously can be ejected with very different velocities. We confirm that the observed wiggling of the jet and counterjet arises from the orbital motion of the jet source in a binary system. Precession is of secondary importance in shaping the jet. We derive an orbital period $\tau_o=114\pm2$ yr and a mass function $m\mu_c^3=0.014\pm0.006$ $M_\odot$. For a mass of the system of $m=0.45\pm0.04$ $M_\odot$ (the value inferred from the disk kinematics) we obtain a mass $m_j=0.31\pm0.04$ $M_\odot$ for the jet source, a mass $m_c=0.14\pm0.03$ $M_\odot$ for the companion, and a binary separation of $a=18.0\pm0.6$ AU. This binary separation coincides with the value required to account for the size of the inner hole observed in the disk, attributed to tidal truncation in a binary system.

HH 223 is the optical counterpart of a larger scale H2 outflow, driven by the protostellar source VLA 2A, in L723. Its poorly collimated and rather chaotic morphology suggested the Integral Field Spectroscopy (IFS) as an appropriate option to map the emission for deriving the physical conditions and the kinematics. Here we present new results based on the IFS observations made with the INTEGRAL system at the WHT. The brightest knots of HH 223 (\sim16 arcsec, 0.02 pc at a distance of 300 pc) were mapped with a single pointing in the spectral range 6200-7700 A. We obtained the emission-line intensity maps for Halpha, [NII] 6584 A and [SII] 6716, 6731 A, and explored the distribution of the excitation and electron density from [NII]/Halpha, [SII]/Halpha, and [SII] 6716/6731 line-ratio maps. Maps of the radial velocity field were obtained. We analysed the 3D-kinematics by combining the knot radial velocities, derived from IFS data, with the knot proper motions derived from multi-epoch, narrow-band images. The intensity maps built from IFS data reproduced well the morphology found in the narrow-band images. We checked the results obtained from previous long-slit observations with those derived from IFS spectra extracted with a similar spatial sampling. At the positions intersected by the slit, the physical conditions and kinematics derived from IFS are compatible with those derived from long-slit data. In contrast, significant discrepancies were found when the results from long-slit data were compared with the ones derived from IFS spectra extracted at positions shifted a few arcsec from those intersected by the slit. This clearly revealed IFS observations as the best choice to get a reliable picture of the HH emission properties.

Aims. We analyze the LBN 978 HII region in order to study the HII/molecular cloud interaction. Methods. We used the IAC-80 telescope to image the region with narrow-band filters and the Very Large Array to obtain a radio continuum map at 3.6 cm. We also used the DSS2 red images and the NRAO VLA Sky Survey at 20 cm. Results. We have discovered a new radial system of dark globules associated with the LBN 978 HII region, containing a group of at least eight cometary bright-rimmed globules with the same morphological type. The brightest source is also detected in the radio continuum. Analysis of optical and radio emission suggest that this object is photoionized by HD 47432, the central star of the LBN 978 HII region.