We present sensitive and high angular resolution ($\sim$0.2-0.3$''$) (sub)millimeter (230 and 345 GHz) continuum and CO(2$-$1)/CO(3$-$2) line archive observations of the disk star system in UX Tauri carried out with ALMA (The Atacama Large Millimeter/Submillimeter Array). These observations reveal the gas and dusty disk surrounding the young star UX Tauri A with a large signal-to-noise ratio ($>$400 in the continuum and $>$50 in the line), and for the first time is detected the molecular gas emission associated with the disk of UX Tauri C (with a size for the disk of $<$56 au). No (sub)millimeter continuum emission is detected at 5$\sigma$-level (0.2 mJy at 0.85 mm) associated with UX Tauri C. For the component UX Tauri C, we estimate a dust disk mass of $\leq$ 0.05 M$_\oplus$. Additionally, we report a strong tidal disk interaction between both disks UX Tauri A/C, separated 360 au in projected distance. The CO line observations reveal marked spiral arms in the disk of UX Tauri A and an extended redshifted stream of gas associated with the UX Tauri C disk. No spiral arms are observed in the dust continuum emission of UX Tauri A. Assuming a Keplerian rotation we estimate the enclosed masses (disk$+$star) from their radial velocities in 1.4 $\pm$ 0.6 M$_\odot$ for UX Tauri A, and 70 $\pm$ 30 / $\sin i$ Jupiter masses for UX Tauri C (the latter coincides with the mass upper limit value for a brown dwarf). The observational evidence presented here lead us to propose that UX Tauri C is having a close approach of a possible wide, evolving and eccentric orbit around the disk of UX Tauri A causing the formation of spiral arms and the stream of molecular gas falling towards UX Tauri C.

We present Atacama Large Millimeter/submillimeter Array (ALMA) and Karl G. Jansky Very Large Array (JVLA) observations of the massive infrared dark cloud NGC 6334S (also known as IRDC G350.56+0.44), located at the southwestern end of the NGC 6334 molecular cloud complex. The H$^{13}$CO$^{+}$ and the NH$_{2}$D lines covered by the ALMA observations at a $\sim$3$^{\prime\prime}$ angular resolution ($\sim$0.02 pc) reveal that the spatially unresolved non-thermal motions are predominantly subsonic and transonic, a condition analogous to that found in low-mass star-forming molecular clouds. The observed supersonic non-thermal velocity dispersions in massive star forming regions, often reported in the literature, might be significantly biased by poor spatial resolutions that broaden the observed line widths due to unresolved motions within the telescope beam. Our 3~mm continuum image resolves 49 dense cores, whose masses range from 0.17 to 14 $M_{\odot}$. The majority of them are resolved with multiple velocity components. Our analyses of these gas velocity components find an anti-correlation between the gas mass and the virial parameter. This implies that the more massive structures tend to be more gravitationally unstable. Finally, we find that the external pressure in the NGC 6334S cloud is important in confining these dense structures, and may play a role in the formation of dense cores, and subsequently, the embedded young stars.

We study the evolution of the virial parameter, $\alpha$, and the star formation activity of a star forming region in a numerical simulation and two observational samples likely to be at different evolutionary stages, the Pipe and the G14.225 clouds. We consider a numerical simulation of turbulence in the warm atomic gas, in which clouds form by a combination of thermal instability and turbulent compressions. In this simulation we select a region with physical properties similar to those of the Pipe in both its global physical properties and the dense cores it contains, and show that this region evolves to become similar to the G14.225 cloud and its substructure in the course of $\sim$1.5 Myr. Within this region, we follow the evolution of $\alpha$ and the star formation activity of the cores over this time interval, during which they grow in mass by accretion and develop further substructure of higher density. We find that the individual cores evolve by exhibiting first a mild decrease in $\alpha$ followed by a rapid increase when star formation begins. We suggest this evolution is due to an early loss of external compressive kinetic energy followed by an increase of the gravitationally driven motions. Nevertheless, collectively, the ensamble of clumps and cores reproduces the recently observed trend of decreasing $\alpha$ for objects of higher mass. We also find that the star formation rate and star formation efficiency increase monotonically as the region evolves. Thus, we propose that the energy balance, $\alpha$, and the star formation activity determine the evolutionary stage of a star forming region, and that a star forming region like the Pipe may evolve into a G14.225-like stage over the course of a few megayears if it is allowed to accrete from its environment. This evolution is consistent with the recently proposed scenario of hierarchical gravitational collapse.

We present a unified description of the scenario of Global Hierarchical Collapse and fragmentation (GHC) in molecular clouds (MCs), owing to the continuous decrease of the average Jeans mass in the contracting cloud. GHC constitutes a regime of collapses within collapses, in which small-scale collapses begin at later times, but occur on shorter timescales than large-scale ones. The difference in timescales allows for most of the clouds' mass to be dispersed by feedback from the first massive stars, maintaining the global star formation rate low. All scales accrete from their parent structures. The main features of GHC are: star-forming MCs are in an essentially pressureless regime, which produces filaments that accrete onto clumps and cores ("hubs"). The filaments constitute the collapse flow from cloud to hub scales and may approach a quasi-stationary state; the molecular and dense mass fractions of the clouds increase over time; the first (low-mass) stars appear several Myr after global contraction began; more massive stars appear after a few Myr in massive hubs resulting from the collapse of larger scales; the minimum fragment mass may extend well into the brown-dwarf regime; Bondi-Hoyle-Lyttleton accretion occurs at the protostellar and core scales, accounting for a near-Salpeter IMF; the extreme anisotropy of the filamentary network explains the difficulty in detecting large-scale infall signatures; the balance between inertial and gravitationally-driven motions in clumps evolves during the contraction; prestellar cores adopt Bonnor-Ebert-like profiles, but are contracting ever since early times when they may appear to be unbound and to require pressure confinement; stellar clusters develop radial age and mass segregation gradients. Finally, we discuss the incompatibility between supersonic turbulence and the observed scalings in the molecular hierarchy.

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

We present a comparative study of the physical properties and the spatial distribution of column density peaks in two Giant Molecular Clouds (GMC), the Pipe Nebula and Orion A, which exemplify opposite cases of star cluster formation stages. The density peaks were extracted from dust extinction maps constructed from Herschel/SPIRE farinfrared images. We compare the distribution functions for dust temperature, mass, equivalent radius and mean volume density of peaks in both clouds, and made a more fair comparison by isolating the less active Tail region in Orion A and by convolving the Pipe Nebula map to simulate placing it at a distance similar to that of the Orion Complex. The peak mass distributions for Orion A, the Tail, and the convolved Pipe, have similar ranges, sharing a maximum near 5 M$_\odot$, and a similar power law drop above 10 M$_\odot$. Despite the clearly distinct evolutive stage of the clouds, there are very important similarities in the physical and spatial distribution properties of the column density peaks, pointing to a scenario where they form as a result of uniform fragmentation of filamentary structures across the various scales of the cloud, with density being the parameter leading the fragmentation, and with clustering being a direct result of thermal fragmentation at different spatial scales. Our work strongly supports the idea that the formation of clusters in GMC could be the result of the primordial organization of pre-stellar material

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

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

The explosive outflows are a newly-discovered family of molecular outflows associated with high-mass star forming regions. Such energetic events are possibly powered by the release of gravitational energy related with the formation of a (proto)stellar merger or a close stellar binary. Here, we present sensitive and high angular resolution observations (0.85$''$) archival CO(J=3-2) observations carried out with the Submillimeter Array (SMA) of the high-mass star forming region G5.89$-$0.39 that reveal the possible presence of an explosive outflow. We find six well-defined and narrow straight filament-like ejections pointing back approximately to the center of an expanding molecular and ionized shell located at the center of this region. These high velocity ($-$120 to $+$100 km s$^{-1}$) filaments follow a Hubble-like velocity law with the radial velocities increasing with the projected distance. The estimated kinematical age of the filaments is about of 1000 yrs, a value similar to the dynamical age found for the expanding ionized shell. G5.89 is the thus the third explosive outflow reported in the galaxy (together with Orion BN-KL and DR21) and argues in favor of the idea that this is a frequent phenomenon. In particular, explosive outflows, in conjunction with runaway stars, demonstrate that dynamical interactions in such groups are a very important ingredient in star formation.

Filaments are ubiquitous structures in molecular clouds and play an important role in the mass assembly of stars. We present results of dynamical stability analyses for filaments in the infrared dark cloud G14.225$-$0.506, where a delayed onset of massive star formation was reported in the two hubs at the convergence of multiple filaments of parsec length. Full-synthesis imaging is performed with the Atacama Large Millimeter/submillimeter Array (ALMA) to map the $\mathrm{N_2H^+} \; (1-0)$ emission in two hub-filament systems with a spatial resolution of $\sim 0.034 \; \mathrm{pc}$. Kinematics are derived from sophisticated spectral fitting algorithm that accounts for line blending, large optical depth, and multiple velocity components. We identify five velocity coherent filaments and derive their velocity gradients with principal component analysis. The mass accretion rates along the filaments are up to $10^{-4} \; \mathrm{M_\odot \, \mathrm{yr^{-1}}}$ and are significant enough to affect the hub dynamics within one free-fall time ($\sim 10^5 \; \mathrm{yr}$). The $\mathrm{N_2H^+}$ filaments are in equilibrium with virial parameter $\alpha_\mathrm{vir} \sim 1.2$. We compare $\alpha_\mathrm{vir}$ measured in the $\mathrm{N_2H^+}$ filaments, $\mathrm{NH_3}$ filaments, $870 \; \mu\mathrm{m}$ dense clumps, and $3 \; \mathrm{mm}$ dense cores. The decreasing trend in $\alpha_\mathrm{vir}$ with decreasing spatial scales persists, suggesting an increasingly important role of gravity at small scales. Meanwhile, $\alpha_\mathrm{vir}$ also decreases with decreasing non-thermal motions. In combination with the absence of high-mass protostars and massive cores, our results are consistent with the global hierarchical collapse scenario.

During the last decades, a great interest has emerged to know if even the most massive stars in our galaxy (namely the spectral O-type stars) are formed in a similar manner as the low- and intermediate-mass stars, that is, through the presence of accreting disks and powerful outflows. Here, using sensitive observations of the Atacama Large Millimeter/Submillimeter Array (ALMA), we report a resolved Keplerian disk (with fifteen synthesized beams across its major axis) surrounding the deeply embedded O-type protostar IRAS16547$-$4247. The disk shows some asymmetries that could arise because of the disk is unstable and fragmenting or maybe because of different excitation conditions within the disk. The enclosed mass estimated from the disk Keplerian radial velocities is 25$\pm$3 M$_\odot$. The molecular disk is at the base of an ionized thermal radio jet and is approximately perpendicular to the jet axis orientation. We additionally find the existence of a binary system of compact dusty objects at the center of the accreting disk, which indicates the possible formation of an O-type star and a companion of lower mass. This is not surprising due to the high binary fraction reported in massive stars. Subtracting the contribution of the dusty disk plus the envelope and the companion, we estimated a mass of 20 M$_\odot$ for the central star.

The filamentary infrared dark cloud L1287 is actively forming a dense cluster of low-mass YSOs at its inner $\sim$0.1 pc region. To help understand the origin of this low-mass YSO cluster, the present work aims at resolving the gas structures and kinematics. We have performed $\sim$1$"$ angular resolution ($\sim$930 AU) SMA observations at $\sim$1.3 mm. From a $\sim$2$"$ resolution 1.3 mm continuum image we identified six dense cores, namely SMA1-6 with masses in the range of $\sim0.4-4$ M$_\odot$. From a $\sim$1$"$ resolution 1.3 mm continuum image, we find a high fragmentation level, with 14 compact millimeter sources within 0.1 pc (two of them associated with the known accretion outburst YSOs RNO 1C and RNO 1B). The dense gas tracer DCN (3--2) traces well the dust continuum emission and shows a clear velocity gradient along the NW-SE direction centered at SMA3. There is another velocity gradient with opposite direction around the most luminous YSO IRAS 00338+6312. The fragmentation within 0.1 pc in L1287 is very high compared to other regions at the same spatial scales. The incoherent motions of dense gas flows are sometimes interpreted by being influenced by (proto)stellar feedback (e.g., outflows), which is not yet ruled out in this particular target source. The directions of the velocity gradients traced by DCN are approximately perpendicular to those of the dominant CO outflow(s). Therefore, we alternatively hypothesize that the velocity gradients revealed by DCN trace the convergence from the $\gtrsim$0.1 pc scales infalling motion towards the rotational motions around the more compact ($\sim0.02$ pc) sources. This global molecular gas converging flow may feed the formation of the dense low-mass YSO cluster. IRAS 00338+6312 is the most likely powering source of the dominant CO outflow. A compact blue-shifted outflow from RNO 1C is also identified.

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

