Source author record

N. Schneider

N. Schneider 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

97works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

97 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

FEEDBACK from the NGC7538 HII region

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

preprint2022arXiv

Self-absorption in [CII], $^{12}$CO, and HI in RCW120. Building up a geometrical and physical model of the region

Revealing the 3D dynamics of HII regions and their associated molecular clouds is important for understanding the longstanding problem as to how stellar feedback affects the density structure and kinematics of the interstellar medium. We employed observations of the HII region RCW 120 in [CII], observed within the SOFIA legacy program FEEDBACK, and the $^{12}$CO and $^{13}$CO (3$\to$2) lines, obtained with APEX. In addition we used HI data from the Southern Galactic Plane Survey. Two radiative transfer models were used to fit the observed data. A line profile analysis with the 1D non-LTE radiative transfer code SimLine proves that the CO emission cannot stem from a spherically symmetric molecular cloud configuration. With a two-layer multicomponent model, we then quantified the amount of warm background and cold foreground gas. There is a deficit of CO emission along the line-of-sight toward the center of the HII region which indicates that the HII region is associated with a flattened molecular cloud. Self-absorption in the CO line may hide signatures of infalling and expanding molecular gas. The [CII] emission arises from an expanding [CII] bubble and from the PDRs. A significant part of [CII] emission is absorbed in a cool (~60-100 K), low-density (<500 cm$^{-3}$) atomic foreground layer with a thickness of a few parsec. We propose that the RCW 120 HII region formed in a flattened molecular cloud and is now bursting out of its parental cloud. The compressed surrounding molecular layer formed a torus around the spherically expanding HII bubble. This scenario can possibly be generalized for other HII bubbles and would explain the observed "flat" structure of molecular clouds associated with HII bubbles. We suggest that the [CII] absorption observed in many star-forming regions is at least partly caused by low-density, cool, HI-envelopes surrounding the molecular clouds.

preprint2022arXiv

Star formation in two irradiated globules around Cygnus OB2

We investigate the young stellar populations associated with DR 18 and ECX 6-21, which are two isolated globules irradiated by the O-type stars of the Cygnus OB2 association. Both are HII regions containing obvious tracers of recent and ongoing star formation. We also study smaller isolated molecular structures in their surroundings. Both globules contain their own embedded populations, with a higher fraction of the less-evolved classes. Masses and temperatures are estimated under the assumption of a common age of 1 Myr, which has been found to appropriately represent the general Cygnus OB2 YSO population but is most probably an overestimate for both globules, especially ECX 6-21. The early-B star responsible for the erosion of DR 18 is found to be part of a small aggregate of intermediate-mass stars still embedded in the cloud, which probably contains a second site of recent star formation, also with intermediate-mass stars. We confirm the two main star forming sites embedded in ECX 6-21 described in previous works, with the southern site being more evolved than the northern site. We also discuss the small globule ECX 6-21-W ($= G79.8+1.2$), and propose that its non thermal radio spectrum is due to synchrotron emission from an embedded jet, whose existence is suggested by our observations. The extreme youth of some of the YSOs suggests that star formation in both globules started after they became externally irradiated. The populations of both globules are not found to be particularly rich, but they contain stars with estimated masses similar or above that of the Sun in numbers that hint at some differences with respect to the star formation process taking place in more quiescent regions where low-mass stars dominate, which deeper observations may confirm.

preprint2022arXiv

The SOFIA FEEDBACK Legacy Survey: Dynamics and mass ejection in the bipolar HII region RCW 36

We present [CII] 158 $μ$m and [OI] 63 $μ$m observations of the bipolar HII region RCW 36 in the Vela C molecular cloud, obtained within the SOFIA legacy project FEEDBACK, which is complemented with APEX $^{12/13}$CO(3-2) and Chandra X-ray (0.5-7 keV) data. This shows that the molecular ring, forming the waist of the bipolar nebula, expands with a velocity of 1 - 1.9 km s$^{-1}$. We also observe an increased linewidth in the ring indicating that turbulence is driven by energy injection from the stellar feedback. The bipolar cavity hosts blue-shifted expanding [CII] shells at 5.2$\pm$0.5$\pm$0.5 km s$^{-1}$ (statistical and systematic uncertainty) which indicates that expansion out of the dense gas happens non-uniformly and that the observed bipolar phase might be relatively short ($\sim$0.2 Myr). The X-ray observations show diffuse emission that traces a hot plasma, created by stellar winds, in and around RCW 36. At least 50 \% of the stellar wind energy is missing in RCW 36. This is likely due to leakage which is clearing even larger cavities around the bipolar RCW 36 region. Lastly, the cavities host high-velocity wings in [CII] which indicates relatively high mass ejection rates ($\sim$5$\times$10$^{-4}$ M$_{\odot}$ yr$^{-1}$). This could be driven by stellar winds and/or radiation pressure, but remains difficult to constrain. This local mass ejection, which can remove all mass within 1 pc of RCW 36 in 1-2 Myr, and the large-scale clearing of ambient gas in the Vela C cloud indicates that stellar feedback plays a significant role in suppressing the star formation efficiency (SFE).

