We have studied four complex organic molecules (COMs), methyl formate ($CH_3OCHO$), dimethyl ether ($CH_3OCH_3$), formamide ($NH_2CHO$), and ethyl cyanide ($C_2H_5CN$), towards a large sample of 39 high-mass star-forming regions representing different evolutionary stages, from early to evolved phases. We aim to identify potential correlations between the molecules and to trace their evolutionary sequence through the star formation process. We analysed spectra obtained at 3, 2, and 0.9 mm with the IRAM-30m telescope. We derived the main physical parameters for each species by fitting the molecular lines. We compared them and evaluated their evolution, also taking several other interstellar environments into account. We report detections in 20 sources, revealing a clear dust absorption effect on column densities. Derived abundances are ~$10^{-10}-10^{-7}$ for $CH_3OCHO$ and $CH_3OCH_3$, ~$10^{-12}-10^{-10}$ for $NH_2CHO$, and ~$10^{-11}-10^{-9}$ for $C_2H_5CN$. The abundances of $CH_3OCHO$, $CH_3OCH_3$, and $C_2H_5CN$ are very strongly correlated (r>0.92) across ~4 orders of magnitude. $CH_3OCHO$ and $CH_3OCH_3$ show the strongest correlations in most parameters, and a nearly constant ratio (~1) over a remarkable ~9 orders of magnitude in luminosity for a wide variety of sources: pre-stellar to evolved cores, low- to high-mass objects, shocks, Galactic clouds, and comets. This indicates that COMs chemistry is likely early developed and then preserved through evolved phases. Moreover, the molecular abundances clearly increase with evolution. We consider $CH_3OCHO$ and $CH_3OCH_3$ to be most likely chemically linked: they could e.g. share a common precursor, or be formed one from the other. We propose a general scenario for all COMs, involving a formation in the cold, earliest phases of star formation and a following increasing desorption with the progressive heating of the evolving core.

We present the study of deuteration of cyanoacetylene (HC$_3$N) towards a sample of 28 high-mass star-forming cores divided into different evolutionary stages, from starless to evolved protostellar cores. We report for the first time the detection of DC$_3$N towards 15 high-mass cores. The abundance ratios of DC$_3$N with respect HC$_3$N range in the interval 0.003$-$0.022, lower than those found in low-mas protostars and dark clouds. No significant trend with the evolutionary stage, or with the kinetic temperature of the region, has been found. We compare the level of deuteration of HC$_3$N with those of other molecules towards the same sample, finding weak correlation with species formed only or predominantly in gas phase (N$_2$H$^+$ and HNC, respectively), and no correlation with species formed only or predominantly on dust grains (CH$_3$OH and NH$_3$, respectively). We also present a single-dish map of DC$_3$N towards the protocluster IRAS 05358+3543, which shows that DC$_3$N traces an extended envelope ($\sim$0.37 pc) and peaks towards two cold condensations separated from the positions of the protostars and the dust continuum. The observations presented in this work suggest that deuteration of HC$_3$N is produced in the gas of the cold outer parts of massive star-forming clumps, giving us an estimate of the deuteration factor prior to the formation of denser gas.

Aims: A unique environment to study how interstellar Complex Organic Molecules (iCOMs) can be formed is the shocked gas along low-mass protostellar outflows, as the dust mantles composition is sputtered into the gas phase. The chemical richness in these environments has been so far studied only in the L1157 blue shifted outflow. Methods: To understand if the L1157-B1 case is unique, we imaged the NGC 1333 IRAS 4A outflows using the NOEMA (NOrthern Extended Millimeter Array) interferometer as part of the IRAM SOLIS (Seeds Of Life in Space) Large Program and compared the observations with the GRAINOBLE+ gas phase astrochemical model. Results: Several iCOMs were detected in the IRAS 4A outflows: methanol (CH$_3$OH), acetaldehyde (CH$_3$CHO), formamide (NH$_2$CHO) and dimethyl ether (CH$_3$OCH$_3$), all sampling upper excitation energy up to $\sim$30 K. We found a significant chemical differentiation between the IRAS 4A1 outflow, showing a richer molecular content, and the IRAS 4A2 one. The CH$_3$OH/CH$_3$CHO abundance ratio is lower by a factor $\sim$4 in the former; furthermore the ratio in both outflows is lower by a factor $\sim$10 with respect to hot corinos values. Conclusions: After L1157-B1, IRAS 4A outflow is now the second outflow to show an evident chemical complexity. Given that CH$_3$OH is a grain surface species, GRAINOBLE+ reproduced our observations assuming acetaldehyde formation in gas phase by the reaction of ethyl radical (CH$_3$CH$_2$) with atomic oxygen. Moreover, the chemical differentiation between the two outflows suggests that the IRAS 4A1 outflow is likely younger than the IRAS 4A2 one. Further investigation is needed to constrain the age of the outflow and observations of even younger shocks are necessary and future spectroscopic studies on CH$_3$CH$_2$ are needed to be able to observe this species and provide strong constraints on the CH$_3$CHO formation.

C-fractionation has been studied from a theoretical point of view with different models of time-dependent chemistry, including both isotope-selective photodissociation and low-temperature isotopic exchange reactions. Recent chemical models predict that the latter may lead to a depletion of $^{13}$C in nitrile-bearing species, with $^{12}$C/$^{13}$C ratios two times higher than the elemental abundance ratio of 68 in the local ISM. Since the carbon isotopic ratio is commonly used to evaluate the $^{14}$N/$^{15}$N ratios with the double-isotope method, it is important to study C-fractionation in detail to avoid incorrect assumptions. In this work we implemented a gas-grain chemical model with new isotopic exchange reactions and investigated their introduction in the context of dense and cold molecular gas. In particular, we investigated the $^{12}$C/$^{13}$C ratios of HNC, HCN, and CN using a grid of models, with temperatures and densities ranging from 10 to 50 K and 2$\times$10$^{3}$ to 2$\times$10$^{7}$ cm$^{-3}$, respectively. We suggest a possible $^{13}$C exchange through the $^{13}$C + C$_{3}$ $\rightarrow$ $^{12}$C +$^{13}$CC$_{2}$ reaction, which does not result in dilution, but rather in $^{13}$C enhancement, for molecules formed starting from atomic carbon. This effect is efficient in a range of time between the formation of CO and its freeze-out on grains. Furthermore, we show that the $^{12}$C/$^{13}$C ratios of nitriles are predicted to be a factor 0.8-1.9 different from the local value of 68 for massive star-forming regions. This result also affects the $^{14}$N/$^{15}$N ratio: a value of 330 obtained with the double-isotope method is predicted to be 260-1150, depending on the physical conditions. Finally, we studied the $^{12}$C/$^{13}$C ratios by varying the cosmic-ray ionization rate: the ratios increase with it because of secondary photons and cosmic-ray reactions.

We present observations of SO and $\rm SO_2$ lines toward the shocked regions along the L1157 chemically rich outflow, taken in the context of the Seeds Of Life In Space IRAM-NOrthern Extended Millimeter Array Large Program, and supported by data from Submillimeter Array and IRAM-30 m telescope at 1.1--3.6 mm wavelengths. We simultaneously analyze, for the first time, all of the brightest shocks in the blueshifted lobe, namely, B0, B1, and B2. We found the following. (1) SO and $\rm SO_2$ may trace different gas, given that the large(-scale) velocity gradient analysis indicates for $\rm SO_2$ a volume density ($\rm 10^5\text{--}10^6\,cm^{-3}$) denser than that of the gas emitting in SO by a factor up to an order of magnitude. (2) Investigating the 0.1 pc scale field of view, we note a tentative gradient along the path of the precessing jet. More specifically, $\rm \chi({SO/SO_2})$ decreases from the B0-B1 shocks to the older B2. (3) At a linear resolution of 500--1400 au, a tentative spatial displacement between the two emitting molecules is detected, with the SO peak closer (with respect to $\rm SO_2$) to the position where the recent jet is impinging on the B1 cavity wall. Our astrochemical modeling shows that the SO and $\rm SO_2$ abundances evolve on timescales less than about 1000 years. Furthermore, the modeling requires high abundances ($2\times10^{-6}$) of both $\rm H_2S/H$ and S/H injected in the gas phase due to the shock occurrence, so pre-frozen OCS only is not enough to reproduce our new observations.

We show the first interferometric maps of the 14N/15N ratio obtained with the Atacama Large Millimeter Array (ALMA) towards the Solar-like forming protocluster OMC-2 FIR4. We observed N2H+, 15NNH+, N15NH+ (1-0), and N2D+ (2-1), from which we derive the isotopic ratios 14N/15N and D/H. The target is one of the closest analogues of the environment in which our Sun may have formed. The ALMA images, having synthesised beam corresponding to ~600 au, show that the emission of the less abundant isotopologues is distributed in several cores of ~10" (i.e. ~0.02~pc or 4000 au) embedded in a more extended N2H+ emission. Overall, our results indicate that: (1) 14N/15N does not change across the region at core scales, and (2) 14N/15N does not depend on temperature variations. Our findings also suggest that the 14N/15N variations found in pristine Solar System objects are likely not inherited from the protocluster stage, and hence their reason has to be found elsewhere.

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.

We present the Forgotten Quadrant Survey (FQS), an ESO large project that used the 12m antenna of the Arizona Radio Observatory to map the Galactic Plane in the range 220\deg$<l<$240\deg and -2.5\deg$<b<$0\deg, both in $^{12}$CO(1-0) and $^{13}$CO(1-0), at a spectral resolution of 0.65 km s$^{-1}$ and 0.26 km s$^{-1}$. Our dataset allows us to easily identify how the molecular dense gas is organised at different spatial scales: from the giant clouds with their denser filamentary networks, down to the clumps and cores that host the newborn stars and to obtain reliable estimates of their key physical parameters. We present the first release of the FQS data and discuss their quality. Spectra with 0.65 km s$^{-1}$ velocity channels have a noise ranging from 0.8 K to 1.3 K for $^{12}$CO(1-0) and from 0.3 K to 0.6 K for $^{13}$CO(1-0). In this paper, we used the $^{12}$CO(1-0) spectral cubes to produce a catalogue of 263 molecular clouds. This is the first selfconsistent, statistical catalogue of molecular clouds of the outer Galaxy, obtained with a subarcminute spatial resolution and therefore able to detect not only the classical giant molecular clouds, but also the small clouds and to resolve the cloud structure at the subparsec scale up to a distance of a few kpc. We found two classes of objects: structures with size above a few parsecs that are typical molecular clouds and may be self-gravitating, and subparsec structures that cannot be in gravitational equilibrium and are likely transient or confined by external pressure. We used the ratio between the Herschel H$_2$ column density and the integrated intensity of the CO lines to calculate the CO conversion factor and we found mean values of (3.3$\pm$1.4)$\times 10^{20}$ cm$^{-2}$(K km s$^{-1})^{-1}$ and (1.2$\pm$0.4)$\times 10^{21}$ cm$^{-2}$(K km s$^{-1})^{-1}$, for $^{12}$CO(1-0) and $^{13}$CO(1-0), respectively.

