Theories of massive star formation predict that massive protostars accrete gas through circumstellar disks. Although several cases have been found already thanks to high-angular resolution interferometry, it remains unknown the internal physical structure of these disks and, in particular, whether they present warps or internal holes as observed in low-mass proto-planetary disks. Here, we report very high angular resolution observations of the H21alpha radio recombination line carried out in Band 9 with the Atacama Large Millimeter/submillimeter Array (beam of 80 mas x 60 mas, or 70 au x 50 au) toward the IRS2 massive young stellar object in the Monoceros R2 star-forming cluster. The H21alpha line shows maser amplification, which allows us to study the kinematics and physical structure of the ionised gas around the massive protostar down to spatial scales of ~1-2 au. Our ALMA images and 3D radiative transfer modelling reveal that the ionized gas around IRS2 is distributed in a Keplerian circumstellar disk and an expanding wind. The H21alpha emission centroids at velocities between -10 and 20 km s-1 deviate from the disk plane, suggesting a warping for the disk. This could be explained by the presence of a secondary object (a stellar companion or a massive planet) within the system. The ionized wind seems to be launched from the disk surface at distances ~11 au from the central star, consistent with magnetically-regulated disk wind models. This suggests a similar wind launching mechanism to that recently found for evolved massive stars such as MWC349A and MWC922.

We discuss the science drivers for ALMA Band 2 which spans the frequency range from 67 to 90 GHz. The key science in this frequency range are the study of the deuterated molecules in cold, dense, quiescent gas and the study of redshifted emission from galaxies in CO and other species. However, Band 2 has a range of other applications which are also presented. The science enabled by a single receiver system which would combine ALMA Bands 2 and 3 covering the frequency range 67 to 116 GHz, as well as the possible doubling of the IF bandwidth of ALMA to 16 GHz, are also considered.

In the past decade, Astrochemistry has witnessed an impressive increase in the number of detections of complex organic molecules. Some of these species are of prebiotic interest such as glycolaldehyde, the simplest sugar, or amino acetonitrile, a possible precursor of glycine. Recently, we have reported the detection of two new nitrogen-bearing complex organics, glycolonitrile and Z-cyanomethanimine, known to be intermediate species in the formation process of ribonucleotides within theories of a primordial ribonucleic acid (RNA)-world for the origin of life. In this paper, we present deep and high-sensitivity observations toward two of the most chemically rich sources in the Galaxy: a Giant Molecular Cloud in the center of the Milky Way (G+0.693-0.027) and a proto-Sun (IRAS16293-2422 B). Our aim is to explore whether the key precursors considered to drive the primordial RNA-world chemistry, are also found in space. Our high-sensitivity observations reveal that urea is present in G+0.693-0.027 with an abundance of about 5x10-11. This is the first detection of this prebiotic species outside a star-forming region. Urea remains undetected toward the proto-Sun IRAS16293-2422 B (upper limit to its abundance of less than 2x10-11). Other precursors of the RNA-world chemical scheme such as glycolaldehyde or cyanamide are abundant in space, but key prebiotic species such as 2- amino-oxazole, glyceraldehyde or dihydroxyacetone are not detected in either source. Future more sensitive observations targeting the brightest transitions of these species will be needed to disentangle whether these large prebiotic organics are certainly present in space.

The SOLIS (Seeds Of Life In Space) IRAM/NOEMA Large Program aims at studying a set of crucial complex organic molecules in a sample of sources, with well-known physical structure, covering the various phases of Solar-type star formation. One representative object of the transition from the prestellar core to the protostar phases has been observed toward the Very Low Luminosity Object (VeLLO) called L1521F. This type of source is important to study to make the link between prestellar cores and Class 0 sources and also to constrain the chemical evolution during the process of star formation. Two frequency windows (81.6-82.6 GHz and 96.65-97.65 GHz) were used to observe the emission from several complex organics toward the L1521F VeLLO. Only 2 transitions of methanol (A+, E2) have been detected in the narrow window centered at 96.7 GHz (with an upper limit on E1) in a very compact emission blob (~7'' corresponding to ~1000au) toward the NE of the L1521F protostar. The CS 2-1 transition is also detected within the WideX bandwidth. Consistently, with what has been found in prestellar cores, the methanol emission appears ~1000au away from the dust peak. The location of the methanol blob coincides with one of the filaments previously reported in the literature. The Tex of the gas inferred from methanol is (10$\pm$2) K, while the H2 gas density (estimated from the detected CS 2-1 emission and previous CS 5-4 ALMA obs.) is a factor >25 higher than the density in the surrounding environment (n(H2) >10$^{7}$ cm$^{-3}$). From its compactness, low excitation temperature and high gas density, we suggest that the methanol emission detected with NOEMA is either a cold and dense shock-induced blob, recently formed ($\leq$ few hundred years) by infalling gas or a cold and dense fragment that may have just been formed as a result of the intense gas dynamics found within the L1521F VeLLO system.

Aims: The physical structure of a shock wave may take a form unique to its shock type, implying that the chemistry of each shock type is unique as well. We aim to investigate the different chemistries of J-type and C-type shocks in order to identify unique molecular tracers of both shock types. We apply these diagnostics to the protostellar outflow L1157 to establish whether the B2 clump could host shocks exhibiting type-specific behaviour. Of particular interest is the L1157-B2 clump, which has been shown to exhibit bright emission in S-bearing species and HNCO. Methods: We simulate, using a parameterised approach, a planar, steady-state J-type shock wave using UCLCHEM. We compute a grid of models using both C-type and J-type shock models to determine the chemical abundance of shock-tracing species as a function of distance through the shock and apply it to the L1157 outflow. We focus on known shock-tracing molecules such as H2O, HCN, and CH3 OH. Results: We find that a range of molecules including H2O and HCN have unique behaviour specific to a J-type shock, but that such differences in behaviour are only evident at low $v_{s}$ and low $n_{H}$. We find that CH3OH is enhanced by shocks and is a reliable probe of the pre-shock gas density. However, we find no difference between its gas-phase abundance in C-type and J-type shocks. Finally, from our application to L1157, we find that the fractional abundances within the B2 region are consistent with both C-type and J-type shock emission.

The size and composition of dust grains are critical in setting the dynamical, physical and chemical evolution of the media in which they are present. Thanks to facilities such as ALMA and in the future the SKA, their thermal emission in the (sub)mm to cm has become a convenient way to trace grain properties. Our aim is to understand the influence of the grain composition and size distribution on the shape of their SED in dense ISM regions such as molecular clouds, prestellar cores, YSOs and protoplanetary discs. Starting from the optical constants defined in the THEMIS model for amorphous hydrogenated carbon and silicate grains in addition to water ice, we define 6 material mixtures representative of the expected dust composition in dense ISM regions. The optical properties of 0.01 micron to 10 cm grains are then calculated with effective medium and Mie theories. The corresponding SEDs are calculated for isolated clouds either externally heated by the ISRF alone or in addition to an internal source. The 3 main outcomes of this study are as follows. First, the dust mass absorption coefficient strongly depends on both its composition and size distribution potentially leading to errors in dust mass estimates by factors up to 3 and 20, respectively. Second, it appears almost impossible to retrieve the grain composition from the (sub)mm to cm SED shape alone as its spectral index for lambda > 3 mm does not depend on composition. Third, using the true dust opacity spectral index to estimate grain sizes may lead to erroneous findings as the observed spectral index can be highly modified by the dust temperature distribution along the line-of-sight, which depends on the specific heating source and on the geometry of the studied region. Based on the interpretation of only the spectral shape of (sub)mm to cm SEDs, the determination of the dust masses, compositions and sizes are highly uncertain.

The main carrier of sulfur in dense clouds, where it is depleted from the gas phase, remains a mystery. Shock waves in young molecular outflows disrupt the ice mantles and allow us to directly probe the material that is ejected into the gas phase. A comprehensive study of sulfur-bearing species towards L1157-B1, a shocked region along a protostellar outflow, has been carried out as part of the IRAM-30m large program ASAI. The dataset contains over 100 lines of CCS, H$_2$CS, OCS, SO, SO$_2$ and isotopologues. The results of these observations are presented, complementing previous studies of sulfur-bearing species in the region. The column densities and fractional abundances of these species are measured and together these species account for 10\% of the cosmic sulfur abundance in the region. The gas properties derived from the observations are also presented, demonstrating that sulfur bearing species trace a wide range of different gas conditions in the region.

In star forming regions in our own Galaxy, the 14N/15N ratio is found to vary from $\sim$ 100 in meteorites, comets and protoplanetary disks up to $\sim$ 1000 in pre-stellar and star forming cores, while in external galaxies the very few single-dish large scale measurements of this ratio lead to values of 100-450. The extent of the contribution of isotopic fractionation to these variations is, to date, unknown. In this paper we present a theoretical chemical study of nitrogen fractionation in external galaxies in order to determine the physical conditions that may lead to a spread of the 14N/15N ratio from the solar value of $\sim$440 and hence evaluate the contribution of chemical reactions in the ISM to nitrogen fractionation. We find that the main cause of ISM enrichment of nitrogen fractionation is high gas densities, aided by high fluxes of cosmic rays.