preprint2020arXiv

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

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

preprint2020arXiv

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

Dense Gas in a Giant Molecular Filament

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

preprint2020arXiv

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

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

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

preprint2020arXiv

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

The PDR structure and kinematics around the compact HII regions S235A and S235C with [CII], [13CII], [OI] and HCO+ line profiles

The aim of the present work is to study structure and gas kinematics in the photodissociation regions (PDRs) around the compact HII regions S235A and S235C. We observe the [CII], [13CII] and [OI] line emission, using SOFIA/upGREAT and complement them by data of HCO+ and CO. We use the [13CII] line to measure the optical depth of the [CII] emission, and find that the [CII] line profiles are influenced by self-absorption, while the [13CII] line remains unaffected by these effects. Hence, for dense PDRs, [13CII] emission is a better tracer of gas kinematics. The optical depth of the [CII] line is up to 10 in S235A. We find an expanding motion of the [CII]-emitting layer of the PDRs into the front molecular layer in both regions. Comparison of the gas and dust columns shows that gas components visible neither in the [CII] nor in low-J CO lines may contribute to the total column across S235A. We test whether the observed properties of the PDRs match the predictions of spherical models of expanding HII region + PDR + molecular cloud. Integrated intensities of the [13CII], [CII] and [OI] lines are well-represented by the model, but the models do not reproduce the double-peaked [CII] line profiles due to an insufficient column density of C+. The model predicts that the [OI] line could be a more reliable tracer of gas kinematics, but the foreground self-absorbing material does not allow using it in the considered regions.

preprint2020arXiv

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

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

preprint2016arXiv

A census of dense cores in the Taurus L1495 cloud from the Herschel Gould Belt Survey

We present a catalogue of dense cores in a $\sim 4^\circ\times2^\circ$ field of the Taurus star-forming region, inclusive of the L1495 cloud, derived from Herschel SPIRE and PACS observations in the 70 $μ$m, 160 $μ$m, 250 $μ$m, 350 $μ$m, and 500 $μ$m continuum bands. Estimates of mean dust temperature and total mass are derived using modified blackbody fits to the spectral energy distributions. We detect 525 starless cores of which $\sim10$-20% are gravitationally bound and therefore presumably prestellar. Our census of unbound objects is $\sim85$% complete for $M>0.015\,M_\odot$ in low density regions ($A_V\stackrel{<}{_\sim}5$ mag), while the bound (prestellar) subset is $\sim85$% complete for $M>0.1\,M_\odot$ overall. The prestellar core mass function (CMF) is consistent with lognormal form, resembling the stellar system initial mass function, as has been reported previously. All of the inferred prestellar cores lie on filamentary structures whose column densities exceed the expected threshold for filamentary collapse, in agreement with previous reports. Unlike the prestellar CMF, the unbound starless CMF is not lognormal, but instead is consistent with a power-law form below $0.3\,M_\odot$ and shows no evidence for a low-mass turnover. It resembles previously reported mass distributions for CO clumps at low masses ($M\stackrel{<}{_\sim}0.3\,M_\odot$). The volume density PDF, however, is accurately lognormal except at high densities. It is consistent with the effects of self-gravity on magnetized supersonic turbulence. The only significant deviation from lognormality is a high-density tail which can be attributed unambiguously to prestellar cores.

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

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

Hi-GAL, the Herschel infrared Galactic Plane Survey: photometric maps and compact source catalogues. First data release for Inner Milky Way: +68°> l > -70°

(Abridged) We present the first public release of high-quality data products (DR1) from Hi-GAL, the {\em Herschel} infrared Galactic Plane Survey. Hi-GAL is the keystone of a suite of continuum Galactic Plane surveys from the near-IR to the radio, and covers five wavebands at 70, 160, 250, 350 and 500 micron, encompassing the peak of the spectral energy distribution of cold dust for 8 < T < 50K. This first Hi-GAL data release covers the inner Milky Way in the longitude range 68° > l > -70° in a |b|<1° latitude strip. Photometric maps have been produced with the ROMAGAL pipeline, that optimally capitalizes on the excellent sensitivity and stability of the bolometer arrays of the {\em Herschel} PACS and SPIRE photometric cameras, to deliver images of exquisite quality and dynamical range, absolutely calibrated with {\em Planck} and {\em IRAS}, and recovering extended emission at all wavelengths and all spatial scales. The compact source catalogues have been generated with the CuTEx algorithm, specifically developed to optimize source detection and extraction in the extreme conditions of intense and spatially varying background that are found in the Galactic Plane in the thermal infrared. Hi-GAL DR1 images will be accessible via a dedicated web-based image cutout service. The DR1 Compact Source Catalogues are delivered as single-band photometric lists containing, in addition to source position, peak and integrated flux and source sizes, a variety of parameters useful to assess the quality and reliability of the extracted sources, caveats and hints to help this assessment are provided. Flux completeness limits in all bands are determined from extensive synthetic source experiments and depend on the specific line of sight along the Galactic Plane. Hi-GAL DR1 catalogues contain 123210, 308509, 280685, 160972 and 85460 compact sources in the five bands, respectively.