To understand how Phosphorus-bearing molecules are formed in star-forming regions, we have analysed ALMA observations of PN and PO towards the massive star-forming region AFGL 5142, combined with a new analysis of the data of the comet 67P/Churyumov-Gerasimenko taken with the ROSINA instrument onboard Rosetta. The ALMA maps show that the emission of PN and PO arises from several spots associated with low-velocity gas with narrow linewidths in the cavity walls of a bipolar outflow. PO is more abundant than PN in most of the spots, with the PO/PN ratio increasing as a function of the distance to the protostar. Our data favor a formation scenario in which shocks sputter phosphorus from the surface of dust grains, and gas-phase photochemistry induced by UV photons from the protostar allows efficient formation of the two species in the cavity walls. Our analysis of the ROSINA data has revealed that PO is the main carrier of P in the comet, with PO/PN>10. Since comets may have delivered a significant amount of prebiotic material to the early Earth, this finding suggests that PO could contribute significantly to the phosphorus reservoir during the dawn of our planet. There is evidence that PO was already in the cometary ices prior to the birth of the Sun, so the chemical budget of the comet might be inherited from the natal environment of the Solar System, which is thought to be a stellar cluster including also massive stars.

Phosphorus-bearing species are an essential key to form life on Earth, however they have barely been detected in the interstellar medium. Since only PN and PO have been identified so far towards star-forming regions, the chemical formation pathways of P-bearing molecules are not easy to constrain and are thus highly debatable. An important factor still missing in the chemical models is the initial elemental abundance of phosphorus, i.e. the depletion level of P. In order to overcome this problem, we study P-bearing species in diffuse/translucent clouds. In these objects phosphorus is mainly in the gas phase and therefore the elemental initial abundance needed in our chemical simulations corresponds to the cosmic one and is thus well constrained. An advanced chemical model with an updated P-chemistry network has been used. Single-pointing observations were performed with the IRAM 30m telescope towards the line of sight to the blazar B0355+508 aiming for the (2-1) transitions of PN, PO, HCP and CP. This line of sight incorporates five diffuse/translucent clouds. The (2-1) transitions of the PN, PO, HCP and CP were not detected. We report detections of the (1-0) lines of $\mathrm{^{13}CO}$, HNC and CN along with a first detection of $\mathrm{C^{34}S}$. We have reproduced the observations of HNC, CN, CS and CO in every cloud by applying typical conditions for diffuse/translucent clouds. According to our best-fit model, the most abundant P-bearing species are HCP and CP ($\sim10^{-10}$), followed by PN, PO and $\mathrm{PH_3}$ ($\sim10^{-11}$). We show that the production of P-bearing species is favoured towards translucent rather than diffuse clouds, where the environment provides a stronger shielding from the interstellar radiation. Based on our improved model, the (1-0) transitions of HCP, CP, PN and PO are expected to be detectable with estimated intensities up to $\sim200$ mK.

Observations of young stellar objects (YSOs) in centimeter bands can probe the continuum emission from growing dust grains, ionized winds, and magnetospheric activity, which are intimately connected to the evolution of protoplanetary disks and the formation of planets. We have carried out sensitive continuum observations toward the Ophiuchus A star-forming region using the Karl G. Jansky Very Large Array (VLA) at 10 GHz over a field-of-view of 6$'$ with a spatial resolution of $\theta_{maj}$ $\times$ $\theta_{min}$ $\sim$ 0.4$''$ $\times$ 0.2$''$. We achieved a 5 $\mu$Jy beam$^{-1}$ root-mean-square noise level at the center of our mosaic field of view. Among the eighteen sources we detected, sixteen are YSOs (three Class 0, five Class I, six Class II, and two Class III) and two are extragalactic candidates. We find that thermal dust emission generally contributes less that 30% of the emission at 10 GHz. The radio emission is dominated by other types of emission such as gyro-synchrotron radiation from active magnetospheres, free-free emission from thermal jets, free-free emission from the outflowing photo-evaporated disk material, and/or synchrotron emission from accelerated cosmic-rays in jet or protostellar surface shocks. These different types of emission could not be clearly disentangled. Our non-detections towards Class II/III disks suggest that extreme UV-driven photoevaporation is insufficient to explain the disk dispersal, assuming that the contribution of UV photoevaporating stellar winds to radio flux does not evolve with time. The sensitivity of our data cannot exclude photoevaporation due to X-ray photons as an efficient mechanism for disk dispersal. Deeper surveys with the Square Kilometre Array will be able to provide strong constraints on disk photoevaporation.

Phosphorus nitride (PN) is the P-bearing species with the highest number of detections in star-forming regions. Multi-line studies of the molecule have shown that the excitation temperature of PN is usually lower than the gas kinetic temperature, suggesting that PN is likely in conditions of sub-thermal excitation. We present an analysis of PN which takes the possible sub-thermal excitation conditions into account in a sample of 24 massive star-forming regions. We observed PN (2-1), (3-2), (4-3), and (6-5) with the IRAM-30m and APEX telescopes and detected PN lines in 15 of them. Together with 9 similar sources detected in PN in previous works, we have analysed the largest sample of star-forming regions to date, made of 33 sources with 24 detections in total (among which 13 are new detections). Hence, we have increased the number of star-forming regions detected in PN by more than a factor 2. Our analysis indicates that the PN lines are indeed sub-thermally excited, but well described by a single excitation temperature. We have compared line profiles and fractional abundances of PN and SiO, a typical shock tracer, and found that almost all objects detected in PN have high-velocity SiO wings. Moreover, the SiO and PN abundances with respect to H2 are correlated over several orders of magnitude, and uncorrelated with gas temperature. This clearly shows that the production of PN is strongly linked to the presence of shocked gas, and rules out alternative scenarios based on thermal evaporation from iced grain mantles.

Magneto-hydrodynamic (MHD) shocks are violent events that inject large amounts of energy in the interstellar medium (ISM) dramatically modifying its physical properties and chemical composition. Indirect evidence for the presence of such shocks has been reported from the especial chemistry detected toward a variety of astrophysical shocked environments. However, the internal physical structure of these shocks remains unresolved since their expected spatial scales are too small to be measured with current instrumentation. Here we report the first detection of a fully spatially resolved, MHD shock toward the Infrared Dark Cloud (IRDC) G034.77-00.55. The shock, probed by Silicon Monoxide (SiO) and observed with the Atacama Large Millimetre/sub-millimetre Array (ALMA), is associated with the collision between the dense molecular gas of the cloud and a molecular gas flow pushed toward the IRDC by the nearby supernova remnant (SNR) W44. The interaction is occurring on sub-parsec spatial scales thanks to the enhanced magnetic field of the SNR, making the dissipation region of the MHD shock large enough to be resolved with ALMA. Our observations suggest that molecular flow-flow collisions can be triggered by stellar feedback, inducing shocked molecular gas densities compatible with those required for massive star formation.

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.

We present ALMA CO(2-1) observations toward a massive infrared dark cloud G28.37+0.07. The ALMA data reveal numerous molecular (CO) outflows with a wide range of sizes throughout the cloud. Sixty-two 1.3 mm continuum cores were identified to be driving molecular outflows. We have determined the position angle in the plane-of-sky of 120 CO outflow lobes and studied their distribution. We find that the distribution of the plane-of-sky outflow position angles peaks at about 100 degree, corresponding to a concentration of outflows with an approximately east-west direction. For most outflows, we have been able to estimate the plane-of-sky angle between the outflow axis and the filament that harbors the protostar that powers the outflow. Statistical tests strongly indicate that the distribution of outflow-filament orientations is consistent with most outflow axes being mostly orthogonal to their parent filament in 3D. Such alignment may result from filament fragmentation or continuous mass transportation from filament to the embedded protostellar core. The latter is suggested by recent numerical studies with moderately strong magnetic fields.

Nitrogen (N) fractionation is used as a tool to search for a link between the chemical history of the Solar System and star-forming regions. A large variation of $^{14}$N/$^{15}$N is observed towards different astrophysical sources, and current chemical models cannot reproduce it. With the advent of high angular resolution radiotelescopes it is now possible to search for N-fractionation at core scales. We present IRAM NOEMA observations of the J=1-0 transition of N$_{2}$H$^{+}$, $^{15}$NNH$^{+}$ and N$^{15}$NNH$^{+}$ towards the high-mass protocluster IRAS 05358+3543. We find $^{14}$N/$^{15}$N ratios that span from $\sim$100 up to $\sim$220 and these values are lower or equal than those observed with single-dish observations towards the same source. Since N-fractionation changes across the studied region, this means that it is regulated by local environmental effects. We find also the possibility, for one of the four cores defined in the protocluster, to have a more abundant $^{15}$NNH$^{+}$ with respect to N$^{15}$NNH$^{+}$. This is another indication that current chemical models may be missing chemical reactions or may not take into account other mechanisms, like photodissociation or grain surface chemistry, that could be important.

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.

The formyl radical HCO has been proposed as the basic precursor of many complex organic molecules such as methanol (CH$_3$OH) or glycolaldehyde (CH$_2$OHCHO). Using ALMA, we have mapped, for the first time at high angular resolution ($\sim$1$^{\prime\prime}$, $\sim$140 au), HCO towards the Solar-type protostellar binary IRAS 16293$-$2422, where numerous complex organic molecules have been previously detected. We also detected several lines of the chemically related species H$_2$CO, CH$_3$OH and CH$_2$OHCHO. The observations revealed compact HCO emission arising from the two protostars. The line profiles also show redshifted absorption produced by foreground material of the circumbinary envelope that is infalling towards the protostars. Additionally, IRAM 30m single-dish data revealed a more extended HCO component arising from the common circumbinary envelope. The comparison between the observed molecular abundances and our chemical model suggests that whereas the extended HCO from the envelope can be formed via gas-phase reactions during the cold collapse of the natal core, the HCO in the hot corinos surrounding the protostars is predominantly formed by the hydrogenation of CO on the surface of dust grains and subsequent thermal desorption during the protostellar phase. The derived abundance of HCO in the dust grains is high enough to produce efficiently more complex species such as H$_2$CO, CH$_3$OH, and CH$_2$OHCHO by surface chemistry. We found that the main formation route of CH$_2$OHCHO is the reaction between HCO and CH$_2$OH.

Interstellar CO2 is an important reservoir of carbon and oxygen, and one of the major constituents of the icy mantles of dust grains, but it is not observable directly in the cold gas because has no permanent dipole moment. Its protonated form, HOCO+, is believed to be a good proxy for gaseous CO2. However, it has been detected in only a few star-forming regions so far, so that its interstellar chemistry is not well understood. We present new detections of HOCO+ lines in 11 high-mass star-forming clumps. Our observations increase by more than three times the number of detections in star-forming regions so far. We have derived beam-averaged abundances relative to H2 in between 0.3 and 3.8 x 10^{-11}. We have compared these values with the abundances of H13CO+, a possible gas-phase precursor of HOCO+, and CH3OH, a product of surface chemistry. We have found a positive correlation with H13CO+, while with CH3OH there is no correlation. We suggest that the gas-phase formation route starting from HCO+ plays an important role in the formation of HOCO+, perhaps more relevant than protonation of CO2 (upon evaporation of this latter from icy dust mantles).

Phosphorus is a crucial element in prebiotic chemistry, especially the P$-$O bond, which is key for the formation of the backbone of the deoxyribonucleic acid. So far, PO had only been detected towards the envelope of evolved stars, and never towards star-forming regions. We report the first detection of PO towards two massive star-forming regions, W51 e1/e2 and W3(OH), using data from the IRAM 30m telescope. PN has also been detected towards the two regions. The abundance ratio PO/PN is 1.8 and 3 for W51 and W3(OH), respectively. Our chemical model indicates that the two molecules are chemically related and are formed via gas-phase ion-molecule and neutral-neutral reactions during the cold collapse. The molecules freeze out onto grains at the end of the collapse and desorb during the warm-up phase once the temperature reaches $\sim$35 K. The observed molecular abundances of 10$^{-10}$ are predicted by the model if a relatively high initial abundance of phosphorus, 5$\times$10$^{-9}$, is assumed.

