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.

Recent high-angular resolution (40 mas) ALMA observations at 1.14 mm resolve a compact (R~200 au) flattened dust structure perpendicular to the HH 80-81 jet emanating from the GGD 27-MM1 high-mass protostar, making it a robust candidate for a true accretion disk. The jet/disk system (HH 80-81 / GGD 27-MM1) resemble those found in association with low- and intermediate-mass protostars. We present radiative transfer models that fit the 1.14 mm ALMA dust image of this disk which allow us to obtain its physical parameters and predict its density and temperature structure. Our results indicate that this accretion disk is compact (Rdisk~170 au) and massive (5Msun), about 20% of the stellar mass of 20 Msun. We estimate the total dynamical mass of the star-disk system from the molecular line emission finding a range between 21 and 30 Msun, which is consistent with our model. We fit the density and temperature structures found by our model with power law functions. These results suggest that accretion disks around massive stars are more massive and hotter than their low-mass siblings, but they still are quite stable. We also compare the temperature distribution in the GGD 27-MM1 disk with that found in low- and intermediate-mass stars and discuss possible implications on the water snow line. We have also carried about a study of the distance based on Gaia DR2 data and the population of young stellar objects (YSOs) in this region, and from the extinction maps. We conclude that the source distance is in within 1.2 and 1.4 kpc, closer than what was derived in previous studies (1.7 kpc).

The compact, non-thermal emission in DoAr21 has been studied with the VLBA to investigate the possibility that the residuals of the astrometry fitting are due to the reflex motion induced by a possible companion. We find that the fitting of VLBA astrometric observations of DoAr21 improves significantly by adding the orbital motions of three companions. We obtain an improved distance to the source of $134.6\pm1.0$ pc, and estimate that the central star, DoAr21, has a mass of about $2.04\pm0.70~M_{\odot}$. We suggest that DoAr21 represents a unique case where two sub-stellar companions, DoAr21b and DoAr21c ($m_{b}\sim 35.6\pm27.2~M_{jup}$ and $m_{c}\sim44.0\pm13.6~M_{jup}$, respectively), have been found associated to a relatively low mass, pre-main sequence star. In addition, we find that this WTTau star is an astrometric double system, having a low-mass star companion, DoAr21$B$ ($m_{B}\sim0.35\pm0.12~M_{\odot}$), in a relatively eccentric orbit. The orbit of this low-mass stellar companion is compact, while the Brown Dwarfs are located in external orbits. DoAr21$c$ has the strongest astrometric signature in the periodogram, while DoAr21$B$ has a weak but significant signature. On the other hand, the astrometric signature of DoAr21$b$ does not appear in the periodogram, however, this Brown Dwarf was directly detected in some of the VLBA observations. The estimated orbital periods of DoAr21$B$, DoAr21$b$ and DoAr21$c$ are $P_{B}\sim92.92\pm0.02$, $P_{b}\sim450.9\pm3.8$ and $P_{c}\sim1013.5\pm25.3$ days, respectively. Since the estimated age of this young star is about 0.4$-$0.8 Myrs, the detected Brown Dwarf companion is among the youngest companions observed to date.

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.

Studying the interaction between AGN jets and lobes and their surrounding environment is important in order to understand how they transfer energy to their environment as well as determining the intrinsic physical properties of the sources themselves. This paper presents broadband VLA polarization and Faraday rotation observations of the radio galaxy Coma A (3C 277.3) from 1 to 4 GHz, including archival VLA observations at 4.9 and 15 GHz. Through broadband polarization model-fitting, we find that an external Faraday screen with a turbulent magnetic field provides an appropriate description to most of the data. By combining the polarization and Faraday rotation results with previous H$\alpha$ observations, we identified the H$\alpha$-emitting gas as the Faraday screen responsible for the observed Faraday depolarization. We were able to derive the magnetic field strength in the H$\alpha$-emitting gas, finding typical field strengths of $\sim1$ $\mu$G, which is consistent with studies of the intra-group medium local to other radio galaxies. However, we find a highly depolarized region of the southern lobe coincident with a H$\alpha$ filament that has a field strength comparable to the equipartition field strength in the radio lobe (i.e. $\gtrsim$36 $\mu$G). This implies that the H$\alpha$ filament is internal to the radio emitting plasma. Such clear examples of internal Faraday depolarization are rare, thus providing another key insight into the evolution of radio galaxies and their ability to provide significant feedback on the local gas that would otherwise cool and form stars.