Aims: We aim to understand the fragmentation as well as the disk formation, outflow generation and chemical processes during high-mass star formation on spatial scales of individual cores. Methods: Using the IRAM Northern Extended Millimeter Array (NOEMA) in combination with the 30m telescope, we have observed in the IRAM large program CORE the 1.37mm continuum and spectral line emission at high angular resolution (~0.4'') for a sample of 20 well-known high-mass star-forming regions with distances below 5.5kpc and luminosities larger than 10^4Lsun. Results: We present the overall survey scope, the selected sample, the observational setup and the main goals of CORE. Scientifically, we concentrate on the mm continuum emission on scales on the order of 1000AU. We detect strong mm continuum emission from all regions, mostly due to the emission from cold dust. The fragmentation properties of the sample are diverse. We see extremes where some regions are dominated by a single high-mass core whereas others fragment into as many as 20 cores. A minimum-spanning-tree analysis finds fragmentation at scales on the order of the thermal Jeans length or smaller suggesting that turbulent fragmentation is less important than thermal gravitational fragmentation. The diversity of highly fragmented versus singular regions can be explained by varying initial density structures and/or different initial magnetic field strengths. Conclusions: The smallest observed separations between cores are found around the angular resolution limit which indicates that further fragmentation likely takes place on even smaller spatial scales. The CORE project with its numerous spectral line detections will address a diverse set of important physical and chemical questions in the field of high-mass star formation.

The formation of the most massive O-type stars is poorly understood. We present a case study of a young massive clump from the ATLASGAL survey, G328.2551-0.5321. It exhibits a bolometric luminosity of 1.3$\times$10$^4$ L$_{\odot}$ corresponding to a current protostellar mass of $\sim$11 and 16 M$_{\odot}$. We analyze high angular-resolution observations with ALMA at $\sim$0.17" corresponding a physical scale of $\sim$400 au in dust continuum and molecular lines. The dust continuum emission reveals a single high-mass protostellar envelope and shows evidence for a marginally resolved continuum source. We detect a rotational line of CH$_3$OH within its $v_{\rm t}$=1 torsionally excited state revealing two bright peaks of emission spatially offset from the dust continuum peak, and exhibiting a distinct velocity component $\pm$4.5 km s$^{-1}$ offset compared to the source $v_{\rm lsr}$. Local thermodynamic equilibrium analysis suggests N(CH$_3$OH)=1.2$-$2$\times$10$^{19}$ cm$^{-2}$, and kinetic temperatures of 160$-$200 K at the position of these peaks. Their velocity shifts correspond well to the expected Keplerian velocity around a central object with 15M$_{\odot}$ consistent with the mass estimate based on the source's bolometric luminosity. We propose a picture where the CH$_3$OH emission peaks trace the accretion shocks around the centrifugal barrier, pinpointing the interaction region between the collapsing envelope and an accretion disk. Because the HC$_3$N $v_{\rm 7}$=1e ($J$=38-37) line shows compact emission, and a velocity pattern consistent with models of Keplerian rotation, we suggest that this could be a new tracer for compact accretion disks around high-mass protostars. The estimated physical properties of the accretion disk suggest a specific angular momentum several times larger than typically observed towards low-mass protostars.

