Source author record

S. Bontemps

S. Bontemps 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

104works
2topics
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

104 published item(s)

preprint2022arXiv

A new phase of massive star formation? A luminous outflow cavity centred on an infrared quiet core

We present APEX, infrared and radio continuum observations of the G345.88-1.10 hub filament system which is a newly discovered star-forming cloud that hosts an unusually bright bipolar infrared nebulosity at its centre. At a distance of 2.26$^{+0.30}_{-0.21}$ kpc, G345.88-1.10 exhibits a network of parsec-long converging filaments. At the junction of these filaments lie four infrared-quiet fragments. The densest fragment (with M=210 M$_{\odot}$, R$_{\rm{eff}}=0.14$ pc) sits at the centre of a wide (opening angle of $\sim$ 90$\pm$15$^{o}$) bipolar nebulosity. $^{12}$CO(2-1) observations show that these infrared-bright nebulosities are spatially associated with a powerful molecular outflow from the central fragment. Negligible radio continuum and no H30$α$ emission is detected towards the cavities, seemingly excluding that ionising radiation drives the evolution of the cavities. Furthermore, radiative transfer simulations are unable to reproduce the observed combination of a low-luminosity ($\lesssim$ 500 L$_{\odot}$) central source and a surrounding high-luminosity ($\sim 4000$ L$_{\odot}$) mid-infrared-bright bipolar cavity. This suggests that radiative heating from a central protostar cannot be responsible for the illumination of the outflow cavities. To our knowledge, this is the first reported object of this type. The rarity of objects like G345.88-1.10 is likely related to a very short phase in the massive star and/or cluster formation process that was so far unidentified. We discuss whether mechanical energy deposition by one episode or successive episodes of powerful mass accretion in a collapsing hub might explain the observations. While promising in some aspects, a fully coherent scenario that explains the presence of a luminous bipolar cavity centred on an infrared-dark fragment remains elusive at this point.

preprint2022arXiv

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

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

preprint2022arXiv

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

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

preprint2021arXiv

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

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

preprint2020arXiv

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

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

preprint2020arXiv

Episodic accretion constrained by a rich cluster of outflows

The accretion history of protostars remains widely mysterious even though it represents one of the best ways to understand the protostellar collapse that leads to the formation of stars. Molecular outflows are here used to characterize the protostellar accretion phase in W43-MM1. The W43-MM1 protocluster host a sufficient number of protostars to statistically investigate molecular outflows in a single, homogeneous region. We used the CO(2-1) and SiO(5-4) line datacubes, taken as part of an ALMA mosaic with a 2000 AU resolution, to search for protostellar outflows, evaluate the influence that the environment has on these outflows' characteristics and put constraints on outflow variability in W43-MM1. We discovered a rich cluster of 46 outflow lobes, driven by 27 protostars with masses of 1-100 Msun. The complex environment inside which these outflow lobes develop has a definite influence on their length, limiting the validity of using outflows' dynamical timescales as a proxy of the ejection timescale in clouds with high dynamics and varying conditions. We performed a detailed study of Position-Velocity (PV) diagrams of outflows that revealed clear events of episodic ejection. The time variability of W43-MM1 outflows is a general trend and is more generally observed than in nearby, low- to intermediate-mass star-forming regions. The typical timescale found between two ejecta, about 500 yr, is consistent with that found in nearby protostars. If ejection episodicity reflects variability in the accretion process, either protostellar accretion is more variable or episodicity is easier to detect in high-mass star-forming regions than in nearby clouds. The timescale found between accretion events could be resulting from disk instabilities, associated with bursts of inflowing gas arising from the dynamical environment of high-mass star-forming cores.

preprint2020arXiv

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

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

preprint2020arXiv

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

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

preprint2020arXiv

Mass segregation and sequential star formation in NGC 2264 revealed by Herschel

The mass segregation of stellar clusters could be primordial rather than dynamical. Despite the abundance of studies of mass segregation for stellar clusters, those for stellar progenitors are still scarce, so the question on the origin and evolution of mass segregation is still open. Our goal is to characterize the structure of the NGC 2264 molecular cloud and compare the populations of clumps and young stellar objects (YSOs) in this region whose rich YSO population has shown evidence of sequential star formation. We separated the Herschel column density map of NGC 2264 in three subregions and compared their cloud power spectra using a multiscale segmentation technique. We identified in the whole NGC 2264 cloud a population of 256 clumps with typical sizes of ~0.1 pc and masses ranging from 0.08 Msun to 53 Msun. Although clumps have been detected all over the cloud, the central subregion of NGC 2264 concentrates most of the massive, bound clumps. The local surface density and the mass segregation ratio indeed indicate a strong degree of mass segregation for the 15 most massive clumps, with a median $Σ_6$ three time that of the whole clumps population and $Λ_{MSR}$ about 8. We showed that this cluster of massive clumps is forming within a high-density cloud ridge, itself formed and probably still fed by the high concentration of gas observed on larger scales in the central subregion. The time sequence obtained from the combined study of the clump and YSO populations in NGC 2264 suggests that the star formation started in the northern subregion, that it is now actively developing at the center and will soon start in the southern subregion. Taken together, the cloud structure and the clump and YSO populations in NGC 2264 argue for a dynamical scenario of star formation.

preprint2020arXiv

Physical properties of the ambient medium and of dense cores in the Perseus star-forming region derived from Herschel Gould Belt Survey observations

(Abridged) In this paper, we present analyses of images taken with the Herschel ESA satellite from 70mu to 500mu. We first constructed column density and dust temperature maps. Next, we identified compact cores in the maps, and characterize the cores using modified blackbody fits to their SEDs: we identified 684 starless cores, of which 199 are bound and potential prestellar cores, and 132 protostars. We also matched the Herschel-identified young stars with GAIA sources to model distance variations across the Perseus cloud. We measure a linear gradient function with right ascension and declination for the entire cloud. From the SED fits, mass and temperature of cores were derived. The core mass function can be modelled with a log-normal distribution that peaks at 0.82~$M_\sun$ suggesting a star formation efficiency of 0.30. The high-mass tail can be modelled with a power law of slope $\sim-2.32$, close to the Salpeter's value. We also identify the filamentary structure of Perseus, confirming that stars form preferentially in filaments. We find that the majority of filaments where star formation is ongoing are transcritical against their own internal gravity because their linear masses are below the critical limit of 16~$M_\sun$pc$^{-1}$ above which we expect filaments to collapse. We find a possible explanation for this result, showing that a filament with a linear mass as low as 8~$M_\sun$pc$^{-1}$ can be already unstable. We confirm a linear relation between star formation efficiency and slope of dust probability density function and a similar relation is also seen with the core formation efficiency. We derive a lifetime for the prestellar core phase of $1.69\pm0.52$~Myr for Perseus but different regions have a wide range in prestellar core fractions, hint that star-formation has started only recently in some clumps. We also derive a free-fall time for prestellar cores of 0.16~Myr.

preprint2020arXiv

Properties of the dense core population in Orion B as seen by the Herschel Gould Belt survey

We present a detailed study of the Orion B clouds (d~400 pc), imaged with the PACS/SPIRE cameras at 70-500 $μ$m by the Herschel Gould Belt survey (HGBS). We release new high-res. maps of column density and dust temperature. In the filamentary sub-regions NGC2023/2024, NGC2068/2071, and L1622, 1768 starless dense cores were identified, ~28-45% of which are self-gravitating prestellar cores. A total of 76 protostellar dense cores were also found. The typical lifetime of the prestellar cores was found to be $t_{\rm pre}=1.7_{-0.6}^{+0.8}$ Myr. The prestellar core mass function (CMF) peaks at ~0.5 $M_\odot$ and is consistent with a power law with log slope -1.27$\pm$0.24 at the high-mass end, compared to the Salpeter slope of -1.35. In this region, we confirm the existence of a transition in prestellar core formation efficiency (CFE) around a fiducial value A_V_bg~7 mag in background visual extinction, similar to the trend observed with Herschel in other clouds. This is not a sharp threshold, but a smooth transition between a regime with very low prestellar CFE at A_V_bg<5 and a regime with higher, roughly constant CFE at A_V_bg$\gtrsim$10. The total mass in the form of prestellar cores represents only ~20% of the dense molecular cloud gas at A_V_bg$\gtrsim$7 mag. About 60-80% of the prestellar cores are closely associated with filaments, and this fraction increases up to >90% when a more complete sample of filamentary structures is considered. Interestingly, the median separation between nearest core neighbors corresponds to the typical inner filament width of ~0.1 pc commonly observed in nearby molecular clouds. Analysis of the CMF observed as a function of background cloud column density shows that the most massive prestellar cores are spatially segregated in the highest column density areas, and suggests that both higher- and lower-mass prestellar cores may form in denser filaments.

preprint2020arXiv

Questioning the spatial origin of complex organic molecules in young protostars with the CALYPSO survey

Complex organic molecules (COMs) have been detected in a few Class 0 protostars but their origin is not well understood. Going beyond studies of individual objects, we want to investigate the origin of COMs in young protostars on a statistical basis. We use the CALYPSO survey performed with the IRAM PdBI to search for COMs at high angular resolution in a sample of 26 solar-type protostars, including 22 Class 0 and four Class I objects. Methanol is detected in 12 sources and tentatively in one source, which represents half of the sample. Eight sources (30%) have detections of at least three COMs. We find a strong chemical differentiation in multiple systems with five systems having one component with at least three COMs detected but the other component devoid of COM emission. The internal luminosity is found to be the source parameter impacting the most the COM chemical composition of the sources, while there is no obvious correlation between the detection of COM emission and that of a disk-like structure. A canonical hot-corino origin may explain the COM emission in four sources, an accretion-shock origin in two or possibly three sources, and an outflow origin in three sources. The CALYPSO sources with COM detections can be classified into three groups on the basis of the abundances of oxygen-bearing molecules, cyanides, and CHO-bearing molecules. These chemical groups correlate neither with the COM origin scenarii, nor with the evolutionary status of the sources if we take the ratio of envelope mass to internal luminosity as an evolutionary tracer. We find strong correlations between molecules that are a priori not related chemically (for instance methanol and methyl cyanide), implying that the existence of a correlation does not imply a chemical link. [abridged]

preprint2020arXiv

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

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

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.

preprint2016arXiv

Characterizing filaments in regions of high-mass star formation: High-resolution submilimeter imaging of the massive star-forming complex NGC 6334 with ArTéMiS

Herschel observations of nearby molecular clouds suggest that interstellar filaments and prestellar cores represent two fundamental steps in the star formation process. The observations support a picture of low-mass star formation according to which ~ 0.1 pc-wide filaments form first in the cold interstellar medium, probably as a result of large-scale compression of interstellar matter by supersonic turbulent flows, and then prestellar cores arise from gravitational fragmentation of the densest filaments. Whether this scenario also applies to regions of high-mass star formation is an open question, in part because Herschel data cannot resolve the inner width of filaments in the nearest regions of massive star formation. We used the bolometer camera ArTeMiS on the APEX telescope to map the central part of the NGC6334 complex at a factor of > 3 higher resolution than Herschel at 350 microns. Combining ArTeMiS data with Herschel data allowed us to study the structure of the main filament of the complex with a resolution of 8" or < 0.07 pc at d ~ 1.7 kpc. Our study confirms that this filament is a very dense, massive linear structure with a line mass ranging from ~ 500 Msun/pc to ~ 2000 Msun/pc over nearly 10 pc. It also demonstrates that its inner width remains as narrow as W ~ 0.15 +- 0.05 pc all along the filament length, within a factor of < 2 of the characteristic 0.1 pc value found with Herschel for lower-mass filaments in the Gould Belt. While it is not completely clear whether the NGC 6334 filament will form massive stars or not in the future, it is two to three orders of magnitude denser than the majority of filaments observed in Gould Belt clouds, and yet has a very similar inner width. This points to a common physical mechanism for setting the filament width and suggests that some important structural properties of nearby clouds also hold in high-mass star forming regions.