Iram 30-m Observations towards eight protostellar outflow sources were taken in the 96-\SI{176}{\giga\hertz} range. Transitions of CH$_3$OH and CH$_3$CHO were detected in seven of them. The integrated emission of the transitions of each species that fell into the observed frequency range were measured and fit using RADEX and LTE models. Column densities and gas properties inferred from this fitting are presented. The ratio of the A and E-type isomers of CH$_3$OH indicate that the methanol observed in these outflows was formed on the grain surface. Both species demonstrate a reduction of terminal velocity in their line profiles in faster outflows, indicating destruction in the post-shock gas phase. This destruction, and a near constant ratio of the CH$_3$OH and CH$_3$CHO column densities imply it is most likely that CH$_3$CHO also forms on the grain surface.

We present the first release of the data and compact-source catalogue for the JCMT Large Program SCUBA-2 Continuum Observations of Pre-protostellar Evolution (SCOPE). SCOPE consists of 850-um continuum observations of 1235 Planck Galactic Cold Clumps (PGCCs) made with the Submillimetre Common-User Bolometer Array 2 on the James Clerk Maxwell Telescope. These data are at an angular resolution of 14.4 arcsec, significantly improving upon the 353-GHz resolution of Planck at 5 arcmin, and allowing for a catalogue of 3528 compact sources in 558 PGCCs. We find that the detected PGCCs have significant sub-structure, with 61 per cent of detected PGCCs having 3 or more compact sources, with filamentary structure also prevalent within the sample. A detection rate of 45 per cent is found across the survey, which is 95 per cent complete to Planck column densities of $N_{H_{2}}$ $>$ 5 $\times$ 10$^{21}$ cm$^{-2}$. By positionally associating the SCOPE compact sources with YSOs, the star formation efficiency, as measured by the ratio of luminosity to mass, in nearby clouds is found to be similar to that in the more distant Galactic Plane, with the column density distributions also indistinguishable from each other.

To probe the initial conditions for high-mass star and cluster formation, we investigate the properties of dense filaments within the infrared dark cloud G035.39-00.33 (IRDC G035.39) in a combined Very Large Array (VLA) and the Green Bank Telescope (GBT) mosaic tracing the NH3 (1,1) and (2,2) emission down to 0.08 pc scales. Using agglomerative hierarchical clustering on multiple line-of-sight velocity component fitting results, we identify seven extended velocity-coherent components in our data, likely representing spatially coherent physical structures, some exhibiting complex gas motions. The velocity gradient magnitude distribution peaks at its mode of 0.35 km/s/pc and has a long tail extending into higher values of 1.5 - 2 km/s/pc, and is generally consistent with those found toward the same cloud in other molecular tracers and with the values found towards nearby low-mass dense cloud cores at the same scales. Contrary to observational and theoretical expectations, we find the non-thermal ammonia line widths to be systematically narrower (by about 20%) than those of N2H+ (1-0) line transition observed with similar resolution. If the observed ordered velocity gradients represent the core envelope solid-body rotation, we estimate the specific angular momentum to be about 2 x 10^21 cm^2/s, similar to the low-mass star-forming cores. Together with the previous finding of subsonic motions in G035.39, our results demonstrate high levels of similarity between kinematics of a high-mass star-forming IRDC and the low-mass star formation regime.

For decades, the detection of phosphorus-bearing molecules in the interstellar medium was restricted to high-mass star-forming regions (as e.g. SgrB2 and Orion KL) and the circumstellar envelopes of evolved stars. However, recent higher-sensitivity observations have revealed that molecules such as PN and PO are present not only toward cold massive cores and low-mass star-forming regions with PO/PN ratios >1, but also toward the Giant Molecular Clouds in the Galactic Center known to be exposed to highly energetic phenomena such as intense UV radiation fields, shock waves and cosmic rays. In this paper, we carry out a comprehensive study of the chemistry of phosphorus-bearing molecules across different astrophysical environments which cover a range of physical conditions (cold molecular dark clouds, warm clouds, hot cores/hot corinos) and are exposed to different physical processes and energetic phenomena (proto-stellar heating, shock waves, intense UV radiation and cosmic rays). We show how the measured PO/PN ratio (either >1 as in e.g. hot molecular cores, or <1 as in UV strongly illuminated environments) can provide constraints on the physical conditions and energetic processing of the source. We propose that the reaction P + OH --> PO + H, not included in previous works, could be an efficient gas-phase PO formation route in shocks. Our modelling provides a template with which to study the detectability of P-bearing species not only in regions in our own Galaxy but also in extragalactic sources.

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

High-mass star forming regions are typically thought to be dominated by supersonic motions. We present combined Very Large Array and Green Bank Telescope (VLA+GBT) observations of NH$_3$ (1,1) and (2,2) in the infrared dark cloud (IRDC) G035.39-00.33, tracing cold and dense gas down to scales of 0.07 pc. We find that, in contrast to previous similar studies of IRDCs, more than a third of the fitted ammonia spectra show subsonic non-thermal motions (mean line width of 0.71 $\mathrm{km~s^{-1}}$), and the sonic Mach number distribution peaks around $\mathcal{M} = 1$. As possible observational and instrumental biases would only broaden the line profiles, our results provide strong upper limits to the actual value of $\mathcal{M}$, further strengthening our findings of narrow line widths. This finding calls for a reevaluation of the role of turbulent dissipation and subsonic regions in massive-star and cluster formation. Based on our findings in G035.39, we further speculate that the coarser spectral resolution used in the previous VLA NH$_3$ studies may have inhibited the detection of subsonic turbulence in IRDCs. The reduced turbulent support suggests that dynamically important magnetic fields of the 1 mG order would be required to support against possible gravitational collapse. Our results offer valuable input into the theories and simulations that aim to recreate the initial conditions of high-mass star and cluster formation.

Infrared dark clouds (IRDCs) are thought to be potential hosts of the elusive early phases of high-mass star formation. Here we conduct an in-depth kinematic analysis of one such IRDC, G034.43+00.24 (Cloud F), using high sensitivity and high spectral resolution IRAM-30m N$_2$H$^+$ ($1-0$) and C$^{18}$O ($1-0$) observations. To disentangle the complex velocity structure within this cloud we use Gaussian decomposition and hierarchical clustering algorithms. We find that four distinct coherent velocity components are present within Cloud F. The properties of these components are compared to those found in a similar IRDC, G035.39-00.33 (Cloud H). We find that the components in both clouds have: high densities (inferred by their identification in N$_2$H$^+$), trans-to-supersonic non-thermal velocity dispersions with Mach numbers of $\sim$ $1.5-4$, a separation in velocity of $\sim$3 km s$^{-1}$, and a mean red-shift of $\sim$ 0.3 km s$^{-1}$ between the N$_2$H$^+$ (dense gas) and C$^{18}$O emission (envelope gas). The latter of these could suggest that these clouds share a common formation scenario. We investigate the kinematics of the larger-scale Cloud F structures, using lower-density-tracing $^{13}$CO ($1-0$) observations. A good correspondence is found between the components identified in the IRAM-30m observations and the most prominent component in the$^{13}$CO data. We find that the IRDC Cloud F is only a small part of a much larger structure, which appears to be an inter-arm filament of the Milky Way.

We present here a systematic search for cyanopolyynes in the shock region L1157-B1 and its associated protostar L1157-mm in the framework of the Large Program "Astrochemical Surveys At IRAM" (ASAI), dedicated to chemical surveys of solar-type star forming regions with the IRAM 30m telescope. Observations of the millimeter windows between 72 and 272 GHz permitted the detection of HC$_3$N and its $^{13}$C isotopologues, and HC$_5$N (for the first time in a protostellar shock region). In the shock, analysis of the line profiles shows that the emission arises from the outflow cavities associated with L1157-B1 and L1157-B2. Molecular abundances and excitation conditions were obtained from analysis of the Spectral Line Energy Distributions under the assumption of Local Thermodynamical Equilibrium or using a radiative transfer code in the Large Velocity Gradient approximation. Towards L1157mm, the HC$_3$N emission arises from the cold envelope ($T_{rot}=10$ K) and a higher-excitation region ($T_{rot}$= $31$ K) of smaller extent around the protostar. We did not find any evidence of $^{13}$C or D fractionation enrichment towards L1157-B1. We obtain a relative abundance ratio HC$_3$N/HC$_5$N of 3.3 in the shocked gas. We find an increase by a factor of 30 of the HC$_3$N abundance between the envelope of L1157-mm and the shock region itself. Altogether, these results are consistent with a scenario in which the bulk of HC$_3$N was produced by means of gas phase reactions in the passage of the shock. This scenario is supported by the predictions of a parametric shock code coupled with the chemical model UCL_CHEM.

Infrared-Dark Clouds (IRDCs) are cold, dense regions of high (optical and infrared) extinction, believed to be the birthplace of high-mass stars and stellar clusters. The physical mechanisms leading to the formation of these IRDCs are not completely understood and it is thus important to study their molecular gas kinematics and chemical content to search for any signature of the IRDCs formation process. Using the 30m-diameter antenna at the Instituto de Radioastronom\'ia Milim\'etrica, we have obtained emission maps of dense gas tracers (H$^{13}$CO$^{+}$ and HN$^{13}$C) and typical shock tracers (SiO and CH$_3$OH) toward three IRDCs, G028.37+00.07, G034.43+00.24 and G034.77-00.55 (clouds C, F and G, respectively). We have studied the molecular gas kinematics in these clouds and, consistent with previous works toward other IRDCs, the clouds show complex gas kinematics with several velocity-coherent sub-structures separated in velocity space by a few km s$^{-1}$. Correlated with these complex kinematic structures, widespread (parsec-scale) emission of SiO and CH$_3$OH is present in all the three clouds. For clouds C and F, known to be actively forming stars, widespread SiO and CH$_3$OH is likely associated with on-going star formation activity. However, for cloud G, which lacks either 8 $\mu$m or 24 $\mu$m sources and 4.5 $\mu$m H$_2$ shock-excited emission, the detected widespread SiO and CH$_3$OH emission may have originated in a large-scale shock interaction, although a scenario involving a population of low-mass stars driving molecular outflows cannot be fully ruled out.