We present subarcsecond 1.3 mm continuum ALMA observations towards the Orion Molecular Cloud 1 South (OMC-1S) region, down to a spatial resolution of 74 AU, which reveal a total of 31 continuum sources. We also present subarcsecond 7 mm continuum VLA observations of the same region, which allow to further study fragmentation down to a spatial resolution of 40 AU. By applying a Mean Surface Density of Companions method we find a characteristic spatial scale at ~560 AU, and we use this spatial scale to define the boundary of 19 `cores' in OMC-1S as groupings of millimeter sources. We find an additional characteristic spatial scale at ~2900 AU, which is the typical scale of the filaments in OMC-1S, suggesting a two-level fragmentation process. We measured the fragmentation level within each core and find a higher fragmentation towards the southern filament. In addition, the cores of the southern filament are also the densest (within 1100 AU) cores in OMC-1S. This is fully consistent with previous studies of fragmentation at spatial scales one order of magnitude larger, and suggests that fragmentation down to 40 AU seems to be governed by thermal Jeans processes in OMC-1S.

The early evolutionary stage of brown dwarfs are not very well characterized, specially during the embedded phase. To gain insight into the dominant formation mechanism of very low-mass objects and brown dwarfs, we conducted deep observations at 870$\mu$m with the LABOCA bolometer at the APEX telescope. Our goal was to identify young sub-mm sources in the Barnard 30 dark cloud. We complemented these data with multi-wavelength observations from the optical to the far-IR and. As a result, we have identified 34 submm sources and a substantial number of possible and probable Barnard 30 members within each individual APEX/LABOCA beam. They can be classified in three distinct groups. First, 15 out of these 34 have a clear optical or IR counterpart to the submm peak and nine of them are potential proto-BDs candidates. Moreover, a substantial number of them could be multiple systems. A second group of 13 sources comprises candidate members with significant infrared excesses located away from the central submm emission. All of them include brown dwarf candidates, some displaying IR excess, but their association with submm emission is unclear. In addition, we have found six starless cores and, based on the total dust mass estimate, three might be pre-substellar (or pre-BDs) cores. Finally, the complete characterization of our APEX/LABOCA sources, focusing on those detected at 24 and/or 70 $\mu$m, indicates that in our sample of 34 submm sources there are, at least: two WTTs, four CTTs, five young stellar objects (YSOs), eight proto-BD candidates (with another three dubious cases), and one Very Low Luminosity object (VeLLO).

In this work we present ALMA continuum observations at 880 $\mu$m of 30 sub-mm cores previously identified with APEX/LABOCA at 870$\mu$m in the Barnard 30 cloud. The main goal is to characterize the youngest and lowest mass population in the cloud. As a result, we report the detection of five (out of 30) spatially unresolved sources with ALMA, with estimated masses between 0.9 and 67 M$_{\rm Jup}$. From these five sources, only two show gas emission. The analysis of multi-wavelength photometry from these two objects, namely B30-LB14 and B30-LB19, is consistent with one Class II- and one Class I low-mass stellar object, respectively. The gas emission is consistent with a rotating disk in the case of B30-LB14, and with an oblate rotating envelope with infall signatures in the case of LB19. The remaining three ALMA detections do not have infrared counterparts and can be classified as either deeply embedded objects or as starless cores if B30 members. In the former case, two of them (LB08 and LB31) show internal luminosity upper limits consistent with Very Low Luminosity objects, while we do not have enough information for LB10. In the starless core scenario, and taking into account the estimated masses from ALMA and the APEX/LABOCA cores, we estimate final masses for the central objects in the substellar domain, so they could be classified as pre-BD core candidates.

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 study the dynamical state of cores by using a simple analytical model, a sample of observational massive cores, and numerical simulations of collapsing massive cores. From the model, we find that, if cores are formed from turbulent compressions, they evolve from small to large column densities, increasing their velocity dispersion as they collapse, while they tend to equipartition between gravity and kinetic energy. From the observational sample, we find that: (a) cores with substantially different column densities in the sample do not follow a Larson-like linewidth-size relation. Instead, cores with higher column densities tend to be located in the upper-left corner of the Larson velocity dispersion-size diagram, a result predicted previously (Ballesteros-Paredes et al. 2011a). (b) The data exhibit cores with overvirial values. Finally, in the simulations we reproduce the behavior depicted by the model and observational sample: cores evolve towards larger velocity dispersions and smaller sizes as they collapse and increase their column density. However, collapsing cores appear to approach overvirial states within a free-fall time. The cause of this apparent excess of kinetic energy is an underestimation of the gravitational energy, due to the assumption that the gravitational energy is given by the energy of an isolated sphere of constant column density. This excess disappears when the gravitational energy is correctly calculated from the actual spatial mass distribution, where inhomogeneities, as well as the potential due to the mass outside of the core, also contribute to the gravitational energy. We conclude that the observed energy budget of cores in surveys is consistent with their non-thermal motions being driven by their self-gravity and in a hierarchical and chaotic collapse.

We present Submillimeter Array (SMA) observations at 345 GHz towards the intermediate/high-mass cluster-forming region NGC 6334 V. From the dust emission we spatially resolve three dense condensations, the brightest one presenting the typical chemistry of a hot core. The magnetic field (derived from the dust polarized emission) shows a bimodal converging pattern towards the hot core. The molecular emission traces two filamentary structures at two different velocities, separated by 2 km/s, converging to the hot core and following the magnetic field distribution. We compare the velocity field and the magnetic field derived from the SMA observations with MHD simulations of star-forming regions dominated by gravity. This comparison allows us to show how the gas falls in from the larger-scale extended dense core (~0.1 pc) of NGC 6334 V towards the higher-density hot core region (~0.02 pc) through two distinctive converging flows dragging the magnetic field, whose strength seems to have been overcome by gravity.

We present a survey for water maser emission toward a sample of 44 low-luminosity young objects, comprising (proto-)brown dwarfs, first hydrostatic cores (FHCs), and other young stellar objects (YSOs) with bolometric luminosities lower than 0.4 L$_\odot$. Water maser emission is a good tracer of energetic processes, such as mass-loss and/or accretion, and is a useful tool to study this processes with very high angular resolution. This type of emission has been confirmed in objects with L$_{\rm bol}\ge 1$ L$_\odot$. Objects with lower luminosities also undergo mass-loss and accretion, and thus, are prospective sites of maser emission. Our sensitive single-dish observations provided a single detection when pointing toward the FHC L1448 IRS 2E. However, follow-up interferometric observations showed water maser emission associated with the nearby YSO L1448 IRS 2 { (a Class 0 protostar of L$_{\rm bol}\simeq 3.6-5.3$ L$_\odot$)}, and did not find any emission toward L1448 IRS 2E. The upper limits for water maser emission determined by our observations are one order of magnitude lower than expected from the correlation between water maser luminosities and bolometric luminosities found for YSOs. This suggests that this correlation does not hold at the lower end of the (sub)stellar mass spectrum. Possible reasons are that the slope of this correlation is steeper at L$_{\rm bol}\le 1$ L$_\odot$, or that there is an absolute luminosity threshold below which water maser emission cannot be produced. Alternatively, if the correlation still stands at low luminosity, the detection rates of masers would be significantly lower than the values obtained in higher-luminosity Class 0 protostars.

We present sensitive 3.0 cm JVLA radio continuum observations of six regions of low-mass star formation that include twelve young brown dwarfs and four young brown dwarf candidates. We detect a total of 49 compact radio sources in the fields observed, of which 24 have no reported counterparts and are considered new detections. Twelve of the radio sources show variability in timescales of weeks to months, suggesting gyrosynchrotron emission produced in active magnetospheres. Only one of the target brown dwarfs, FU Tau A, was detected. However, we detected radio emission associated with two of the brown dwarf candidates, WL 20S and CHLT 2. The radio flux densities of the sources associated with these brown dwarf candidates are more than an order of magnitude larger than expected for a brown dwarf and suggest a revision of their classification. In contrast, FU Tau A falls on the well-known correlation between radio luminosity and bolometric luminosity, suggesting that the emission comes from a thermal jet and that this brown dwarf seems to be forming as a scaled-down version of low-mass stars.

The early evolution of massive cluster progenitors is poorly understood. We investigate the fragmentation properties from 0.3 pc to 0.06 pc scales of a homogenous sample of infrared-quiet massive clumps within 4.5 kpc selected from the ATLASGAL survey. Using the ALMA 7m array we detect compact dust continuum emission towards all targets, and find that fragmentation, at these scales, is limited. The mass distribution of the fragments uncovers a large fraction of cores above 40 $M_\odot$, corresponding to massive dense cores (MDCs) with masses up to ~400 $M_\odot$. 77 % of the clumps contain at most 3 MDCs per clump, and we also reveal single clumps/MDCs. The most massive cores are formed within the more massive clumps, and a high concentration of mass on small scales reveals a high core formation efficiency. The mass of MDCs highly exceeds the local thermal Jeans-mass, and observational evidence is lacking for a sufficiently high level of turbulence or strong enough magnetic fields to keep the most massive MDCs in equilibrium. If already collapsing, the observed fragmentation properties with a high core formation efficiency are consistent with the collapse setting in at parsec scales.

With the recent recognition of a second, distinctive class of molecular outflows, namely the explosive ones not directly connected to the accretion-ejection process in the star formation, a juxtaposition of the morphological and kinematic properties of both classes is warranted. By applying the same method used in Zapata et al. (2009), and using $^{12}$CO(J=2-1) archival data from the Submillimeter Array (SMA), we contrast two well known explosive objects, Orion KL and DR21, to HH211 and DG Tau B, two flows representative of classical low-mass protostellar outflows. At the moment there are only two well established cases of explosive outflows, but with the full availability of ALMA we expect that more examples will be found in the near future. Main results are the largely different spatial distributions of the explosive flows, consisting of numerous narrow straight filament-like ejections with different orientations and in almost an isotropic configuration, the red with respect to the blueshifted components of the flows (maximally separated in protostellar, largely overlapping in explosive outflows), the very well-defined Hubble flow-like increase of velocity with distance from the origin in the explosive filaments versus the mostly non-organized CO velocity field in protostellar objects, and huge inequalities in mass, momentum and energy of the two classes, at least for the case of low-mass flows. Finally, all the molecular filaments in the explosive outflows point back to approximately a central position i.e. the place where its "exciting source" was located, contrary to the bulk of the molecular material within the protostellar outflows.

We report on subarcsecond observations of complex organic molecules (COMs) in the high-mass protostar IRAS20126+4104 with the Plateau de Bure Interferometer in its most extended configurations. In addition to the simple molecules SO, HNCO and H2-13CO, we detect emission from CH3CN, CH3OH, HCOOH, HCOOCH3, CH3OCH3, CH3CH2CN, CH3COCH3, NH2CN, and (CH2OH)2. SO and HNCO present a X-shaped morphology consistent with tracing the outflow cavity walls. Most of the COMs have their peak emission at the putative position of the protostar, but also show an extension towards the south(east), coinciding with an H2 knot from the jet at about 800-1000 au from the protostar. This is especially clear in the case of H2-13CO and CH3OCH3. We fitted the spectra at representative positions for the disc and the outflow, and found that the abundances of most COMs are comparable at both positions, suggesting that COMs are enhanced in shocks as a result of the passage of the outflow. By coupling a parametric shock model to a large gas-grain chemical network including COMs, we find that the observed COMs should survive in the gas phase for about 2000 yr, comparable to the shock lifetime estimated from the water masers at the outflow position. Overall, our data indicate that COMs in IRAS20126+4104 may arise not only from the disc, but also from dense and hot regions associated with the outflow.

We have performed a dense core survey toward the Infrared Dark Cloud G14.225-0.506 at 3 mm continuum emission with the Atacama Large Millimeter/Submillimeter Array (ALMA). This survey covers the two hub-filament systems with an angular resolution of $\sim3$\arcsec ($\sim0.03$ pc). We identified 48 dense cores. Twenty out of the 48 cores are protostellar due to their association with young stellar objects (YSOs) and/or X-ray point-sources, while the other 28 cores are likely prestellar and unrelated with known IR or X-ray emission. Using APEX 870 $\mu$m continuum emission, we also identified the 18 clumps hosting these cores. Through virial analysis using the ALMA N$_2$H$^+$ and VLA/Effelsberg NH$_3$ molecular line data, we found a decreasing trend in the virial parameter with decreasing scales from filaments to clumps, and then to cores. The virial parameters of $0.1-1.3$ in cores, indicate that cores are likely undergoing dynamical collapse. The cumulative Core Mass Function (CMF) for the prestellar cores candidates has a power law index of $\alpha=1.6$, with masses ranging from 1.5 to 22 $M_\odot$. We find no massive prestellar or protostellar cores. Previous studies suggest that massive O-tpye stars have not been produced yet in this region. Therefore, high-mass stars should be formed in the prestellar cores by accreting a significant amount of gas from the surrounding medium. Another possibility is that low-mass YSOs become massive by accreting from their parent cores that are fed by filaments. These two possibilities might be consistent with the scenario of global hierarchical collapse.

Pre-stellar cores within molecular clouds provide the very initial conditions in which stars are formed. We use the IRAM 30m telescope and the PdBI to study the chemical and physical properties of the starless core FeSt 1-457 (Core 109), in the Pipe nebula. We fit the hyperfine structure of the N2H+(1-0) IRAM 30m data. This allow us to measure with high precision the velocity field, line widths and opacity and derive the excitation temperature and column density in the core. We use a modified Bonnor-Ebert sphere model adding a temperature gradient towards the center to fit the 1.2 mm continuum emission and visual extinction maps. Using this model, we estimate the abundances of the N2H+ and the rest of molecular lines detected in the 30 GHz wide line survey performed at 3 mm with IRAM 30m using ARTIST software. The core presents a rich chemistry with emission from early (C3H2, HCN, CS) and late-time molecules (e.g., N2H+), with a clear chemical spatial differentiation for nitrogen, oxygen and sulphurated molecules. For most of the molecules detected (HCN, HCO+, CH3OH, CS, SO, 13CO and C18O), abundances are best fitted with three values, presenting a clear decrease of abundance of at least 1 or 2 orders of magnitude towards the center of the core. The Bonnor-Ebert analysis indicates the core is gravitationally unstable and the magnetic field is not strong enough to avoid the collapse. Depletion of molecules onto the dust grains occurs at the interior of the core, where dust grain growth and dust depolarization also occurs. This suggests that these properties may be related. On the other hand, some molecules exhibit asymmetries in their integrated emission maps, which appear to be correlated with a previously reported submillimetre polarization asymmetry. These asymmetries could be due to a stronger interstellar radiation field in the western side of the core.

We present sensitive (rms-noises $\sim$ 4 -- 25 $\mu$Jy) and high angular resolution ($\sim$1--2$"$) 8.9 GHz (3.3 cm) Karl G. Jansky Very Large Array (JVLA) radio continuum observations of 10 presumed transitional disks associated with young low-mass stars. We report the detection of radio continuum emission in 5 out of the 10 objects (RXJ1615, UX Tau A, LkCa15, RXJ1633, and SR24s). In the case of LkCa15, the centimeter emission is extended, and has a similar morphology to that of the transitional disk observed at mm wavelengths with an inner depression. For these five detections, we construct the Spectral Energy Distributions (SEDs) from the centimeter to submillimeter wavelengths, and find that they can be well fitted with a single (RXJ 1633 and UX Tau A) or a two component power-law (LkCa 15, RXJ 1615, and SR24s). For the cases where a single power-law fits well the data, the centimeter emission is likely produced by optically thin dust with large grains i.e. centimeter-size pebbles) present in the transitional disks. For the cases where a double power-law fits the data, the centimeter emission might be produced by the combination of photoevaporation and/or a free-free jet. We conclude that RXJ1633 and UX Tau A are excellent examples of transitional disks where the structure of the emission from centimeter/millimeter pebbles can be studied. In the other cases, some other physical emitting mechanisms are also important in the centimeter regime.

Giant molecular clouds contain supersonic turbulence that can locally heat small fractions of gas to over 100 K. We run shock models for low-velocity, C-type shocks propagating into gas with densities between 10^3 and 10^5 cm^(-3) and find that CO lines are the most important cooling lines. Comparison to photodissociation region (PDR) models indicates that mid-J CO lines (J = 8-7 and higher) should be dominated by emission from shocked gas. In Papers I and II we presented CO J = 3-2, 8-7, and 9-8 observations toward four primarily quiescent clumps within infrared dark clouds. Here we fit PDR models to the combined spectral line energy distributions and show that the PDR models that best fit the low-J CO emission underpredict the mid-J CO emission by orders of magnitude, strongly hinting at a hot gas component within these clumps. The low-J CO data clearly show that the integrated intensities of both the CO J = 8-7 and 9-8 lines are anomalously high, such that the line ratio can be used to characterize the hot gas component. Shock models are reasonably consistent with the observed mid-J CO emission, with models with densities near 10^(4.5) cm^(-3) providing the best agreement. Where this mid-J CO is detected, the mean volume filling factor of the hot gas is 0.1%. Much of the observed mid-J CO emission, however, is also associated with known protostars and may be due to protostellar feedback.

Phosphorus is a crucial element for the development of life, but so far P-bearing molecules have been detected only in a few astrophysical objects, hence its interstellar chemistry is almost totally unknown. Here we show new detections of phosphorus nitride in a sample of dense cores in different evolutionary stages of the intermediate- and high-mass star formation process: starless, with protostellar objects, and with ultracompact HII regions. All detected PN line widths are smaller than ~5 km/s , and they arise from regions associated with kinetic temperatures smaller than 100 K. Because the few previous detections reported in the literature are associated with warmer and more turbulent sources, the results of this work show that PN can arise from relatively quiescent and cold gas. This information is challenging for theoretical models that invoke either high desorption temperatures or grain sputtering from shocks to release phosphorus into the gas phase. Derived column densities are of the order of 10^{11-12} cm^{-2}, marginally lower than the values derived in the few high-mass star forming regions detected so far. New constraints on the column density upper limit of phosphorus monoxide, the fundamental unit of biologically relevant molecules, are also given.

The Kleinmann-Low nebula in Orion (Orion-KL) is the nearest example of a high-mass star-forming environment. For the first time, we complemented 1.3 mm Submillimeter Array (SMA) interferometric line survey with IRAM 30 m single-dish observations of the Orion-KL region. Covering a 4 GHz bandwidth in total, this survey contains over 160 emission lines from 20 species (25 isotopologues), including 11 complex organic molecules (COMs). At a spatial resolution of 1200 AU, the continuum substructures are resolved. Extracting the spectra from individual substructures and providing the intensity-integrated distribution map for each species, we studied the small-scale chemical variations in this region. Our main results are: (1) We identify lines from the low-abundance COMs CH3COCH3 and CH3CH2OH, as well as tentatively detect CH3CHO and long carbon-chains C6H and HC7N. (2) We find that while most COMs are segregated by type, peaking either towards the hot core (e.g., N-bearing species) or the compact ridge (e.g., O-bearing species like HCOOCH3 and CH3OCH3), while the distributions of others do not follow this segregated structure (e.g., CH3CH2OH, CH3OH, CH3COCH3). (3) We find a second velocity component of HNCO, SO2, 34SO2, and SO lines, which may be associated with a strong shock event in the low-velocity outflow. (4) Temperatures and molecular abundances show large gradients between central condensations and the outflow regions, illustrating a transition between hot molecular core and shock-chemistry dominated regimes. Our observations of spatially resolved chemical variations in Orion-KL provide the nearest reference source for hot molecular core and outflow chemistry, which will be an important example for interpreting the chemistry of more distant HMSFRs.

We present observations of the 1.3 mm continuum emission toward hub-N and hub-S of the infrared dark cloud G14.225-0.506 carried out with the Submillimeter Array, together with observations of the dust emission at 870 and 350 microns obtained with APEX and CSO telescopes. The large scale dust emission of both hubs consists of a single peaked clump elongated in the direction of the associated filament. At small scales, the SMA images reveal that both hubs fragment into several dust condensations. The fragmentation level was assessed under the same conditions and we found that hub-N presents 4 fragments while hub-S is more fragmented, with 13 fragments identified. We studied the density structure by means of a simultaneous fit of the radial intensity profile at 870 and 350 microns and the spectral energy distribution adopting a Plummer-like function to describe the density structure. The parameters inferred from the model are remarkably similar in both hubs, suggesting that density structure could not be responsible in determining the fragmentation level. We estimated several physical parameters such as the level of turbulence and the magnetic field strength, and we found no significant differences between these hubs. The Jeans analysis indicates that the observed fragmentation is more consistent with thermal Jeans fragmentation compared with a scenario that turbulent support is included. The lower fragmentation level observed in hub-N could be explained in terms of stronger UV radiation effects from a nearby HII region, evolutionary effects, and/or stronger magnetic fields at small scales, a scenario that should be further investigated.

Infrared dark clouds are kinematically complex molecular structures in the interstellar medium that can host sites of massive star formation. We present 4 square arcminute maps of the 12CO, 13CO, and C18O J = 3 to 2 lines from selected locations within the C and F (G028.37+00.07 and G034.43+00.24) infrared dark clouds (IRDCs), as well as single pointing observations of the 13CO and C18O J = 2 to 1 lines towards three cores within these clouds. We derive CO gas temperatures throughout the maps and find that CO is significantly frozen out within these IRDCs. We find that the CO depletion tends to be the highest near column density peaks, with maximum depletion factors between 5 and 9 in IRDC F and between 16 and 31 in IRDC C. We also detect multiple velocity components and complex kinematic structure in both IRDCs. Therefore, the kinematics of IRDCs seem to point to dynamically evolving structures yielding dense cores with considerable depletion factors.

The origin of massive stars is a fundamental open issue in modern astrophysics. Pre-ALMA interferometric studies reveal precursors to early B to late O type stars with collapsing envelopes of 15-20 M$_\odot$ on 1000-3000 AU size-scales. To search for more massive envelopes we selected the most massive nearby young clumps from the ATLASGAL survey to study their protostellar content with ALMA. Our first results using the intermediate scales revealed by the ALMA ACA array providing 3-5" angular resolution, corresponding to $\sim$0.05-0.1 pc size-scales, reveals a sample of compact objects. These massive dense cores are on average two-times more massive than previous studies of similar types of objects. We expect that once the full survey is completed, it will provide a comprehensive view on the origin of the most massive stars.

Context. Chamaeleon II molecular cloud is an active star forming region that offers an excellent opportunity for studying the formation of brown dwarfs in the southern hemisphere. Aims. Our aims are to identify a population of pre- and proto- brown dwarfs (5 sigma mass limit threshold of ~0.015 Msun) and provide information on the formation mechanisms of substellar objects. Methods. We performed high sensitivity observations at 870 microns using the LABOCA bolometer at the APEX telescope towards an active star forming region in Chamaeleon II. The data are complemented with an extensive multiwavelength catalogue of sources from the optical to the far-infrared to study the nature of the LABOCA detections. Results. We detect fifteen cores at 870 microns, and eleven of them show masses in the substellar regime. The most intense objects in the surveyed field correspond to the submillimeter counterparts of the well known young stellar objects DK Cha and IRAS 12500-7658. We identify a possible proto-brown dwarf candidate (ChaII-APEX-L) with IRAC emission at 3.6 and 4.5 microns. Conclusions. Our analysis indicates that most of the spatially resolved cores are transient, and that the point-like starless cores in the sub-stellar regime (with masses between 0.016 Msun and 0.066 Msun) could be pre-brown dwarfs cores gravitationally unstable if they have radii smaller than 220 AU to 907 AU (1.2" to 5" at 178 pc) respectively for different masses. ALMA observations will be the key to reveal the energetic state of these pre-brown dwarfs candidates.

IRAS 22134+5834 was observed in the centimeter with (E)VLA, 3~mm with CARMA, 2~mm with PdBI, and 1.3~mm with SMA, to study the continuum emission as well as the molecular lines, that trace different physical conditions of the gas to study the influence of massive YSOs on nearby starless cores, and the possible implications in the clustered star formation process. The multi-wavelength centimeter continuum observations revealed two radio sources within the cluster, VLA1 and VLA2. VLA1 is considered to be an optically thin UCHII region with a size of 0.01~pc and sits at the edge of the near-infrared (NIR) cluster. The flux of ionizing photons of the VLA1 corresponds to a B1 ZAMS star. VLA2 is associated with an infrared point source and has a negative spectral index. We resolved six millimeter continuum cores at 2~mm, MM2 is associated with the UCHII region VLA1, and other dense cores are distributed around the UCH{\sc ii} region. Two high-mass starless clumps (HMSC), HMSC-E (east) and HMSC-W (west), are detected around the NIR cluster with N$_2$H$^+$(1--0) and NH$_3$ emission, and show different physical and chemical properties. Two N$_2$D$^+$ cores are detected on an NH$_3$ filament close to the UCHII region, with a projected separation of $\sim$8000~AU at the assumed distance of 2.6~kpc. The kinematic properties of the molecular line emission confirm the expansion of the UCHII region and that the molecular cloud around the near infrared (NIR) cluster is also expanding. Our multi-wavelength study has revealed different generations of star formation in IRAS 22134+5834. The formed intermediate- to massive stars show strong impact on nearby starless clumps. We propose that while the stellar wind from the UCHII region and the NIR cluster drives the large scale bubble, the starless clumps and HMPOs formed at the edge of the cluster.

We combine previously published interferometric and single-dish data of relatively nearby massive dense cores that are actively forming stars to test whether their `fragmentation level' is controlled by turbulent or thermal support. We find no clear correlation between the fragmentation level and velocity dispersion, nor between the observed number of fragments and the number of fragments expected when the gravitationally unstable mass is calculated including various prescriptions for `turbulent support'. On the other hand, the best correlation is found for the case of pure thermal Jeans fragmentation, for which we infer a core formation efficiency around 13 per cent, consistent with previous works. We conclude that the dominant factor determining the fragmentation level of star-forming massive dense cores at 0.1 pc scale seems to be thermal Jeans fragmentation.