preprint2016arXiv

Dense gas in low-metallicity galaxies

Stars form out of the densest parts of molecular clouds. Far-IR emission can be used to estimate the Star Formation Rate (SFR) and high dipole moment molecules, typically HCN, trace the dense gas. A strong correlation exists between HCN and Far-IR emission, with the ratio being nearly constant, over a large range of physical scales. A few recent observations have found HCN to be weak with respect to the Far-IR and CO in subsolar metallicity (low-Z) objects. We present observations of the Local Group galaxies M33, IC10, and NGC6822 with the IRAM 30meter and NRO 45m telescopes, greatly improving the sample of low-Z galaxies observed. HCN, HCO$^+$, CS, C$_2$H, and HNC have been detected. Compared to solar metallicity galaxies, the Nitrogen-bearing species are weak (HCN, HNC) or not detected (CN, HNCO, N$_2$H$^+$) relative to Far-IR or CO emission. HCO$^+$ and C$_2$H emission is normal with respect to CO and Far-IR. While $^{13}$CO is the usual factor 10 weaker than $^{12}$CO, C$^{18}$O emission was not detected down to very low levels. Including earlier data, we find that the HCN/HCO$^+$ ratio varies with metallicity (O/H) and attribute this to the sharply decreasing Nitrogen abundance. The dense gas fraction, traced by the HCN/CO and HCO$^+$/CO ratios, follows the SFR but in the low-Z objects the HCO$^+$ is much easier to measure. Combined with larger and smaller scale measurements, the HCO$^+$ line appears to be an excellent tracer of dense gas and varies linearly with the SFR for both low and high metallicities.

preprint2016arXiv

First image of the L1157 molecular jet by the CALYPSO IRAM-PdBI survey

Fast jets are thought to be a crucial ingredient of star formation because they might extract angular momentum from the disk and thus allow mass accretion onto the star. However, it is unclear whether jets are ubiquitous, and likewise, their contribution to mass and angular momentum extraction during protostar formation remains an open question. Our aim is to investigate the ejection process in the low-mass Class 0 protostar L1157. This source is associated with a spectacular bipolar outflow, and the recent detection of high-velocity SiO suggests the occurrence of a jet. Observations of CO 2-1 and SiO 5-4 at 0.8" resolution were obtained with the IRAM Plateau de Bure Interferometer as part of the CALYPSO large program. The jet and outflow structure were fit with a precession model. We derived the column density of CO and SiO, as well as the jet mass-loss rate and mechanical luminosity. High-velocity CO and SiO emission resolve for the first time the first 200 au of the outflow-driving molecular jet. The jet is strongly asymmetric, with the blue lobe 0.65 times slower than the red lobe. This suggests that the large-scale asymmetry of the outflow is directly linked to the jet velocity and that the asymmetry in the launching mechanism has been at work for the past 1800 yr. Velocity asymmetries are common in T Tauri stars, which suggests that the jet formation mechanism from Class 0 to Class II stages might be similar. Our model simultaneously fits the inner jet and the clumpy 0.2 pc scale outflow by assuming that the jet precesses counter-clockwise on a cone inclined by 73 degree to the line of sight with an opening angle of 8 degree on a period of 1640 yr. The estimated jet mass flux and mechanical luminosity are 7.7e-7 Msun/yr, and 0.9 Lsun, indicating that the jet could extract at least 25% of the gravitational energy released by the forming star.

preprint2016arXiv

Globules and Pillars in Cygnus X I. Herschel Far-infrared imaging of the Cyg OB2 environment

The radiative feedback of massive stars on molecular clouds creates pillars, globules and other features at the interface between the HII region and molecular cloud. We present here Herschel observations between 70 and 500 micron of the immediate environment of the Cygnus OB2 association, performed within the HOBYS program. All structures were detected based on their appearance at 70 micron, and have been classified as pillars, globules, evaporating gasous globules (EGGs), proplyd-like objects, and condensations. From the 70 and 160 micron flux maps, we derive the local FUV field on the PDR surfaces. In parallel, we use a census of the O-stars to estimate the overall FUV-field, that is 10^3-10^4 G_0 close to the central OB cluster (within 10 pc) and decreases down to a few tens G_0, in a distance of 50 pc. From a SED fit to the four longest Herschel wavelengths, we determine column density and temperature maps and derive masses, volume densities and surface densities for these structures. We find that the morphological classification corresponds to distinct physical properties. Pillars and globules have the longest estimated photoevaporation lifetimes, a few 10^6 yr, while all other features should survive less than that. These lifetimes are consistent with that found in simulations of turbulent, UV-illuminated clouds. We propose a tentative evolutionary scheme in which pillars can evolve into globules, which in turn then evolve into EGGs, condensations and proplyd-like objects.

preprint2016arXiv

Probing the CO and methanol snow lines in young protostars. Results from the CALYPSO IRAM-PdBI survey

Context. "Snow lines", marking regions where abundant volatiles freeze out onto the surface of dust grains, play an important role for planet growth and bulk composition in protoplanetary disks. They can already be observed in the envelopes of the much younger, low-mass Class 0 protostars that are still in their early phase of heavy accretion. Aims. We aim at using the information on the sublimation regions of different kinds of ices to understand the chemistry of the envelope, its temperature and density structure, and the history of the accretion process. Methods. As part of the CALYPSO IRAM Large Program, we have obtained observations of C$^{18}$O, N$_2$H$^+$ and CH$_3$OH towards nearby Class 0 protostars with the IRAM Plateau de Bure interferometer at sub-arcsecond resolution. For four of these sources we have modeled the emission using a chemical code coupled with a radiative transfer module. Results. We observe an anti-correlation of C$^{18}$O and N$_2$H$^+$ in NGC 1333-IRAS4A, NGC 1333-IRAS4B, L1157, and L1448C, with N$_2$H$^+$ forming a ring around the centrally peaked C$^{18}$O emission due to N$_2$H$^+$ being chemically destroyed by CO. The emission regions of models and observations match for a CO binding energy of 1200 K, which is higher than the binding energy of pure CO ices ($\sim$855 K). Furthermore, we find very low CO abundances inside the snow lines in our sources, about an order of magnitude lower than the total CO abundance observed in the gas on large scales in molecular clouds before depletion sets in. Conclusions. The high CO binding energy may hint at CO being frozen out in a polar ice environment like amorphous water ice or in non-polar CO$_2$-rich ice. The low CO abundances are comparable to values found in protoplanetary disks, which may indicate an evolutionary scenario where these low values are already established in the protostellar phase. (Abbr. Version)

preprint2016arXiv

The Herschel-HIFI view of mid-IR quiet massive protostellar objects

We present Herschel/HIFI observations (WISH KP) of 14 water lines in a small sample of galactic massive protostellar objects: NGC6334I(N), DR21(OH), IRAS16272-4837, and IRAS05358+3543. We analyze the gas dynamics from the line profiles. Through modeling of the observations using RATRAN, we estimate outflow, infall, turbulent velocities, molecular abundances, and investigate any correlation with the evolutionary status of each source. The molecular line profiles exhibit a broad component coming from the shocks along the cavity walls associated with the protostars, and an infalling (or expansion for IRAS05358+3543) and passively heated envelope component, with highly supersonic turbulence likely increasing with the distance from the center. Accretion rates between 6.3 10^{-5} and 5.6 10^{-4} \msun yr^{-1} are derived from the infall observed in three of our sources. The outer water abundance is estimated to be at the typical value of a few 10^{-8} while the inner abundance varies from 1.7 10^{-6} to 1.4 10^{-4} with respect to H2 depending on the source. We confirm that regions of massive star formation are highly turbulent and that the turbulence likely increases in the envelope with the distance to the star. The inner abundances are lower than the expected 10^{-4} perhaps because our observed lines do not probe deep enough into the inner envelope, or because photodissociation through protostellar UV photons is more efficient than expected. We show that the higher the infall/expansion velocity in the protostellar envelope, the higher is the inner abundance, maybe indicating that larger infall/expansion velocities generate shocks that will sputter water from the ice mantles of dust grains in the inner region. High-velocity water must be formed in the gas-phase from shocked material.

preprint2016arXiv

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

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

preprint2015arXiv

A census of dense cores in the Aquila cloud complex: SPIRE/PACS observations from the Herschel Gould Belt survey

We present and discuss the results of the Herschel Gould Belt survey observations in a ~11 deg^2 area of the Aquila molecular cloud complex at d~260 pc, imaged with the SPIRE/PACS cameras from 70 to 500 micron. We identify a complete sample of starless dense cores and embedded protostars in this region, and analyze their global properties and spatial distributions. We find a total of 651 starless cores, ~60% of which are gravitationally bound prestellar cores, and they will likely form stars in the future. We also detect 58 protostellar cores. The core mass function (CMF) derived for the prestellar cores is very similar in shape to the stellar initial mass function (IMF), supporting the earlier view that there is a close physical link between the IMF and the CMF. The global shift in mass scale observed between the CMF and the IMF is consistent with a typical star formation efficiency of ~40%. By comparing the numbers of starless cores to the number of young stellar objects, we estimate that the lifetime of prestellar cores is ~1 Myr. We find a strong correlation between the spatial distribution of prestellar cores and the densest filaments. About 90% of the Herschel-identified prestellar cores are located above a background column density corresponding to A_V~7, and ~75% of them lie within filamentary structures with supercritical masses per unit length >~16 M_sun/pc. These findings support a picture wherein the cores making up the peak of the CMF (and probably responsible for the base of the IMF) result primarily from the gravitational fragmentation of marginally supercritical filaments. Given that filaments appear to dominate the mass budget of dense gas at A_V>7, our findings also suggest that the physics of prestellar core formation within filaments is responsible for a characteristic "efficiency" SFR/M_dense ~5+-2 x 10^-8 yr^-1 for the star formation process in dense gas.

preprint2015arXiv

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

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

preprint2015arXiv

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

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

preprint2015arXiv

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

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

preprint2015arXiv

First detection of CF+ towards a high-mass protostar

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

preprint2015arXiv

From forced collapse to H ii region expansion in Mon R2: Envelope density structure and age determination with Herschel

The surroundings of HII regions can have a profound influence on their development, morphology, and evolution. This paper explores the effect of the environment on H II regions in the MonR2 molecular cloud. We aim to investigate the density structure of envelopes surrounding HII regions and to determine their collapse and ionisation expansion ages. The Mon R2 molecular cloud is an ideal target since it hosts an H II region association. Column density and temperature images derived from Herschel data were used together to model the structure of HII bubbles and their surrounding envelopes. The resulting observational constraints were used to follow the development of the Mon R2 ionised regions with analytical calculations and numerical simulations. The four hot bubbles associated with H II regions are surrounded by dense, cold, and neutral gas envelopes. The radial density profiles are reminiscent of those of low-mass protostellar envelopes. The inner parts of envelopes of all four HII regions could be free-falling because they display shallow density profiles. As for their outer parts, the two compact HII regions show a density profile, which is typical of the equilibrium structure of an isothermal sphere. In contrast, the central UCHii region shows a steeper outer profile, that could be interpreted as material being forced to collapse. The size of the heated bubbles, the spectral type of the irradiating stars, and the mean initial neutral gas density are used to estimate the ionisation expansion time, texp, 0.1Myr,for the dense UCHII and compact HII regions and 0.35 Myr for the extended HII region. The envelope transition radii between the shallow and steeper density profiles are used to estimate the time elapsed since the formation of the first proto stellar embryo, Tinf : 1Myr, for the ultra-compact, 1.5 / 3Myr for the compact, and greater than 6Myr for the extended HII regions.