Proxima Centauri, the star closest to our Sun, is known to host at least one terrestrial planet candidate in a temperate orbit. Here we report the ALMA detection of the star at 1.3 mm wavelength and the discovery of a belt of dust orbiting around it at distances ranging between 1 and 4 au, approximately. Given the low luminosity of the Proxima Centauri star, we estimate a characteristic temperature of about 40 K for this dust, which might constitute the dust component of a small-scale analog to our solar system Kuiper belt. The estimated total mass, including dust and bodies up to 50 km in size, is of the order of 0.01 Earth masses, which is similar to that of the solar Kuiper belt. Our data also show a hint of warmer dust closer to the star. We also find signs of two additional features that might be associated with the Proxima Centauri system, which, however, still require further observations to be confirmed: an outer extremely cold (about 10 K) belt around the star at about 30 au, whose orbital plane is tilted about 45 degrees with respect to the plane of the sky; and additionally, we marginally detect a compact 1.3 mm emission source at a projected distance of about 1.2 arcsec from the star, whose nature is still unknown.

Complex organic molecules have been observed for decades in the interstellar medium. Some of them might be considered as small bricks of the macromolecules at the base of terrestrial life. It is hence particularly important to understand organic chemistry in Solar-like star forming regions. In this article, we present a new observational project: SOLIS (Seeds Of Life In Space). This is a Large Project at the IRAM-NOEMA interferometer, and its scope is to image the emission of several crucial organic molecules in a sample of Solar-like star forming regions in different evolutionary stage and environments. Here, we report the first SOLIS results, obtained from analysing the spectra of different regions of the Class 0 source NGC1333-IRAS4A, the protocluster OMC-2 FIR4, and the shock site L1157-B1. The different regions were identified based on the images of formamide (NH2CHO) and cyanodiacetylene (HC5N) lines. We discuss the observed large diversity in the molecular and organic content, both on large (3000-10000 au) and relatively small (300-1000 au) scales. Finally, we derive upper limits to the methoxy fractional abundance in the three observed regions of the same order of magnitude of that measured in few cold prestellar objects, namely ~10^-12-10^-11 with respect to H2 molecules.

The interstellar delivery of carbon atoms locked into molecules might be one of the key ingredients for the emergence of life. Cyanopolyynes are carbon chains delimited at their two extremities by an atom of hydrogen and a cyano group, so that they might be excellent reservoirs of carbon. The simplest member, HC3N, is ubiquitous in the galactic interstellar medium and found also in external galaxies. Thus, understanding the growth of cyanopolyynes in regions forming stars similar to our Sun, and what affects it, is particularly relevant. In the framework of the IRAM/NOEMA Large Program SOLIS (Seeds Of Life In Space), we have obtained a map of two cyanopolyynes, HC3N and HC5N, in the protocluster OMC2-FIR4. Because our Sun is thought to be born in a rich cluster, OMC2-FIR4 is one of the closest and best known representatives of the environment in which the Sun may have been born. We find a HC3N/HC5N abundance ratio across the source in the range ~ 1 - 30, with the smallest values (< 10) in FIR5 and in the Eastern region of FIR4. The ratios < 10 can be reproduced by chemical models only if: (1) the cosmic-ray ionisation rate $\zeta$ is ~ $4 \times 10^{-14}$ s$^{-1}$; (2) the gaseous elemental ratio C/O is close to unity; (3) O and C are largely depleted. The large $\zeta$ is comparable to that measured in FIR4 by previous works and was interpreted as due to a flux of energetic (> 10 MeV) particles from embedded sources. We suggest that these sources could lie East of FIR4 and FIR5. A temperature gradient across FIR4, with T decreasing by about 10 K, could also explain the observed change in the HC3N/HC5N line ratio, without the need of a cosmic ray ionisation rate gradient. However, even in this case, a high constant cosmic-ray ionisation rate (of the order of $10^{-14}$ s$^{-1}$) is necessary to reproduce the observations.

Aims. Infrared dark clouds represent the earliest stages of high-mass star formation. Detailed observations of their physical conditions on all physical scales are required to improve our understanding of their role in fueling star formation. Methods. We investigate the large-scale structure of the IRDC G035.39-00.33, probing the dense gas with the classical ammonia thermometer. This allows us to put reliable constraints on the temperature of the extended, pc-scale dense gas reservoir and to probe the magnitude of its non-thermal motions. Available far-infrared observations can be used in tandem with the observed ammonia emission to estimate the total gas mass contained in G035.39-00.33. Results. We identify a main velocity component as a prominent filament, manifested as an ammonia emission intensity ridge spanning more than 6 pc, consistent with the previous studies on the Northern part of the cloud. A number of additional line-of-sight components are found, and a large scale, linear velocity gradient of ~0.2 km s$^{-1}$ pc$^{-1}$ is found along the ridge of the IRDC. In contrast to the dust temperature map, an ammonia-derived kinetic temperature map, presented for the entirety of the cloud, reveals local temperature enhancements towards the massive protostellar cores. We show that without properly accounting for the line of sight contamination, the dust temperature is 2-3 K larger than the gas temperature measured with NH$_3$. Conclusions. While both the large scale kinematics and temperature structure are consistent with that of starless dark filaments, the kinetic gas temperature profile on smaller scales is suggestive of tracing the heating mechanism coincident with the locations of massive protostellar cores.

Nitrogen is one of the most abundant elements in the Universe and its 14N/15N isotopic ratio has the potential to provide information about the initial environment in which our Sun formed. Recent findings suggest that the Solar System may have formed in a massive cluster since the presence of short-lived radioisotopes in meteorites can only be explained by the influence of a supernova. The aim of this project is to determine the 14N/15N ratio towards a sample of cold, massive dense cores at the initial stages in their evolution. We have observed the J=1-0 transitions of HCN, H13CN, HC15N, HN13C and H15NC toward a sample of 22 cores in 4 Infrared-Dark Clouds (IRDCs). IRDCs are believed to be the precursors of high-mass stars and star clusters. Assuming LTE and a temperature of 15K, the column densities of HCN, H13CN, HC15N, HN13C and H15NC are calculated and their 14N/15N ratio is determined for each core. The 14N/15N ratio measured in our sample of IRDC cores range between ~70 and >763 in HCN and between ~161 and ~541 in HNC. They are consistent with the terrestrial atmosphere (TA) and protosolar nebula (PSN) values, and with the ratios measured in low-mass pre-stellar cores. However, the 14N/15N ratios measured in cores C1, C3, F1, F2 and G2 do not agree with the results from similar studies toward the same massive cores using nitrogen bearing molecules with nitrile functional group (-CN) and nitrogen hydrides (-NH) although the ratio spread covers a similar range. Amongst the 4 IRDCs we measured relatively low 14N/15N ratios towards IRDC G which are comparable to those measured in small cosmomaterials and protoplanetary disks. The low average gas density of this cloud suggests that the gas density, rather than the gas temperature, may be the dominant parameter influencing the initial nitrogen isotopic composition in young PSN.

The formation of deuterated molecules is favoured at low temperatures and high densities. Therefore, the deuteration fraction D$_{frac}$ is expected to be enhanced in cold, dense prestellar cores and to decrease after protostellar birth. Previous studies have shown that the deuterated forms of species such as N2H+ (formed in the gas phase) and CH3OH (formed on grain surfaces) can be used as evolutionary indicators and to constrain their dominant formation processes and time-scales. Formaldehyde (H2CO) and its deuterated forms can be produced both in the gas phase and on grain surfaces. However, the relative importance of these two chemical pathways is unclear. Comparison of the deuteration fraction of H2CO with respect to that of N2H+, NH3 and CH3OH can help us to understand its formation processes and time-scales. With the new SEPIA Band 5 receiver on APEX, we have observed the J=3-2 rotational lines of HDCO and D2CO at 193 GHz and 175 GHz toward three massive star forming regions hosting objects at different evolutionary stages: two High-mass Starless Cores (HMSC), two High-mass Protostellar Objects (HMPOs), and one Ultracompact HII region (UCHII). By using previously obtained H2CO J=3-2 data, the deuteration fractions HDCO/H2CO and D2CO/HDCO are estimated. Our observations show that singly-deuterated H2CO is detected toward all sources and that the deuteration fraction of H2CO increases from the HMSC to the HMPO phase and then sharply decreases in the latest evolutionary stage (UCHII). The doubly-deuterated form of H2CO is detected only in the earlier evolutionary stages with D2CO/H2CO showing a pattern that is qualitatively consistent with that of HDCO/H2CO, within current uncertainties. Our initial results show that H2CO may display a similar D$_{frac}$ pattern as that of CH3OH in massive young stellar objects. This finding suggests that solid state reactions dominate its formation.

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.