We present new high angular resolution interferometer observations of the v=0 J=14-13 and 15-14 SiS lines towards IRC+10216, carried out with CARMA and ALMA. The maps, with angular resolutions of ~0.25"and 0.55", reveal (1) an extended, roughly uniform, and weak emission with a size of ~0.5", (2) a component elongated approximately along the East-West direction peaking at ~0.13" and 0.17" at both sides of the central star, and (3) two blue- and red-shifted compact components peaking around 0.07" to the NW of the star. We have modeled the emission with a 3D radiation transfer code finding that the observations cannot be explained only by thermal emission. Several maser clumps and one arc-shaped maser feature arranged from 5 to 20R* from the central star, in addition to a thin shell-like maser structure at ~13R* are required to explain the observations. This maser emitting set of structures accounts for 75% of the total emission while the other 25% is produced by thermally excited molecules. About 60% of the maser emission comes from the extended emission and the rest from the set of clumps and the arc. The analysis of a time monitoring of these and other SiS and 29SiS lines carried out with the IRAM 30m telescope from 2015 to present suggests that the intensity of some spectral components of the maser emission strongly depends on the stellar pulsation while other components show a mild variability. This monitoring evidences a significant phase lag of ~0.2 between the maser and NIR light-curves.

Here we present deep (16 mumJy), very high (40 mas) angular resolution 1.14 mm, polarimetric, Atacama Large Millimeter/submillimeter Array (ALMA) observations towards the massive protostar driving the HH 80-81 radio jet. The observations clearly resolve the disk oriented perpendicular to the radio jet, with a radius of ~0.171 arcsec (~291 au at 1.7 kpc distance). The continuum brightness temperature, the intensity profile, and the polarization properties clearly indicate that the disk is optically thick for a radius of R<170 au. The linear polarization of the dust emission is detected almost all along the disk and its properties suggest that dust polarization is produced mainly by self-scattering. However, the polarization pattern presents a clear differentiation between the inner (optically thick) part of the disk and the outer (optically thin) region of the disk, with a sharp transition that occurs at a radius of 0.1 arcsec (~170 au). The polarization characteristics of the inner disk suggest that dust settling has not occurred yet with a maximum dust grain size between 50 and 500 mum. The outer part of the disk has a clear azimuthal pattern but with a significantly higher polarization fraction compared to the inner disk. This pattern is broadly consistent with self-scattering of a radiation field that is beamed radially outward, as expected in the optically thin outer region, although contribution from non-spherical grains aligned with respect to the radiative flux cannot be excluded.

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.

We present subarcsecond angular resolution observations carried out with the Submillimeter Array (SMA) at 880 $\mu$m centered at the B0-type protostar GGD27~MM1, the driving source of the parsec scale HH 80-81 jet. We constrain its polarized continuum emission to be $\lesssim0.8\%$ at this wavelength. Its submm spectrum is dominated by sulfur-bearing species tracing a rotating disk--like structure (SO and SO$_2$ isotopologues mainly), but also shows HCN-bearing and CH$_3$OH lines, which trace the disk and the outflow cavity walls excavated by the HH 80-81 jet. The presence of many sulfurated lines could indicate the presence of shocked gas at the disk's centrifugal barrier or that MM1 is a hot core at an evolved stage. The resolved SO$_2$ emission traces very well the disk kinematics and we fit the SMA observations using a thin-disk Keplerian model, which gives the inclination (47$^{\circ}$), the inner ($\lesssim170$ AU) and outer ($\sim950-1300$~AU) radii and the disk's rotation velocity (3.4 km s$^{-1}$ at a putative radius of 1700 AU). We roughly estimate a protostellar dynamical mass of 4-18\msun. MM2 and WMC cores show, comparatively, an almost empty spectra suggesting that they are associated with extended emission detected in previous low-angular resolution observations, and therefore indicating youth (MM2) or the presence of a less massive object (WMC).

