Source author record

T. Csengeri

T. Csengeri appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

64works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

64 published item(s)

preprint2022arXiv

ALMA-IMF III -- Investigating the origin of stellar masses: Top-heavy core mass function in the W43-MM2&MM3 mini-starburst

The ALMA-IMF Large Program observed the W43-MM2-MM3 ridge, whose 1.3mm and 3mm ALMA 12m array continuum images reach a 2500au spatial resolution. We used both the best-sensitivity and the line-free ALMA-IMF images, reduced the noise with the multi-resolution segmentation technique MnGSeg, and derived the most complete and most robust core catalog possible. Using two different extraction software packages, getsf and GExt2D, we identified 200 compact sources, whose 100 common sources have on average fluxes consistent to within 30%. We filtered sources with non-negligible free-free contamination and corrected fluxes from line contamination, resulting in a W43-MM2-MM3 catalog of 205 getsf cores. With a median deconvolved FWHM size of 3400au, core masses range from 0.1Msun to 70Msun and the getsf catalog is 90% complete down to 0.8Msun. The high-mass end of the core mass function (CMF) of W43-MM2-MM3 is top-heavy compared to the canonical IMF. Fitting the cumulative CMF with a single power law of the form N(>logM)\propto M^a, we measured a=-0.95\pm0.04, compared to the canonical a=-1.35 Salpeter IMF slope. The slope of the CMF is robust with respect to map processing, extraction software package, and reasonable variations in the assumptions taken to estimate core masses. We explore several assumptions on how cores transfer their mass to stars and sub-fragment to predict the IMF resulting from the W43-MM2-MM3 CMF. In stark contrast to the commonly accepted paradigm, our result argues against the universality of the CMF shape. More robust functions of the star-formation efficiency and core sub-fragmentation are required to better predict the resulting IMF, here suggested to remain top-heavy at the end of the star-formation phase. If confirmed, the IMFs emerging from starburst events could inherit their top-heavy shape from their parental CMFs, challenging the IMF universality.

preprint2022arXiv

SOFIA/GREAT observations of OD and OH rotational lines towards high-mass star forming regions

Only recently, OD, the deuterated isotopolog of hydroxyl, OH, has become accessible in the interstellar medium; spectral lines from both species have been observed in the supra-Terahertz and far infrared regime. Here we study rotational lines of OD and OH towards 13 Galactic high-mass star forming regions, with the aim to constrain the OD abundance and infer the deuterium fractionation of OH in their molecular envelopes. We used the Stratospheric Observatory for Infrared Astronomy (SOFIA) to observe the $^2Π_{3/2}$ $J=5/2-3/2$ ground-state transition of OD at 1.3 THz ($215~μ$m) and the rotationally excited OH line at 1.84 THz ($163~μ$m). We also used published high-spectral-resolution SOFIA data of the OH ground-state transition at 2.51 THz ($119.3~μ$m). Our results show that absorption from the $^2Π_{3/2}$ OD $J=5/2-3/2$ ground-state transition is prevalent in the dense clumps surrounding active sites of high-mass star formation. We performed detailed radiative transfer modelling to investigate the OD abundance profile in the inner envelope for a large fraction of our sample. Our modelling suggests that part of the absorption arises from the denser inner parts, while the bulk of it as seen with SOFIA originates in the outer, cold layers of the envelope for which our constraints on the molecular abundance suggest a strong enhancement in deuterium fractionation. We find a weak negative correlation between the OD abundance and the bolometric luminosity to mass ratio, an evolutionary indicator, suggesting a slow decrease of OD abundance with time. A comparison with HDO shows a similarly high deuterium fractionation for the two species in the cold envelopes, which is of the order of 0.48% for the best studied source, G34.26+0.15. Our results are consistent with chemical models that favour rapid exchange reactions to form OD in the dense cold gas.

preprint2022arXiv

Sulphur-rich cold gas around the hot core precursor G328.2551-0.5321. An APEX unbiased spectral survey of the 2 mm, 1.2 mm, and 0.8 mm atmospheric windows

During star formation, the dense gas undergoes significant chemical evolution leading to the emergence of a rich variety of molecules associated with hot cores and hot corinos. The physical and chemical conditions are poorly constrained; the early phases of emerging hot cores in particular represent an unexplored territory. We provide here a full molecular inventory of a massive protostellar core that is proposed to be a precursor of a hot core. We performed an unbiased spectral survey towards the hot core precursor associated with clump G328.2551-0.5321 between 159GHz and 374GHz. To identify the spectral lines, we used rotational diagrams and radiative transfer modelling assuming LTE. We detected 39 species and 26 isotopologues, and were able to distinguish a warm and compact inner region, a colder more extended envelope, and the kinematic signatures of the accretion shocks that have previously been observed with ALMA. We associate most of the emission of the small molecules with the cold gas, while the molecular emission of the warm gas is enriched by complex organic molecules (COMs). We find a high abundance of S-bearing molecules in the cold gas phase suggesting a low sulphur depletion, with a factor of > 1%. We identify nine COMs in the warm gas, four in the cold gas, and four towards the accretion shocks. The high abundances of S-bearing species originating from the undisturbed gas may suggest a contribution from shocked gas at the outflow cavity walls. The molecular composition of the warm gas is similar to that of both hot cores and hot corinos, but the molecular abundances are closer to the values found towards hot corinos than to values found towards hot cores. Considering the compactness of the warm region and its moderate temperature, we suggest that thermal desorption has not been completed towards this object yet, representing an early phase of the emergence of hot cores.

preprint2022arXiv

The SEDIGISM survey: Molecular cloud morphology. I. Classification and star formation

We present one of the very first extensive classifications of a large sample of molecular clouds based on their morphology. This is achieved using a recently published catalogue of 10663 clouds obtained from the first data release of the SEDIGISM survey. The clouds are classified into four different morphologies by visual inspection and using an automated algorithm -- J plots. The visual inspection also serves as a test for the J plots algorithm, as this is the first time it has been used on molecular gas. Generally, it has been found that the structure of molecular clouds is highly filamentary and our observations indeed verify that most of our molecular clouds are elongated structures. Based on our visual classification of the 10663 SEDIGISM clouds, 15% are ring-like, 57% are elongated, 15% are concentrated and 10% are clumpy clouds. The remaining clouds do not belong to any of these morphology classes and are termed unclassified. We compare the SEDIGISM molecular clouds with structures identified through other surveys, i.e. ATLASGAL elongated structures and the bubbles from Milky Way Project (MWP). We find that many of the ATLASGAL and MWP structures are velocity coherent. ATLASGAL elongated structures overlap with ~ 21% of the SEDIGISM elongated structures (elongated and clumpy clouds) and MWP bubbles overlap with ~ 25% of the SEDIGISM ring-like clouds. We also analyse the star-formation associated with different cloud morphologies using two different techniques. The first technique examines star formation efficiency (SFE) and the dense gas fraction (DGF), based on SEDIGISM clouds and ATLASGAL clumps data. The second technique uses the high-mass star formation (HMSF) threshold for molecular clouds. The results indicate that clouds with ring-like and clumpy morphologies show a higher degree of star formation.

preprint2022arXiv

The SEDIGISM survey: Molecular cloud morphology. II. Integrated source properties

The Structure, Excitation, and Dynamics of the Inner Galactic InterStellar Medium (SEDIGISM) survey has produced high (spatial and spectral) resolution $^{13}$CO (2-1) maps of the Milky Way. It has allowed us to investigate the molecular interstellar medium in the inner Galaxy at an unprecedented level of detail and characterise it into molecular clouds. In a previous paper, we have classified the SEDIGISM clouds into four morphologies. However, how the properties of the clouds vary for these four morphologies is not well understood. Here, we use the morphological classification of SEDIGISM clouds to find connections between the cloud morphologies, their integrated properties, and their location on scaling relation diagrams. We observe that ring-like clouds show the most peculiar properties, having, on average, higher masses, sizes, aspect ratios and velocity dispersions compared to other morphologies. We speculate that this is related to the physical mechanisms that regulate their formation and evolution, for example, turbulence from stellar feedback can often results in the creation of bubble-like structures. We also see a trend of morphology with virial parameter whereby ring-like, elongated, clumpy and concentrated clouds have virial parameters in a decreasing order. Our findings provide a foundation for a better understanding of the molecular cloud behaviour based on their measurable properties.

preprint2021arXiv

ALMA-IMF I -- Investigating the origin of stellar masses: Introduction to the Large Program and first results

The ALMA-IMF Large Program imaged a total noncontiguous area of 53pc2, covering 15 extreme, nearby protoclusters of the Milky Way. They were selected to span relevant early protocluster evolutionary stages. Our 1.3mm and 3mm observations provide continuum images that are homogeneously sensitive to point-like cores with masses of 0.2 and 0.6Msun, respectively, with a matched spatial resolution of 2000au. We also detect lines that probe the protocluster structure, kinematics, chemistry, and feedback over scales from clouds to filaments to cores. We classify ALMA-IMF protoclusters as Young, Intermediate, or Evolved based on the amount of dense gas in the cloud that has potentially been impacted by HII regions. The ALMA-IMF catalog contains 700 cores that span a mass range of 0.15-250Msun at a typical size of 2100au. We show that this core sample has no significant distance bias and can be used to build core mass functions at similar physical scales. Significant gas motions, which we highlight here in the G353.41 region, are traced down to core scales and can be used to look for inflowing gas streamers and to quantify the impact of the possible associated core mass growth on the shape of the CMF with time. Our first analysis does not reveal any significant evolution of the matter concentration from clouds to cores or from the youngest to more evolved protoclusters, indicating that cloud dynamical evolution and stellar feedback have for the moment only had a slight effect on the structure of high-density gas in our sample. Furthermore, the first-look analysis of the line richness toward bright cores indicates that the survey encompasses several tens of hot cores, of which we highlight the most massive in the G351.77 cloud. Their homogeneous characterization can be used to constrain the emerging molecular complexity in protostars of high to intermediate masses.

preprint2021arXiv

The SEDIGISM survey: a search for molecular outflows

Context. The formation processes of massive stars are still unclear but a picture is emerging involving accretion disks and molecular outflows in what appears to be a scaled-up version of low-mass star formation. A census of outflow activity towards massive star-forming clumps in various evolutionary stages has the potential to shed light on massive star formation (MSF). Aims. We conducted an outflow survey towards ATLASGAL clumps using SEDIGISM data and aimed to obtain a large sample of clumps exhibiting outflows in different evolutionary stages. Methods. We identify the high-velocity wings of the 13CO lines toward ATLASGAL clumps by (1) extracting the simultaneously observed 13CO and C18O spectra from SEDIGISM, and (2) subtracting Gaussian fits to the scaled C18O from the 13CO, line after considering opacity broadening. Results. We have detected high-velocity gas towards 1192 clumps out of a total sample of 2052, giving an overall detection rate of 58%. Outflow activity has been detected in the earliest quiescent clumps (i.e., 70$μ$m weak), to the most evolved HII region stages i.e., 8$μ$m bright with MSF tracers. The detection rate increases as a function of evolution (quiescent=51%, protostellar=47%, YSO=57%, UCHII regions=76%). Conclusion. Our sample is the largest outflow sample identified so far. The high-detection rate from this large sample is consistent with previous results and supports that outflows are a ubiquitous feature of MSF. The lower detection rate in early evolutionary stages may be due to that outflows in the early stages are weak and difficult to detect. We obtain a statistically significant sample of outflow clumps for every evolutionary stage, especially for outflow clumps in the 70$μ$m dark stage. This suggests that the absence of 70$μ$m emission is not a robust indicator of starless/pre-stellar cores.

preprint2021arXiv

