Source author record

H. Beuther

H. Beuther 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

96works
5topics
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

96 published item(s)

preprint2026arXiv

The converging gas flow around the infrared dark cloud G28.3

Aims: The G28.37+0.07 star-forming region is a prototypical infrared dark cloud (IRDC) located at the interface of a converging gas flow. This study characterizes the properties of this dynamic gas flow. Methods: Combining data from the Northern Extended Millimeter Array (NOEMA) with single-dish data from the IRAM30m observatory, we mapped large spatial scales (~81pc^2) at high angular resolution (7.0''x2.6'' corresponding ~2.3x10^4au or ~0.1pc) down to core scales. The spectral setup in the 3mm band covers many spectral lines as well as the continuum emission. Results: The data reveal the proposed west-east converging gas flow in all observed dense gas tracers. We estimate a mass-flow rate along that flow around 10^-3M_sun/yr. Comparing these west-east flow rates to infall rates toward sources along the line of sight, the gas flow rates are roughly a factor of 25 greater than than those along the line of sight. This confirms the dominance of longitudinal motions along the converging gas flow in G28.37. For comparison, in the main north-south IRDC formed by the west-east converging gas flow, infall rates along the line of sight are about an order of magnitude greater than those along the west-east flow. In addition to the kinematic analysis, a comparison of CH_3CN-derived gas temperatures with Herschel-derived dust temperatures typically show higher gas temperatures toward high-density sources. We discuss whether mechanical heating from the conversion of the flow's kinetic energy into thermal energy may explain some of the observed temperature differences. Conclusions: The differences between flow rates along the converging flow, perpendicular to it, and toward the sources at the IRDC center indicate that at the interfaces of converging gas flows - where most of the active star formation takes place - originally more directed gas flows can convert into multidirectional infall motions.

preprint2025arXiv

ALMAGAL V. Relations between the core populations and the parent clump physical properties

Context. The fragmentation of massive molecular clumps into smaller, potentially star-forming cores plays a key role in the processes of high-mass star formation. The ALMAGAL project offers high-resolution data to investigate these processes across various evolutionary stages in the Galactic plane. Aims. This study aims at correlating the fragmentation properties of massive clumps, obtained from ALMA observations, with their global physical parameters (e.g., mass, surface density, and temperature) and evolutionary indicators (such as luminosity-to-mass ratio and bolometric temperature) obtained from Herschel observations. It seeks to assess whether the cores evolve in number and mass in tandem with their host clumps, and to determine the possible factors influencing the formation of massive cores (M > 24M_\odot). Methods. We analyzed the masses of 6348 fragments, estimated from 1.4 mm continuum data for 1007 ALMAGAL clumps. Leveraging this unprecedentedly large data set, we evaluated statistical relationships between clump parameters, estimated over about 0.1 pc scales, and fragment properties, corresponding to scales of a few 1000 au, while accounting for potential biases related to distance and observational resolution. Our results were further compared with predictions from numerical simulations. Results. The fragmentation level correlates preferentially with clump surface density, supporting a scenario of density-driven fragmentation, whereas it does not show any clear dependence on total clump mass. Both the mass of the most massive core and the core formation efficiency show a broad range and increase on average by an order of magnitude in the intervals spanned by evolutionary indicators such as clump dust temperature and the luminosity-to-mass ratio. This suggests that core growth continues throughout the clump evolution, favoring clump-fed over core-fed theoretical scenarios.

preprint2025arXiv

From theory to observation: understanding filamentary flows in high-mass star-forming clusters

Here we use data from multi-scale galactic MHD simulations to observe filaments and star forming clumps on 10's of pc scales and investigate flow rate relationships along, and onto filaments as well as flows towards the clumps. Using the FilFinderPPV identification technique, we identify the prominent filamentary structures in each data cube. Each filament and its corresponding clump are analysed by calculating flow rates along each filament towards the clump, onto each filament from increasing distances, and radially around each clump. This analysis is conducted for two cubes, one feedback dominated region, and one with less feedback. Looking at the face-on inclination of the simulations (0 degrees), we observe different trends depending on the environmental conditions (more or less feedback). The median flow rate in the region with more feedback is 8.9$\times$10$^{-5}$ M$_{sun}\mathrm{yr}^{-1}$ and we see that flow rates along the filaments toward the clumps generally decrease in these regions. In the region with less feedback we have a median flow rate of 2.9$\times$10$^{-4}$ M$_{sun}\mathrm{yr}^{-1}$ and when looking along the filaments here we see the values either increase or remain constant. We find that the flow rates from the environments onto the primary filaments are of an order of magnitude sufficient to sustain the flow rates along these filaments. When discussing the effects of galactic and filamentary inclination, we also observe that viewing the filaments from different galactic inclinations can reveal the presence of feeder structures (smaller filamentary structures aiding in the flow of material). The method used to estimate these flow rates, which has been previously applied to observational data, produced results consistent with those obtained from the simulations themselves, providing high confidence in the flow rate calculation method.

preprint2022arXiv

FEEDBACK from the NGC7538 HII region

Context: How do expanding HII regions interact with their environmental cloud? This is one of the central questions driving the SOFIA legacy program FEEDBACK. Here, we present a case study toward the prototypical H{\sc ii} region NGC7538. Methods: With SOFIA we mapped an area of ~210'^2 around NGC7538 in the [CII] line at 1.9THz. Complementary observed atomic carbon [CI] and high-J CO(8-7) data as well as archival NIR/FIR, cm continuum, CO(3-2) and HI data are folded into the analysis. Results: While the overall [CII] morphology follows the general ionized gas, the channel maps show multiple bubble-like structures with sizes on the order of ~80-100" (~1.0-1.28pc). While at least one of them may be an individual feedback bubble driven by the main exciting sources of the region, the other bubble-morphologies may also be due to the intrinsically porous structure of the HII region. An analysis of the expansion velocities around 10km s^{-1} indicates that thermal expansion is not sufficient but that wind-driving from the central O-stars is required. The most blue-shifted [CII] component has barely any molecular or atomic counterparts. At the interface to the molecular cloud, we find a typical photon-dominated region (PDR) with a bar-shape. Ionized, atomic and molecular carbon show a layered structure in this PDR. The carbon in the PDR is dominated by its ionized form with atomic and molecular masses of ~0.45+-0.1M_{\odot} and ~1.2+-0.1M_{\odot}, respectively, compared to the ionized carbon in the range of 3.6-9.7M_{\odot}. Conclusions: The NGC7538 HII region exhibits a diverse set of sub-structures that interact with each other as well as with the adjacent cloud. Compared to other recent [CII] observations of HII regions (e.g., Orion Veil, RCW120, RCW49), bubble-shape morphologies revealed in [CII] emission, indicative of expanding shells, are recurring structures of PDRs.

preprint2022arXiv

On the accuracy of HI observations in molecular clouds -- More cold HI than thought?

We present a study of the cold atomic hydrogen (HI) content of molecular clouds simulated within the SILCC-Zoom project. We produce synthetic observations of HI at 21 cm including HI self-absorption (HISA) and observational effects. We find that HI column densities, $N_\textrm{HI}$, of $\gtrsim$10$^{22}$ cm$^{-2}$ are frequently reached in molecular clouds with HI temperatures as low as $\sim$10 K. Hence, HISA observations assuming a fixed HI temperature tend to underestimate the amount of cold HI in molecular clouds by a factor of 3 - 10 and produce an artificial upper limit of $N_\textrm{HI}$ around 10$^{21}$ cm$^{-2}$. We thus argue that the cold HI mass in clouds could be a factor of a few higher than previously estimated. Also $N_\textrm{HI}$-PDFs obtained from HISA observations might be subject to observational biases and should be considered with caution. The underestimation of cold HI in HISA observations is due to both the large HI temperature variations and the effect of noise in regions of high optical depth. We find optical depths of cold HI around 1 - 10 making optical depth corrections essential. We show that the high HI column densities ($\gtrsim$10$^{22}$ cm$^{-2}$) can in parts be attributed to the occurrence of up to 10 individual HI-H$_2$ transitions along the line of sight. This is also reflected in the spectra, necessitating Gaussian decomposition algorithms for their analysis. However, also for a single HI-H$_2$ transition, $N_\textrm{HI}$ frequently exceeds 10$^{21}$ cm$^{-2}$, challenging 1D, semi-analytical models. This is due to non-equilibrium chemistry effects and the fact that HI-H$_2$ transition regions usually do not possess a 1-dimensional geometry. Finally, we show that the HI gas is moderately supersonic with Mach numbers of a few. The corresponding non-thermal velocity dispersion can be determined via HISA observations within a factor of $\sim$2.

preprint2022arXiv

The Cygnus Allscale Survey of Chemistry and Dynamical Environments: CASCADE: Overview and first results toward DR20 from the Max Planck IRAM Observatory program (MIOP)

Context: While star formation on large molecular cloud scales and on small core and disk scales has been investigated intensely over the past decades, the connection of the large-scale interstellar material with the densest small-scale cores has been a largely neglected field. Methods: Using NOEMA and the IRAM 30\,m telescope, we mapped large areas (640\,arcmin$^2$) of the archetypical star formation complex Cygnus X at 3.6\,mm wavelengths in line and continuum emission. Results: The scope and outline of The Cygnus Allscale Survey of Chemistry and Dynamical Environments (CASCADE) is presented. We then focus on the first observed subregion in Cygnus X, namely the DR20 star formation site, which comprises sources in a range of evolutionary stages from cold pristine gas clumps to more evolved ultracompact H{\sc ii} regions. The data covering cloud to cores scales at a linear spatial resolution of $<5000$\,au reveal several kinematic cloud components that are likely part of several large-scale flows onto the central cores. The temperature structure of the region is investigated by means of the HCN/HNC intensity ratio and compared to dust-derived temperatures. We find that the deuterated DCO$^+$ emission is almost exclusively located toward regions at low temperatures below 20\,K. Investigating the slopes of spatial power spectra of dense gas tracer intensity distributions (HCO$^+$, H$^{13}$CO$^+$, and N$_2$H$^+$), we find comparatively flat slopes between $-2.9$ and $-2.6$, consistent with high Mach numbers and/or active star formation in DR20.

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 sharp ALMA view of infall and outflow in the massive protocluster G31.41+0.31

Context. To better understand the formation of high-mass stars, it is fundamental to investigate how matter accretes onto young massive stars, how it is ejected, and how all this differs from the low-mass case. The massive protocluster G31.41+0.31 is the ideal target to study all these processes because observations at millimeter and centimeter wavelengths have resolved the emission of the Main core into at least four massive dust continuum sources, named A, B, C, and D, within 1" or 0.018 pc, and have identified signatures of infall and several outflows associated with the core. Aims. We study the interplay between infall and outflow in G31.41+0.31 by investigating at a spatial resolution of a few 100 au their properties and their possible impact on the core. Methods. We carried out molecular line observations of typical high-density tracers, such as CH3CN or H2CO, and shock and outflow tracers, such as SiO, with ALMA at 1.4 mm that achieved an angular resolution of 0.09" (340 au). Results. The observations have revealed inverse P-Cygni profiles in CH3CN and H2CO toward the four sources embedded in the Main core, suggesting that all of them are undergoing collapse. The infall rates, estimated from the red-shifted absorption are on the order of 1E-2 Msun/yr. The individual infall rates imply that the accretion timescale of the Main core is an order of magnitude smaller than its rotation timescale. This confirms that rotating toroids such as the G31 Main core are non-equilibrium, transient collapsing structures that need to be constantly replenished with fresh material from a large-scale reservoir. For sources B, C, and D, the infall could be accelerating inside the sources, while for source A, the presence of a second emission component complicates the interpretation. The SiO observations have revealed the presence of at least six outflows in the G31.41+0.31 star-forming region, ...

preprint2021arXiv

A relation between the radial velocity dispersion of young clusters and their age: Evidence for hardening as the formation scenario of massive close binaries

The majority of massive stars ($>8$ $\rm{M_{\odot}}$) in OB associations are found in close binary systems. Nonetheless, the formation mechanism of these close massive binaries is not understood yet. Using literature data, we measured the radial-velocity dispersion ($σ_\mathrm{RV}$) as a proxy for the close binary fraction in ten OB associations in the Galaxy and the Large Magellanic Cloud, spanning an age range from 1 to 6 Myrs. We find a positive trend of this dispersion with the cluster's age, which is consistent with binary hardening. Assuming a universal binary fraction of $f_\mathrm{bin}$ = 0.7, we converted the $σ_\mathrm{RV}$ behavior to an evolution of the minimum orbital period $P_\mathrm{cutoff}$ from $\sim$9.5 years at 1 Myr to $\sim$1.4 days for the oldest clusters in our sample at $\sim$6 Myr. Our results suggest that binaries are formed at larger separations, and they harden in around 1 to 2 Myrs to produce the period distribution observed in few million year-old OB binaries. Such an inward migration may either be driven by an interaction with a remnant accretion disk or with other young stellar objects present in the system. Our findings constitute the first empirical evidence in favor of migration as a scenario for the formation of massive close binaries.

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

[CII] 158 μm self-absorption and optical depth effects

Context. The [CII] 158 μm far-infrared (FIR) fine-structure line is one of the most important cooling lines of the star-forming interstellar medium (ISM). High spectral resolution observations have shown complex structures in the line profiles of the [CII] emission. Aims. Our aim is to determine whether the complex profiles observed in [^{12}CII] are due to individual velocity components along the line-of-sight or to self-absorption based on a comparison of the [^{12}CII] and isotopic [^{13}CII] line profiles. Methods. Deep integrations with the SOFIA/upGREAT 7-pixel array receiver in M43, Horsehead~PDR, Monoceros~R2, and M17~SW allow for the detection of optically thin [^{13}CII] emission lines, along with the [^{12}CII] emission lines, with a high signal-to-noise ratio. We first derived the [^{12}CII] optical depth and the [CII] column density from a single component model. However, the complex line profiles observed require a double layer model with an emitting background and an absorbing foreground. A multi-component velocity fit allows us to derive the physical conditions of the [CII] gas: column density and excitation temperature. Results. We find moderate to high [^{12}CII] optical depths in all four sources and self-absorption of [^{12}CII] in Mon R2 and M17 SW. The high column density of the warm background emission corresponds to an equivalent Av of up to 41 mag. The foreground absorption requires substantial column densities of cold and dense [CII] gas, with an equivalent Av ranging up to about 13 mag. Conclusions. The column density of the warm background material requires multiple photon-dominated region (PDR) surfaces stacked along the line of sight and in velocity. The substantial column density of dense and cold foreground [CII] gas detected in absorption cannot be explained with any known scenario and we can only speculate on its origins

preprint2020arXiv

A Detailed View of the Circumstellar Environment and Disk of the Forming O-star AFGL 4176