preprint2015arXiv

Herschel Observations of the W3 GMC (II): Clues to the Formation of Clusters of High-Mass Stars

The W3 GMC is a prime target for investigating the formation of high-mass stars and clusters. This second study of W3 within the HOBYS Key Program provides a comparative analysis of subfields within W3 to further constrain the processes leading to the observed structures and stellar population. Probability density functions (PDFs) and cumulative mass distributions (CMDs) were created from dust column density maps, quantified as extinction Av. The shape of the PDF, typically represented with a lognormal function at low Av "breaking" to a power-law tail at high Av, is influenced by various processes including turbulence and self-gravity. The breaks can also be identified, often more readily, in the CMDs. The PDF break from lognormal (Av(SF)= 6-10 mag) appears to shift to higher Av by stellar feedback, so that high-mass star-forming regions tend to have higher PDF breaks. A second break at Av > 50 mag traces structures formed or influenced by a dynamic process. Because such a process has been suggested to drive high-mass star formation in W3, this second break might then identify regions with potential for hosting high-mass stars/clusters. Stellar feedback appears to be a major mechanism driving the local evolution and state of regions within W3. A high initial star formation efficiency in a dense medium could result in a self-enhancing process, leading to more compression and favourable star-formation conditions (e.g., colliding flows), a richer stellar content, and massive stars. This scenario would be compatible with the "convergent constructive feedback" model introduced in our previous Herschel study.

preprint2015arXiv

SOFIA/FORCAST Observations of Warm Dust in S106: A Fragmented Environment

We present mid-IR (19 - 37 microns) imaging observations of S106 from SOFIA/FORCAST, complemented with IR observations from Spitzer/IRAC (3.6 - 8.0 microns), IRTF/MIRLIN (11.3 and 12.5 microns), and Herschel/PACS (70 and 160 microns). We use these observations, observations in the literature, and radiation transfer modeling to study the heating and composition of the warm (~ 100 K) dust in the region. The dust is heated radiatively by the source S106 IR, with little contributions from grain-electron collisions and Ly-alpha radiation. The dust luminosity is >~ (9.02 +/- 1.01) x 10^4 L_sun, consistent with heating by a mid- to late-type O star. We find a temperature gradient (~ 75 - 107 K) in the lobes, which is consistent with a dusty equatorial geometry around S106 IR. Furthermore, the SOFIA observations resolve several cool (~ 65 - 70 K) lanes and pockets of warmer (~ 75 - 90 K) dust in the ionization shadow, indicating that the environment is fragmented. We model the dust mass as a composition of amorphous silicates, amorphous carbon, big grains, very small grains, and PAHs. We present the relative abundances of each grain component for several locations in S106.

preprint2015arXiv

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

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

preprint2015arXiv

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

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

preprint2015arXiv

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

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

preprint2015arXiv

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

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

preprint2014arXiv

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

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

preprint2014arXiv

Class 0 Protostars in the Perseus Molecular Cloud: A Correlation Between the Youngest Protostars and the Dense Gas Distribution

We use PACS and SPIRE continuum data at 160 um, 250 um, 350 um, and 500 um from the Herschel Gould Belt Survey to sample seven clumps in Perseus: B1, B1-E, B5, IC348, L1448, L1455, and NGC1333. Additionally, we identify and characterize the embedded Class 0 protostars using detections of compact Herschel sources at 70 um as well as archival Spitzer catalogues and SCUBA 850 um photometric data. We identify 28 candidate Class 0 protostars, four of which are newly discovered sources not identified with Spitzer. We find that the star formation efficiency of clumps, as traced by Class 0 protostars, correlates strongly with the flatness of their respective column density distributions at high values. This correlation suggests that the fraction of high column density material in a clump reflects only its youngest protostellar population rather than its entire source population. We propose that feedback from either the formation or evolution of protostars changes the local density structure of clumps.

preprint2014arXiv

First results from the CALYPSO IRAM-PdBI survey. I. Kinematics of the inner envelope of NGC1333-IRAS2A

The structure and kinematics of Class 0 protostars on scales of a few hundred AU is poorly known. Recent observations have revealed the presence of Keplerian disks with a diameter of 150-180 AU in L1527-IRS and VLA1623A, but it is not clear if such disks are common in Class 0 protostars. Here we present high-angular-resolution observations of two methanol lines in NGC1333-IRAS2A. We argue that these lines probe the inner envelope, and we use them to study the kinematics of this region. Our observations suggest the presence of a marginal velocity gradient normal to the direction of the outflow. However, the position velocity diagrams along the gradient direction appear inconsistent with a Keplerian disk. Instead, we suggest that the emission originates from the infalling and perhaps slowly rotating envelope, around a central protostar of 0.1-0.2 M$_\odot$. If a disk is present, it is smaller than the disk of L1527-IRS, perhaps suggesting that NGC1333-IRAS2A is younger.

preprint2014arXiv

First results from the CALYPSO IRAM-PdBI survey. II. Resolving the hot corino in the Class 0 protostar NGC 1333-IRAS2A

We investigate the origin of complex organic molecules (COMs) in the gas phase around the low-mass Class~0 protostar NGC1333-IRAS2A, to determine if the COM emission lines trace an embedded disk, shocks from the protostellar jet, or the warm inner parts of the protostellar envelope. In the framework of the CALYPSO (Continuum And Lines in Young ProtoStellar Objects) IRAM Plateau de Bure survey, we obtained large bandwidth spectra at sub-arcsecond resolution towards NGC 1333-IRAS2A. We identify the emission lines towards the central protostar and perform Gaussian fits to constrain the size of the emitting region for each of these lines, tracing various physical conditions and scales. The emission of numerous COMs such as methanol, ethylene glycol, and methyl formate is spatially resolved by our observations. This allows us to measure, for the first time, the size of the COM emission inside the protostellar envelope, finding that it originates from a region of radius 40-100 AU, centered on the NGC 1333-IRAS2A protostellar object. Our analysis shows no preferential elongation of the COM emission along the jet axis, and therefore does not support the hypothesis that COM emission arises from shocked envelope material at the base of the jet. Down to similar sizes, the dust continuum emission is well reproduced with a single envelope model, and therefore does not favor the hypothesis that COM emission arises from the thermal sublimation of grains embedded in a circumstellar disk. Finally, the typical scale $\sim$60 AU observed for COM emission is consistent with the size of the inner envelope where $T_{\rm{dust}} > 100$ K is expected. Our data therefore strongly suggest that the COM emission traces the hot corino in IRAS2A, i.e., the warm inner envelope material where the icy mantles of dust grains evaporate because they are passively heated by the central protostellar object.

preprint2014arXiv

Fragmentation of massive dense cores down to ~1000 AU: Relation between fragmentation and density structure

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

preprint2014arXiv

Herschel's view of the large-scale structure in the Chamaeleon dark clouds

The Chamaeleon molecular cloud complex is one of the nearest star-forming sites encompassing three molecular clouds with a different star-formation history, from quiescent (Cha III) to actively forming stars (Cha II), and reaching the end of star-formation (Cha I). To charactize its large-scale structure, we derived column density and temperature maps using PACS and SPIRE observations from the Herschel Gould Belt Survey, and applied several tools, such as filament tracing, power-spectra, Δ-variance, and probability distribution functions of column density (PDFs), to derive physical properties. The column density maps reveal a different morphological appearance for the three clouds, with a ridge-like structure for Cha I, a clump-dominated regime for Cha II, and an intricate filamentary network for Cha III. The filament width is measured to be around 0.12\pm0.04 pc in the three clouds, and the filaments found to be gravitationally unstable in Cha I and II, but mostly subcritical in Cha III. Faint filaments (striations) are prominent in Cha I showing a preferred alignment with the large-scale magnetic field. The PDFs of all regions show a lognormal distribution at low column densities. For higher densities, the PDF of Cha I shows a turnover indicative of an extended higher density component, culminating with a power-law tail. Cha II shows a power-law tail with a slope characteristic of gravity. The PDF of Cha III can be best fit by a single lognormal. The turbulence properties of the three regions are found to be similar, pointing towards a scenario where the clouds are impacted by large-scale processes. The magnetic field could possibly play an important role for the star-formation efficiency in the Chamaeleon clouds if proven that it can effectively channel material on Cha I, and possibly Cha II, but probably less efficiently on the quiescent Cha III cloud.

preprint2014arXiv

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

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

preprint2014arXiv

Properties of Starless and Prestellar Cores in Taurus Revealed by Herschel SPIRE/PACS Imaging

The density and temperature structures of dense cores in the L1495 cloud of the Taurus star-forming region are investigated using Herschel SPIRE and PACS images in the 70 $μ$m, 160 $μ$m, 250 $μ$m, 350 $μ$m and 500 $μ$m continuum bands. A sample consisting of 20 cores, selected using spectral and spatial criteria, is analysed using a new maximum likelihood technique, COREFIT, which takes full account of the instrumental point spread functions. We obtain central dust temperatures, $T_0$, in the range 6-12 K and find that, in the majority of cases, the radial density falloff at large radial distances is consistent with the $r^{-2}$ variation expected for Bonnor-Ebert spheres. Two of our cores exhibit a significantly steeper falloff, however, and since both appear to be gravitationally unstable, such behaviour may have implications for collapse models. We find a strong negative correlation between $T_0$ and peak column density, as expected if the dust is heated predominantly by the interstellar radiation field. At the temperatures we estimate for the core centres, carbon-bearing molecules freeze out as ice mantles on dust grains, and this behaviour is supported here by the lack of correspondence between our estimated core locations and the previously-published positions of H$^{13}$CO$^+$ peaks. On this basis, our observations suggest a sublimation-zone radius typically $\sim 10^4$ AU. Comparison with previously-published N$_2$H$^+$ data at 8400 AU resolution, however, shows no evidence for N$_2$H$^+$ depletion at that resolution.

preprint2014arXiv

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

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

preprint2014arXiv

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

The IRAM M33 CO(2-1) Survey - A complete census of the molecular gas out to 7 kpc