High-mass stars shape the interstellar medium in galaxies, and yet, largely because the initial conditions are poorly constrained, we do not know how they form. One possibility is that high-mass stars and star clusters form at the junction of filamentary networks, referred to as 'hubs'. In this letter we present the complex anatomy of a protocluster hub within an Infrared Dark Cloud (IRDC), G035.39-00.33, believed to be in an early phase of its evolution. We use high-angular resolution ($\{\theta_{\rm maj}, \theta_{\rm min}\}=\{1.4{\rm arcsec}, 0.8{\rm arcsec}\}\sim\{0.02{\rm pc}, 0.01{\rm pc}\}$) and high-sensitivity ($0.2$mJy beam$^{-1}$; $\sim0.2$M$_{\odot}$) 1.07 mm dust continuum observations from the Atacama Large Millimeter Array (ALMA) to identify a network of narrow, $0.028\,\pm\,0.005$pc wide, filamentary structures. These are a factor of $\gtrsim3$ narrower than the proposed 'quasi-universal' $\sim0.1$pc width of interstellar filaments. Additionally, 28 compact objects are reported, spanning a mass range $0.3{\rm M_{\odot}}<M_{\rm c}<10.4{\rm M_{\odot}}$. This indicates that at least some low-mass objects are forming coevally with more massive counterparts. Comparing to the popular 'bead-on-a-string' analogy, the protocluster hub is poorly represented by a monolithic clump embedded within a single filament. Instead, it comprises multiple intra-hub filaments, each of which retains its integrity as an independent structure and possesses its own embedded core population.

The detection of complex organic molecules (COMs) toward cold sources such as pre-stellar cores (with T<10 K), has challenged our understanding of the formation processes of COMs in the interstellar medium. Recent modelling on COM chemistry at low temperatures has provided new insight into these processes predicting that COM formation depends strongly on parameters such as visual extinction and the level of CO freeze out. We report deep observations of COMs toward two positions in the L1544 pre-stellar core: the dense, highly-extinguished continuum peak with Av>=30 mag within the inner 2700 au; and a low-density shell with average Av~7.5-8 mag located at 4000 au from the core's center and bright in CH3OH. Our observations show that CH3O, CH3OCH3 and CH3CHO are more abundant (by factors ~2-10) toward the low-density shell than toward the continuum peak. Other COMs such as CH3OCHO, c-C3H2O, HCCCHO, CH2CHCN and HCCNC show slight enhancements (by factors <=3) but the associated uncertainties are large. This suggests that COMs are actively formed and already present in the low-density shells of pre-stellar cores. The modelling of the chemistry of O-bearing COMs in L1544 indicates that these species are enhanced in this shell because i) CO starts freezing out onto dust grains driving an active surface chemistry; ii) the visual extinction is sufficiently high to prevent the UV photo-dissociation of COMs by the external interstellar radiation field; and iii) the density is still moderate to prevent severe depletion of COMs onto grains.

We present 3.7 arcsec (~0.05 pc) resolution 3.2 mm dust continuum observations from the IRAM PdBI, with the aim of studying the structure and fragmentation of the filamentary Infrared Dark Cloud G035.39-00.33. The continuum emission is segmented into a series of 13 quasi-regularly spaced (~0.18pc) cores, following the major axis of the IRDC. We compare the spatial distribution of the cores with that predicted by theoretical work describing the fragmentation of hydrodynamic fluid cylinders, finding a significant (a factor of ~8) discrepancy between the two. Our observations are consistent with the picture emerging from kinematic studies of molecular clouds suggesting that the cores are harboured within a complex network of independent sub-filaments. This result emphasizes the importance of considering the underlying physical structure, and potentially, dynamically important magnetic fields, in any fragmentation analysis. The identified cores exhibit a range in (peak) beam-averaged column density ($3.6{\rm x}10^{23}{\rm cm}^{-2}<N_{H,c}<8.0{\rm x}10^{23}{\rm cm}^{-2}$), mass ($8.1M_{\odot}<M_{c}<26.1M_{\odot}$), and number density ($6.1{\rm x}10^{5}{\rm cm}^{-3}<n_{H, c, eq}<14.7{\rm x}10^{5}{\rm cm}^{-3}$). Two of these cores, dark in the mid-infrared, centrally-concentrated, monolithic (with no traceable substructure at our PdBI resolution), and with estimated masses of the order ~20-25$M_{\odot}$, are good candidates for the progenitors of intermediate-to-high-mass stars. Virial parameters span a range $0.2<\alpha_{\rm vir}<1.3$. Without additional support, possibly from dynamically important magnetic fields with strengths of the order 230$\mu$G<B<670$\mu$G, the cores are susceptible to gravitational collapse. These results may imply a multilayered fragmentation process, which incorporates the formation of sub-filaments, embedded cores, and the possibility of further fragmentation.

Sulfur-bearing molecules are highly reactive in the gas phase of the ISM. However, the form in which most of the sulfur is locked onto interstellar dust grains is unknown. By taking advantage of the short time-scales of shocks in young molecular outflows, one could track back the main form of sulfur in the ices. In this paper, six transitions of H$_2$S and its isotopologues in the L1157-B1 bowshock have been detected using data from the Herschel-CHESS survey and the IRAM-30m ASAI large program. These detections are used to calculate the properties of H$_2$S gas in L1157-B1 through use of a rotation diagram and to explore the possible carriers of sulfur on the grains. The isotopologue detections allow the first calculation of the H$_2$S deuteration fraction in an outflow from a low mass protostar. The fractional abundance of H$_2$S in the region is found to be 6.0$\times$10$^{-7}$ and the deuteration fraction is 2$\times$10$^{-2}$. In order to investigate the form of sulfur on the grains, a chemical model is run with four different networks, each with different branching ratios for the freeze out of sulfur bearing species into molecules such as OCS and H$_2$S. It is found that the model best fits the data when at least half of each sulfur bearing species hydrogenates when freezing. We therefore conclude that a significant fraction of sulfur in L1157-B1 is likely to be locked in H$_2$S on the grains.

As part of the Large Program ASAI (Astrochemical Surveys At IRAM), we have used the IRAM 30m telescope to lead a systematic search for the emission of rotational transitions of P-bearing species between 80 and 350 GHz towards L1157-B1, a shock position in the solar-type star forming region L1157. We report the detection of several transitions of PN and, for the first time, of prebiotic molecule PO. None of these species are detected towards the driving protostar of the outflow L1157-mm. Analysis of the line profiles shows that PN arises from the outflow cavity, where SiO, a strong shock tracer, is produced. Radiative transfer analysis yields an abundance of 2.5e-9 and 0.9e-9 for PO and PN, respectively. These results imply a strong depletion (approx 100) of Phosphorus in the quiescent cloud gas. Shock modelling shows that atomic N plays a major role in the chemistry of PO and PN. The relative abundance of PO and PN brings constraints both on the duration of the pre-shock phase, which has to be about 1 Myr, and on the shock parameters. The maximum temperature in the shock has to be larger than 4000K, which implies a shock velocity of 40 km/s.

Investigations into the substructure of massive star forming regions are essential for understanding the observed relationships between core mass distributions and mass distributions in stellar clusters, differentiating between proposed mechanisms of massive star formation. We study the substructure in the two largest fragments (i.e. cores) MM1 and MM2, in the infrared dark cloud complex SDC13. As MM1 appears to be in a later stage of evolution than MM2, comparing their substructure provides an insight in to the early evolution of massive clumps. We report the results of high resolution SMA dust continuum observations towards MM1 and MM2. Combining these data with Herschel observations, we carry out RADMC-3D radiative transfer modelling to characterise the observed substructure. SMA continuum data indicates 4 sub-fragments in the SDC13 region. The nature of the second brightest sub-fragment (B) is uncertain as it does not appear as prominent at the lower MAMBO resolution or at radio wavelengths. Statistical analysis indicates that it is unlikely to be a background source, an AGB star, or the free-free emission of a HII region. It is plausible that B is a runaway object ejected from MM1. MM1, which is actively forming stars, consists of two sub-fragments A and C. This is confirmed by 70 micron Herschel data. While MM1 and MM2 appear quite similar in previous low resolution observations, at high resolution, the sub-fragment at the centre of MM2 (D) is much fainter than sub-fragment at the centre of MM1 (A). RADMC-3D models of MM1 and MM2 are able to reproduce these results, modelling MM2 with a steeper density profile and higher mass than is required for MM1. The relatively steep density profile of MM2 depends on a significant temperature decrease in its centre, justified by the lack of star formation in MM2. A final stellar population for MM1 was extrapolated, indicating a ...

Previous studies of the initial conditions of massive star formation have mainly targeted Infrared-Dark Clouds (IRDCs) toward the inner Galaxy. This is due to the fact that IRDCs were first detected in absorption against the bright mid-IR background, requiring a favourable location to be observed. By selection, IRDCs represent only a fraction of the Galactic clouds capable of forming massive stars and star clusters. Due to their low dust temperatures, IRDCs are bright in the far-IR and millimeter and thus, observations at these wavelengths have the potential to provide a complete sample of star-forming massive clouds across the Galaxy. Our aim is to identify the clouds at the initial conditions of massive star formation across the Galaxy and compare their physical properties as a function of their Galactic location. We have examined the physical properties of a homogeneous galactic cold core sample obtained with the Planck satellite across the Galactic Plane. With the use of Herschel Hi-GAL observations, we have characterized the internal structure of them. By using background-subtracted Herschel images, we have derived the H2 column density and dust temperature maps for 48 Planck clumps. Their basic physical parameters have been calculated and analyzed as a function of location within the Galaxy. These properties have also been compared with the empirical relation for massive star formation derived by Kauffmann & Pillai (2010). Most of the Planck clumps contain signs of star formation. About 25% of them are massive enough to form high mass stars. Planck clumps toward the Galactic center region show higher peak column densities and higher average dust temperatures than those of the clumps in the outer Galaxy. Although we only have seven clumps without associated YSOs, the Hi-GAL data show no apparent differences in the properties of Planck cold clumps with and without star formation.