LRLL 54361 is an infrared source located in the star forming region IC 348 SW. Remarkably, its infrared luminosity increases by a factor of 10 during roughly one week every 25.34 days. To understand the origin of these remarkable periodic variations, we obtained sensitive 3.3 cm JVLA radio continuum observations of LRLL 54361 and its surroundings in six different epochs: three of them during the IR-on state and three during the IR-off state. The radio source associated with LRLL 54361 remained steady and did not show a correlation with the IR variations. We suggest that the IR is tracing the results of fast (with a timescale of days) pulsed accretion from an unseen binary companion, while the radio traces an ionized outflow with an extent of $\sim$100 AU that smooths out the variability over a period of order a year. The average flux density measured in these 2014 observations, 27$\pm$5 $\mu$Jy, is about a factor of two less than that measured about 1.5 years before, $53\pm$11 $\mu$Jy, suggesting that variability in the radio is present, but over larger timescales than in the IR. We discuss other sources in the field, in particular two infrared/X-ray stars that show rapidly varying gyrosynchrotron emission.

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 Square Kilometre Array (SKA) is called to revolutionise essentially all areas of Astrophysics. With a collecting area of about a square kilometre, the SKA will be a transformational instrument, and its scientific potential will go beyond the interests of astronomers. Its technological challenges and huge cost requires a multinational effort, and Europe has recognised this by putting the SKA on the roadmap of the European Strategy Forum for Research Infrastructures (ESFRI). The Spanish SKA White Book is the result of the coordinated effort of 120 astronomers from 40 different research centers. The book shows the enormous scientific interest of the Spanish astronomical community in the SKA and warrants an optimum scientific exploitation of the SKA by Spanish researchers, if Spain enters the SKA project.

We report on the first measurements of the isotopic ratio 14N/15N in N2H+ toward a statistically significant sample of high-mass star forming cores. The sources belong to the three main evolutionary categories of the high-mass star formation process: high-mass starless cores, high-mass protostellar objects, and ultracompact HII regions. Simultaneous measurements of 14N/15N in CN have been made. The 14N/15N ratios derived from N2H+ show a large spread (from ~180 up to ~1300), while those derived from CN are in between the value measured in the terrestrial atmosphere (~270) and that of the proto-Solar nebula (~440) for the large majority of the sources within the errors. However, this different spread might be due to the fact that the sources detected in the N2H+ isotopologues are more than those detected in the CN ones. The 14N/15N ratio does not change significantly with the source evolutionary stage, which indicates that time seems to be irrelevant for the fractionation of nitrogen. We also find a possible anticorrelation between the 14N/15N (as derived from N2H+) and the H/D isotopic ratios. This suggests that 15N enrichment could not be linked to the parameters that cause D enrichment, in agreement with the prediction by recent chemical models. These models, however, are not able to reproduce the observed large spread in 14N/15N, pointing out that some important routes of nitrogen fractionation could be still missing in the models.

There is increasing evidence that episodic accretion is a common phenomenon in Young Stellar Objects (YSOs). Recently, the source HOPS 383 in Orion was reported to have a $\times 35$ mid-infrared -- and bolometric -- luminosity increase between 2004 and 2008, constituting the first clear example of a class 0 YSO (a protostar) with a large accretion burst. The usual assumption that in YSOs accretion and ejection follow each other in time needs to be tested. Radio jets at centimeter wavelengths are often the only way of tracing the jets from embedded protostars. We searched the Very Large Array archive for the available observations of the radio counterpart of HOPS 383. The data show that the radio flux of HOPS 383 varies only mildly from January 1998 to December 2014, staying at the level of $\sim 200$ to 300 $\mu$Jy in the X band ($\sim 9$ GHz), with a typical uncertainty of 10 to 20 $\mu$Jy in each measurement. We interpret the absence of a radio burst as suggesting that accretion and ejection enhancements do not follow each other in time, at least not within timescales shorter than a few years. Time monitoring of more objects and specific predictions from simulations are needed to clarify the details of the connection betwen accretion and jets/winds in YSOs.

We observed with the JVLA at 3.6 and 1.3 cm a sample of 11 proto-brown dwarf candidates in Taurus in a search for thermal radio jets driven by the most embedded brown dwarfs. We detected for the first time four thermal radio jets in proto-brown dwarf candidates. We compiled data from UKIDSS, 2MASS, Spitzer, WISE and Herschel to build the Spectral Energy Distribution (SED) of the objects in our sample, which are similar to typical Class~I SEDs of Young Stellar Objects (YSOs). The four proto-brown dwarf candidates driving thermal radio jets also roughly follow the well-known trend of centimeter luminosity against bolometric luminosity determined for YSOs, assuming they belong to Taurus, although they present some excess of radio emission compared to the known relation for YSOs. Nonetheless, we are able to reproduce the flux densities of the radio jets modeling the centimeter emission of the thermal radio jets using the same type of models applied to YSOs, but with corresponding smaller stellar wind velocities and mass-loss rates, and exploring different possible geometries of the wind or outflow from the star. Moreover, we also find that the modeled mass outflow rates for the bolometric luminosities of our objects agree reasonably well with the trends found between the mass outflow rates and bolometric luminosities of YSOs, which indicates that, despite the "excess" centimeter emission, the intrinsic properties of proto-brown dwarfs are consistent with a continuation of those of very low mass stars to a lower mass range. Overall, our study favors the formation of brown dwarfs as a scaled-down version of low-mass stars.

Water maser emission at 22 GHz is a useful probe to study the transition between the nearly spherical mass-loss in the AGB to a collimated one in the post-AGB phase. In their turn, collimated jets in the post-AGB phase could determine the shape of planetary nebulae (PNe) once photoionization starts. We intend to find new cases of post-AGB stars and PNe with water maser emission, including water fountains or water-maser-emitting PNe. We observed water maser emission in a sample of 133 objects, with a significant fraction being post-AGB and young PN candidate sources with strong obscuration. We detected this emission in 15 of them, of which seven are reported here for the first time. We identified three water fountain candidates: IRAS 17291-2147, with a total velocity spread of ~96 km/s in its water maser components and two sources (IRAS 17021-3109 and IRAS 17348-2906) that show water maser emission outside the velocity range covered by OH masers. We have also identified IRAS 17393-2727 as a possible new water-maser-emitting PN. The detection rate is higher in obscured objects (14%) than in those with optical counterparts (7%), consistent with previous results. Water maser emission seems to be common in objects that are bipolar in the near-IR (43% detection rate). The water maser spectra of water fountain candidates like IRAS 17291-2147 show significantly less maser components than others (e.g., IRAS 18113-2503). We speculate that most post-AGBs may show water maser emission with wide enough velocity spread (> 100 km/s) when observed with enough sensitivity and/or for long enough periods of time. Therefore, it may be necessary to single out a special group of "water fountains", probably defined by their high maser luminosities. We also suggest that the presence of both water and OH masers in a PN is a better tracer of its youth, rather than the presence of just one of these species.

We present high angular resolution ($\sim$ 0.3$"$) submillimeter continuum (0.85 mm) and line observations of the O-type protostar IRAS 16547$-$4247 carried out with the Atacama Large Millimeter/Submillimeter Array (ALMA). In the 0.85 mm continuum band, the observations revealed two compact sources (with a separation of 2$"$), one of them associated with IRAS 16547$-$4247, and the other one to the west. Both sources are well resolved angularly, revealing a clumpy structure. On the other hand, the line observations revealed a rich variety of molecular species related to both continuum sources. In particular, we found a large number of S-bearing molecules, such as the rare molecule methyl mercaptan (CH$_3$SH). At scales larger than 10,000 AU, molecules (e.g., SO$_2$ or OCS) mostly with low excitation temperatures in the upper states (E$_k$ $\lesssim$ 300 K) are present in both millimeter continuum sources, and show a southeast-northwest velocity gradient of 7 km s$^{-1}$ over 3$"$ (165 km s$^{-1}$ pc$^{-1}$). We suggest that this gradient probably is produced by the thermal (free-free) jet emerging from this object with a similar orientation at the base. At much smaller scales (about 1000 AU), molecules with high excitation temperatures (E$_k$ $\gtrsim$ 500 K) are tracing a rotating structure elongated perpendicular to the orientation of the thermal jet, which we interpret as a candidate disk surrounding IRAS 16547$-$4247. The dynamical mass corresponding to the velocity gradient of the candidate to disk is about 20 M$_\odot$, which is consistent with the bolometric luminosity of IRAS 16547$-$4247.

An ever growing number of observational and theoretical evidence suggests that the deuterated fraction (column density ratio between a species containing D and its hydrogenated counterpart, Dfrac) is an evolutionary indicator both in the low- and the high-mass star formation process. However, the role of surface chemistry in these studies has not been quantified from an observational point of view. In order to compare how the deuterated fractions of species formed only in the gas and partially or uniquely on grain surfaces evolve with time, we observed rotational transitions of CH3OH, 13CH3OH, CH2DOH, CH3OD at 3 and 1.3~mm, and of NH2D at 3~mm with the IRAM-30m telescope, and the inversion transitions (1,1) and (2,2) of NH3 with the GBT, towards most of the cores already observed by Fontani et al.~(2011, 2014) in N2H+, N2D+, HNC, DNC. NH2D is detected in all but two cores, regardless of the evolutionary stage. Dfrac(NH3) is on average above 0.1, and does not change significantly from the earliest to the most evolved phases, although the highest average value is found in the protostellar phase (~0.3). Few lines of CH2DOH and CH3OD are clearly detected, and only towards protostellar cores or externally heated starless cores. This work clearly confirms an expected different evolutionary trend of the species formed exclusively in the gas (N2D+ and N2H+) and those formed partially (NH2D and NH3) or totally (CH2DOH and CH3OH) on grain mantles. The study also reinforces the idea that Dfrac(N2H+) is the best tracer of massive starless cores, while high values of Dfrac(CH3OH) seem rather good tracers of the early protostellar phases, at which the evaporation/sputtering of the grain mantles is most efficient.