The SEDIGISM survey: The influence of spiral arms on the molecular gas distribution of the inner Milky Way

The morphology of the Milky Way is still a matter of debate. In order to shed light on uncertainties surrounding the structure of the Galaxy, in this paper, we study the imprint of spiral arms on the distribution and properties of its molecular gas. To do so, we take full advantage of the SEDIGISM survey that observed a large area of the inner Galaxy in the $^{13}$CO(2-1) line at an angular resolution of 28". We analyse the influences of the spiral arms by considering the features of the molecular gas emission as a whole across the longitude-velocity map built from the full survey. Additionally, we examine the properties of the molecular clouds in the spiral arms compared to the properties of their counterparts in the inter-arm regions. Through flux and luminosity probability distribution functions, we find that the molecular gas emission associated with the spiral arms does not differ significantly from the emission between the arms. On average, spiral arms show masses per unit length of $\sim10^5-10^6$ M$_{\odot} $kpc$^{-1}$. This is similar to values inferred from data sets in which emission distributions were segmented into molecular clouds. By examining the cloud distribution across the Galactic plane, we infer that the molecular mass in the spiral arms is a factor of 1.5 higher than that of the inter-arm medium, similar to what is found for other spiral galaxies in the local Universe. We observe that only the distributions of cloud mass surface densities and aspect ratio in the spiral arms show significant differences compared to those of the inter-arm medium; other observed differences appear instead to be driven by a distance bias. By comparing our results with simulations and observations of nearby galaxies, we conclude that the measured quantities would classify the Milky Way as a flocculent spiral galaxy, rather than as a grand-design one.

preprint2020arXiv

Dense gas formation in the Musca filament due to the dissipation of a supersonic converging flow

Observations with the Herschel Space Telescope have established that most of the star forming gas is organised in interstellar filaments, a finding that is supported by numerical simulations of the supersonic interstellar medium (ISM) where dense filamentary structures are ubiquitous. We aim to understand the formation of these dense structures by performing observations covering the $^{12}$CO(4-3), $^{12}$CO(3-2), and various CO(2-1) isotopologue lines of the Musca filament, using the APEX telescope. The observed CO intensities and line ratios cannot be explained by PDR (photodissociation region) emission because of the low ambient far-UV field that is strongly constrained by the non-detections of the [C II] line at 158 $μ$m and the [O I] line at 63 $μ$m, observed with the upGREAT receiver on SOFIA, as well as a weak [C I] 609 $μ$m line detected with APEX. We propose that the observations are consistent with a scenario in which shock excitation gives rise to warm and dense gas close to the highest column density regions in the Musca filament. Using shock models, we find that the CO observations can be consistent with excitation by J-type low-velocity shocks. A qualitative comparison of the observed CO spectra with synthetic observations of dynamic filament formation simulations shows a good agreement with the signature of a filament accretion shock that forms a cold and dense filament from a converging flow. The Musca filament is thus found to be dense molecular post-shock gas. Filament accretion shocks that dissipate the supersonic kinetic energy of converging flows in the ISM may thus play a prominent role in the evolution of cold and dense filamentary structures.

preprint2020arXiv

FEEDBACK: a SOFIA Legacy Program to Study Stellar Feedback in Regions of Massive Star Formation

FEEDBACK is a SOFIA legacy program dedicated to study the interaction of massive stars with their environment. It performs a survey of 11 galactic high mass star forming regions in the 158 $μ$m (1.9 THz) line of CII and the 63 $μ$m (4.7 THz) line of OI. We employ the 14 pixel LFA and 7 pixel HFA upGREAT instrument to spectrally resolve (0.24 MHz) these FIR structure lines. With an observing time of 96h, we will cover $\sim$6700 arcmin$^2$ at 14.1$''$ angular resolution for the CII line and 6.3$''$ for the OI line. The observations started in spring 2019 (Cycle 7). Our aim is to understand the dynamics in regions dominated by different feedback processes from massive stars such as stellar winds, thermal expansion, and radiation pressure, and to quantify the mechanical energy injection and radiative heating efficiency. The CII line provides the kinematics of the gas and is one of the dominant cooling lines of gas for low to moderate densities and UV fields. The OI line traces warm and high-density gas, excited in photodissociations regions with a strong UV field or by shocks. The source sample spans a broad range in stellar characteristics from single OB stars, to small groups of O stars, to rich young stellar clusters, to ministarburst complexes. It contains well-known targets such as Aquila, the Cygnus X region, M16, M17, NGC7538, NGC6334, Vela, and W43 as well as a selection of HII region bubbles, namely RCW49, RCW79, and RCW120. These CII maps, together with the less explored OI 63 $μ$m line, provide an outstanding database for the community. They will be made publically available and will trigger further studies and follow-up observations.

preprint2020arXiv

Formation of the Musca filament: Evidence for asymmetries in the accretion flow due to a cloud-cloud collision

Context. Dense molecular filaments are ubiquituous in the interstellar medium, yet their internal physical conditions and formation mechanism remain debated. Aims. We study the kinematics and physical conditions in the Musca filament and the Chamaeleon-Musca complex to constrain the physics of filament formation. Methods. We produced CO(2-1) isotopologue maps with the APEX telescope that cut through the Musca filament. We further study a NANTEN2 $^{12}$CO(1-0) map of the Musca cloud and the HI emission of the Chamaeleon-Musca complex. Results. The Musca cloud contains multiple velocity components. Radiative transfer modelling of the CO emission indicates that the Musca filament consists of a cold ($\sim$10 K), dense (n$_{H_2}\sim$10$^4$ cm$^{-3}$) crest, which is best described with a cylindrical geometry. Connected to the crest, a separate gas component at T$\sim$15 K and n$_{H_2}\sim$10$^3$ cm$^{-3}$ is found, the so-called strands. The filament crest has a transverse velocity gradient that is linked to the kinematics of the nearby ambient cloud. Studying the large scale kinematics, we find coherence of the asymmetric kinematics from the 50 pc HI cloud down to the Musca filament. We also report a strong [C$^{18}$O]/[$^{13}$CO] abundance drop by an order of magnitude from the filament crest to the strands over a distance $<$ 0.2 pc in a weak far-ultraviolet (FUV) field. Conclusions. The dense Musca filament crest is a long-lived (several crossing times), dynamic structure that can form stars in the near future because of continuous mass accretion. This mass accretion appears to be triggered by a HI cloud-cloud collision, which bends the magnetic field around dense filaments. This bending of the magnetic field is then responsible for the observed asymmetric accretion scenario of the Musca filament, which is, for instance, seen as a V-shape in the position-velocity (PV) diagram.

preprint2020arXiv

The accretion history of high-mass stars: An ArTéMiS pilot study of Infrared Dark Clouds

The mass growth of protostars is a central element to the determination of fundamental stellar population properties such as the initial mass function. Constraining the accretion history of individual protostars is therefore an important aspect of star formation research. The goal of the study presented here is to determine whether high-mass (proto)stars gain their mass from a compact (<0.1pc) fixed-mass reservoir of gas, often referred to as dense cores, in which they are embedded, or whether the mass growth of high-mass stars is governed by the dynamical evolution of the parsec-scale clump that typically surrounds them. To achieve this goal, we performed a 350micron continuum mapping of 11 infrared dark clouds, along side some of their neighbouring clumps, with the ArTéMiS camera on APEX. By identifying about 200 compact ArTéMiS sources, and matching them with Herschel Hi-GAL 70micron sources, we have been able to produce mass vs. temperature diagrams. We compare the nature (i.e. starless or protostellar) and location of the ArTéMiS sources in these diagrams with modelled evolutionary tracks of both core-fed and clump-fed accretion scenarios. We argue that the latter provide a better agreement with the observed distribution of high-mass star-forming cores. However, a robust and definitive conclusion on the question of the accretion history of high-mass stars requires larger number statistics.

preprint2020arXiv

The role of Galactic HII regions in the formation of filaments. High-resolution submilimeter imaging of RCW 120 with ArTéMiS

Massive stars and their associated ionized (HII) regions could play a key role in the formation and evolution of filaments that host star formation. However, the properties of filaments that interact with H regions are still poorly known. To investigate the impact of HII regions on the formation of filaments, we imaged the Galactic HII region RCW 120 and its surroundings where active star formation takes place and where the role of ionization feedback on the star formation process has already been studied. We used the ArTéMiS camera on the APEX telescope and combined the ArTéMiS data at 350 and 450 microns with Herschel-SPIRE/HOBYS. We studied the dense gas distribution around RCW 120 with a resolution of 8 arcsec (0.05 pc at a distance of 1.34 kpc). Our study allows us to trace the median radial intensity profile of the dense shell of RCW 120. This profile is asymmetric, indicating a clear compression from the HII region on the inner part of the shell. The profile is observed to be similarly asymmetric on both lateral sides of the shell, indicating a homogeneous compression over the surface. On the contrary, the profile analysis of a radial filament associated with the shell, but located outside of it, reveals a symmetric profile, suggesting that the compression from the ionized region is limited to the dense shell. The mean intensity profile of the internal part of the shell is well fitted by a Plummer like profile with a deconvolved Gaussian FWHM of 0.09 pc, as observed for filaments in low-mass star-forming regions. This study suggests that compression exerted by HII regions may play a key role in the formation of filaments and may further act on their hosted star formation. ArTéMiS data also suggest that RCW 120 might be a 3D ring, rather than a spherical structure

preprint2016arXiv

ATLASGAL -- A Galaxy-wide sample of dense filamentary structures

[Abridged] Aims. We study the properties of filamentary structures from the ATLASGAL survey. Methods. We use the DisPerSE algorithm to identify spatially coherent structures located across the inner-Galaxy (300 < l < 60 and |b| < 1.5). Results. We have determined distances, masses and physical sizes for 241 of the filamentary structures. We find a median distance of 3.8 kpc, a mean mass of a few 10^3 m_sun, a mean length of ~6pc and a mass-to-length ratio of (M/L) ~200-2000M_sun/ pc. We also find that these filamentary structures are tightly correlated with the spiral arms in longitude and velocity, and that their semi-major axis is preferentially aligned parallel to the Galactic mid-plane and therefore with the direction of large-scale Galactic magnetic field. We find many examples where the dense filaments identified in ATLASGAL are associated with larger scale filamentary structures (~100 pc), and argue that this is likely to be common, and as such these may indicate a connection between large-scale Galactic dynamics and star formation. Conclusions. We have produced a large and Galaxy-wide catalogue of dense filamentary structures that are representative of a particular size and mass range not previously well studied in the literature. Analyses of the properties and distribution of these filaments reveals that they are correlated with the spiral arms and make a significant contribution to star formation in the Galaxy. Massive star formation is ongoing within ~20% of the filaments and is strongly correlated with the filaments with the largest mass-to- length ratios. The luminosity of the embedded sources has a similar distribution to the Galactic-wide samples of young massive stars and can therefore be considered to be representative.

preprint2016arXiv

New debris disks in nearby young moving groups

A significant fraction of nearby young moving group members harbor circumstellar debris dust disks. Due to their proximity and youth, these disks are attractive targets for studying the early evolution of debris dust and planetesimal belts. Here we present 70 and 160$μ$m observations of 31 systems in the $β$ Pic moving group, and in the Tucana-Horologium, Columba, Carina and Argus associations, using the Herschel Space Observatory. None of these stars were observed at far-infrared wavelengths before. Our Herschel measurements were complemented by photometry from the WISE satellite for the whole sample, and by submillimeter/millimeter continuum data for one source, HD 48370. We identified six stars with infrared excess, four of them are new discoveries. By combining our new findings with results from the literature, we examined the incidence and general characteristics of debris disks around Sun-like members of the selected groups. With their dust temperatures of <45 K the newly identified disks around HD 38397, HD 48370, HD 160305, and BD-20 951 represent the coldest population within this sample. For HD 38397 and HD 48370, the emission is resolved in the 70$μ$m PACS images, the estimated radius of these disks is ~90 au. Together with the well-known disk around HD 61005, these three systems represent the highest mass end of the known debris disk population around young G-type members of the selected groups. In terms of dust content, they resemble the hypothesized debris disk of the ancient Solar System.