We present high resolution (0.2", 1000 AU) 1.3 mm ALMA observations of massive infrared dark cloud clump, G028.37+00.07-C1, thought to harbor the early stages of massive star formation. Using $\rm N_2D^+$(3-2) we resolve the previously identified C1-S core, separating the bulk of its emission from two nearby protostellar sources. C1-S is thus identified as a massive ($\sim50\:M_\odot$), compact ($\sim0.1\:$pc diameter) starless core, e.g., with no signs of outflow activity. Being highly deuterated, this is a promising candidate for a pre-stellar core on the verge of collapse. An analysis of its dynamical state indicates a sub-virial velocity dispersion compared to a trans-Alfv\'enic turbulent core model. However, virial equilibrium could be achieved with sub-Alfv\'enic conditions involving $\sim2\:$mG magnetic field strengths.

The formation of asteroids, comets and planets occurs in the interior of protoplanetary disks during the early phase of star formation. Consequently, the chemical composition of the disk might shape the properties of the emerging planetary system. In this context, it is crucial to understand whether and what organic molecules are synthesized in the disk. In this Letter, we report the first detection of formic acid (HCOOH) towards the TW Hydrae protoplanetary disk. The observations of the trans-HCOOH 6$_{(1,6)-5(1,5)}$ transition were carried out at 129~GHz with ALMA. We measured a disk-averaged gas-phase t-HCOOH column density of $\sim$ (2-4)$\times$10$^{12}$~cm$^{-2}$, namely as large as that of methanol. HCOOH is the first organic molecules containing two oxygen atoms detected in a protoplanetary disk, a proof that organic chemistry is very active even though difficult to observe in these objects. Specifically, this simplest acid stands as the basis for synthesis of more complex carboxylic acids used by life on Earth.

We report new interferometric images of cyclopropenylidene, c-C$_3$H$_2$, towards the young protocluster OMC-2 FIR\,4. The observations were performed at 82 and 85 GHz with the NOrthern Extended Millimeter Array (NOEMA) as part of the project Seeds Of Life In Space (SOLIS). In addition, IRAM-30m data observations were used to investigate the physical structure of OMC-2 FIR\,4. We find that the c-C$_3$H$_2$ gas emits from the same region where previous SOLIS observations showed bright HC$_5$N emission. From a non-LTE analysis of the IRAM-30m data, the c-C$_3$H$_2$ gas has an average temperature of $\sim$40K, a H$_2$ density of $\sim$3$\times$10$^{5}$~cm$^{-3}$, and a c-C$_3$H$_2$ abundance relative to H$_2$ of ($7\pm1$)$\times$10$^{-12}$. In addition, the NOEMA observations provide no sign of significant c-C$_3$H$_2$ excitation temperature gradients across the region (about 3-4 beams), with T$_{ex}$ in the range 8$\pm$3 up to 16$\pm$7K. We thus infer that our observations are inconsistent with a physical interaction of the OMC-2 FIR\,4 envelope with the outflow arising from OMC-2 FIR\,3, as claimed by previous studies. The comparison of the measured c-C$_3$H$_2$ abundance with the predictions from an astrochemical PDR model indicates that OMC-2 FIR\,4 is irradiated by a FUV field $\sim$1000 times larger than the interstellar one, and by a flux of ionising particles $\sim$4000 times larger than the canonical value of $1\times10^{-17}$~s$^{-1}$ from the Galaxy cosmic rays, which is consistent with our previous HC$_5$N observations. This provides an important and independent confirmation of other studies that one or more sources inside the OMC-2 FIR\,4 region emit energetic ($\geq10$~MeV) particles.

The fractionation of nitrogen (N) in star-forming regions is a poorly understood process. To put more stringent observational constraints on the N-fractionation, we have observed with the IRAM-30m telescope a large sample of 66 cores in massive star-forming regions. We targeted the (1-0) rotational transition of HN$^{13}$C, HC$^{15}$N, H$^{13}$CN and HC$^{15}$N, and derived the $^{14}$N/$^{15}$N ratio for both HCN and HNC. We have completed this sample with that already observed by Colzi et al. (2018), and thus analysed a total sample of 87 sources. The $^{14}$N/$^{15}$N ratios are distributed around the Proto-Solar Nebula value with a lower limit near the terrestrial atmosphere value ($\sim$272). We have also derived the $^{14}$N/$^{15}$N ratio as a function of the Galactocentric distance and deduced a linear trend based on unprecedented statistics. The Galactocentric dependences that we have found are consistent, in the slope, with past works but we have found a new local $^{14}$N/$^{15}$N value of $\sim$400, i.e. closer to the Prosolar Nebula value. A second analysis was done, and a parabolic Galactocentric trend was found. Comparison with Galactic chemical evolution models shows that the slope until 8 kpc is consistent with the linear analysis, while the flattening trend above 8 kpc is well reproduced by the parabolic analysis.

Fragmentation of massive dense molecular clouds is the starting point in the formation of rich clusters and massive stars. Theory and numerical simulations indicate that the population of the fragments (number, mass, diameter, separation) resulting from the gravitational collapse of such clumps is probably regulated by the balance between the magnetic field and the other competitors of self-gravity, in particular turbulence and protostellar feedback. We have observed 11 massive, dense and young star-forming clumps with the Atacama Large Millimeter Array (ALMA) in the thermal dust continuum emission at $\sim 1$~mm with an angular resolution of 0.25 arcseconds with the aim of determining their population of fragments. We find fragments on sub-arcsecond scales in 8 out of the 11 sources. The ALMA images indicate two different fragmentation modes: a dominant fragment surrounded by companions with much smaller mass and size, and many ($\geq 8$) fragments with a gradual change in masses and sizes. On average, the largest number of fragments is found towards the warmer and more massive clumps. Also, the warmer clumps tend to form fragments with larger mass and size. To understand the role of the different physical parameters to regulate the final population of the fragments, we have simulated the collapse of a massive clump of $100$ and $300$ M$_{\odot}$ having different magnetic support. The simulations indicate that: (1) fragmentation is inhibited when the initial turbulence is low, independent of the other physical parameters. (2) a filamentary distribution of the fragments is favoured in a highly magnetised clump. We conclude that the clumps that show many fragments distributed in a filamentary-like structure are likely characterised by a strong magnetic field, while the others are possible also in a weaker magnetic field.

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.

Stars are born from dense cores in molecular clouds. Observationally, it is crucial to capture the formation of cores in order to understand the necessary conditions and rate of the star formation process. The {\it Atacama Large Mm/sub-mm Array} (ALMA) is extremely powerful for identifying dense gas structures, including cores, at mm wavelengths via their dust continuum emission. Here we use ALMA to carry out a survey of dense gas and cores in the central region of the massive ($\sim10^5\:M_\odot$) Infrared Dark Cloud (IRDC) G28.37+0.07. The observation consists of a mosaic of 86 pointings of the 12m-array and produces an unprecedented view of the densest structures of this IRDC. In this first paper about this data set, we focus on a comparison between the 1.3 mm continuum emission and a mid-infrared (MIR) extinction map of the IRDC. This allows estimation of the "dense gas" detection probability function (DPF), i.e., as a function of the local mass surface density, $\Sigma$, for various choices of thresholds of mm continuum emission to define "dense gas". We then estimate the dense gas mass fraction, $f_{\rm dg}$, in the central region of the IRDC and, via extrapolation with the DPF and the known $\Sigma$ probability distribution function, to the larger-scale surrounding regions, finding values of about 5\% to 15\% for the fiducial choice of threshold. We argue that this observed dense gas is a good tracer of the protostellar core population and, in this context, estimate a star formation efficiency per free-fall time in the central IRDC region of $\epsilon_{\rm ff}\sim$10\%, with approximately a factor of two systematic uncertainties.

We present multi-transition observations of PN towards a sample of nine massive dense cores in different evolutionary stages. Using transitions with different excitation conditions, we have found for the first time that the excitation temperatures of PN are in the range 5-30 K. To investigate the main chemical route for the PN formation (surface-chemistry vs. gas-phase chemistry), and the dominant desorption mechanism (thermal vs. shock), we have compared our results with those obtained from molecules tracing different chemical and physical conditions (SiO, SO, CH3OH, and N2H+). We have found that the PN line profiles are very well correlated with those of SiO and SO in six out of the nine targets, which indicates that PN may be released by sputtering of dust grains due to shocks. This finding is corroborated by a faint but statistically significant positive trend between the PN abundance and those of SiO and SO. However, in three objects the PN lines have no hints of high velocity wings, which indicates an alternative origin of PN. Overall, our results indicate that the origin of PN is not unique, as it can be formed in protostellar shocks, but also in colder and more quiescent gas through alternative pathways.

Towards the pre-stellar core L1544, the methanol (CH$_3$OH) emission forms an asymmetric ring around the core centre, where CH$_3$OH is mostly in solid form, with a clear peak 4000~au to the north-east of the dust continuum peak. As part of the NOEMA Large Project SOLIS (Seeds of Life in Space), the CH$_3$OH peak has been spatially resolved to study its kinematics and physical structure and to investigate the cause behind the local enhancement. We find that methanol emission is distributed in a ridge parallel to the main axis of the dense core. The centroid velocity increases by about 0.2~km~s$^{-1}$ and the velocity dispersion increases from subsonic to transonic towards the central zone of the core, where the velocity field also shows complex structure. This could be indication of gentle accretion of material onto the core or interaction of two filaments, producing a slow shock. We measure the rotational temperature and show that methanol is in local thermodynamic equilibrium (LTE) only close to the dust peak, where it is significantly depleted. The CH$_3$OH column density, $N_{tot}({\rm CH_3OH})$, profile has been derived with non-LTE radiative transfer modelling and compared with chemical models of a static core. The measured $N_{tot}({\rm CH_3OH})$ profile is consistent with model predictions, but the total column densities are one order of magnitude lower than those predicted by models, suggesting that the efficiency of reactive desorption or atomic hydrogen tunnelling adopted in the model may be overestimated; or that an evolutionary model is needed to better reproduce methanol abundance.

Although chemical models predict that the deuterium fractionation in N$_2$H$^+$ is a good evolutionary tracer in the star formation process, the fractionation of nitrogen is still a poorly understood process. Recent models have questioned the similar evolutionary trend expected for the two fractionation mechanisms in N$_2$H$^+$, based on a classical scenario in which ion-neutral reactions occurring in cold gas should have caused an enhancement of the abundance of N$_2$D$^+$, $^{15}$NNH$^+$, and N$^{15}$NH$^+$. In the framework of the ASAI IRAM-30m large program, we have investigated the fractionation of deuterium and $^{15}$N in N$_2$H$^+$ in the best known representatives of the different evolutionary stages of the Sun-like star formation process. The goal is to ultimately confirm (or deny) the classical "ion-neutral reactions" scenario that predicts a similar trend for D and $^{15}$N fractionation. We do not find any evolutionary trend of the $^{14}$N/$^{15}$N ratio from both the $^{15}$NNH$^+$ and N$^{15}$NH$^+$ isotopologues. Therefore, our findings confirm that, during the formation of a Sun-like star, the core evolution is irrelevant in the fractionation of $^{15}$N. The independence of the $^{14}$N/$^{15}$N ratio with time, found also in high-mass star-forming cores, indicates that the enrichment in $^{15}$N revealed in comets and protoplanetary disks is unlikely to happen at core scales. Nevertheless, we have firmly confirmed the evolutionary trend expected for the H/D ratio, with the N$_2$H$^+$/N$_2$D$^+$ ratio decreasing before the pre-stellar core phase, and increasing monotonically during the protostellar phase. We have also confirmed clearly that the two fractionation mechanisms are not related.

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.

Phosphorus (P) is one of the essential elements for life due to its central role in biochemical processes. Recent searches have shown that P-bearing molecules (in particular PN and PO) are present in star-forming regions, although their formation routes remain poorly understood. In this Letter, we report observations of PN and PO towards seven molecular clouds located in the Galactic Center, which are characterized by different types of chemistry. PN is detected in five out of seven sources, whose chemistry is thought to be shock-dominated. The two sources with PN non-detections correspond to clouds exposed to intense UV/X-rays/cosmic-ray radiation. PO is detected only towards the cloud G+0.693$-$0.03, with a PO/PN abundance ratio of $\sim$1.5. We conclude that P-bearing molecules likely form in shocked gas as a result of dust grain sputtering, while are destroyed by intense UV/X-ray/cosmic ray radiation.

We present a complete census of the 14N/15N isotopic ratio in the most abundant N-bearing molecules towards the cold envelope of the protocluster OMC-2 FIR4, the best known Sun progenitor. To this scope, we analysed the unbiased spectral survey obtained with the IRAM-30m telescope at 3mm, 2mm and 1mm. We detected several lines of CN, HCN, HNC, HC3N, N2H+, and their respective 13C and 15N isotopologues. The lines relative fluxes are compatible with LTE conditions and moderate line opacities have been corrected via a Population Diagram method or theoretical relative intensity ratios of the hyperfine structures. The five species lead to very similar 14N/15N isotopic ratios, without any systematic difference between amine and nitrile bearing species as previously found in other protostellar sources. The weighted average of the 14N/15N isotopic ratio is 270 +/- 30. This 14N/15N value is remarkably consistent with the [250-350] range measured for the local galactic ratio but significantly differs from the ratio measured in comets (around 140). High-angular resolution observations are needed to examine whether this discrepancy is maintained at smaller scales. In addition, using the CN, HCN and HC3N lines, we derived a 12C/13C isotopic ratio of 50 +/- 5.

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.

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 ratio between the two stable isotopes of nitrogen, $^{14}$N and $^{15}$N, is well measured in the terrestrial atmosphere ($\sim 272$), and in the pre-Solar nebula ($\sim 441$). Interestingly, some pristine Solar System materials show enrichments in $^{15}$N with respect to the pre-Solar nebula value. However, it is not yet clear if, and how, these enrichments are linked to the past chemical history, due to the limited number of measurements in dense star-forming regions. In this respect, dense cores believed to be precursors of clusters containing also intermediate- and high-mass stars are important targets, as the Solar System was probably born within a rich stellar cluster. In this work, we show the results of IRAM-30m observations of the J=1-0 rotational transition of the molecules HCN and HNC, and their $^{15}$N-bearing counterparts, towards 27 intermediate/high-mass dense cores divided in three evolutionary categories: high-mass starless cores, high-mass protostellar objects, and ultra-compact HII regions. We have also observed the DNC(2-1) rotational transition, in order to search for a relation between the isotopic ratios D/H and $^{14}$N/$^{15}$N. We derive average $^{14}$N/$^{15}$N ratios of $359\pm16$ in HCN and of $438\pm21$ in HNC, with a dispersion of about 150-200. We find no trend of the $^{14}$N/$^{15}$N ratio with the evolutionary stage. This result agrees with what found from N$_{2}$H$^{+}$ and its isotopologues in the same sources, although the $^{14}$N/$^{15}$N ratios from N$_{2}$H$^{+}$ show a dispersion larger than that in HCN/HNC. Moreover, we have found no correlation between D/H and $^{14}$N/$^{15}$N in HNC. These findings indicate that: (1) the chemical evolution does not seem to play a role in the fractionation of nitrogen; (2) the fractionation of hydrogen and nitrogen in these objects are not related.

Context: Modern versions of the Miller-Urey experiment claim that formamide (NH$_2$CHO) could be the starting point for the formation of metabolic and genetic macromolecules. Intriguingly, formamide is indeed observed in regions forming Solar-type stars as well as in external galaxies. Aims: How NH$_2$CHO is formed has been a puzzle for decades: our goal is to contribute to the hotly debated question of whether formamide is mostly formed via gas-phase or grain surface chemistry. Methods: We used the NOEMA interferometer to image NH$_2$CHO towards the L1157-B1 blue-shifted shock, a well known interstellar laboratory, to study how the components of dust mantles and cores released into the gas phase triggers the formation of formamide. Results: We report the first spatially resolved image (size $\sim$ 9", $\sim$ 2300 AU) of formamide emission in a shocked region around a Sun-like protostar: the line profiles are blueshifted and have a FWHM $\simeq$ 5 km s$^{-1}$. A column density of $N_{\rm NH_2CHO}$ = 8 $\times$ 10$^{12}$ cm$^{-1}$, and an abundance (with respect to H-nuclei) of 4 $\times$ 10$^{-9}$ are derived. We show a spatial segregation of formamide with respect to other organic species. Our observations, coupled with a chemical modelling analysis, indicate that the formamide observed in L1157-B1 is formed by gas-phase chemical process, and not on grain surfaces as previously suggested. Conclusions: The SOLIS interferometric observations of formamide provide direct evidence that this potentially crucial brick of life is efficiently formed in the gas-phase around Sun-like protostars.

L1157-B1 is the brightest shocked region of the large-scale molecular outflow, considered the prototype of chemically rich outflows, being the ideal laboratory to study how shocks affect the molecular gas. Several deuterated molecules have been previously detected with the IRAM 30m, most of them formed on grain mantles and then released into the gas phase due to the shock. We aim to observationally investigate the role of the different chemical processes at work that lead to formation the of DCN and test the predictions of the chemical models for its formation. We performed high-angular resolution observations with NOEMA of the DCN(2-1) and H13CN(2-1) lines to compute the deuterated fraction, Dfrac(HCN). We detected emission of DCN(2-1) and H13CN(2-1) arising from L1157-B1 shock. Dfrac(HCN) is ~4x10$^{-3}$ and given the uncertainties, we did not find significant variations across the bow-shock. Contrary to HDCO, whose emission delineates the region of impact between the jet and the ambient material, DCN is more widespread and not limited to the impact region. This is consistent with the idea that gas-phase chemistry is playing a major role in the deuteration of HCN in the head of the bow-shock, where HDCO is undetected as it is a product of grain-surface chemistry. The spectra of DCN and H13CN match the spectral signature of the outflow cavity walls, suggesting that their emission result from shocked gas. The analysis of the time dependent gas-grain chemical model UCL-CHEM coupled with a C-type shock model shows that the observed Dfrac(HCN) is reached during the post-shock phase, matching the dynamical timescale of the shock. Our results indicate that the presence of DCN in L1157-B1 is a combination of gas-phase chemistry that produces the widespread DCN emission, dominating in the head of the bow-shock, and sputtering from grain mantles toward the jet impact region.

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.

The shock L1157-B1 driven by the low-mass protostar L1157-mm is an unique environment to investigate the chemical enrichment due to molecules released from dust grains. IRAM-30m and Plateau de Bure Interferometer observations allow a census of Si-bearing molecules in L1157-B1. We detect SiO and its isotopologues and, for the first time in a shock, SiS. The strong gradient of the [SiO/SiS] abundance ratio across the shock (from >=180 to ~25) points to a different chemical origin of the two species. SiO peaks where the jet impacts the cavity walls ([SiO/H2] ~ 1e-6), indicating that SiO is directly released from grains or rapidly formed from released Si in the strong shock occurring at this location. In contrast, SiS is only detected at the head of the cavity opened by previous ejection events ([SiS/H2] ~ 2e-8). This suggests that SiS is not directly released from the grain cores but instead should be formed through slow gas-phase processes using part of the released silicon. This finding shows that Si-bearing molecules can be useful to distinguish regions where grains or gas-phase chemistry dominates.

We present the results of formaldehyde and methanol deuteration measurements towards the Class I low-mass protostar SVS13-A, in the framework of the IRAM 30-m ASAI (Astrochemical Surveys At IRAM) project. We detected emission lines of formaldehyde, methanol, and their deuterated forms (HDCO, D2CO, CHD2OH, CH3OD) with Eup up to 276 K. The formaldehyde analysis indicates Tkin = 15 - 30 K, n (H2) >= 10^6 cm^-3, and a size of about 1200 AU suggesting an origin in the protostellar envelope. For methanol we find two components: (i) a high temperature (Tkin = 80 K) and very dense (> 10^8 cm^-3}) gas from a hot corino (radius about 35 AU), and (ii) a colder Tkin <= 70 K) and more extended (radius about 350 AU) region. The deuterium fractionation is 9 10^-2 for HDCO, 4 10^-3 for D2CO, and 2 - 7 10^-3 for CH2DOH, up to two orders of magnitude lower than the values measured in Class 0 sources. We derive also formaldehyde deuteration in the outflow: 4 10^-3, in agreement with what found in the L1157-B1 protostellar shock. Finally, we estimate [CH2DOH]/[CH3OD] about 2. The decrease of deuteration in the Class I source SVS13-A with respect to Class 0 sources can be explained by gas-phase processes. Alternatively, a lower deuteration could be the effect of a gradual collapse of less deuterated external shells of the protostellar evelope. The present measurements fill in the gap between prestellar cores and protoplanetary disks in the context of organics deuteration measurements.