We present a detailed analysis of the disk and circumstellar environment of the forming O-type star AFGL 4176 mm1, placing results from the Atacama Large Millimeter/submillimeter Array (ALMA) into context with multiwavelength data. With ALMA, we detect seventeen 1.2 mm continuum sources within 5$''$ (21,000 au) of AFGL 4176 mm1. We find that mm1 has a spectral index of 3.4$\pm$0.2 across the ALMA band, with $>$87% of its 1.2 mm continuum emission from dust. The source mm2, projected 4200 au from mm1, may be a companion or a blueshifted knot in a jet. We also explore the morphological differences between the molecular lines detected with ALMA, finding 203 lines from 25 molecules, which we categorize into several morphological types. Our results show that AFGL 4176 mm1 provides an example of a forming O-star with a large and chemically complex disk, which is mainly traced by nitrogen-bearing molecules. Lines that show strong emission on the blueshifted side of the disk are predominantly oxygen-bearing, which we suggest are tracing a disk accretion shock. The molecules C$^{34}$S, H$_2$CS and CH$_{3}$CCN trace a slow wide-angle wind or dense structures in the outflow cavity walls. With the Australia Telescope Compact Array (ATCA), we detect a compact continuum source ($<$2000 $\times$ 760 au) at 1.2 cm, associated with mm1, of which $>$96% is from ionized gas. The ATCA NH$_3$(1,1) and (2,2) emission traces a large-scale (r$\sim$0.5 pc) rotating toroid with the disk source mm1 in the blueshifted part of this structure offset to the northwest.

preprint2020arXiv

Autonomous Gaussian decomposition of the Galactic Ring Survey. II. The Galactic distribution of 13CO

Knowledge about the distribution of CO emission in the Milky Way is essential to understand the impact of Galactic environment on the formation and evolution of structures in the interstellar medium. However, currently our insight about the fraction of CO in spiral arm and interarm regions is still limited by large uncertainties in assumed rotation curve models or distance determination techniques. In this work we use a Bayesian approach to obtain the current best assessment of the distribution of 13CO from the Galactic Ring Survey. We performed two different distance estimates that either included or excluded a model for Galactic features. We also include a prior for the solution of the kinematic distance ambiguity that was determined from a compilation of literature distances and an assumed size-linewidth relationship. We find that the fraction of 13CO emission associated with spiral arm features varies from 76% to 84% between the two distance runs. The vertical distribution of the gas is concentrated around the Galactic midplane showing FWHM values of ~75 pc. We do not find any significant difference between gas emission properties associated with spiral arm and interarm features. In particular the distribution of velocity dispersion values of gas emission in spurs and spiral arms is very similar. We detect a trend of higher velocity dispersion values with increasing heliocentric distance, which we attribute to beam averaging effects caused by differences in spatial resolution. We argue that the true distribution of the gas emission is likely more similar to a combination of our two distance results, and highlight the importance of using complementary distance estimations to safeguard against the pitfalls of any single approach. We conclude that the methodology presented in this work is a promising way to determine distances to gas emission features in Galactic plane surveys.

preprint2020arXiv

Characterization of unresolved and unclassified sources detected in radio continuum surveys of the Galactic plane

The continuum emission from 1 to 2 GHz of The HI/OH/Recombination line survey of the inner Milky Way (THOR) at $\lesssim$18" resolution covers $\sim 132$ square degrees of the Galactic plane and detects 10387 sources. Similarly, the first data release of the Global View of Star Formation in the Milky Way (GLOSTAR) surveys covers $\sim 16$ square degrees of the Galactic plane from 4-8 GHz at 18" resolution and detects 1575 sources. However, a large fraction of the unresolved discrete sources detected in these radio continuum surveys of the Galactic plane remain unclassified. Here, we study the Euclidean-normalized differential source counts of unclassified and unresolved sources detected in these surveys and compare them with simulated extragalactic radio source populations as well as previously established source counts. We find that the differential source counts for THOR and GLOSTAR surveys are in excellent agreement with both simulation and previous observations. We also estimate the angular two-point correlation function of unclassified and unresolved sources detected in THOR survey. We find a higher clustering amplitude in comparison with the Faint Images of the Radio Sky at Twenty-cm (FIRST) survey up to the angular separation of $5^{\circ}$. The decrease in angular correlation with increasing flux cut and the excellent agreement of clustering pattern of sources above 1 mJy with high $z$ samples ($z >0.5$) of the FIRST survey indicates that these sources might be high $z$ extragalactic compact objects. The similar pattern of one-point and two-point statistics of unclassified and compact sources with extragalactic surveys and simulations confirms the extragalactic origin of these sources.

preprint2020arXiv

Cloud formation in the atomic and molecular phase: HI self absorption (HISA) towards a Giant Molecular Filament

Molecular clouds form from the atomic phase of the interstellar medium. However, characterizing the transition between the atomic and the molecular interstellar medium (ISM) is a difficult observational task. Here we address cloud formation processes by combining HSIA with molecular line data. One scenario proposed by numerical simulations is that the column density probability density functions (N-PDF) evolves from a log-normal shape at early times to a power-law-like shape at later times. In this paper, we study the cold atomic component of the giant molecular filament GMF38a (d=3.4 kpc, length$\sim230$ pc). We identify an extended HISA feature, which is partly correlated with the 13CO emission. The peak velocities of the HISA and 13CO observations agree well on the eastern side of the filament, whereas a velocity offset of approximately 4 km s$^{-1}$ is found on the western side. The sonic Mach number we derive from the linewidth measurements shows that a large fraction of the HISA, which is ascribed to the cold neutral medium (CNM), is at subsonic and transonic velocities. The column density of the CNM is on the order of 10$^{20}$ to 10$^{21}$ cm$^{-2}$. The column density of molecular hydrogen is an order of magnitude higher. The N-PDFs from HISA (CNM), HI emission (WNM+CNM), and 13CO (molecular component) are well described by log-normal functions, which is in agreement with turbulent motions being the main driver of cloud dynamics. The N-PDF of the molecular component also shows a power law in the high column-density region, indicating self-gravity. We suggest that we are witnessing two different evolutionary stages within the filament. The eastern subregion seems to be forming a molecular cloud out of the atomic gas, whereas the western subregion already shows high column density peaks, active star formation and evidence of related feedback processes.

preprint2020arXiv

Complex organic molecules in low-mass protostars on solar system scales -- I. Oxygen-bearing species

Complex organic molecules (COMs) are thought to form on icy dust grains in the earliest phase of star formation. The evolution of these COMs from the youngest Class 0/I protostellar phases toward the more evolved Class II phase is still not fully understood. Since planet formation seems to start early, and mature disks are too cold for characteristic COM emission lines, studying the inventory of COMs on solar system scales in the Class 0/I stage is relevant. ALMA Band 3 (3 mm) and Band 6 (1 mm) observations are obtained of seven Class 0 protostars in the Perseus and Serpens star-forming regions. By modeling the inner protostellar region using 'LTE' models, the excitation temperature and column densities are determined for several O-bearing COMs. B1-c, B1-bS, and Serpens S68N show COM emission, i.e, three out of the seven sources. No clear correlation seems to exist between the occurrence of COMs and source luminosity. The abundances of several COMs with respect to CH3OH are remarkably similar for the three COM-rich sources, and to IRAS 16293-2422B and HH 212. For other COMs the abundances differ by up to an order of magnitude, indicating that local source conditions are case determining. B1-c hosts a cold ($T_{ex}\approx60$ K), more extended component of COM emission with a column density of typically a few % of the warm/hot ($T_{ex}\sim 200$ K), central component. A D/H ratio of 1-3 % is derived based on the CH2DOH/CH3OH ratio suggesting a temperature of $\sim$15~K during the formation of methanol. This ratio is consistent with other low-mass protostars. Future mid-infrared facilities such as JWST/MIRI will be essential to directly observe COM ices. Combining this with a larger sample of COM-rich sources with ALMA will allow for directly linking ice and gas-phase abundances in order to constrain the routes that produce and maintain chemical complexity during the star formation process.

preprint2020arXiv

Dense Gas in a Giant Molecular Filament

Recent surveys of the Galactic plane in the dust continuum and CO emission lines reveal that large ($\gtrsim 50$~pc) and massive ($\gtrsim 10^5$~$M_\odot$) filaments, know as giant molecular filaments (GMFs), may be linked to galactic dynamics and trace the mid-plane of the gravitational potential in the Milky Way. We have imaged one entire GMF located at $l\sim$52--54$^\circ$ longitude, GMF54 ($\sim$68~pc long), in the empirical dense gas tracers using the HCN(1--0), HNC(1--0), HCO$^+$(1--0) lines, and their $^{13}$C isotopologue transitions, as well as the N$_2$H$^+$(1--0) line. We study the dense gas distribution, the column density probability density functions (N-PDFs) and the line ratios within the GMF. The dense gas molecular transitions follow the extended structure of the filament with area filling factors between 0.06 and 0.28 with respect to $^{13}$CO(1--0). We constructed the N-PDFs of H$_2$ for each of the dense gas tracers based on their column densities and assumed uniform abundance. The N-PDFs of the dense gas tracers appear curved in log-log representation, and the HCO$^+$ N-PDF has the largest log-normal width and flattest power-law slope index. Studying the N-PDFs for sub-regions of GMF54, we found an evolutionary trend in the N-PDFs that high-mass star forming and Photon-Dominate Regions (PDRs) have flatter power-law indices. The integrated intensity ratios of the molecular lines in GMF54 are comparable to those in nearby galaxies. In particular, the N$_2$H$^+$/$^{13}$CO ratio, which traces the dense gas fraction, has similar values in GMF54 and all nearby galaxies except ULIRGs.

preprint2020arXiv

Dynamical cloud formation traced by atomic and molecular gas

Context: Atomic and molecular cloud formation is a dynamical process. However, kinematic signatures of these processes are still observationally poorly constrained. Methods: Targeting the cloud-scale environment of the prototypical infrared dark cloud G28.3, we employ spectral line imaging observations of the two atomic lines HI and [CI] as well as molecular lines observations in 13CO in the 1--0 and 3--2 transitions. The analysis comprises investigations of the kinematic properties of the different tracers, estimates of the mass flow rates, velocity structure functions, a Histogram of Oriented Gradients (HOG) study as well as comparisons to simulations. Results: The central IRDC is embedded in a more diffuse envelope of cold neutral medium (CNM) traced by HI self-absorption (HISA) and molecular gas. The spectral line data as well as the HOG and structure function analysis indicate a possible kinematic decoupling of the HI from the other gas compounds. Spectral analysis and position-velocity diagrams reveal two velocity components that converge at the position of the IRDC. Estimated mass flow rates appear rather constant from the cloud edge toward the center. The velocity structure function analysis is consistent with gas flows being dominated by the formation of hierarchical structures. Conclusions: The observations and analysis are consistent with a picture where the IRDC G28 is formed at the center of two converging gas flows. While the approximately constant mass flow rates are consistent with a self-similar, gravitationally driven collapse of the cloud, external compression by, e.g., spiral arm shocks or supernovae explosions cannot be excluded yet. Future investigations should aim at differentiating the origin of such converging gas flows.

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

Large-scale Map of Millimeter-wavelength Hydrogen Radio Recombination Lines around a Young Massive Star Cluster

We report the first map of large-scale (10 pc in length) emission of millimeter-wavelength hydrogen recombination lines (mm-RRLs) toward the giant H II region around the W43-Main young massive star cluster (YMC). Our mm-RRL data come from the IRAM 30 m telescope and are analyzed together with radio continuum and cm-RRL data from the Karl G. Jansky Very Large Array and HCO$^{+}$ 1-0 line emission data from the IRAM 30 m. The mm-RRLs reveal an expanding wind-blown ionized gas shell with an electron density ~70-1500 cm$^{-3}$ driven by the WR/OB cluster, which produces a total Ly$α$ photon flux of 1.5 x 10$^{50}$ s$^{-1}$. This shell is interacting with the dense neutral molecular gas in the W43-Main dense cloud. Combining the high spectral and angular resolution mm-RRL and cm-RRL cubes, we derive the two-dimensional relative distributions of dynamical and pressure broadening of the ionized gas emission and find that the RRL line shapes are dominated by pressure broadening (4-55 km s$^{-1}$) near the YMC and by dynamical broadening (8-36 km s$^{-1}$) near the shell's edge. Ionized gas clumps hosting ultra-compact H II regions found at the edge of the shell suggest that large-scale ionized gas motion triggers the formation of new star generation near the periphery of the shell.

preprint2020arXiv

Multi-wavelength modelling of the circumstellar environment of the massive proto-star AFGL 2591 VLA 3

We have studied the dust density, temperature and velocity distributions of the archetypal massive young stellar object (MYSO) AFGL 2591. Given its high luminosity ($L=2 \times 10^5$ L$_\odot$) and distance ($d=3.3$ kpc), AFGL 2591 has one of the highest $\sqrt{L}/d$ ratio, giving better resolved dust emission than any other MYSO. As such, this paper provides a template on how to use resolved multi-wavelength data and radiative transfer to obtain a well-constrained 2-D axi-symmetric analytic rotating infall model. We show for the first time that the resolved dust continuum emission from Herschel 70 $μ$m observations is extended along the outflow direction, whose origin is explained in part from warm dust in the outflow cavity walls. However, the model can only explain the kinematic features from CH$_3$CN observations with unrealistically low stellar masses ($<15$ M$_\odot$), indicating that additional physical processes may be playing a role in slowing down the envelope rotation. As part of our 3-step continuum and line fitting, we have identified model parameters that can be further constrained by specific observations. High-resolution mm visibilities were fitted to obtain the disc mass (6 M$_\odot$) and radius (2200 au). A combination of SED and near-IR observations were used to estimate the luminosity and envelope mass together with the outflow cavity inclination and opening angles.

preprint2020arXiv

Spiral arms and instability within the AFGL 4176 mm1 disc

We present high-resolution (30 mas or 130 au at 4.2 kpc) Atacama Large Millimeter/submillimeter Array observations at 1.2 mm of the disc around the forming O-type star AFGL 4176 mm1. The disc (AFGL 4176 mm1-main) has a radius of ~1000 au and contains significant structure, most notably a spiral arm on its redshifted side. We fitted the observed spiral with logarithmic and Archimedean spiral models. We find that both models can describe its structure, but the Archimedean spiral with a varying pitch angle fits its morphology marginally better. As well as signatures of rotation across the disc, we observe gas arcs in CH$_3$CN that connect to other millimetre continuum sources in the field, supporting the picture of interactions within a small cluster around AFGL 4176 mm1-main. Using local thermodynamic equilibrium modelling of the CH$_3$CN K-ladder, we determine the temperature and velocity field across the disc, and thus produce a map of the Toomre stability parameter. Our results indicate that the outer disc is gravitationally unstable and has already fragmented or is likely to fragment in the future, possibly producing further companions. These observations provide evidence that disc fragmentation is one possible pathway towards explaining the high fraction of multiple systems around high-mass stars.

preprint2020arXiv

Ubiquitous velocity fluctuations throughout the molecular interstellar medium

The density structure of the interstellar medium (ISM) determines where stars form and release energy, momentum, and heavy elements, driving galaxy evolution. Density variations are seeded and amplified by gas motion, but the exact nature of this motion is unknown across spatial scale and galactic environment. Although dense star-forming gas likely emerges from a combination of instabilities, convergent flows, and turbulence, establishing the precise origin is challenging because it requires quantifying gas motion over many orders of magnitude in spatial scale. Here we measure the motion of molecular gas in the Milky Way and in nearby galaxy NGC 4321, assembling observations that span an unprecedented spatial dynamic range ($10^{-1}{-}10^3$ pc). We detect ubiquitous velocity fluctuations across all spatial scales and galactic environments. Statistical analysis of these fluctuations indicates how star-forming gas is assembled. We discover oscillatory gas flows with wavelengths ranging from $0.3{-}400$ pc. These flows are coupled to regularly-spaced density enhancements that likely form via gravitational instabilities. We also identify stochastic and scale-free velocity and density fluctuations, consistent with the structure generated in turbulent flows. Our results demonstrate that ISM structure cannot be considered in isolation. Instead, its formation and evolution is controlled by nested, interdependent flows of matter covering many orders of magnitude in spatial scale.