In order to study the ISM and the interplay between the atomic and molecular components in a low-metallicity environment, we present a complete high angular and spectral resolution map and data cube of the 12CO(2-1) emission from the Local Group galaxy M33. Its metallicity is roughly half-solar, such that we can compare its ISM with that of the Milky Way with the main changes being the metallicity and the gas mass fraction. The data have a 12" angular resolution (50pc) with a spectral resolution of 2.6 km/s and a mean noise level of 20 mK per channel in antenna temperature. A radial cut along the major axis was also observed in the 12CO(1-0) line. The CO data cube and integrated intensity map are optimal when using HI data to define the baseline window and the velocities over which the CO emission is integrated. Great care was taken when building these maps, testing different windowing and baseline options and investigating the effect of error beam pickup. The total CO(2-1) luminosity is 2.8e7 K km/s pc2, following the spiral arms in the inner disk. There is no clear variation in the CO(2-1/1-0) intensity ratio with radius and the average value is roughly 0.8. The total molecular gas mass is estimated, using a N(H2)/Ico(1-0)=4e20 cm-2/(K km/s) conversion factor, to be 3.1e8 Msol. The CO spectra in the cube were shifted to zero velocity by subtracting the velocity of the HI peak from the CO spectra. Hence, the velocity dispersion between the atomic and molecular components is extremely low, independently justifying the use of the HI line in building our maps. Stacking the spectra in concentric rings shows that the CO linewidth and possibly the CO-HI velocity dispersion decrease in the outer disk. Using the CO(2-1) emission to trace the molecular gas, the probability distribution function of the H2 column density shows an excess at high column density above a log normal distribution.

preprint2013arXiv

CO outflows from high-mass Class 0 protostars in Cygnus-X

As natural consequences of the accretion process, outflows are one of the few (indirect) tracers of accretion. We used CO(2-1) PdBI observations towards 6 MDCs in Cygnus-X, containing 9 high-mass cores, to investigate what the accretion process and origin of the material feeding the precursors of high-mass stars are. We compared our sample to low-mass objects from the literature and developed simple evolutionary models to reproduce the observables. We find that 8/9 high-mass cores drive clear individual outflows as true equivalents of Class 0 protostars in the high-mass regime. The remaining core has only a tentative outflow detection. It could be amongst the first examples of a true individual high-mass prestellar core. We find that the momentum flux of high-mass objects has a linear relation to the envelope mass, as a scale-up of the relations found for low-mass protostars. This suggests a fundamental proportionality between accretion rates and mass reservoir, suggesting identical collapse timescales for all masses. We conclude that if the pre-collapse evolution is quasi-static, the fragmentation scale (which is similar for all masses) would limit the size of the initial mass reservoirs, leading to shorter free-fall times for higher mass stars. However, as we find identical collapse timescales for all masses, a significant turbulent/magnetic support is needed to slow down the collapse of the more massive envelopes in a quasi-static view. With this support still to be discovered, and with indications of large dynamics in pre-collapse gas for high-mass star formation, we propose that such an identical collapse timescale implies that the initial densities, which should set the duration of the collapse, should be similar for all masses. This suggests that the mass that incorporates massive stars has to have been accreted in a dynamical way from larger scales than those of low-mass stars.

preprint2013arXiv

First results from the Herschel Gould Belt Survey in Taurus

The whole of the Taurus region (a total area of 52 sq. deg.) has been observed by the Herschel SPIRE and PACS instruments at wavelengths of 70, 160, 250, 350 and 500 μm as part of the Herschel Gould Belt Survey. In this paper we present the first results from the part of the Taurus region that includes the Barnard 18 and L1536 clouds. A new source-finding routine, the Cardiff Source-finding AlgoRithm (CSAR), is introduced, which is loosely based on CLUMPFIND, but that also generates a structure tree, or dendrogram, which can be used to interpret hierarchical clump structure in a complex region. Sources were extracted from the data using the hierarchical version of CSAR and plotted on a mass-size diagram. We found a hierarchy of objects with sizes in the range 0.024-2.7 pc. Previous studies showed that gravitationally bound prestellar cores and unbound starless clumps appeared in different places on the mass-size diagram. However, it was unclear whether this was due to a lack of instrumental dynamic range or whether they were actually two distinct populations. The excellent sensitivity of Herschel shows that our sources fill the gap in the mass-size plane between starless and pre-stellar cores, and gives the first clear supporting observational evidence for the theory that unbound clumps and (gravitationally bound) prestellar cores are all part of the same population, and hence presumably part of the same evolutionary sequence (c.f. Simpson et al. 2011).

preprint2013arXiv

Global collapse of molecular clouds as a formation mechanism for the most massive stars

The relative importance of primordial molecular cloud fragmentation versus large-scale accretion still remains to be assessed in the context of massive core/star formation. Studying the kinematics of the dense gas surrounding massive-star progenitors can tell us the extent to which large-scale flow of material impacts the growth in mass of star-forming cores. Here we present a comprehensive dataset of the 5500(+/-800) Msun infrared dark cloud SDC335.579-0.272 (hereafter SDC335) which exhibits a network of cold, dense, parsec-long filaments. Atacama Large Millimeter Array (ALMA) Cycle 0 observations reveal two massive star-forming cores, MM1 and MM2, sitting at the centre of SDC335 where the filaments intersect. With a gas mass of 545(+770,-385) Msun contained within a source diameter of 0.05pc, MM1 is one of the most massive, compact protostellar cores ever observed in the Galaxy. As a whole, SDC335 could potentially form an OB cluster similar to the Trapezium cluster in Orion. ALMA and Mopra single-dish observations of the SDC335 dense gas furthermore reveal that the kinematics of this hub-filament system are consistent with a global collapse of the cloud. These molecular-line data point towards an infall velocity V_{inf} =0.7(+/-0.2) km/s, and a total mass infall rate \dot{M}_{inf} = 2.5(+/-1.0) x 10^{-3} Msun/yr towards the central pc-size region of SDC335. This infall rate brings 750(+/-300) Msun of gas to the centre of the cloud per free-fall time (t_{ff}=3x10^5 yr). This is enough to double the mass already present in the central pc-size region in 3.5(+2.2,-1.0) x t_{ff}. These values suggest that the global collapse of SDC335 over the past million year resulted in the formation of an early O-type star progenitor at the centre of the cloud's gravitational potential well.

preprint2013arXiv

Herschel Observations of the W3 GMC: Clues to the Formation of Clusters of High-Mass Stars

The W3 GMC is a prime target for the study of the early stages of high-mass star formation. We have used Herschel data from the HOBYS key program to produce and analyze column density and temperature maps. Two preliminary catalogs were produced by extracting sources from the column density map and from Herschel maps convolved to the 500 micron resolution. Herschel reveals that among the compact sources (FWHM<0.45 pc), W3 East, W3 West, and W3 (OH) are the most massive and luminous and have the highest column density. Considering the unique properties of W3 East and W3 West, the only clumps with on-going high-mass star formation, we suggest a 'convergent constructive feedback' scenario to account for the formation of a cluster with decreasing age and increasing system/source mass toward the innermost regions. This process, which relies on feedback by high-mass stars to ensure the availability of material during cluster formation, could also lead to the creation of an environment suitable for the formation of Trapezium-like systems. In common with other scenarios proposed in other HOBYS studies, our results indicate that an active/dynamic process aiding in the accumulation, compression, and confinement of material is a critical feature of the high-mass star/cluster formation, distinguishing it from classical low-mass star formation. The environmental conditions and availability of triggers determine the form in which this process occurs, implying that high-mass star/cluster formation could arise from a range of scenarios: from large scale convergence of turbulent flows, to convergent constructive feedback or mergers of filaments.

preprint2013arXiv

Herschel Reveals Massive Cold Clumps in NGC 7538

We present the first overview of the Herschel observations of the nearby high-mass star-forming region NGC 7538, taken as part of the Herschel imaging study of OB Young Stellar objects (HOBYS) Key Programme. These PACS and SPIRE maps cover an approximate area of one square degree at five submillimeter and far-infrared wavebands. We have identified 780 dense sources and classified 224 of those. With the intention of investigating the existence of cold massive starless or class 0-like clumps that would have the potential to form intermediate- to high-mass stars, we further isolate 13 clumps as the most likely candidates for followup studies. These 13 clumps have masses in excess of 40 M_sun and temperatures below 15 K. They range in size from 0.4 pc to 2.5 pc and have densities between 3x10^3 cm^-3 to 4x10^4 cm^-3. Spectral energy distributions are then used to characterize their energetics and evolutionary state through a luminosity-mass diagram. NGC 7538 has a highly filamentary structure, previously unseen in the dust continuum of existing submillimeter surveys. We report the most complete imaging to date of a large, evacuated ring of material in NGC 7538 which is bordered by many cool sources.

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

Low-velocity shocks traced by extended SiO emission along the W43 ridges: witnessing the formation of young massive clusters

The formation of high-mass stars is tightly linked to that of their parental clouds. We here focus on the high-density parts of W43, a molecular cloud undergoing an efficient event of formation. The cloud structure is studied with a column density image derived from Herschel continuum maps obtained at 70, 160, 250, 350, and 500 micron. We identify two high-column density filamentary clouds, quoted as the W43-MM1 and W43-MM2 ridges, which both account for 1.5x10^4 Msun gas mass above 10^23 cm-2 and within areas of 5 and 14pc^2, respectively. We used the N_2H^+ 1--0 line to confirm that the W43-MM1 and W43-MM2 ridges are structures coherent in velocity and gravitationally bound, despite their large velocity dispersion and ~5 kms line widths. The most intriguing result of the W43 large program is the bright wide-spread SiO 2--1 emission: 1--11 K kms$ stretching an area of ~28 pc^2. Concentrated toward the W43-MM1 and W43-MM2 ridges and their immediate surroundings, it leads to a total luminosity of L_SiO 2-1 ~4 10^4 K kms kpc^2pc^2. We measured a steep relation between the luminosity and velocity extent of the SiO~2--1 lines and propose to use it to distinguish the low-velocity shocks observed here from the more classical high-velocity ones associated with outflows of high-mass young stellar objects. We used state-of-the-art shock models to demonstrate that low-velocity (<10 kms^-1) shocks with a small amount (10%) of Si atoms initially in gas phase or in grain mantles can explain the observed SiO column density in W43. The spatial and velocity overlaps between the ridges high-density gas (n_H2>10^4-10^5 cm^-3) and the shocked SiO gas suggests that ridges could be forming via colliding flows driven by gravity and accompanied by low-velocity shocks. This mechanism may be the initial conditions for the formation of young massive clusters in these ridges.

preprint2013arXiv

MALT90: The Millimetre Astronomy Legacy Team 90 GHz Survey

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

preprint2013arXiv

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

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

preprint2013arXiv

Probing the Turbulent Ambipolar Diffusion Scale in Molecular Clouds with Spectroscopy

We estimate the turbulent ambipolar diffusion length scale and magnetic field strength in the massive dense cores CygX-N03 and CygX-N53, located in the Cygnus-X star-forming region. The method we use requires comparing the velocity dispersions in the spectral line profiles of the coexistent ion and neutral pair H13CN and H13CO+ (J=1-0) at different length scales. We fit Kolmogorov-type power laws to the lower envelopes of the velocity dispersion spectra of the two species. This allows to calculate the turbulent ambipolar diffusion scale, which in turn determines the plane-of-the-sky magnetic field strength. We find turbulent ambipolar diffusion length scales of 3.8+-0.1 mpc and 21.2+-0.4 mpc, and magnetic field strengths of 0.33 mG and 0.76 mG for CygX-N03 and CygX-N53, respectively. These magnetic field values have uncertainties of a factor of a few. Despite a lower signal-to-noise ratio of the data in CygX-N53 than in CygX-N03, and the caveat that its stronger field might stem in part from projection effects, the difference in field strengths suggests different fragmentation activities of the two cores. Even though the quality of our data, obtained with the IRAM Plateau de Bure Interferometer (PdBI), is somewhat inferior to previous single-dish data, we demonstrate that this method is suited also for observations at high spatial resolution.

preprint2013arXiv