We carry out an ALMA $\rm N_2D^+$(3-2) and 1.3~mm continuum survey towards 32 high mass surface density regions in seven Infrared Dark Clouds with the aim of finding massive starless cores, which may be the initial conditions for the formation of massive stars. Cores showing strong $\rm N_2D^+$(3-2) emission are expected to be highly deuterated and indicative of early, potentially pre-stellar stages of star formation. We also present maps of these regions in ancillary line tracers, including C$^{18}$O(2-1), DCN(3-2) and DCO$^+$(3-2). Over 100 $\rm N_2D^+$ cores are identified with our newly developed core-finding algorithm based on connected structures in position-velocity space. The most massive core has $\sim70\:M_\odot$ (potentially $\sim170\:M_\odot$) and so may be representative of the initial conditions for massive star formation. The existence and dynamical properties of such cores constrain massive star formation theories. We measure the line widths and thus velocity dispersion of six of the cores with strongest $\rm N_2D^+$(3-2) line emission, finding results that are generally consistent with virial equilibrium of pressure confined cores.

The detection of organic molecules with increasing complexity and potential biological relevance is opening the possibility to understand the formation of the building blocks of life in the interstellar medium. One of the families of molecules with astrobiological interest are the esters, whose simplest member, methyl formate, is rather abundant in star-forming regions. The next step in the chemical complexity of esters is ethyl formate, C$_2$H$_5$OCHO. Only two detections of this species have been reported so far, which strongly limits our understanding of how complex molecules are formed in the interstellar medium. We have searched for ethyl formate towards the W51 e2 hot molecular core, one of the most chemically rich sources in the Galaxy and one of the most promising regions to study prebiotic chemistry, especially after the recent discovery of the P$-$O bond, key in the formation of DNA. We have analyzed a spectral line survey towards the W51 e2 hot molecular core, which covers 44 GHz in the 1, 2 and 3 mm bands, carried out with the IRAM 30m telescope. We report the detection of the trans and gauche conformers of ethyl formate. A Local Thermodynamic Equilibrium analysis indicates that the excitation temperature is 78$\pm$10 K and that the two conformers have similar source-averaged column densities of (2.0$\pm$0.3)$\times$10$^{16}$ cm$^{-2}$ and an abundance of $\sim$10$^{-8}$. We compare the observed molecular abundances of ethyl formate with different competing chemical models based on grain surface and gas-phase chemistry. We propose that grain-surface chemistry may have a dominant role in the formation of ethyl formate (and other complex organic molecules) in hot molecular cores, rather than reactions in the gas phase.

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.