preprint2019arXiv

Feedback in W49A diagnosed with radio recombination lines and models

We present images of radio recombination lines (RRLs) at wavelengths around 17 cm from the star-forming region W49A to determine the kinematics of ionized gas in the THOR survey (The HI/OH/Recombination line survey of the inner Milky Way) at an angular resolution of 16.8"x13.8". The distribution of ionized gas appears to be affected by feedback processes from the star clusters in W49A. The velocity structure of the RRLs shows a complex behavior with respect to the molecular gas. We find a shell-like distribution of ionized gas as traced by RRL emission surrounding the central cluster of OB stars in W49A. We describe the evolution of the shell with the recent feedback model code WARPFIELD that includes the important physical processes and has previously been applied to the 30 Doradus region in the Large Magellanic Cloud. The cloud structure and dynamics of W49A are in agreement with a feedback-driven shell that is re-collapsing. The shell may have triggered star formation in other parts of W49A. We suggest that W49A is a potential candidate for star formation regulated by feedback-driven and re-collapsing shells.

preprint2019arXiv

Probing the cold magnetized Universe with SPICA-POL (B-BOP)

SPICA, the cryogenic infrared space telescope recently pre-selected for a `Phase A' concept study as one of the three remaining candidates for ESA's fifth medium class (M5) mission, is foreseen to include a far-infrared polarimetric imager (SPICA-POL, now called B-BOP), which would offer a unique opportunity to resolve major issues in our understanding of the nearby, cold magnetized Universe. This paper presents an overview of the main science drivers for B-BOP, including high dynamic range polarimetric imaging of the cold interstellar medium (ISM) in both our Milky Way and nearby galaxies. Thanks to a cooled telescope, B-BOP will deliver wide-field 100-350 micron images of linearly polarized dust emission in Stokes Q and U with a resolution, signal-to-noise ratio, and both intensity and spatial dynamic ranges comparable to those achieved by Herschel images of the cold ISM in total intensity (Stokes I). The B-BOP 200 micron images will also have a factor ~30 higher resolution than Planck polarization data. This will make B-BOP a unique tool for characterizing the statistical properties of the magnetized interstellar medium and probing the role of magnetic fields in the formation and evolution of the interstellar web of dusty molecular filaments giving birth to most stars in our Galaxy. B-BOP will also be a powerful instrument for studying the magnetism of nearby galaxies and testing galactic dynamo models, constraining the physics of dust grain alignment, informing the problem of the interaction of cosmic rays with molecular clouds, tracing magnetic fields in the inner layers of protoplanetary disks, and monitoring accretion bursts in embedded protostars.

preprint2019arXiv

Strong excess Faraday rotation on the Inside of the Sagittarius spiral arm

We present first results for Faraday rotation of compact polarized sources (1 to 2 GHz continuum) in The HI/OH/Recombination line (THOR) survey of the inner Galaxy. In the Galactic longitude range 39 degr < l < 52 degr, we find rotation measures in the range -310 rad/m2 < RM < +4219 rad/m2, with the highest values concentrated within a degree of l = 48 degrees at the Sagittarius arm tangent. Most of the high RMs arise in diffuse plasma, along lines of sight that do not intersect HII regions. For l > 49 degr, RM drops off rapidly, while at l < 47 degr, the mean RM is higher with a larger standard deviation than at l > 49 degr. We attribute the RM structure to the compressed diffuse Warm Ionized Medium in the spiral arm, upstream of the major star formation regions. The Sagittarius arm acts as a significant Faraday screen inside the Galaxy. This has implications for models of the Galactic magnetic field and the expected amount of Faraday rotation of Fast Radio Bursts from their host galaxies. We emphasize the importance of sensitivity to high Faraday depth in future polarization surveys.

preprint2019arXiv

The HI/OH/Recombination line survey of the inner Milky Way (THOR): data release 2 and HI overview

With the $Karl~G.~Jansky$ Very Large Array (VLA) in C-configuration, we observed a large portion of the first Galactic quadrant ($l=14.0-67.4^\circ $ and $\lvert b \rvert \leq 1.25^\circ $) achieving an angular resolution of $\leq 40^{\prime\prime}$. At $L$ Band, the WIDAR correlator at the VLA was set to cover the 21~cm HI line, four OH transitions, a series of H$nα$ radio recombination lines (RRLs; $n=151$ to 186), and eight 128~MHz wide continuum spectral windows (SPWs) simultaneously. The HI emission shows clear filamentary substructures at negative velocities with low velocity crowding. The emission at positive velocities is more smeared-out likely due to higher spatial and velocity crowding of structures at the positive velocities. Comparing to the spiral arm model of the Milky Way, the atomic gas follows the Sagittarius and Perseus Arm well but with significant material in the inter-arm regions. With the C-configuration-only HI+continuum data, we produced a HI optical depth map of the THOR areal coverage from 228 absorption spectra with the nearest-neighbor method. With this $τ$ map, we corrected the HI emission for optical depth and the derived column density is 38% higher than the column density with optically thin assumption. The total HI mass with optical depth correction in the survey region is 4.7$\times10^8~M_\odot$, 31% more than the mass derived assuming the emission is optically thin. If we apply this 31% correction to the whole Milky Way, the total atomic gas mass would be 9.4-10.5$\times 10^9~M_\odot$. Comparing the HI with existing CO data, we find a significant increase in the atomic-to-molecular gas ratio from the spiral arms to the inter-arm regions.

preprint2019arXiv

The young stellar content of the giant HII regions M8, G333.6-0.2, and NGC6357 with VLT/KMOS

Context: The identification and characterisation of populations of young massive stars in (giant) HII regions provides important constraints on i) the formation process of massive stars and their early feedback on the environment, and ii) the initial conditions for population synthesis models predicting the evolution of ensembles of stars. Aims: We identify and characterise the stellar populations of the following young giant HII regions: M8, G333.6-0.2, and NGC6357. Methods: We acquired H- and K-band spectra of around 200 stars using The K-band KMOS on the ESO Very Large Telescope. The targets for M8 and NGC6357 were selected from the MYStIX project, which combines X-ray observations with near-infrared and mid-infrared data. For G333.6-0.2, the sample selection is based on the near-infrared colours combined with X-ray data. We introduce an automatic spectral classification method in order to obtain temperatures and luminosities for the observed stars. We analyse the stellar populations using their photometric, astrometric, and spectroscopic properties and compared the position of the stars in the Hertzprung-Russell diagram with stellar evolution models to constrain their ages and mass ranges. Results: We confirm the presence of candidate ionising sources in the three regions and report new ones, including the first spectroscopically identified O stars in G333.6-0.2. In M8 and NGC6357, two populations are identified: i) OB main-sequence stars ($M > 5~\rm{M_{\odot}}$) and ii) pre-main sequence stars ($M\approx0.5-5~\rm{M_{\odot}}$). The ages of the clusters are $\sim$1-3~Myr, $< 3$~Myr, and $\sim$0.5-3~Myr for M8, G333.6-0.2, and NGC6357, respectively. We show that MYStIX selected targets have $>$ 90\% probability of being members of the HII region, whereas a selection based on near infrared (NIR) colours leads to a membership probability of only $\sim$70\%.

preprint2016arXiv

Are Infrared Dark Clouds Really Quiescent?

The dense, cold regions where high-mass stars form are poorly characterised, yet they represent an ideal opportunity to learn more about the initial conditions of high-mass star formation (HMSF), since high-mass starless cores (HMSCs) lack the violent feedback seen at later evolutionary stages. We present continuum maps obtained from the Submillimeter Array (SMA) interferometry at 1.1 mm for four infrared dark clouds (IRDCs, G28.34S, IRDC 18530, IRDC 18306, and IRDC 18308). We also present 1 mm/3 mm line surveys using IRAM 30 m single-dish observations. Our results are: (1) At a spatial resolution of 10^4 AU, the 1.1 mm SMA observations resolve each source into several fragments. The mass of each fragment is on average >10 Msun, which exceeds the predicted thermal Jeans mass of the whole clump by a factor of up to 30, indicating that thermal pressure does not dominate the fragmentation process. Our measured velocity dispersions in the 30 m lines imply that non-thermal motions provides the extra support against gravity in the fragments. (2) Both non-detection of high-J transitions and the hyperfine multiplet fit of N2H+(1-0), C2H(1-0), HCN(1-0), and H13CN(1-0) indicate that our sources are cold and young. However, obvious detection of SiO and the asymmetric line profile of HCO+(1-0) in G28.34S indicate a potential protostellar object and probable infall motion. (3) With a large number of N-bearing species, the existence of carbon rings and molecular ions, and the anti-correlated spatial distributions between N2H+/NH2D and CO, our large-scale high-mass clumps exhibit similar chemical features as small-scale low-mass prestellar objects. This study of a small sample of IRDCs illustrates that thermal Jeans instability alone cannot explain the fragmentation of the clump into cold (~15 K), dense (>10^5 cm-3) cores and that these IRDCs are not completely quiescent.

preprint2016arXiv

Continuum sources from the THOR survey between 1 and 2 GHz

We carried out a large program with the Karl G. Jansky Very Large Array (VLA): "THOR: The HI, OH, Recombination line survey of the Milky Way". We observed a significant portion of the Galactic plane in the first quadrant of the Milky Way in the 21cm HI line, 4 OH transitions, 19 radio recombination lines, and continuum from 1 to 2 GHz. In this paper we present a catalog of the continuum sources in the first half of the survey (l=14.0-37.9deg and l=47.1-51.2deg, |b|<1.1deg) at a spatial resolution of 10-25", with a spatially varying noise level of ~0.3-1 mJy/beam. The catalog contains ~4400 sources. Around 1200 of these are spatially resolved, and ~1000 are possible artifacts, given their low signal-to-noise ratios. Since the spatial distribution of the unresolved objects is evenly distributed and not confined to the Galactic plane, most of them are extragalactic. Thanks to the broad bandwidth of the observations from 1 to 2 GHz, we are able to determine a reliable spectral index for ~1800 sources. The spectral index distribution reveals a double-peaked profile with maxima at spectral indices of alpha = -1 and alpha = 0 , corresponding to steep declining and flat spectra, respectively. This allows us to distinguish between thermal and non-thermal emission, which can be used to determine the nature of each source. We examine the spectral index of ~300 known HII regions, for which we find thermal emission with spectral indices around alpha = 0. In contrast, supernova remnants (SNR) show non-thermal emission with alpha = -0.5 and extragalactic objects generally have a steeper spectral index of alpha = -1. Using the spectral index information of the THOR survey, we investigate potential SNR candidates. We classify the radiation of four SNR candidates as non-thermal, and for the first time, we provide strong evidence for the SNR origin of these candidates.

preprint2016arXiv

Giant molecular filaments in the Milky Way II: The fourth Galactic quadrant

Filamentary structures are common morphological features of the cold, molecular interstellar medium (ISM). Recent studies have discovered massive, hundred-parsec-scale filaments that may be connected to the large-scale, Galactic spiral arm structure. Addressing the nature of these Giant Molecular Filaments (GMFs) requires a census of their occurrence and properties. We perform a systematic search of GMFs in the fourth Galactic quadrant and determine their basic physical properties. We identify GMFs based on their dust extinction signatures in near- and mid-infrared and velocity structure probed by ^{13}CO line emission. We use the ^{13}CO line emission and ATLASGAL dust emission data to estimate the total and dense gas masses of the GMFs. We combine our sample with an earlier sample from literature and study the Galactic environment of the GMFs. We identify nine GMFs in the fourth Galactic quadrant; six are located in the Centaurus spiral arm and three in inter-arm regions. Combining this sample with an earlier study using the same identification criteria in the first Galactic quadrant results in 16 GMFs, nine of which are located within spiral arms. The GMFs have sizes of 80-160 pc and ^{13}CO-derived masses between 5-90 x 10^{4} Msun. Their dense gas mass fractions are between 1.5-37%, being higher in the GMFs connected to spiral arms. We also compare the different GMF-identification methods and find that emission and extinction based techniques overlap only partially, highlighting the need to use both to achieve a complete census.

preprint2016arXiv

Inferring the evolutionary stages of the internal structures of NGC 7538 S and IRS1 from chemistry

To unambiguously diagnose the detailed physical mechanisms and the evolutionary status of high-mass star-forming regions (HMSFRs), we have performed 0.4" ($\sim1000$ AU) resolution observations towards NGC7538S and IRS1, using the Plateau de Bure Interferometer (PdBI). This paper presents a joint analysis of the 1.37 mm continuum emission and the line intensity of 15 molecular species. Assuming local thermal equilibrium, we derived molecular column densities and molecular abundances for each internal gas substructure which is spatially resolved. These derived quantities are compared with a suite of 1-D gas-grain models. NGC7538S is resolved into at least three dense gas condensations. Despite the comparable continuum intensity of these condensations, their differing molecular line emission is suggestive of an overall chemical evolutionary trend from the northeast to the southeast. Since these condensations are embedded within the same parent gas core, their differing chemical properties are the most likely due to the different warm-up histories, rather than the different dynamic timescales. Despite remaining spatially unresolved, in IRS1 we detect abundant complex organic molecules, indicating that IRS1 is the most chemically evolved hot molecular core presented here. We observe a continuum that is dominated by absorption features with at least three strong emission lines, potentially from $\rm CH_3OH$. The $\rm CH_3OH$ lines which are purely in emission have higher excitation than the ones showing purely absorption. Potential reasons for that difference are discussed. This is the first comprehensive comparison of observations of the high-mass cores NGC7538S and IRS1 and a chemical model. We have found that different chemical evolutionary stages can coexist in the same natal gas core. Our achievement illustrates the strength of chemical analysis for understanding HMSFRs.

preprint2016arXiv

Outflow detection in a 70 $μ$m dark high-mass core

We present observations towards a high-mass ($\rm >40\,M_{\odot}$), low luminosity ($\rm <10\,L_{\odot}$) $\rm 70\,μ$m dark molecular core G 28.34 S-A at 3.4 mm, using the IRAM 30 m telescope and the NOEMA interferometer. We report the detection of $\rm SiO$ $J=\rm 2\rightarrow1$ line emission, which is spatially resolved in this source at a linear resolution of $\sim$0.1 pc, while the 3.4 mm continuum image does not resolve any internal sub-structures. The SiO emission exhibits two W-E oriented lobes centring on the continuum peak. Corresponding to the red-shifted and blue-shifted gas with velocities up to $\rm 40\,km\,s^{-1}$ relative to the quiescent cloud, these lobes clearly indicate the presence of a strong bipolar outflow from this $\rm 70\,μ$m dark core, a source previously considered as one of the best candidates of "starless" core. Our SiO detection is consistent with ALMA archival data of $\rm SiO$ $J=\rm 5\rightarrow4$, whose high-velocity blue-shifted gas reveals a more compact lobe spatially closer to the dust center. This outflow indicates that the central source may be in an early evolutionary stage of forming a high-mass protostar. We also find that the low-velocity components (in the range of $\rm V_{lsr}$$\rm_{-5}^{+3}\,km\,s^{-1}$) have an extended, NW-SE oriented distribution. Discussing the possible accretion scenarios of the outflow-powering young stellar object, we argue that the molecular line emission and the molecular outflows may provide a better indication of the accretion history when forming young stellar object, than that from a snapshot observations of the present bolometric luminosity. This is particularly significant for the cases of episodic accretion, which may occur during the collapse of the parent molecular core.

preprint2016arXiv

Resolving the chemical substructure of Orion-KL

The Kleinmann-Low nebula in Orion (Orion-KL) is the nearest example of a high-mass star-forming environment. For the first time, we complemented 1.3 mm Submillimeter Array (SMA) interferometric line survey with IRAM 30 m single-dish observations of the Orion-KL region. Covering a 4 GHz bandwidth in total, this survey contains over 160 emission lines from 20 species (25 isotopologues), including 11 complex organic molecules (COMs). At a spatial resolution of 1200 AU, the continuum substructures are resolved. Extracting the spectra from individual substructures and providing the intensity-integrated distribution map for each species, we studied the small-scale chemical variations in this region. Our main results are: (1) We identify lines from the low-abundance COMs CH3COCH3 and CH3CH2OH, as well as tentatively detect CH3CHO and long carbon-chains C6H and HC7N. (2) We find that while most COMs are segregated by type, peaking either towards the hot core (e.g., N-bearing species) or the compact ridge (e.g., O-bearing species like HCOOCH3 and CH3OCH3), while the distributions of others do not follow this segregated structure (e.g., CH3CH2OH, CH3OH, CH3COCH3). (3) We find a second velocity component of HNCO, SO2, 34SO2, and SO lines, which may be associated with a strong shock event in the low-velocity outflow. (4) Temperatures and molecular abundances show large gradients between central condensations and the outflow regions, illustrating a transition between hot molecular core and shock-chemistry dominated regimes. Our observations of spatially resolved chemical variations in Orion-KL provide the nearest reference source for hot molecular core and outflow chemistry, which will be an important example for interpreting the chemistry of more distant HMSFRs.

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.

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

Chemical evolution in the early phases of massive star formation II: Deuteration

The chemical evolution in high-mass star-forming regions is still poorly constrained. Studying the evolution of deuterated molecules allows to differentiate between subsequent stages of high-mass star formation regions due to the strong temperature dependence of deuterium isotopic fractionation. We observed a sample of 59 sources including 19 infrared dark clouds, 20 high-mass protostellar objects, 11 hot molecular cores and 9 ultra-compact HII regions in the (3-2) transitions of the four deuterated molecules, DCN, DNC, DCO+ and N2D+ as well as their non-deuterated counterpart. The overall detection fraction of DCN, DNC and DCO+ is high and exceeds 50% for most of the stages. N2D+ was only detected in a few infrared dark clouds and high-mass protostellar objects. It can be related to problems in the bandpass at the frequency of the transition and to low abundances in the more evolved, warmer stages. We find median D/H ratios of ~0.02 for DCN, ~0.005 for DNC, ~0.0025 for DCO+ and ~0.02 for N2D+. While the D/H ratios of DNC, DCO+ and N2D+ decrease with time, DCN/HCN peaks at the hot molecular core stage. We only found weak correlations of the D/H ratios for N2D+ with the luminosity of the central source and the FWHM of the line, and no correlation with the H2 column density. In combination with a previously observed set of 14 other molecules (Paper I) we fitted the calculated column densities with an elaborate 1D physico-chemical model with time-dependent D-chemistry including ortho- and para-H2 states. Good overall fits to the observed data have been obtained the model. It is one of the first times that observations and modeling have been combined to derive chemically based best-fit models for the evolution of high-mass star formation including deuteration.

preprint2015arXiv

Filament Fragmentation in High-Mass Star Formation

Aims: We resolve the length-scales for filament formation and fragmentation (res. <=0.1pc), in particular the Jeans length and cylinder fragmentation scale. Methods: We observed the prototypical high-mass star-forming filament IRDC18223 with the Plateau de Bure Interferometer (PdBI) in the 3.2mm continuum and N2H+(1-0) line emission in a ten field mosaic at a spatial resolution of ~4'' (~14000AU). Results: The dust continuum emission resolves the filament into a chain of at least 12 relatively regularly spaced cores. The mean separation between cores is ~0.40(+-0.18)pc. While this is approximately consistent with the fragmentation of an infinite, isothermal, gravitationally bound gas cylinder, a high mass-to-length ratio of M/l~1000M_sun/pc requires additional turbulent and/or magnetic support against radial collapse of the filament. The N2H+(1-0) data reveal a velocity gradient perpendicular to the main filament. Although rotation of the filament cannot be excluded, the data are also consistent with the main filament being comprised of several velocity-coherent sub-filaments. Furthermore, this velocity gradient perpendicular to the filament resembles recent results toward Serpens south that are interpreted as signatures of filament formation within magnetized and turbulent sheet-like structures. Lower-density gas tracers ([CI] and C18O) reveal a similar red/blueshifted velocity structure on scales around 60'' east and west of the IRDC18223 filament. This may tentatively be interpreted as a signature of the large-scale cloud and the smaller-scale filament being kinematically coupled. We do not identify a velocity gradient along the axis of the filament. This may either be due to no significant gas flows along the filamentary axis, but it may partly also be caused by a low inclination angle of the filament with respect to the plane of the sky that could minimize such signature.

preprint2015arXiv

Hierarchical fragmentation and collapse signatures in a high-mass starless region

Aims: Understanding the fragmentation and collapse properties of the dense gas during the onset of high-mass star formation. Methods: We observed the massive (~800M_sun) starless gas clump IRDC18310-4 with the Plateau de Bure Interferometer (PdBI) at sub-arcsecond resolution in the 1.07mm continuum andN2H+(3-2) line emission. Results: Zooming from a single-dish low-resolution map to previous 3mm PdBI data, and now the new 1.07mm continuum observations, the sub-structures hierarchically fragment on the increasingly smaller spatial scales. While the fragment separations may still be roughly consistent with pure thermal Jeans fragmentation, the derived core masses are almost two orders of magnitude larger than the typical Jeans mass at the given densities and temperatures. However, the data can be reconciled with models using non-homogeneous initial density structures, turbulence and/or magnetic fields. While most sub-cores remain (far-)infrared dark even at 70mum, we identify weak 70mum emission toward one core with a comparably low luminosity of ~16L_sun, re-enforcing the general youth of the region. The spectral line data always exhibit multiple spectral components toward each core with comparably small line widths for the individual components (in the 0.3 to 1.0km/s regime). Based on single-dish C18O(2-1) data we estimate a low virial-to-gas-mass ratio <=0.25. We discuss that the likely origin of these spectral properties may be the global collapse of the original gas clump that results in multiple spectral components along each line of sight. Even within this dynamic picture the individual collapsing gas cores appear to have very low levels of internal turbulence.

preprint2015arXiv

Kinematic and Thermal Structure at the onset of high-mass star formation

We want to understand the kinematic and thermal properties of young massive gas clumps prior to and at the earliest evolutionary stages of high-mass star formation. Do we find signatures of gravitational collapse? Do we find temperature gradients in the vicinity or absence of infrared emission sources? Do we find coherent velocity structures toward the center of the dense and cold gas clumps? To determine kinematics and gas temperatures, we used ammonia, because it is known to be a good tracer and thermometer of dense gas. We observed the NH$_3$(1,1) and (2,2) lines within seven very young high-mass star-forming regions with the VLA and the Effelsberg 100m telescope. This allows us to study velocity structures, linewidths, and gas temperatures at high spatial resolution of 3-5$"$, corresponding to $\sim$0.05 pc. We find on average cold gas clumps with temperatures in the range between 10 K and 30 K. The observations do not reveal a clear correlation between infrared emission peaks and ammonia temperature peaks. We report an upper limit for the linewidth of $\sim$1.3 km s$^{-1}$, at the spectral resolution limit of our VLA observation. This indicates a relatively low level of turbulence on the scale of the observations. Velocity gradients are present in almost all regions with typical velocity differences of 1 to 2 km s$^{-1}$ and gradients of 5 to 10 km s$^{-1}$ pc$^{-1}$. These velocity gradients are smooth in most cases, but there is one exceptional source (ISOSS23053), for which we find several velocity components with a steep velocity gradient toward the clump centers that is larger than 30 km s$^{-1}$ pc$^{-1}$. This steep velocity gradient is consistent with recent models of cloud collapse. Furthermore, we report a spatial correlation of ammonia and cold dust, but we also find decreasing ammonia emission close to infrared emission sources.

preprint2015arXiv

Relationship between the column density distribution and evolutionary class of molecular clouds as viewed by ATLASGAL

We present the first study of the relationship between the column density distribution of molecular clouds within nearby Galactic spiral arms and their evolutionary status as measured from their stellar content. We analyze a sample of 195 molecular clouds located at distances below 5.5 kpc, identified from the ATLASGAL 870 micron data. We define three evolutionary classes within this sample: starless clumps, star-forming clouds with associated young stellar objects, and clouds associated with HII regions. We find that the N(H2) probability density functions (N-PDFs) of these three classes of objects are clearly different: the N-PDFs of starless clumps are narrowest and close to log-normal in shape, while star-forming clouds and HII regions exhibit a power-law shape over a wide range of column densities and log-normal-like components only at low column densities. We use the N-PDFs to estimate the evolutionary time-scales of the three classes of objects based on a simple analytic model from literature. Finally, we show that the integral of the N-PDFs, the dense gas mass fraction, depends on the total mass of the regions as measured by ATLASGAL: more massive clouds contain greater relative amounts of dense gas across all evolutionary classes.

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 disk-outflow system in the S255IR area of high mass star formation

We report the results of our observations of the S255IR area with the SMA at 1.3 mm in the very extended configuration and at 0.8 mm in the compact configuration as well as with the IRAM-30m at 0.8 mm. The best achieved angular resolution is about 0.4 arcsec. The dust continuum emission and several tens of molecular spectral lines are observed. The majority of the lines is detected only towards the S255IR-SMA1 clump, which represents a rotating structure (probably disk) around the young massive star. The achieved angular resolution is still insufficient for conclusions about Keplerian or non-Keplerian character of the rotation. The temperature of the molecular gas reaches 130-180 K. The size of the clump is about 500 AU. The clump is strongly fragmented as follows from the low beam filling factor. The mass of the hot gas is significantly lower than the mass of the central star. A strong DCN emission near the center of the hot core most probably indicates a presence of a relatively cold ($\lesssim 80$ K) and rather massive clump there. High velocity emission is observed in the CO line as well as in lines of high density tracers HCN, HCO+, CS and other molecules. The outflow morphology obtained from combination of the SMA and IRAM-30m data is significantly different from that derived from the SMA data alone. The CO emission detected with the SMA traces only one boundary of the outflow. The outflow is most probably driven by jet bow shocks created by episodic ejections from the center. We detected a dense high velocity clump associated apparently with one of the bow shocks. The outflow strongly affects the chemical composition of the surrounding medium.

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 Musca cloud: A 6 pc-long velocity-coherent, sonic filament

Filaments play a central role in the molecular clouds' evolution, but their internal dynamical properties remain poorly characterized. To further explore the physical state of these structures, we have investigated the kinematic properties of the Musca cloud. We have sampled the main axis of this filamentary cloud in $^{13}$CO and C$^{18}$O (2--1) lines using APEX observations. The different line profiles in Musca shows that this cloud presents a continuous and quiescent velocity field along its $\sim$6.5 pc of length. With an internal gas kinematics dominated by thermal motions (i.e., $σ_{NT}/c_s\lesssim1$) and large-scale velocity gradients, these results reveal Musca as the longest velocity-coherent, sonic-like object identified so far in the ISM. The transonic properties of Musca present a clear departure from the predicted supersonic velocity dispersions expected in the Larson's velocity dispersion-size relationship, and constitute the first observational evidence of a filament fully decoupled from the turbulent regime over multi-parsec scales.

preprint2015arXiv

THOR - The HI, OH, Recombination Line Survey of the Milky Way - The pilot study: HI observations of the giant molecular cloud W43

To study the atomic, molecular and ionized emission of Giant Molecular Clouds (GMCs), we have initiated a Large Program with the VLA: 'THOR - The HI, OH, Recombination Line survey of the Milky Way'. We map the 21cm HI line, 4 OH lines, 19 H_alpha recombination lines and the continuum from 1 to 2 GHz of a significant fraction of the Milky Way (l=15-67deg, |b|<1deg) at ~20" resolution. In this paper, we focus on the HI emission from the W43 star-formation complex. Classically, the HI 21cm line is treated as optically thin with column densities calculated under this assumption. This might give reasonable results for regions of low-mass star-formation, however, it is not sufficient to describe GMCs. We analyzed strong continuum sources to measure the optical depth, and thus correct the HI 21cm emission for optical depth effects and weak diffuse continuum emission. Hence, we are able to measure the HI mass of W43 more accurately and our analysis reveals a lower limit of M=6.6x10^6 M_sun, which is a factor of 2.4 larger than the mass estimated with the assumption of optically thin emission. The HI column densities are as high as N(HI)~150 M_sun/pc^2 ~ 1.9x10^22 cm^-2, which is an order of magnitude higher than for low mass star formation regions. This result challenges theoretical models that predict a threshold for the HI column density of ~10 M_sun/pc^2, at which the formation of molecular hydrogen should set in. By assuming an elliptical layered structure for W43, we estimate the particle density profiles. While at the cloud edge atomic and molecular hydrogen are well mixed, the center of the cloud is strongly dominated by molecular hydrogen. We do not identify a sharp transition between hydrogen in atomic and molecular form. Our results are an important characterization of the atomic to molecular hydrogen transition in an extreme environment and challenges current theoretical models.

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-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

Carbon in different phases ([CII], [CI], and CO) in infrared dark clouds: Cloud formation signatures and carbon gas fractions

Context: How do molecular clouds form out of the atomic phase? And what are the relative fractions of carbon in the ionized, atomic and molecular phase? These are questions at the heart of cloud and star formation. Methods: Using multiple observatories from Herschel and SOFIA to APEX and the IRAM 30m telescope, we mapped the ionized, atomic and molecular carbon ([CII]@1900GHz, [CI]@492GHz and C18O(2-1)@220GHz) at high spatial resolution (12"-25") in four young massive infrared dark clouds (IRDCs). Results: The three carbon phases were successfully mapped in all four regions, only in one source the [CII] line remained a non-detection. Both the molecular and atomic phases trace the dense structures well, with [CI] also tracing material at lower column densities. [CII] exhibits diverse morphologies in our sample, from compact to diffuse structures probing the cloud environment. In at least two out of the four regions, we find kinematic signatures strongly indicating that the dense gas filaments have formed out of a dynamically active and turbulent atomic/molecular cloud, potentially from converging gas flows. The atomic-to-molecular carbon gas mass ratios are low between 7% and 12% with the lowest values found toward the most quiescent region. In the three regions where [CII] is detected, its mass is always higher by a factor of a few than that of the atomic carbon. The ionized carbon emission depends as well on the radiation field, however, we also find strong [CII] emission in a region without significant external sources, indicating that other processes, e.g., energetic gas flows can contribute to the [CII] excitation as well.

preprint2014arXiv

Chemical evolution in the early phases of massive star formation. I

Understanding the chemical evolution of young (high-mass) star-forming regions is a central topic in star formation research. Chemistry is employed as a unique tool 1) to investigate the underlying physical processes and 2) to characterize the evolution of the chemical composition. We observed a sample of 59 high-mass star-forming regions at different evolutionary stages varying from the early starless phase of infrared dark clouds to high-mass protostellar objects to hot molecular cores and, finally, ultra-compact HII regions at 1mm and 3mm with the IRAM 30m telescope. We determined their large-scale chemical abundances and found that the chemical composition evolves along with the evolutionary stages. On average, the molecular abundances increase with time. We modeled the chemical evolution, using a 1D physical model where density and temperature vary from stage to stage coupled with an advanced gas-grain chemical model and derived the best-fit chi^2 values of all relevant parameters. A satisfying overall agreement between observed and modeled column densities for most of the molecules was obtained. With the best-fit model we also derived a chemical age for each stage, which gives the timescales for the transformation between two consecutive stages. The best-fit chemical ages are ~10,000 years for the IRDC stage, ~60,000 years for the HMPO stage, ~40,000 years for the HMC stage, and ~10,000 years for the UCHII stage. The total chemical timescale for the entire evolutionary sequence of the high-mass star formation process is on the order of 10^5 years, which is consistent with theoretical estimates. Furthermore, based on the approach of a multiple-line survey of unresolved data, we were able to constrain an intuitive and reasonable physical and chemical model. The results of this study can be used as chemical templates for the different evolutionary stages in high-mass star formation.

preprint2014arXiv

Kinematic structure of massive star-forming regions - I. Accretion along filaments

The mid- and far-infrared view on high-mass star formation, in particular with the results from the Herschel space observatory, has shed light on many aspects of massive star formation. However, these continuum studies lack kinematic information. We study the kinematics of the molecular gas in high-mass star-forming regions. We complemented the PACS and SPIRE far-infrared data of 16 high-mass star-forming regions from the Herschel key project EPoS with N2H+ molecular line data from the MOPRA and Nobeyama 45m telescope. Using the full N2H+ hyperfine structure, we produced column density, velocity, and linewidth maps. These were correlated with PACS 70micron images and PACS point sources. In addition, we searched for velocity gradients. For several regions, the data suggest that the linewidth on the scale of clumps is dominated by outflows or unresolved velocity gradients. IRDC18454 and G11.11 show two velocity components along several lines of sight. We find that all regions with a diameter larger than 1pc show either velocity gradients or fragment into independent structures with distinct velocities. The velocity profiles of three regions with a smooth gradient are consistent with gas flows along the filament, suggesting accretion flows onto the densest regions. We show that the kinematics of several regions have a significant and complex velocity structure. For three filaments, we suggest that gas flows toward the more massive clumps are present.

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

The dynamics and star-forming potential of the massive Galactic centre cloud G0.253+0.016

The massive infrared dark cloud G0.253+0.016 projected 45pc from the Galactic centre contains ~10^5Msun of dense gas whilst being mostly devoid of observed star-formation tracers. To scrutinise the physical properties, dynamics and structure of this cloud with reference to its star-forming potential, we have carried out a concerted SMA and IRAM 30m study of this cloud in dust continuum, CO isotopologues, shock tracing molecules, as well as H$_2$CO to trace the gas temperature. We detect and characterise the dust cores within G0.253+0.016 at ~1.3 mm and find that the kinetic temperature of the gas is >320K on size-scales of ~0.15 pc. Analysis of the position-velocity diagrams of our observed lines show broad linewidths and strong shock emission in the south of the cloud, indicating that G0.253+0.016 is colliding with another cloud at v(LSR)~70 km/s. We confirm via an analysis of the observed dynamics in the CMZ that it is an elongated structure, orientated with Sgr B2 closer to the Sun than Sgr A*, however our results suggest that the actual geometry may be more complex than an elliptical ring. We find that the column density PDF of G0.253+0.016 is log-normal with no discernible power-law tail, consistent with little star formation, and that its width can be explained in the framework of theory predicting the density structure of clouds created by supersonic, magnetised turbulence. We also present the delta-variance spectrum of this region, and show it is consistent with that expected for clouds with no star formation. Using G0.253+0.016 as a test-bed of the conditions required for star formation in a different physical environment to that of nearby clouds, we also conclude that there is not one column density threshold for star formation, but instead this value is dependant on the local physical conditions. [Abbrv.]

preprint2014arXiv

The Earliest Phases of Star formation (EPoS): Temperature, density, and kinematic structure of the star-forming core CB 17

Context: The initial conditions for the gravitational collapse of molecular cloud cores and the subsequent birth of stars are still not well constrained. The characteristic cold temperatures (about 10 K) in such regions require observations at sub-millimetre and longer wavelengths. The Herschel Space Observatory and complementary ground-based observations presented in this paper have the unprecedented potential to reveal the structure and kinematics of a prototypical core region at the onset of stellar birth. Aims: This paper aims to determine the density, temperature, and velocity structure of the star-forming Bok globule CB 17. This isolated region is known to host (at least) two sources at different evolutionary stages: a dense core, SMM1, and a Class I protostar, IRS. Methods: We modeled the cold dust emission maps from 100 micron to 1.2 mm with both a modified blackbody technique to determine the optical depth-weighted line-of-sight temperature and column density and a ray-tracing technique to determine the core temperature and volume density structure. Furthermore, we analysed the kinematics of CB17 using the high-density gas tracer N2H+. Results: From the ray-tracing analysis, we find a temperature in the centre of SMM1 of 10.6 K, a flat density profile with radius 9500 au, and a central volume density of n(H) = 2.3x10^5 cm-3. The velocity structure of the N2H+ observations reveal global rotation with a velocity gradient of 4.3 km/s/pc. Superposed on this rotation signature we find a more complex velocity field, which may be indicative of differential motions within the dense core. Conclusions: SMM is a core in an early evolutionary stage at the verge of being bound, but the question of whether it is a starless or a protostellar core remains unanswered.

preprint2014arXiv

The formation of the W43 complex: constraining its atomic-to-molecular transition and searching for colliding clouds

Numerical simulations have explored the possibility to form molecular clouds through either a quasi-static, self-gravitating mechanism or the collision of gas streams or lower-density clouds. They also quantitatively predict the distribution of matter at the transition from atomic to molecular gases. We aim to observationally test these models by studying the environment of W43, a molecular cloud complex near the tip of the Galactic long bar. Using Galaxy-wide HI and 12CO surveys we searched for gas flowing toward the W43 molecular cloud complex. We also estimated the HI and H2 mass surface densities to constrain the transition from atomic to molecular gas around and within W43. We found 3 cloud ensembles within the position-velocity diagrams of 12CO and HI gases. They are separated by 20km/s along the line of sight and extend into the 13CO velocity structure of W43. Since their velocity gradients are consistent with free-fall, they could be nearby clouds attracted by, and streaming toward, the W43 10^7Msun potential well. We show that the HI surface density, Sigma_HI=45-85Msun/pc2, does not reach any threshold level but increases when entering the 130pc-wide molecular complex previously defined. This suggests that an equilibrium between H2 formation and photodissociation has not yet been reached. The H2-to-HI ratio measured over the W43 region and its surroundings, R_H2~3.5, is high, indicating that most of the gas is already in molecular form in W43 and in structures several hundreds of parsecs downstream along the Scutum-Centaurus arm. The W43 molecular cloud complex may have formed, and in fact may still be accreting mass from the agglomeration of clouds. Already in the molecular-dominated regime, most of these clouds are streaming from the Scutum-Centaurus arm. This is in clear disagreement with quasi-static and steady-state models of molecular cloud formation.

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

Characterization of Infrared Dark Clouds -- NH$_3$ Observations of an Absorption-contrast Selected IRDC Sample

Despite increasing research in massive star formation, little is known about its earliest stages. Infrared Dark Clouds (IRDCs) are cold, dense and massive enough to harbour the sites of future high-mass star formation. But up to now, mainly small samples have been observed and analysed. To understand the physical conditions during the early stages of high-mass star formation, it is necessary to learn more about the physical conditions and stability in relatively unevolved IRDCs. Thus, for characterising IRDCs studies of large samples are needed. We investigate a complete sample of 218 northern hemisphere high-contrast IRDCs using the ammonia (1,1)- and (2,2)-inversion transitions. We detected ammonia (1,1)-inversion transition lines in 109 of our IRDC candidates. Using the data we were able to study the physical conditions within the star-forming regions statistically. We compared them with the conditions in more evolved regions which have been observed in the same fashion as our sample sources. Our results show that IRDCs have, on average, rotation temperatures of 15 K, are turbulent (with line width FWHMs around 2 km s$^{-1}$), have ammonia column densities on the order of $10^{14}$ cm$^{-2}$ and molecular hydrogen column densities on the order of $10^{22}$ cm$^{-2}$. Their virial masses are between 100 and a few 1000 M$_\odot$. The comparison of bulk kinetic and potential energies indicate that the sources are close to virial equilibrium. IRDCs are on average cooler and less turbulent than a comparison sample of high-mass protostellar objects, and have lower ammonia column densities. Virial parameters indicate that the majority of IRDCs are currently stable, but are expected to collapse in the future.

preprint2013arXiv

Deep near-infrared imaging of W3 Main: constraints on stellar cluster formation

Embedded clusters like W3 Main are complex and dynamically evolving systems that represent an important phase of the star formation process. We aim at the characterization of the entire stellar content of W3 Main in a statistical sense to identify possible differences in evolutionary phase of the stellar populations and find clues about the formation mechanism of this massive embedded cluster. Methods. Deep JHKs imaging is used to derive the disk fraction, Ks-band luminosity functions and mass functions for several subregions in W3 Main. A two dimensional completeness analysis using artificial star experiments is applied as a crucial ingredient to assess realistic completeness limits for our photometry. We find an overall disk fraction of 7.7 $\pm$ 2.3%, radially varying from 9.4 $\pm$ 3.0 % in the central 1 pc to 5.6 $\pm$ 2.2 % in the outer parts of W3 Main. The mass functions derived for three subregions are consistent with a Kroupa and Chabrier mass function. The mass function of IRSN3 is complete down to 0.14 Msun and shows a break at M $\sim$ 0.5 Msun. We interpret the higher disk fraction in the center as evidence for a younger age of the cluster center. We find that the evolutionary sequence observed in the low-mass stellar population is consistent with the observed age spread among the massive stars. An analysis of the mass function variations does not show evidence for mass segregation. W3 Main is currently still actively forming stars, showing that the ionizing feedback of OB stars is confined to small areas ($\sim$ 0.5 pc). The FUV feedback might be influencing large regions of the cluster as suggested by the low overall disk fraction.

preprint2013arXiv

Fragmentation and dynamical collapse of the starless high-mass star-forming region IRDC18310-4

Aims: We study the fragmentation and dynamical properties of a massive starless gas clump at the onset of high-mass star formation. Methods: Based on Herschel continuum data we identify a massive gas clump that remains far-infrared dark up to 100mum wavelengths. The fragmentation and dynamical properties are investigated by means of Plateau de Bure Interferometer and Nobeyama 45m single-dish spectral line and continuum observations. Results: The massive gas reservoir fragments at spatial scales of ~18000AU in four cores. Comparing the spatial extent of this high-mass region with intermediate- to low-mass starless cores from the literature, we find that linear sizes do not vary significantly over the whole mass regime. However, the high-mass regions squeeze much more gas into these similar volumes and hence have orders of magnitude larger densities. The fragmentation properties of the presented low-to high-mass regions are consistent with gravitational instable Jeans fragmentation. Furthermore, we find multiple velocity components associated with the resolved cores. Recent radiative transfer hydrodynamic simulations of the dynamic collapse of massive gas clumps also result in multiple velocity components along the line of sight because of the clumpy structure of the regions. This result is supported by a ratio between viral and total gas mass for the whole region <1. Conclusions: This apparently still starless high-mass gas clump exhibits clear signatures of early fragmentation and dynamic collapse prior to the formation of an embedded heating source. A comparison with regions of lower mass reveals that the linear size of star-forming regions does not necessarily have to vary much for different masses, however, the mass reservoirs and gas densities are orders of magnitude enhanced for high-mass regions compared to their lower-mass siblings.

preprint2013arXiv

Fragmentation, infall, and outflow around the showcase massive protostar NGC7538 IRS1 at 500 AU resolution

Aims: Revealing the fragmentation, infall, and outflow processes in the immediate environment around massive young stellar objects is crucial for understanding the formation of the most massive stars. Methods: With this goal in mind we present the so far highest spatial-resolution thermal submm line and continuum observations toward the young high-mass protostar NGC7538 IRS1. Using the Plateau de Bure Interferometer in its most extended configuration at 843mum wavelength, we achieved a spatial resolution of 0.2"x0.17", corresponding to ~500AU at a distance of 2.7\,kpc. Results: For the first time, we have observed the fragmentation of the dense inner core of this region with at least three subsources within the inner 3000 AU. The outflow exhibits blue- and red-shifted emission on both sides of the central source indicating that the current orientation has to be close to the line-of-sight, which differs from other recent models. We observe rotational signatures in northeast-southwest direction; however, even on scales of 500 AU, we do not identify any Keplerian rotation signatures. This implies that during the early evolutionary stages any stable Keplerian inner disk has to be very small (<=500 AU). The high-energy line HCN(4-3)v2=1 (E_u/k=1050K) is detected over an extent of approximately 3000 AU. In addition to this, the detection of red-shifted absorption from this line toward the central dust continuum peak position allows us to estimate infall rates of ~1.8x10^(-3)Msun/yr on the smallest spatial scales. Although all that gas will not necessarily be accreted onto the central protostar, nevertheless, such inner core infall rates are among the best proxies of the actual accretion rates one can derive during the early embedded star formation phase. These data are consistent with collapse simulations and the observed high multiplicity of massive stars.

preprint2013arXiv

Large scale IRAM 30m CO-observations in the giant molecular cloud complex W43

We aim to give a full description of the distribution and location of dense molecular clouds in the giant molecular cloud complex W43. It has previously been identified as one of the most massive star-forming regions in our Galaxy. To trace the moderately dense molecular clouds in the W43 region, we initiated an IRAM 30m large program, named W43-HERO, covering a large dynamic range of scales (from 0.3 to 140 pc). We obtained on-the-fly-maps in 13CO (2-1) and C18O (2-1) with a high spectral resolution of 0.1 km/s and a spatial resolution of 12". These maps cover an area of ~1.5 square degrees and include the two main clouds of W43, as well as the lower density gas surrounding them. A comparison with Galactic models and previous distance calculations confirms the location of W43 near the tangential point of the Scutum arm at a distance from the Sun of approximately 6 kpc. The resulting intensity cubes of the observed region are separated into sub-cubes, centered on single clouds which are then analyzed in detail. The optical depth, excitation temperature, and H2 column density maps are derived out of the 13CO and C18O data. These results are then compared with those derived from Herschel dust maps. The mass of a typical cloud is several 10^4 solar masses while the total mass in the dense molecular gas (>100 cm^-3) in W43 is found to be about 1.9e6 solar masses. Probability distribution functions obtained from column density maps derived from molecular line data and Herschel imaging show a log-normal distribution for low column densities and a power-law tail for high densities. A flatter slope for the molecular line data PDF may imply that those selectively show the gravitationally collapsing gas.

preprint2013arXiv

The earliest phases of star formation - A Herschel key project. The thermal structure of low-mass molecular cloud cores

The temperature and density structure of molecular cloud cores are the most important physical quantities that determine the course of the protostellar collapse and the properties of the stars they form. Nevertheless, density profiles often rely either on the simplifying assumption of isothermality or on observationally poorly constrained model temperature profiles. With the aim of better constraining the initial physical conditions in molecular cloud cores at the onset of protostellar collapse, we initiated the Guaranteed Time Key Project (GTKP) "The Earliest Phases of Star Formation" (EPoS) with the Herschel satellite. This paper gives an overview of the low-mass sources in the EPoS project, including all observations, the analysis method, and the initial results of the survey. We study the thermal dust emission of 12 previously well-characterized, isolated, nearby globules using FIR and submm continuum maps at up to eight wavelengths between 100 micron and 1.2 mm. Our sample contains both globules with starless cores and embedded protostars at different early evolutionary stages. The dust emission maps are used to extract spatially resolved SEDs, which are then fit independently with modified blackbody curves to obtain line-of-sight-averaged dust temperature and column density maps. We find that the thermal structure of all globules is dominated by external heating from the interstellar radiation field and moderate shielding by thin extended halos. All globules have warm outer envelopes (14-20 K) and colder dense interiors (8-12 K). The protostars embedded in some of the globules raise the local temperature of the dense cores only within radii out to about 5000 AU, but do not significantly affect the overall thermal balance of the globules.

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

Fragmentation in the Massive Star-Forming Region IRAS 19410+2336

The Core Mass Functions (CMFs) of low-mass star-forming regions are found to resemble the shape of the Initial Mass Function (IMF). A similar result is observed for the dust clumps in high-mass star forming regions, although at spatial scales of clusters that do not resolve the substructure found in them. The region IRAS 19410+2336 is one exception, having been observed at spatial scales on the order of ~2500AU, resolving the clump substructure into individual cores. We mapped that region with the PdBI in the 1.4mm and 3mm continuum and several transitions of H2CO and CH3CN. The H2CO transitions were also observed with the IRAM 30m Telescope. We detected 26 continuum sources at 1.4mm with a spatial resolution down to ~2200 AU, distributed in two protoclusters. With the lines emission we derived the temperature structure of the region, ranging from 35 to 90K. With them we calculated the core masses of the detected sources, ranging from ~0.7 to ~8 M_sun. These masses were strongly (~90%) affected by the interferometer spatial filtering. Considering only the detected dense cores we derived a CMF with a power-law index b=-2.3+-0.2. We resolve the Jeans length of the protoclusters by one order of magnitude, and only find little velocity dispersion between the different subsources. Since we cannot unambiguously differentiate protostellar and prestellar cores, the derived CMF is not prestellar. Also, because of the large missing flux, we cannot establish a firm link between the CMF and the IMF. This implies that future high-mass CMF studies will need to complement the interferometer continuum data with the short spacing data, a task suitable for ALMA. We note that the method of extracting temperatures using H2CO lines becomes less applicable when reaching the dense core scales of the interferometric observations because most of the H2CO appears to originate in the envelope structure.

preprint2012arXiv

Multiple episodes of star formation in the CN15/16/17 molecular complex

We have started a campaign to identify massive star clusters inside bright molecular bubbles towards the Galactic Center. The CN15/16/17 molecular complex is the first example of our study. The region is characterized by the presence of two young clusters, DB10 and DB11, visible in the NIR, an ultra-compact HII region identified in the radio, several young stellar objects visible in the MIR, a bright diffuse nebulosity at 8μm coming from PAHs and sub-mm continuum emission revealing the presence of cold dust. Given its position on the sky (l=0.58, b=-0.85) and its kinematic distance of ~7.5 kpc, the region was thought to be a very massive site of star formation in proximity of the CMZ. The cluster DB11 was estimated to be as massive as 10^4 M_sun. However the region's properties were known only through photometry and its kinematic distance was very uncertain given its location at the tangential point. We aimed at better characterizing the region and assess whether it could be a site of massive star formation located close to the Galactic Center. We have obtained NTT/SofI JHKs photometry and long slit K band spectroscopy of the brightest members. We have additionally collected data in the radio, sub-mm and mid infrared, resulting in a quite different picture of the region. We have confirmed the presence of massive early B type stars and have derived a spectro-photometric distance of ~1.2 kpc, much smaller than the kinematic distance. Adopting this distance we obtain clusters masses of M(DB10) ~ 170 M_sun and M(DB11) ~ 275 M_sun. This is consistent with the absence of any O star, confirmed by the excitation/ionization status of the nebula. No HeI diffuse emission is detected in our spectroscopic observations at 2.113μm, which would be expected if the region was hosting more massive stars. Radio continuum measurements are also consistent with the region hosting at most early B stars.

preprint2012arXiv

Search for starless clumps in the ATLASGAL survey

In this study, we present an unbiased sample of the earliest stages of massive star formation across 20 square-degree of the sky. Within the region 10deg < l < 20deg and |b| < 1deg, we search the ATLASGAL survey at 870 micron for dense gas condensations. These clumps are carefully examined for indications of ongoing star formation using YSOs from the GLIMPSE source catalog as well as sources in the 24 micron MIPSGAL images, to search for starless clumps. We calculate the column densities as well as the kinematic distances and masses for sources where the v_lsr is known from spectroscopic observations. Within the given region, we identify 210 starless clumps with peak column densities > 1 x 10e23 cm^(-2). In particular, we identify potential starless clumps on the other side of the Galaxy. The sizes of the clumps range between 0.1 pc and 3 pc with masses between a few tens of solar masses up to several ten thousands of solar masses. Most of them may form massive stars, but in the 20 square-degree we only find 14 regions massive enough to form stars more massive than 20 solar masses and 3 regions with the potential to form stars more massive than 40 40 solar masses. The slope of the high-mass tail of the clump mass function for clumps on the near side of the Galaxy is 2.2 and, therefore, Salpeter-like. We estimate the lifetime of the most massive starless clumps to be 60000 yr. The sample offers a uniform selection of starless clumps. In the large area surveyed, we only find a few potential precursors of stars in the excess of 40 solar masses. It appears that the lifetime of these clumps is somewhat shorter than their free-fall times, although both values agree within the errors. In addition, these are ideal objects for detailed studies and follow-up observations.

preprint2012arXiv

The Earliest Phases of Star formation observed with Herschel (EPoS): The dust temperature and density distributions of B68

(Abriged) In the framework of the Herschel GTKP "The earliest phases of star formation", we have imaged B68 between 100 and 500 um. Ancillary (sub)mm data, spectral line maps of the 12/13CO(2-1) transitions as well as a NIR extinction map were added to the analysis. We employed a ray-tracing algorithm to derive the 2D mid-plane dust temperature and volume density distribution without suffering from LoS averaging effects of simple SED fitting procedures. Additional 3D radiative transfer calculations were employed to investigate the connection between the external irradiation and the peculiar crescent shaped morphology found in the FIR maps. For the first time, we spatially resolve the dust temperature and density distribution of B68. We find T_dust dropping from 16.7 K at the edge to 8.2 K in the centre, which is about 4 K lower than the result of the simple SED fitting approach. N_H peaks at 4.3x10^22 cm^-2 and n_H at 3.4x10^5 cm^-3 in the centre. B68 has a mass of 3.1 M_sun of material with A_K > 0.2 mag for an assumed distance of 150 pc. We detect a compact source in the southeastern trunk, which is also seen in extinction and CO. We find the radial density distribution from the edge of the inner plateau outward to be n_H ~ r^-3.5. Such a steep profile can arise from either or both of the following: external irradiation with a significant UV contribution or the fragmentation of filamentary structures. Our 3D radiative transfer model of an externally irradiated core by an anisotropic ISRF reproduces the crescent morphology. Our CO observations show that B68 is part of a chain of globules in both space and velocity, which may indicate that it was once part of a filament which dispersed. We also resolve a new compact source in the SE trunk and find that it is slightly shifted in centroid velocity from B68, lending qualitative support to core collision scenarios.

preprint2012arXiv

Triggered/sequential star formation? A multi-phase ISM study around the prominent IRDC G18.93-0.03

G18.93-0.03 is a prominent dust complex within an 0.8deg long filament, with the molecular clump G18.93/m being IR dark from near IR wavelength up to 160mu. Spitzer composite images show an IR bubble spatially associated with G18.93. We use GRS 13CO and IRAM 30m H13CO+ data to disentangle the spatial structure of the region. From ATLASGAL submm data we calculate the gas mass, while we use the H13CO+ line width to estimate its virial mass. Using HERSCHEL data we produce temperature maps from fitting the SED. With the MAGPIS 20cm and SuperCOSMOS Halpha data we trace the ionized gas, and the VGPS HI survey provides information on the atomic hydrogen gas. We show that the bubble is spatially associated with G18.93, located at a kinematic near distance of 3.6kpc. With 280Msun, the most massive clump within G18.93 is G18.93/m. The virial analysis shows that it may be gravitationally bound and has neither Spitzer young stellar objects nor mid-IR point sources within. Fitting the SED reveals a temperature distribution that decreases towards its center, but heating from the ionizing source puts it above the general ISM temperature. We find that the bubble is filled by HII gas, ionized by an O8.5 star. Between the ionizing source and the IR dark clump G18.93/m we find a layered structure, from ionized to atomic to molecular hydrogen, revealing a PDR. Furthermore, we identify an additional velocity component within the bubble's 8mu emission rim at the edge of the infrared dark cloud and speculate that it might be shock induced by the expanding HII region. While the elevated temperature allows for the build-up of larger fragments, and the shock induced velocity component may lead to additional turbulent support, we do not find conclusive evidence that the massive clump G18.93/m is prone to collapse because of the expanding HII region.

preprint2011arXiv

Age spread in W3 Main: LBT/LUCI near-infrared spectroscopy of the massive stellar content

We present near-infrared multi-object spectroscopy and JHKs imaging of the massive stellar content of the Galactic star-forming region W3 Main, obtained with LUCI at the Large Binocular Telescope. We confirm 15 OB stars in W3 Main and derive spectral types between O5V and B4V from their absorption line spectra. Three massive Young Stellar Objects are identified by their emission line spectra and near-infrared excess. The color-color diagram of the detected sources allows a detailed investigation of the slope of the near-infrared extinction law towards W3 Main. Analysis of the Hertzsprung Russell diagram suggests that the Nishiyama extinction law fits the stellar population of W3 Main best (E(J-H)/E(H-Ks) = 1.76 and R_(Ks) = 1.44). From our spectrophotometric analysis of the massive stars and the nature of their surrounding HII regions we derive the evolutionary sequence of W3 Main and we find evidence of an age spread of at least 2-3 Myr. While the most massive star (IRS2) is already evolved, indications for high-mass pre-main-sequence evolution is found for another star (IRS N1), deeply embedded in an ultra compact HII region, in line with the different evolutionary phases observed in the corresponding HII regions. We derive a stellar mass of W3 Main of (4 +- 1) 10^3 Msun, by extrapolating from the number of OB stars using a Kroupa IMF and correcting for our spectroscopic incompleteness. We have detected the photospheres of OB stars from the more evolved diffuse HII region to the much younger UCHII regions, suggesting that these stars have finished their formation and cleared away their circumstellar disks very fast. Only in the hyper-compact HII region (IRS5), the early type stars seem to be still surrounded by circumstellar material.

preprint2011arXiv

Dissecting a hot molecular core: The case of G31.41+0.31

We made a detailed observational analysis of a well known hot molecular core lying in the high-mass star-forming region G31.41+0.31. This core is believed to contain deeply embedded massive stars and presents a velocity gradient that has been interpreted either as rotation or as expansion, depending on the authors. Our aim was to shed light on this question and possibly prepare the ground for higher resolution ALMA observations which could directly detect circumstellar disks around the embedded massive stars. Observations at sub-arcsecond resolution were performed with the Submillimeter Array in methyl cyanide, a typical hot molecular core tracer, and 12CO and 13CO, well known outflow tracers. We also obtained sensitive continuum maps at 1.3 mm. Our findings confirm the existence of a sharp velocity gradient across the core, but cannot confirm the existence of a bipolar outflow perpendicular to it. The improved angular resolution and sampling of the uv plane allow us to attain higher quality channel maps of the CH3CN lines with respect to previous studies and thus significantly improve our knowledge of the structure and kinematics of the hot molecular core. While no conclusive argument can rule out any of the two interpretations (rotation or expansion) proposed to explain the velocity gradient observed in the core, in our opinion the observational evidence collected so far indicates the rotating toroid as the most likely scenario. The outflow hypothesis appears less plausible, because the dynamical time scale is too short compared to that needed to form species such as CH3CN, and the mass loss and momentum rates estimated from our measurements appear too high.

preprint2011arXiv

Galactic structure based on the ATLASGAL 870mum survey

The ATLASGAL 870mum continuum survey conducted with the APEX telescope is the first survey covering the whole inner Galactic plane (60deg>l>-60deg & b<1.5deg) in submm continuum emission tracing the cold dust of dense and young star-forming regions. Here, we present the overall distribution of sources within our Galactic disk. The submm continuum emission is confined to a narrow range around the galactic plane, but shifted on average by ~0.07deg below the plane. Source number counts show strong enhancements toward the Galactic center, the spiral arms and toward prominent star-forming regions. Comparing the distribution of ATLASGAL dust continuum emission to that of young intermediate- to high-mass young stellar objects (YSOs) derived from Spitzer data, we find similarities as well as differences. In particular, the distribution of submm dust continuum emission is significantly more confined to the plane than the YSO distribution (FWHM of 0.7 and 1.1deg, corresponding to mean physical scale heights of approximately 46 and 80pc, respectively). While this difference may partly be caused by the large extinction from the dense submm cores, gradual dispersal of stellar distributions after their birth could also contribute to this effect. Compared to other tracers of Galactic structure, the ATLASGAL data are strongly confined to a narrow latitude strip around the Galactic plane.

preprint2011arXiv

High-mass star formation at high luminosities: W31 at >10^6 L_sun

Context: High-mass star formation has been a very active field over the last decade, however, most studies targeted regions of luminosities between 10^4 and 10^5 L_sun. Methods: We selected the W31 star-forming complex with a total luminosity of ~6x10^6 L_sun for a multi-wavelength spectral line and continuum study covering wavelengths from the near- and mid-infrared via (sub)mm wavelength observations to radio data in the cm regime. Results: While the overall structure of the multi-wavelength continuum data resembles each other well, there are several intriguing differences. The 24mum emission stemming largely from small dust grains follows tightly the spatial structure of the cm emission tracing the ionized free-free emission. Hence warm dust resides in regions that are spatially associated with the ionized hot gas (~10^4 K) of the HII regions. Furthermore, we find several evolutionary stages within the same complexes, ranging from infrared-observable clusters, via deeply embedded regions associated with active star formation traced by 24\,$μ$m and cm emission, to at least one high-mass gas clump devoid of any such signature. The 13CO(2-1) and C18O(2-1) spectral line observations reveal a large kinematic breadth in the entire region with a total velocity range of approximately 90 km/s. While the average virial mass ratio for W31 is close to unity, the line width analysis indicates large-scale evolutionary differences between the southern and northern sub-regions (G10.2-0.3 and G10.3-0.1) of the whole W31 complex. The clump mass function - tracing cluster scales and not scales of individual stars - derived from the 875mum continuum data has a slope of 1.5+-0.3, consistent with previous cloud mass functions.

preprint2011arXiv

Looking for outflow and infall signatures in high mass star forming regions

(Context) Many physical parameters change with time in star forming regions. Here we attempt to correlate changes in infall and outflow motions in high mass star forming regions with evolutionary stage using JCMT observations. (Aims) From a sample of 45 high mass star forming regions in three phases of evolution, we investigate the presence of established infall and outflow tracers to determine whether there are any trends attributable to the age of the source. (Methods) We obtained JCMT observations of HCO+/H13CO+ J=4-3 to trace large scale infall, and SiO J=8-7 to trace recent outflow activity. We compare the infall and outflow detections to the evolutionary stage of the host source (high mass protostellar objects, hypercompact HII regions and ultracompact HII regions). We also note that the integrated intensity of SiO varies with the full width at half maximum of the H13CO+. (Results) We find a surprising lack of SiO detections in the middle stage (Hypercompact HII regions), which may be due to an observational bias. When SiO is detected, we find that the integrated intensity of the line increases with evolutionary stage. We also note that all of the sources with infall signatures onto Ultracompact HII regions have corresponding outflow signatures as well.

preprint2011arXiv

Mass reservoirs surrounding massive infrared dark clouds: A view by near-infrared dust extinction

Context: Infrared Dark Clouds (IRDCs) harbor progenitors of high-mass stars. Little is known of the parental molecular clouds of the IRDCs. Aims: We demonstrate the feasibility of the near-infrared (NIR) dust extinction mapping in tracing the parental molecular clouds of IRDCs at the distances of D = 2.5 - 8 kpc. Methods: We derive NIR extinction maps for 10 prominent IRDC complexes using a color-excess mapping technique and NIR data from the UKIDSS/Galactic Plane Survey. We compare the resulting maps to the 13CO emission line data, to the 8 μm dust opacity data, and to the millimeter dust emission data. We derive distances for the clouds by comparing the observed NIR source densities to the Besancon stellar distribution model and compare them to the kinematic distance estimates. Results: The NIR extinction maps provide a view to the IRDC complexes over the dynamical range of Av = 2 - 40 mag, in spatial resolution of 30". The NIR extinction data correlate well with the 13CO data and probe a similar gas component, but also extend to higher column densities. The NIR data reveal a wealth of extended structures surrounding the dense gas traced by the 8 μm shadowing features and sub-mm dust emission, showing that the clouds contain typically > 10 times more mass than traced by those tracers. The IRDC complexes of our sample contain relatively high amount of high-column density material, and their cumulative column density distributions resemble active nearby star-forming clouds like Orion rather than less active clouds like California. Conclusions: NIR dust extinction data provide a new powerful tool to probe the mass distribution of the parental molecular clouds of IRDCs up to the distances of D = 8 kpc. This encourages for deeper NIR observations of IRDCs, because the sensitivity and resolution of the data can be directly enhanced with dedicated observations.

preprint2011arXiv

Molecular outflows and hot molecular core in G24.78+0.08 at sub-arcsecond angular resolution

Context. This study is part of a large project to study the physics of accretion and molecular outflows towards a selected sample of high-mass star-forming regions that show evidence of infall and rotation from previous studies. Aims. We wish to make a thorough study at high-angular resolution of the structure and kinematics of the HMCs and corresponding molecular outflows in the high-mass star-forming region G24.78+0.08. Methods. We carried out SMA and IRAM PdBI observations at 1.3 and 1.4 mm, respectively, of dust and of typical high-density and molecular outflow tracers with resolutions of <1". Complementary IRAM 30-m 12CO and 13CO observations were carried out to recover the short spacing information of the molecular outflows. Results. The millimeter continuum emission towards cores G24 A1 and A2 has been resolved into 3 and 2 cores, respectively, and named A1, A1b, A1c, A2, and A2b. All these cores are aligned in a southeast-northwest direction coincident with that of the molecular outflows detected in the region, which suggests a preferential direction for star formation in this region. The masses of the cores range from 7 to 22 Msun, and the rotational temperatures from 128 to 180 K. The high-density tracers have revealed the existence of 2 velocity components towards A1, one of them peaks close to the position of the millimeter continuum peak and of the HC HII region, and is associated with the velocity gradient seen in CH3CN towards this core, while the other one peaks southwest of core A1 and is not associated with any millimeter continuum emission peak. The position-velocity plots along outflow A and the 13CO averaged blueshifted and redshifted emission indicate that this outflow is driven by core A2.

preprint2011arXiv

Probing the evolution of molecular cloud structure II: From chaos to confinement

We present an analysis of the large-scale molecular cloud structure and of the stability of clumpy structures in nearby molecular clouds. In our recent work, we identified a structural transition in molecular clouds by studying the probability distributions of gas column densities in them. In this paper, we further examine the nature of this transition. The transition takes place at the visual extinction of A_V^tail = 2-4 mag, or equivalently, at Σ^tail = 40-80 Ms pc^{-2}. The clumps identified above this limit have wide ranges of masses and sizes, but a remarkably constant mean volume density of n = 10^3 cm^{-3}. This is 5-10 times larger than the density of the medium surrounding the clumps. By examining the stability of the clumps, we show that they are gravitationally unbound entities, and that the external pressure from the parental molecular cloud is a significant source of confining pressure for them. Then, the structural transition at A_V^tail may be linked to a transition between this population and the surrounding medium. The star formation rates in the clouds correlate strongly with the total mass in the clumps, i.e, with the mass above A_V^tail, dropping abruptly below that threshold. These results imply that the formation of pressure confined clumps introduces a prerequisite for star formation. Furthermore, they give a physically motivated explanation for the recently reported relation between the star formation rates and the amount of dense material in molecular clouds. Likewise, they give rise to a natural threshold for star formation at A_V^tail.

preprint2011arXiv

The Coalsack near and far

Context: The large Coalsack dark cloud is one of the rare molecular clouds without clear signs of star formation. Aims: We investigate the dynamical properties of the gas within the Coalsack. Methods: The two highest extinction regions were mapped with the APEX telescope in 13CO(2-1) comprising a region of ~1 square degree. Results: In addition to the well-known, nearby gas component around -4km/s, we identified additional molecular gas components - in particular a second extended molecular cloud at a velocity of ~-30km/s and an estimated distance of ~3.1kpc - that dominate the column density and visual extinction distributions in the northeastern part of the Coalsack. Although comprising ~2600Msun, the mass of this distant cloud is distributed over an extent of ~73pc, much larger than typical high-mass infrared dark clouds. Its filamentary structure is consistent with a compressible gaseous self-gravitating cylinder, and its low mass per length indicates that it may be stable against gravitational collapse. We find barely any mid-infrared emission in archival MSX data, which is indicative of almost no star-formation activity in the near and far cloud complexes. The nearby clouds have narrow, almost thermal velocity dispersions with median values between 0.2 and 0.4km/s, which is also consistent with low star-formation activity. Only Tapia's Globule 2 exhibits a velocity dispersion increase toward the extinction peak and peak-velocity gradients over the core, which is indicative of a state of elevated dynamical properties. Conclusions: The Coalsack is not one single coherent structure, but consists of several cloud complexes nearby as well as at several kpc distance. All studied clouds appear as starless low-turbulence regions that may not even collapse in the future. Only one globule exhibits more dynamical signatures and is a good candidate for present/future star formation.

preprint2011arXiv

The onset of high-mass star formation in the direct vicinity of the galactic mini-starburst W43

Context: The earliest stages of high-mass star formation are still poorly characterized. It is also unknown whether the initial conditions vary with environment. Aims: We want to investigate the youngest massive gas clumps in the environment of extremely active star formation. Methods: We selected the IRDC18454 complex associated with the W43 Galactic mini-starburst, and observed it in the continuum emission between 70mu and 1.2mm, and in spectral line emission of N2H+ and 13CO with single-dish instruments and an interferometer. Results: The multi-wavelength continuum study allows us to identify clumps that are infrared dark even at 70mu and hence the best candidates to be genuine high-mass starless gas clumps. The spectral energy distributions reveal elevated temperatures and luminosities compared to more quiescent environments. Furthermore, we identify a temperature gradient from the W43 mini-starburst toward the starless clumps. We discuss whether the radiation impact of the nearby mini-starburst changes the fragmentation properties of the gas clumps and by that maybe favors more high-mass star formation in such an environment. The spectral line data reveal two different velocity components of the gas at 100 and 50km/s. While chance projection is a possibility to explain these components, the projected associations of the emission sources as well as the prominent location at the Galactic bar -- spiral arm interface also allow the possibility that these two components may be spatially associated and even interacting. Conclusions: High-mass starless gas clumps can exist in the close environment of very active star formation without being destroyed. The impact of the active star formation sites may even allow for more high-mass stars to form in these 2nd generation gas clumps.

preprint2010arXiv

Chemistry in Infrared Dark Clouds

Massive stars play an important role in shaping the structure of galaxies. Infrared dark clouds (IRDCs), with their low temperatures and high densities, have been identified as the potential birthplaces of massive stars. In order to understand the formation processes of massive stars the physical and chemical conditions in infrared dark clouds have to be characterized. The goal of this paper is to investigate the chemical composition of a sample of southern infrared dark clouds. One important aspect of the observations is to check, if the molecular abuncances in IRDCs are similar to the low-mass pre-stellar cores, or whether they show signatures of more evolved evolutionary stages. We performed observations toward 15 IRDCs in the frequency range between 86 and 93 GHz using the 22-m Mopra radio telescope. We detect HNC, HCO$^+$ and HNC emission in all clouds and N$_2$H$^+$ in all IRDCs except one. In some clouds we detect SiO emission. Complicated shapes of the HCO$^+$ emission line profile are found in all IRDCs. Both signatures indicates the presence of infall and outflow motions and beginning of star formation activity, at least in some parts of the IRDCs. Where possible, we calculate molecular abundances and make a comparison with previously obtained values for low-mass pre-stellar cores and high-mass protostellar objects (HMPOs). We show a tendency for IRDCs to have molecular abundances similar to low-mass pre-stellar cores rather than to HMPOs abundances on the scale of our single-dish observations.

preprint2010arXiv

Different Evolutionary Stages in the Massive Star Forming Region S255 Complex

To understand evolutionary and environmental effects during the formation of high-mass stars, we observed three regions of massive star formation at different evolutionary stages that reside in the same natal molecular cloud. Methods. The three regions S255IR, S255N and S255S were observed at 1.3 mm with the Submillimeter Array (SMA) and followup short spacing information was obtained with the IRAM 30m telescope. Near infrared (NIR) H + K-band spectra and continuum observations were taken for S255IR with VLT-SINFONI to study the different stellar populations in this region. The combination of millimeter (mm) and near infrared data allow us to characterize different stellar populations within the young forming cluster in detail. While we find multiple mm continuum sources toward all regions, their outflow, disk and chemical properties vary considerably. The most evolved source S255IR exhibits a collimated bipolar outflow visible in CO and H2 emission, the outflows from the youngest region S255S are still small and rather confined in the regions of the mm continuum peaks. Also the chemistry toward S255IR is most evolved exhibiting strong emission from complex molecules, while much fewer molecular lines are detected in S255N, and in S255S we detect only CO isotopologues and SO lines. Also, rotational structures are found toward S255N and S255IR. Furthermore, a comparison of the NIR SINFONI and mm data from S255IR clearly reveal two different (proto) stellar populations with an estimated age difference of approximately 1 Myr. A multi-wavelength spectroscopy and mapping study reveals different evolutionary phases of the star formation regions. We propose the triggered outside-in collapse star formation scenario for the bigger picture and the fragmentation scenario for S255IR.

preprint2010arXiv

Disk and outflow signatures in Orion-KL: The power of high-resolution thermal infrared spectroscopy

We used the CRIRES spectrograph on the VLT to study the ro-vibrational 12CO/13CO, the Pfund beta and H2 emission between 4.59 and 4.72mu wavelengths toward the BN object, the disk candidate source n, and a proposed dust density enhancement IRC3. We detected CO absorption and emission features toward all three targets. Toward the BN object, the data partly confirm the results obtained more than 25 years ago by Scoville et al., however, we also identify several new features. While the blue-shifted absorption is likely due to outflowing gas, toward the BN object we detect CO in emission extending in diameter to ~3300AU. Although at the observational spectral resolution limit, the 13CO line width of that feature increases with energy levels, consistent with a disk origin. If one attributes the extended CO emission also to a disk origin, its extent is consistent with other massive disk candidates in the literature. For source n, we also find the blue-shifted CO absorption likely from an outflow. However, it also exhibits a narrower range of redshifted CO absorption and adjacent weak CO emission, consistent with infalling motions. We do not spatially resolve the emission for source n. For both sources we conduct a Boltzmann analysis of the 13CO absorption features and find temperatures between 100 and 160K, and H2 column densities of the order a few times 10^23cm^-2. The observational signatures from IRC3 are very different with only weak absorption against a much weaker continuum source. However, the CO emission is extended and shows wedge-like position velocity signatures consistent with jet-entrainment of molecular gas, potentially associated with the Orion-KL outflow system. We also present and discuss the Pfund beta and H2 emission in the region.

preprint2010arXiv

From high-mass starless cores to high-mass protostellar objects

Aims: Our aim is to understand the evolutionary sequence of high-mass star formation from the earliest evolutionary stage of high-mass starless cores, via high-mass cores with embedded low- to intermediate-mass objects, to finally high-mass protostellar objects. Methods: Herschel far-infrared PACS and SPIRE observations are combined with existing data at longer and shorter wavelengths to characterize the spectral and physical evolution of massive star-forming regions. Results: The new Herschel images spectacularly show the evolution of the youngest and cold high-mass star-forming regions from mid-infrared shadows on the Wien-side of the spectral energy distribution (SED), via structures almost lost in the background emission around 100mum, to strong emission sources at the Rayleigh-Jeans tail. Fits of the SEDs for four exemplary regions covering evolutionary stages from high-mass starless cores to high-mass protostellar objects reveal that the youngest regions can be fitted by single-component black-bodies with temperatures on the order of 17K. More evolved regions show mid-infrared excess emission from an additional warmer component, which however barely contributes to the total luminosities for the youngest regions. Exceptionally low values of the ratio between bolometric and submm luminosity additionally support the youth of the infrared-dark sources. Conclusions: The Herschel observations reveal the spectral and physical properties of young high-mass star-forming regions in detail. The data clearly outline the evolutionary sequence in the images and SEDs. Future work on larger samples as well as incorporating full radiative transfer calculations will characterize the physical nature at the onset of massive star formation in even more depth.

preprint2010arXiv

Magnetic field structure in a high-mass outflow/disk system

To characterize the magnetic field structure of the outflow and core region within a prototypical high-mass star-forming region, we analyzed polarized CO(3-2) - for the first time observed with the Submillimeter Array -- as well as 880mum submm continuum emission from the high-mass outflow/disk system IRAS18089-1732. Both emission features with polarization degrees at a few percent level indicate that the magnetic field structure is largely aligned with the outflow/jet orientation from the small core scales to the larger outflow scales. Although quantitative estimates are crude, the analysis indicates that turbulent energy dominates over magnetic energy. The data also suggest a magnetic field strength increase from the lower-density envelope to the higher-density core.

preprint2010arXiv

Mopra line survey mapping of NGC6334I and I(N) at 3mm

A 5'x5' region encompassing NGC6334I and I(N) has been mapped at a wavelength of 3mm (from 83.5 to 115.5GHz) with the Mopra telescope at an angular resolution between 33 arcsec and 36 arcsec. This investigation has made use of the recently installed 3mm MMIC receiver and the Mopra Spectrometer (MOPS) with broadband capabilities permitting total coverage of the entire frequency range with just five different observations. In total, the spatial distribution of nineteen different molecules, ions and radicals, along with additional selected isotopologues have been studied. Whilst most species trace the sites of star formation, CH_3CN appears to be most closely associated with NGC6334I and I(N). Both CN and C_2H appear to be widespread, tracing gas that is not associated with active star formation. Both N_2H^+ and HC_3N closely resemble dust continuum emission, showing they are reliable tracers of dense material, as well as the youngest stages of high mass star formation. Hot (E_u/k>100K) thermal CH_3OH emission is preferentially found towards NGC6334I, contrasting with I(N), where only cold (E_u/k<22K) thermal CH_3OH emission is found.

preprint2010arXiv

Sequential Star Formation in RCW 34: A Spectroscopic Census of the Stellar Content of High-mass Star-forming Regions

We present VLT/SINFONI integral field spectroscopy of RCW 34 along with Spitzer/IRAC photometry of the surroundings. RCW 34 consists of three different regions. A large bubble has been detected on the IRAC images in which a cluster of intermediate- and low-mass class II objects is found. At the northern edge of this bubble, an HII region is located, ionized by 3 OB stars. Intermediate mass stars (2 - 3 Msun) are detected of G- and K- spectral type. These stars are still in the pre-main sequence (PMS) phase. North of the HII region, a photon-dominated region is present, marking the edge of a dense molecular cloud traced by H2 emission. Several class 0/I objects are associated with this cloud, indicating that star formation is still taking place. The distance to RCW 34 is revised to 2.5 +- 0.2 kpc and an age estimate of 2 - 1 Myrs is derived from the properties of the PMS stars inside the HII region. The most likely scenario for the formation of the three regions is that star formation propagates from South to North. First the bubble is formed, produced by intermediate- and low-mass stars only, after that, the HII region is formed from a dense core at the edge of the molecular cloud, resulting in the expansion as a champagne flow. More recently, star formation occurred in the rest of the molecular cloud. Two different formation scenarios are possible: (a) The bubble with the cluster of low- and intermediate mass stars triggered the formation of the O star at the edge of the molecular cloud which in turn induces the current star-formation in the molecular cloud. (b) An external triggering is responsible for the star-formation propagating from South to North. [abridged]

preprint2010arXiv

The structured environments of embedded star-forming cores. PACS and SPIRE mapping of the enigmatic outflow source UYSO 1

The intermediate-mass star-forming core UYSO 1 has previously been found to exhibit intriguing features. While deeply embedded and previously only identified by means of its (sub-)millimeter emission, it drives two powerful, dynamically young, molecular outflows. Although the process of star formation has obviously started, the chemical composition is still pristine. We present Herschel PACS and SPIRE continuum data of this presumably very young region. The now complete coverage of the spectral energy peak allows us to precisely constrain the elevated temperature of 26 - 28 K for the main bulge of gas associated with UYSO1, which is located at the interface between the hot HII region Sh 2-297 and the cold dark nebula LDN 1657A. Furthermore, the data identify cooler compact far-infrared sources of just a few solar masses, hidden in this neighbouring dark cloud.

preprint2009arXiv

ATLASGAL - The APEX Telescope Large Area Survey of the Galaxy at 870 microns

(Abridged) Studying continuum emission from interstellar dust is essential to locating and characterizing the highest density regions in the interstellar medium. In particular, the early stages of massive star formation remain poorly understood. Our goal is to produce a large-scale, systematic database of massive pre- and proto-stellar clumps in the Galaxy, to understand how and under what conditions star formation takes place. A well characterized sample of star-forming sites will deliver an evolutionary sequence and a mass function of high-mass, star-forming clumps. This systematic survey at submm wavelengths also represents a preparatory work for Herschel and ALMA. The APEX telescope is ideally located to observe the inner Milky Way. The Large APEX Bolometer Camera (LABOCA) is a 295-element bolometer array observing at 870 microns, with a beam size of 19". Taking advantage of its large field of view (11.4') and excellent sensitivity, we started an unbiased survey of the Galactic Plane, with a noise level of 50-70 mJy/beam: the APEX Telescope Large Area Survey of the Galaxy (ATLASGAL). As a first step, we covered 95 sq. deg. These data reveal 6000 compact sources brighter than 0.25 Jy, as well as extended structures, many of them filamentary. About two thirds of the compact sources have no bright infrared counterpart, and some of them are likely to correspond to the precursors of (high-mass) proto-stars or proto-clusters. Other compact sources harbor hot cores, compact HII regions or young embedded clusters. Assuming a typical distance of 5 kpc, most sources are clumps smaller than 1 pc with masses from a few 10 to a few 100 M_sun. In this introductory paper, we show preliminary results from these ongoing observations, and discuss the perspectives of the survey.

preprint2009arXiv

Digging into NGC 6334I(N): Multiwavelength Imaging of a Massive Protostellar Cluster

We present a high-resolution, multi-wavelength study of the massive protostellar cluster NGC 6334I(N) that combines new spectral line data from the Submillimeter Array (SMA) and VLA with a reanalysis of archival VLA continuum data, 2MASS and Spitzer images. As shown previously, the brightest 1.3 mm source SMA1 contains substructure at subarcsecond resolution, and we report the first detection of SMA1b at 3.6 cm along with a new spatial component at 7 mm (SMA1d). We find SMA1 (aggregate of sources a, b, c, and d) and SMA4 to be comprised of free-free and dust components, while SMA6 shows only dust emission. Our 1.5" resolution 1.3 mm molecular line images reveal substantial hot-core line emission toward SMA1 and to a lesser degree SMA2. We find CH3OH rotation temperatures of 165\pm 9 K and 145\pm 12 K for SMA1 and SMA2, respectively. We estimate a diameter of 1400 AU for the SMA1 hot core emission, encompassing both SMA1b and SMA1d, and speculate that these sources comprise a >800 AU separation binary that may explain the previously-suggested precession of the outflow emanating from the SMA1 region. The LSR velocities of SMA1, SMA2, and SMA4 all differ by 1-2 km/s. Outflow activity from SMA1, SMA2, SMA4, and SMA6 is observed in several molecules including SiO(5--4) and IRAC 4.5 micron emission; 24 micron emission from SMA4 is also detected. Eleven water maser groups are detected, eight of which coincide with SMA1, SMA2, SMA4, and SMA6. We also detect a total of 83 Class I CH3OH 44GHz maser spots which likely result from the combined activity of many outflows. Our observations paint the portrait of multiple young hot cores in a protocluster prior to the stage where its members become visible in the near-infrared.

preprint2009arXiv

Hot high-mass accretion disk candidates

To better understand the physical properties of accretion disks in high-mass star formation, we present a study of a 12 high-mass accretion disk candidates observed at high spatial resolution with the Australia Telescope Compact Array (ATCA) in the NH3 (4,4) and (5,5) lines. Almost all sources were detected in NH3, directly associated with CH3OH Class II maser emission. From the remaining eleven sources, six show clear signatures of rotation and/or infall motions. These signatures vary from velocity gradients perpendicular to the outflows, to infall signatures in absorption against ultracompact HII regions, to more spherical infall signatures in emission. Although our spatial resolution is ~1000AU, we do not find clear Keplerian signatures in any of the sources. Furthermore, we also do not find flattened structures. In contrast to this, in several of the sources with rotational signatures, the spatial structure is approximately spherical with sizes exceeding 10^4 AU, showing considerable clumpy sub-structure at even smaller scales. This implies that on average typical Keplerian accretion disks -- if they exist as expected -- should be confined to regions usually smaller than 1000AU. It is likely that these disks are fed by the larger-scale rotating envelope structure we observe here. Furthermore, we do detect 1.25cm continuum emission in most fields of view.

preprint2009arXiv

Probing the evolution of molecular cloud structure: From quiescence to birth

Aims: We derive the probability density functions (PDFs) of column density for a complete sample of prominent molecular cloud complexes closer than 200 pc. Methods: We derive near-infrared dust extinction maps for 23 molecular cloud complexes, using the "nicest" colour excess mapping technique and data from the 2MASS archive. The extinction maps are then used to examine the column density PDFs in the clouds. Results: The column density PDFs of most molecular clouds are well-fitted by log-normal functions at low column densities (0.5 mag < A_v < 3-5 mag). However, at higher column densities prominent, power-law-like wings are common. In particular, we identify a trend among the PDFs: active star-forming clouds always have prominent non-log-normal wings. In contrast, clouds without active star formation resemble log-normals over the whole observed column density range, or show only low excess of higher column densities. This trend is also reflected in the cumulative PDFs, showing that the fraction of high column density material is significantly larger in star-forming clouds. These observations are in agreement with an evolutionary trend where turbulent motions are the main cloud-shaping mechanism for quiescent clouds, but the density enhancements induced by them quickly become dominated by gravity (and other mechanisms) which is strongly reflected by the shape of the column density PDFs. The dominant role of the turbulence is restricted to the very early stages of molecular cloud evolution, comparable to the onset of active star formation in the clouds.

preprint2009arXiv

Rotational Structure and Outflow in the Infrared Dark Cloud 18223-3

We examine an Infrared Dark Cloud (IRDC) at high spatial resolution as a means to study rotation, outflow, and infall at the onset of massive star formation. Submillimeter Array observations combined with IRAM 30 meter data in 12CO(2--1) reveal the outflow orientation in the IRDC 18223-3 region, and PdBI 3 mm observations confirm this orientation in other molecular species. The implication of the outflow's presence is that an accretion disk is feeding it, so using high density tracers such as C18O, N2H+, and CH3OH, we looked for indications of a velocity gradient perpendicular to the outflow direction. Surprisingly, this gradient turns out to be most apparent in CH3OH. The large size (28,000 AU) of the flattened rotating object detected indicates that this velocity gradient cannot be due solely to a disk, but rather from inward spiraling gas within which a Keplerian disk likely exists. From the outflow parameters, we derive properties of the source such as an outflow dynamical age of ~37,000 years, outflow mass of ~13 M_sun, and outflow energy of ~1.7 x 10^46 erg. While the outflow mass and energy are clearly consistent with a high-mass star forming region, the outflow dynamical age indicates a slightly more evolved evolutionary stage than previous spectral energy distribution (SED) modeling indicates. The calculated outflow properties reveal that this is truly a massive star in the making. We also present a model of the observed methanol velocity gradient. The rotational signatures can be modeled via rotationally infalling gas. These data present evidence for one of the youngest known outflow/infall/disk systems in massive star formation. A tentative evolutionary picture for massive disks is discussed.

preprint2009arXiv

Three intermediate-mass YSOs with different properties emerging from the same natal cloud in IRAS 00117+6412

We observed with the VLA, PdBI, and SMA the centimeter and millimeter continuum, N2H+(1-0), and CO(2-1) emission associated with a dusty cloud harboring a nascent cluster with intermediate-mass protostars. At centimeter wavelengths we found a strong source, tracing a UCHII region, at the eastern edge of the dusty cloud, with a shell-like structure, and with the near-infrared counterpart falling in the center of the shell. This is presumably the most massive source of the forming cluster. About 15'' to the west of the UCHII region and well embedded in the dusty cloud, we detected a strong millimeter source, MM1, associated with centimeter and near-infrared emission. MM1 seems to be driving a prominent high-velocity CO bipolar outflow, and is embedded in a ridge of dense gas traced by N2H+. We estimated that MM1 is an intermediate-mass source in the Class 0/I phase. About 15'' to the south of MM1, and still more deeply embedded in the dusty cloud, we detected a compact millimeter source, MM2, with neither centimeter nor near-infrared emission, but with water maser emission. MM2 is associated with a clump of N2H+, whose kinematics reveal a clear velocity gradient and additionally we found signposts of infall motions. MM2, being deeply embedded within the dusty cloud, with an associated water maser but no hints of CO outflow emission, is an intriguing object, presumably of intermediate mass. In conclusion, the UCHII region is found at the border of a dusty cloud which is currently undergoing active star formation. Two intermediate-mass protostars in the dusty cloud seem to have formed after the UCHII region and have different properties related to the outflow phenomenon.

preprint2007arXiv

The Protostar in the massive Infrared Dark Cloud IRDC18223-3

At the onset of high-mass star formation, accreting protostars are deeply embedded in massive cores made of gas and dust. Their spectral energy distribution is still dominated by the cold dust and rises steeply from near-to far-infrared wavelengths. The young massive star-forming region IRDC18223-3 is a prototypical Infrared-Dark-Cloud with a compact mm continuum core that shows no protostellar emission below 8mum. However, based on outflow tracers, early star formation activity was previously inferred for this region. Here, we present recent Spitzer observations from the MIPSGAL survey that identify the central protostellar object for the first time at 24 and 70mum. Combining the mid- to far-infrared data with previous mm continuum observations and the upper limits below 8mum, one can infer physical properties of the central source. At least two components with constant gas mass M and dust temperature T are necessary: one cold component (~15K and ~576M_sun) that contains most of the mass and luminosity, and one warmer component (>=51K and >=0.01M_sun) to explain the 24mum data. The integrated luminosity of ~177L_sun can be used to constrain additional parameters of the embedded protostar from the turbulent core accretion model for massive star formation. The data of IRDC18223-3 are consistent with a massive gas core harboring a low-mass protostellar seed of still less than half a solar mass with high accretion rates of the order 10^-4M_sun/yr. In the framework of this model, the embedded protostar is destined to become a massive star at the end of its formation processes.

preprint2005arXiv

Submillimeter Array 440$μ$m/690GHz line and continuum observations of Orion-KL

Submillimeter Array observations of Orion-KL at 1'' resolution in the 440mu/690GHz band reveal new insights about the continuum and line emission of the region. The 440mu continuum flux density measurement from source I allows us to differentiate among the various proposed physical models: Source I can be well modeled by a ``normal'' protostellar SED consisting of a proton-electron free-free emission component at low frequencies and a strong dust component in the submillimeter bands. Furthermore, we find that the protostellar object SMA1 is clearly distinct from the hot core. The non-detection of SMA1 at cm and infrared wavelengths suggests that it may be one of the youngest sources in the entire Orion-KL region. The molecular line maps show emission mainly from the sources I, SMA1 and the hot core peak position. An analysis of the CH$3CN(37_K-36_K) K-ladder (K=0...3) indicates a warm gas component of the order 600+-200K. In addition, we detect a large fraction (~58%) of unidentified lines and discuss the difficulties of line identifications at these frequencies.

preprint2002arXiv

IRAS 05358+3543: Multiple outflows at the earliest stages of massive star formation

We present a high-angular-resolution molecular line and millimeter continuum study of the massive star formation site IRAS 05358+3543. The most remarkable feature is a highly collimated (collimation factor ~10) and massive (>10 M_sun) bipolar outflow of 1 pc length, which is part of a quadrupolar outflow system. The three observed molecular outflows forming the IRAS 05358+3543 outflow system resemble, in structure and collimation, those typical of low-mass star-forming regions. They might therefore, just like low-mass outflows, be explained by shock entrainment models of jets. We estimate a mass accretion rate of 10^{-4) M_sun/yr, sufficient to overcome the radiative pressure of the central object and to build up a massive star, lending further support to the hypothesis that massive star formation occurs similarly to low-mass star formation, only with higher accretion rates and energetics.