Two Mass Distributions in the L 1641 Molecular Clouds: The Herschel connection of Dense Cores and Filaments in Orion A

We present the Herschel Gould Belt survey maps of the L1641 molecular clouds in Orion A. We extracted both the filaments and dense cores in the region. We identified which of dense sources are proto- or pre-stellar, and studied their association with the identified filaments. We find that although most (71%) of the pre-stellar sources are located on filaments there is still a significant fraction of sources not associated with such structures. We find that these two populations (on and off the identified filaments) have distinctly different mass distributions. The mass distribution of the sources on the filaments is found to peak at 4 Solar masses and drives the shape of the CMF at higher masses, which we fit with a power law of the form dN/dlogM \propto M^{-1.4+/-0.4}. The mass distribution of the sources off the filaments, on the other hand, peaks at 0.8 Solar masses and leads to a flattening of the CMF at masses lower than ~4 Solar masses. We postulate that this difference between the mass distributions is due to the higher proportion of gas that is available in the filaments, rather than in the diffuse cloud.

preprint2013arXiv

What determines the density structure of molecular clouds ? A case study of Orion B with Herschel

A key parameter to the description of all star formation processes is the density structure of the gas. In this letter, we make use of probability distribution functions (PDFs) of Herschel column density maps of Orion B, Aquila, and Polaris, obtained with the Herschel Gould Belt survey (HGBS). We aim to understand which physical processes influence the PDF shape, and with which signatures. The PDFs of Orion B (Aquila) show a lognormal distribution for low column densities until Av 3 (6), and a power-law tail for high column densities, consistent with a rho r^-2 profile for the equivalent spherical density distribution. The PDF of Orion B is broadened by external compression due to the nearby OB stellar aggregates. The PDF of a quiescent subregion of the non-star-forming Polaris cloud is nearly lognormal, indicating that supersonic turbulence governs the density distribution. But we also observe a deviation from the lognormal shape at Av>1 for a subregion in Polaris that includes a prominent filament. We conclude that (i) the point where the PDF deviates from the lognormal form does not trace a universal Av-threshold for star formation, (ii) statistical density fluctuations, intermittency and magnetic fields can cause excess from the lognormal PDF at an early cloud formation stage, (iii) core formation and/or global collapse of filaments and a non-isothermal gas distribution lead to a power-law tail, and (iv) external compression broadens the column density PDF, consistent with numerical simulations.

preprint2012arXiv

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

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

preprint2012arXiv

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

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

preprint2012arXiv

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

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

preprint2012arXiv

Herschel Observations of a Potential Core Forming Clump: Perseus B1-E

We present continuum observations of the Perseus B1-E region from the Herschel Gould Belt Survey. These Herschel data reveal a loose grouping of substructures at 160 - 500 micron not seen in previous submillimetre observations. We measure temperature and column density from these data and select the nine densest and coolest substructures for follow-up spectral line observations with the Green Bank Telescope. We find that the B1-E clump has a mass of ~ 100 solar masses and appears to be gravitationally bound. Furthermore, of the nine substructures examined here, one substructure (B1-E2) appears to be itself bound. The substructures are typically less than a Jeans length from their nearest neighbour and thus, may interact on a timescale of ~ 1 Myr. We propose that B1-E may be forming a first generation of dense cores, which could provide important constraints on the initial conditions of prestellar core formation. Our results suggest that B1-E may be influenced by a strong, localized magnetic field, but further observations are still required.

preprint2012arXiv

Herschel view of the Taurus B211/3 filament and striations: Evidence of filamentary growth?

We present first results from the Herschel Gould Belt survey for the B211/L1495 region in the Taurus molecular cloud. Thanks to their high sensitivity and dynamic range, the Herschel images reveal the structure of the dense, star-forming filament B211 with unprecedented detail, along with the presence of striations perpendicular to the filament and generally oriented along the magnetic field direction as traced by optical polarization vectors. Based on the column density and dust temperature maps derived from the Herschel data, we find that the radial density profile of the B211 filament approaches a power-law behavior ρ {\propto} r^(-2.0{\pm}0.4) at large radii and that the temperature profile exhibits a marked drop at small radii. The observed density and temperature profiles of the B211 filament are in good agreement with a theoretical model of a cylindrical filament undergoing gravitational contraction with a polytropic equation of state: P {\propto} ρ^γ and T {\propto} ρ^(γ-1), with γ=0.97{\pm}0.01<1 (i.e. not strictly isothermal). The morphology of the column density map, where some of the perpendicular striations are apparently connected to the B211 filament, further suggests that the material may be accreting along the striations onto the main filament. The typical velocities expected for the infalling material in this picture are ~0.5-1 km/s, which are consistent with the existing kinematical constraints from previous CO observations.

preprint2012arXiv

Interstellar matter and star formation in W5-E - A Herschel view

W5-E has been observed with the Herschel-PACS and -SPIRE photometers, at 100, 160, 250, 350, and 500 microns. The dust temperature map shows a rather uniform temperature, in the range 17.5-20 K in the dense condensations or filaments, 21-22 K in the photodissociation regions, and 24-31 K in the direction of the ionized regions. The column densities are rather low, everywhere lower than 10^23 cm-2, and of the order of a few 10^21 cm-2 in the PDRs. About 8000 solar masses of neutral material surrounds the ionized region, which is low with respect to the volume of this HII region; we suggest that the exciting stars of the W5-E, W5-W, Sh~201, A and B HII regions formed along a dense filament or sheet rather than inside a more spherical cloud. Fifty point sources have been detected at 100 microns. Most of them are Class 0/I YSOs. The SEDs of their envelopes have been fitted using a modified blackbody model. These envelopes are cold, with a mean temperature of 15.7+-1.8K. Their masses are in the range 1.3-47 solar masses. Eleven of these point sources are candidate Class 0 YSOs. Twelve of these point sources are possibly at the origin of bipolar outflows detected in this region. None of the YSOs contain a massive central object, but a few may form a massive star as they have both a massive envelope and also a high envelope accretion rate. Most of the Class 0/I YSOs are observed in the direction of high column density material, for example in the direction of the massive condensations present at the waist of the bipolar Sh 201 HII region or enclosed by the bright-rimmed cloud BRC14. The overdensity of Class 0/I YSOs on the borders of the HII regions strongly suggests that triggered star formation is at work in this region but, due to insufficient resolution, the exact processes at the origin of the triggering are difficult to determine.

preprint2012arXiv

Recent star formation in the Lupus clouds as seen by Herschel

We present a study of the star formation histories of the Lupus I, III, and IV clouds using the Herschel 70-500 micron maps obtained by the Herschel Gould Belt Survey Key Project. By combining the new Herschel data with the existing Spitzer catalog we obtained an unprecedented census of prestellar sources and young stellar objects in the Lupus clouds, which allowed us to study the overall star formation rate (SFR) and efficiency (SFE). The high SFE of Lupus III, its decreasing SFR, and its large number of pre-main sequence stars with respect to proto- and prestellar sources, suggest that Lupus III is the most evolved cloud, and after having experienced a major star formation event in the past, is now approaching the end of its current star-forming cycle. Lupus I is currently undergoing a large star formation event, apparent by the increasing SFR, the large number of prestellar objects with respect to more evolved objects, and the high percentage of material at high extinction (e.g., above A_V=8 mag). Also Lupus IV has an increasing SFR; however, the relative number of prestellar sources is much lower, suggesting that its star formation has not yet reached its peak.

preprint2012arXiv

Resolving the Vela C ridge with P-ArTeMiS and Herschel

We present APEX/P-ArTéMiS 450μm continuum observations of RCW 36 and the adjacent ridge, a high-mass high-column density filamentary structure at the centre of the Vela C molecular cloud. These observations, at higher resolution than Herschel's SPIRE camera, reveal clear fragmentation of the central star-forming ridge. Combined with PACS far-infrared and SPIRE sub-millimetre observations from the Herschel HOBYS project we build a high resolution column density map of the region mapped with P-ArTéMiS. We extract the radial density profile of the Vela C ridge which with a ~ 0.1pc central width is consistent with that measured for low-mass star-forming filaments in the Herschel Gould Belt survey. Direct comparison with Serpens South, of the Gould Belt Aquila complex, reveals many similarities between the two regions. Despite likely different formation mechanisms and histories, the Vela C ridge and Serpens South filament share common characteristics, including their filament central widths.

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 Dust Properties of Bubble HII Regions as seen by Herschel

Because of their relatively simple morphology, "bubble" HII regions have been instrumental to our understanding of star formation triggered by HII regions. With the far-infrared (FIR) spectral coverage of the Herschel satellite, we can access the wavelengths where these regions emit the majority of their energy through their dust emission. At Herschel wavelengths 70 micron to 500 micron, the emission associated with HII regions is dominated by the cool dust in their photodissociation regions (PDRs). We find average dust temperatures of 26K along the PDRs, with little variation between the HII regions in the sample, while local filaments and infrared dark clouds average 19K and 15K respectively. Higher temperatures lead to higher values of the Jeans mass, which may affect future star formation. The mass of the material in the PDR, collected through the expansion of the HII region, is between ~300 and ~10,000 Solar masses for the HII regions studied here. These masses are in rough agreement with the expected masses swept up during the expansion of the \hii\ regions. Approximately 20% of the total FIR emission is from the direction of the bubble central regions. This suggests that we are detecting emission from the "near-side" and "far-side" PDRs along the line of sight and that bubbles are three-dimensional structures. We find only weak support for a relationship between dust temperature and beta, of a form similar to that caused by noise and calibration uncertainties alone.

preprint2012arXiv

The Herschel view of the on-going star formation in the Vela-C molecular cloud

As part of the Herschel guaranteed time key program 'HOBYS', we present the photometric survey of the star forming region Vela-C, one of the nearest sites of low-to-high-mass star formation in the Galactic plane. Vela-C has been observed with PACS and SPIRE in parallel mode between 70 um and 500 um over an area of about 3 square degrees. A photometric catalogue has been extracted from the detections in each band, using a threshold of 5 sigma over the local background. Out of this catalogue we have selected a robust sub-sample of 268 sources, of which 75% are cloud clumps and 25% are cores. Their Spectral Energy Distributions (SEDs) have been fitted with a modified black body function. We classify 48 sources as protostellar and 218 as starless. For two further sources, we do not provide a secure classification, but suggest they are Class 0 protostars. From SED fitting we have derived key physical parameters. Protostellar sources are in general warmer and more compact than starless sources. Both these evidences can be ascribed to the presence of an internal source(s) of moderate heating, which also causes a temperature gradient and hence a more peaked intensity distribution. Moreover, the reduced dimensions of protostellar sources may indicate that they will not fragment further. A virial analysis of the starless sources gives an upper limit of 90% for the sources gravitationally bound and therefore prestellar. We fit a power law N(logM) prop M^-1.1 to the linear portion of the mass distribution of prestellar sources. This is in between that typical of CO clumps and those of cores in nearby star-forming regions. We interpret this as a result of the inhomogeneity of our sample, which is composed of comparable fractions of clumps and cores.

preprint2012arXiv

The M16 molecular complex under the influence of NGC6611. Herschel's perspective of the heating effect on the Eagle Nebula