preprint2016arXiv

The Milky Way Project and ATLASGAL: The distribution and physical properties of cold clumps near infrared bubbles

We present a statistical study of the distribution and physical properties of cold dense material in and around the inner Galactic Plane near infrared bubbles as catalogued by the Milky Way Project citizen scientists. Using data from the ATLASGAL 870 um survey, we show that 48 +/- 2% of all cold clumps in the studied survey region (|l| <= 65 degrees, |b| <= 1 degree) are found in close proximity to a bubble, and 25 +/- 2% appear directly projected towards a bubble rim. A two-point correlation analysis confirms the strong correlation of massive cold clumps with expanding bubbles. It shows an overdensity of clumps along bubble rims that grows with increasing bubble size, which shows how interstellar medium material is reordered on large scales by bubble expansion around regions of massive star formation. The highest column density clumps appear resistent to the expansion, remaining overdense towards the bubbles' interior rather than being swept up by the expanding edge. Spectroscopic observations in ammonia show that cold dust clumps near bubbles appear to be denser, hotter and more turbulent than those in the field, offering circumstantial evidence that bubble-associated clumps are more likely to be forming stars. These observed differences in physical conditions persist for beyond the region of the bubble rims.

preprint2016arXiv

The physical and chemical structure of Sagittarius B2, I. Three-dimensional thermal dust and free-free continuum modeling on 100 au to 45 pc scales

We model the dust and free-free continuum emission in the high-mass star-forming region Sagittarius B2 in order to reconstruct the three-dimensional density and dust temperature distribution, as a crucial input to follow-up studies of the gas velocity field and molecular abundances. We employ the three-dimensional radiative transfer program RADMC-3D to calculate the dust temperature self-consistently, provided a given initial density distribution. This density distribution of the entire cloud complex is then recursively reconstructed based on available continuum maps, including both single-dish and high-resolution interferometric maps covering a wide frequency range (40 GHz - 4 THz). The model covers spatial scales from 45 pc down to 100 au, i.e. a spatial dynamic range of 10^5. We find that the density distribution of Sagittarius B2 can be reasonably well fitted by applying a superposition of spherical cores with Plummer-like density profiles. In order to reproduce the spectral energy distribution, we position Sgr B2(N) along the line of sight behind the plane containing Sgr B2(M). We find that the entire cloud complex comprises a total gas mass of 8.0 x 10^6 Msun within a diameter of 45 pc, corresponding to an averaged gas density of 170 Msun/pc^3. We estimate stellar masses of 2400 Msun and 20700 Msun and luminosities of 1.8 x 10^6 Lsun and 1.2 x 10^7 Lsun for Sgr B2(N) and Sgr B2(M), respectively. We report H_2 column densities of 2.9 x 10^24 cm^-2 for Sgr B2(N) and 2.5 x 10^24 cm^-2 for Sgr B2(M) in a 40" beam. For Sgr B2(S), we derive a stellar mass of 1100 Msun, a luminosity of 6.6 x 10^5 Lsun and a H_2 column density of 2.2 x 10^24 cm^-2 in a 40" beam. We calculate a star formation efficiency of 5% for Sgr B2(N) and 50% for Sgr B2(M), indicating that most of the gas content in Sgr B2(M) has already been converted to stars or dispersed.

preprint2016arXiv

The rate and latency of star formation in dense, massive clumps in the Milky Way

Newborn stars form within the localized, high density regions of molecular clouds. The sequence and rate at which stars form in dense clumps and the dependence on local and global environments are key factors in developing descriptions of stellar production in galaxies. We seek to observationally constrain the rate and latency of star formation in dense massive clumps that are distributed throughout the Galaxy and to compare these results to proposed prescriptions for stellar production. A sample of 24 micron-based Class~I protostars are linked to dust clumps that are embedded within molecular clouds selected from the APEX Telescope Large Area Survey of the Galaxy. We determine the fraction of star-forming clumps, f*, that imposes a constraint on the latency of star formation in units of a clump's lifetime. Protostellar masses are estimated from models of circumstellar environments of young stellar objects from which star formation rates are derived. Physical properties of the clumps are calculated from 870 micron dust continuum emission and NH_3 line emission. Linear correlations are identified between the star formation rate surface density, Sigma_{SFR}, and the quantities Sigma_{H2}/tau_{ff} and Sigma_{H2}/tau_{cross}, suggesting that star formation is regulated at the local scales of molecular clouds. The measured fraction of star forming clumps is 23%. Accounting for star formation within clumps that are excluded from our sample due to 24 micron saturation, this fraction can be as high as 31%. Dense, massive clumps form primarily low mass (< 1-2 msun) stars with emergent 24 micron fluxes below our sensitivity limit or are incapable of forming any stars for the initial 70% of their lifetimes. The low fraction of star forming clumps in the Galactic center relative to those located in the disk of the Milky Way is verified.

preprint2016arXiv

The upGREAT 1.9 THz multi-pixel high resolution spectrometer for the SOFIA Observatory

We present a new multi-pixel high resolution (R >10^7) spectrometer for the Stratospheric Observatory for Far-Infrared Astronomy (SOFIA). The receiver uses 2 x 7-pixel subarrays in orthogonal polarization, each in an hexagonal array around a central pixel. We present the first results for this new instrument after commissioning campaigns in May and December 2015 and after science observations performed in May 2016 . The receiver is designed to ultimately cover the full 1.8-2.5 THz frequency range but in its first implementation, the observing range was limited to observations of the [CII] line at 1.9 THz in 2015 and extended to 1.83-2.07 THz in 2016. The instrument sensitivities are state-of-the-art and the first scientific observations performed shortly after the commissioning confirm that the time efficiency for large scale imaging is improved by more than an order of magnitude as compared to single pixel receivers. An example of large scale mapping around the Horsehead Nebula is presented here illustrating this improvement. The array has been added to SOFIA's instrument suite already for ongoing observing cycle 4.

preprint2016arXiv

Understanding star formation in molecular clouds III. Probability distribution functions of molecular lines in Cygnus X

Column density (N) PDFs serve as a powerful tool to characterize the physical processes that influence the structure of molecular clouds. Star-forming clouds can best be characterized by lognormal PDFs for the lower N range and a power-law tail for higher N, commonly attributed to turbulence and self-gravity and/or pressure, respectively. We report here on PDFs obtained from observations of 12CO, 13CO, C18O, CS, and N2H+ in the Cygnus X North region and compare to a PDF derived from dust observations with the Herschel satellite. The PDF of 12CO is lognormal for Av~1-30, but is cut for higher Av due to optical depth effects. The PDFs of C18O and 13CO are mostly lognormal up for Av~1-15, followed by excess up to Av~40. Above that value, all CO PDFs drop, most likely due to depletion. The high density tracers CS and N2H+ exhibit only a power law distribution between Av~15 and 400, respectively. The PDF from dust is lognormal for Av~3-15 and has a power-law tail up to Av~500. Absolute values for the molecular line column densities are, however, rather uncertain due to abundance and excitation temperature variations. Taken the dust PDF face value, we 'calibrate' the molecular line PDF of CS to the one of the dust and determined an abundance [CS]/[H2] of 10^-9. The slopes of the power-law tails of the CS, N2H+, and dust PDFs are consistent with free-fall collapse of filaments and clumps. A quasi static configuration of filaments and clumps can possibly also account for the observed N-PDFs, as long as they have a sufficiently condensed density structure and external ram pressure by gas accretion is provided. The somehow flatter slopes of N2H+ and CS can reflect an abundance change and/or subthermal excitation at low column densities.

preprint2015arXiv

ATLASGAL - Kinematic distances and the dense gas mass distribution of the inner Galaxy

The formation of high mass stars and clusters occurs in giant molecular clouds. Objects in evolved stages of massive star formation such as protostars, hot molecular cores, and ultracompact HII regions have been studied in more detail than earlier, colder objects. With this in mind, the APEX Telescope Large Area Survey of the whole inner Galactic plane at 870 micron (ATLASGAL) has been carried out to provide a global view of cold dust and star formation at submillimetre wavelengths. To derive kinematic distances to a large sample of ATLASGAL clumps we divided them into groups of sources, which are located close together, mostly within a radius of 2 pc, and have velocities in a similar range with a median velocity dispersion of ~ 1 km/s. Using NH3, N2H+ and CS velocities we calculate near and far kinematic distances to 296 groups of ATLASGAL sources in the first quadrant and 393 groups in the fourth quadrant. We analyse HI self-absorption and HI absorption to resolve the kinematic distance ambiguity. We obtain a scale height of ~ 28+/-2 pc and displacement below the Galactic midplane of ~ -7+/-1 pc. Within distances from 2 to 18 kpc ATLASGAL clumps have a broad range of gas masses with a median of 1050 solar masses and a wide distribution of radii with a median of 0.4 pc. Their distribution in galactocentric radii is correlated with spiral arms. Using a statistically significant ATLASGAL sample we derive a power-law exponent of -2.2+/-0.1 of the clump mass function. This is consistent with the slope derived for clusters and with that of the stellar initial mass function. Examining the power-law index for different galactocentric distances and various source samples shows that it is independent of environment and evolutionary phase. Fitting the mass-size relationship by a power law gives a slope of 1.76+/-0.01 for cold sources such as IRDCs and warm clumps associated with HII regions.

preprint2015arXiv

ATLASGAL-selected massive clumps in the inner Galaxy, II: Characterisation of different evolutionary stages and their SiO emission

The processes leading to the birth of high-mass stars are poorly understood. We characterise here a sample of 430 massive clumps from the ATLASGAL survey, which are representative of different evolutionary stages. To establish a census of molecular tracers of their evolution we performed an unbiased spectral line survey covering the 3-mm atmospheric window between 84-117 GHz with the IRAM 30m. A smaller sample of 128 clumps has been observed in the SiO (5-4) transition with the APEX telescope to complement the SiO (2-1) line and probe the excitation conditions of the emitting gas, which is the main focus of the current study. We report a high detection rate of >75% of the SiO (2-1) line and a >90% detection rate from the dedicated follow-ups in the (5-4) transition. The SiO (2-1) line with broad line profiles and high detection rates, is a powerful probe of star formation activity, while the ubiquitous detection of SiO in all evolutionary stages suggests a continuous star formation process in massive clumps. We find a large fraction of infrared-quiet clumps to exhibit SiO emission, the majority of them only showing a low-velocity component (FWHM~5-6 km/s) centred at the rest velocity of the clump. In the current picture, where this is attributed to low-velocity shocks from cloud-cloud collisions, this can be used to pinpoint the youngest, thus, likely prestellar massive structures. Based on the line ratio of the (5-4) to the (2-1) line, our study reveals a trend of changing excitation conditions that lead to brighter emission in the (5-4) line towards more evolved sources. Our analysis delivers a more robust estimate of SiO column density and abundance than previous studies and questions the decrease of jet activity in massive clumps as a function of age.

preprint2015arXiv

Detection of two power-law tails in the probability distribution functions of massive GMCs

