STATCONT is a python-based tool designed to determine the continuum emission level in spectral data, in particular for sources with a line-rich spectrum. The tool inspects the intensity distribution ...of a given spectrum and automatically determines the continuum level by using different statistical approaches. The different methods included in STATCONT are tested against synthetic data. We conclude that the sigma-clipping algorithm provides the most accurate continuum level determination, together with information on the uncertainty in its determination. This uncertainty can be used to correct the final continuum emission level, resulting in the here called ‘corrected sigma-clipping method’ or c-SCM. The c-SCM has been tested against more than 750 different synthetic spectra reproducing typical conditions found towards astronomical sources. The continuum level is determined with a discrepancy of less than 1% in 50% of the cases, and less than 5% in 90% of the cases, provided at least 10% of the channels are line free. The main products of STATCONT are the continuum emission level, together with a conservative value of its uncertainty, and datacubes containing only spectral line emission, i.e., continuum-subtracted datacubes. STATCONT also includes the option to estimate the spectral index, when different files covering different frequency ranges are provided.
Context. The discovery of amino acids in meteorites fallen to Earth and the detection of glycine, the simplest of them, in samples returned from a comet to Earth strongly suggest that the chemistry ...of the interstellar medium is capable of producing such complex organic molecules and that they may be widespread in our Galaxy. Aims. Our goal is to investigate the degree of chemical complexity that can be reached in the interstellar medium, in particular in dense star-forming regions. Methods. We performed an unbiased, spectral line survey toward Sgr B2(N) and (M), two regions where high-mass stars are formed, with the IRAM 30 m telescope in the 3 mm atmospheric transmission window. Partial surveys at 2 and 1.3 mm were performed in parallel. The spectra were analyzed with a simple radiative transfer model that assumes local thermodynamic equilibrium but takes optical depth effects into account. Results. About 3675 and 945 spectral lines with a peak signal-to-noise ratio higher than 4 are detected at 3 mm toward Sgr B2(N) and (M), i.e. about 102 and 26 lines per GHz, respectively. This represents an increase by about a factor of two over previous surveys of Sgr B2. About 70% and 47% of the lines detected toward Sgr B2(N) and (M) are identified and assigned to 56 and 46 distinct molecules as well as to 66 and 54 less abundant isotopologues of these molecules, respectively. In addition, we report the detection of transitions from 59 and 24 catalog entries corresponding to vibrationally or torsionally excited states of some of these molecules, respectively, up to a vibration energy of 1400 cm-1 (2000 K). Excitation temperatures and column densities were derived for each species but should be used with caution. The rotation temperatures of the detected complex molecules typically range from ~50 to 200 K. Among the detected molecules, aminoacetonitrile, n-propyl cyanide, and ethyl formate were reported for the first time in space based on this survey, as were five rare isotopologues of vinyl cyanide, cyanoacetylene, and hydrogen cyanide. We also report the detection of transitions from within twelve new vibrationally or torsionally excited states of known molecules. Absorption features produced by diffuse clouds along the line of sight are detected in transitions with low rotation quantum numbers of many simple molecules and are modeled with ~30–40 velocity components with typical linewidths of ~3–5 km s-1. Conclusions. Although the large number of unidentified lines may still allow future identification of new molecules, we expect most of these lines to belong to vibrationally or torsionally excited states or to rare isotopologues of known molecules for which spectroscopic predictions are currently missing. Significant progress in extending the inventory of complex organic molecules in Sgr B2(N) and deriving tighter constraints on their location, origin, and abundance is expected in the near future thanks to an ongoing spectral line survey at 3 mm with ALMA in its cycles 0 and 1. The present single-dish survey will serve as a solid basis for the line identification and analysis of such an interferometric survey.
Aims. We aim to characterise certain physical properties of high-mass star-forming sites in the NGC 6334 molecular cloud, such as the core mass function (CMF), spatial distribution of cores, and mass ...segregation. Methods. We used the Atacama Large Millimeter/sub-millimeter Array (ALMA) to image the embedded clusters NGC 6334-I and NGC 6334-I(N) in the continuum emission at 87.6 GHz. We achieved a spatial resolution of 1300 au, enough to resolve different compact cores and fragments, and to study the properties of the clusters. Results. We detected 142 compact sources distributed over the whole surveyed area. The ALMA compact sources are clustered in different regions. We used different machine-learning algorithms to identify four main clusters: NGC 6334-I, NGC 6334-I(N), NGC 6334-I(NW), and NGC 6334-E. The typical separations between cluster members range from 4000 au to 12 000 au. These separations, together with the core masses (0.1–100 M⊙), are in agreement with the fragmentation being controlled by turbulence at scales of 0.1 pc. We find that the CMFs show an apparent excess of high-mass cores compared to the stellar initial mass function. We evaluated the effects of temperature and unresolved multiplicity on the derived slope of the CMF. Based on this, we conclude that the excess of high-mass cores might be spurious and due to inaccurate temperature determinations and/or resolution limitations. We searched for evidence of mass segregation in the clusters and we find that clusters NGC 6334-I and NGC 6334-I(N) show hints of segregation with the most massive cores located in the centre of the clusters. Conclusions. We searched for correlations between the physical properties of the four embedded clusters and their evolutionary stage (based on the presence of H II regions and infrared sources). NGC 6334-E appears as the most evolved cluster, already harbouring a well-developed H II region. NGC 6334-I is the second-most evolved cluster with an ultra-compact H II region. NGC 6334-I(N) contains the largest population of dust cores distributed in two filamentary structures and no dominant H II region. Finally, NGC 6334-I(NW) is a cluster of mainly low-mass dust cores with no clear signs of massive cores or H II regions. We find a larger separation between cluster members in the more evolved clusters favouring the role of gas expulsion and stellar ejection with evolution. The mass segregation, seen in the NGC 6334-I and NGC 6334-I(N) clusters, suggests a primordial origin for NGC 6334-I(N). In contrast, the segregation in NGC 6334-I might be due to dynamical effects. Finally, the lack of massive cores in the most evolved cluster suggests that the gas reservoir is already exhausted, while the less evolved clusters still have a large gas reservoir along with the presence of massive cores. In general, the fragmentation process of NGC 6334 at large scales (from filament to clump, i.e. at about 1 pc) is likely governed by turbulent pressure, while at smaller scales (scale of cores and sub-fragments, i.e. a few hundred au) thermal pressure starts to be more significant.
Aims. We describe the assignment of a previously unidentified interstellar absorption line to ArH+ and discuss its relevance in the context of hydride absorption in diffuse gas with a low H2 ...fraction. The confidence of the assignment to ArH+ is discussed, and the column densities are determined toward several lines of sight. The results are then discussed in the framework of chemical models, with the aim of explaining the observed column densities. Methods. We fitted the spectral lines with multiple velocity components, and determined column densities from the line-to-continuum ratio. The column densities of ArH+ were compared to those of other species, tracing interstellar medium (ISM) components with different H2 abundances. We constructed chemical models that take UV radiation and cosmic ray ionization into account. Results. Thanks to the detection of two isotopologues, 36ArH+ and 38ArH+, we are confident about the carrier assignment to ArH+. NeH+ is not detected with a limit of NeH+/ArH+ ≤ 0.1. The derived column densities agree well with the predictions of chemical models. ArH+ is a unique tracer of gas with a fractional H2 abundance of 10-4 − 10-3 and shows little correlation to H2O+, which traces gas with a fractional H2 abundance of ≈0.1. Conclusions. A careful analysis of variations in the ArH+, OH+, H2O+, and HF column densities promises to be a faithful tracer of the distribution of the H2 fractional abundance by providing unique information on a poorly known phase in the cycle of interstellar matter and on its transition from atomic diffuse gas to dense molecular gas traced by CO emission. Abundances of these species put strong observational constraints upon magnetohydrodynamical (MHD)simulations of the interstellar medium, and potentially could evolve into a tool characterizing the ISM. Paradoxically, the ArH+ molecule is a better tracer of almost purely atomic hydrogen gas than Hi itself, since Hi can also be present in gas with a significant molecular content, but ArH+ singles out gas that is >99.9% atomic.
Aims. Star formation involves the collapse of gas from the scale of giant molecular clouds down to dense cores. Our aim is to trace the velocities in the filamentary, massive star-forming region NGC ...6334 and to explain its dynamics. Methods. The main filament was mapped with the single-dish telescope APEX in HCO+ (J = 3–2) at 267.6 GHz to trace the dense gas. In order to reproduce the position−velocity diagram, we use a 3D radiative transfer code and create a model of a cylinder that undergoes a gravitational collapse toward its center. Results. We find a velocity gradient in the filament from the end toward its center, with the highest masses being found at both ends. Similar velocities have been predicted by recent calculations of the gravitational collapse of a sheet or cylinder of gas, and the observed velocities are consistent with these predictions. The 3D structure is revealed by taking the inclination and curvature of the filament into account. The free-fall collapse timescale of the filamentary molecular cloud is estimated to be ~1 Myr.
As part of our effort to search for circumstellar disks around high-mass stellar objects, we observed the well-known core G31.41 +0.31 with ALMA at 1.4 mm with an angular resolution of ~0.′′22 (~1700 ...au). The dust continuum emission has been resolved into two cores namely Main and NE. The Main core, which has the stronger emission and is the more chemically rich, has a diameter of ~5300 au, and is associated with two free-free continuum sources. The Main core looks featureless and homogeneous in dust continuum emission and does not present any hint of fragmentation. Each transition of CH3CN and CH3OCHO, both ground and vibrationally excited, as well as those of CH3CN isotopologues, shows a clear velocity gradient along the NE–SW direction, with velocity linearly increasing with distance from the center, consistent with solid-body rotation. However, when comparing the velocity field of transitions with different upper level energies, the rotation velocity increases with increasing energy of the transition, which suggests that the rotation speeds up toward the center. Spectral lines towardtoward the dust continuum peak show an inverse P-Cygni profile that supports the existence of infall in the core. The infall velocity increases with the energy of the transition suggesting that the infall is accelerating toward the center of the core, consistent with gravitational collapse. Despite the monolithic appearance of the Main core, the presence of red-shifted absorption, the existence of two embedded free-free sources at the center, and the rotational spin-up are consistent with an unstable core undergoing fragmentation with infall and differential rotation due to conservation of angular momentum. Therefore, the most likely explanation for the monolithic morphology is that the large opacity of the dust emission prevents the detection of any inhomogeneity in the core.
Aims.Our current knowledge of high-mass star formation is mainly based on follow-up studies of bright sources found by IRAS, and is thus biased against its earliest phases, inconspicuous at infrared ...wavelengths. We therefore started searching, in an unbiased way and in the closest high-mass star-forming complexes, for the high-mass analogs of low-mass pre-stellar cores and class 0 protostars. Methods.We have made an extensive 1.2 mm continuum mosaicing study of the Cygnus X molecular cloud complex using the MAMBO cameras at the IRAM 30 m telescope. The $\sim 3^{\circ^2}$ imaged areas cover all the high-column density (AV ≥ 15 mag) clouds of this nearby (~1.7 kpc) cloud complex actively forming OB stars. We then compared our millimeter maps with mid-infrared images, and have made SiO(2-1) follow-up observations of the best candidate progenitors of high-mass stars. Results.Our complete study of Cygnus X with ~0.09 pc resolution provides, for the first time, an unbiased census of massive young stellar objects. We discover 129 massive dense cores (FWHM size ~0.1 pc, $M_{1.2~mm}$ = 4-950 $\mbox{$M_\odot$}$, volume-averaged density ~105 cm-3), among which ~42 are probable precursors of high-mass stars. A large fraction of the Cygnus X dense cores (2/3 of the sample) remain undetected by the MSX satellite, regardless of the mass range considered. Among the most massive (≥40 $\mbox{$M_\odot$}$) cores, infrared-quiet objects are driving powerful outflows traced by SiO emission. Our study qualifies 17 cores as good candidates for hosting massive infrared-quiet protostars, while up to 25 cores potentially host high-luminosity infrared protostars. We fail to discover the high-mass analogs of pre-stellar dense cores (~0.1 pc, > 104 cm-3) in Cygnus X, but find several massive starless clumps (~ 0.8 pc, 7 $\times$ 103 cm-3) that might be gravitationally bound. Conclusions.Since our sample is derived from a single molecular complex and covers every embedded phase of high-mass star formation, it gives the first statistical estimates of their lifetime. In contrast to what is found for low-mass class 0 and class I phases, the infrared-quiet protostellar phase of high-mass stars may last as long as their better-known high-luminosity infrared phase. The statistical lifetimes of high-mass protostars and pre-stellar cores (~ 3 $\times$ 104 yr and < 103 yr) in Cygnus X are one and two order(s) of magnitude smaller, respectively, than what is found in nearby, low-mass star-forming regions. We therefore propose that high-mass pre-stellar and protostellar cores are in a highly dynamic state, as expected in a molecular cloud where turbulent processes dominate.
We have mapped 20 spectral lines in the Central Molecular Zone (CMZ) around the Galactic Centre, emitting from 85.3 to 93.3 GHz. This work used the 22-m Mopra radio telescope in Australia, equipped ...with the 8-GHz bandwidth University of New South Wales-Mopra Spectrometer (UNSW-MOPS) digital filter bank, obtaining ∼2 km s−1 spectral and ∼40 arcsec spatial resolution. The lines measured include emission from the c-C3H2, CH3CCH, HOCO+, SO, H13CN, H13CO+, SO, H13NC, C2H, HNCO, HCN, HCO+, HNC, HC3N, 13CS and N2H+ molecules. The area covered is Galactic longitude −0°.7 to 1°.8 and latitude −0°.3 to 0°.2, including the bright dust cores around Sgr A, Sgr B2, Sgr C and G1.6−0.025. We present images from this study and conduct a principal component analysis on the integrated emission from the brightest eight lines. This is dominated by the first component, showing that the large-scale distribution of all molecules is very similar. We examine the line ratios and optical depths in selected apertures around the bright dust cores, as well as for the complete mapped region of the CMZ. We highlight the behaviour of the bright HCN, HNC and HCO+ line emission, together with that from the 13C isotopologues of these species, and compare the behaviour with that found in extragalactic sources where the emission is unresolved spatially. We also find that the isotopologue line ratios (e.g. HCO+/H13CO+) rise significantly with increasing redshifted velocity in some locations. Line luminosities are also calculated and compared to that of CO, as well as to line luminosities determined for external galaxies.
Context. Numerical simulations have explored the possibility of forming molecular clouds through either a quasi-static, self-gravitating mechanism or the collision of gas streams or lower density ...clouds. They also quantitatively predict the distribution of matter at the transition from atomic to molecular gases. Aims. We aim to observationally test these models by studying the environment of W43, a molecular cloud complex recently identified near the tip of the Galactic long bar. Methods. Using Galaxy-wide H i and 12CO 1–0 surveys, we searched for gas flowing toward the W43 molecular cloud complex. We also estimated the H i and H2 mass surface densities to constrain the transition from atomic to molecular gas around and within W43. Results. We found three cloud ensembles within the position-velocity diagrams of 12CO and H i gases. They are separated by ~20 km s-1 along the line of sight and extend into the 13CO velocity structure of W43. Since their velocity gradients are consistent with free fall, they could be nearby clouds attracted by and streaming toward the W43 ~107 M⊙ potential well. We show that the H i surface density, ΣH i = 45−85 M⊙pc-2, does not reach any threshold level but increases when entering the 130 pc-wide molecular complex previously defined. This suggests that an equilibrium between H2 formation and photodissociation has not yet been reached. The H2-to-H i ratio measured over the W43 region and its surroundings, RH2 ~ 3.5±23, is high, indicating that most of the gas is already in molecular form in W43 and in structures several hundred parsecs downstream along the Scutum-Centaurus arm. Conclusions. The W43 molecular cloud complex may have formed and, in fact may still be accreting mass from the agglomeration of clouds. Already in the molecular-dominated regime, most of these clouds are streaming from the Scutum-Centaurus arm. This clearly disagrees with quasi-static and steady-state models of molecular cloud formation.