Massive stars, multiple stellar systems and clusters are born from the gravitational collapse of massive dense gaseous clumps, and the way these systems form strongly depends on how the parent clump fragments into cores during collapse. Numerical simulations show that magnetic fields may be the key ingredient in regulating fragmentation. Here we present ALMA observations at ~0.25'' resolution of the thermal dust continuum emission at ~278 GHz towards a turbulent, dense, and massive clump, IRAS 16061-5048c1, in a very early evolutionary stage. The ALMA image shows that the clump has fragmented into many cores along a filamentary structure. We find that the number, the total mass and the spatial distribution of the fragments are consistent with fragmentation dominated by a strong magnetic field. Our observations support the theoretical prediction that the magnetic field plays a dominant role in the fragmentation process of massive turbulent clump.

We study the molecular abundance and spatial distribution of the simplest sugar alcohol, ethylene glycol (EG), the simplest sugar glycoladehyde (GA), and other chemically related complex organic species towards the massive star-forming region G31.41+0.31. We have analyzed multiple single dish and interferometric data, and obtained excitation temperatures and column densities using an LTE analysis. We have reported for the first time the presence of EG towards G31.41+0.31, and we have also detected multiple transitions of other complex organic molecules such as GA, methyl formate (MF), dimethyl ether (DME) and ethanol (ET). The high angular resolution images show that the EG emission is very compact, peaking towards the maximum of the continuum. These observations suggest that low abundance complex organic molecules, like EG or GA, are good probes of the gas located closer to the forming stars. Our analysis confirms that EG is more abundant than GA in G31.41+0.31, as previously observed in other interstellar regions. Comparing different star-forming regions we find evidence of an increase of the EG/GA abundance ratio with the luminosity of the source. The DME/MF and EG/ET ratios are nearly constant with luminosity. We have also found that the abundance ratios of pairs of isomers GA/MF and ET/DME decrease with the luminosity of the sources. The most likely explanation for the behavior of the EG/GA ratio is that these molecules are formed by different chemical formation routes not directly linked; although warm-up timescales effects and different formation and destruction efficiencies in the gas phase cannot be ruled out. The most likely formation route of EG is by combination of two CH$_{2}$OH radicals on dust grains. We also favor that GA is formed via the solid-phase dimerization of the formyl radical HCO, and a chemical link between MF and DME.

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.

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

In the context of the ASAI (Astrochemical Surveys At IRAM) project, we carried out an unbiased spectral survey in the millimeter window towards the well known low-mass Class I source SVS13-A. The high sensitivity reached (3-12 mK) allowed us to detect at least 6 HDO broad (FWHM ~ 4-5 km/s) emission lines with upper level energies up to Eu = 837 K. A non-LTE LVG analysis implies the presence of very hot (150-260 K) and dense (> 3 10^7 cm-3) gas inside a small radius ($\sim$ 25 AU) around the star, supporting, for the first time, the occurrence of a hot corino around a Class I protostar. The temperature is higher than expected for water molecules are sublimated from the icy dust mantles (~ 100 K). Although we cannot exclude we are observig the effects of shocks and/or winds at such small scales, this could imply that the observed HDO emission is tracing the water abundance jump expected at temperatures ~ 220-250 K, when the activation barrier of the gas phase reactions leading to the formation of water can be overcome. We derive X(HDO) ~ 3 10-6, and a H2O deuteration > 1.5 10-2, suggesting that water deuteration does not decrease as the protostar evolves from the Class 0 to the Class I stage.

Phosphorus is a crucial element in biochemistry, especially the P-O bond, which is key for the formation of the backbone of the deoxyribonucleic acid. So far, PO has only been detected towards the envelope of evolved stars, and never towards star-forming regions. We report the first detections of PO towards two massive star-forming regions, W51 e1/e2 and W3(OH), using data from the IRAM 30m telescope. PN has also been detected towards the two regions. The abundance ratio PO/PN is 1.8 and 3 for W51 and W3(OH), respectively. Our chemical model indicates that the two molecules are chemically related and are formed via gas-phase ion-molecule and neutral-neutral reactions during the cold collapse. The molecules freeze out onto grains at the end of the collapse and desorb during the warm-up phase once the temperature reaches 35 K. Similar abundances of the two species are expected during a period of 5x10^{4} yr at the early stages of the warm-up phase, when the temperature is in the range 35-90 K. The observed molecular abundances of 10^{-10} are predicted by the model if a relatively high initial abundance of 5x10^{-9} of depleted phosphorus is assumed.

We study deuterium fractionation in two massive starless/early-stage cores C1-N and C1-S in Infrared Dark Cloud (IRDC) G028.37+00.07, first identified by Tan et al. (2013) with ALMA. Line emission from multiple transitions of $\rm N_2H^+$ and $\rm N_2D^+$ were observed with the ALMA, CARMA, SMA, JCMT, NRO 45m and IRAM 30m telescopes. By simultaneously fitting the spectra, we estimate the excitation conditions and deuterium fraction, $D_{\rm frac}^{\rm N_2H^+} \equiv [\rm N_2D^+]/[N_2H^+]$, with values of $D_{\rm frac}^{\rm N_2H^+} \simeq 0.2$--$0.7$, several orders of magnitude above the cosmic [D]/[H] ratio. Additional observations of o-H$_2$D$^+$ are also presented that help constrain the ortho-to-para ratio of $\rm H_2$, which is a key quantity affecting the degree of deuteration. We then present chemodynamical modeling of the two cores, exploring especially the implications for the collapse rate relative to free-fall, $\alpha_{\rm ff}$. In order to reach the high level of observed deuteration of $\rm N_2H^+$, we find that the most likely evolutionary history of the cores involves collapse at a relatively slow rate, $\lesssim1/10$th of free-fall.

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

We present ALMA follow-up observations of two massive, early-stage core candidates, C1-N & C1-S, in Infrared Dark Cloud (IRDC) G028.37+00.07, which were previously identified by their N2D+(3-2) emission and show high levels of deuteration of this species. The cores are also dark at far infrared wavelengths up to ~100 microns. We detect 12CO(2-1) from a narrow, highly-collimated bipolar outflow that is being launched from near the center of the C1-S core, which is also the location of the peak 1.3mm dust continuum emission. This protostar, C1-Sa, has associated dense gas traced by C18O(2-1) and DCN(3-2), from which we estimate it has a radial velocity that is near the center of the range exhibited by the C1-S massive core. A second outflow-driving source is also detected within the projected boundary of C1-S, but appears to be at a different radial velocity. After considering properties of the outflows, we conclude C1-Sa is a promising candidate for an early-stage massive protostar and as such it shows that these early phases of massive star formation can involve highly ordered outflow, and thus accretion, processes, similar to models developed to explain low-mass protostars.

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.

Aims: Using the unprecedented combination of high resolution and sensitivity offered by ALMA, we aim to investigate whether and how hot corinos, circumstellar disks, and ejected gas are related in young solar-mass protostars. Methods: We observed CH$_3$CHO and deuterated water (HDO) high-excitation ($E_{\rm u}$ up to 335 K) lines towards the Sun-like protostar HH212--MM1. Results: For the first time, we have obtained images of CH$_3$CHO and HDO emission in the inner $\simeq$ 100 AU of HH212. The multifrequency line analysis allows us to contrain the density ($\geq$ 10$^{7}$ cm$^{-3}$), temperature ($\simeq$ 100 K), and CH$_3$CHO abundance ($\simeq$ 0.2--2 $\times$ 10$^{-9}$) of the emitting region. The HDO profile is asymmetric at low velocities ($\leq$ 2 km s$^{-1}$ from $V_{\rm sys}$). If the HDO line is optically thick, this points to an extremely small ($\sim$ 20--40 AU) and dense ($\ge$ 10$^{9}$ cm$^{-3}$) emitting region. Conclusions: We report the first detection of a hot corino in Orion. The HDO asymmetric profile indicates a contribution of outflowing gas from the compact central region, possibly associated with a dense disk wind.

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

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

To understand massive star formation requires study of its initial conditions. Two massive starless core candidates, C1-N & C1-S, have been detected in IRDC G028.37+00.07 in $\rm N_2D^+$(3-2) with $ALMA$. From their line widths, either the cores are subvirial and are thus young structures on the verge of near free-fall collapse, or they are threaded by $\sim1$ mG $B$-fields that help support them in near virial equilibrium and potentially have older ages. We modeled the deuteration rate of $\rm N_2H^+$ to constrain collapse rates of the cores. First, to measure their current deuterium fraction, $D_{\rm frac}^{\rm N_2H^+}$ $\equiv [\rm N_2D^+]/[N_2H^+]$, we observed multiple transitions of $\rm N_2H^+$ and $\rm N_2D^+$ with $CARMA$, $SMA$, $JCMT$, $NRO~45m$ and $IRAM~30m$, to complement the $ALMA$ data. For both cores we derived $D_{\rm frac}^{\rm N_2H^+}\sim0.3$, several orders of magnitude above the cosmic [D]/[H] ratio. We then carried out chemodynamical modeling, exploring how collapse rate relative to free-fall, $\alpha_{\rm ff}$, affects the level of $D_{\rm frac}^{\rm N_2H^+}$ that is achieved from a given initial condition. To reach the observed $D_{\rm frac}^{\rm N_2H^+}$, most models require slow collapse with $\alpha_{\rm ff}\sim0.1$, i.e., $\sim1/10$th of free-fall. This makes it more likely that the cores have been able to reach a near virial equilibrium state and we predict that strong $B$-fields will eventually be detected. The methods developed here will be useful for measurement of the pre-stellar core mass function.

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

The Premiale Project "Science and Technology in Italy for the upgraded ALMA Observatory - iALMA" has the goal of strengthening the scientific, technological and industrial Italian contribution to the Atacama Large Millimeter/submillimeter Array (ALMA), the largest ground based international infrastructure for the study of the Universe in the microwave. One of the main objectives of the Science Working Group (SWG) inside iALMA, the Work Package 1, is to develop the Italian contribution to the Science Case for the ALMA Band 2 or Band 2+3 receiver. ALMA Band 2 receiver spans from ~67 GHz (bounded by an opaque line complex of ozone lines) up to 90 GHz which overlaps with the lower frequency end of ALMA Band 3. Receiver technology has advanced since the original definition of the ALMA frequency bands. It is now feasible to produce a single receiver which could cover the whole frequency range from 67 GHz to 116 GHz, encompassing Band 2 and Band 3 in a single receiver cartridge, a so called Band 2+3 system. In addition, upgrades of the ALMA system are now foreseen that should double the bandwidth to 16 GHz. The science drivers discussed below therefore also discuss the advantages of these two enhancements over the originally foreseen Band 2 system.