We report on Submillimeter Array observations of the 870 micron continuum and CO(3-2), 13CO(2-1) and C18O(2-1) line emission of a faint object, SMM2E, near the driving source of the HH797 outflow in the IC348 cluster. The continuum emission shows an unresolved source for which we estimate a mass of gas and dust of 30 Mjup, and the CO(3-2) line reveals a compact bipolar outflow centred on SMM2E, and barely seen also in 13CO(2-1). In addition, C18O(2-1) emission reveals hints of a possible rotating envelope/disk perpendicular to the outflow, for which we infer a dynamical mass of ~16 Mjup. In order to further constrain the accreted mass of the object, we gathered data from Spitzer, Herschel, and new and archive submillimetre observations, and built the Spectral Energy Distribution (SED). The SED can be fitted with one single modified black-body from 70 micron down to 2.1 cm, using a dust temperature of ~24 K, a dust emissivity index of 0.8, and an envelope mass of ~35 Mjup. The bolometric luminosity is 0.10 Lsun, and the bolometric temperature is 35 K. Thus, SMM2E is comparable to the known Class 0 objects in the stellar domain. An estimate of the final mass indicates that SMM2E will most likely remain substellar, and the SMM2E outflow force matches the trend with luminosity known for young stellar objects. Thus, SMM2E constitutes an excellent example of a Class 0 proto-brown dwarf candidate which forms as a scaled-down version of low-mass stars. Finally, SMM2E seems to be part of a wide (~2400 AU) multiple system of Class 0 sources.

We present sensitive 2.1 and 3.3 cm JVLA radio continuum observations of the region IC 348 SW. We detect a total of 10 compact radio sources in the region, of which seven are first reported here. One of the sources is associated with the remarkable periodic time-variable infrared source LRLL 54361, opening the possibility of monitoring this object at radio wavelengths. Four of the sources appear to be powering outflows in the region, including HH 211 and HH 797. In the case of the rotating outflow HH 797 we detect at its center a double radio source, separated by $\sim3"$. Two of the sources are associated with infrared stars that possibly have gyrosynchrotron emission produced in active magnetospheres. Finally, three of the sources are interpreted as background objects.

We carried out a 12CO(3-2) survey of 52 southern stars with a wide range of IR excesses (LIR/L*) using the single dish telescopes APEX and ASTE. The main aims were (1) to characterize the evolution of molecular gas in circumstellar disks using LIR/L* values as a proxy of disk dust evolution, and (2) to identify new gas-rich disk systems suitable for detailed study with ALMA. About 60% of the sample (31 systems) have LIR/L* > 0.01 typical of T-Tauri or Herbig AeBe stars, and the rest (21 systems) have LIR/L* < 0.01 typical of debris disks. We detect CO(3-2) emission from 20 systems, and 18 (90%) of these have LIR/L* > 0.01. However, the spectra of only four of the newly detected systems appear free of contamination from background or foreground emission from molecular clouds. These include the early-type stars HD 104237 (A4/5V, 116 pc) and HD 98922 (A2 III, 507 pc, as determined in this work), where our observations reveal the presence of CO-rich circumstellar disks for the first time. Of the other detected sources, many could harbor gaseous circumstellar disks, but our data are inconclusive. For these two newly discovered gas-rich disks, we present radiative transfer models that simultaneously reproduce their spectral energy distributions and the 12CO(3-2) line profiles. For both of these systems, the data are fit well by geometrically flat disks, placing them in the small class of non-flaring disks with significant molecular gas reservoirs.

NGC 7129 FIRS 2 (hereafter FIRS 2) is an intermediate-mass (2 to 8 Msun) protostar located at a distance of 1250 pc. High spatial resolution observations are required to resolve the hot core at its center. We present a molecular survey from 218200 MHz to 221800 MHz carried out with the IRAM Plateau de Bure Interferometer. These observations were complemented with a long integration single-dish spectrum taken with the IRAM 30m telescope. We used a Local Thermodynamic Equilibrium (LTE) single temperature code to model the whole dataset. The interferometric spectrum is crowded with a total of ~300 lines from which a few dozens remain unidentified yet. The spectrum has been modeled with a total of 20 species and their isomers, isotopologues and deuterated compounds. Complex molecules like methyl formate (CH3OCHO), ethanol (CH3CH2OH), glycolaldehyde (CH2OHCHO), acetone (CH3COCH3), dimethyl ether (CH3OCH3), ethyl cyanide (CH3CH2CN) and the aGg' conformer of ethylene glycol (aGg'-(CH2OH)_2) are among the detected species. The detection of vibrationally excited lines of CH3CN, CH3OCHO, CH3OH, OCS, HC3N and CH3CHO proves the existence of gas and dust at high temperatures. In fact, the gas kinetic temperature estimated from the vibrational lines of CH3CN, ~405 K, is similar to that measured in massive hot cores. Our data allow an extensive comparison of the chemistry in FIRS~2 and the Orion hot core. We find a quite similar chemistry in FIRS 2 and Orion. Most of the studied fractional molecular abundances agree within a factor of 5. Larger differences are only found for the deuterated compounds D2CO and CH2DOH and a few molecules (CH3CH2CN, SO2, HNCO and CH3CHO). Since the physical conditions are similar in both hot cores, only different initial conditions (warmer pre-collapse phase in the case of Orion) and/or different crossing time of the gas in the hot core can explain this behavior.

We present line and continuum observations made with the Submillimeter Array (SMA) of the young stellar object Barnard 5 - IRS1 located in the Perseus molecular cloud. Our $^{12}$CO(2-1) line observations resolve the high-velocity bipolar northeast-southwest outflow associated with this source. We find that the outflowing gas shows different structures at three different velocity regimes, in both lobes, resulting in a spider-like morphology. In addition to the low-velocity, cone-like (wide-angle) lobes that have previously been observed, we report the presence of intermediate-velocity parabolic shells emerging very close to the Class I protostar, as well as high velocity molecular bullets that appear to be associated to the optical/IR jet emanating from this source. These compact high-velocity features reach radial velocities of about 50 km s$^{-1}$ away from the cloud velocity. We interpret the peculiar spider-like morphology is a result of the molecular material being entrained by a wind with both a collimated jet-like component and a wide-angle component. We suggest the outflow is in a transitional evolutionary phase between a mostly jet-driven flow and an outflow in which the entrainment is dominated by the wide-angle wind component. We also detect 1300 $\mu$m continuum emission at the position of the protostar, which likely arises from the dusty envelope and disk surrounding the protostar. Finally, we report the detection of $^{13}$CO(2-1) and SO(6$_5$-5$_4$) emission arising from the outflow and the location of the young stellar object.

In order to study the temperature distribution of the extended gas within the Orion Kleinmann-Low nebula, we have mapped the emission by methyl cyanide (CH3CN) in its J=6_K-5_K, J=12_K-11_K, J=13_K-12_K, and J=14_K-13_K transitions at an average angular resolution of ~10 arcsec (22 arcsec for the 6_K-5_K lines), as part of a new 2D line survey of this region using the IRAM 30m telescope. These fully sampled maps show extended emission from warm gas to the northeast of IRc2 and the distinct kinematic signatures of the hot core and compact ridge source components. We have constructed population diagrams for the four sets of K-ladder emission lines at each position in the maps and have derived rotational excitation temperatures and total beam-averaged column densities from the fitted slopes. In addition, we have fitted LVG model spectra to the observations to determine best-fit physical parameters at each map position, yielding the distribution of kinetic temperatures across the region. The resulting temperature maps reveal a region of hot (T > 350 K) material surrounding the northeastern edge of the hot core, whereas the column density distribution is more uniform and peaks near the position of IRc2. We attribute this region of hot gas to shock heating caused by the impact of outflowing material from active star formation in the region, as indicated by the presence of broad CH3CN lines. This scenario is consistent with predictions from C-shock chemical models that suggest that gas-phase methyl cyanide survives in the post-shock gas and can be somewhat enhanced due to sputtering of grain mantles in the passing shock front.

In order to shed light on the main physical processes controlling fragmentation of massive dense cores, we present a uniform study of the density structure of 19 massive dense cores, selected to be at similar evolutionary stages, for which their relative fragmentation level was assessed in a previous work. We inferred the density structure of the 19 cores through a simultaneous fit of the radial intensity profiles at 450 and 850 micron (or 1.2 mm in two cases) and the Spectral Energy Distribution, assuming spherical symmetry and that the density and temperature of the cores decrease with radius following power-laws. We find a weak (inverse) trend of fragmentation level and density power-law index, with steeper density profiles tending to show lower fragmentation, and vice versa. In addition, we find a trend of fragmentation increasing with density within a given radius, which arises from a combination of flat density profile and high central density and is consistent with Jeans fragmentation. We considered the effects of rotational-to-gravitational energy ratio, non-thermal velocity dispersion, and turbulence mode on the density structure of the cores, and found that compressive turbulence seems to yield higher central densities. Finally, a possible explanation for the origin of cores with concentrated density profiles, which are the cores showing no fragmentation, could be related with a strong magnetic field, consistent with the outcome of radiation magnetohydrodynamic simulations.