We present the first resolved observations of the 1.3mm polarized emission from the disk-like structure surrounding the high-mass protostar Cepheus A HW2. These CARMA data partially resolve the dust polarization, suggesting an uniform morphology of polarization vectors with an average position angle of 57 degrees and an average polarization fraction of 2.0%. The distribution of the polarization vectors can be attributed to (1) the direct emission of magnetically aligned grains of dust by a uniform magnetic field, or (2) the pattern produced by the scattering of an inclined disk. We show that both models can explain the observations, and perhaps a combination of the two mechanisms produce the polarized emission. A third model including a toroidal magnetic field does not match the observations. Assuming scattering is the polarization mechanism, these observations suggest that during the first few 10000 years of high-mass star formation, grain sizes can grow from 1 to several 10s micron.

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 present high angular resolution ($\sim$0.94$"$) $^{12}$CO(1-0) Atacama Large Millimeter/Submillimeter Array (ALMA) observations obtained during the 2014 long baseline campaign from the eruptive bipolar flow from the multiple XZ Tau stellar system discovered by the Hubble Space Telescope (HST). These observations reveal, for the first time, the kinematics of the molecular flow. The kinematics of the different ejections close to XZ Tau reveal a rotating and expanding structure with a southeast-northwest velocity gradient. The youngest eruptive bubbles unveiled in the optical HST images are inside of this molecular expanding structure. Additionally, we report a very compact and collimated bipolar outflow emanating from XZ Tau A, which indicates that the eruptive outflow is indeed originating from this object. The mass (3 $\times$ 10$^{-7}$ $M_\odot$) and energetics (E$_{kin}$ $=$ 3 $\times$ 10$^{37}$ ergs) for the collimated outflow are comparable with those found in molecular outflows associated with young brown dwarfs.

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.

Three optical telescopes located at the San Pedro M\'artir National Observatory were used for the first time to obtain multi-filter defocused photometry of the transiting extrasolar planets WASP-39b and WASP-43b. We observed WASP-39b with the 2.12m telescope in the U filter for the first time, and additional observations were carried out in the R and I filters using the 0.84m telescope. WASP-43b was observed in VRI with the same instrument, and in the i filter with the robotic 1.50m telescope. We reduced the data using different pipelines and performed aperture photometry with the help of custom routines, in order to obtain the light curves. The fit of the light curves (1.5--2.5mmag rms), and of the period analysis, allowed a revision of the orbital and physical parameters, revealing for WASP-39b a period ($4.0552947 \pm 9.65 \times 10^{-7}$ days) which is $3.084 \pm 0.774$ seconds larger than previously reported. Moreover, we find for WASP-43b a planet/star radius ($0.1738 \pm 0.0033$) which is $0.01637 \pm 0.00371$ larger in the i filter with respect to previous works, and that should be confirmed with additional observations. Finally, we confirm no evidence of constant period variations in WASP-43b.

We present low spectral resolution molecular interferometric observations at 1.2 mm obtained with the Combined Array for Research in Millimetre-wave Astronomy (CARMA) towards the C-rich AGB star IRC+10216. We have mapped the emission of several lines of SiS, H13CN, SiO, and SiC2 in the ground and first excited vibrational states with a high angular resolution of 0.25 arcsec. These observations have allowed us to partially resolve the emission of the envelope at distances from the star <50 stellar radii (R*), where the stellar wind is mainly accelerated. The structure of the molecular emission has been modelled with a 3D radiation transfer code. The emission of line SiS(v=0,J=14-13) is best reproduced with a set of maser emitting arcs arranged between 5 and 20 R*. The abundance of H13CN with respect to H2 decreases from 8e-7 at 1-5 R* to 3e-7 at 20 R*. The SiO observations are explained with an abundance <2e-8 in the shell-like region between 1 and 5 R*. At this point, the SiO abundance sharply increases up to (2-3)e-7. The vibrational temperature of SiO increases by a factor of 2 due North-East between 20 and 50 R*. SiC2 is formed at the stellar surface with an abundance of 8e-7 decreasing down to 8e-8 at 20 R* probably due to depletion on to dust grains. Several asymmetries are found in the abundance distributions of H13CN, SiO, and SiC2 which define three remarkable directions (North-East, South-Southwest, and South-East) in the explored region of the envelope. There are some differences between the red- and blue-shifted emissions of these molecules suggesting the existence of additional asymmetries in their abundance distributions along the line-of-sight.

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.

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 present Combined Array for Research in Millimeter-wave Astronomy (CARMA) 3.5 mm observations and SubMillimeter Array (SMA) 870 microns observations toward the high-mass star-forming region IRAS 1816-2048, the core of the HH 80/81/80N system. Molecular emission from HCN, HCO+ and SiO is tracing two molecular outflows (the so-called Northeast and Northwest outflows). These outflows have their origin in a region close to the position of MM2, a millimeter source known to harbor a couple protostars. We estimate for the first time the physical characteristics of these molecular outflows, which are similar to those of 1000-5000 Lo protostars, suggesting that MM2 harbors high-mass protostars. High-angular resolution CO observations show an additional outflow due southeast. We identify for the first time its driving source, MM2(E), and see evidence of precession. All three outflows have a monopolar appearance, but we link the NW and SE lobes, explaining their asymmetric shape as a consequence of possible deflection.

We present (sub)millimeter line and continuum observations in a mosaicing mode of the massive star forming region Cepheus A East made with the Submillimeter Array (SMA). Our mosaic covers a total area of about 3$'$ $\times$ 12$'$ centered in the HW 2/3 region. For the first time, this observational study encloses a high angular resolution ($\sim$ 3$''$) together with a large scale mapping of Cepheus A East. We report compact and high velocity $^{12}$CO(2-1) emission associated with the multiple east-west bright H$_2$ condensations present in the region. Blueshifted and redshifted gas emission is found towards the east as well as west of HW 2/3. The observations suggest the presence of multiple large-scale east-west outflows that seems to be powered at smaller scales by radio sources associated with the young stars HW2, HW3c and HW3d. A kinematical study of part of the data suggests that the molecular outflow powered by HW2 is precesing with time as recently reported. Our data reveal five periodic ejections of material separated approximately every 10$^\circ$ as projected in the plane of the sky. The most recent ejections appear to move toward the plane of the sky. An energetic explosive event as the one that occurred in Orion BN/KL or DR21 does not explain the kinematics, and the dynamical times of the multiple ejections found here. The continuum observations only revealed a strong millimeter source associated with the HW 2/3 region. High angular resolution observations allow us to resolve this extended dusty object in only two compact sources (with spatial sizes of approximately 300 AU) associated with HW2 and HW3c. Finally, the bright optical/X-Ray HH 168 -- GDD37 object might be produced by strong shocks related with the outflow from HW3c.

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

In this {\it Letter}, we present sensitive millimeter SiO (J=5-4; $\nu$=0) line observations of the outflow arising from the enigmatic object Orion Source I made with the Atacama Large Millimeter/Submillimeter Array (ALMA). The observations reveal that at scales of a few thousand AU, the outflow has a marked "butterfly" morphology along a northeast-southwest axis. However, contrary to what is found in the SiO and H$_2$O maser observations at scales of tens of AU, the blueshifted radial velocities of the moving gas are found to the northwest, while the redshifted velocities are in the southeast. The ALMA observations are complemented with SiO (J=8-7; $\nu$=0) maps (with a similar spatial resolution) obtained with the Submillimeter Array (SMA). These observations also show a similar morphology and velocity structure in this outflow. We discuss some possibilities to explain these differences at small and large scales across the flow.

VLBI multi-epoch water maser observations are a powerful tool to study the dense, warm shocked gas very close to massive protostars. The very high-angular resolution of these observations allow us to measure the proper motions of the masers in a few weeks, and together with the radial velocity, to determine their full kinematics. In this paper we present a summary of the main observational results obtained toward the massive star-forming regions of Cepheus A and W75N, among them: (i) the identification of different centers of high-mass star formation activity at scales of 100 AU; (ii) the discovery of new phenomena associated with the early stages of high-mass protostellar evolution (e.g., isotropic gas ejections); and (iii) the identification of the simultaneous presence of a wide-angle outflow and a highly collimated jet in the massive object Cep A HW2, similar to what is observed in some low-mass protostars. Some of the implications of these results in the study of high-mass star formation are discussed.

VLBI multi-epoch water maser observations are a powerful tool to study the gas very close to the central engine responsible for the phenomena associated with the early evolution of massive protostars. In this paper we present a summary of the main observational results obtained toward the massive star-forming regions of Cepheus A and W75N. These observations revealed unexpected phenomena in the earliest stages of evolution of massive objects (e.g., non-collimated "short-lived" pulsed ejections in different massive protostars), and provided new insights in the study of the dynamic scenario of the formation of high-mass stars (e.g., simul- taneous presence of a jet and wide-angle outflow in the massive object Cep A HW2, similar to what is observed in low-mass protostars). In addition, with these observations it has been possible to identify new, previously unseen centers of high-mass star formation through outflow activity.

We present several molecular line emission arcsec and subarcsec observations obtained with the Submillimeter Array (SMA) in the direction of the massive protostar IRAS 18162-2048, the exciting source of HH 80-81. The data clearly indicates the presence of a compact (radius~425-850 AU) SO2 structure, enveloping the more compact (radius~150 AU) 1.4 millimeter dust emission (reported in a previous paper). The emission spatially coincides with the position of the prominent thermal radio jet which terminates at the HH 80-81 and HH 80N Herbig-Haro objects. Furthermore, the molecular emission is elongated in the direction perpendicular to the axis of the thermal radio jet, suggesting a disk-like structure. We derive a total dynamic mass (disk-like structure and protostar) of 11-15 msun. The SO2 spectral line data also allow us to constrain the structure temperature between 120-160 K and the volume density > 2x10^9 cm-3. We also find that such a rotating flattened system could be unstable due to gravitational disturbances. The data from C17O line emission show a dense core within this star-forming region. Additionally, the H2CO and the SO emissions appear clumpy and trace the disk-like structure, a possible interaction between a molecular core and the outflows, and in part, the cavity walls excavated by the thermal radio jet.

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

The GGD27 complex includes the HH 80-81-80N system, which is one of the most powerful molecular outflows associated with a high mass star-forming region observed up to now. This outflow is powered by the star associated with the source IRAS 18162-2048. Here we report the detection of continuum emission at sub-arcsec/arcsec resolution with the Submillimeter Array at 1.36mm and 456microns, respectively. We detected dust emission arising from two compact cores, MM1 and MM2, separated by about 7" (~12000AU in projected distance). MM1 spatially coincides with the powerful thermal radio continuum jet that powers the very extended molecular outflow, while MM2 is associated with the protostar that drives the compact molecular outflow recently found in this region. High angular resolution obervations at 1.36mm show that MM1 is unresolved and that MM2 splits into two subcomponents separated by ~1". The mass of MM1 is about 4Msun and it has a size of <300AU. This is consistent with MM1 being associated with a massive and dense (n(H2)>10^9cm-3) circumstellar dusty disk surrounding a high-mass protostar, which has not developed yet a compact HII region. On the other hand, the masses of the two separate components of MM2 are about 2Msun each. One of these components is a compact core with an intermediate-mass young protostar inside and the other component is probably a pre-stellar core. MM1 is the brigthest source at 1.36mm, while MM2 dominates the emission at 456microns. These are the only (sub)millimeter sources detected in the SMA observations. Hence, it seems that both sources may contribute significantly to the bolometric luminosity of the region. Finally, we argue that the characteristics of these two sources indicate that MM2 is probably in an earlier evolutionary stage than MM1.

We present five epochs of VLBI water maser observations around the massive protostar Cepheus A HW2 with 0.4 mas (0.3 AU) resolution. The main goal of these observations was to follow the evolution of the remarkable water maser linear/arcuate structures found in earlier VLBI observations. Comparing the data of our new epochs of observation with those observed five years before, we find that at "large" scales of > 1" (700 AU) the main regions of maser emission persist, implying that both the surrounding medium and the exciting sources of the masers have been relatively stable during that time span. However, at smaller scales of < 0.1" (70 AU) we see large changes in the maser structures, particularly in the expanding arcuate structures R4 and R5. R4 traces a nearly elliptical patchy ring of ~ 70 mas size (50 AU) with expanding motions of ~ 5 mas/yr (15 km/s). This structure is probably driven by the wind of a still unidentified YSO located at the centre of the ring (~ 0.18" south of HW2). On the other hand, the R5 expanding bubble structure (driven by the wind of a previously identified YSO located ~ 0.6" south of HW2) is currently dissipating in the circumstellar medium and losing its previous degree of symmetry, indicating a very short-lived event. In addition, our results reveal, at scales of ~ 1" (700 AU), the simultaneous presence of a relatively slow (~ 10-70 km/s) wide-angle outflow (opening angle of ~ 102 deg, traced by the masers, and the fast (~ 500~km/s) highly collimated radio jet associated with HW2 (opening angle of ~ 18 deg, previously observed with the VLA. This simultaneous presence of a wide-angle outflow and a highly collimated jet associated with a massive protostar is similar to what is found in some low-mass YSOs. The implications of these results in the study of the formation of high-mass stars are discussed.

We have done a statistical analysis of Very Long Baseline Array (VLBA) data of water masers in the star-forming regions (SFRs) Cepheus A and W75 N, using correlation functions to study the spatial clustering and Doppler-velocity distribution of these masers. Two-point spatial correlation functions show a characteristic scale size for clusters of water maser spots < or ~1 AU, similar to the values found in other SFRs. This suggests that the scale for water maser excitation tends to be < or ~1 AU. Velocity correlation functions show power-law dependences with indices that can be explained by regular velocity fields, such as expansion and/or rotation. These velocity fields are similar to those indicated by the water maser proper-motion measurements; therefore, the velocity correlation functions appear to reveal the organized motion of water maser spots on scales larger than 1 AU.

Context. Mathematical optimization can be used as a computational tool to obtain the optimal solution to a given problem in a systematic and efficient way. For example, in twice-differentiable functions and problems with no constraints, the optimization consists of finding the points where the gradient of the objective function is zero and using the Hessian matrix to classify the type of each point. Sometimes, however it is impossible to compute these derivatives and other type of techniques must be employed such as the steepest descent/ascent method and more sophisticated methods such as those based on the evolutionary algorithms. Aims. We present a simple algorithm based on the idea of genetic algorithms (GA) for optimization. We refer to this algorithm as AGA (Asexual Genetic Algorithm) and apply it to two kinds of problems: the maximization of a function where classical methods fail and model fitting in astronomy. For the latter case, we minimize the chi-square function to estimate the parameters in two examples: the orbits of exoplanets by taking a set of radial velocity data, and the spectral energy distribution (SED) observed towards a YSO (Young Stellar Object). Methods. The algorithm AGA may also be called genetic, although it differs from standard genetic algorithms in two main aspects: a) the initial population is not encoded, and b) the new generations are constructed by asexual reproduction. Results. Applying our algorithm in optimizing some complicated functions, we find the global maxima within a few iterations. For model fitting to the orbits of exoplanets and the SED of a YSO, we estimate the parameters and their associated errors.