This booklet contains a collection of contributions to the meeting of the JEts and Disks at INAF (JEDI) group, which took place at the Capodimonte Observatory during 9-10 April 2015. Scope of the meeting was to bring together the JEDI researchers of the Italian Istituto Nazionale di Astrofisica (INAF) working in the field of circumstellar disks and jets in young stars, to discuss together the different agents affecting the structure and the evolution of disks, namely accretion, jets and winds. More information on the JEDI group and its activities can be found at \texttt{http://www.oa-roma.inaf.it/irgroup/JEDI}.

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

Infrared dark clouds (IRDCs) are dense, molecular structures in the interstellar medium that can harbour sites of high-mass star formation. IRDCs contain supersonic turbulence, which is expected to generate shocks that locally heat pockets of gas within the clouds. We present observations of the CO J = 8-7, 9-8, and 10-9 transitions, taken with the Herschel Space Observatory, towards four dense, starless clumps within IRDCs (C1 in G028.37+00.07, F1 and F2 in G034.43+0007, and G2 in G034.77-0.55). We detect the CO J = 8-7 and 9-8 transitions towards three of the clumps (C1, F1, and F2) at intensity levels greater than expected from photodissociation region (PDR) models. The average ratio of the 8-7 to 9-8 lines is also found to be between 1.6 and 2.6 in the three clumps with detections, significantly smaller than expected from PDR models. These low line ratios and large line intensities strongly suggest that the C1, F1, and F2 clumps contain a hot gas component not accounted for by standard PDR models. Such a hot gas component could be generated by turbulence dissipating in low velocity shocks.

A major goal of the Atacama Large Millimeter/submillimeter Array (ALMA) is to make accurate images with resolutions of tens of milliarcseconds, which at submillimeter (submm) wavelengths requires baselines up to ~15 km. To develop and test this capability, a Long Baseline Campaign (LBC) was carried out from September to late November 2014, culminating in end-to-end observations, calibrations, and imaging of selected Science Verification (SV) targets. This paper presents an overview of the campaign and its main results, including an investigation of the short-term coherence properties and systematic phase errors over the long baselines at the ALMA site, a summary of the SV targets and observations, and recommendations for science observing strategies at long baselines. Deep ALMA images of the quasar 3C138 at 97 and 241 GHz are also compared to VLA 43 GHz results, demonstrating an agreement at a level of a few percent. As a result of the extensive program of LBC testing, the highly successful SV imaging at long baselines achieved angular resolutions as fine as 19 mas at ~350 GHz. Observing with ALMA on baselines of up to 15 km is now possible, and opens up new parameter space for submm astronomy.

Molecular complexity builds up at each step of the Sun-like star formation process, starting from simple molecules and ending up in large polyatomic species. Complex organic molecules (COMs; such as methyl formate, HCOOCH$_3$, dymethyl ether, CH$_3$OCH$_3$, formamide, NH$_2$CHO, or glycoaldehyde, HCOCH$_2$OH) are formed in all the components of the star formation recipe (e.g. pre-stellar cores, hot-corinos, circumstellar disks, shocks induced by fast jets), due to ice grain mantle sublimation or sputtering as well as gas-phase reactions. Understanding in great detail the involved processes is likely the only way to predict the ultimate molecular complexity reached in the ISM, as the detection of large molecules is increasingly more difficult with the increase of the number of atoms constituting them. Thanks to the recent spectacular progress of astronomical observations, due to the Herschel (sub-mm and IR), IRAM and SMA (mm and sub-mm), and NRAO (cm) telescopes, an enormous activity is being developed in the field of Astrochemistry, extending from astronomical observatories to chemical laboratories. We are involved in several observational projects providing unbiased spectral surveys (in the 80-300 and 500-2000 GHz ranges) with unprecedented sensitivity of templates of dense cores and protostars. Forests of COM lines have been detected. In this chapter we will focus on the chemistry of both cold prestellar cores and hot shocked regions, (i) reviewing results and open questions provided by mm-FIR observations, and (ii) showing the need of carrying on the observations of COMs at lower frequencies, where SKA will operate. We will also emphasize the importance of analysing the spectra by the light of the experimental studies performed by our team, who is investigating the chemical effects induced by ionising radiation bombarding astrophysically relevant ices.

The formation of complex organic molecules (COMs) in protostellar environments is a hotly debated topic. In particular, the relative importance of the gas phase processes as compared to a direct formation of COMs on the dust grain surfaces is so far unknown. We report here the first high-resolution images of acetaldehyde (CH$_3$CHO) emission towards the chemically rich protostellar shock L1157-B1, obtained at 2 mm with the IRAM Plateau de Bure interferometer. Six blueshifted CH$_3$CHO lines with $E_{\rm u}$ = 26-35 K have been detected. The acetaldehyde spatial distribution follows the young ($\sim$ 2000 yr) outflow cavity produced by the impact of the jet with the ambient medium, indicating that this COM is closely associated with the region enriched by iced species evaporated from dust mantles and released into the gas phase. A high CH$_3$CHO relative abundance, 2-3 $\times$ 10$^{-8}$, is inferred, similarly to what found in hot-corinos. Astrochemical modelling indicates that gas phase reactions can produce the observed quantity of acetaldehyde only if a large fraction of carbon, of the order of 0.1%, is locked into iced hydrocarbons.

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

The enormous radiative and mechanical luminosities of massive stars impact a vast range of scales and processes, from the reionization of the universe, to the evolution of galaxies, to the regulation of the interstellar medium, to the formation of star clusters, and even to the formation of planets around stars in such clusters. Two main classes of massive star formation theory are under active study, Core Accretion and Competitive Accretion. In Core Accretion, the initial conditions are self-gravitating, centrally concentrated cores that condense with a range of masses from the surrounding, fragmenting clump environment. They then undergo relatively ordered collapse via a central disk to form a single star or a small-N multiple. In this case, the pre-stellar core mass function has a similar form to the stellar initial mass function. In Competitive Accretion, the material that forms a massive star is drawn more chaotically from a wider region of the clump without passing through a phase of being in a massive, coherent core. In this case, massive star formation must proceed hand in hand with star cluster formation. If stellar densities become very high near the cluster center, then collisions between stars may also help to form the most massive stars. We review recent theoretical and observational progress towards understanding massive star formation, considering physical and chemical processes, comparisons with low and intermediate-mass stars, and connections to star cluster formation.

In the low-mass regime, it is found that the gas-phase abundances of C-bearing molecules in cold starless cores rapidly decrease with increasing density, as the molecules form mantles on dust grains. We study CO depletion in 102 massive clumps selected from the ATLASGAL 870 micron survey, and investigate its correlation with evolutionary stage and with the physical parameters of the sources. Moreover, we study the gradients in [12C]/[13C] and [18O]/[17O] isotopic ratios across the inner Galaxy, and the virial stability of the clumps. We use low-J emission lines of CO isotopologues and the dust continuum emission to infer the depletion factor fD. RATRAN one-dimensional models were also used to determine fD and to investigate the presence of depletion above a density threshold. The isotopic ratios and optical depth were derived with a Bayesian approach. We find a significant number of clumps with a large fD, up to ~20. Larger values are found for colder clumps, thus for earlier evolutionary phases. For massive clumps in the earliest stages of evolution we estimate the radius of the region where CO depletion is important to be a few tenths of a pc. Clumps are found with total masses derived from dust continuum emission up to ~20 times higher than the virial mass, especially among the less evolved sources. These large values may in part be explained by the presence of depletion: if the CO emission comes mainly from the low-density outer layers, the molecules may be subthermally excited, leading to an overestimate of the dust masses. CO depletion in high-mass clumps seems to behave as in the low-mass regime, with less evolved clumps showing larger values for the depletion than their more evolved counterparts, and increasing for denser sources. The C and O isotopic ratios are consistent with previous determinations, and show a large intrinsic scatter.

We present the results of a line identification analysis using data from the IRAM Plateau de Bure Inferferometer, focusing on six massive star-forming hot cores: G31.41+0.31, G29.96-0.02, G19.61-0.23, G10.62-0.38, G24.78+0.08A1 and G24.78+0.08A2. We identify several transitions of vibrationally excited methyl formate (HCOOCH$_3$) for the first time in these objects as well as transitions of other complex molecules, including ethyl cyanide (C$_2$H$_5$CN), and isocyanic acid (HNCO). We also postulate a detection of one transition of glycolaldehyde (CH$_2$(OH)CHO) in two new hot cores. We find G29.96-0.02, G19.61-0.23, G24.78+0.08A1 and 24.78+0.08A2 to be chemically very similar. G31.41+0.31, however, is chemically different: it manifests a larger chemical inventory and has significantly larger column densities. We suggest that it may represent a different evolutionary stage to the other hot cores in the sample, or it may surround a star with a higher mass. We derive column densities for methyl formate in G31.41+0.31, using the rotation diagram method, of $\times$10$^{17}$ cm$^{-2}$ and a T$_{rot}$ of $\sim$170 K. For G29.96-0.02, G24.78+0.08A1 and G24.78+0.08A2, glycolaldehyde, methyl formate and methyl cyanide all seem to trace the same material and peak at roughly the same position towards the dust emission peak. For G31.41+0.31, however, glycolaldehyde shows a different distribution to methyl formate and methyl cyanide and seems to trace the densest, most compact inner part of hot cores.

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.

Chemical models predict that the deuterated fraction (the column density ratio between a molecule containing D and its counterpart containing H) of N2H+, Dfrac(N2H+), is high in massive pre-protostellar cores and rapidly drops of an order of magnitude after the protostar birth, while that of HNC, Dfrac(HNC), remains constant for much longer. We tested these predictions by deriving Dfrac(HNC) in 22 high-mass star forming cores divided in three different evolutionary stages, from high-mass starless core candidates (HMSCs, 8) to high-mass protostellar objects (HMPOs, 7) to Ultracompact HII regions (UCHIIs, 7). For all of them, Dfrac (N2H+) was already determined through IRAM-30m Telescope observations, which confirmed the theoretical rapid decrease of Dfrac(N2H+) after protostar birth (Fontani et al. 2011). Therefore our comparative study is not affected by biases introduced by the source selection. We have found average Dfrac(HNC) of 0.012, 0.009 and 0.008 in HMSCs, HMPOs and UCHIIs, respectively, with no statistically significant differences among the three evolutionary groups. These findings confirm the predictions of the chemical models, and indicate that large values of Dfrac(N2H+) are more suitable than large values of Dfrac(HNC) to identify cores on the verge of forming high-mass stars, likewise what found in the low-mass regime.

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

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.

How do stars that are more massive than the Sun form, and thus how is the stellar initial mass function (IMF) established? Such intermediate- and high-mass stars may be born from relatively massive pre-stellar gas cores, which are more massive than the thermal Jeans mass. The Turbulent Core Accretion model invokes such cores as being in approximate virial equilibrium and in approximate pressure equilibrium with their surrounding clump medium. Their internal pressure is provided by a combination of turbulence and magnetic fields. Alternatively, the Competitive Accretion model requires strongly sub-virial initial conditions that then lead to extensive fragmentation to the thermal Jeans scale, with intermediate- and high-mass stars later forming by competitive Bondi-Hoyle accretion. To test these models, we have identified four prime examples of massive (~100Msun) clumps from mid-infrared extinction mapping of infrared dark clouds (IRDCs). Fontani et al. found high deuteration fractions of N2H+ in these objects, which are consistent with them being starless. Here we present ALMA observations of these four clumps that probe the N2D+(3-2) line at 2.3" resolution. We find six N2D+ cores and determine their dynamical state. Their observed velocity dispersions and sizes are broadly consistent with the predictions of the Turbulent Core model of self-gravitating, magnetized (with Alfven Mach number m_A~1) and virialized cores that are bounded by the high pressures of their surrounding clumps. However, in the most massive cores, with masses up to ~60Msun, our results suggest that moderately enhanced magnetic fields (so that m_A~0.3) may be needed for the structures to be in virial and pressure equilibrium. Magnetically regulated core formation may thus be important in controlling the formation of massive cores, inhibiting their fragmentation, and thus helping to establish the stellar IMF.

We present the first detection of N2H+ towards a low-mass protostellar outflow, namely the L1157-B1 shock, at about 0.1 pc from the protostellar cocoon. The detection was obtained with the IRAM 30-m antenna. We observed emission at 93 GHz due to the J = 1-0 hyperfine lines. The analysis of the emission coupled with the HIFI CHESS multiline CO observations leads to the conclusion that the observed N2H+(1-0) line originates from the dense (> 10^5 cm-3) gas associated with the large (20-25 arcsec) cavities opened by the protostellar wind. We find a N2H+ column density of few 10^12 cm-2 corresponding to an abundance of (2-8) 10^-9. The N2H+ abundance can be matched by a model of quiescent gas evolved for more than 10^4 yr, i.e. for more than the shock kinematical age (about 2000 yr). Modelling of C-shocks confirms that the abundance of N2H+ is not increased by the passage of the shock. In summary, N2H+ is a fossil record of the pre-shock gas, formed when the density of the gas was around 10^4 cm-3, and then further compressed and accelerated by the shock.

The formation of the first virialized structures in overdensities dates back to ~9 Gyr ago, i.e. in the redshift range z ~ 1.4 - 1.6. Some models of structure formation predict that the star formation activity in clusters was high at that epoch, implying large reservoirs of cold molecular gas. Aiming at finding a trace of this expected high molecular gas content in primeval clusters, we searched for the 12CO(2-1) line emission in the most luminous active galactic nucleus (AGN) of the cluster around the radio galaxy 7C 1756+6520 at z ~ 1.4, one of the farthest spectroscopic confirmed clusters. This AGN, called AGN.1317, is located in the neighbourhood of the central radio galaxy at a projected distance of ~780 kpc. The IRAM Plateau de Bure Interferometer was used to investigate the molecular gas quantity in AGN.1317, observing the 12CO(2-1) emission line. We detect CO emission in an AGN belonging to a galaxy cluster at z ~ 1.4. We measured a molecular gas mass of 1.1 x 10^10 Msun, comparable to that found in submillimeter galaxies. In optical images, AGN.1317 does not seem to be part of a galaxy interaction or merger.We also derived the nearly instantaneous star formation rate (SFR) from Halpha flux obtaining a SFR ~65 Msun/yr. This suggests that AGN.1317 is actively forming stars and will exhaust its reservoir of cold gas in ~0.2-1.0 Gyr.

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

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

OMC-2 FIR 4 is one of the closest known young intermediate-mass protoclusters, located at a distance of 420 pc in Orion. This region is one of the few where the complete 500-2000 GHz spectrum has been observed with the heterodyne spectrometer HIFI on board the Herschel satellite, and unbiased spectral surveys at 0.8, 1, 2 and 3 mm have been obtained with the JCMT and IRAM 30-m telescopes. In order to investigate the morphology of this region, we used the IRAM Plateau de Bure Interferometer to image OMC-2 FIR 4 in the 2-mm continuum emission, as well as in DCO+(2-1), DCN(2-1), C34S(3-2), and several CH3OH lines. In addition, we analysed observations of the NH3(1,1) and (2,2) inversion transitions made with the Very Large Array of the NRAO. The resulting maps have an angular resolution which allows us to resolve structures of 5", equivalent to 2000 AU. Our observations reveal three spatially resolved sources within OMC-2 FIR 4, of one or several solar masses each, with hints of further unresolved substructure within them. Two of these sources have elongated shapes and are associated with dust continuum emission peaks, thus likely containing at least one molecular core each. One of them also displays radio continuum emission, which may be attributed to a young B3-B4 star that dominates the overall luminosity output of the region. The third source identified displays a DCO+(2-1) emission peak, and weak dust continuum emission. Its higher abundance of DCO+ relative to the other two regions suggests a lower temperature and therefore its possible association with either a younger low-mass protostar or a starless core. It may alternatively be part of the colder envelope of OMC-2 FIR 4. Our interferometric observations evidence the complexity of this region, where multiple cores, chemical differentiation and an ionised region all coexist within an area of only 10000 AU.

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

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

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 have mapped in the 2.7 mm continuum and 12CO with the PdBI the IR-dark "tail" that crosses the IC 1396N globule from south to north, and is the most extincted part of this cloud. These observations have allowed us to distinguish all possible associations of molecular hydrogen emission features by revealing the presence of two well-collimated low-mass protostellar outflows at the northern part of the globule. The outflows are located almost in the plane of the sky and are colliding with each other towards the position of a strong 2.12 microns H2 line emission feature.

Infrared-dark high-mass clumps are among the most promising objects to study the initial conditions of the formation process of high-mass stars and rich stellar clusters. In this work, we have observed the (3-2) rotational transition of C18O with the APEX telescope, and the (1,1) and (2,2) inversion transitions of NH3 with the Australia Telescope Compact Array in 21 infrared-dark clouds already mapped in the 1.2 mm continuum, with the aim of measuring basic chemical and physical parameters such as the CO depletion factor (fD), the gas kinetic temperature and the gas mass. In particular, the C18O (3-2) line allows us to derive fD in gas at densities higher than that traced by the (1-0) and (2-1) lines, typically used in previous works. We have detected NH3 and C18O in all targets. The clumps possess mass, H2 column and surface densities consistent with being potentially the birthplace of high-mass stars. We have measured fD in between 5 and 78, with a mean value of 32 and a median of 29. These values are, to our knowledge, larger than the typical CO depletion factors measured towards infrared-dark clouds and high-mass dense cores, and are comparable to or larger than the values measured in low-mass pre-stellar cores close to the onset of the gravitational collapse. This result suggests that the earliest phases of the high-mass star and stellar cluster formation process are characterised by fD larger than in low-mass pre-stellar cores. Thirteen out of 21 clumps are undetected in the 24 {\mu}m Spitzer images, and have slightly lower kinetic temperatures, masses and H2 column densities with respect to the eight Spitzer-bright sources. This could indicate that the Spitzer-dark clumps are either less evolved or are going to form less massive objects.

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

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

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

IGRJ17361-4441 is a newly discovered INTEGRAL hard X-ray transient, located in the globular cluster NGC6388. We report here the results of the X-ray and radio observations performed with Swift, INTEGRAL, RXTE, and the Australia Telescope Compact Array (ATCA) after the discovery of the source on 2011 August 11. In the X-ray domain, IGRJ17361-4441 showed virtually constant flux and spectral parameters up to 18 days from the onset of the outburst. The broad-band (0.5-100 keV) spectrum of the source could be reasonably well described by using an absorbed power-law component with a high energy cut-off (N_H\simeq0.8x10^(22) cm^(-2), {\Gamma}\simeq0.7-1.0, and E_cut\simeq25 keV) and displayed some evidence of a soft component below \sim2 keV. No coherent timing features were found in the RXTE data. The ATCA observation did not detect significant radio emission from IGRJ17361-4441, and provided the most stringent upper limit (rms 14.1 {\mu}Jy at 5.5 GHz) to date on the presence of any radio source close to the NGC6388 center of gravity. The improved position of IGRJ17361-4441 in outburst determined from a recent target of opportunity observation with Chandra, together with the X-ray flux and radio upper limits measured in the direction of the source, argue against its association with the putative intermediate-mass black hole residing in the globular cluster and with the general hypothesis that the INTEGRAL source is a black hole candidate. IGRJ17361-4441 might be more likely a new X-ray binary hosting an accreting neutron star. The ATCA radio non-detection also permits us to derive an upper limit to the mass of the suspected intermediate massive black hole in NGC6388 of <600 M\odot. This is a factor of 2.5 lower than the limit reported previously.

We analyze the long-term evolution of the fluxes of six active galactic nuclei (AGN) - 0923+392, 3C 111, 3C 273, 3C 345, 3C 454.3, and 3C 84 - in the frequency range 80-267 GHz using archival calibration data of the IRAM Plateau de Bure Interferometer. Our dataset spans a long timeline of ~14 years with 974-3027 flux measurements per source. We find strong (factors ~2-8) flux variability on timescales of years for all sources. The flux density distributions of five out of six sources show clear signatures of bi- or even multimodality. Our sources show mostly steep (alpha~0.5-1), variable spectral indices that indicate outflow dominated emission; the variability is most probably due to optical depth variations. The power spectra globally correspond to red-noise spectra with five sources being located between the cases of white and flicker noise and one source (3C 111) being closer to the case of random walk noise. For three sources the low-frequency ends of their power spectra appear to be upscaled in spectral power by factors ~2-3 with respect to the overall powerlaws. In two sources, 3C 454.3 and 3C 84, the 1.3-mm emission preceeds the 3-mm emission by ~55 and ~300 days, respectively, probably due to (combinations of) optical depth and emission region geometry effects. We conclude that the source emission cannot be described by uniform stochastic emission processes; instead, a distinction of "quiescent" and (maybe multiple) "flare" states of the source emission appears to be necessary.

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.

Among the tracers of the earliest phases in the massive star formation process, methanol masers have gained increasing importance. The phenomenological distinction between Class I and II methanol masers is based on their spatial association with objects such as jets, cores, and ultracompact HII regions, but is also believed to correspond to different pumping mechanisms: radiation for Class II masers, collisions for Class I masers. In this work, we have surveyed a large sample of massive star forming regions - 296 objects divided into two groups named 'High' and 'Low' according to their [25-12] and [60-12] IRAS colours - in Class I and II methanol masers. Previous studies indicate that the High sources are likely more evolved. Therefore, the sample can be used to assess the existence of a sequence for the occurrence of Class I and II methanol masers during the evolution of a massive star forming region. We observed the 6 GHz (Class II) CH3OH maser with the Effelsberg 100-m telescope, and the 44 GHz and 95 GHz (Class I) CH3OH masers with the Nobeyama 45-m telescope. We have detected: 55 sources in the Class II line (12 new detections); 27 sources in the 44 GHz Class I line (17 new detections); 11 sources in the 95 GHz Class I line (all except one are new detections). Our statistical analysis shows that the ratio between the detection rates of Class II and Class I methanol masers is basically the same in High and Low sources. Therefore, both masers are equally associated with each evolutionary phase. In contrast, all maser species have about 3 times higher detection rates in High than in Low sources. This might indicate that the phenomena that originate all masers become progressively more active with time, during the earliest evolutionary phases of a high-mass star forming region.

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.

To better understand the initial conditions of the high-mass star formation process, it is crucial to study at high-angular resolution the morphology, the kinematics, and eventually the interactions of the coldest condensations associated with intermediate-/high-mass star forming regions. The paper studies the cold condensations in the intermediate-/high-mass proto-cluster IRAS 05345+3157, focusing the attention on the interaction with the other objects in the cluster. We have performed millimeter high-angular resolution observations, both in the continuum and several molecular lines, with the PdBI and the SMA. In a recent paper, we have already published part of these data. The main finding of that work was the detection of two cold and dense gaseous condensations, called N and S (masses ~2 and ~9 M_sun), characterised by high values of the deuterium fractionation (~0.1 in both cores). In this paper, we present a full report of the observations, and a complete analysis of the data obtained. The millimeter maps reveal the presence of 3 cores inside the interferometers primary beam, called C1-a, C1-b and C2. None of them are associated with cores N and S. C1-b is very likely associated with a newly formed early-B ZAMS star embedded inside a hot-core, while C1-a is more likely associated with a class 0 intermediate-mass protostar. The nature of C2 is unclear. Both C1-a and C1-b are good candidates as driving sources of a powerful CO outflow, which strongly interacts with N and S, as demonstrated by the velocity gradient across both condensations. Our major conclusion is that the chemical properties of these pre-stellar cores are similar to those observed in low-mass isolated ones, while the kinematics is dominated by the turbulence triggered by the CO outflow and can influece their evolution.

We have observed the deuterated gas in the high-mass star formation region IRAS 05345+3157 at high-angular resolution, in order to determine the morphology and the nature of such gas. We have mapped the N2H+ (1-0) line with the Plateau de Bure Interferometer, and the N2D+ (3-2) and N2H+ (3-2) lines with the Submillimeter Array. The N2D+ (3-2) integrated emission is concentrated in two condensations, with masses of 2-3 and 9 M_sun and diameters of 0.05 and 0.09 pc, respectively. The high deuterium fractionation (0.1) and the line parameters in the N2D+ condensations indicate that they are likely low- to intermediate-mass pre-stellar cores, even though other scenarios are possible.

We have observed several emission lines of two Nitrogen-bearing (C2H5CN and C2H3CN) and two Oxygen-bearing (CH3OCH3 and HCOOCH3) molecules towards a sample of well-known hot molecular cores (HMCs) in order to check whether the chemical differentiation seen in the Orion-HMC and W3(H_2O) between O- and N-bearing molecules is a general property of HMCs. With the IRAM-30m telescope we have observed 12 HMCs in 21 bands, centered at frequencies from 86250 to 258280 MHz. The rotational temperatures obtained range from ~100 to ~150 K in these HMCs. Single Gaussian fits performed to unblended lines show a marginal difference in the line peak velocities of the C2H5CN and CH3OCH3 lines, indicating a possible spatial separation between the region traced by the two molecules. On the other hand, neither the linewidths nor the rotational temperatures and column densities confirm such a result. By comparing the abundance ratio of the pair C2H5CN/C2H3CN with the predictions of theoretical models, we derive that the age of our cores ranges between 3.7 and 5.9x10^{4} yrs. The abundances of C2H5CN and C2H3CN are strongly correlated, as expected from theory which predicts that C2H3CN is formed through gas phase reactions involving C2H5CN. A correlation is also found between the abundances of C2H3CN and CH3OCH3, and C2H5CN and CH3OCH3. In all tracers the fractional abundances increase with the H_2 column density while they are not correlated with the gas temperature.

We have measured the deuterium fractionation (through the column density ratio N(N2D+)/N(N2H+)) and the CO depletion factor (ratio between expected and observed CO abundance) in a sample of 10 high-mass protostellar candidates, in order to understand whether the earliest evolutionary stages of high-mass stars have chemical characteristics similar to those of low-mass ones. The observations were carried out with the IRAM-30m telescope and the JCMT. We have detected N2D+ emission in 7 of the 10 sources of our sample, and found an average value N(N2D+)/N(N2H+)~0.015. This value is 3 orders of magnitude larger than the interstellar D/H ratio, indicating the presence of cold and dense gas, in which the physical-chemical conditions are similar to those observed in low-mass pre-stellar cores. Also, the integrated CO depletion factors show that in the majority of the sources the expected CO abundances are larger than the observed values, with a median ratio of 3.2. In principle, the cold gas that gives origin to the N2D+ emission can be the remnant of the massive molecular core in which the high-mass (proto-)star was born, not yet heated up by the central object. If so, our results indicate that the chemical properties of the clouds in which high-mass stars are born are similar to their low-mass counterparts. Alternatively, this cold gas can be located into one (or more) starless core (cores) near the protostellar object. Due to the poor angular resolution of our data, we cannot decide which is the correct scenario.

In an ongoing effort to identify and study high-mass protostellar candidates we have observed in various tracers a sample of 235 sources selected from the IRAS Point Source Catalog, mostly with dec < -30 deg, with the SEST antenna at millimeter wavelengths. The sample contains 142 Low sources and 93 High, which are believed to be in different evolutionary stages. Both sub-samples have been studied in detail by comparing their physical properties and morphologies. Massive dust clumps have been detected in all but 8 regions, with usually more than one clump per region. The dust emission shows a variety of complex morphologies, sometimes with multiple clumps forming filaments or clusters. The mean clump has a linear size of ~0.5 pc, a mass of ~320 Msolar for a dust temperature Td=30 K, an H_2 density of 9.5E5 cm-3, and a surface density of 0.4 g cm-2. The median values are 0.4 pc, 102 Msolar, 4E4 cm-3, and 0.14 g cm-2, respectively. The mean value of the luminosity-to-mass ratio, L/M ~99 Lsolar/Msolar, suggests that the sources are in a young, pre-ultracompact HII phase. We have compared the millimeter continuum maps with images of the mid-IR MSX emission, and have discovered 95 massive millimeter clumps non-MSX emitters, either diffuse or point-like, that are potential prestellar or precluster cores. The physical properties of these clumps are similar to those of the others, apart from the mass that is ~3 times lower than for clumps with MSX counterpart. Such a difference could be due to the potential prestellar clumps having a lower dust temperature. The mass spectrum of the clumps with masses above M ~100 Msolar is best fitted with a power-law dN/dM proportional to M-alpha with alpha=2.1, consistent with the Salpeter (1955) stellar IMF, with alpha=2.35.

(Abridged) We have observed CS and C17O lines, and 1.2 mm cont. emission towards a sample of 130 high-mass protostellar candidates with DEC<-30 deg. This is the first step of the southern extension of a project started more than a decade ago aimed at the identification of massive protostellar candidates. We selected from the IRAS PSC 429 sources which potentially are compact molecular clouds. The sample is divided into two groups: the 298 sources with [25-12]>0.57 and [60-12]>1.30 we call 'High' sources, the remaining 131 we call 'Low' sources. In this paper, we check the association with dense gas and dust in 130 'Low' sources. We find a detection rate of ca. 85% in CS, demonstrating a tight association with dense molecular clumps. Among the sources detected in CS, ca. 76% have also been detected in C17O and ca. 93% in the 1.2 mm cont. Mm-cont. maps show the presence of clumps with diameters 0.2-2 pc and masses from a few Msun to 10^5 Msun; H2 volume densities lie between ca. 10^{4.5} and 10^{5.5} cm^{-3}. The L(bol) are 10^3-10^6 Lsun, consistent with embedded high-mass objects. Based on our results and those found in the literature for other samples, we conclude that our sources are massive objects probably in a stage prior to the formation of an HII region. We propose a scenario in which 'High' and 'Low' sources are both made of a massive clump hosting a high-mass protostellar candidate and a nearby stellar cluster. The difference might be due to the fact that the IRAS 12mu flux, the best discriminant between the two groups, is dominated by the emission from the cluster in 'Lows' and from the massive protostellar object in 'Highs'.

We present observations of continuum and molecular lines towards the protostar candidates IRAS 21307+5049 and IRAS 22172+5549. Single-dish maps in the 12CO(2--1), C18O(2--1), HCO+(1--0) lines and sub-mm continuum are compared with interferometric maps in the 12CO(1--0) line and 3 mm continuum, and with mid- and near-infrared images. A plausible interpretation of our data based on the continuum maps and spectral energy distributions is that two components are present: a compact molecular core, responsible for the continuum emission at wavelengths longer than ~25 \mum, and a cluster of stars located close to the center of the core, but not spatially coincident with it, responsible for the emission at shorter wavelengths. The core is approximately located at the center of the associated molecular outflow, detected for both sources in the 12CO(1--0) and (2--1) lines. The cores have masses of ~50 M_{\odot}, and luminosities of ~10^{3} L_{\odot}. The outflows parameters are consistent with those typically found in high-mass young stellar objects. Our results support the hypothesis that in these sources the luminosity is dominated by accretion rather than by nuclear burning. We conclude that the sources embedded inside the cores are likely protostars with mass ~5-8 M_{\odot}.

We present the results of a multi-line and continuum study towards the source IRAS 23385+6053,performed with the IRAM-30m telescope, the Plateau de Bure Interferometer, the Very Large Array Interferometer and the James Clerk Maxwell Telescope. The new results confirm our earlier findings, namely that IRAS 23385+6053 is a good candidate high-mass protostellar object, precursor of an ultracompact H$_{II}$ region. The source is roughly composed of two regions: a molecular core $\sim0.03\div0.04$ pc in size, with a temperature of $\sim40$ K and an H$_{2}$ volume density of the order of 10$^{7}$ cm$^{-3}$, and an extended halo of diameter $\leq$0.4 pc, with an average kinetic temperature of $\sim 15$ K and H$_{2}$ volume density of the order of 10$^{5}$ cm$^{-3}$. The core temperature is much smaller than what is typically found in molecular cores of the same diameter surrounding massive ZAMS stars. We deduce that the core luminosity is between 150 and $1.6\times10^{4}L_{\odot}$, and we believe that the upper limit is near the ``true'' source luminosity. Moreover, by comparing the H$_{2}$ volume density obtained at different radii from the IRAS source, we find that the halo has a density profile of the type $n_{\rm H_{2}}\propto r^{-2.3}$. This suggests that the source is gravitationally unstable. Finally, we demonstrate that the temperature at the core surface is consistent with a core luminosity of $10^3 L_{\odot}$ and conclude that we might be observing a protostar still accreting material from its parental cloud, whose mass at present is $\sim 6 M_{\odot}$.

We have used the IRAM 30-m and FCRAO 14-m telescopes to observe the molecular clumps associated with 12 ultracompact (UC) HII regions in the J=6-5, 8-7 and 13-12 rotational transitions of methyl-acetylene (CH3C2H). Under the assumption of LTE and optically thin emission, we have derived temperature estimates ranging from 30 to 56 K. We estimate that the clumps have diameters of 0.2-1.6 pc, H_2 densities of 10^5-10^6 {cm^{-3}}, and masses of 10^2-2 10^4 M_\odot. We compare these values with those obtained by other authors from different molecular tracers and find that the H_2 density and the temperature inside the clumps vary respectively like n_{H_2} ~ R^{-2.6} and T ~ R^{-0.5}, with R distance from the centre. We also find that the virial masses of the clumps are ~3 times less than those derived from the CH3C2H column densities: we show that a plausible explanation is that magnetic fields play an important role to stabilise the clumps, which are on the verge of gravitational collapse. Finally, we show that the CH3C2H line width increases for decreasing distance from the clump centre: this effect is consistent with infall in the inner regions of the clumps. We conclude that the clumps around UC HII regions are likely to be transient (~10^(5) yr) entities, remnants of isothermal spheres currently undergoing gravitational collapse: the high mass accretion rates (~10^{-2} M_\odot yr^{-1}) lead to massive star formation at the centre of such clumps.