We present Herschel images from the HOBYS key program of the Eagle Nebula (M16) in the far-infrared and sub-millimetre, using the PACS and SPIRE cameras at 70μm, 160μm, 250μm, 350μm, 500μm. M16, home to the Pillars of Creation, is largely under the influence of the nearby NGC6611 high-mass star cluster. The Herschel images reveal a clear dust temperature gradient running away from the centre of the cavity carved by the OB cluster. We investigate the heating effect of NGC6611 on the entire M16 star-forming complex seen by Herschel including the diffuse cloud environment and the dense filamentary structures identified in this region. In addition, we interpret the three-dimensional geometry of M16 with respect to the nebula, its surrounding environment, and the NGC6611 cavity. The dust temperature and column density maps reveal a prominent eastern filament running north-south and away from the high-mass star-forming central region and the NGC6611 cluster, as well as a northern filament which extends around and away from the cluster. The dust temperature in each of these filaments decreases with increasing distance from the NGC6611 cluster, indicating a heating penetration depth of \sim 10 pc in each direction in 3 - 6 \times 10^{22} cm-2 column density filaments. We show that in high-mass star-forming regions OB clusters impact the temperature of future star-forming sites, modifying the initial conditions for collapse and effecting the evolutionary criteria of protostars developed from spectral energy distributions. Possible scenarios for the origin of the morphology seen in this region are discussed, including a western equivalent to the eastern filament, which was destroyed by the creation of the OB cluster and its subsequent winds and radiation.

preprint2012arXiv

The massive protostar W43-MM1 as seen by Herschel-HIFI water spectra: high turbulence and accretion luminosity

We present Herschel/HIFI observations of fourteen water lines in W43-MM1, a massive protostellar object in the luminous star cluster-forming region W43. We analyze the gas dynamics from the line profiles using Herschel-HIFI observations (WISH-KP) of fourteen far-IR water lines (H2O, H217O, H218O), CS(11-10), and C18O(9-8) lines, and using our modeling of the continuum spectral energy distribution. As for lower mass protostellar objects, the molecular line profiles are a mix of emission and absorption, and can be decomposed into 'medium', and 'broad' velocity components. The broad component is the outflow associated with protostars of all masses. Our modeling shows that the remainder of the water profiles can be well fitted by an infalling and passively heated envelope, with highly supersonic turbulence varying from 2.2 km/s in the inner region to 3.5 km/s in the outer envelope. Also, W43-MM1 has a high accretion rate, between 4.0 x 10^{-4} and 4.0 x 10^{-2} \msun /yr, derived from the fast (0.4-2.9 km/s) infall observed. We estimate a lower mass limit of gaseous water of 0.11 \msun and total water luminosity of 1.5 \lsun (in the 14 lines presented here). The central hot core is detected with a water abundance of 1.4 x 10^{-4} while the water abundance for the outer envelope is 8 x10^{-8}. The latter value is higher than in other sources, most likely related to the high turbulence and the micro-shocks created by its dissipation. Examining water lines of various energies, we find that the turbulent velocity increases with the distance to the center. While not in clear disagreement with the competitive accretion scenario, this behavior is predicted by the turbulent core model. Moreover, the estimated accretion rate is high enough to overcome the expected radiation pressure.

preprint2012arXiv

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

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

preprint2011arXiv

Characterizing interstellar filaments with Herschel in IC5146

We provide a first look at the results of the Herschel Gould Belt survey toward the IC5146 molecular cloud and present a preliminary analysis of the filamentary structure in this region. The column density map, derived from our 70-500 micron Herschel data, reveals a complex network of filaments, and confirms that these filaments are the main birth sites of prestellar cores. We analyze the column density profiles of 27 filaments and show that the underlying radial density profiles fall off as r^{-1.5} to r^{-2.5} at large radii. Our main result is that the filaments seem to be characterized by a narrow distribution of widths having a median value of 0.10 +- 0.03 pc, which is in stark contrast to a much broader distribution of central Jeans lengths. This characteristic width of ~0.1 pc corresponds to within a factor of ~2 to the sonic scale below which interstellar turbulence becomes subsonic in diffuse gas, supporting the argument that the filaments may form as a result of the dissipation of large-scale turbulence.

preprint2011arXiv

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

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

preprint2011arXiv

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

The end of star formation in Chamaeleon I ? A LABOCA census of starless and protostellar cores

Chamaeleon I is the most active region in terms of star formation in the Chamaeleon molecular cloud complex. Its population of prestellar and protostellar cores is not known and a controversy exists concerning its history of star formation. Our goal is to characterize the earliest stages of star formation in this cloud. We used the bolometer array LABOCA at APEX to map the cloud in dust continuum emission at 870 micron. The detected sources are analysed by carefully taking into account the spatial filtering inherent in the data reduction process. A search for associations with YSOs is performed using Spitzer data and the SIMBAD database. Most of the detected 870 micron emission is distributed in 5 filaments. We identify 59 starless cores, one candidate first hydrostatic core, and 21 sources associated with more evolved YSOs. The starless cores are only found above a visual extinction threshold of 5 mag. They are less dense than those detected in other nearby molecular clouds by a factor of a few on average. The core mass distribution is consistent with the IMF at the high-mass end but is overpopulated at the low-mass end. In addition, at most 17% of the cores have a mass larger than the critical Bonnor-Ebert mass. Both results suggest that a large fraction of the starless cores may not be prestellar. Based on the census of prestellar cores, Class 0 protostars, and more evolved YSOs, we conclude that the star formation rate has decreased with time in this cloud. The low fraction of candidate prestellar cores among the population of starless cores, the small number of Class 0 protostars, the high global star formation efficiency, the decrease of the star formation rate with time, and the low mass per unit length of the detected filaments all suggest that we may be witnessing the end of the star formation process in Cha I {abridged}.

preprint2011arXiv

The Herschel view of massive star formation in G035.39--00.33: Dense and cold filament of W48 undergoing a mini-starburst

The filament IRDC G035.39--00.33 in the W48 molecular complex is one of the darkest infrared clouds observed by \textit{Spitzer}. It has been observed by the PACS (70 and 160\,$\micron$) and SPIRE (250, 350, and 500\,$\micron$) cameras of the \textit{Herschel} Space Observatory as part of the W48 molecular cloud complex in the framework of the HOBYS key programme. The observations reveal a sample of 28 compact sources (deconvolved FWHM sizes $<$0.3 pc) complete down to $\sim$$5 \msun$ in G035.39--00.33 and its surroundings. Among them, 13 compact sources are massive dense cores with masses $>$$20 \msun$. The cloud characteristics we derive from the analysis of their spectral energy distributions are masses of $20-50 \msun$, sizes of 0.1--0.2 pc, and average densities of $2-20 \times 10^{5} \cmc$, which make these massive dense cores excellent candidates to form intermediate- to high-mass stars. Most of the massive dense cores are located inside the G035.39--00.33 ridge and host IR-quiet high-mass protostars. The large number of protostars found in this filament suggests that we are witnessing a mini-burst of star formation with an efficiency of $\sim$15% and a rate density of $\sim$$40 \msun\,$yr$^{-1}\,$kpc$^{-2}$ within $\sim$8 pc$^2$, a large area covering the full ridge. Part of the extended SiO emission observed towards G035.39--00.33 is not associated with obvious protostars and may originate from low-velocity shocks within converging flows, as advocated by previous studies.

preprint2011arXiv

W43: the closest molecular complex of the Galactic Bar?

We used a large database extracted from Galaxy-wide surveys of H {\scriptsize I}, $^{13}$CO 1-0, $8 \micron$ and $870 \micron$ continuum to trace diffuse atomic gas, low- to medium-density molecular gas, high-density molecular gas, and star formation activity which we complemented by dedicated $^{12}$CO 2--1, 3--2 observations of the W43 region. From the detailed 3D (space-space-velocity) analysis of the molecular and atomic cloud tracers through the region and despite its wide velocity range (\emph{FWHM}$\sim22.3 \kms$ around $\vlsr\sim 95.9 \kms$), we identified W43 as a large (equivalent diameter $\sim 140$ pc) and coherent complex of molecular clouds which is surrounded by an atomic gas envelope (equivalent diameter $\sim 290$ pc). We measured the total mass of this newly-identified molecular complex ($M_{\tiny total}\sim 7.1 \times 10^6 \msun$), the mass contained in dense $870 \micron$ clumps ($<5$ pc dense cloud structures, $M_{\tiny clumps}\sim 8.4 \times 10^5 \msun$) and conclude that W43 is particularly massive and concentrated. The distance we assume for the W43 complex is 6 kpc from the Sun, which may place it at the meeting point of the Scutum-Centaurus (or Scutum-Crux) Galactic arm and the Bar, a dynamically complex region where high-velocity streams could easily collide. The star formation rate of W43 is suggested not to be steady but it is increasing from $\, \sim 0.01 \msun\, {yr}^{-1}$ (measured from its $8\,μ$m luminosity) to $\sim 0.1 \msun\, {yr}^{-1}$ (measured from its molecular content). From the global properties of W43, we claim it is an extreme molecular complex in the Milky Way and it could even be forming starburst clusters in the near future.

preprint2011arXiv

Will the starless cores in Chamaeleon I and III turn prestellar?

The nearby Chamaeleon molecular cloud complex is a good laboratory to study the process of low-mass star formation since it consists of three clouds with very different properties. Cha III does not show any sign of star formation, while star formation has been very active in Cha I and may already be finishing. Our goal is to determine whether star formation can proceed in Cha III, and to compare the results to our recent survey of Cha I. We used the Large APEX Bolometer Array (LABOCA) to map Cha III in dust continuum emission at 870 micron. 29 sources are extracted from the map, all of them being starless. The starless cores are found down to a visual extinction of 1.9 mag, in marked contrast with other molecular clouds, including Cha I. Apart from this difference, the Cha III starless cores share very similar properties with those found in Cha I. At most two sources have a mass larger than the critical Bonnor-Ebert mass, which suggests that the fraction of prestellar cores is very low, even lower than in Cha I. Only the most massive sources are candidate prestellar cores, in agreement with the correlation found earlier in the Pipe nebula. The mass distribution of the 85 starless cores of Cha I and III that are not candidate prestellar cores is consistent with a single power law down to the 90% completeness limit, with an exponent close to the Salpeter value. A fraction of the starless cores in Cha I and III may still grow in mass and become gravitationally unstable. Based on predictions of numerical simulations of turbulent molecular clouds, we estimate that at most 50% and 20% of the starless cores of Cha I and III, respectively, may form stars. The LABOCA survey reveals that Cha III, and Cha I to some extent too, is a prime target to study the formation of prestellar cores, and thus the onset of star formation. (abridged).

preprint2010arXiv

A Herschel study of the properties of starless cores in the Polaris Flare dark cloud region using PACS and SPIRE

The Polaris Flare cloud region contains a great deal of extended emission. It is at high declination and high Galactic latitude. It was previously seen strongly in IRAS Cirrus emission at 100 microns. We have detected it with both PACS and SPIRE on Herschel. We see filamentary and low-level structure. We identify the five densest cores within this structure. We present the results of a temperature, mass and density analysis of these cores. We compare their observed masses to their virial masses, and see that in all cases the observed masses lie close to the lower end of the range of estimated virial masses. Therefore, we cannot say whether they are gravitationally bound prestellar cores. Nevertheless, these are the best candidates to be potentialprestellar cores in the Polaris cloud region.

preprint2010arXiv

Clouds, filaments and protostars: the Herschel Hi-GAL Milky Way