We present new 7 mm and archive 1.3 cm high angular resolution observations of the HL/XZ Tau region made with the VLA. At 7 mm, the emission from HL Tau seems to be arising in a clumpy disk with radius of order 25 AU. The 1.3 cm emission from XZ Tau shows the emission from a binary system with 0"3 (42 AU) separation, known from previous optical/IR observations. However, at 7 mm, the southern radio component resolves into a binary with 0"09 (13 AU) separation, suggesting that XZ Tau is actually a triple star system. We suggest that the remarkable ejection of gas from the XZ Tau system observed with the HST may be related to a periastron passage of this newly discovered close binary system.

S140 IRS1 is a remarkable source where the radio source at the center of the main bipolar molecular outflow in the region is elongated perpendicular to the axis of the outflow, an orientation opposite to that expected if the radio source is a thermal jet exciting the outflow. We present results of 1.3 cm continuum and H2O maser emission observations made with the VLA in its A configuration toward this region. In addition, we also present results of continuum observations at 7 mm and re-analyse observations at 2, 3.5 and 6 cm (previously published). IRS 1A is detected at all wavelengths, showing an elongated structure. Three water maser spots are detected along the major axis of the radio source IRS 1A. We have also detected a new continuum source at 3.5 cm (IRS 1C) located ~0.6'' northeast of IRS 1A. The presence of these two YSOs (IRS 1A and 1C) could explain the existence of the two bipolar molecular outflows observed in the region. In addition, we have also detected three continuum clumps (IRS 1B, 1D and 1E) located along the major axis of IRS 1A. We discuss two possible models to explain the nature of IRS 1A: a thermal jet and an equatorial wind.