We report the novel detection of complex high-column density tails in the probability distribution functions (PDFs) for three high-mass star-forming regions (CepOB3, MonR2, NGC6334), obtained from dust emission observed with Herschel. The low column density range can be fit with a lognormal distribution. A first power-law tail starts above an extinction (Av) of ~6-14. It has a slope of alpha=1.3-2 for the rho~r^-alpha profile for an equivalent density distribution (spherical or cylindrical geometry), and is thus consistent with free-fall gravitational collapse. Above Av~40, 60, and 140, we detect an excess that can be fitted by a flatter power law tail with alpha>2. It correlates with the central regions of the cloud (ridges/hubs) of size ~1 pc and densities above 10^4 cm^-3. This excess may be caused by physical processes that slow down collapse and reduce the flow of mass towards higher densities. Possible are: 1. rotation, which introduces an angular momentum barrier, 2. increasing optical depth and weaker cooling, 3. magnetic fields, 4. geometrical effects, and 5. protostellar feedback. The excess/second power-law tail is closely linked to high-mass star-formation though it does not imply a universal column density threshold for the formation of (high-mass) stars.

preprint2015arXiv

Discovery of molecular gas around HD 131835 in an APEX molecular line survey of bright debris disks

Debris disks are considered to be gas-poor, but recent observations revealed molecular or atomic gas in several 10-40 Myr old systems. We used the APEX and IRAM 30m radiotelescopes to search for CO gas in 20 bright debris disks. In one case, around the 16 Myr old A-type star HD 131835, we discovered a new gas-bearing debris disk, where the CO 3-2 transition was successfully detected. No other individual system exhibited a measurable CO signal. Our Herschel Space Observatory far-infrared images of HD 131835 marginally resolved the disk both at 70 and 100$μ$m, with a characteristic radius of ~170 au. While in stellar properties HD 131835 resembles $β$ Pic, its dust disk properties are similar to those of the most massive young debris disks. With the detection of gas in HD 131835 the number of known debris disks with CO content has increased to four, all of them encircling young ($\leq$40 Myr) A-type stars. Based on statistics within 125 pc, we suggest that the presence of detectable amount of gas in the most massive debris disks around young A-type stars is a common phenomenon. Our current data cannot conclude on the origin of gas in HD 131835. If the gas is secondary, arising from the disruption of planetesimals, then HD 131835 is a comparably young and in terms of its disk more massive analogue of the $β$ Pic system. However, it is also possible that this system similarly to HD 21997 possesses a hybrid disk, where the gas material is predominantly primordial, while the dust grains are mostly derived from planetesimals.

preprint2015arXiv

First detection of CF+ towards a high-mass protostar

We report the first detection of the J = 1 - 0 (102.6 GHz) rotational lines of CF+ (fluoromethylidynium ion) towards CygX-N63, a young and massive protostar of the Cygnus X region. This detection occurred as part of an unbiased spectral survey of this object in the 0.8-3 mm range, performed with the IRAM 30m telescope. The data were analyzed using a local thermodynamical equilibrium model (LTE model) and a population diagram in order to derive the column density. The line velocity (-4 km s-1) and line width (1.6 km s-1) indicate an origin from the collapsing envelope of the protostar. We obtain a CF+ column density of 4.10e11 cm-2. The CF+ ion is thought to be a good tracer for C+ and assuming a ratio of 10e-6 for CF+/C+, we derive a total number of C+ of 1.2x10e53 within the beam. There is no evidence of carbon ionization caused by an exterior source of UV photons suggesting that the protostar itself is the source of ionization. Ionization from the protostellar photosphere is not efficient enough. In contrast, X-ray ionization from the accretion shock(s) and UV ionization from outflow shocks could provide a large enough ionizing power to explain our CF+ detection. Surprisingly, CF+ has been detected towards a cold, massive protostar with no sign of an external photon dissociation region (PDR), which means that the only possibility is the existence of a significant inner source of C+. This is an important result that opens interesting perspectives to study the early development of ionized regions and to approach the issue of the evolution of the inner regions of collapsing envelopes of massive protostars. The existence of high energy radiations early in the evolution of massive protostars also has important implications for chemical evolution of dense collapsing gas and could trigger peculiar chemistry and early formation of a hot core.

preprint2015arXiv

Infall through the evolution of high-mass star-forming clumps

With the GREAT receiver at the Stratospheric Observatory for Infrared Astronomy (SOFIA), nine massive molecular clumps have been observed in the ammonia $3_{2+}- 2_{2-}$ line at 1.8~THz in a search for signatures of infall. The sources were selected from the ATLASGAL submillimeter dust continuum survey of our Galaxy. Clumps with high masses covering a range of evolutionary stages based on their infrared properties were chosen. The ammonia line was detected in all sources, leading to five new detections and one confirmation of a previous detection of redshifted absorption in front of their strong THz continuum as a probe of infall in the clumps. These detections include two clumps embedded in infrared dark clouds. The measured velocity shifts of the absorptions compared to optically thin \CSEO\ (3--2) emission are 0.3--2.8~km/s, corresponding to fractions of 3\%\ to 30\% of the free-fall velocities of the clumps. The ammonia infall signature is compared with complementary data of different transitions of HCN, HNC, CS, and HCO$^+$, which are often used to probe infall via their blue-skewed line profiles. The best agreement with the ammonia results is found for the HCO$^+$ (4--3) transitions, but the latter is still strongly blended with emission from associated outflows. This outflow signature is far less prominent in the THz ammonia lines, which confirms it as a powerful probe of infall in molecular clumps. Infall rates in the range from 0.3 to 16~$10^{-3}\,M_\odot/$yr were derived with a tentative correlation with the virial parameters of the clumps. The new observations show that infall on clump scales is ubiquitous through a wide range of evolutionary stages, from $L/M$ covering about ten to several hundreds.

preprint2015arXiv

Infrared dark clouds on the far side of the Galaxy

Context: Infrared dark clouds are the coldest and densest portions of giant molecular clouds. The most massive ones represent some of the most likely birthplaces for the next generation of massive stars in the Milky Way. Because a strong mid-IR background is needed to make them appear in absorption, they are usually assumed to be nearby. Aims: We use THz absorption spectroscopy to solve the distance ambiguity associated with kinematic distances for the IR-dark clouds in the TOP100 ATLASGAL sample, a flux-limited selection of massive clumps in different evolutionary phases of star formation. Methods: The para-H2O ground state transition at 1113.343 GHz, observed with Herschel/HIFI, was used to investigate the occurrence of foreground absorption along the line of sight directly towards infrared-dark clouds. Additional consistency checks were performed using MALT90 and HiGAL archival data and targeted Mopra and APEX spectroscopic observations. Results: We report the first discovery of five IRDCs in the TOP100 lying conclusively at the far kinematic distance, showing that the mere presence of low-contrast mid-IR absorption is not sufficient to unequivocally resolve the near/far ambiguity in favour of the former. All IRDCs are massive and actively forming high-mass stars; four of them also show infall signatures. Conclusions: We give a first estimate of the fraction of dark sources at the far distance (~11% in the TOP100) and describe their appearance and properties. The assumption that all dark clouds lie at the near distance may lead, in some cases, to underestimating masses, sizes, and luminosities, possibly causing clouds to be missed that will form very massive stars and clusters.

preprint2015arXiv

The ATLASGAL survey: distribution of cold dust in the Galactic plane. Combination with Planck data

Sensitive ground-based submillimeter surveys, such as ATLASGAL, provide a global view on the distribution of cold dense gas in the Galactic plane. Here we use the 353 GHz maps from the Planck/HFI instrument to complement the ground-based APEX/LABOCA observations with information on larger angular scales. The resulting maps reveal the distribution of cold dust in the inner Galaxy with a larger spatial dynamic range. We find examples of elongated structures extending over angular scales of 0.5 degree. Corresponding to >30 pc structures in projection at a distance of 3 kpc, these dust lanes are very extended and show large aspect ratios. Furthermore, we assess the fraction of dense gas ($f_{\rm DG}$), and estimate 2-5% (above A$_{\rm{v}}>$7 mag) on average in the Galactic plane. PDFs of the column density reveal the typically observed log-normal distribution for low- and exhibit an excess at high column densities. As a reference for extragalactic studies, we show the line-of-sight integrated N-PDF of the inner Galaxy, and derive a contribution of this excess to the total column density of $\sim2.2$%, above $N_{\rm H_2} = 2.92\times10^{22}$ cm$^{-2}$. Taking the total flux density, we provide an independent estimate of the mass of molecular gas in the inner Galaxy of $\sim1\times10^9\,M_{\odot}$, which is consistent with previous estimates using CO emission. From the mass and $f_{\rm DG}$ we estimate a Galactic SFR of $\dot M = 1.3\,M_{\odot}$ yr$^{-1}$. While the distribution of diffuse gas is homogenous in the inner Galaxy, the CMZ stands out with a higher dense gas fraction. The low star formation efficiency of the Milky Way is well explained by the low $f_{\rm DG}$ in the Galactic ISM, while the high $f_{\rm DG}$ towards the CMZ, despite its low star formation activity, suggests that, in that particular region of our Galaxy, high-density gas is not the bottleneck for star formation.

preprint2015arXiv

The first Galaxy scale hunt for the youngest high-mass protostars

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

preprint2015arXiv

The RMS Survey: Ammonia mapping of the environment of massive young stellar objects

We present the results of ammonia observations towards 66 massive star forming regions identified by the Red MSX source survey. We have used the Green Bank Telescope and the K-band focal plane array to map the ammonia NH3 (1,1) and (2,2) inversion emission at a resolution of 30 arcsec in 8 arcmin regions towards the positions of embedded massive star formation. We have identified a total of 115 distinct clumps, approximately two-thirds of which are associated with an embedded massive young stellar object or compact HII region, while the others are classified as quiescent. There is a strong spatial correlation between the peak NH3 emission and the presence of embedded objects. We derive the spatial distribution of the kinetic gas temperatures, line widths, and NH$_3$ column densities from these maps, and by combining these data with dust emission maps we estimate clump masses, H$_2$ column densities and ammonia abundances. The clumps have typical masses of ~1000 Msun and radii ~0.5 pc, line widths of ~2 km/s and kinetic temperatures of ~16-20 K. We find no significant difference between the sizes and masses of the star forming and quiescent subsamples; however, the distribution maps reveal the presence of temperature and line width gradients peaking towards the centre for the star forming clumps while the quiescent clumps show relatively uniform temperatures and line widths throughout. Virial analysis suggests that the vast majority of clumps are gravitationally bound and are likely to be in a state of global free fall in the absence of strong magnetic fields. The similarities between the properties of the two subsamples suggest that the quiescent clumps are also likely to form massive stars in the future, and therefore provide a excellent opportunity to study the initial conditions of massive pre-stellar and protostellar clumps.

preprint2015arXiv

Understanding star formation in molecular clouds I. Effects of line-of-sight contamination on the column density structure

Column-density maps of molecular clouds are one of the most important observables in the context of molecular cloud- and star-formation (SF) studies. With the Herschel satellite it is now possible to determine the column density from dust emission. We use observations and simulations to demonstrate how LOS contamination affects the column density probability distribution function (PDF). We apply a first-order approximation (removing a constant level) to the molecular clouds of Auriga, Maddalena, Carina and NGC3603. In perfect agreement with the simulations, we find that the PDFs become broader, the peak shifts to lower column densities, and the power-law tail of the PDF flattens after correction. All PDFs have a lognormal part for low column densities with a peak at Av~2, a deviation point (DP) from the lognormal at Av(DP)~4-5, and a power-law tail for higher column densities. Assuming a density distribution rho~r^-alpha, the slopes of the power-law tails correspond to alpha(PDF)=1.8, 1.75, and 2.5 for Auriga, Carina, and NGC3603 (alpha~1.5-2 is consistent gravitational collapse). We find that low-mass and high-mass SF clouds display differences in the overall column density structure. Massive clouds assemble more gas in smaller cloud volumes than low-mass SF ones. However, for both cloud types, the transition of the PDF from lognormal shape into power-law tail is found at the same column density (at Av~4-5 mag). Low-mass and high-mass SF clouds then have the same low column density distribution, most likely dominated by supersonic turbulence. At higher column densities, collapse and external pressure can form the power-law tail. The relative importance of the two processes can vary between clouds and thus lead to the observed differences in PDF and column density structure.