Infrared Dark Clouds (IRDCs) are cold, dense regions that are usually found within Giant Molecular Clouds (GMCs). Ongoing star formation within IRDCs is typically still deeply embedded within the surrounding molecular gas. Characterising the properties of relatively quiescent IRDCs may therefore help us to understand the earliest phases of the star formation process. Studies of local molecular clouds have revealed that deuterated species are enhanced in the earliest phases of star formation. In this paper we test this towards IRDC G035.39-00.33. We present an 80 arcsec by 140 arcsec map of the J=2-1 transition of N2D+, obtained with the IRAM-30m telescope. We find that N2D+ is widespread throughout G035.39-00.33. Complementary observations of N2H+(1-0) are used to estimate the deuterium fraction, N(N2D+)/N(N2H+). We report a mean deuterium fraction of 0.04+-0.01, with a maximum of 0.09+-0.02. The mean deuterium fraction is ~3 orders of magnitude greater than the interstellar [D]/[H] ratio. High angular resolution observations are required to exclude beam dilution effects of compact deuterated cores. Using chemical modelling, we find that the average observed values of the deuterium fraction are in agreement with an equilibrium deuterium fraction, given the general properties of the cloud. This implies that the IRDC is at least ~3Myr old, which is ~8 times longer than the mean free-fall time of the observed deuterated region.

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 present the first high-angular resolution study of the MonR2 star-forming complex carried out with the Submillimeter Array at (sub-)millimeter wavelengths. We image the continuum and molecular line emission toward the young stellar objects in MonR2 at 0.85mm and 1.3mm, with resolutions ranging from 0.5" to ~3". While free-free emission dominates the IRS1 and IRS2 continuum, dust thermal emission prevails for IRS3 and IRS5, giving envelope masses of ~0.1-0.3 M_Sun. IRS5 splits into at least two sub-arcsecond scale sources, IRS5B and the more massive IRS5A. Our 12CO(2-1) images reveal 11 previously unknown molecular outflows in the MonR2 clump. Comparing these outflows with known IR sources in the IRS5 and IRS3 subclusters allows for tentative identification of driving stars. Line images of molecular species such as CH3CN or CH3OH show that, besides IRS3 (a well-known hot molecular core), IRS5 is also a chemically active source in the region. The gas excitation temperature derived from CH3CN lines toward IRS5 is 144 \pm 15 K, indicating a deeply embedded protostar at the hot-core evolutionary stage. SED fitting of IRS5 gives a mass of ~7 M_Sun and a luminosity of 300 L_Sun for the central source. The derived physical properties of the CO outflows suggest that they contribute to the turbulent support of the MonR2 complex and to the gas velocity dispersion in the clump's center. The detection of a large number of CO outflows widespread across the region supports the competitive accretion scenario as origin of the MonR2 star cluster.

Sufficiently massive clumps of molecular gas collapse under self-gravity and fragment to spawn a cluster of stars that have a range of masses. We investigate observationally the early stages of formation of a stellar cluster in a massive filamentary infrared dark cloud, G28.34+0.06 P1, in the 1.3mm continuum and spectral line emission using the ALMA. Sensitive continuum data reveal further fragmentation in five dusty cores at a resolution of several 10^3 AU. Spectral line emission from C18O, CH3OH, 13CS, H2CO and N2D+ are detected for the first time toward these dense cores. We found that three cores are chemically more evolved as compared with the other two; interestingly though, all of them are associated with collimated outflows as suggested by evidence from the CO, SiO, CH3OH, H2CO and SO emissions. The parsec-scale kinematics in NH3 exhibit velocity gradients along the filament, consistent with accretion flows toward the clumps and cores. The moderate luminosity and the chemical signatures indicate that the five cores harbor low- to intermediate-mass protostars that likely become massive ones at the end of the accretion. Despite the fact that the mass limit reached by the 1\sigma dust continuum sensitivity is 30 times lower than the thermal Jeans mass, there is a lack of a distributed low-mass protostellar population in the clump. Our observations indicate that in a protocluster, low-mass stars form at a later stage after the birth of more massive protostars.

Glycine (NH2CH2COOH) is the simplest amino acid relevant for life. Its detection in the interstellar medium is key to understand the formation mechanisms of pre-biotic molecules and their subsequent delivery onto planetary systems. Glycine has extensively been searched for toward hot molecular cores, although these studies did not yield any firm detection. In contrast to hot cores, low-mass star forming regions, and in particular their earliest stages represented by cold pre-stellar cores, may be better suited for the detection of glycine as well as more relevant for the study of pre-biotic chemistry in young Solar System analogs. We present 1D spherically symmetric radiative transfer calculations of the glycine emission expected to arise from the low-mass pre-stellar core L1544. Water vapour has recently been reported toward this core, indicating that a small fraction of the grain mantles in L1544 (~0.5%) has been injected into the gas phase. Assuming that glycine is photo-desorbed together with water in L1544, and considering a solid abundance of glycine on ices of ~1E-4 with respect to water, our calculations reveal that several glycine lines between 67 GHz and 80 GHz have peak intensities larger than 10 mK. These results show for the first time that glycine could reach detectable levels in cold objects such as L1544. This opens up the possibility to detect glycine, and other pre-biotic species, at the coldest and earliest stages in the formation of Solar-type systems with near-future instrumentation such as the Band 2 receivers of ALMA.

Infrared Dark Clouds (IRDCs) are unique laboratories to study the initial conditions of high-mass star and star cluster formation. We present high-sensitivity and high-angular resolution IRAM PdBI observations of N2H+ (1-0) towards IRDC G035.39-00.33. It is found that G035.39-00.33 is a highly complex environment, consisting of several mildly supersonic filaments (sigma_NT/c_s ~1.5), separated in velocity by <1 km s^-1 . Where multiple spectral components are evident, moment analysis overestimates the non-thermal contribution to the line-width by a factor ~2. Large-scale velocity gradients evident in previous single-dish maps may be explained by the presence of substructure now evident in the interferometric maps. Whilst global velocity gradients are small (<0.7 km s^-1 pc^-1), there is evidence for dynamic processes on local scales (~1.5-2.5 km s^-1 pc^-1 ). Systematic trends in velocity gradient are observed towards several continuum peaks. This suggests that the kinematics are influenced by dense (and in some cases, starless) cores. These trends are interpreted as either infalling material, with accretion rates ~(7 \pm 4)x10^-5 M_sun yr^-1 , or expanding shells with momentum ~24 \pm 12 M_sun km s^-1 . These observations highlight the importance of high-sensitivity and high-spectral resolution data in disentangling the complex kinematic and physical structure of massive star forming regions.

Some theories of dense molecular cloud formation involve dynamical environments driven by converging atomic flows or collisions between preexisting molecular clouds. The determination of the dynamics and physical conditions of the gas in clouds at the early stages of their evolution is essential to establish the dynamical imprints of such collisions, and to infer the processes involved in their formation. We present multi-transition 13CO and C18O maps toward the IRDC G035.39-00.33, believed to be at the earliest stages of evolution. The 13CO and C18O gas is distributed in three filaments (Filaments 1, 2 and 3), where the most massive cores are preferentially found at the intersecting regions between them. The filaments have a similar kinematic structure with smooth velocity gradients of ~0.4-0.8 km s-1 pc-1. Several scenarios are proposed to explain these gradients, including cloud rotation, gas accretion along the filaments, global gravitational collapse, and unresolved sub-filament structures. These results are complemented by HCO+, HNC, H13CO+ and HN13C single-pointing data to search for gas infall signatures. The 13CO and C18O gas motions are supersonic across G035.39-00.33, with the emission showing broader linewidths toward the edges of the IRDC. This could be due to energy dissipation at the densest regions in the cloud. The average H2 densities are ~5000-7000 cm-3, with Filaments 2 and 3 being denser and more massive than Filament 1. The C18O data unveils three regions with high CO depletion factors (f_D~5-12), similar to those found in massive starless cores.

We have studied the young low-mass pre-main sequence (PMS) stellar population associated with the massive star-forming region DR 21 by using archival X-ray Chandra observations and by complementing them with existing optical and IR surveys. The Chandra observations have revealed for the first time a new highly extincted population of PMS low-mass stars previously missed in observations at other wavelengths. The X-ray population exhibits three main stellar density peaks, coincident with the massive star-forming regions, being the DR 21 core the main peak. The cross-correlated X-ray/IR sample exhibits a radial "Spokes-like" stellar filamentary structure that extends from the DR 21 core towards the northeast. The near IR data reveal a centrally peaked structure for the extinction, which exhibits its maximum in the DR 21 core and gradually decreases with the distance to the N-S cloud axis and to the cluster center. We find evidence of a global mass segregation in the full low-mass stellar cluster, and of an stellar age segregation, with the youngest stars still embedded in the N-S cloud, and more evolved stars more spatially distributed. The results are consistent with the scenario where an elongated overall potential well created by the full low-mass stellar cluster funnels gas through filaments feeding stellar formation. Besides the full gravitational well, smaller-scale local potential wells created by dense stellar sub-clusters of low-mass stars are privileged in the competition for the gas of the common reservoir, allowing the formation of massive stars. We also discuss the possibility that a stellar collision in the very dense stellar cluster revealed by Chandra in the DR 21 core is the origin of the large-scale and highly-energetic outflow arising from this region.