We report SMA 335 GHz continuum observations with angular resolution of ~0.''3, together with VLA ammonia observations with ~1'' resolution toward Cep A HW 2. We find that the flattened disk structure of the dust emission observed by Patel et al. is preserved at the 0.''3 scale, showing an elongated structure of ~$0.''6 size (450 AU) peaking on HW 2. In addition, two ammonia cores are observed, one associated with a hot-core previously reported, and an elongated core with a double peak separated by ~1.''3 and with signs of heating at the inner edges of the gas facing HW 2. The double-peaked ammonia structure, as well as the double-peaked CH3CN structure reported previously (and proposed to be two independent hot-cores), surround both the dust emission as well as the double-peaked SO2 disk structure found by Jimenez-Serra et al. All these results argue against the interpretation of the elongated dust-gas structure as due to a chance-superposition of different cores; instead, they imply that it is physically related to the central massive object within a disk-protostar-jet system.

The processes leading to the birth of low-mass stars such as our Sun have been well studied, but the formation of high-mass (> 8 x Sun's mass) stars has heretofore remained poorly understood. Recent observational studies suggest that high-mass stars may form in essentially the same way as low-mass stars, namely via an accretion process, instead of via merging of several low-mass (< 8 Msun) stars. However, there is as yet no conclusive evidence. Here, we report the discovery of a flattened disk-like structure observed at submillimeter wavelengths, centered on a massive 15 Msun protostar in the Cepheus-A region. The disk, with a radius of about 330 astronomical units (AU) and a mass of 1 to 8 Msun, is detected in dust continuum as well as in molecular line emission. Its perpendicular orientation to, and spatial coincidence with the central embedded powerful bipolar radio jet, provides the best evidence yet that massive stars form via disk accretion in direct analogy to the formation of low-mass stars.

We report centimeter continuum observations with the VLA and the VLBA as well as mid-infrared observations with COMICS/SUBARU toward the components of the YLW 15 very young binary system, VLA 1 and VLA 2. The centimeter emission of the two components traces partially thick free-free emission, likely due to collimated, ionized winds. VLA 1 is an embedded protostar, undetected in the near-IR, and possibly in the Class 0 to Class I transition and powering a Herbig-Haro outflow. Its mid-IR emission appears slightly resolved with a diameter of ~16 AU, possibly tracing circumstellar material from both the envelope and the disk. VLA 2 is a typical Class I object, unresolved in the mid-IR, and is the responsible of the strong X-ray emission associated with YLW 15. The expected centimeter ''peri-stellar'' emission associated with the X-ray emission is not detected with the VLBA at 6 cm likely due to the high optical depth of the free-free emission. Strikingly, the near to mid-IR properties of YLW 15 suggest that VLA 1 is a more embedded YSO, or alternatively, less luminous than VLA 2, whereas orbital proper motions of this binary system indicate that VLA 1 is more massive than VLA 2. This result is apparently against the expected evolutionary scenario, where one expects that the more massive YSO in a binary system is the more evolved and luminous YSO. Finally, the nearby source YLW 16A is detected with the VLA, its position coincides with reported near-IR and X-ray sources.

Spatial and spectral profiles of O VI emission behind a shock wave on the northern edge of the Cygnus Loop were obtained with the FUSE satellite. The velocity width of the narrowest O VI profile places a tight constraint on the electron-ion and ion-ion thermal equilibration in this 350 $\rm km s^{-1}$ collisionless shock. Unlike faster shocks in SN1006 and in the heliosphere, this shock brings oxygen ions and protons to within a factor of 2.5 of the same temperature. Comparison with other shocks suggests that shock speed, rather than Alfv\'{e}n Mach number, may control the degree of thermal equilibration. We combine the O VI observations with a low resolution far UV spectrum from HUT, an H$\alpha$ image and ROSAT PSPC X-ray data to constrain the pre-shock density and the structure along the line-of-sight. As part of this effort, we model the effects of resonance scattering of O VI photons within the shocked gas and compute time-dependent ionization models of the X-ray emissivity. Resonance scattering affects the O VI intensities at the factor of 2 level, and the soft spectrum of the X-ray rim can be mostly attributed to departures from ionization equilibrium. The pre-shock density is about twice the canonical value for the Cygnus Loop X-ray emitting shocks.

We present multi-frequency VLA continuum observations towards 8 star forming regions with molecular and optical outflows: L1489, HH 68-69, HH 94-95, NGC 2264D, L1681B, L778, MWC 1080 and V645 Cyg. We detect three thermal radio jets, L1489, YLW 16A in L1681B and NGC 2264D VLA 7, associated with molecular and/or HH outflows. The L1489 and NGC 2264D VLA 7 thermal radio jets appear elongated in the direction of the larger scale outflow. We report the first tentative detection of a non-thermal radio jet, L778 VLA 5, associated with a low mass Class I protostar and powering a molecular outflow. For HH68-69, HH 94-95 and the molecular outflow in NGC 2264D we could not identify a candidate of the exciting source of these outflows. The radio emission associated with V645 Cyg is quite extended, ~0.1 pc, and time variable. We detect three radio sources in the MWC 1080 that could be associated with YSOs.

We present a long-slit spectrum of the Herbig-Haro object HH 111, obtained with SITS on the Hubble Space Telescope. This spectrum has a spectral resolution of approx. 50 km/s, so it gives a good picture of the kinematical properties of the observed object at very high angular resolution. We find that some of the knots along the jet are associated with sudden drops in the radial velocity (modulus), confirming that the emission from the knots is formed in shocks. We interprt the observed position-velocity diagrams in terms of a model of a jet from a variable source, and we attempt to carry out a reconstruction of the ejection variability necessary for reproducing the observed kinematical structure.

We present images and spectra of the Cepheus E (Cep E) region at both optical and infrared wavelengths. Only the brightest region of the southern lobe of the Cep E outflow reveals optical emission, suggesting that the extinction close to the outflow source plays an important r\^ole in the observed difference between the optical and IR morphologies. Cep E is a unique object since it provides a link between the spectroscopic properties of the optical Herbig-Haro (HH) objects and those of deeply embedded outflows.