Io's atmosphere is predominately SO2 sustained by a combination of volcanic outgassing and sublimation. The loss from the atmosphere is the main mass source for Jupiter's large magnetosphere. Previous studies attributed various transient phenomena in Io's environment and Jupiter's magnetosphere to a sudden change in the mass loss from the atmosphere supposedly triggered by a change in volcanic activity. Since the gas in volcanic plumes does not escape directly, such causal correlation would require a transient volcano-induced change in atmospheric abundance, which has never been observed so far. Here we report four observations of atmospheric SO2 and NaCl obtained with the IRAM NOEMA interferometer. These observations are compared to measurements of volcanic hot spots and Io's neutral and plasma environment. We find a stable NaCl column density in Io's atmosphere. The SO2 column density derived for December 2016 is about 30% lower compared to the period of March to April 2017. This increase in SO2 from December 2016 to March 2017 might be related to increasing volcanic activity observed at several sites in spring 2017, but the stability of the volcanic trace gas NaCl and resulting decrease in NaCl/SO2 ratio do not support this interpretation. Observed dimmings in both the sulfur ion torus and Na neutral cloud suggest rather a decrease in mass loading in the period of increasing SO2 abundance. The dimming Na brightness and stable atmospheric NaCl furthermore dispute an earlier suggested positive correlation of the sodium cloud and the hot spot activity at Loki Patara, which considerably increased in this period. The environment of Io overall appears to be in a quiescent state, preventing further conclusions. Only Jupiter's aurora morphology underwent several short-term changes, which are apparently unrelated to Io's quiescent environment or the relatively stable atmosphere.

This paper was withdrawed from the ApJ after the comments from the referee, please be carefully.

High-mass stars and star clusters commonly form within hub-filament systems. Monoceros R2, harbors one of the closest such systems, making it an excellent target for case studies. We investigate the morphology, stability and dynamical properties of the hub-filament system on basis of 13CO and C18O observations obtained with the IRAM-30m telescope and H2 column density maps derived from Herschel dust emission observations. We identified the filamentary network and characterized the individual filaments as either main (converging into the hub) or secondary (converging to a main filament) filaments. The main filaments have line masses of 30-100 Msun/pc and show signs of fragmentation. The secondary filaments have line masses of 12-60 Msun/pc and show fragmentation only sporadically. In the context of Ostriker's hydrostatic filament model, the main filaments are thermally super-critical. If non-thermal motions are included, most of them are trans-critical. Most of the secondary filaments are roughly trans-critical regardless of whether non-thermal motions are included or not. From the main filaments, we estimate a mass accretion rate of 10(-4)-10(-3) Msun/pc into the hub. The secondary filaments accrete into the main filaments with a rate of 0.1-0.4x10(-4) Msun/pc. The main filaments extend into the hub. Their velocity gradients increase towards the hub, suggesting acceleration of the gas. We estimate that with the observed infall velocity, the mass-doubling time of the hub is ~2.5 Myr, ten times larger than the free-fall time, suggesting a dynamically old region. These timescales are comparable with the chemical age of the HII region. Inside the hub, the main filaments show a ring- or a spiral-like morphology that exhibits rotation and infall motions. One possible explanation for the morphology is that gas is falling into the central cluster following a spiral-like pattern.

We present the highest angular resolution (20x15mas - 44x33au) Atacama Large Millimeter/sub-millimeter Array (ALMA) observations currently possible of the proto-O-star G17.64+0.16 in Band 6. The Cycle 5 observations with baselines out to 16km probes scales <50au and reveal the rotating disc around G17.64+0.16, a massive forming O-type star. The disc has a ring-like enhancement in the dust emission, especially visible as arc structures to the north and south. The Keplerian kinematics are most prominently seen in the vibrationally excited water line, H2O (Eu=3461.9K). The mass of the central source found by modelling the Keplerian rotation is consistent with 45+/-10Mo. The H30alpha (231.9GHz) radio-recombination line and the SiO (5-4) molecular line were detected at up to the 10 sigma$ level. The estimated disc mass is 0.6-2.6Mo under the optically thin assumption. Analysis of the Toomre Q parameter, in the optically thin regime, indicates that the disc stability is highly dependent on temperature. The disc currently appears stable for temperatures >150K, this does not preclude that the substructures formed earlier through disc fragmentation.

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.