To distinguish between the different theories proposed to explain massive star formation, it is crucial to establish the distribution, the extinction, and the density of low-mass stars in massive star-forming regions. We analyzed deep X-ray observations of the Orion massive star-forming region using the Chandra Orion Ultradeep Project (COUP) catalog. We found that pre-main sequence (PMS) low-mass stars cluster toward the three massive star-forming regions: the Trapezium Cluster (TC), the Orion Hot Core (OHC), and OMC1-S. We derived low-mass stellar densities of 10^{5} stars pc^{-3} in the TC and OMC1-S, and of 10^{6} stars pc^{-3} in the OHC. The close association between the low-mass star clusters with massive star cradles supports the role of these clusters in the formation of massive stars. The X-ray observations show for the first time in the TC that low-mass stars with intermediate extinction are clustered toward the position of the most massive star, which is surrounded by a ring of non-extincted low-mass stars. Our analysis suggests that at least two basic ingredients are needed in massive star formation: the presence of dense gas and a cluster of low-mass stars. The scenario that better explains our findings assumes high fragmentation in the parental core, accretion at subcore scales that forms a low-mass stellar cluster, and subsequent competitive accretion.

Outflows arising from very young stars affect their surroundings and influence the star formation in the parental core. Multiple molecular outflows and Herbig-Haro (HH) objects have been observed in Orion, many of them originating from the embedded massive star-forming region known as OMC1-S. The detection of the outflow driving sources is commonly difficult, because they are still hidden behind large extinction, preventing their direct observation at optical and even near and mid-IR wavelengths. With the aim of improving the identification of the driving sources of the multiple outflows detected in OMC1-S, we used the catalog provided by deep X-ray observations, which have unveiled the very embedded population of pre-main sequence stars. We compared the position of stars observed by the Chandra Orion Ultra Deep project (COUP) in OMC1-S with the morphology of the molecular outflows and the directions of measured proper motions of HH optical objects. We find that 6 out of 7 molecular outflows reported in OMC1-S (detection rate of 86 %) have an extincted X-ray COUP star located at the expected position of the driving source. In several cases, X-rays detected the possible driving sources for the first time. This clustered embedded population revealed by Chandra is very young, with an estimated average age of few 10^{5} yr. It is also likely responsible for the multiple HH objects, which are the optical correspondence of flows arising from the cloud. We show that the molecular outflows driven by the members of the OMC1-S cluster can account for the observed turbulence at core-scales and regulate the star formation efficiency. We discuss the effects of outflow feedback in the formation of massive stars, concluding that the injected turbulence in OMC1-S is compatible with a competitive accretion scenario.

We have observed the IRDC G028.23-00.19 at 3.3 mm using CARMA. In its center, the IRDC hosts one of the most massive (~1520 Msun) quiescent, cold (12 K) clumps known (MM1). The low temperature, high NH_2D abundance, narrow molecular line widths, and absence of embedded IR sources (from 3.6 to 70 mu) indicate that the clump is likely prestellar. Strong SiO emission with broad line widths (6-9 km s^-1) and high abundances (0.8-4 x 10^-9) is detected in the northern and southern regions of the IRDC, unassociated with MM1. We suggest that SiO is released to the gas phase from the dust grains through shocks produced by outflows from undetected intermediate-mass stars or clusters of low-mass stars deeply embedded in the IRDC. A weaker SiO component with narrow line widths (~2 km s^-1) and low abundances (4.3 x 10^-11) is detected in the center-west region, consistent with either a "subcloud-subcloud" collision or an unresolved population of a few low-mass stars. We report widespread CH_3OH emission throughout the whole IRDC and the first detection of extended narrow methanol emission (~2 km s^-1) in a cold, massive prestellar clump (MM1). We suggest that the most likely mechanism releasing methanol into the gas phase in such a cold region is the exothermicity of grain-surface reactions. HN^13C reveals that the IRDC is actually composed of two distinct substructures ("subclouds") separated in velocity space by ~1.4 km s^-1. The narrow SiO component arises where the subclouds overlap. The spatial distribution of C_2H resembles that of NH_2D, which suggests that C_2H also traces cold gas in this IRDC.

We present 450 {\mu}m polarimetric observations of the M17 molecular cloud obtained with the SHARP polarimeter at the Caltech Submillimeter Observatory. Across the observed region, the magnetic field orientation is consistent with previous submillimeter and far-infrared polarization measurements. Our observations are centered on a region of the molecular cloud that has been compressed by stellar winds from a cluster of OB stars. We have compared these new data with previous 350 {\mu}m polarimetry and find an anti-correlation between the 450 to 350 {\mu}m polarization magnitude ratio and the ratio of 21 cm to 450 {\mu}m intensity. The polarization ratio is lower near the east end of the studied region where the cloud is exposed to stellar winds and radiation. At the west end of the region, the polarization ratio is higher. We interpret the varying polarization spectrum as evidence supporting the radiative alignment torque (RAT) model for grain alignment, implying higher alignment efficiency in the region that is exposed to a higher anisotropic radiation field.

To distinguish between the different theories proposed to explain massive star formation, it is crucial to establish the distribution, the extinction, and the density of low-mass stars in massive star-forming regions. We analyze deep X-ray observations of the Orion massive star-forming region using the Chandra Orion Ultradeep Project (COUP) catalog. We studied the stellar distribution as a function of extinction, with cells of 0.03 pc x 0.03 pc, the typical size of protostellar cores. We derived stellar density maps and calculated cluster stellar densities. We found that low-mass stars cluster toward the three massive star-forming regions: the Trapezium Cluster (TC), the Orion Hot Core (OHC), and OMC1-S. We derived low-mass stellar densities of 10^{5} stars pc^{-3} in the TC and OMC1-S, and of 10^{6} stars pc^{-3} in the OHC. The close association between the low-mass star clusters with massive star cradles supports the role of these clusters in the formation of massive stars. The X-ray observations show for the first time in the TC that low-mass stars with intermediate extinction are clustered toward the position of the most massive star, which is surrounded by a ring of non-extincted low-mass stars. This 'envelope-core' structure is also supported by infrared and optical observations. Our analysis suggests that at least two basic ingredients are needed in massive star formation: the presence of dense gas and a cluster of low-mass stars. The scenario that better explains our findings assumes high fragmentation in the parental core, accretion at subcore scales that forms a low-mass stellar cluster, and subsequent competitive accretion. Finally, although coalescence does not seem a common mechanism for building up massive stars, we show that a single stellar merger may have occurred in the evolution of the OHC cluster, favored by the presence of disks, binaries, and gas accretion.

We report the detection of a new radio recombination line (RRL) maser object toward the IRS2 source in the MonR2 ultracompact HII region. The continuum emission at 1.3mm and 0.85mm and the H30a and H26a lines were observed with the Submillimeter Array (SMA) at angular resolutions of about 0.5"-3". The SMA observations show that the MonR2-IRS2 source is very compact and remains unresolved at spatial scales <=400AU. Its continuum power spectrum at millimeter wavelengths is almost flat (alpha=-0.16, with S_nu proportional to nu^alpha), indicating that this source is dominated by optically thin free-free emission. The H30a and H26a RRL emission is also compact and peaks toward the position of the MonR2-IRS2 source. The measured RRL profiles are double-peaked with the H26a line showing a clear asymmetry in its spectrum. Since the derived line-to-continuum flux ratios (80 and 180kms-1 for H30a and H26a, respectively) exceed the LTE predictions, the RRLs toward MonR2-IRS2 are affected by maser amplification. The amplification factors are however smaller than those found toward the emission line star MWC349A, indicating that MonR2-IRS2 is a weakly amplified maser. Radiative transfer modelling of the RRL emission toward this source shows that the RRL masers arise from a dense and collimated jet embedded in a cylindrical ionized wind, oriented nearly along the direction of the line-of-sight. High-angular resolution observations at sub-millimeter wavelengths are needed to unveil weakly amplified RRL masers in very young massive stars.

Infrared Dark Clouds (IRDCs) host the initial conditions under which massive stars and stellar clusters form. We have obtained high sensitivity and high spectral resolution observations with the IRAM 30m antenna, which allowed us to perform detailed analysis of the kinematics within one IRDC, G035.39-00.33. We focus on the 1-0 and 3-2 transitions of N2H+, C18O (1-0), and make comparison with SiO (2-1) observations and extinction mapping. Three interacting filaments of gas are found. We report large-scale velocity coherence throughout the cloud, evidenced through small velocity gradients and relatively narrow line widths. This suggests that the merging of these filaments is somewhat "gentle", possibly regulated by magnetic fields. This merging of filaments may be responsible for the weak parsec-scale SiO emission detected by Jimenez-Serra et al. 2010, via grain mantle vaporization. A systematic velocity shift between the N2H+ (1-0) and C18O (1-0) gas throughout the cloud of 0.18 +/- 0.04 kms^{-1} is also found, consistent with a scenario of collisions between filaments which is still ongoing. The N2H+ (1-0) is extended throughout the IRDC and it does not only trace dense cores, as found in nearby low-mass star-forming regions. The average H2 number density across the IRDC is ~ 5 x 10^4 cm^{-3}, at least one order of magnitude larger than in nearby molecular clouds where low-mass stars are forming. A temperature gradient perpendicular to the filament is found. From our study, we conclude that G035.39-00.33 (clearly seen in the extinction map and in N2H+) has been formed via the collision between two relatively quiescent filaments with average densities of ~ 5 x 10^3 cm^{-3}, moving with relative velocities of ~ 5 kms^{-1}. The accumulation of material at the merging points started > 1 Myr ago and it is still ongoing.