preprint2016arXiv

Photoevaporation and close encounters: how the environment around Cygnus OB2 affects the evolution of protoplanetary disks

In our Galaxy, star formation occurs in a variety of environments, with a large fraction of stars formed in clusters hosting massive stars. OB stars have an important feedback on the evolution of protoplanetary disks around nearby young stars and likely on the process of planet formation occurring in them. The nearby massive association Cygnus OB2 is an outstanding laboratory to study this feedback. It is the closest massive association to our Sun, and hosts hundreds of massive stars and thousands of low mass members. In this paper, we analyze the spatial variation of the disk fraction in Cygnus OB2 and we study its correlation with the local values of Far and Extreme ultraviolet radiation fields and the local stellar surface density. We present definitive evidence that disks are more rapidly dissipated in the regions of the association characterized by intense local UV field and large stellar density. In particular, the FUV radiation dominates disks dissipation timescales in the proximity (i.e. within 0.5 pc) of the O stars. In the rest of the association, EUV photons potentially induce a significant mass loss from the irradiated disks across the entire association, but the efficiency of this process is reduced at increasing distances from the massive stars due to absorption by the intervening intracluster material. We find that disk dissipation due to close stellar encounters is negligible in Cygnus OB2, and likely to have affected 1% or fewer of the stellar population. Disk dissipation is instead dominated by photoevaporation. We also compare our results to what has been found in other young clusters with different massive populations, concluding that massive associations like Cygnus OB2 are potentially hostile to protoplanetary disks, but that the environments where disks can safely evolve in planetary systems are likely quite common in our Galaxy.

preprint2016arXiv

Tracing extended low-velocity shocks through SiO emission - Case study of the W43-MM1 ridge

Previous literature suggests that the densest structures in the interstellar medium form through colliding flows, but patent evidence of this process is still missing. Recent literature proposes using SiO line emission to trace low-velocity shocks associated with cloud formation through collision. In this paper we investigate the bright and extended SiO(2-1) emission observed along the ~5 pc-long W43-MM1 ridge to determine its origin. We used high angular resolution images of the SiO(2-1) and HCN(1-0) emission lines obtained with the IRAM plateau de Bure (PdBI) interferometer and combined with data from the IRAM 30 m radiotelescope. These data were complemented by a Herschel column density map of the region. We performed spectral analysis of SiO and HCN emission line profiles to identify protostellar outflows and spatially disentangle two velocity components associated with low- and high-velocity shocks. Then, we compared the low-velocity shock component to a dedicated grid of one-dimensional (1D) radiative shock models. We find that the SiO emission originates from a mixture of high-velocity shocks caused by bipolar outflows and low-velocity shocks. Using SiO and HCN emission lines, we extract seven bipolar outflows associated with massive dense cores previously identified within the W43-MM1 mini-starburst cluster. Comparing observations with dedicated Paris-Durham shock models constrains the velocity of the low-velocity shock component from 7 to 12km/s. The SiO arising from low-velocity shocks spreads along the complete length of the ridge. Its contribution represents at least 45% and up to 100% of the total SiO emission depending on the area considered. The low-velocity component of SiO is most likely associated with the ridge formation through colliding flows or cloud-cloud collision.

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

A possible link between the power spectrum of interstellar filaments and the origin of the prestellar core mass function

Two major features of the prestellar CMF are: 1) a broad peak below 1 Msun, presumably corresponding to a mean gravitational fragmentation scale, and 2) a characteristic power-law slope, very similar to the Salpeter slope of the stellar initial mass function (IMF) at the high-mass end. While recent Herschel observations have shown that the peak of the prestellar CMF is close to the thermal Jeans mass in marginally supercritical filaments, the origin of the power-law tail of the CMF/IMF at the high-mass end is less clear. Inutsuka (2001) proposed a theoretical scenario in which the origin of the power-law tail can be understood as resulting from the growth of an initial spectrum of density perturbations seeded along the long axis of filaments by interstellar turbulence. Here, we report the statistical properties of the line-mass fluctuations of filaments in nearby molecular clouds observed with Herschel using a 1-D power spectrum analysis. The observed filament power spectra were fitted by a power-law function $(P_{true}(s) \propto s^α)$ after removing the effect of beam convolution at small scales. A Gaussian-like distribution of power-spectrum slopes was found centered at -1.6, close to that of the one-dimensional velocity power spectrum generated by subsonic Kolomogorov turbulence (-1.67). An empirical correlation, $P^{0.5}(s_0) \propto <N_{\rm H_2}>^{1.4 \pm 0.1} $, was also found between the amplitude of each filament power spectrum $P(s_0)$ and the mean column density along the filament $<N_{\rm H_2}>$. Finally, the dispersion of line-mass fluctuations along each filament $σ_{\rm M_{line}}$ was found to scale with the physical length $L$ of the filament, roughly as $σ_{M_{line}} \propto L^{0.7}$. Overall, our results are consistent with the suggestion that the bulk of the CMF/IMF results from the gravitational fragmentation of filaments.

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