Understanding the matter cycle in the interstellar medium of galaxies from the assembly of clouds to star formation and stellar feedback remains an important and exciting field in comtemporary astrophysics. Many open questions regarding cloud and structure formation, the role of turbulence, and the relative importance of the various feedback processes can only be addressed with observations of spectrally resolved lines. We here stress the importance of two specific sets of lines: the finestructure lines of atomic carbon as a tracer of the dark molecular gas and mid-J CO lines as tracers of the warm, active molecular gas in regions of turbulence dissipation and feedback. The observations must cover a wide range of environments (i.e., physical conditions), which will be achieved by large scale surveys of Galactic molecular clouds, the Galactic Center, the Magellanic clouds, and nearby galaxies. To date, such surveys are completely missing and thus constitute an important science opportunity for the next decade and beyond. For the successful interpretation of the observations, it will be essential to combine them with results from (chemical) modelling and simulations of the interstellar medium.

We present an initial overview of the filamentary structure in the Orion A molecular cloud utilizing a high angular and velocity resolution C$^{18}$O(1-0) emission map that was recently produced as part of the CARMA-NRO Orion Survey. The main goal of this study is to build a credible method to study varying widths of filaments which has previously been linked to star formation in molecular clouds. Due to the diverse star forming activities taking place throughout its $\sim$20 pc length, together with its proximity of 388 pc, the Orion A molecular cloud provides an excellent laboratory for such an experiment to be carried out with high resolution and high sensitivity. Using the widely-known structure identification algorithm, DisPerSE, on a 3-dimensional (PPV) C$^{18}$O cube, we identified 625 relatively short (the longest being 1.74 pc) filaments over the entire cloud. We study the distribution of filament widths using FilChaP, a python package that we have developed and made publicly available. We find that the filaments identified in a 2 square degree PPV cube do not overlap spatially, except for the complex OMC-4 region that shows distinct velocity components along the line of sight. The filament widths vary between 0.02 and 0.3 pc depending on the amount of substructure that a filament possesses. The more substructure a filament has, the larger is its width. We also find that despite this variation, the filament width shows no anticorrelation with the central column density which is in agreement with previous Herschel observations.

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.

We model the emission of methyl cyanide (CH3CN) lines towards the massive hot molecular core SgrB2(M). We aim at reconstructing the CH3CN abundance field and investigating the gas temperature distribution as well as the velocity field. SgrB2(M) was observed with the ALMA in a spectral line survey from 211 to 275 GHz. This frequency range includes several transitions of CH3CN (including isotopologues and vibrationally excited states). We employ the three-dimensional radiative transfer toolbox Pandora in order to retrieve the velocity and abundance field by modeling different CH3CN lines. For this purpose, we base our model on the results of a previous study that determined the physical structure of SgrB2(M), i.e.\ the distribution of dust dense cores, ionized regions and heating sources. The morphology of the CH3CN emission can be reproduced by a molecular density field that consists of a superposition of cores with modified Plummer-like density profiles. The averaged relative abundance of CH3CN with respect to H2 ranges from 4x10^{-11} to 2x10^{-8} in the northern part of SgrB2(M) and from 2x10^{-10} to 5x10^{-7} in the southern part. In general, we find that the relative abundance of CH3CN is lower at the center of the very dense and hot cores, causing the general morphology of the CH3CN emission to be shifted with respect to the dust continuum emission. The dust temperature calculated by the radiative transfer simulation based on the available luminosity reaches values up to 900 K. However, in some regions vibrationally excited transitions of CH3CN are underestimated by the model, indicating that the predicted gas temperature, which is assumed to be equal to the dust temperature, is partly underestimated. The determination of the velocity component along the line of sight reveals that a velocity gradient from the north to the south exists in SgrB2(M).

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.