preprint2015arXiv

Understanding star formation in molecular clouds II. Signatures of gravitational collapse of IRDCs

We analyse column density and temperature maps derived from Herschel dust continuum observations of a sample of massive infrared dark clouds (G11.11-0.12, G18.82-0.28, G28.37+0.07, G28.53-0.25). We disentangle the velocity structure of the clouds using 13CO 1-0 and 12CO 3-2 data, showing that these IRDCs are the densest regions in massive giant molecular clouds and not isolated features. The probability distribution function (PDF) of column densities for all clouds have a power-law distribution over all (high) column densities, regardless of the evolutionary stage of the cloud: G11.11-0.12, G18.82-0.28, and G28.37+0.07 contain (proto)-stars, while G28.53-0.25 shows no signs of star formation. This is in contrast to the purely log-normal PDFs reported for near/mid-IR extinction maps. We only find a log-normal distribution for lower column densities, if we perform PDFs of the column density maps of the whole GMC in which the IRDCs are embedded. By comparing the PDF slope and the radial column density profile, we attribute the power law to the effect of large-scale gravitational collapse and to local free-fall collapse of pre- and protostellar cores. Independent from the PDF analysis, we find infall signatures in the spectral profiles of 12CO for G28.37+0.07 and G11.11-0.12, supporting the scenario of gravitational collapse. IRDCs are the densest regions within GMCs, which may be the progenitors of massive stars or clusters. At least some of the IRDCs are probably the same features as ridges (high column density regions with N>1e23 cm-2 over small areas), which were defined for nearby IR-bright GMCs. Because IRDCs are only confined to the densest (gravity dominated) cloud regions, the PDF constructed from this kind of a clipped image does not represent the (turbulence dominated) low column density regime of the cloud.

preprint2014arXiv

ATLASGAL --- Complete compact source catalogue: 280\degr\ $ <\ell <$ 60\degr

The APEX Telescope Large Area Survey of the Galaxy (ATLASGAL) is the largest and most sensitive systematic survey of the inner Galactic plane in the submillimetre wavelength regime. The observations were carried out with the Large APEX Bolometer Camera (LABOCA), an array of 295 bolometers observing at 870\,$μ$m (345 GHz). Aim: In this research note we present the compact source catalogue for the 280\degr\ $ <\ell <$ 330\degr\ and 21\degr\ $ <\ell <$ 60\degr\ regions of this survey. Method: The construction of this catalogue was made with the source extraction routine \sex\ using the same input parameters and procedures used to analyse the inner Galaxy region presented in an earlier publication (i.e., 330\degr\ $ <\ell <$ 21\degr). Results: We have identified 3523 compact sources and present a catalogue of their properties. When combined with the regions already published this provides a comprehensive and unbiased database of ~10163 massive, dense clumps located across the inner Galaxy.

preprint2014arXiv

ATLASGAL --- towards a complete sample of massive star forming clumps

By matching infrared-selected, massive young stellar objects (MYSOs) and compact HII regions in the RMS survey to massive clumps found in the submillimetre ATLASGAL survey, we have identified ~1000 embedded young massive stars between 280\degr < $\ell$ < 350\degr and 10degr < $\ell$ < 60\degr with |b|<1.5degr. Combined with an existing sample of radio-selected methanol masers and compact HII regions, the result is a catalogue of ~1700 massive stars embedded within ~1300 clumps located across the inner Galaxy, containing three observationally distinct subsamples, methanol-maser, MYSO and HII-region associations, covering the most important tracers of massive star formation, thought to represent key stages of evolution. We find that massive star formation is strongly correlated with the regions of highest column density in spherical, centrally condensed clumps. We find no significant differences between the three samples in clump structure or the relative location of the embedded stars, which suggests that the structure of a clump is set before the onset of star formation, and changes little as the embedded object evolves towards the main sequence. There is a strong linear correlation between clump mass and bolometric luminosity, with the most massive stars forming in the most massive clumps. We find that the MYSO and HII-region subsamples are likely to cover a similar range of evolutionary stages and that the majority are near the end of their main accretion phase. We find few infrared-bright MYSOs associated with the most massive clumps, probably due to very short pre-main sequence lifetimes in the most luminous sources.

preprint2014arXiv

ATLASGAL-selected massive clumps in the inner Galaxy: I. CO depletion and isotopic ratios

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.

preprint2014arXiv

Ionization compression impact on dense gas distribution and star formation, Probability density functions around H ii regions as seen by Herschel

Ionization feedback should impact the probability distribution function (PDF) of the column density around the ionized gas. We aim to quantify this effect and discuss its potential link to the Core and Initial Mass Function (CMF/IMF). We used in a systematic way Herschel column density maps of several regions observed within the HOBYS key program: M16, the Rosette and Vela C molecular cloud, and the RCW 120 H ii region. We fitted the column density PDFs of all clouds with two lognormal distributions, since they present a double-peak or enlarged shape in the PDF. Our interpretation is that the lowest part of the column density distribution describes the turbulent molecular gas while the second peak corresponds to a compression zone induced by the expansion of the ionized gas into the turbulent molecular cloud. The condensations at the edge of the ionized gas have a steep compressed radial profile, sometimes recognizable in the flattening of the power-law tail. This could lead to an unambiguous criterion able to disentangle triggered from pre-existing star formation. In the context of the gravo-turbulent scenario for the origin of the CMF/IMF, the double peaked/enlarged shape of the PDF may impact the formation of objects at both the low-mass and the high-mass end of the CMF/IMF. In particular a broader PDF is required by the gravo-turbulent scenario to fit properly the IMF with a reasonable initial Mach number for the molecular cloud. Since other physical processes (e.g. the equation of state and the variations among the core properties) have already been suggested to broaden the PDF, the relative importance of the different effects remains an open question.

preprint2014arXiv

Methanol Maser Associated Outflows: Detection statistics and properties

We have selected the positions of 54 6.7GHz methanol masers from the Methanol Multibeam Survey catalogue, covering a range of longitudes between $20^{\circ}$ and $34^{\circ}$ of the Galactic Plane. These positions were mapped in the J=3-2 transition of both the $\rm{^{13}CO}$ and $\rm{C^{18}O}$ lines. A total of 58 $\rm{^{13}CO}$ emission peaks are found in the vicinity of these maser positions. We search for outflows around all $\rm{^{13}CO}$ peaks, and find evidence for high-velocity gas in all cases, spatially resolving the red and blue outflow lobes in 55 cases. Of these sources, 44 have resolved kinematic distances, and are closely associated with the 6.7GHz masers, a sub-set referred to as Methanol Maser Associated Outflows (MMAOs). We calculate the masses of the clumps associated with each peak using 870 $\rm{μm}$ continuum emission from the ATLASGAL survey. A strong correlation is seen between the clump mass and both outflow mass and mechanical force, lending support to models in which accretion is strongly linked to outflow. We find that the scaling law between outflow activity and clump masses observed for low-mass objects, is also followed by the MMAOs in this study, indicating a commonality in the formation processes of low-mass and high-mass stars.

preprint2014arXiv

SiO emission from low- and high-velocity shocks in Cygnus-X massive dense clumps

We used PdBI observations of SiO (2-1) to investigate the morphology and profile of the SiO emission within several massive dense clumps (MDCs) in Cygnus-X. We find that most molecular outflows are detected in both SiO and CO, although there are some cases of CO outflows with no SiO counterpart. We find a significant amount of narrow line SiO emission that appears to be unrelated to outflows. The fraction of the total SiO luminosity that is not associated with outflows is highly variable in the different MDCs (from 10% to 90%); this might be a problem when extrapolating outflow properties from SiO luminosities without resolving individual outflows. The extent of the narrow SiO emission varies from rather compact (~ 0.03 pc) to widespread (~0.2 pc), and its kinematics often differs from those found by other high-density tracers such as H13CO+. We find that the least centrally concentrated clumps with the least massive protostellar cores have the most widespread narrow SiO emission. In line with previous evidence of SiO emission associated with low-velocity shocks, we propose an evolutionary picture to explain the existence and distribution of narrow SiO line profiles. In this scenario, the least centrally condensed MDCs are at an early stage where the SiO emission traces shocks from the large-scale collapse of material onto the MDC (e.g. CygX-N40). As the MDC collapses, the SiO emission becomes more confined to the close surroundings of cores, tracing the post-shock material from the infalling MDC against the dense cores (e.g. CygX-N3, N12, and N48). At later stages, when single massive protostars are formed, the SiO luminosity is largely dominated by powerful outflows, and the weaker narrow component shows perhaps the last remnants of the initial collapse (e.g. CygX-N53 and N63).

preprint2014arXiv

SiO excitation from dense shocks in the earliest stages of massive star formation

Molecular outflows are a direct consequence of accretion, and therefore they represent one of the best tracers of accretion processes in the still poorly understood early phases of high-mass star formation. Previous studies suggested that the SiO abundance decreases with the evolution of a massive young stellar object probably because of a decay of jet activity, as witnessed in low-mass star-forming regions. We investigate the SiO excitation conditions and its abundance in outflows from a sample of massive young stellar objects through observations of the SiO(8-7) and CO(4-3) lines with the APEX telescope. Through a non-LTE analysis, we find that the excitation conditions of SiO increase with the velocity of the emitting gas. We also compute the SiO abundance through the SiO and CO integrated intensities at high velocities. For the sources in our sample we find no significant variation of the SiO abundance with evolution for a bolometric luminosity-to-mass ratio of between 4 and 50 $L_\odot/M_\odot$. We also find a weak increase of the SiO(8-7) luminosity with the bolometric luminosity-to-mass ratio. We speculate that this might be explained with an increase of density in the gas traced by SiO. We find that the densities constrained by the SiO observations require the use of shock models that include grain-grain processing. For the first time, such models are compared and found to be compatible with SiO observations. A pre-shock density of $10^5\, $cm$^{-3}$ is globally inferred from these comparisons. Shocks with a velocity higher than 25 km s$^{-1}$ are invoked for the objects in our sample where the SiO is observed with a corresponding velocity dispersion. Our comparison of shock models with observations suggests that sputtering of silicon-bearing material (corresponding to less than 10% of the total silicon abundance) from the grain mantles is occurring.

preprint2014arXiv

The almost ubiquitous association of 6.7 GHz methanol masers with dust

We report the results of 870-$μ$m continuum observations, using the Large APEX Bolometer Camera (LABOCA), towards 77 class-II, 6.7-GHz methanol masers identified by the Methanol Multibeam (MMB) survey to map the thermal emission from cool dust towards these objects. These data complement a study of 630 methanol masers associated with compact dense clumps identified from the ATLASGAL survey. Compact dust emission is detected towards 70 sources, which implies a dust-association rate of 99% for the full MMB catalogue. Evaluation of the derived dust and maser properties leads us to conclude that the combined sample represents a single population tracing the same phenomenon. We find median clump masses of a few 10$^3$ M$\odot$ and that all but a handful of sources satisfy the mass-size criterion required for massive star formation. This study provides the strongest evidence of the almost ubiquitous association of methanol masers with massive, star-forming clumps. The fraction of methanol-maser associated clumps is a factor of ~2 lower in the outer Galaxy than the inner Galaxy, possibly a result of the lower metallicity environment of the former. We find no difference in the clump-mass and maser-luminosity distributions of the inner and outer Galaxy. The maser-pumping and clump-formation mechanisms are therefore likely to be relatively invariant to Galactic location. Finally, we use the ratio of maser luminosity and clump mass to investigate the hypothesis that the maser luminosity is a good indicator of the evolutionary stage of the embedded source, however, we find no evidence to support this.