Filaments in the Lupus molecular clouds

We have studied the filaments extracted from the column density maps of the nearby Lupus 1, 3, and 4 molecular clouds, derived from photometric maps observed with the Herschel satellite. Filaments in the Lupus clouds have quite low column densities, with a median value of $\sim$1.5$\times$10$^{21}$ cm$^{-2}$ and most have masses per unit length lower than the maximum critical value for radial gravitational collapse. Indeed, no evidence of filament contraction has been seen in the gas kinematics. We find that some filaments, that on average are thermally subcritical, contain dense cores that may eventually form stars. This is an indication that in the low column density regime, the critical condition for the formation of stars may be reached only locally and this condition is not a global property of the filament. Finally, in Lupus we find multiple observational evidences of the key role that the magnetic field plays in forming filaments, and determining their confinement and dynamical evolution.

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

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

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

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

A Herschel and BIMA study of the sequential star formation near the W48A HII region

We present the results of Herschel HOBYS photometric mapping combined with BIMA observations and additional archival data, and perform an in-depth study of the evolutionary phases of the star-forming clumps in W 48A and their surroundings. Age estimates for the compact sources were derived from bolometric luminosities and envelope masses, which were obtained from the dust continuum emission, and agree within an order of magnitude with age estimates from molecular line and radio data. The clumps in W 48A are linearly aligned by age (east-old to west-young): we find a ultra compact (UC) HII region, a young stellar object (YSO) with class II methanol maser emission, a YSO with a massive outflow, and finally the NH_2D prestellar cores from Pillai et al. This remarkable positioning reflects the (star) formation history of the region. We find that it is unlikely that the star formation in the W 48A molecular cloud was triggered by the UCHII region and discuss the Aquila supershell expansion as a mayor influence on the evolution of W 48A. We conclude that the combination of Herschel continuum data with interferometric molecular line and radio continuum data is important to derive trustworthy age estimates and interpret the origin of large scale structures through kinematic information.

preprint2014arXiv

Age, size, and position of H ii regions in the Galaxy. Expansion of ionized gas in turbulent molecular clouds

This work aims at improving the current understanding of the interaction between H ii regions and turbulent molecular clouds. We propose a new method to determine the age of a large sample of OB associations by investigating the development of their associated H ii regions in the surrounding turbulent medium. Using analytical solutions, one-dimensional (1D), and three-dimensional (3D) simulations, we constrained the expansion of the ionized bubble depending on the turbulent level of the parent molecular cloud. A grid of 1D simulations was then computed in order to build isochrone curves for H ii regions in a pressure-size diagram. This grid of models allowed to date large sample of OB associations and was used on the H ii Region Discovery Survey (HRDS). Analytical solutions and numerical simulations showed that the expansion of H ii regions is slowed down by the turbulence up to the point where the pressure of the ionized gas is in a quasi-equilibrium with the turbulent ram pressure. Based on this result, we built a grid of 1D models of the expansion of H ii regions in a profile based on Larson laws. The 3D turbulence is taken into account by an effective 1D temperature profile. The ages estimated by the isochrones of this grid agree well with literature values of well-known regions such as Rosette, RCW 36, RCW 79, and M16. We thus propose that this method can be used to give ages of young OB associations through the Galaxy such as the HRDS survey and also in nearby extra-galactic sources.

preprint2014arXiv

Characterizing the structure of diffuse emission in Hi-GAL maps