We present a multiwavelength study of the filamentary infrared dark cloud (IRDC) G333.73+0.37. The region contains two distinct mid-infrared sources S1 and S2 connected by dark lanes of gas and dust. Cold dust emission from the IRDC is detected at seven wavelength bands and we have identified 10 high density clumps in the region. The physical properties of the clumps such as temperature: 14.3-22.3 K and mass: 87-1530 M_sun are determined by fitting a modified blackbody to the spectral energy distribution of each clump between 160 micron and 1.2 mm. The total mass of the IRDC is estimated to be $~4700 M_sun. The molecular line emission towards S1 reveals signatures of protostellar activity. Low frequency radio emission at 1300 and 610 MHz is detected towards S1 (shell-like) and S2 (compact morphology), confirming the presence of newly formed massive stars in the IRDC. Photometric analysis of near and mid-infrared point sources unveil the young stellar object population associated with the cloud. Fragmentation analysis indicates that the filament is supercritical. We observe a velocity gradient along the filament, that is likely to be associated with accretion flows within the filament rather than rotation. Based on various age estimates obtained for objects in different evolutionary stages, we attempt to set a limit to the current age of this cloud.

STATCONT is a python-based tool designed to determine the continuum emission level in spectral data, in particular for sources with a line-rich spectrum. The tool inspects the intensity distribution of a given spectrum and automatically determines the continuum level by using different statistical approaches. The different methods included in STATCONT are tested against synthetic data. We conclude that the sigma-clipping algorithm provides the most accurate continuum level determination, together with information on the uncertainty in its determination. This uncertainty can be used to correct the final continuum emission level, resulting in the here called `corrected sigma-clipping method' or c-SCM. The c-SCM has been tested against more than 750 different synthetic spectra reproducing typical conditions found towards astronomical sources. The continuum level is determined with a discrepancy of less than 1% in 50% of the cases, and less than 5% in 90% of the cases, provided at least 10% of the channels are line free. The main products of STATCONT are the continuum emission level, together with a conservative value of its uncertainty, and datacubes containing only spectral line emission, i.e., continuum-subtracted datacubes. STATCONT also includes the option to estimate the spectral index, when different files covering different frequency ranges are provided.

The high-mass star forming sites SgrB2(M) and SgrB2(N) have been the target of numerous studies, revealing e.g. a rich chemistry. We want to characterize their physical and chemical structure using ALMA high-angular resolution observations at mm wavelengths, reaching spatial scales of about 4000 au, and covering the whole band 6 (from 211 to 275 GHz). In order to determine the continuum emission in line-rich sources, we use a new statistical method: STATCONT. We detect 27 continuum sources in SgrB2(M) and 20 in SgrB2(N). We study the continuum emission across the ALMA band 6, and compare it with previous SMA 345 GHz and VLA 40 GHz observations, to study the nature of the sources detected. The brightest sources are dominated by (partially optically thick) dust emission, while there is an important degree of contamination from ionized gas free-free emission in weaker sources. While the total mass in SgrB2(M) is distributed in many fragments, most of the mass in SgrB2(N) arises from a single object, with filamentary-like structures converging towards the center. There seems to be a lack of low-mass dense cores in both regions. We determine H2 volume densities for the cores of about 10^5-10^7 Msun pc^-3, one to two orders of magnitude higher than the stellar densities of super star clusters. In general, SgrB2(N) is chemically richer than SgrB2(M). There seems to be a correlation between the chemical richness and the mass of the fragments, with more massive clumps being more chemically rich. Both SgrB2(N) and SgrB2(M) harbour a cluster of hot molecular cores. We compare the continuum images with predictions from a detailed 3D radiative transfer model that reproduces the structure of SgrB2 from 45 pc down to 100 au. This dataset, together with ongoing projects in the range 5 to 200 GHz, better constrain the 3D structure of SgrB2, and allow us to understand its physical and chemical structure.

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.

The CO+ reactive ion is thought to be a tracer of the boundary between a HII region and the hot molecular gas. In this study, we present the spatial distribution of the CO+ rotational emission toward the Mon R2 star-forming region. The CO+ emission presents a clumpy ring-like morphology, arising from a narrow dense layer around the HII region. We compare the CO+ distribution with other species present in photon-dominated regions (PDR), such as [CII] 158 mm, H2 S(3) rotational line at 9.3 mm, polycyclic aromatic hydrocarbons (PAHs) and HCO+. We find that the CO+ emission is spatially coincident with the PAHs and [CII] emission. This confirms that the CO+ emission arises from a narrow dense layer of the HI/H2 interface. We have determined the CO+ fractional abundance, relative to C+ toward three positions. The abundances range from 0.1 to 1.9x10^(-10) and are in good agreement with previous chemical model, which predicts that the production of CO+ in PDRs only occurs in dense regions with high UV fields. The CO+ linewidth is larger than those found in molecular gas tracers, and their central velocity are blue-shifted with respect to the molecular gas velocity. We interpret this as a hint that the CO+ is probing photo-evaporating clump surfaces.

Almost 200 different species have been detected in the interstellar medium (ISM) during the last decades, revealing not only simple species but complex molecules with more than 6 atoms. Other exotic compounds, like the weakly-bound dimer (H2)2, have also been detected in astronomical sources like Jupiter. We aim at detecting for the first time the CO-H2 van der Waals complex in the ISM, which if detected can be a sensitive indicator for low temperatures. We use the IRAM30m telescope, located in Pico Veleta (Spain), to search for the CO-H2 complex in a cold, dense core in TMC-1C (with a temperature of 10 K). All the brightest CO-H2 transitions in the 3 mm (80-110 GHz) band have been observed with a spectral resolution of 0.5-0.7 km/s, reaching a rms noise level of 2 mK. The simultaneous observation of a broad frequency band, 16 GHz, has allowed us to conduct a serendipitous spectral line survey. No lines belonging to the CO-H2 complex have been detected. We have set up a new, more stringent upper limit for its abundance to be [CO-H2]/[CO] = 5x10^{-6}, while we expect the abundance of the complex to be in the range 10^{-8}-10^{-3}. The spectral line survey has allowed us to detect 75 lines associated with 41 different species (including isotopologues). We detect a number of complex organic species, e.g. methyl cyanide (CH3CN), methanol (CH3OH), propyne (CH3CCH) and ketene (CH2CO), associated with cold gas (excitation temperatures about 7 K), confirming the presence of these complex species not only in warm objects but also in cold regimes.

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.

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 describe the assignment of a previously unidentified interstellar absorption line to ArH$^+$ and discuss its relevance in the context of hydride absorption in diffuse gas with a low H$_2$ fraction. The column densities along several lines of sight are determined and discussd in the framework of chemical models. The column densities of ArH$^+$ are compared to those of other species, tracing interstellar medium (ISM) components with different H$_2$ abundances. Chemical models are constructed, taking UV radiation and cosmic ray ionization into account. Due to the detection of two isotopologues, $^{36}$ArH$^+$ and $^{38}$ArH$^+$, we are confident about the carrier assignment to ArH$^+$. NeH$^+$ is not detected with a limit of [NeH$^+$]/[ArH$^+$] $\le$ 0.1. The derived column densities agree well with the predictions of chemical models. ArH$^+$ is a unique tracer of gas with a fractional H$_2$ abundance of $10^{-4}- 10^{-3}$ and shows little correlation with H$_2$O$^+$, which traces gas with a fractional H$_2$ abundance of $\approx $0.1. A careful analysis of variations in the ArH$^+$, OH$^+$, H$_2$O$^+$ and HF column densities promises to be a faithful tracer of the distribution of the H$_2$ fractional abundance, providing unique information on a poorly known phase in the cycle of interstellar matter, its transition from atomic diffuse gas to dense molecular gas traced by CO emission. Abundances of these species put strong observational constraints upon magnetohydrodynamical (MHD) simulations of the interstellar medium, and potentially could evolve into a tool to characterize the ISM. Paradoxically, the ArH$^+$ molecule is a better tracer of \new{almost} purely atomic hydrogen gas than H{\sc i} itself, since H{\sc i} can also be present in gas with a significant molecular content, but ArH$^+$ singles out gas that is $>99.9$\% atomic.

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.

(Abridged) Aims. To investigate the first stages of the process of high-mass star formation, we selected a sample of massive clumps previously observed with the SEST at 1.2 mm and with the ATNF ATCA at 1.3 cm. We want to characterize the physical conditions in such sources, and test whether their properties depend on the evolutionary stage of the clump. Methods. With ATCA we observed the selected sources in the NH3(1,1) and (2,2) transitions and in the 22 GHz H2O maser line. Ammonia lines are a good temperature probe that allow us to accurately determine the mass and the column-, volume-, and surface densities of the clumps. We also collected all data available to construct the spectral energy distribution of the individual clumps and to determine if star formation is already occurring, through observations of its most common signposts, thus putting constraints on the evolutionary stage of the source. We fitted the spectral energy distribution between 1.2 mm and 70 microns with a modified black body to derive the dust temperature and independently determine the mass. Results. The clumps are cold (T~10-30 K), massive (M~10^2-10^3 Mo), and dense (n(H2)>~10^5 cm^-3) and they have high column densities (N(H2)~10^23 cm^-2). All clumps appear to be potentially able to form high-mass stars. The most massive clumps appear to be gravitationally unstable, if the only sources of support against collapse are turbulence and thermal pressure, which possibly indicates that the magnetic field is important in stabilizing them. Conclusions. After investigating how the average properties depend on the evolutionary phase of the source, we find that the temperature and central density progressively increase with time. Sources likely hosting a ZAMS star show a steeper radial dependence of the volume density and tend to be more compact than starless clumps.

Theoretical models suggest that massive stars form via disk-mediated accretion, with bipolar outflows playing a fundamental role. A recent study toward massive molecular outflows has revealed a decrease of the SiO line intensity as the object evolves. The present study aims at characterizing the variation of the molecular outflow properties with time, and at studying the SiO excitation conditions in outflows associated with massive YSOs. We used the IRAM30m telescope to map 14 massive star-forming regions in the SiO(2-1), SiO(5-4) and HCO+(1-0) outflow lines, and in several dense gas and hot core tracers. Hi-GAL data was used to improve the spectral energy distributions and the L/M ratio, which is believed to be a good indicator of the evolutionary stage of the YSO. We detect SiO and HCO+ outflow emission in all the sources, and bipolar structures in six of them. The outflow parameters are similar to those found toward other massive YSOs. We find an increase of the HCO+ outflow energetics as the object evolve, and a decrease of the SiO abundance with time, from 10^(-8) to 10^(-9). The SiO(5-4) to (2-1) line ratio is found to be low at the ambient gas velocity, and increases as we move to high velocities, indicating that the excitation conditions of the SiO change with the velocity of the gas (with larger densities and/or temperatures for the high-velocity gas component). The properties of the SiO and HCO+ outflow emission suggest a scenario in which SiO is largely enhanced in the first evolutionary stages, probably due to strong shocks produced by the protostellar jet. As the object evolves, the power of the jet would decrease and so does the SiO abundance. During this process, however, the material surrounding the protostar would have been been swept up by the jet, and the outflow activity, traced by entrained molecular material (HCO+), would increase with time.

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 ALMA observations of continuum and molecular line emission with 0.4" resolution towards the high-mass star forming region G35.20-0.74 N. Two dense cores are detected in typical hot-core tracers, such as CH3CN, which reveal velocity gradients. In one of these cores, the velocity field can be fitted with an almost edge-on Keplerian disk rotating about a central mass of 18 Msun. This finding is consistent with the results of a recent study of the CO first overtone bandhead emission at 2.3mum towards G35.20-0.74 N. The disk radius and mass are >2500 au and 3 Msun. To reconcile the observed bolometric luminosity (3x10^4 Lsun) with the estimated stellar mass of 18 Msun, we propose that the latter is the total mass of a binary system.

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 ATCA observations of the H2O maser line and radio continuum at 18.0GHz and 22.8GHz, toward a sample of 192 massive star forming regions containing several clumps already imaged at 1.2mm. The main aim of this study is to investigate the water maser and centimeter continuum emission (likely tracing thermal free-free emission) in sources at different evolutionary stages, using the evolutionary classifications proposed by Palla et al (1991) and Molinari et al (2008). We used the recently comissioned CABB backend at ATCA obtaining images with 20arcsec resolution in the 1.3cm continuum and H2O maser emission, in all targets. For the evolutionary analysis of the sources we used the millimeter continuum emission from Beltran et al (2006) and the infrared emission from the MSX Point Source Catalogue. We detect centimeter continuum emission in 88% of the observed fields with a typical rms noise level of 0.45mJy/beam. Most of the fields show a single radio continuum source, while in 20% of them we identify multiple components. A total of 214 centimeter continuum sources have been identified, likely tracing optically thin HII regions, with physical parameters typical of both extended and compact HII regions. Water maser emission was detected in 41% of the regions, resulting in a total of 85 distinct components. The low angular (20arcsec) and spectral (14km/s) resolutions do not allow a proper analysis of the water maser emission, but suffice to investigate its association with the continuum sources. We have also studied the detection rate of HII regions in the two types of IRAS sources defined by Palla et (1991) on the basis of the IRAS colours: High and Low. No significant differences are found, with large detection rates (>90%) for both High and Low sources. We classify the millimeter and infrared sources in our fields in three evolutionary stages following the scheme presented by ...

The Pipe nebula is a massive, nearby, filamentary dark molecular cloud with a low star-formation efficiency threaded by a uniform magnetic field perpendicular to its main axis. It harbors more than a hundred, mostly quiescent, very chemically young starless cores. The cloud is, therefore, a good laboratory to study the earliest stages of the star-formation process. We aim to investigate the primordial conditions and the relation among physical, chemical, and magnetic properties in the evolution of low-mass starless cores. We used the IRAM 30-m telescope to map the 1.2 mm dust continuum emission of five new starless cores, which are in good agreement with previous visual extinction maps. For the sample of nine cores, which includes the four cores studied in a previous work, we derived a Av to NH2 factor of (1.27$\pm$0.12)$\times10^{-21}$ mag cm$^{2}$ and a background visual extinction of ~6.7 mag possibly arising from the cloud material. We derived an average core diameter of ~0.08 pc, density of ~10$^5$ cm$^{-3}$, and mass of ~1.7 Msun. Several trends seem to exist related to increasing core density: (i) diameter seems to shrink, (ii) mass seems to increase, and (iii) chemistry tends to be richer. No correlation is found between the direction of the surrounding diffuse medium magnetic field and the projected orientation of the cores, suggesting that large scale magnetic fields seem to play a secondary role in shaping the cores. The full abstract is available in the pdf.

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.

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.

Hypercompact (HC) HII regions are, by nature, very young HII regions, associated with the earliest stages of massive star formation. They may represent the transition phase as an early B-type star grows into an O-type star. Unfortunately, so few HCHII regions are presently known that their general attributes and defining characteristics are based on small number statistics. A larger sample is needed for detailed studies and good statistics. Class II methanol masers are one of the best indicators of the early stages of massive star formation. Using the Arecibo Methanol Maser Galactic Plane Survey - the most sensitive blind survey for 6.7 GHz methanol masers to date - we selected 24 HCHII region candidates. We made EVLA continuum observations at 3.6 and 1.3 cm to search for HCHII regions associated with these masers. We identified six potential HCHII regions in our sample based on the presence of optically thick free-free emission. Overall, we find that 30% of the methanol masers have an associated centimeter radio continuum source (separation less than 0.1 pc), which is in general agreement with previous studies.

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.

The Pipe Nebula is a massive, nearby dark molecular cloud with a low star-formation efficiency which makes it a good laboratory to study the very early stages of the star formation process. The Pipe Nebula is largely filamentary, and appears to be threaded by a uniform magnetic field at scales of few parsecs, perpendicular to its main axis. The field is only locally perturbed in a few regions, such as the only active cluster forming core B59. The aim of this study is to investigate primordial conditions in low-mass pre-stellar cores and how they relate to the local magnetic field in the cloud. We used the IRAM 30-m telescope to carry out a continuum and molecular survey at 3 and 1 mm of early- and late-time molecules toward four selected starless cores inside the Pipe Nebula. We found that the dust continuum emission maps trace better the densest regions than previous 2MASS extinction maps, while 2MASS extinction maps trace better the diffuse gas. The properties of the cores derived from dust emission show average radii of ~0.09 pc, densities of ~1.3x10^5 cm^-3, and core masses of ~2.5 M_sun. Our results confirm that the Pipe Nebula starless cores studied are in a very early evolutionary stage, and present a very young chemistry with different properties that allow us to propose an evolutionary sequence. All of the cores present early-time molecular emission, with CS detections toward all the sample. Two of them, Cores 40 and 109, present strong late-time molecular emission. There seems to be a correlation between the chemical evolutionary stage of the cores and the local magnetic properties that suggests that the evolution of the cores is ruled by a local competition between the magnetic energy and other mechanisms, such as turbulence.

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.

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.