The initial conditions of massive star and star cluster formation are expected to be cold, dense and high column density regions of the interstellar medium, which can reveal themselves via near, mid and even far-infrared absorption as Infrared Dark Clouds (IRDCs). Elucidating the dynamical state of IRDCs thus constrains theoretical models of these complex processes. In particular, it is important to assess whether IRDCs have reached virial equilibrium, where the internal pressure balances that due to the self-gravitating weight of the cloud plus the pressure of the external environmental. We study this question for the filamentary IRDC G035.39-00.33 by deriving mass from combined NIR & MIR extinction maps and velocity dispersion from C18O (1-0) & (2-1) line emission. In contrast to our previous moderately super-virial results based on 13CO emission and MIR-only extinction mapping, with improved mass measurements we now find that the filament is consistent with being in virial equilibrium, at least in its central parsec-wide region where ~1000 M_Sun snakes along several parsecs. This equilibrium state does not require large-scale net support or confinement by magnetic fields.

We present high angular resolution observations (0.5"x0.3") carried out with the Submillimeter Array (SMA) toward the AFGL2591 high-mass star forming region. Our SMA images reveal a clear chemical segregation within the AFGL2591 VLA 3 hot core, where different molecular species (Type I, II and III) appear distributed in three concentric shells. This is the first time that such a chemical segregation is ever reported at linear scales <3000 AU within a hot core. While Type I species (H2S and 13CS) peak at the AFGL2591 VLA 3 protostar, Type II molecules (HC3N, OCS, SO and SO2) show a double-peaked structure circumventing the continuum peak. Type III species, represented by CH3OH, form a ring-like structure surrounding the continuum emission. The excitation temperatures of SO2, HC3N and CH3OH (185+-11 K, 150+-20 K and 124+-12 K, respectively) show a temperature gradient within the AFGL2591 VLA 3 envelope, consistent with previous observations and modeling of the source. By combining the H2S, SO2 and CH3OH images, representative of the three concentric shells, we find that the global kinematics of the molecular gas follow Keplerian-like rotation around a 40 Mo-star. The chemical segregation observed toward AFGL2591 VLA 3 is explained by the combination of molecular UV photo-dissociation and a high-temperature (~1000 K) gas-phase chemistry within the low extinction innermost region in the AFGL2591 VLA 3 hot core.

We investigate the origin and nature of the profiles of water and ammonia observ ed toward the L1157-B1 clump as part of the HIFI CHESS survey (Ceccarelli et al. 2010) using a new code coupling a gas-grain chemical model with a parametric shock model. Fir st results from the unbiased survey (Lefloch et al. 2010, Codella et al. 2010) reveal different molecular components at different excitation conditions coexisting in the B1 bow shock structure, with NH$_3$, H$_2$CO and CH$_3$OH emitting only at relatively low outflow velocities whereas H$_2$O shows bright emission at high velocities. Our model suggests that these differences are purely chemical and can be explained by the presence of a C-type shock whose maximum temperature must be close to 4000 K along the B1 clump.

The detection of narrow SiO thermal emission toward young outflows has been proposed to be a signature of the magnetic precursor of C-shocks. Recent modeling of the SiO emission across C-shocks predicts variations in the SiO line intensity and line shape at the precursor and intermediate-velocity regimes in only few years. We present high-angular resolution (3.8"x3.3") images of the thermal SiO J=2-1 emission toward the L1448-mm outflow in two epochs (November 2004-February 2005, March-April 2009). Several SiO condensations have appeared at intermediate velocities (20-50 km/s) toward the red-shifted lobe of the outflow since 2005. Toward one of the condensations (clump D), systematic differences of the dirty beams between 2005 and 2009 could be responsible for the SiO variability. At higher velocities (50-80 km/s), SiO could also have experienced changes in its intensity. We propose that the SiO variability toward L1448-mm is due to a real SiO enhancement by young C-shocks at the internal working surface between the jet and the ambient gas. For the precursor regime (5.2-9.2 km/s), several narrow and faint SiO components are detected. Narrow SiO tends to be compact, transient and shows elongated (bow-shock) morphologies perpendicular to the jet. We speculate that these features are associated with the precursor of C-shocks appearing at the interface of the new SiO components seen at intermediate velocities.

Infrared Dark Clouds (IRDCs) are cold, high mass surface density and high density structures, likely to be representative of the initial conditions for massive star and star cluster formation. CO emission from IRDCs has the potential to be useful for tracing their dynamics, but may be affected by depleted gas phase abundances due to freeze-out onto dust grains. Here we analyze C18O J=1-0 and J=2-1 emission line data, taken with the IRAM 30m telescope, of the highly filamentary IRDC G035.39.-0033. We derive the excitation temperature as a function of position and velocity, with typical values of ~7K, and thus derive total mass surface densities, Sigma_C18O, assuming standard gas phase abundances and accounting for optical depth in the line, which can reach values of ~1. The mass surface densities reach values of ~0.07 g/cm^2. We compare these results to the mass surface densities derived from mid-infrared (MIR) extinction mapping, Sigma_SMF, by Butler & Tan, which are expected to be insensitive to the dust temperatures in the cloud. With a significance of >10sigma, we find Sigma_C18O/Sigma_SMF decreases by about a factor of 5 as Sigma increases from ~0.02 to ~0.2 g/cm^2, which we interpret as evidence for CO depletion. Several hundred solar masses are being affected, making this one of the most massive clouds in which CO depletion has been observed directly. We present a map of the depletion factor in the filament and discuss implications for the formation of the IRDC.

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 present the first detection of the H40a, H34a and H31a radio recombination lines (RRLs) at millimeter wavelengths toward the high-velocity, ionized jet in the Cepheus A HW2 star forming region. From our single-dish and interferometric observations, we find that the measured RRLs show extremely broad asymmetric line profiles with zero-intensity linewidths of ~1100 kms-1. From the linewidths, we estimate a terminal velocity for the ionized gas in the jet of >500 kms-1, consistent with that obtained from the proper motions of the HW2 radio jet. The total integrated line-to-continuum flux ratios of the H40a, H34a and H31a lines are 43, 229 and 280 kms-1, clearly deviating from LTE predictions. These ratios are very similar to those observed for the RRL maser toward MWC349A, suggesting that the intensities of the RRLs toward HW2 are affected by maser emission. Our radiative transfer modeling of the RRLs shows that their asymmetric profiles could be explained by maser emission arising from a bi-conical radio jet with a semi-aperture angle of 18 deg, electron density distribution varying as r^(-2.11) and turbulent and expanding wind velocities of 60 and 500 kms-1.

We present high-sensitivity 2'x4' maps of the J=2-1 rotational lines of SiO, CO, 13CO and C18O, observed toward the filamentary Infrared Dark Cloud (IRDC) G035.39-00.33. Single-pointing spectra of the SiO J=2-1 and J=3-2 lines toward several regions in the filament, are also reported. The SiO images reveal that SiO is widespread along the IRDC (size >2 pc), showing two different components: one bright and compact arising from three condensations (N, E and S), and the other weak and extended along the filament. While the first component shows broad lines (linewidths of ~4-7 kms-1) in both SiO J=2-1 and SiO J=3-2, the second one is only detected in SiO J=2-1 and has narrow lines (~0.8 kms-1). The maps of CO and its isotopologues show that low-density filaments are intersecting the IRDC and appear to merge toward the densest portion of the cloud. This resembles the molecular structures predicted by flow-driven, shock-induced and magnetically-regulated cloud formation models. As in outflows associated with low-mass star formation, the excitation temperatures and fractional abundances of SiO toward N, E and S, increase with velocity from ~6 to 40 K, and from ~1E-10 to >1E-8 respectively, over a velocity range of ~7 kms-1. Since 8 micron sources, 24 micron sources and/or extended 4.5 micron emission are detected in N, E and S, broad SiO is likely produced in outflows associated with high-mass protostars. The excitation temperatures and fractional abundances of the narrow SiO lines, however, are very low (~9 K and ~1E-11, respectively), and consistent with the processing of interstellar grains by the passage of a shock with vs~12 kms-1. This emission could be generated i) by a large-scale shock, perhaps remnant of the IRDC formation process; ii) by decelerated or recently processed gas in large-scale outflows driven by 8 micron and 24 micron sources; or iii) by an undetected and widespread population of lower mass protostars. High-angular resolution observations are needed to disentangle between these three scenarios.

We present high angular resolution PdBI images (beam of ~0.33'') of the J=27-26 line from several vibrational levels (v7=1 and v6=1) of HC3N toward Cepheus A HW2. These images reveal the two main heating sources in the cluster: one centered in the disk collimating the HW2 radio jet (the HW2 disk), and the other associated with a hot core 0.3'' northeast HW2 (the HC). This is the first time that vibrationally excited emission of HC3N is spatially resolved in a disk. The kinematics of this emission shows that the HW2 disk rotates following a Keplerian law. We derive the temperature profiles in the two objects from the excitation of HC3N along the HW2 disk and the HC. These profiles reveal that both objects are centrally heated and show temperature gradients. The inner and hotter regions have temperatures of 350+-30K and 270+-20K for the HW2 disk and the HC, respectively. In the cooler and outer regions, the temperature drops to 250+-30K in the HW2 disk, and to 220+-15K in the HC. The estimated luminosity of the heating source of the HW2 disk is ~2.2E4Lo, and 3000Lo for the HC. The most massive protostar in the HW2 region is the powering source of the HW2 radio jet. We discuss the formation of multiple systems in this cluster. The proximity of the HC to HW2 suggest that these sources likely form a binary system of B stars, explaining the observed precession of the HW2 radio jet.