We present a study of the structure of the Galactic interstellar medium through the Delta-variance technique, related to the power spectrum and the fractal properties of infrared/sub-mm maps. Through this method, it is possible to provide quantitative parameters which are useful to characterize different morphological and physical conditions, and to better constrain the theoretical models. In this respect, the Herschel Infrared Galactic Plane Survey carried out at five photometric bands from 70 to 500 μm constitutes an unique database for applying statistical tools to a variety of regions across the Milky Way. In this paper, we derive a robust estimate of the power-law portion of the power spectrum of four contiguous 2°x2° Hi-GAL tiles located in the third Galactic quadrant (217° < l < 225°, -2° < b < 0°). The low level of confusion along the line of sight testified by CO observations makes this region an ideal case. We find very different values of the power spectrum slope from tile to tile but also from wavelength to wavelength (2 < β< 3), with similarities between fields attributable to components located at the same distance. Thanks to the comparison with models of turbulence, an explanation of the determined slopes in terms of the fractal geometry is also provided, and possible relations with the underlying physics are investigated. In particular, an anti-correlation between ISM fractal dimension and star formation efficiency is found for the two main distance components observed in these fields. A possible link between the fractal properties of the diffuse emission and the resulting clump mass function is discussed.

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

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

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

Reconstructing the density and temperature structure of prestellar cores from $Herschel$ data: A case study for B68 and L1689B

Utilizing multi-wavelength dust emission maps acquired with $Herschel$, we reconstruct local volume density and dust temperature profiles for the prestellar cores B68 and L1689B using inverse-Abel transform based technique. We present intrinsic radial dust temperature profiles of starless cores directly from dust continuum emission maps disentangling the effect of temperature variations along the line of sight which was previously limited to the radiative transfer calculations. The reconstructed dust temperature profiles show a significant drop in core center, a flat inner part, and a rising outward trend until the background cloud temperature is reached. The central beam-averaged dust temperatures obtained for B68 and L1689B are 9.3 $\pm$ 0.5 K and 9.8 $\pm$0.5 K, respectively, which are lower than the temperatures of 11.3 K and 11.6 K obtained from direct SED fitting. The best mass estimates derived by integrating the volume density profiles of B68 and L1689B are 1.6 M_sol and 11 M_sol, respectively. Comparing our results for B68 with the near-infrared extinction studies, we find that the dust opacity law adopted by the HGBS project, $κ_λ =0.1(λ/300 μm)^{-2}$, agrees to within 50% with the dust extinction constraints

preprint2013arXiv

The first Hi-GAL observations of the outer Galaxy: a look to star formation in the third Galactic quadrant in the longitude range 216.5 < l < 225.5

We present the first Herschel PACS and SPIRE photometric observations in a portion of the outer Galaxy ($216.5^{\circ} \lesssim \ell \lesssim 225.5^{\circ}$ and $-2^{\circ} \lesssim b \lesssim 0^{\circ}$) as a part of the Hi-GAL survey. The maps between 70 and 500 $μ$m, the derived column density and temperature maps, and the compact source catalog are presented. NANTEN CO(1-0) line observations are used to derive cloud kinematics and distances, so that we can estimate distance-dependent physical parameters of the compact sources (cores and clumps) having a reliable spectral energy distribution, that we separate in 255 proto-stellar and 688 starless. Both typologies are found in association with all the distance components observed in the field, up to $\sim 5.8$ kpc, testifying the presence of star formation beyond the Perseus arm at these longitudes. Selecting the starless gravitationally bound sources we identify 590 pre-stellar candidates. Several sources of both proto- and pre-stellar nature are found to exceed the minimum requirement for being compatible with massive star formation, based on the mass-radius relation. For the pre-stellar sources belonging to the Local arm ($d\lesssim1.5$ kpc) we study the mass function, whose high-mass end shows a power-law $N(\log M) \propto M^{-1.0 \pm 0.2}$. Finally, we use a luminosity vs mass diagram to infer the evolutionary status of the sources, finding that most of the proto-stellar are in the early accretion phase (with some cases compatible with a Class I stage), while for pre-stellar sources, in general, accretion has not started yet.

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.

preprint2013arXiv

Young stellar clusters in the Rosette molecular cloud. Arguments against triggered star formation

The Rosette complex is a well studied region of the galactic plane which presents the apparent characteristics of a triggered star forming region. This is however still debated as no strong evidence corroborates this statement. We focus on characterizing the young stellar population in the Rosette complex to improve our understanding of the processes that regulate the star formation in this region. We propose an original method that relies on the joint analysis of the star color and density in the near-infrared. It leads to mapping the molecular cloud spatial distribution and detecting the embedded clusters with their characterization in terms of member number and age estimation. We have identified 13 clusters, 2 of which are new discoveries, and we estimate that the total number of young stellar objects in the Rosette ranges between 4000 and 8000 members. We find that the age distribution of the young clusters is not consistent with a general triggered scenario for the star formation in this molecular cloud. This study proves that the Rosette complex evolution is not governed by the influence of its OB star population. It suggests that the simple morphological appearance of an active region is not sufficient to conclude much about the triggering role in the star formation process. Our method of constraining the cluster properties using UKIDSS and WISE data has proven efficient, and studies of other regions of the galactic plane would definitely benefit from this approach.

preprint2012arXiv