We present the first results from the science demonstration phase for the Hi-GAL survey, the Herschel key-project that will map the inner Galactic Plane of the Milky Way in 5 bands. We outline our data reduction strategy and present some science highlights on the two observed 2° x 2° tiles approximately centered at l=30° and l=59°. The two regions are extremely rich in intense and highly structured extended emission which shows a widespread organization in filaments. Source SEDs can be built for hundreds of objects in the two fields, and physical parameters can be extracted, for a good fraction of them where the distance could be estimated. The compact sources (which we will call 'cores' in the following) are found for the most part to be associated with the filaments, and the relationship to the local beam-averaged column density of the filament itself shows that a core seems to appear when a threshold around A_V of about 1 is exceeded for the regions in the l=59° field; a A_V value between 5 and 10 is found for the l=30° field, likely due to the relatively larger distances of the sources. This outlines an exciting scenario where diffuse clouds first collapse into filaments, which later fragment to cores where the column density has reached a critical level. In spite of core L/M ratios being well in excess of a few for many sources, we find core surface densities between 0.03 and 0.5 g cm-2. Our results are in good agreement with recent MHD numerical simulations of filaments forming from large-scale converging flows.

preprint2010arXiv

Dynamic star formation in the massive DR21 filament

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

preprint2010arXiv

Filamentary structures and compact objects in the Aquila and Polaris clouds observed by Herschel

Our PACS and SPIRE images of the Aquila Rift and part of the Polaris Flare regions, taken during the science demonstration phase of Herschel discovered fascinating, omnipresent filamentary structures that appear to be physically related to compact cores. We briefly describe a new multi-scale, multi-wavelength source extraction method used to detect objects and measure their parameters in our Herschel images. All of the extracted starless cores (541 in Aquila and 302 in Polaris) appear to form in the long and very narrow filaments. With its combination of the far-IR resolution and sensitivity, Herschel directly reveals the filaments in which the dense cores are embedded; the filaments are resolved and have deconvolved widths of 35 arcsec in Aquila and 59 arcsec in Polaris (9000 AU in both regions). Our first results of observations with Herschel enable us to suggest that in general dense cores may originate in a process of fragmentation of complex networks of long, thin filaments, likely formed as a result of an interplay between gravity, interstellar turbulence, and magnetic fields. To unravel the roles of the processes, one has to obtain additional kinematic and polarization information; these follow-up observations are planned.

preprint2010arXiv

From filamentary clouds to prestellar cores to the stellar IMF: Initial highlights from the Herschel Gould Belt survey

We summarize the first results from the Gould Belt survey, obtained toward the Aquila Rift and Polaris Flare regions during the 'science demonstration phase' of Herschel. Our 70-500 micron images taken in parallel mode with the SPIRE and PACS cameras reveal a wealth of filamentary structure, as well as numerous dense cores embedded in the filaments. Between ~ 350 and 500 prestellar cores and ~ 45-60 Class 0 protostars can be identified in the Aquila field, while ~ unbound starless cores and no protostars are observed in the Polaris field. The prestellar core mass function (CMF) derived for the Aquila region bears a strong resemblance to the stellar initial mass function (IMF), already confirming the close connection between the CMF and the IMF with much better statistics than earlier studies. Comparing and contrasting our Herschel results in Aquila and Polaris, we propose an observationally-driven scenario for core formation according to which complex networks of long, thin filaments form first within molecular clouds, and then the densest filaments fragment into a number of prestellar cores via gravitational instability.

preprint2010arXiv

Gas dynamics in Massive Dense Cores in Cygnus-X

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

preprint2010arXiv

Herschel observations of embedded protostellar clusters in the Rosette Molecular Cloud

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

preprint2010arXiv

Herschel-HIFI observations of high-J CO lines in the NGC 1333 low-mass star-forming region

Herschel-HIFI observations of high-J lines (up to J_u=10) of 12CO, 13CO and C18O are presented toward three deeply embedded low-mass protostars, NGC 1333 IRAS 2A, IRAS 4A, and IRAS 4B, obtained as part of the Water In Star-forming regions with Herschel (WISH) key program. The spectrally-resolved HIFI data are complemented by ground-based observations of lower-J CO and isotopologue lines. The 12CO 10-9 profiles are dominated by broad (FWHM 25-30 km s^-1) emission. Radiative transfer models are used to constrain the temperature of this shocked gas to 100-200 K. Several CO and 13CO line profiles also reveal a medium-broad component (FWHM 5-10 km s^-1), seen prominently in H2O lines. Column densities for both components are presented, providing a reference for determining abundances of other molecules in the same gas. The narrow C18O 9-8 lines probe the warmer part of the quiescent envelope. Their intensities require a jump in the CO abundance at an evaporation temperature around 25 K, thus providing new direct evidence for a CO ice evaporation zone around low-mass protostars.

preprint2010arXiv

HERSCHEL-HIFI spectroscopy of the intermediate mass protostar NGC7129 FIRS 2

HERSCHEL-HIFI observations of water from the intermediate mass protostar NGC7129 FIRS 2 provide a powerful diagnostic of the physical conditions in this star formation environment. Six spectral settings, covering four H216O and two H218O lines, were observed and all but one H218O line were detected. The four H2 16 O lines discussed here share a similar morphology: a narrower, \approx 6 km/s, component centered slightly redward of the systemic velocity of NGC7129 FIRS 2 and a much broader, \approx 25 km/s component centered blueward and likely associated with powerful outflows. The narrower components are consistent with emission from water arising in the envelope around the intermediate mass protostar, and the abundance of H2O is constrained to \approx 10-7 for the outer envelope. Additionally, the presence of a narrow self-absorption component for the lowest energy lines is likely due to self-absorption from colder water in the outer envelope. The broader component, where the H2O/CO relative abundance is found to be \approx 0.2, appears to be tracing the same energetic region that produces strong CO emission at high J.

preprint2010arXiv

Herschel-PACS spectroscopy of the intermediate mass protostar NGC7129 FIRS 2

Aims: We present preliminary results of the first Herschel spectroscopic observations of NGC7129 FIRS2, an intermediate mass star-forming region. We attempt to interpret the observations in the framework of an in-falling spherical envelope. Methods: The PACS instrument was used in line spectroscopy mode (R=1000-5000) with 15 spectral bands between 63 and 185 microns. This provided good detections of 26 spectral lines seen in emission, including lines of H2O, CO, OH, O I, and C II. Results: Most of the detected lines, particularly those of H2O and CO, are substantially stronger than predicted by the spherical envelope models, typically by several orders of magnitude. In this paper we focus on what can be learned from the detected CO emission lines. Conclusions: It is unlikely that the much stronger than expected line emission arises in the (spherical) envelope of the YSO. The region hot enough to produce such high excitation lines within such an envelope is too small to produce the amount of emission observed. Virtually all of this high excitation emission must arise in structures such as as along the walls of the outflow cavity with the emission produced by a combination of UV photon heating and/or non-dissociative shocks.

preprint2010arXiv

Herschel-SPIRE spectroscopy of G29.96-0.02: fitting the full SED

We use the SPIRE Fourier-Transform Spectrometer (FTS) on-board the ESA Herschel Space Telescope to analyse the submillimetre spectrum of the Ultra-compact HII region G29.96-0.02. Spectral lines from species including 13CO, CO, [CI], and [NII] are detected. A sparse map of the [NII] emission shows at least one other HII region neighbouring the clump containing the UCHII. The FTS spectra are combined with ISO SWS and LWS spectra and fluxes from the literature to present a detailed spectrum of the source spanning three orders of magnitude in wavelength. The quality of the spectrum longwards of 100 μm allows us to fit a single temperature greybody with temperature 80.3\pm0.6K and dust emissivity index 1.73\pm0.02, an accuracy rarely obtained with previous instruments. We estimate a mass of 1500 Msol for the clump containing the HII region. The clump's bolometeric luminosity of 4 x 10^6 Lsol is comparable to, or slightly greater than, the known O-star powering the UCHII region.

preprint2010arXiv

Herschel-SPIRE spectroscopy of the DR21 molecular cloud core

We present far-infrared spectra and maps of the DR21 molecular cloud core between 196 and 671 microns, using the Herschel-SPIRE spectrometer. Nineteen molecular lines originating from CO, 13CO, HCO+ and H2O, plus lines of [N II] and [CI] were recorded, including several transitions not previously detected. The CO lines are excited in warm gas with Tkin ~ 125 K and nH2 ~ 7 x 10^4 cm-3, CO column density N(CO) ~ 3.5 x 10^18 cm^-2 and a filling factor of ~ 12%, and appear to trace gas associated with an outflow. The rotational temperature analysis incorporating observations from ground-based telescopes reveals an additional lower excitation CO compoment which has a temperature ~ 78 K and N(CO) ~ 4.5 x 10^21 cm^-2. Astronomy & Astrophysics HERSCHEL special Issue, in press.

preprint2010arXiv

Molecular and Atomic Gas in the Local Group Galaxy M33

We present high resolution large scale observations of the molecular and atomic gas in the Local Group Galaxy M33. The observations were carried out using the HERA at the 30m IRAM telescope in the CO(2-1) line achieving a resolution of 12"x2.6 km/s, enabling individual GMCs to be resolved. The observed region mainly along the major axis out to a radius of 8.5 kpc, and covers the strip observed with HIFI/PACS Spectrometers as part of the HERM33ES Herschel key program. The achieved sensitivity in main beam temperature is 20-50 mK at 2.6 km/s velocity resolution. The CO(2-1) luminosity of the observed region is 1.7\pm0.1x10^7 Kkm/s pc^2, corresponding to H2 masses of 1.9x10^8 Msun (including He), calculated with a NH2/ICO twice the Galactic value due to the half-solar metallicity of M33. HI 21 cm VLA archive observations were reduced and the mosaic was imaged and cleaned using the multi-scale task in CASA, yielding a series of datacubes with resolutions ranging from 5" to 25". The HI mass within a radius of 8.5 kpc is estimated to be 1.4x10^9 Msun. The azimuthally averaged CO surface brightness decreases exponentially with a scale length of 1.9\pm0.1 kpc whereas the atomic gas surface density is constant at Sigma_HI=6\pm2 Msun/pc^2 deprojected to face-on. The central kiloparsec H_2 surface density is Sigma_H2=8.5\pm0.2 Msun/pc^2. The star formation rate per unit molecular gas (SF Efficiency, the rate of transformation of molecular gas into stars), as traced by the ratio of CO to Halpha and FIR brightness, is constant with radius. The SFE appears 2-4 times greater than of large spiral galaxies. A morphological comparison of molecular and atomic gas with tracers of star formation shows good agreement between these maps both in terms of peaks and holes. A few exceptions are noted. Several spectra, including those of a molecular cloud situated more than 8 kpc from the galaxy center, are presented.

preprint2010arXiv

Origin of the hot gas in low-mass protostars: Herschel-PACS spectroscopy of HH 46