preprint2014arXiv

The ATLASGAL survey: a catalog of dust condensations in the Galactic plane

The formation processes and the evolutionary stages of high-mass stars are poorly understood compared to low-mass stars. Large-scale surveys are needed to provide an unbiased census of high column density sites which can potentially host precursors to high-mass stars. Here we use the ATLASGAL survey covering 420 sq. degree of the Galactic plane at 870 $μ$m; and use the MRE-GLC method to identify the population of embedded sources throughout the inner Galaxy. We identify in total 10861 compact sub-millimeter sources with fluxes above 5 sigma. Completeness tests show that our catalogue is 97% complete above 5 sigma and >99% complete above 7$σ$. We correlate this sample with mid-infrared point source catalogues (MSX at 21.3 $μ$m and WISE at 22 $μ$m) and determine a lower limit of ~33% that are associated with embedded protostellar objects. We note that the proportion of clumps associated with mid-infrared sources increases with increasing flux density, achieving a rather constant fraction of ~75% of all clumps with fluxes over 5 Jy/beam being associated with star-formation. Examining the source counts as a function of Galactic longitude we are able to identify the most prominent star forming regions in the Galaxy. From the fraction of the likely massive quiescent clumps (~25%) we estimate a formation time-scale of ~7.5+/-2.5 $\times$ 10$^4$yr for the deeply embedded phase before the emergence of luminous YSOs. Such a short duration for the formation of high-mass stars in massive clumps clearly proves that the earliest phases have to be dynamic with supersonic motions.

preprint2014arXiv

Water emission from the high-mass star-forming region IRAS 17233-3606. High water abundances at high velocities

We investigate the physical and chemical processes at work during the formation of a massive protostar based on the observation of water in an outflow from a very young object previously detected in H2 and SiO in the IRAS 17233-3606 region. We estimated the abundance of water to understand its chemistry, and to constrain the mass of the emitting outflow. We present new observations of shocked water obtained with the HIFI receiver onboard Herschel. We detected water at high velocities in a range similar to SiO. We self-consistently fitted these observations along with previous SiO data through a state-of-the-art, one-dimensional, stationary C-shock model. We found that a single model can explain the SiO and H2O emission in the red and blue wings of the spectra. Remarkably, one common area, similar to that found for H2 emission, fits both the SiO and H2O emission regions. This shock model subsequently allowed us to assess the shocked water column density, N(H2O)=1.2x10^{18} cm^{-2}, mass, M(H2O)=12.5 M_earth, and its maximum fractional abundance with respect to the total density, x(H2O)=1.4x10^{-4}. The corresponding water abundance in fractional column density units ranges between 2.5x10^{-5} and 1.2x10^{-5}, in agreement with recent results obtained in outflows from low- and high-mass young stellar objects.

preprint2013arXiv

A resolved debris disk around the candidate planet-hosting star HD95086

Recently, a new planet candidate was discovered on direct images around the young (10-17 Myr) A-type star HD95086. The strong infrared excess of the system indicates that, similarly to HR8799, ß Pic, and Fomalhaut, the star harbors a circumstellar disk. Aiming to study the structure and gas content of the HD95086 disk, and to investigate its possible interaction with the newly discovered planet, here we present new optical, infrared and millimeter observations. We detected no CO emission, excluding the possibility of an evolved gaseous primordial disk. Simple blackbody modeling of the spectral energy distribution suggests the presence of two spatially separate dust belts at radial distances of 6 and 64 AU. Our resolved images obtained with the Herschel Space Observatory reveal a characteristic disk size of ~6.0x5.4 arcsec (540x490 AU) and disk inclination of ~25 degree. Assuming the same inclination for the planet candidate's orbit, its re-projected radial distance from the star is 62 AU, very close to the blackbody radius of the outer cold dust ring. The structure of the planetary system at HD95086 resembles the one around HR8799. Both systems harbor a warm inner dust belt and a broad colder outer disk and giant planet(s) between the two dusty regions. Modelling implies that the candidate planet can dynamically excite the motion of planetesimals even out to 270 AU via their secular perturbation if its orbital eccentricity is larger than about 0.4. Our analysis adds a new example to the three known systems where directly imaged planet(s) and debris disks co-exist.

preprint2013arXiv

ALMA continuum observations of a 30 Myr old gaseous debris disk around HD 21997

Circumstellar disks around stars older than 10 Myr are expected to be gas-poor. There are, however, two examples of old (30-40 Myr) debris-like disks containing a detectable amount of cold CO gas. Here we present ALMA and Herschel Space Observatory observations of one of these disks, around HD 21997, and study the distribution and origin of the dust and its connection to the gas. Our ALMA continuum images at 886um clearly resolve a broad ring of emission within a diameter of ~4.5 arcsec, adding HD 21997 to the dozen debris disks resolved at (sub)millimeter wavelengths. Modeling the morphology of the ALMA image with a radiative transfer code suggests inner and outer radii of ~55 and ~150 AU, and a dust mass of 0.09 M_Earth. Our data and modeling hints at an extended cold outskirt of the ring. Comparison with the morphology of the CO gas in the disk reveals an inner dust-free hole where gas nevertheless can be detected. Based on dust grain lifetimes, we propose that the dust content of this gaseous disk is of secondary origin produced by planetesimals. Since the gas component is probably primordial, HD 21997 is one of the first known examples of a hybrid circumstellar disk, a so-far little studied late phase of circumstellar disk evolution.

preprint2013arXiv

ALMA observations of the molecular gas in the debris disk of the 30 Myr old star HD 21997

The 30 Myr old A3-type star HD 21997 is one of the two known debris dust disks having a measurable amount of cold molecular gas. With the goal of understanding the physical state, origin, and evolution of the gas in young debris disks, we obtained CO line observations with the Atacama Large Millimeter/submillimeter Array (ALMA). Here we report on the detection of 12CO and 13CO in the J=2-1 and J=3-2 transitions and C18O in the J=2-1 line. The gas exhibits a Keplerian velocity curve, one of the few direct measurements of Keplerian rotation in young debris disks. The measured CO brightness distribution could be reproduced by a simple star+disk system, whose parameters are r_in < 26 AU, r_out = 138 +/- 20 AU, M_*=1.8 +0.5 -0.2 M_Sun, and i = 32.6 +/- 3.1 degrees. The total CO mass, as calculated from the optically thin C18O line, is about (4-8) x 10^-2 M_Earth, while the CO line ratios suggest a radiation temperature on the order of 6-9 K. Comparing our results with those obtained for the dust component of the HD 21997 disk from the ALMA continuum observations by Moór et al., we conclude that comparable amounts of CO gas and dust are present in the disk. Interestingly, the gas and dust in the HD 21997 system are not co-located, indicating a dust-free inner gas disk within 55 AU of the star. We explore two possible scenarios for the origin of the gas. A secondary origin, which involves gas production from colliding or active planetesimals, would require unreasonably high gas production rates and would not explain why the gas and dust are not co-located. We propose that HD 21997 is a hybrid system where secondary debris dust and primordial gas coexist. HD 21997, whose age exceeds both the model predictions for disk clearing and the ages of the oldest T Tauri-like or transitional gas disks in the literature, may be a key object linking the primordial and the debris phases of disk evolution.

preprint2013arXiv

ATLASGAL --- Environments of 6.7\,GHz methanol masers

Using the 870-$μ$m APEX Telescope Large Area Survey of the Galaxy (ATLASGAL), we have identified 577 submillimetre continuum sources with masers from the methanol multibeam (MMB) survey in the region $280\degr < \ell < 20\degr$; $|\,b\,| < 1.5\degr$. 94\,per\,cent of methanol masers in the region are associated with sub-millimetre dust emission. We estimate masses for ~450 maser-associated sources and find that methanol masers are preferentially associated with massive clumps. These clumps are centrally condensed, with envelope structures that appear to be scale-free, the mean maser position being offset from the peak column density by 0.0 \pm 4". Assuming a Kroupa initial mass function and a star-formation efficiency of ~30\,per\,cent, we find that over two thirds of the clumps are likely to form clusters with masses >20\,\msun. Furthermore, almost all clumps satisfy the empirical mass-size criterion for massive star formation. Bolometric luminosities taken from the literature for ~100 clumps range between ~100 and 10$^6$\,\lsun. This confirms the link between methanol masers and massive young stars for 90\,per\,cent of our sample. The Galactic distribution of sources suggests that the star-formation efficiency is significantly reduced in the Galactic-centre region, compared to the rest of the survey area, where it is broadly constant, and shows a significant drop in the massive star-formation rate density in the outer Galaxy. We find no enhancement in source counts towards the southern Scutum-Centaurus arm tangent at $\ell ~ 315\degr$, which suggests that this arm is not actively forming stars.

preprint2013arXiv

ATLASGAL --- properties of compact HII regions and their natal clumps

We present a complete sample of molecular clumps containing compact and ultra-compact (UC) HII regions between \ell=10\degr and 60\degr\ and $|b|<1\degr, identified by combining the the ATLASGAL submm and CORNISH radio continuum surveys with visual examination of archival infrared data. Our sample is complete to optically thin, compact and UCHII regions driven by a zero age main sequence star of spectral type B0 or earlier embedded within a 1,000 Msun clump. In total we identify 213 compact and UCHII regions, associated with 170 clumps. Unambiguous kinematic distances are derived for these clumps and used to estimate their masses and physical sizes, as well as the Lyman continuum fluxes and sizes of their embedded HII regions. We find a clear lower envelope for the surface density of molecular clumps hosting massive star formation of 0.05 g cm^{-2}, which is consistent with a similar sample of clumps associated with 6.7 GHz masers. The mass of the most massive embedded stars is closely correlated with the mass of their natal clump. Young B stars appear to be significantly more luminous in the ultraviolet than predicted by current stellar atmosphere models. The properties of clumps associated with compact and UCHII regions are very similar to those associated with 6.7 GHz methanol masers and we speculate that there is little evolution in the structure of the molecular clumps between these two phases. Finally, we identify a significant peak in the surface density of compact and UCHII regions associated with the W49A star-forming complex, noting that this complex is truly one of the most massive and intense regions of star formation in the Galaxy.

preprint2013arXiv

MALT90: The Millimetre Astronomy Legacy Team 90 GHz Survey