G79.29+0.46 is a nebula created by a Luminous Blue Variable (LBV) star candidate characterized by two almost circular concentric shells. In order to investigate whether the shells are interacting with the infrared dark cloud (IRDC) G79.3+0.3 located at the southwestern border of the inner shell, we conducted Jansky Very Large Array observations of NH3(1,1), (2,2) and c-C3H2, and combined them with previous Effelsberg data. The overall NH3 emission consists of one main clump, named G79A, elongated following the shape of the IRDC, plus two fainter and smaller cores to the north, which spatially match the inner infrared shell. We analysed the NH3 spectra at each position with detected emission and inferred linewidth, rotational temperature, column density and abundance maps, and find that: i) the linewidth of NH3(1,1) in the northern cores is 0.5 km/s, slightly larger than in their surroundings; ii) the NH3 abundance is enhanced by almost one order of magnitude towards the northwestern side of G79A; iii) there is one `hot slab' at the interface between the inner infrared shell and the NH3 peak of G79A; iv) the western and southern edges of G79A present chemical differentiation, with c-C3H2 tracing more external layers than NH3, similar to what is found in PDRs. Overall, the kinematics and physical conditions of G79A are consistent with both shock-induced and UV radiation-induced chemistry driven by the LBV star. Therefore, the IRDC is not likely associated with the star-forming region DR15, but located farther away, near G79.29+0.46 at 1.4 kpc.

The surroundings of Luminous Blue Variable (LBV) stars are excellent laboratories to study the effects of their high UV radiation, powerful winds, and strong ejection events onto the surrounding gas and dust. The LBV G79.29+0.46 powered two concentric infrared rings which may interact with the infrared dark cloud (IRDC) G79.3+0.3. The Effelsberg 100m telescope was used to observe the NH_3 (1,1), (2,2) emission surrounding G79.29+0.46 and the IRDC. In addition, we observed particular positions in the (3,3) transition toward the strongest region of the IRDC. We report here the first coherent shell-like structure of dense NH_3 gas associated with an evolved massive star. The shell, two or three orders of magnitude more tenuous than the IRDC, is well traced in both ammonia lines, and surrounds the ionized nebula. The NH_3 emission in the IRDC is characterized by a low and uniform rotational temperature (T_rot ~ 10 K) and moderately high opacities in the (1,1) line. The rest of the observed field is spotted by warm or hot zones (T_rot > 30 K) and characterized by optically thin emission of the (1,1) line. The NH_3 abundances are about 10^{-8} in the IRDC, and 10^{-10}-10^{-9} elsewhere. The warm temperatures and low abundances of NH_3 in the shell suggest that the gas is being heated and photo-dissociated by the intense UV field of the LBV star. An outstanding region is found to the south-west (SW) of the LBV star within the IRDC. The NH_3 (3,3) emission at the centre of the SW region reveals two velocity components tracing gas at temperatures > 30K. The northern edge of the SW region agrees with the border of the ring nebula and a region of continuum enhancement; here, the opacity of the (1,1) line and the NH_3 abundance do not decrease as expected in a typical clump of an isolated cold dark cloud. This strongly suggests some kind of interaction between the ring nebula and the IRDC.

We present CO 3-2, SiO 8-7, C34S 7-6, and 878 mum dust continuum subarcsecond angular resolution observations with the SMA toward IRAS 16293-2422 (I16293). The C34S emission traces the 878 mum dust continuum well, and clearly shows a smooth velocity gradient along the major axis of component I16293A. The CO shows emission at moderate high velocities arising from two bipolar outflows, which appear to be perpendicular with respect to each other. The high sensitivity and higher angular resolution of these observations allows us to pinpoint well the origin of these two outflows at the center of component I16293A. Interestingly, the most compact outflow appears to point toward I16293B. Our data show that the previously reported monopolar blueshifted CO outflow associated with component I16293B seems to be part of the compact outflow arising from component I16293A. In addition, the SiO emission is also tracing this compact outflow: on one hand, the SiO emission appears to have a jet-like morphology along the southern redshifted lobe; on the other hand, the SiO emission associated with the blueshifted northern lobe traces a well defined arc on the border of component I16293B facing I16293A. The blueshifted CO lobe of the compact outflow splits into two lobes around the position of this SiO arc. All these results lead us to propose that the compact outflow from component I16293A is impacting on the circumstellar gas around component I16293B, possibly being diverged as a consequence of the interaction.

We present a study of cyanoacetylene (HC3N) and cyanodiacetylene (HC5N) in Orion KL, through observations from two line surveys performed with the IRAM 30m telescope and the HIFI instrument on board the Herschel telescope. The frequency ranges covered are 80-280 GHz and 480-1906 GHz. We model the observed lines of HC3N, HC5N, their isotopologues (including DC3N), and vibrational modes, using a non-LTE radiative transfer code. To investigate the chemical origin of HC3N and DC3N in Orion KL, we use a time-dependent chemical model. We detect 40 lines of the ground state of HC3N and 68 lines of its 13C isotopologues. We also detect 297 lines of six vibrational modes of this molecule (nu_7, 2nu_7, 3nu_7, nu_6, nu_5, and nu_6+nu_7) and 35 rotational lines of the ground state of HC5N. We report the first tentative detection of DC3N in a giant molecular cloud with a DC3N/HC3N abundance ratio of 0.015. We provide column densities and isotopic and molecular abundances. We also perform a 2x2" map around Orion IRc2 and we present maps of HC3N lines and maps of lines of the HC3N vibrational modes nu_6 and nu_7. In addition, a comparison of our results for HC3N with those in other clouds allows us to derive correlations between the column density, the FWHM, the mass, and the luminosity of the clouds. The high column densities of HC3N obtained in the hot core, make this molecule an excellent tracer of hot and dense gas. In addition, the large frequency range covered reveals the need to consider a temperature and density gradient in the hot core in order to obtain better line fits. The high D/H ratio (comparable to that obtained in cold clouds) that we derive suggests a deuterium enrichment. Our chemical models indicate that the possible deuterated HC3N present in Orion KL is formed during the gas-phase. This fact provides new hints concerning the processes leading to deuteration.

We present a study of the sulfur-bearing species detected in a line confusion-limited survey towards Orion KL performed with the IRAM 30m telescope in the range 80-281 GHz. The study is part of an analysis of the line survey divided into families of molecules. Our aim is to derive accurate physical conditions and molecular abundances in the different components of Orion KL from observed SO and SO2 lines. First we assumed LTE conditions obtain rotational temperatures. We then used a radiative transfer model, assuming either LVG or LTE excitation to derive column densities of these molecules in the different components of Orion KL. We have detected 68 lines of SO, 34SO, 33SO, and S18O and 653 lines of SO2, 34SO2, 33SO2, SO18O and SO2 v2=1. We provide column densities for all of them and also upper limits for the column densities of S17O, 36SO, 34S18O, SO17O and 34SO2 v2=1 and for several undetected sulfur-bearing species. In addition, we present 2'x2' maps around Orion IRc2 of SO2 transitions with energies from 19 to 131 K and also maps with four transitions of SO, 34SO and 34SO2. We observe an elongation of the gas along the NE-SW direction. An unexpected emission peak appears at 20.5 km/s in most lines of SO and SO2. A study of the spatial distribution of this emission feature shows that it is a new component ~5" in diameter, which lies ~4" west of IRc2. We suggest the emission from this feature is related to shocks associated to the BN object. The highest column densities for SO and SO2 are found in the high-velocity plateau (a region dominated by shocks) and in the hot core. These values are up to three orders of magnitude higher than the results for the ridge components. We also find high column densities for their isotopologues in both components. Therefore, we conclude that SO and SO2 are good tracers, not only of regions affected by shocks, but also of regions with warm dense gas.

We aim at characterising dense cores in the clustered environments associated with massive star-forming regions. For this, we present an uniform analysis of VLA NH3(1,1) and (2,2) observations towards a sample of 15 massive star-forming regions, where we identify a total of 73 cores, classify them as protostellar, quiescent starless, or perturbed starless, and derive some physical properties. The average sizes and ammonia column densities are 0.06 pc and 10^15 cm^-2, respectively, with no significant differences between the starless and protostellar cores, while the linewidth and rotational temperature of quiescent starless cores are smaller, 1.0 km/s and 16 K, than those of protostellar (1.8 km/s, 21 K), and perturbed starless (1.4 km/s, 19 K) cores. Such linewidths and temperatures for these quiescent starless cores in the surroundings of massive stars are still significantly larger than the typical values measured in starless cores of low-mass star-forming regions, implying an important non-thermal component. We confirm at high angular resolutions the correlations previously found with single-dish telescopes between the linewidth, the temperature of the cores, and the bolometric luminosity. In addition, we find a correlation between the temperature of each core and the incident flux from the most massive star in the cluster, suggesting that the large temperatures measured in the starless cores of our sample could be due to heating from the nearby massive star. A simple virial equilibrium analysis seems to suggest a scenario of a self-similar, self-graviting, turbulent, virialised hierarchy of structures from clumps (0.1-10 pc) to cores (0.05 pc). A closer inspection of the dynamical state taking into account external pressure effects, reveal that relatively strong magnetic field support may be needed to stabilise the cores, or that they are unstable and thus on the verge of collapse.

We report on the long term X-ray monitoring of the outburst decay of the low magnetic field magnetar SGR 0418+5729, using all the available X-ray data obtained with RXTE, SWIFT, Chandra, and XMM-Newton observations, from the discovery of the source in June 2009, up to August 2012. The timing analysis allowed us to obtain the first measurement of the period derivative of SGR 0418+5729: \dot{P}=4(1)x10^{-15} s/s, significant at ~3.5 sigma confidence level. This leads to a surface dipolar magnetic field of B_dip ~6x 10^{12} G. This measurement confirms SGR 0418+5729 as the lowest magnetic field magnetar. Following the flux and spectral evolution from the beginning of the outburst up to ~1200 days, we observe a gradual cooling of the tiny hot spot responsible for the X-ray emission, from a temperature of ~0.9 to 0.3 keV. Simultaneously, the X-ray flux decreased by about 3 orders of magnitude: from about 1.4x10^{-11} to 1.2x10^{-14} erg/s/cm^2 . Deep radio, millimeter, optical and gamma-ray observations did not detect the source counterpart, implying stringent limits on its multi-band emission, as well as constraints on the presence of a fossil disk. By modeling the magneto-thermal secular evolution of SGR 0418+5729, we infer a realistic age of ~550 kyr, and a dipolar magnetic field at birth of ~10^{14} G. The outburst characteristics suggest the presence of a thin twisted bundle with a small heated spot at its base. The bundle untwisted in the first few months following the outburst, while the hot spot decreases in temperature and size. We estimate the outburst rate of low magnetic field magnetars to be about one per year per galaxy, and we briefly discuss the consequences of such result in several other astrophysical contexts.

We present sub-arcsecond observations toward the massive star forming region G75.78+0.34. We used the Very Large Array to study the centimeter continuum and H2O and CH3OH maser emission, and the Owens Valley Radio Observatory and Submillimeter Array to study the millimeter continuum and recombination lines (H40alpha and H30alpha). We found radio continuum emission at all wavelengths, coming from three components: (1) a cometary ultracompact (UC) HII region with an electron density 3.7x10^4 cm^{-3}, excited by a B0 type star, and with no associated dust emission; (2) an almost unresolved UCHII region (EAST), located 6 arcsec to the east of the cometary UCHII region, with an electron density 1.3x10^5 cm^{-3}, and associated with a compact dust clump detected at millimeter and mid-infrared wavelengths; and (3) a compact source (CORE), located 2 arcsec to the southwest of the cometary arc, with a flux density increasing with frequency, and embedded in a dust condensation of 30 Msun. The CORE source is resolved into two compact and unresolved sources which can be well-fit by two homogeneous hypercompact HII regions each one photo-ionized by a B0.5 ZAMS star, or by free-free radiation from shock-ionized gas resulting from the interaction of a jet/outflow system with the surrounding environment. The spatial distribution and kinematics of water masers close to the CORE-N and S sources, together with excess emission at 4.5 mum and the detected dust emission, suggest that the CORE source is a massive protostar driving a jet/outflow.

We present the results of combined NH3(1,1) and (2,2) line emission observed with the Very Large Array and the Effelsberg 100m telescope of the Infrared Dark Cloud G14.225-0.506. The NH3 emission reveals a network of filaments constituting two hub-filament systems. Hubs are associated with gas of rotational temperature Trot \sim 25 K, non-thermal velocity dispersion ~1.1 km/s, and exhibit signs of star formation, while filaments appear to be more quiescent (Trot \sim 11 K, non-thermal velocity dispersion ~0.6 km/s). Filaments are parallel in projection and distributed mainly along two directions, at PA \sim 10 deg and 60 deg, and appear to be coherent in velocity. The averaged projected separation between adjacent filaments is between 0.5 pc and 1pc, and the mean width of filaments is 0.12 pc. Cores within filaments are separated by ~0.33 pc, which is consistent with the predicted fragmentation of an isothermal gas cylinder due to the 'sausage'-type instability. The network of parallel filaments observed in G14.225-0.506 is consistent with the gravitational instability of a thin gas layer threaded by magnetic fields. Overall, our data suggest that magnetic fields might play an important role in the alignment of filaments, and polarization measurements in the entire cloud would lend further support to this scenario.

In order to study the fragmentation of massive dense cores, which constitute the cluster cradles, we observed with the PdBI in the most extended configuration the continuum at 1.3 mm and the CO(2-1) emission of four massive cores. We detect dust condensations down to ~0.3 Msun and separate millimeter sources down to 0.4" or ~1000 AU, comparable to the sensitivities and separations reached in optical/infrared studies of clusters. The CO(2-1) high angular resolution images reveal high-velocity knots usually aligned with previously known outflow directions. This, in combination with additional cores from the literature observed at similar mass sensitivity and spatial resolution, allowed us to build a sample of 18 protoclusters with luminosities spanning 3 orders of magnitude. Among the 18 regions, ~30% show no signs of fragmentation, while 50% split up into ~4 millimeter sources. We compiled a list of properties for the 18 massive dense cores, such as bolometric luminosity, total mass, and mean density, and found no correlation of any of these parameters with the fragmentation level. In order to investigate the combined effects of magnetic field, radiative feedback and turbulence in the fragmentation process, we compared our observations to radiation magneto-hydrodynamic simulations, and found that the low-fragmented regions are well reproduced in the magnetized core case, while the highly-fragmented regions are consistent with cores where turbulence dominates over the magnetic field. Overall, our study suggests that the fragmentation in massive dense cores could be determined by the initial magnetic field/turbulence balance in each particular core.

The theory of massive star formation currently suffers from a scarce observational base of massive young stellar objects to compare with. In this paper, we present OVRO CO(1-0), 13CO(1-0), C18O(1-0), and 2.6 mm continuum images of the infrared source IRAS 19520+2759 together with complementary single-dish observations of CS(1-0), obtained with the 34m antenna DSS-54 at the Madrid Deep Space Communications Complex, as well as archive images at different wavelengths. As a result from our work, IRAS 19520+2759, with a controversial nature in the past, is firmly established as a massive young stellar object associated with a strong and compact millimetre source and driving a collimated outflow. In addition, a second fainter millimetre source is discovered about 4 arcsec to the south, which is also driving an outflow. Furthermore, the two millimetre sources are associated with C18O clumps elongated perpendicularly to the outflows, which may be related to rotating toroids. The masses of gas and dust of the millimetre sources are estimated to be around 100 and 50 Msun. MM1, the dominant source at all wavelengths, with a total luminosity of (1-2)x10^5 Lsun at 9 kpc, is however not associated with 6 cm emission down to a rms noise level of 0.1 mJy. We propose that IRAS 19520+2759 could be an example of the recent theoretical prediction of `bloated' or `swollen' star, i.e., a massive young stellar object whose radius has increased due to effects of accretion at a high-mass accretion rate.

In an attempt to study whether the formation of brown dwarfs (BDs) takes place as a scaled-down version of low-mass stars, we conducted IRAM30m/MAMBO-II observations at 1.2 mm in a sample of 12 proto-BD candidates selected from Spitzer/IRAC data in the B213-L1495 clouds in Taurus. Subsequent observations with the CSO at 350 micron, VLA at 3.6 and 6 cm, and IRAM30m/EMIR in the 12CO(1-0), 13CO(1-0), and N2H+(1-0) transitions were carried out toward the two most promising Spitzer/IRAC source(s), J042118 and J041757. J042118 is associated with a compact (<10 arcsec or <1400 AU) and faint source at 350 micron, while J041757 is associated with a partially resolved (~16 arcsec or ~2000 AU) and stronger source emitting at centimetre wavelengths with a flat spectral index. The corresponding masses of the dust condensations are ~1 and ~5 Mjup for J042118 and J041757, respectively. In addition, about 40 arcsec to the northeast of J041757 we detect a strong and extended submillimetre source, J041757-NE, which is not associated with NIR/FIR emission down to our detection limits, but is clearly detected in 13CO and N2H+ at ~7 km/s, and for which we estimated a total mass of ~100 Mjup, close to the mass required to be gravitationally bound. In summary, our observational strategy has allowed us to find in B213-L1495 two proto-BD candidates and one pre-substellar core candidate, whose properties seem to be consistent with a scaled-down version of low-mass stars.