'Water in Star-forming regions with Herschel' (WISH) is a Herschel Key Programme aimed at understanding the physical and chemical structure of young stellar objects (YSOs) with a focus on water and related species. The low-mass protostar HH 46 was observed with the Photodetector Array Camera and Spectrometer (PACS) on the Herschel Space Observatory to measure emission in H2O, CO, OH, [OI], and [CII] lines located between 63 and 186 um. The excitation and spatial distribution of emission can disentangle the different heating mechanisms of YSOs, with better spatial resolution and sensitivity than previously possible. Far-IR line emission is detected at the position of the protostar and along the outflow axis. The OH emission is concentrated at the central position, CO emission is bright at the central position and along the outflow, and H2O emission is concentrated in the outflow. In addition, [OI] emission is seen in low-velocity gas, assumed to be related to the envelope, and is also seen shifted up to 170 km/s in both the red- and blue-shifted jets. Envelope models are constructed based on previous observational constraints. They indicate that passive heating of a spherical envelope by the protostellar luminosity cannot explain the high-excitation molecular gas detected with PACS, including CO lines with upper levels at >2500 K above the ground state. Instead, warm CO and H2O emission is probably produced in the walls of an outflow-carved cavity in the envelope, which are heated by UV photons and non-dissociative C-type shocks. The bright OH and [OI] emission is attributed to J-type shocks in dense gas close to the protostar. In the scenario described here, the combined cooling by far-IR lines within the central spatial pixel is estimated to be 2 \times 10-2 L_sun, with 60-80% attributed to J- and C-type shocks produced by interactions between the jet and the envelope.

preprint2010arXiv

Physical properties of the Sh2-104 HII region as seen by Herschel

Context: Sh2-104 is a Galactic H ii region with a bubble morphology, detected at optical and radio wavelengths. It is considered the first observational confirmation of the collect-and-collapse model of triggered star-formation. Aims: We aim to analyze the dust and gas properties of the Sh2-104 region to better constrain its effect on local future generations of stars. In addition, we investigate the relationship between the dust emissivity index β and the dust temperature, T_dust. Methods: Using Herschel PACS and SPIRE images at 100, 160, 250, 350 and 500 μm we determine T_dust and β throughout Sh2-104, fitting the spectral energy distributions (SEDs) obtained from aperture photometry. With the SPIRE Fourier transform spectrometer (FTS) we obtained spectra at different positions in the Sh2-104 region. We detect J-ladders of CO and 13CO, with which we derive the gas temperature and column density. We also detect proxies of ionizing flux as the [NII] 3P1-3P0 and [CI] 3P2-3P1 transitions. Results: We find an average value of β ~ 1.5 throughout Sh2-104, as well as a T dust difference between the photodissociation region (PDR, ~ 25 K) and the interior (~ 40 K) of the bubble. We recover the anti-correlation between β and dust temperature reported numerous times in the literature. The relative isotopologue abundances of CO appear to be enhanced above the standard ISM values, but the obtained value is very preliminary and is still affected by large uncertainties.

preprint2010arXiv

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

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

preprint2010arXiv

Star formation triggered by the Galactic HII region RCW 120: First results from the Herschel Space Observatory

By means of different physical mechanisms, the expansion of HII regions can promote the formation of new stars of all masses. RCW 120 is a nearby Galactic HII region where triggered star formation occurs. This region is well-studied - there being a wealth of existing data - and is nearby. However, it is surrounded by dense regions for which far infrared data is essential to obtain an unbiased view of the star formation process and in particular to establish whether very young protostars are present. We attempt to identify all Young Stellar Objects (YSOs), especially those previously undetected at shorter wavelengths, to derive their physical properties and obtain insight into the star formation history in this region. We use Herschel-PACS and -SPIRE images to determine the distribution of YSOs observed in the field. We use a spectral energy distribution fitting tool to derive the YSOs physical properties. Herschel-PACS and -SPIRE images confirm the existence of a young source and allow us to determine its nature as a high-mass (8-10 MSun) Class 0 object (whose emission is dominated by a massive envelope) towards the massive condensation 1 observed at (sub)-millimeter wavelengths. This source was not detected at 24 micron and only barely seen in the MISPGAL 70 micron data. Several other red sources are detected at Herschel wavelengths and coincide with the peaks of the millimeter condensations. SED fitting results for the brightest Herschel sources indicate that, apart from the massive Class 0 that forms in condensation 1, young low mass stars are forming around RCW 120. The YSOs observed on the borders of RCW 120 are younger than its ionizing star, which has an age of about 2.5 Myr.

preprint2010arXiv

The Aquila prestellar core population revealed by Herschel

The origin and possible universality of the stellar initial mass function (IMF) is a major issue in astrophysics. One of the main objectives of the Herschel Gould Belt Survey is to clarify the link between the prestellar core mass function (CMF) and the IMF. We present and discuss the core mass function derived from Herschel data for the large population of prestellar cores discovered with SPIRE and PACS in the Aquila Rift cloud complex at d ~ 260 pc. We detect a total of 541 starless cores in the entire ~11 deg^2 area of the field imaged at 70-500 micron with SPIRE/PACS. Most of these cores appear to be gravitationally bound, and thus prestellar in nature. Our Herschel results confirm that the shape of the prestellar CMF resembles the stellar IMF, with much higher quality statistics than earlier submillimeter continuum ground-based surveys.

preprint2010arXiv

The Herschel first look at protostars in the Aquila Rift

As part of the science demonstration phase of the Herschel mission of the Gould Belt Key Program, the Aquila Rift molecular complex has been observed. The complete ~ 3.3deg x 3.3deg imaging with SPIRE 250/350/500 micron and PACS 70/160 micron allows a deep investigation of embedded protostellar phases, probing of the dust emission from warm inner regions at 70 and 160 micron to the bulk of the cold envelopes between 250 and 500 micron. We used a systematic detection technique operating simultaneously on all Herschel bands to build a sample of protostars. Spectral energy distributions are derived to measure luminosities and envelope masses, and to place the protostars in an M_env - L_bol evolutionary diagram. The spatial distribution of protostars indicates three star-forming sites in Aquila, with W40/Sh2-64 HII region by far the richest. Most of the detected protostars are newly discovered. For a reduced area around the Serpens South cluster, we could compare the Herschel census of protostars with Spitzer results. The Herschel protostars are younger than in Spitzer with 7 Class 0 YSOs newly revealed by Herschel. For the entire Aquila field, we find a total of ~ 45-60 Class 0 YSOs discovered by Herschel. This confirms the global statistics of several hundred Class~0 YSOs that should be found in the whole Gould Belt survey.

preprint2010arXiv

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

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

preprint2010arXiv

The link between molecular cloud structure and turbulence

We aim to better understand how the spatial structure of molecular clouds is governed by turbulence. For that, we study the large-scale spatial distribution of low density molecular gas and search for characteristic length scales. We employ a 35 square degrees 13CO 1-0 molecular line survey of Cygnus X and visual extinction (A_V) maps of 17 Galactic clouds to analyse the spatial structure using the Delta-variance method. This sample contains a large variety of different molecular cloud types with different star forming activity. The Delta-variance spectra obtained from the A_V maps show differences between low-mass star-forming (SF) clouds and massive giant molecular clouds (GMC) in terms of shape of the spectrum and its power-law exponent beta. Low-mass SF clouds have a double-peak structure with characteristic size scales around 1 pc (though with a large scatter around this value) and 4 pc. GMCs show no characteristic scale in the A_V-maps, which can partly be ascribed to a distance effect due to a larger line-of-sight (LOS) confusion. The Delta-variance for Cygnus, determined from the 13CO survey, shows characteristic scales at 4 pc and 40 pc, either reflecting the filament structure and large-scale turbulence forcing or - for the 4 pc scale - the scale below which the 13CO 1-0 line becomes optically thick. Though there are different processes that can introduce characteristic scales, i.e. geometry, decaying turbulence the transition scale from supersonic to subsonic turbulence (the sonic scale), line-of-sight effects and energy injection due to expanding supernova shells, outflows, HII-regions, and the relative contribution of these effects strongly varies from cloud to cloud, it is remarkable that the resulting turbulent structure of molecular clouds shows similar characteristics.

preprint2010arXiv

The Mass Distributions of Starless and Protostellar Cores in Gould Belt Clouds

Using data from the SCUBA Legacy Catalogue (850 um) and Spitzer Space Telescope (3.6 - 70 um), we explore dense cores in the Ophiuchus, Taurus, Perseus, Serpens, and Orion molecular clouds. We develop a new method to discriminate submillimeter cores found by SCUBA as starless or protostellar, using point source photometry from Spitzer wide field surveys. First, we identify infrared sources with red colors associated with embedded young stellar objects (YSOs). Second, we compare the positions of these YSO-candidates to our submillimeter cores. With these identifications, we construct new, self-consistent starless and protostellar core mass functions (CMFs) for the five clouds. We find best fit slopes to the high-mass end of the CMFs of -1.26 +/- 0.20, -1.22 +/- 0.06, -0.95 +/- 0.20, and -1.67 +/- 0.72 for Ophiuchus, Taurus, Perseus, and Orion, respectively. Broadly, these slopes are each consistent with the -1.35 power-law slope of the Salpeter IMF at higher masses, but suggest some differences. We examine a variety of trends between these CMF shapes and their parent cloud properties, potentially finding a correlation between the high-mass slope and core temperature. We also find a trend between core mass and effective size, but we are very limited by sensitivity. We make similar comparisons between core mass and size with visual extinction (for A_V >= 3) and find no obvious trends. We also predict the numbers and mass distributions of cores that future surveys with SCUBA-2 may detect in each of these clouds.

preprint2010arXiv

Variations in H2O+/H2O ratios toward massive star-forming regions

Early results from the Herschel Space Observatory revealed the water cation H2O+ to be an abundant ingredient of the interstellar medium. Here we present new observations of the H2O and H2O+ lines at 1113.3 and 1115.2 GHz using the Herschel Space Observatory toward a sample of high-mass star-forming regions to observationally study the relation between H2O and H2O+ . Nine out of ten sources show absorption from H2O+ in a range of environments: the molecular clumps surrounding the forming and newly formed massive stars, bright high-velocity outflows associated with the massive protostars, and unrelated low-density clouds along the line of sight. Column densities per velocity component of H2 O+ are found in the range of 10^12 to a few 10^13 cm-2 . The highest N(H2O+) column densities are found in the outflows of the sources. The ratios of H2O+/H2O are determined in a range from 0.01 to a few and are found to differ strongly between the observed environments with much lower ratios in the massive (proto)cluster envelopes (0.01-0.1) than in outflows and diffuse clouds. Remarkably, even for source components detected in H2O in emission, H2O+ is still seen in absorption.

preprint2010arXiv

Water cooling of shocks in protostellar outflows: Herschel-PACS map of L1157

In the framework of the Water in Star-forming regions with Herschel (WISH) key program, maps in water lines of several outflows from young stars are being obtained, to study the water production in shocks and its role in the outflow cooling. This paper reports the first results of this program, presenting a PACS map of the o-H2O 179 um transition obtained toward the young outflow L1157. The 179 um map is compared with those of other important shock tracers, and with previous single-pointing ISO, SWAS, and Odin water observations of the same source that allow us to constrain the water abundance and total cooling. Strong H2O peaks are localized on both shocked emission knots and the central source position. The H2O 179 um emission is spatially correlated with emission from H2 rotational lines, excited in shocks leading to a significant enhancement of the water abundance. Water emission peaks along the outflow also correlate with peaks of other shock-produced molecular species, such as SiO and NH3. A strong H2O peak is also observed at the location of the proto-star, where none of the other molecules have significant emission. The absolute 179 um intensity and its intensity ratio to the H2O 557 GHz line previously observed with Odin/SWAS indicate that the water emission originates in warm compact clumps, spatially unresolved by PACS, having a H2O abundance of the order of 10^-4. This testifies that the clumps have been heated for a time long enough to allow the conversion of almost all the available gas-phase oxygen into water. The total water cooling is ~10^-1 Lo, about 40% of the cooling due to H2 and 23% of the total energy released in shocks along the L1157 outflow.

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

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

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