We present the first results of the expected variations of the molecular line emission arising from material recently affected by C-shocks (shock precursors). Our parametric model of the structure of C-shocks has been coupled with a radiative transfer code to calculate the molecular excitation and line profiles of shock tracers such as SiO, and of ion and neutral molecules such as H13CO+ and HN13C, as the shock propagates through the unperturbed medium. Our results show that the SiO emission arising from the early stage of the magnetic precursor typically has very narrow line profiles slightly shifted in velocity with respect to the ambient cloud. This narrow emission is generated in the region where the bulk of the ion fluid has already slipped to larger velocities in the precursor as observed toward the young L1448-mm outflow. This strongly suggests that the detection of narrow SiO emission and of an ion enhancement in young shocks, is produced by the magnetic precursor of C-shocks. In addition, our model shows that the different velocity components observed toward this outflow can be explained by the coexistence of different shocks at different evolutionary stages, within the same beam of the single-dish observations.

Context: The detection of narrow SiO lines toward the young shocks of the L1448-mm outflow has been interpreted as a signature of the magnetic precursor of C-shocks. In contrast with the low SiO abundances (<10E-12) in the ambient gas, the narrow SiO emission at almost ambient velocities reveals enhanced SiO abundances of 10E-11. This enhancement has been proposed to be produced by the sputtering of the grain mantles at the first stages of C-shocks. However, modelling of the sputtering of grains has usually averaged the SiO abundances over the dissipation region of C-shocks, which cannot explain the recent observations. Aims: To model the evolution of the gas phase abundances of SiO, CH3OH and H2O, produced by the sputtering of grains as the shock propagates through the ambient gas. Methods: We propose a parametric model to describe the physical structure of C-shocks as a function of time. Using the known sputtering yields for water mantles (with minor constituents like silicon and CH3OH) and olivine cores by collisions with H2, He, C, O, Si, Fe and CO, we follow the evolution of the abundances of silicon, CH3OH and H2O ejected from grains. Results: The evolution of these abundances shows that CO seems to be the most efficient sputtering agent in low velocity shocks. The velocity threshold for the sputtering of silicon from the grain mantles is reduced by 5-10 km s-1 by CO compared to other models. The sputtering by CO can generate SiO abundances of 10E-11 at the early stages of low velocity shocks, consistent with those observed in the magnetic precursor of L1448-mm. Our model also satisfactorily reproduce the progressive enhancement of SiO, CH3OH and H2O observed in this outflow by the coexistence of two shocks with vs=30 and 60kms-1 within the same region.

We present VLA and PdBI subarcsecond images (0.15"-0.6") of the radiocontinuum emission at 7 mm and of the SO2 J=19_{2,18}-18_{3,15} and J=27_{8,20}-28_{7,21} lines toward the Cepheus A HW2 region. The SO2 images reveal the presence of a hot core internally heated by an intermediate mass protostar, and a circumstellar rotating disk around the HW2 radio jet with size 600AUx100AU and mass of 1M_sun. Keplerian rotation for the disk velocity gradient of 5 kms-1 requires a 9 M_sun central star, which cannot explain the total luminosity observed in the region. This may indicate that the disk does not rotate with a Keplerian law due to the extreme youth of this object. Our high sensitivity radiocontinuum image at 7 mm shows in addition to the ionized jet, an extended emission to the west (and marginally to the south) of the HW2 jet, filling the south-west cavity of the HW2 disk. From the morphology and location of this free-free continuum emission at centimeter and millimeter wavelengths (spectral index of 0.4-1.5), we propose that the disk is photoevaporating due to the UV radiation from the central star. All this indicates that the Cepheus A HW2 region harbors a cluster of massive stars. Disk accretion seems to be the most plausible way to form massive stars in moderate density/luminosity clusters.

Due to its relatively small distance (725 pc), the Cepheus A East star-forming region is an ideal laboratory to study massive star formation processes. Based on its morphology, it has been suggested that the flattened molecular gas distribution around the YSO HW2 may be a 350-AU-radius massive protostellar disk. Goal of our work is to ascertain the nature of this structure. We have employed the Plateau de Bure Interferometer to acquire (sub-)arcsecond-resolution imaging of high-density and shock tracers, such as methyl cyanide (CH3CN) and silicon monoxide (SiO), towards the HW2 position. On the 1-arcsecond (about 725 AU) scale, the flattened distribution of molecular gas around HW2 appears to be due to the projected superposition, on the plane of the sky, of at least three protostellar objects, of which at least one is powering a molecular outflow at a small angle with respect to the line of sight. The presence of a protostellar disk around HW2 is not ruled out, but such structure is likely to be detected on a smaller spatial scale, or using different molecular tracers.

Magnetic precursors of C-shocks accelerate, compress and heat molecular ions, modifying the kinematics and the physical conditions of the ion fluid with respect to the neutral one. Electron densities are also expected to be significantly enhanced in shock precursors. In this Letter, we present observations of strongly polar ion and neutral molecules such as SiO, H13CO+, HN13C and H13CN, which reveal the electron density enhancements associated with the precursor of the young L1448-mm outflow. While in the ambient gas the excitation of the ions and neutrals is explained by collisional excitation by H2 with a single density of 10E5cmE-3, H13CO+ shows an over excitation in the shock precursor component that requires H2 densities of a factor of >10 larger than those derived from the neutral species. This over excitation in H13CO+ can be explained if we consider an additional excitation by collisions with electrons and an electron density enhancement in the precursor stage by a factor of 500, i.e. a fractional ionization of 5x10E-5. These results show that multiline observations can be used to study the evolution of the ion and electron fluids at the first stages of the C-shock interaction.

We present the discovery of the first molecular hot core associated with an intermediate mass protostar in the CepA HW2 region. The hot condensation was detected from single dish and interferometric observations of several high excitation rotational lines (from 100 to 880K above the ground state) of SO2 in the ground vibrational state and of HC3N in the vibrationally excited states v7=1 and v7=2. The kinetic temperature derived from both molecules is 160K. The high-angular resolution observations (1.25'' x 0.99'') of the SO2 J=28(7,21)-29(6,24) line (488K above the ground state) show that the hot gas is concentrated in a compact condensation with a size of 0.6''(430AU), located 0.4'' (300AU) east from the radio-jet HW2. The total SO2 column density in the hot condensation is 10E18cm-2, with a H2 column density ranging from 10E23 to 6 x 10E24cm-2. The H2 density and the SO2 fractional abundance must be larger than 10E7cm-3 and 2 x 10E-7 respectively. The most likely alternatives for the nature of the hot and very dense condensation are discussed. From the large column densities of hot gas, the detection of the HC3N vibrationally excited lines and the large SO2 abundance, we favor the interpretation of a hot core heated by an intermediate mass protostar of 10E3 Lo. This indicates that the CepA HW2 region contains a cluster of very young stars.

The recent detection of shock-precursors toward the very young L1448-mm outflow offers us the possibility to study the grain chemistry during the first stages of the shock evolution, constraining the molecules ejected from grains and the species formed in gas phase. Observations of key molecules in the grain chemistry such as SiO, CH3OH, SO, CS, H2S, OCS, and SO2 toward this outflow are presented. The line profiles and the derived abundances show three distinct velocity regimes that trace the shock evolution: the preshock, the shock-precursor and the postshock gas. The SiO, CH3OH, SO, and CS abundances are enhanced with respect to the quiescent gas by 1 order of magnitude in the shock-precursor component, and by 3 orders of magnitude in the postshock gas. The derived SiO and CH3OH abundances are consistent with the recent ejection of these molecules from grains. Since H2S is only enhanced in the shock-precursor component, and OCS and SO2 are undetected, SO and CS are the most abundant sulfur-bearing species in the grain mantles of L1448-mm. The ejection of mainly SO and CS rather than H2S or OCS from grains, suggests that the sulfur chemistry will depend on the chemical "history" of the grain mantles in outflows and hot cores.

We present the detection of the SiO v=0 J=2-1 and J=3-2 lines, and of the HCO 101-000 J=3/2-1/2 F=2-1 line at ambient velocities towards the molecular outflows in L1448-mm and L1448-IRS3. This is the first detection of HCO in a dark cloud. We have also measured lines of H13CO+, H13CN, HN13C, CH3OH, and N2H+. While the HCO and the SiO lines have the narrowest profiles with linewidths of 0.5 km s-1, the other lines have widths of 1km s-1. Towards L1448-mm, all lines except those of SiO and HCO, show two distinct velocity components centered at 4.7 and 5.2 km s-1. HCO is only observed in the 4.7 km s-1 cloud, and SiO in the 5.2 km s-1 component. The SiO abundance is 10E-11 in the 5.2 km s-1 clouds, one order of magnitude larger than in the 4.7 km s-1 component and in other dark clouds. The HCO abundance is 10E-11, similar to that predicted by the ion-molecule reactions models for the quiescent gas in dark clouds. The large change in the SiO/HCO abundance ratio (>150) from the 4.7 to the 5.2 km s-1 component, and the distribution and kinematics of the SiO emission towards L1448-mm suggest that the ambient SiO is associated with the molecular outflows. We propose that the narrow linewidths and the abundances of SiO in the ambient gas are produced by the interaction of the magnetic and/or radiative precursors of the shocks with the clumpy pre-shocked ambient gas.