We report on high-sensitivity and high-angular resolution archival Submillimeter Array (SMA) observations of the large ($\sim$15000 AU) putative circumstellar disk associated with the O-type protostar NGC 7538 S. Observations of the continuum resolve this putative circumstellar disk into five compact sources, with sizes $\sim$ 3000 AU and masses $\sim10 M_\odot$. This confirm the results of recent millimeter observations made with CARMA/BIMA towards this object. However, we find that from most of these compact sources eject collimated bipolar outflows, revealed by our silicon monoxide (SiO {\it J}=5-4) observations and confirm that these sources have a (proto)stellar nature. All outflows are perpendicular to the large and rotating dusty structure. We propose therefore that, rather than being a single massive circumstellar disk, NGC 7538 S could be instead a large and massive contracting or rotating filament that is fragmenting at scales of 0.1 to 0.01 pc to form several B-type stars, via the standard process involving outflows and disks. As in recent high spatial resolution studies of dusty filaments, our observations also suggest that thermal pressure does not seem to be sufficient to support the filament, so that either additional support needs to be invoked, or else the filament must be in the process of collapsing. An SPH numerical simulation of the formation of a molecular cloud by converging warm neutral medium flows produces contracting filaments whose dimensions and spacings between the stars forming within them, as well as their column densities, strongly resemble those observed in the filament reported here.

The intermediate- to high-mass star-forming region IRAS 20343+4129 is an excellent laboratory to study the influence of high- and intermediate-mass young stellar objects on nearby starless dense cores, and investigate for possible implications in the clustered star formation process. We present 3 mm observations of continuum and rotational transitions of several molecular species (C2H, c-C3H2, N2H+, NH2D) obtained with the Combined Array for Research in Millimetre-wave Astronomy, as well as 1.3 cm continuum and NH3 observations carried out with the Very Large Array, to reveal the properties of the dense gas. We confirm undoubtedly previous claims of an expanding cavity created by an ultracompact HII region associated with a young B2 zero-age main sequence (ZAMS) star. The dense gas surrounding the cavity is distributed in a filament that seems squeezed in between the cavity and a collimated outflow associated with an intermediate-mass protostar. We have identified 5 millimeter continuum condensations in the filament. All of them show column densities consistent with potentially being the birthplace of intermediate- to high-mass objects. These cores appear different from those observed in low-mass clustered environments in sereval observational aspects (kinematics, temperature, chemical gradients), indicating a strong influence of the most massive and evolved members of the protocluster. We suggest a possible scenario in which the B2 ZAMS star driving the cavity has compressed the surrounding gas, perturbed its properties and induced the star formation in its immediate surroundings.

We have mapped at small spatial scales the temperature and the velocity field in the protocluster associated with IRAS 05345+3157, which contains both intermediate-/high-mass protostellar candidates and starless condensations, and is thus an excellent location to investigate the role of massive protostars on protocluster evolution. We observed the ammonia (1,1) and (2,2) inversion transitions with the VLA. Ammonia is the best thermometer for dense and cold gas, and the observed transitions have critical densities able to trace the kinematics of the intracluster gaseous medium. The ammonia emission is extended and distributed in two filamentary structures. The starless condensations are colder than the star-forming cores, but the gas temperature across the whole protocluster is higher (by a factor of ~1.3-1.5) than that measured typically in both infrared dark clouds and low-mass protoclusters. The non-thermal contribution to the observed line broadening is at least a factor of 2 larger than the expected thermal broadening even in starless condensations, contrary to the close-to-thermal line widths measured in low-mass quiescent dense cores. The NH3-to-N2H+ abundance ratio is greatly enhanced (a factor of 10) in the pre--stellar core candidates, probably due to freeze-out of most molecular species heavier than He. The more massive and evolved objects likely play a dominant role in the physical properties and kinematics of the protocluster. The high level of turbulence and the fact that the measured core masses are larger than the expected thermal Jeans masses indicate that turbulence likely was an important factor in the initial fragmentation of the parental clump.

Observations with the Plateau de Bure Interferometer in the most extended configuration toward two intermediate-mass star-forming regions, IRAS22198+6336 and AFGL5142, reveal the presence of several complex organic molecules at ~500 AU scales, confirming the presence of hot cores in both regions. The hot cores are not rich in CN-bearing molecules, as often seen in massive hot cores, and are mainly traced by CH3CH2OH, (CH2OH)2, CH3COCH3, and CH3OH, with additionally CH3CHO, CH3OD and HCOOD for IRAS22198+6336, and C6H, and O13CS for AFGL5142. The emission of complex molecules is resolved down to sizes of ~300 and ~600 AU, for IRAS22198+6336 and AFGL5142, respectively, and most likely is tracing protostellar disks rather than flattened envelopes or toroids as usually found. This is specially clear for the case of IRAS22198+6336, where we detect a velocity gradient for all the mapped molecules perpendicular to the most chemically rich outflow of the region, yielding a dynamic mass >4 Msun. As for AFGL5142, the hot core emission is resolved into two elongated cores separated 1800 AU. A detailed comparison of the complex molecule peaks to the new CO(2-1) data and H2O maser data from literature suggests that also for AFGL5142 the complex molecules are mainly associated with disks, except for a faint and extended molecular emission found to the west, which is possibly produced in the interface between one of the outflows and the dense surrounding gas.

Theory predicts, and observations confirm, that the column density ratio of a molecule containing D to its counterpart containing H can be used as an evolutionary tracer in the low-mass star formation process. Since it remains unclear if the high-mass star formation process is a scaled-up version of the low-mass one, we investigated whether the relation between deuteration and evolution can be applied to the high-mass regime. With the IRAM-30m telescope, we observed rotational transitions of N2D+ and N2H+ and derived the deuterated fraction in 27 cores within massive star-forming regions understood to represent different evolutionary stages of the massive-star formation process. Results. Our results clearly indicate that the abundance of N2D+ is higher at the pre-stellar/cluster stage, then drops during the formation of the protostellar object(s) as in the low-mass regime, remaining relatively constant during the ultra-compact HII region phase. The objects with the highest fractional abundance of N2D+ are starless cores with properties very similar to typical pre-stellar cores of lower mass. The abundance of N2D+ is lower in objects with higher gas temperatures as in the low-mass case but does not seem to depend on gas turbulence. Our results indicate that the N2D+-to-N2H+ column density ratio can be used as an evolutionary indicator in both low- and high-mass star formation, and that the physical conditions influencing the abundance of deuterated species likely evolve similarly during the processes that lead to the formation of both low- and high-mass stars.

We aim at investigating with high angular resolution the NH3/N2H+ ratio toward the high-mass star-forming region AFGL 5142 in order to study whether this ratio behaves similarly to the low-mass case, for which the ratio decreases from starless cores to cores associated with YSOs. CARMA was used to observe the 3.2 mm continuum and N2H+(1-0) emission. We used NH3(1,1) and (2,2), HCO+(1-0) and H13CO+(1-0) data from the literature and we performed a time-dependent chemical modeling of the region. The 3.2 mm continuum emission reveals a dust condensation of ~23 Msun associated with the massive YSOs, deeply embedded in the strongest NH3 core (hereafter central core). The N2H+ emission reveals two main cores, the western and eastern core, located to the west and to the east of the mm condensation, and surrounded by a more extended and complex structure of ~0.5 pc. Toward the central core the N2H+ emission drops significantly, indicating a clear chemical differentiation in the region. We found low values of the NH3/N2H+ ratio ~50-100 toward the western/eastern cores, and high values up to 1000 in the central core. The chemical model indicates that density, and in particular temperature, are key parameters in determining the NH3/N2H+ ratio. The high density and temperature reached in the central core allow molecules like CO to evaporate from grain mantles. The CO desorption causes a significant destruction of N2H+, favoring the formation of HCO+. This result is supported by our observations, which show that N2H+ and HCO+ are anticorrelated in the central core. The observed values of the NH3/N2H+ ratio in the central core can be reproduced by our model for times t~4.5-5.3x10^5 yr (central) and t~10^4-3x10^6 yr (western/eastern). The NH3/N2H+ ratio in AFGL 5142 does not follow the same trend as in regions of low-mass star formation mainly due to the high temperature reached in hot cores.

We present new SMA and PdBI observations of the intermediate-mass object IRAS 22198+6336 in the millimeter continuum and in several molecular line transitions. The millimeter continuum emission reveals a strong and compact source with a mass of around 5 Msun and with properties of Class 0 objects. CO emission shows an outflow with a quadrupolar morphology centered on the position of the dust condensation. The CO outflow emission seems to come from two distinct outflows, one of them associated with SiO outflow emission. A large set of molecular lines has been detected toward a compact dense core clearly coincident with the compact millimeter source, and showing a velocity gradient perpendicular to the outflow traced by CO and SiO. The chemically rich spectrum and the rotational temperatures derived from CH$_3$CN and CH$_3$OH (100-150 K) indicate that IRAS 22198+6336 is harbouring one the few intermediate-mass hot cores known at present.

The deuterium fractionation, Dfrac, has been proposed as an evolutionary indicator in pre-protostellar and protostellar cores of low-mass star-forming regions. We investigate Dfrac, with high angular resolution, in the cluster environment surrounding the UCHII region IRAS 20293+3952. We performed high angular resolution observations with the IRAM Plateau de Bure Interferometer (PdBI) of the ortho-NH2D 1_{11}-1_{01} line at 85.926 GHz and compared them with previously reported VLA NH3 data. We detected strong NH2D emission toward the pre-protostellar cores identified in NH3 and dust emission, all located in the vicinity of the UCHII region IRAS 20293+3952. We found high values of Dfrac~0.1-0.8 in all the pre-protostellar cores and low values, Dfrac<0.1, associated with young stellar objects. The high values of Dfrac in pre-protostellar cores could be indicative of evolution, although outflow interactions and UV radiation could also play a role.

Aiming to perform a study of the warm dust and gas in the luminous blue variable star G79.29+0.46 and its associated nebula, we present infrared Spitzer imaging and spectroscopy, and new CO J=2-->1 and 4-->3 maps obtained with the IRAM 30m radio telescope and with the Submillimeter Telescope, respectively. We have analyzed the nebula detecting multiple shells of dust and gas connected to the star. Using Infrared Spectrograph-Spitzer spectra, we have compared the properties of the central object, the nebula, and their surroundings. These spectra show a rich variety of solid-state features (amorphous silicates, polycyclic aromatic hydrocarbons, and CO2 ices) and narrow emission lines, superimposed on a thermal continuum. We have also analyzed the physical conditions of the nebula, which point to the existence of a photo-dissociation region.

HH 211 is a highly collimated jet originating from a nearby young Class 0 protostar. Here is a follow-up study of the jet with our previous observations at unprecedented resolution up to ~ 0.3" in SiO (J=8-7), CO (J=3-2), and SO (N_J=8_9-7_8). SiO, CO, and SO can all be a good tracer of the HH 211 jet, tracing the internal shocks in the jet. Although the emissions of these molecules show roughly the same morphology of the jet, there are detailed differences. In particular, the CO emission traces the jet closer to the source than the SiO and SO emissions. In addition, in the better resolved internal shocks, both the CO and SO emission are seen slightly ahead of the SiO emission. The jet is clearly seen on both sides of the source with more than one cycle of wiggle. The wiggle is reflection-symmetric about the source and can be reasonably fitted by an orbiting source jet model. The best-fit parameters suggest that the source itself could be a very low-mass protobinary with a total mass of ~ 60 M_Jup and a binary separation of ~ 4.6 AU. The abundances of SiO and SO in the gas phase are found to be highly enhanced in the jet as compared to the quiescent molecular clouds, even close to within 300 AU from the source where the dynamical time scale is <10 yrs. The abundance enhancements of these molecules are closely related to the internal shocks. The detected SiO is either the consequence of the release of Si-bearing material from dust grains or of its formation via gas chemistry in the shocks. The SO, on the other hand, seems to form via gas chemistry in the shocks.

We observed with the VLA, PdBI, and SMA the centimeter and millimeter continuum, N2H+(1-0), and CO(2-1) emission associated with a dusty cloud harboring a nascent cluster with intermediate-mass protostars. At centimeter wavelengths we found a strong source, tracing a UCHII region, at the eastern edge of the dusty cloud, with a shell-like structure, and with the near-infrared counterpart falling in the center of the shell. This is presumably the most massive source of the forming cluster. About 15'' to the west of the UCHII region and well embedded in the dusty cloud, we detected a strong millimeter source, MM1, associated with centimeter and near-infrared emission. MM1 seems to be driving a prominent high-velocity CO bipolar outflow, and is embedded in a ridge of dense gas traced by N2H+. We estimated that MM1 is an intermediate-mass source in the Class 0/I phase. About 15'' to the south of MM1, and still more deeply embedded in the dusty cloud, we detected a compact millimeter source, MM2, with neither centimeter nor near-infrared emission, but with water maser emission. MM2 is associated with a clump of N2H+, whose kinematics reveal a clear velocity gradient and additionally we found signposts of infall motions. MM2, being deeply embedded within the dusty cloud, with an associated water maser but no hints of CO outflow emission, is an intriguing object, presumably of intermediate mass. In conclusion, the UCHII region is found at the border of a dusty cloud which is currently undergoing active star formation. Two intermediate-mass protostars in the dusty cloud seem to have formed after the UCHII region and have different properties related to the outflow phenomenon.