The Millimetre Astronomy Legacy Team 90 GHz (MALT90) survey aims to characterise the physical and chemical evolution of high-mass star-forming clumps. Exploiting the unique broad frequency range and on-the-fly mapping capabilities of the Australia Telescope National Facility Mopra 22 m single-dish telescope, MALT90 has obtained 3' x 3' maps toward ~2000 dense molecular clumps identified in the ATLASGAL 870 um Galactic plane survey. The clumps were selected to host the early stages of high-mass star formation and to span the complete range in their evolutionary states (from prestellar, to protostellar, and on to HII regions and photodissociation regions). Because MALT90 mapped 16 lines simultaneously with excellent spatial (38") and spectral (0.11 km/s) resolution, the data reveal a wealth of information about the clump's morphologies, chemistry, and kinematics. In this paper we outline the survey strategy, observing mode, data reduction procedure, and highlight some early science results. All MALT90 raw and processed data products are available to the community. With its unprecedented large sample of clumps, MALT90 is the largest survey of its type ever conducted and an excellent resource for identifying interesting candidates for high resolution studies with ALMA.

preprint2013arXiv

Pillars and globules at the edges of H ii regions, Confronting Herschel observations and numerical simulations

Pillars and globules are present in many high-mass star-forming regions, such as the Eagle nebula (M16) and the Rosette molecular cloud, and understanding their origin will help characterize triggered star formation. The formation mechanisms of these structures are still being debated. Recent numerical simulations have shown how pillars can arise from the collapse of the shell in on itself and how globules can be formed from the interplay of the turbulent molecular cloud and the ionization from massive stars. The goal here is to test this scenario through recent observations of two massive star-forming regions, M16 and Rosette. The column density structure of the interface between molecular clouds and H ii regions was characterized using column density maps obtained from far-infrared imaging of the Herschel HOBYS key programme. Then, the DisPerSe algorithm was used on these maps to detect the compressed layers around the ionized gas and pillars in different evolutionary states. Finally, their velocity structure was investigated using CO data, and all observational signatures were tested against some distinct diagnostics established from simulations. The column density profiles have revealed the importance of compression at the edge of the ionized gas. The velocity properties of the structures, i.e. pillars and globules, are very close to what we predict from the numerical simulations. We have identified a good candidate of a nascent pillar in the Rosette molecular cloud that presents the velocity pattern of the shell collapsing on itself, induced by a high local curvature. Globules have a bulk velocity dispersion that indicates the importance of the initial turbulence in their formation, as proposed from numerical simulations. Altogether, this study re-enforces the picture of pillar formation by shell collapse and globule formation by the ionization of highly turbulent clouds.

preprint2012arXiv

ATLASGAL - Compact source catalogue: 330 < l < 21 degrees

Context. The APEX Telescope Large Area Survey of the Galaxy (ATLASGAL) is the first systematic survey of the inner Galactic plane in the sub-millimetre. The observations were carried out with the Large APEX Bolometer Camera (LABOCA), an array of 295 bolometers observing at 870 microns (345 GHz). Aims. Here we present a first version of the compact source catalogue extracted from this survey. This catalogue provides an unbiased database of dusty clumps in the inner Galaxy. Methods. The construction of this catalogue was made using the source extraction routine SExtractor. We have cross-associated the obtained sources with the IRAS and MSX catalogues, in order to constrain their nature. Results. We have detected 6639 compact sources in the range from 330 < l < 21 degrees and |b| < 1.5 degrees. The catalogue has a 99% completeness for sources with a peak flux above 6 sigma, which corresponds to a flux density of ~0.4 Jy/beam. The parameters extracted for sources with peak fluxes below the 6 sigma completeness threshold should be used with caution. Tests on simulated data find the uncertainty in the flux measurement to be ~12%, however, in more complex regions the flux values can be overestimated by a factor of 2 due to the additional background emission. Using a search radius of 30" we found that 40% of ATLASGAL compact sources are associated with an IRAS or MSX point source, but, ~50% are found to be associated with MSX 21 microns fluxes above the local background level, which is probably a lower limit to the actual number of sources associated with star formation. Conclusions. Although infrared emission is found towards the majority of the clumps detected, this catalogue is still likely to include a significant number of clumps that are devoid of star formation activity and therefore excellent candidates for objects in the coldest, earliest stages of (high-mass) star formation.

preprint2012arXiv

Cluster-formation in the Rosette molecular cloud at the junctions of filaments

For many years feedback processes generated by OB-stars in molecular clouds, including expanding ionization fronts, stellar winds, or UV-radiation, have been proposed to trigger subsequent star formation. However, hydrodynamic models including radiation and gravity show that UV-illumination has little or no impact on the global dynamical evolution of the cloud. The Rosette molecular cloud, irradiated by the NGC2244 cluster, is a template region for triggered star-formation, and we investigated its spatial and density structure by applying a curvelet analysis, a filament-tracing algorithm (DisPerSE), and probability density functions (PDFs) on Herschel column density maps, obtained within the HOBYS key program. The analysis reveals not only the filamentary structure of the cloud but also that all known infrared clusters except one lie at junctions of filaments, as predicted by turbulence simulations. The PDFs of sub-regions in the cloud show systematic differences. The two UV-exposed regions have a double-peaked PDF we interprete as caused by shock compression. The deviations of the PDF from the log-normal shape typically associated with low- and high-mass star-forming regions at Av~3-4m and 8-10m, respectively, are found here within the very same cloud. This shows that there is no fundamental difference in the density structure of low- and high-mass star-forming regions. We conclude that star-formation in Rosette - and probably in high-mass star-forming clouds in general - is not globally triggered by the impact of UV-radiation. Moreover, star formation takes place in filaments that arose from the primordial turbulent structure built up during the formation of the cloud. Clusters form at filament mergers, but star formation can be locally induced in the direct interaction zone between an expanding HII--region and the molecular cloud.

preprint2012arXiv

Globules and pillars seen in the [CII] 158 micron line with SOFIA

Molecular globules and pillars are spectacular features, found only in the interface region between a molecular cloud and an HII-region. Impacting Far-ultraviolet (FUV) radiation creates photon dominated regions (PDRs) on their surfaces that can be traced by typical cooling lines. With the GREAT receiver onboard SOFIA we mapped and spectrally resolved the [CII] 158 micron atomic fine-structure line and the highly excited 12CO J=11-10 molecular line from three objects in Cygnus X (a pillar, a globule, and a strong IRAS source). We focus here on the globule and compare our data with existing Spitzer data and recent Herschel Open-Time PACS data. Extended [CII] emission and more compact CO-emission was found in the globule. We ascribe this emission mainly to an internal PDR, created by a possibly embedded star-cluster with at least one early B-star. However, external PDR emission caused by the excitation by the Cyg OB2 association cannot be fully excluded. The velocity-resolved [CII] emission traces the emission of PDR surfaces, possible rotation of the globule, and high-velocity outflowing gas. The globule shows a velocity shift of ~2 km/s with respect to the expanding HII-region, which can be understood as the residual turbulence of the molecular cloud from which the globule arose. This scenario is compatible with recent numerical simulations that emphazise the effect of turbulence. It is remarkable that an isolated globule shows these strong dynamical features traced by the [CII]-line, but it demands more observational studies to verify if there is indeed an embedded cluster of B-stars.

preprint2012arXiv

SOFIA observations of far-infrared hydroxyl emission toward classical ultracompact HII/OH maser regions

The hydroxyl radical (OH) is found in various environments within the interstellar medium (ISM) of the Milky Way and external galaxies, mostly either in diffuse interstellar clouds or in the warm, dense environments of newly formed low-mass and high-mass stars, i.e, in the dense shells of compact and ultracompact HII regions (UCHIIRs). Until today, most studies of interstellar OH involved the molecule's radio wavelength hyperfine structure (hfs) transitions. These lines are generally not in LTE and either masing or over-cooling complicates their interpretation. In the past, observations of transitions between different rotational levels of OH, which are at far-infrared wavelengths, have suffered from limited spectral and angular resolution. Since these lines have critical densities many orders of magnitude higher than the radio wavelength ground state hfs lines and are emitted from levels with more than 100 K above the ground state, when observed in emission, they probe very dense and warm material. We probe the warm and dense molecular material surrounding the UCHIIR/OH maser sources W3(OH), G10.62-0.39 and NGC 7538 IRS1 by studying the $^2Π_{1/2}, J = {3/2} - {1/2}$ rotational transition of OH in emission and, toward the last source also the molecule's $^2Π_{3/2}, J = 5/2 - 3/2$ ground-state transition in absorption. We used the Stratospheric Observatory for Infrared Astronomy (SOFIA) to observe these OH lines, which are near 1.84 THz ($163 μ$m) and 2.51 THz ($119.3 μ$m). We clearly detect the OH lines, some of which are blended with each other. Employing non-LTE radiative transfer calculations we predict line intensities using models of a low OH abundance envelope versus a compact, high-abundance source corresponding to the origin of the radio OH lines.

preprint2012arXiv

The chemistry of C3 & Carbon Chain Molecules in DR21(OH)

(Abridged) We have observed velocity resolved spectra of four ro-vibrational far-infrared transitions of C3 between the vibrational ground state and the low-energy nu2 bending mode at frequencies between 1654--1897 GHz using HIFI on board Herschel, in DR21(OH), a high mass star forming region. Several transitions of CCH and c-C3H2 have also been observed with HIFI and the IRAM 30m telescope. A gas and grain warm-up model was used to identify the primary C3 forming reactions in DR21(OH). We have detected C3 in absorption in four far-infrared transitions, P(4), P(10), Q(2) and Q(4). The continuum sources MM1 and MM2 in DR21(OH) though spatially unresolved, are sufficiently separated in velocity to be identified in the C3 spectra. All C3 transitions are detected from the embedded source MM2 and the surrounding envelope, whereas only Q(4) & P(4) are detected toward the hot core MM1. The abundance of C3 in the envelope and MM2 is \sim6x10^{-10} and \sim3x10^{-9} respectively. For CCH and c-C3H2 we only detect emission from the envelope and MM1. The observed CCH, C3, and c-C3H2 abundances are most consistent with a chemical model with n(H2)\sim5x10^{6} cm^-3 post-warm-up dust temperature, T_max =30 K and a time of \sim0.7-3 Myr. Post warm-up gas phase chemistry of CH4 released from the grain at t\sim 0.2 Myr and lasting for 1 Myr can explain the observed C3 abundance in the envelope of DR21(OH) and no mechanism involving photodestruction of PAH molecules is required. The chemistry in the envelope is similar to the warm carbon chain chemistry (WCCC) found in lukewarm corinos. The observed lower C3 abundance in MM1 as compared to MM2 and the envelope could be indicative of destruction of C3 in the more evolved MM1. The timescale for the chemistry derived for the envelope is consistent with the dynamical timescale of 2 Myr derived for DR21(OH) in other studies.

preprint2012arXiv

The spine of the swan: A Herschel study of the DR21 ridge and filaments in Cygnus X

In order to characterise the cloud structures responsible for the formation of high-mass stars, we present Herschel observations of the DR21 environment. Maps of the column density and dust temperature unveil the structure of the DR21 ridge and several connected filaments. The ridge has column densities larger than 1e23/cm^2 over a region of 2.3 pc^2. It shows substructured column density profiles and branching into two major filaments in the north. The masses in the studied filaments range between 130 and 1400 Msun whereas the mass in the ridge is 15000 Msun. The accretion of these filaments onto the DR21 ridge, suggested by a previous molecular line study, could provide a continuous mass inflow to the ridge. In contrast to the striations seen in e.g., the Taurus region, these filaments are gravitationally unstable and form cores and protostars. These cores formed in the filaments potentially fall into the ridge. Both inflow and collisions of cores could be important to drive the observed high-mass star formation. The evolutionary gradient of star formation running from DR21 in the south to the northern branching is traced by decreasing dust temperature. This evolution and the ridge structure can be explained by two main filamentary components of the ridge that merged first in the south.

preprint2011arXiv

Convergent flows and low-velocity shocks in DR21(OH)

DR21(OH) is a pc-scale massive, 7000 Msun clump hosting three massive dense cores (MDCs) at an early stage of their evolution. We present a high angular-resolution mosaic, covering 70" by 100", with the IRAM PdBI at 3 mm to trace the dust continuum emission and the N2H+ (J=1-0) and CH3CN (J=5-4) molecular emission. The cold, dense gas traced by the compact emission in N2H+ is associated with the three MDCs and shows several velocity components towards each MDC. These velocity components reveal local shears in the velocity fields which are best interpreted as convergent flows. Moreover, we report the detection of weak extended emission from CH3CN at the position of the N2H+ velocity shears. We propose that this extended CH3CN emission is tracing warm gas associated with the low-velocity shocks expected at the location of convergence of the flows where velocity shears are observed. This is the first detection of low-velocity shocks associated with small (sub-parsec) scale convergent flows which are proposed to be at the origin of the densest structures and of the formation of (high-mass) stars. In addition, we propose that MDCs may be active sites of star-formation for more than a crossing time as they continuously receive material from larger scale flows as suggested by the global picture of dynamical, gravity driven evolution of massive clumps which is favored by the present observations.

preprint2011arXiv

Molecular gas in young debris disks

Gas-rich primordial disks and tenuous gas-poor debris disks are usually considered as two distinct evolutionary phases of the circumstellar matter. Interestingly, the debris disk around the young main-sequence star 49 Ceti possesses a substantial amount of molecular gas, and possibly represents the missing link between the two phases. Motivated to understand the evolution of the gas component in circumstellar disks via finding more 49 Ceti-like systems, we carried out a CO J=3-2 survey with Atacama Pathfinder EXperiment, targeting 20 infrared-luminous debris disks. These systems fill the gap between primordial and old tenuous debris disks in terms of fractional luminosity. Here we report on the discovery of a second 49 Ceti-like disk around the 30 Myr old A3-type star HD21997, a member of the Columba Association. This system was also detected in the CO(2-1) transition, and the reliable age determination makes it an even clearer example of an old gas-bearing disk than 49 Ceti. While the fractional luminosities of HD21997 and 49 Ceti are not particularly high, these objects seem to harbor the most extended disks within our sample. The double-peaked profiles of HD21997 were reproduced by a Keplerian disk model combined with the LIME radiative transfer code. Based on their similarities, 49 Ceti and HD21997 may be the first representatives of a so far undefined new class of relatively old (>~8 Myr), gaseous dust disks. From our results, neither primordial origin nor steady secondary production from icy planetesimals can unequivocally explain the presence of CO gas in the disk of HD21997.

preprint2010arXiv

Dynamic star formation in the massive DR21 filament

The formation of massive stars is a highly complex process in which it is not clear whether the star-forming gas is in global gravitational collapse or in an equilibrium state, supported by turbulence. By studying one of the most massive and dense star-forming regions in the Galaxy at a distance of less than 3 kpc, the filament containing the well-known sources DR21 and DR21(OH), we expect to find observational signatures that allow to discriminate between the two views. We use molecular line data from our 13CO 1-0, CS 2-1, and N2H+ 1-0 survey of the Cygnus X region obtained with the FCRAO and high-angular resolution observations of CO, CS, HCO+, N2H+, and H2CO, obtained with the IRAM 30m telescope. We observe a complex velocity field and velocity dispersion in the DR21 filament in which regions of highest column-density, i.e. dense cores, have a lower velocity dispersion than the surrounding gas and velocity gradients that are not (only) due to rotation. Infall signatures in optically thick line profiles of HCO+ and 12CO are observed along and across the whole DR21 filament. From modelling the observed spectra, we obtain a typical infall speed of 0.6 km/s and mass accretion rates of the order of a few 10^-3 Msun/yr for the two main clumps constituting the filament. These massive (4900 and 3300 Msun) clumps are both gravitationally contracting. All observed kinematic features in the DR21 filament can be explained if it is formed by the convergence of flows at large scales and is now in a state of global gravitational collapse. Whether this convergence of flows originated from self-gravity at larger scales or from other processes can not be settled with the present study. The observed velocity field and velocity dispersion are consistent with results from (magneto)-hydrodynamic simulations where the cores lie at the stagnation points of convergent turbulent flows.

preprint2010arXiv

Gas dynamics in Massive Dense Cores in Cygnus-X

We study the kinematic properties of dense gas surrounding massive protostars recognized by Bontemps et a. (2010) in a sample of five Massive Dense Cores in Cygnus-X. We investigate whether turbulent support plays a major role in stabilizing the core against fragmentation into Jeans-mass objects or alternatively, the observed kinematics could indicate a high level of dynamics. We present IRAM 30m single-dish (HCO+ and H13CO+) and IRAM PdBI high angular-resolution observations of dense gas tracers (H13CO+ and H13CN) to reveal the kinematics of molecular gas at scales from 0.03 to 0.1 pc. Radiative transfer modeling shows that H13CO+ is depleted within the envelopes of massive protostars and traces the bulk of material surrounding the protostars rather than their inner envelopes. H13CN shows a better correspondence with the peak of the continuum emission, possibly due to abundance anomalies and specific chemistry in the close vicinity of massive protostars. Analyzing the line-widths we show that the observed line-dispersion of H13CO+ at the scale of MDCs is smaller than expected from the quasi-static, turbulent-core model. At large-scales, global organized bulk motions are identified for 3 of the MDCs. At small-scales, several spectral components are identified in all MDCs showing filamentary structures and intrinsic velocity gradients towards the continuum peaks. The dynamics of these flows show diversity among the sample and we link this to the specific fragmentation properties of the MDCs. No clear evidence is found for a turbulence regulated, equilibrium scenario within the sample of MDCs. We propose a picture in which MDCs are not in equilibrium and their dynamics is governed by small-scale converging flows, which may initiate star-formation via their shears.

preprint2010arXiv

Herschel observations of embedded protostellar clusters in the Rosette Molecular Cloud

The Herschel OB young stellar objects survey (HOBYS) has observed the Rosette molecular cloud, providing an unprecedented view of its star formation activity. These new far-infrared data reveal a population of compact young stellar objects whose physical properties we aim to characterise. We compiled a sample of protostars and their spectral energy distributions that covers the near-infrared to submillimetre wavelength range. These were used to constrain key properties in the protostellar evolution, bolometric luminosity, and envelope mass and to build an evolutionary diagram. Several clusters are distinguished including the cloud centre, the embedded clusters in the vicinity of luminous infrared sources, and the interaction region. The analysed protostellar population in Rosette ranges from 0.1 to about 15 Msun with luminosities between 1 and 150 Lsun, which extends the evolutionary diagram from low-mass protostars into the high-mass regime. Some sources lack counterparts at near- to mid-infrared wavelengths, indicating extreme youth. The central cluster and the Phelps & Lada 7 cluster appear less evolved than the remainder of the analysed protostellar population. For the central cluster, we find indications that about 25% of the protostars classified as Class I from near- to mid-infrared data are actually candidate Class 0 objects. As a showcase for protostellar evolution, we analysed four protostars of low- to intermediate-mass in a single dense core, and they represent different evolutionary stages from Class 0 to Class I. Their mid- to far-infrared spectral slopes flatten towards the Class I stage, and the 160 to 70um flux ratio is greatest for the presumed Class 0 source. This shows that the Herschel observations characterise the earliest stages of protostellar evolution in detail.

preprint2010arXiv

Small-scale structure in the Rosette molecular cloud revealed by Herschel

We present a preliminary analysis of the small-scale structure found in new 70-520 micron continuum maps of the Rosette molecular cloud (RMC), obtained with the SPIRE and PACS instruments of the Herschel Space Observatory. We find 473 clumps within the RMC using a new structure identification algorithm, with sizes up to ~1.0 pc in diameter. A comparison with recent Spitzer maps reveals that 371 clumps are "starless" (without an associated young stellar object), while 102 are "protostellar." Using the respective values of dust temperature, we determine the clumps have masses (M_C) over the range -0.75 <= log (M_C/M_sun) <= 2.50. Linear fits to the high-mass tails of the resulting clump mass spectra (CMS) have slopes that are consistent with those found for high-mass clumps identified in CO emission by other groups.

preprint2010arXiv

Structure and evolution of debris disks around F-type stars: I. Observations, database and basic evolutionary aspects

Although photometric and spectroscopic surveys with the Spitzer Space Telescope increased remarkably the number of well studied debris disks around A-type and Sun-like stars, detailed analyzes of debris disks around F-type stars remained less frequent. Using the MIPS camera and the IRS spectrograph we searched for debris dust around 82 F-type stars with Spitzer. We found 27 stars that harbor debris disks, nine of which are new discoveries. The dust distribution around two of our stars, HD 50571 and HD 170773, was found to be marginally extended on the 70um MIPS images. Combining the MIPS and IRS measurements with additional infrared and submillimeter data, we achieved excellent spectral coverage for most of our debris systems. We have modeled the excess emission of 22 debris disks using a single temperature dust ring model and of 5 debris systems with two-temperature models. The latter systems may contain two dust rings around the star. In accordance with the expected trends, the fractional luminosity of the disks declines with time, exhibiting a decay rate consistent with the range of model predictions. We found the distribution of radial dust distances as a function of age to be consistent with the predictions of both the self stirred and the planetary stirred disk evolution models. A more comprehensive investigation of the evolution of debris disks around F-type stars, partly based on the presented data set, will be the subject of an upcoming paper.

preprint2010arXiv

The Herschel view of star formation in the Rosette molecular cloud under the influence of NGC 2244

The Rosette molecular cloud is promoted as the archetype of a triggered star-formation site. This is mainly due to its morphology, because the central OB cluster NGC 2244 has blown a circular-shaped cavity into the cloud and the expanding HII-region now interacts with the cloud. Studying the spatial distribution of the different evolutionary states of all star-forming sites in Rosette and investigating possible gradients of the dust temperature will help to test the 'triggered star-formation' scenario in Rosette. We use continuum data obtained with the PACS (70 and 160 micron) and SPIRE instruments (250, 350, 500 micron) of the Herschel telescope during the Science Demonstration Phase of HOBYS. Three-color images of Rosette impressively show how the molecular gas is heated by the radiative impact of the NGC 2244 cluster. A clear negative temperature gradient and a positive density gradient (running from the HII-region/molecular cloud interface into the cloud) are detected. Studying the spatial distribution of the most massive dense cores (size scale 0.05 to 0.3 pc), we find an age-sequence (from more evolved to younger) with increasing distance to the cluster NGC 2244. No clear gradient is found for the clump (size-scale up to 1 pc) distribution. The existence of temperature and density gradients and the observed age-sequence imply that star formation in Rosette may indeed be influenced by the radiative impact of the central NGC 2244 cluster. A more complete overview of the prestellar and protostellar population in Rosette is required to obtain a firmer result.

preprint2009arXiv

Fragmentation and mass segregation in the massive dense cores of Cygnus X

We present Plateau de Bure interferometer observations obtained in continuum at 1.3 and 3.5 mm towards the six most massive and young (IR-quiet) dense cores in Cygnus X. Located at only 1.7 kpc, the Cygnus X region offers the opportunity of reaching small enough scales (of the order of 1700 AU at 1.3 mm) to separate individual collapsing objects. The cores are sub-fragmented with a total of 23 fragments inside 5 cores. Only the most compact core, CygX-N63, could actually be a single massive protostar with an envelope mass as large as 60 Msun. The fragments in the other cores have sizes and separations similar to low-mass pre-stellar and proto-stellar condensations in nearby protoclusters, and are probably of the same nature. A total of 9 out of these 23 protostellar objects are found to be probable precursors of OB stars with envelope masses ranging from 6 to 23 Msun. The level of fragmentation is globally higher than in the turbulence regulated, monolithic collapse scenario, but is not as high as expected in a pure gravo-turbulent scenario where the distribution of mass is dominated by low-mass protostars/stars. Here, the fractions of the total core masses in the high-mass fragments are reaching values as high as 28, 44, and 100 % in CygX-N12, CygX-N53, and CygX-N63, respectively, much higher than what an IMF-like mass distribution would predict. The increase of the fragmentation efficiency as a function of density in the cores is proposed to be due to the increasing importance of self-gravity leading to gravitational collapse at the scale of the dense cores. At the same time, the cores tend to fragment into a few massive protostars within their central regions. We are therefore probably witnessing here the primordial mass segregation of clusters in formation.