3D simulations of globules and pillars formation around HII regions: turbulence and shock curvature

We investigate the interplay between the ionization radiation from massive stars and the turbulence inside the surrounding molecular gas thanks to 3D numerical simulations. We used the 3D hydrodynamical code HERACLES to model an initial turbulent medium that is ionized and heated by an ionizing source. Three different simulations are performed with different mean Mach numbers (1, 2 and 4). A non-equilibrium model for the ionization and the associated thermal processes was used. This revealed to be crucial when turbulent ram pressure is of the same order as the ionized-gas pressure. The density structures initiated by the turbulence cause local curvatures of the dense shell formed by the ionization compression. When the curvature of the shell is sufficient, the shell collapse on itself to form a pillar while a smaller curvature leads to the formation of dense clumps that are accelerated with the shell and therefore remain in the shell during the simulation. When the turbulent ram pressure of the cold gas is sufficient to balance the ionized-gas pressure, some dense-gas bubbles have enough kinetic energy to penetrate inside the ionized medium, forming cometary globules. This suggests a direct relation in the observations between the presence of globules and the relative importance of the turbulence compared to the ionized-gas pressure. The probability density functions present a double peak structure when the turbulence is low relative to the ionized-gas pressure. This could be used in observations as an indication of the turbulence inside molecular clouds.

preprint2012arXiv

A multi-scale, multi-wavelength source extraction method: getsources

We present a multi-scale, multi-wavelength source extraction algorithm called getsources. Although it has been designed primarily for use in the far-infrared surveys of Galactic star-forming regions with Herschel, the method can be applied to many other astronomical images. Instead of the traditional approach of extracting sources in the observed images, the new method analyzes fine spatial decompositions of original images across a wide range of scales and across all wavebands. It cleans those single-scale images of noise and background, and constructs wavelength-independent single-scale detection images that preserve information in both spatial and wavelength dimensions. Sources are detected in the combined detection images by following the evolution of their segmentation masks across all spatial scales. Measurements of the source properties are done in the original background-subtracted images at each wavelength; the background is estimated by interpolation under the source footprints and overlapping sources are deblended in an iterative procedure. In addition to the main catalog of sources, various catalogs and images are produced that aid scientific exploitation of the extraction results. We illustrate the performance of getsources on Herschel images by extracting sources in sub-fields of the Aquila and Rosette star-forming regions. The source extraction code and validation images with a reference extraction catalog are freely available.

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 observations of B1-bS and B1-bN: two first hydrostatic core candidates in the Perseus star-forming cloud

We report far-IR Herschel observations obtained between 70 $μ$m and 500 $μ$m of two star-forming dusty condensations, B1-bS and B1-bN, in the B1 region of the Perseus star-forming cloud. In the Western part of the Perseus cloud, B1-bS is the only source detected in all of the 6 PACS and SPIRE photometric bands without being visible in the Spitzer map at 24 $μ$m. B1-bN is clearly detected between 100 $μ$m and 250 $μ$m. We have fitted the spectral energy distributions of these sources to derive their physical properties, and find that a simple greybody model fails to reproduce the observed SEDs. At least a two-component model, consisting of a central source surrounded by a dusty envelope is required. The properties derived from the fit, however, suggest that the central source is not a Class 0 object. We then conclude that while B1-bS and B1-bN appear to be more evolved than a pre-stellar core, the best-fit models suggest that their central objects are younger than a Class 0 source. Hence, they may be good candidates to be examples of the first hydrostatic core phase. The projected distance between B1-bS and B1-bN is a few Jeans lengths. If their physical separation is close to this value, this pair would allow the mutual interactions between two forming stars at a very early stage of their evolution to be studied.

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

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

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 Pipe Nebula as seen with Herschel: Formation of filamentary structures by large-scale compression ?

A growing body of evidence indicates that the formation of filaments in interstellar clouds is a key component of the star formation process. In this paper, we present new Herschel PACS and SPIRE observations of the B59 and Stem regions in the Pipe Nebula complex, revealing a rich, organized network of filaments. The asymmetric column density profiles observed for several filaments, along with the bow-like edge of B59, indicates that the Pipe Nebula is being compressed from its western side, most likely by the winds from the nearby Sco OB2 association. We suggest that this compressive flow has contributed to the formation of some of the observed filamentary structures. In B59, the only region of the entire Pipe complex showing star formation activity, the same compressive flow has likely enhanced the initial column density of the clump, allowing it to become globally gravitationally unstable. Although more speculative, we propose that gravity has also been responsible for shaping the converging filamentary pattern observed in B59. While the question of the relative impact of large-scale compression and gravity remains open in B59, large-scale compression appears to be a plausible mechanism for the initial formation of filamentary structures in the rest of the complex

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.

preprint2012arXiv

Worldwide site comparison for submillimetre astronomy