We aim at studying with high angular resolution a dense core associated with a low-luminosity IRAS source, IRAS 00213+6530, in order to investigate whether low mass star formation is really taking place in isolation. We performed observations at 1.2mm with the IRAM 30m telescope, VLA observations at 6cm, 3.6cm, 1.3cm, 7mm, and H2O maser and NH3 lines, and observations with the NASA 70m antenna in CCS and H2O maser. The cm and mm continuum emission, together with the near infrared data from the 2MASS allowed us to identify 3 YSOs, IRS1, VLA8A, and VLA8B, with different radio and infrared properties, and which seem to be in different evolutionary stages. The NH3 emission consists of three clouds. Two of these, MM1 and MM2, are associated with dust emission, while the southern cloud is only detected in NH3. The YSOs are embedded in MM1, where we found evidence of line broadening and temperature enhancements. On the other hand, the southern cloud and MM2 appear to be quiescent and starless. We modeled the radial intensity profile at 1.2mm of MM1. The model fits reasonably well the data, but it underestimates the intensity at small projected distances from the 1.2mm peak, probably due to the presence of multiple YSOs embedded in the envelope. There is a differentiation in the relative NH3 abundance with low values, ~2x10^-8, toward MM1, and high values, up to 10^-6, toward the southern cloud and MM2, suggesting that these clouds could be in a young evolutionary stage. IRAS 00213+6530 is harboring a multiple system of low-mass protostars, indicating that star formation in this cloud is taking place in groups, rather than in isolation. The low-mass YSOs found in IRAS 00213+6530 are in different evolutionary stages suggesting that star formation is taking place in different episodes.

HH 211 is a nearby young protostellar system with a highly collimated jet. We have mapped it in 352 GHz continuum, SiO (J=8-7), and HCO+ (J=4-3) emission at up to ~ 0.2" resolution with the Submillimeter Array (SMA). The continuum source is now resolved into two sources, SMM1 and SMM2, with a separation of ~ 84 AU. SMM1 is seen at the center of the jet, probably tracing a (inner) dusty disk around the protostar driving the jet. SMM2 is seen to the southwest of SMM1 and may trace an envelope-disk around a small binary companion. A flattened envelope-disk is seen in HCO+ around SMM1 with a radius of ~ 80 AU perpendicular to the jet axis. Its velocity structure is consistent with a rotation motion and can be fitted with a Keplerian law that yields a mass of ~ 50+-15 Jupiter mass (a mass of a brown dwarf) for the protostar. Thus, the protostar could be the lowest mass source known to have a collimated jet and a rotating flattened envelope-disk. A small-scale (~ 200 AU) low-speed (~ 2 km/s) outflow is seen in HCOP+ around the jet axis extending from the envelope-disk. It seems to rotate in the same direction as the envelope-disk and may carry away part of the angular momentum from the envelope-disk. The jet is seen in SiO close to ~ 100 AU from SMM1. It is seen with a "C-shaped" bending. It has a transverse width of <~ 40 AU and a velocity of ~ 170+-60 km/s. A possible velocity gradient is seen consistently across its innermost pair of knots, with ~ 0.5 km/s at ~ 10 AU, consistent with the sense of rotation of the envelope-disk. If this gradient is an upper limit of the true rotational gradient of the jet, then the jet carries away a very small amount of angular momentum of ~ 5 AU km/s and thus must be launched from the very inner edge of the disk near the corotation radius.

Sensitive and high angular resolution ($\sim$ 0.4\arcsec) SO$_2$[22$_{2,20}$ $\to$ 22$_{1,21}$] and SiO[5$\to$4] line and 1.3 and 7 mm continuum observations made with the Submillimeter Array (SMA) and the Very Large Array (VLA) towards the young massive cluster W51 IRS2 are presented. We report the presence of a large (of about 3000 AU) and massive (40 M$_\odot$) dusty circumstellar disk and a hot gas molecular ring around a high-mass protostar or a compact small stellar system associated with W51 North. The simultaneous observations of the silicon monoxide molecule, an outflow gas tracer, further revealed a massive (200 M$_\odot$) and collimated ($\sim14^\circ$) outflow nearly perpendicular to the dusty and molecular structures suggesting thus the presence of a single very massive protostar with a bolometric luminosity of more than 10$^5$ L$_\odot$. A molecular hybrid LTE model of a Keplerian and infalling ring with an inner cavity and a central stellar mass of more than 60 M$_\odot$ agrees well with the SO$_2$[22$_{2,20}$ $\to$ 22$_{1,21}$] line observations. Finally, these results suggest that mechanisms, such as mergers of low- and intermediate- mass stars, might be not necessary for forming very massive stars.

We present the results of millimeter and centimeter continuum observations, made with the IRAM 30m telescope and the VLA, toward a sample of 11 luminous IRAS sources classified as high-mass protostellar object candidates. We find 1.2 mm emission for all (but one) regions likely tracing the dust core in which the massive young stellar object is forming, for which we estimate masses ranging from 10 to 140 Msun. For all the sources, but one, we detect centimeter emission associated with the IRAS source, being compact or ultracompact HII region candidates, with early B-type stars as ionizing stars. The 7 mm emission is partially resolved for the four sources observed at this wavelength, with contribution of dust emission at 7 mm ranging from negligible to 44%. By combining our data with infrared surveys we fitted the spectral energy distribution of the sources. Finally, we find a correlation between the degree of disruption of the natal cloud, estimated from the fraction of dust emission associtaed with the centimeter source relative to the total amount of dust in its surroundings, and the size of the centimeter source. From this correlation, we establish an evolutionary sequence which is consistent with the evolutionary stage expected from maser/outflow/dense gas emission and with the infrared excess.

We present high angular resolution ($\sim$ 3$''$) and sensitive 1.3 mm continuum, cyanogen (CN) and vinyl cyanide (C$_2$H$_3$CN) line observations made with the Submillimeter Array (SMA) toward one of most highly obscured objects of the W51 IRS2 region, W51 North. We find that the CN line exhibits a pronounced inverse P-Cygni profile indicating that the molecular gas is infalling inwards this object with a mass accretion rate between 4 and 7 $\times$ 10$^{-2}$ M$_\odot$ yr$^{-1}$. The C$_2$H$_3$CN traces an east-west rotating molecular envelope that surrounds either a single obscured (proto)star with a kinematic mass of 40 M$_{\odot}$ or a small central cluster of B-type stars and that is associated with a compact high velocity bipolar outflow traced by H$_2$O masers and SiO molecular emission. We thus confirm that the W51 North region is part of the growing list of young massive star forming regions that have been associated with infalling motions and with large mass accretion rates ($\sim$ 10$^{-2}$ -- 10$^{-4}$), strengthening the evidence for massive stars forming with very high accretion rates sufficient to quench the formation of an UCHII region.

In order to elucidate the nature of the brightest infrared sources associated with IRAS 20343+4129, IRS1 and IRS3, we observed with the Submillimeter Array (SMA) the 1.3 mm continuum and CO(2-1) emission of the region. Faint millimeter dust continuum emission was detected toward IRS1, and we derived an associated gas mass of ~0.8 Msun. The IRS1 spectral energy distribution agrees with IRS1 being an intermediate-mass Class I source of about 1000 Lsun, whose circumstellar material is producing the observed large infrared excess. We have discovered a high-velocity CO bipolar outflow in the east-west direction, which is clearly associated with IRS1, and the outflow parameters are similar to those of intermediate-mass young stellar objects. Associated with the blue large scale CO outflow lobe, detected with single-dish observations, we only found two elongated low-velocity structures on either side of IRS3. The large-scale outflow lobe is almost completely resolved out by the SMA. Our detected low-velocity CO structures are coincident with elongated H2 emission features. The strongest millimeter continuum condensations in the region are found on either side of IRS3, where the infrared emission is extremely weak, and the CO and H2 elongated structures follow the border of the millimeter continuum emission that is facing IRS3. All these results suggest that the dust is associated with the walls of an expanding cavity driven by IRS3, estimated to be a B2 star. Within and beyond the expanding cavity, the millimeter continuum sources can be sites of future low-mass star formation.

We have mapped the protostellar jet HH 211 in 342 GHz continuum, SiO ($J=8-7$), and CO ($J=3-2$) emission at $\sim$ \arcs{1} resolution with the Submillimeter Array (SMA). Thermal dust emission is seen in continuum at the center of the jet, tracing an envelope and a possible optically thick compact disk (with a size $<$ 130 AU) around the protostar. A knotty jet is seen in CO and SiO as in \H2{}, but extending closer to the protostar. It consists of a chain of knots on each side of the protostar, with an interknot spacing of $\sim$ \arcs{2}$-$\arcs{3} or 600$-$900 AU and the innermost pair of knots at only $\sim$ \arcsa{1}{7} or 535 AU from the protostar. These knots likely trace unresolved internal (bow) shocks (i.e., working surfaces) in the jet, with a velocity range up to $\sim$ 25 \vkm{}. The two-sided mass-loss rate of the jet is estimated to be $\sim (0.7-2.8)\times 10^{-6}$ \solarmass{} yr$^{-1}$. The jet is episodic, precessing, and bending. A velocity gradient is seen consistently across two bright SiO knots (BK3 and RK2) perpendicular to the jet axis, with $\sim$ 1.5$\pm$0.8 \vkm{} at $\sim$ 30$\pm$15 AU, suggesting a presence of a jet rotation. The launching radius of the jet, derived from the potential jet rotation, is $\sim$ 0.15$-$0.06 AU in the inner disk.

We aim at studying the cluster environment surrounding the UCHII region in IRAS 20293+3952, a region in the first stages of formation of a cluster around a high-mass star. BIMA and VLA were used to observe the 3 mm continuum, N2H+(1-0), NH3(1,1) and (2,2), and CH3OH(2-1) emission of the surroundings of the UCHII region. We studied the kinematics of the region and computed the rotational temperature and column density maps by fitting the hyperfine structure of N2H+ and NH3. The dense gas traced by N2H+ and NH3 shows two different clouds, a main cloud to the east of the UCHII region, of 0.5 pc and 250 Msun, and a western cloud, of 0.15 pc and 30 Msun. The dust emission reveals two strong components in the northern side of the main cloud, associated with Young Stellar Objects (YSOs) driving molecular outflows, and two fainter components in the southern side, with no signs of star forming activity. Regarding the CH3OH, we found strong emission in a fork-like structure associated with outflow B, as well as emission associated with outflow A. The rotational temperature is higher in the northern side of the main cloud, around 22 K, where there are most of the YSOs, than in the southern side, around 16 K. There is strong chemical differentiation in the region, since we determined low values of the NH3/N2H+ ratio, around 50, associated with YSOs in the north of the main cloud, and high values, up to 300, associated with cores with no detected YSOs, in the south of the main cloud. Such a chemical differentiation is likely due to abundance/depletion effects. Finally, interaction between the different sources in the region is important, mainly through their UV radiation and outflows (abridged).

We observed the HH 211 jet in the submillimeter continuum and the CO(3-2) and SiO(8-7) transitions with the Submillimeter Array. The continuum source detected at the center of the outflow shows an elongated morphology, perpendicular to the direction of the outflow axis. The high-velocity emission of both molecules shows a knotty and highly collimated structure. The SiO(8-7) emission at the base of the outflow, close to the driving source, spans a wide range of velocities, from -20 up to 40 km s^{-1}. This suggests that a wide-angle wind may be the driving mechanism of the HH 211 outflow. For distances greater than 5" (1500 AU) from the driving source, emission from both transitions follows a Hubble-law behavior, with SiO(8-7) reaching higher velocities than CO(3-2), and being located upstream of the CO(3-2) knots. This indicates that the SiO(8-7) emission is likely tracing entrained gas very close to the primary jet, while the CO(3-2) is tracing less dense entrained gas. From the SiO(5-4) data of Hirano et al. we find that the SiO(8-7)/SiO(5-4) brightness temperature ratio along the jet decreases for knots far from the driving source. This is consistent with the density decreasing along the jet, from (3-10)x10^6 cm^{-3} at 500 AU to (0.8-4)x10^6 cm^{-3} at 5000 AU from the driving source.

We present centimeter and millimeter observations of gas and dust around IRAS 21391+5802, an intermediate-mass source embedded in the core of IC 1396N. Continuum observations from 3.6 cm to 1.2 mm are used to study the embedded objects and overall distribution of the dust, while molecular line observations of CO, CS, and CH3OH are used to probe the structure and chemistry of the outflows in the region. The continuum emission at centimeter and millimeter wavelengths has been resolved into three sources separated about 15 arcsec from each other, and with one of them, BIMA 2, associated with IRAS 21391+5802. The dust emission around this source shows a very extended envelope, which accounts for most of the circumstellar mass of 5.1 Msun. This source is powering a strong molecular outflow, elongated in the E--W direction, which presents a complex structure and kinematics. While at high outflow velocities the outflow is clearly bipolar, at low outflow velocities the blueshifted and redshifted emission are highly overlapping, and the strongest emission shows a V-shaped morphology. The outflow as traced by CS and CH3OH exhibits two well differentiated and clumpy lobes, with two prominent northern blueshifted and redshifted clumps. The curved shape of the clumps and the spectral shape at these positions are consistent with shocked material. In addition, CS and CH3OH are strongly enhanced toward these positions with respect to typical quiescent material abundances in other star-forming regions.