The most important limitation for ground-based submillimetre (submm) astronomy is the broad-band absorption of the total water vapour in the atmosphere above an observation site, often expressed as the Precipitable Water Vapour (PWV). A long-term statistic on the PWV is thus mandatory to characterize the quality of an existing or potential site for observational submm-astronomy. In this study we present a three-year statistic (2008-2010) of the PWV for ground-based telescope sites all around the world and for stratospheric altitudes relevant for SOFIA (Stratospheric Observatory for Far-infrared astronomy). The submm-transmission is calculated for typical PWVs using an atmospheric model. We present the absolute PWV values for each site sorted by year and time percentage. The PWV corresponding to the first decile (10%) and the quartiles (25%, 50%, 75%) are calculated and transmission curves between 150 μm and 3 mm for these values are shown. The Antarctic and South-American sites present very good conditions for submillimetre astronomy. The 350 μm and 450 μm atmospheric windows are open all year long whereas the 200 μm atmospheric window opens reasonably for 25 % of the time in Antarctica and the extremely high-altitude sites in Chile. Potential interesting new facilities are Macon in Argentinia and Summit in Greenland that show similar conditions as for example Mauna Kea (Hawaii). For SOFIA, we present in more detail transmission curves for different altitudes (11 to 14 km), PWV values, and higher frequencies (up to 5 THz). Though the atmosphere at these altitude is generally very transparent, the absorption at very high frequencies becomes more important, partly caused by minor species. In conclusion, the method presented in this paper could identify sites on Earth with a great potential for submillimetre astronomy, and guide future site testing campaigns in situ.

preprint2011arXiv

3D simulations of pillars formation around HII regions: the importance of shock curvature

Radiative feedback from massive stars is a key process to understand how HII regions may enhance or inhibit star formation in pillars and globules at the interface with molecular clouds. We aim to contribute to model the interactions between ionization and gas clouds to better understand the processes at work. We study in detail the impact of modulations on the cloud-HII region interface and density modulations inside the cloud. We run three-dimensional hydrodynamical simulations based on Euler equations coupled with gravity using the HERACLES code. We implement a method to solve ionization/recombination equations and we take into account typical heating and cooling processes at work in the interstellar medium and due to ionization/recombination physics. UV radiation creates a dense shell compressed between an ionization front and a shock ahead. Interface modulations produce a curved shock that collapses on itself leading to stable growing pillar-like structures. The narrower the initial interface modulation, the longer the resulting pillar. We interpret pillars resulting from density modulations in terms of the ability of these density modula- tions to curve the shock ahead the ionization front. The shock curvature is a key process to understand the formation of structures at the edge of HII regions. Interface and density modulations at the edge of the cloud have a direct impact on the morphology of the dense shell during its formation. Deeper in the cloud, structures have less influence due to the high densities reached by the shell during its expansion.

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

Site testing for submillimetre astronomy at Dome C, Antarctica

Over the past few years a major effort has been put into the exploration of potential sites for the deployment of submillimetre astronomical facilities. Amongst the most important sites are Dome C and Dome A on the Antarctic Plateau, and the Chajnantor area in Chile. In this context, we report on measurements of the sky opacity at 200 um over a period of three years at the French-Italian station, Concordia, at Dome C, Antarctica. We also present some solutions to the challenges of operating in the harsh polar environ- ment. Dome C offers exceptional conditions in terms of absolute atmospheric transmission and stability for submillimetre astron- omy. Over the austral winter the PWV exhibits long periods during which it is stable and at a very low level (0.1 to 0.3 mm). Higher values (0.2 to 0.8 mm) of PWV are observed during the short summer period. Based on observations over three years, a transmission of around 50% at 350 um is achieved for 75% of the time. The 200-um window opens with a typical transmission of 10% to 15% for 25% of the time. Dome C is one of the best accessible sites on Earth for submillimetre astronomy. Observations at 350 or 450 μm are possible all year round, and the 200-um window opens long enough and with a sufficient transparency to be useful. Although the polar environment severely constrains hardware design, a permanent observatory with appropriate technical capabilities is feasible. Because of the very good astronomical conditions, high angular resolution and time series (multi-year) observations at Dome C with a medium size single dish telescope would enable unique studies to be conducted, some of which are not otherwise feasible even from space.

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.

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

Discovery of an extremely bright submillimeter galaxy at z=3.93

Serendipitously we have discovered a rare, bright submillimeter galaxy (SMG) with a flux density of 30 +/- 2 mJy at lambda=1.2mm, using MAMBO2 at the IRAM 30-meter millimeter telescope. Although no optical counterpart is known for MM18423+5938, we were able to measure the redshift z=3.92960 +/- 0.00013 from the detection of CO lines using the IRAM Eight MIxer Receiver (EMIR). In addition, by collecting all available photometric data in the far-infrared and radio to constrain its spectral energy distribution, we derive the FIR luminosity 4.8 10^14/m Lsol and mass 6.0 10^9/m Msol for its dust, allowing for a magnification factor m caused by a probable gravitational lens. The corresponding star-formation rate is 8.3 10^4/m Msol/yr. The detection of three lines of the CO rotational ladder, and a significant upper limit for a fourth CO line, allow us to estimate an H2 mass of between 1.9 10^11/m Msol and 1.1 10^12/m Msol. The two lines CI(3p1-3p0) and CI(3p2-3p1) were clearly detected and yield a [CI]/[H2] number abundance between 1.4 10^-5 and 8.0 10^-5. Upper limits are presented for emission lines of HCN, HCO^+, HNC, H_2O and other molecules observed. The moderate excitation of the CO lines is indicative of an extended starburst, and excludes the dominance of an AGN in heating this high-redshift SMG.

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-SPIRE observations of the Polaris flare : structure of the diffuse interstellar medium at the sub-parsec scale

We present a power spectrum analysis of the Herschel-SPIRE observations of the Polaris flare, a high Galactic latitude cirrus cloud midway between the diffuse and molecular phases. The SPIRE images of the Polaris flare reveal for the first time the structure of the diffuse interstellar medium down to 0.01 parsec over a 10 square degrees region. These exceptional observations highlight the highly filamentary and clumpy structure of the interstellar medium even in diffuse regions of the map. The power spectrum analysis shows that the structure of the interstellar medium is well described by a single power law with an exponent of -2.7 +- 0.1 at all scales from 30" to 8 degrees. That the power spectrum slope of the dust emission is constant down to the SPIRE angular resolution is an indication that the inertial range of turbulence extends down to the 0.01 pc scale. The power spectrum analysis also allows the identification of a Poissonian component at sub-arcminute scales in agreement with predictions of the cosmic infrared background level at SPIRE wavelengths. Finally, the comparison of the SPIRE and IRAS 100 micron data of the Polaris flare clearly assesses the capability of SPIRE in maping diffuse emission over large areas.

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

HIFI observations of warm gas in DR21: Shock versus radiative heating

The molecular gas in the DR21 massive star formation region is known to be affected by the strong UV field from the central star cluster and by a fast outflow creating a bright shock. The relative contribution of both heating mechanisms is the matter of a long debate. By better sampling the excitation ladder of various tracers we provide a quantitative distinction between the different heating mechanisms. HIFI observations of mid-J transitions of CO and HCO+ isotopes allow us to bridge the gap in excitation energies between observations from the ground, characterizing the cooler gas, and existing ISO LWS spectra, constraining the properties of the hot gas. Comparing the detailed line profiles allows to identify the physical structure of the different components. In spite of the known shock-excitation of H2 and the clearly visible strong outflow, we find that the emission of all lines up to > 2 THz can be explained by purely radiative heating of the material. However, the new Herschel/HIFI observations reveal two types of excitation conditions. We find hot and dense clumps close to the central cluster, probably dynamically affected by the outflow, and a more widespread distribution of cooler, but nevertheless dense, molecular clumps.

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

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 Herschel-SPIRE instrument and its in-flight performance

The Spectral and Photometric Imaging Receiver (SPIRE), is the Herschel Space Observatory`s submillimetre camera and spectrometer. It contains a three-band imaging photometer operating at 250, 350 and 500 microns, and an imaging Fourier Transform Spectrometer (FTS) which covers simultaneously its whole operating range of 194-671 microns (447-1550 GHz). The SPIRE detectors are arrays of feedhorn-coupled bolometers cooled to 0.3 K. The photometer has a field of view of 4' x 8', observed simultaneously in the three spectral bands. Its main operating mode is scan-mapping, whereby the field of view is scanned across the sky to achieve full spatial sampling and to cover large areas if desired. The spectrometer has an approximately circular field of view with a diameter of 2.6'. The spectral resolution can be adjusted between 1.2 and 25 GHz by changing the stroke length of the FTS scan mirror. Its main operating mode involves a fixed telescope pointing with multiple scans of the FTS mirror to acquire spectral data. For extended source measurements, multiple position offsets are implemented by means of an internal beam steering mirror to achieve the desired spatial sampling and by rastering of the telescope pointing to map areas larger than the field of view. The SPIRE instrument consists of a cold focal plane unit located inside the Herschel cryostat and warm electronics units, located on the spacecraft Service Module, for instrument control and data handling. Science data are transmitted to Earth with no on-board data compression, and processed by automatic pipelines to produce calibrated science products. The in-flight performance of the instrument matches or exceeds predictions based on pre-launch testing and modelling: the photometer sensitivity is comparable to or slightly better than estimated pre-launch, and the spectrometer sensitivity is also better by a factor of 1.5-2.

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.

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.