Source author record

P. Schilke

P. Schilke 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

70works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

70 published item(s)

preprint2025arXiv

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

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

preprint2022arXiv

HyGAL: Characterizing the Galactic ISM with observations of hydrides and other small molecules -- I. Survey description and a first look toward W3(OH), W3 IRS5 and NGC 7538 IRS1

The HyGAL SOFIA legacy program surveys six hydride molecules -- ArH+, OH+, H2O+, SH, OH, and CH -- and two atomic constituents -- C+ and O -- within the diffuse interstellar medium (ISM) by means of absorption-line spectroscopy toward 25 bright Galactic background continuum sources. This detailed spectroscopic study is designed to exploit the unique value of specific hydrides as tracers and probes of different phases of the ISM, as demonstrated by recent studies with the Herschel Space Observatory. The observations performed under the HyGAL program will allow us to address several questions related to the lifecycle of molecular material in the ISM and the physical processes that impact its phase transition, such as: (1) What is the distribution function of the H2 fraction in the ISM? (2) How does the ionization rate due to low-energy cosmic-rays vary within the Galaxy? (3) What is the nature of interstellar turbulence, and what mechanisms lead to its dissipation? This overview discusses the observing strategy, synergies with ancillary and archival observations, the data reduction and analysis schemes adopted; and presents the first results obtained toward three of the survey targets, W3(OH), W3IRS5 and NGC7538IRS1. Robust measurements of the column densities of these hydrides -- obtained through widespread observations of absorption lines-- help address the questions raised, and there is a timely synergy between these observations and the development of theoretical models, particularly pertaining to the formation of H2 within the turbulent ISM. The provision of enhanced HyGAL data products will therefore serve as a legacy for future ISM studies.

preprint2022arXiv

Importance of source structure on complex organics emission. I. Observations of CH$_3$OH from low-mass to high-mass protostars

Complex organic molecules (COMs) are often observed toward embedded Class 0 and I protostars. However, not all Class 0 and I protostars exhibit COMs emission. In this work, variations in methanol (CH$_3$OH) emission are studied to test if absence of CH$_3$OH emission can be linked to source properties. Combining both new and archival observations with ALMA and sources from the literature, a sample of 184 low-mass and high-mass protostars is investigated. The warm (T > 100 K) gaseous CH$_3$OH mass, $M_{\rm CH_3OH}$, is determined for each source using primarily optically thin isotopologues. On average, Class I protostellar systems seem to have less warm $M_{\rm CH_3OH}$ ($<10^{-10}$ M$_\odot$) than younger Class 0 sources ($\sim10^{-7}$ M$_\odot$). High-mass sources in our sample show higher warm $M_{\rm CH_3OH}$ up to $10^{-7}-10^{-3}$ M$_\odot$. To take into account the effect of the source's overall mass on $M_{\rm CH_3OH}$, a normalized CH$_3$OH mass is defined as $M_{\rm CH_3OH}/M_{\rm dust,0}$, where $M_{\rm dust,0}$ is the cold + warm dust mass within a fixed radius. Excluding upper limits, a simple power-law fit to the normalized warm CH$_3$OH masses results in $M_{\rm CH_3OH}/M_{\rm dust,0}\propto L_{\rm bol}^{0.70\pm0.05}$. This is in good agreement with a simple hot core toy model which predicts that the normalized $M_{\rm CH_3OH}$ increases with $L_{\rm bol}^{0.75}$ due to the snowline moving outward. Sources for which the size of the disk is equivalent or smaller than the estimated 100 K radius agree well with the best-fit power-law model, whereas sources with significantly larger disks show up to two orders of magnitude lower normalized warm CH$_3$OH masses. Based on the latter results, we suggest that source structure such as a disk can result in colder gas and thus fewer COMs in the gas phase. Additionally, optically thick dust can hide the emission of COMs.

preprint2022arXiv

The sharp ALMA view of infall and outflow in the massive protocluster G31.41+0.31

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

preprint2020arXiv

A search for radio jets from massive young stellar objects. Association of radio jets with H2O and CH3OH masers

Recent theoretical and observational studies debate the similarities between the formation process of high-mass (>8 Msun) and low-mass stars. The formation of low-mass star formation is directly associated with the presence of disks and jets. According to this scenario, radio jets are expected to be common in high-mass star-forming regions. We aim to increase the number of known radio jets in high-mass star forming regions by searching for radio jet candidates at radio continuum wavelengths. We have used the Karl G. Jansky Very Large Array (VLA) to observe 18 high-mass star-forming regions in the C band (6 cm, ~1.0 arcsec resolution) and K band (1.3 cm, ~0.3 arcsec resolution). We have searched for radio jet candidates by studying the association of radio continuum sources with shock activity signposts. We have identified 7 as the most probable radio jets. The radio luminosity of the radio jet candidates is correlated with the bolometric luminosity and the outflow momentum rate. About 7-36% of the radio jet candidates are associated with non-thermal emission. The radio jet candidates associated with 6.7 GHz CH3OH maser emission are preferentially thermal winds and jets, while a considerable fraction of radio jet candidates associated with H2O masers show non-thermal emission, likely due to strong shocks. Our sample of 18 regions is divided in 8 less evolved, infrared-dark regions and 10 more evolved, infrared-bright regions. We have found that ~71% of the identified radio jet candidates are located in the more evolved regions. Similarly, 25% of the less evolved regions harbor one of the most probable radio jets, while up to 50% of the more evolved regions contain one of these radio jet candidates. This suggests that the detection of radio jets in high-mass star forming regions is larger in slightly more evolved regions.

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

The GUAPOS project: G31.41+0.31 Unbiased ALMA sPectral Observational Survey -- I. Isomers of C$_{2}$H$_{4}$O$_{2}$

Understanding the degree of chemical complexity that can be reached in star-forming regions, together with the identification of precursors of the building blocks of life in the interstellar medium, is one of the goals of astrochemistry. Unbiased spectral surveys with large bandwidth and high spectral resolution are thus needed, to resolve line blending in chemically rich sources and identify complex organic molecules. This kind of observations has been successfully carried out, mainly towards the Galactic Center, a region that shows peculiar environmental conditions. We present an unbiased spectral survey at 3mm of one of the most chemically rich hot molecular cores located outside the Galactic Center, in the high-mass star-forming region G31.41+0.31. In this first paper, we present the survey and discuss the detection of the 3 isomers of C$_{2}$H$_{4}$O$_{2}$: methyl formate, glycolaldehyde and acetic acid. Observations were carried out with ALMA and cover the entire Band 3 from 84 to 116 GHz with an angular resolution of $1.2^{''}$x$1.2^{''}$ and a spectral resolution of $\sim0.488$ MHz. The transitions of the 3 molecules have been analyzed with the software XCLASS. All three isomers were detected and methyl formate and acetic acid abundances in G31 are the highest detected up to now, if compared to sources in literature. The size of the emission varies among the three isomers with acetic acid showing the most compact emission while methyl formate the most extended. The comparison with chemical models suggests the necessity of grain-surface routes for the formation of methyl formate in G31, while for glycolaldehyde both scenarios could be feasible. Proposed grain-surface reaction for acetic acid is not able to reproduce the observed abundance in this work, while gas-phase scenario should be further tested due to large uncertainties.

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

Molecular gas kinematics within the central 250 pc of the Milky Way

Using spectral-line observations of HNCO, N2H+, and HNC, we investigate the kinematics of dense gas in the central ~250 pc of the Galaxy. We present SCOUSE (Semi-automated multi-COmponent Universal Spectral-line fitting Engine), a line fitting algorithm designed to analyse large volumes of spectral-line data efficiently and systematically. Unlike techniques which do not account for complex line profiles, SCOUSE accurately describes the {l, b, v_LSR} distribution of CMZ gas, which is asymmetric about Sgr A* in both position and velocity. Velocity dispersions range from 2.6 km/s<σ<53.1 km/s. A median dispersion of 9.8 km/s, translates to a Mach number, M_3D>28. The gas is distributed throughout several "streams", with projected lengths ~100-250 pc. We link the streams to individual clouds and sub-regions, including Sgr C, the 20 and 50 km/s clouds, the dust ridge, and Sgr B2. Shell-like emission features can be explained by the projection of independent molecular clouds in Sgr C and the newly identified conical profile of Sgr B2 in {l ,b, v_LSR} space. These features have previously invoked supernova-driven shells and cloud-cloud collisions as explanations. We instead caution against structure identification in velocity-integrated emission maps. Three geometries describing the 3-D structure of the CMZ are investigated: i) two spiral arms; ii) a closed elliptical orbit; iii) an open stream. While two spiral arms and an open stream qualitatively reproduce the gas distribution, the most recent parameterisation of the closed elliptical orbit does not. Finally, we discuss how proper motion measurements of masers can distinguish between these geometries, and suggest that this effort should be focused on the 20 km/s and 50 km/s clouds and Sgr C.

preprint2016arXiv

OH+ and H2O+ absorption toward PKS1830-211

We report the detection of OH+ and H2O+ in the z=0.89 absorber toward the lensed quasar PKS1830-211. The abundance ratio of OH+ and H2O+ is used to quantify the molecular hydrogen fraction (fH2) and the cosmic-ray ionization rate of atomic hydrogen (zH) along two lines of sight, located at ~2 kpc and ~4 kpc to either side of the absorber's center. The molecular fraction decreases outwards, from ~0.04 to ~0.02, comparable to values measured in the Milky Way at similar galactocentric radii. For zH, we find values of ~2x10^-14 s^-1 and ~3x10^-15 s^-1, respectively, which are slightly higher than in the Milky Way at comparable galactocentric radii, possibly due to a higher average star formation activity in the z=0.89 absorber. The ALMA observations of OH+, H2O+, and other hydrides toward PKS1830-211 reveal the multi-phase composition of the absorbing gas. Taking the column density ratios along the southwest and northeast lines of sight as a proxy of molecular fraction, we classify the species ArH+, OH+, H2Cl+, H2O+, CH, and HF as tracing gases increasingly more molecular. Incidentally, our data allow us to improve the accuracy of H2O+ rest frequencies and thus refine the spectroscopic parameters.

preprint2016arXiv

The CO-H2 van der Waals complex and complex organic molecules in cold molecular clouds: a TMC-1C survey

Almost 200 different species have been detected in the interstellar medium (ISM) during the last decades, revealing not only simple species but complex molecules with more than 6 atoms. Other exotic compounds, like the weakly-bound dimer (H2)2, have also been detected in astronomical sources like Jupiter. We aim at detecting for the first time the CO-H2 van der Waals complex in the ISM, which if detected can be a sensitive indicator for low temperatures. We use the IRAM30m telescope, located in Pico Veleta (Spain), to search for the CO-H2 complex in a cold, dense core in TMC-1C (with a temperature of 10 K). All the brightest CO-H2 transitions in the 3 mm (80-110 GHz) band have been observed with a spectral resolution of 0.5-0.7 km/s, reaching a rms noise level of 2 mK. The simultaneous observation of a broad frequency band, 16 GHz, has allowed us to conduct a serendipitous spectral line survey. No lines belonging to the CO-H2 complex have been detected. We have set up a new, more stringent upper limit for its abundance to be [CO-H2]/[CO] = 5x10^{-6}, while we expect the abundance of the complex to be in the range 10^{-8}-10^{-3}. The spectral line survey has allowed us to detect 75 lines associated with 41 different species (including isotopologues). We detect a number of complex organic species, e.g. methyl cyanide (CH3CN), methanol (CH3OH), propyne (CH3CCH) and ketene (CH2CO), associated with cold gas (excitation temperatures about 7 K), confirming the presence of these complex species not only in warm objects but also in cold regimes.

preprint2016arXiv

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

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

preprint2016arXiv

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.

preprint2015arXiv

A study of the C$_3$H$_2$ isomers and isotopologues: first interstellar detection of HDCCC

The partially deuterated linear isomer HDCCC of the ubiquitous cyclic carbene ($c$-C$_3$H$_2$) was observed in the starless cores TMC-1C and L1544 at 96.9 GHz, and a confirming line was observed in TMC-1 at 19.38 GHz. To aid the identification in these narrow line sources, four centimetre-wave rotational transitions (two in the previously reported $K_a =0$ ladder, and two new ones in the $K_a =1$ ladder), and 23 transitions in the millimetre band between 96 and 272 GHz were measured in high-resolution laboratory spectra. Nine spectroscopic constants in a standard asymmetric top Hamiltonian allow the principal transitions of astronomical interest in the $K_a \le 3$ rotational ladders to be calculated to within 0.1 km s$^{-1}$ in radial velocity up to 400 GHz. Conclusive evidence for the identification of the two astronomical lines of HDCCC was provided by the $V_{\rm{LSR}}$ which is the same as that of the normal isotopic species (H$_2$CCC) in the three narrow line sources. In these sources, deuterium fractionation in singly substituted H$_2$CCC (HDCCC/H$_2$CCC $\sim4\%\text{-}19\%$) is comparable to that in $c$-C$_3$H$_2$ ($c$-C$_3$H$_2$/$c$-C$_3$HD $\sim5\%\text{-}17\%$), and similarly in doubly deuterated $c$-C$_3$H$_2$ ($c$-C$_3$D$_2$/$c$-C$_3$HD $\sim3\%\text{-}17\%$), implying that the efficiency of the deuteration processes in the H$_2$CCC and $c$-C$_3$H$_2$ isomers are comparable in dark clouds.

preprint2015arXiv

An Overview of the 2014 ALMA Long Baseline Campaign

A major goal of the Atacama Large Millimeter/submillimeter Array (ALMA) is to make accurate images with resolutions of tens of milliarcseconds, which at submillimeter (submm) wavelengths requires baselines up to ~15 km. To develop and test this capability, a Long Baseline Campaign (LBC) was carried out from September to late November 2014, culminating in end-to-end observations, calibrations, and imaging of selected Science Verification (SV) targets. This paper presents an overview of the campaign and its main results, including an investigation of the short-term coherence properties and systematic phase errors over the long baselines at the ALMA site, a summary of the SV targets and observations, and recommendations for science observing strategies at long baselines. Deep ALMA images of the quasar 3C138 at 97 and 241 GHz are also compared to VLA 43 GHz results, demonstrating an agreement at a level of a few percent. As a result of the extensive program of LBC testing, the highly successful SV imaging at long baselines achieved angular resolutions as fine as 19 mas at ~350 GHz. Observing with ALMA on baselines of up to 15 km is now possible, and opens up new parameter space for submm astronomy.

preprint2015arXiv

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

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

preprint2015arXiv

Sulphur-bearing molecules in diffuse molecular clouds: new results from SOFIA/GREAT and the IRAM 30 m telescope

We have observed five sulphur-bearing molecules in foreground diffuse molecular clouds lying along the sight-lines to five bright continuum sources. We have used the GREAT instrument on SOFIA to observe the 1383 GHz $^2Π_{3/2} J=5/2-3/2$ transitions of SH towards the star-forming regions W31C, G29.96-0.02, G34.3+0.1, W49N and W51, detecting foreground absorption towards all five sources; and the EMIR receivers on the IRAM 30m telescope at Pico Veleta to detect the H$_2$S 1(10)-1(01), CS J=2-1 and SO 3(2)-2(1) transitions. In nine foreground absorption components detected towards these sources, the inferred column densities of the four detected molecules showed relatively constant ratios, with N(SH)/N(H$_2$S) in the range 1.1 - 3.0, N(CS)/N(H$_2$S) in the range 0.32 - 0.61, and N(SO)/N(H$_2$S) in the range 0.08 - 0.30. The observed SH/H$_2$ ratios - in the range (0.5-2.6) $\times 10^{-8}$ - indicate that SH (and other sulphur-bearing molecules) account for << 1% of the gas-phase sulphur nuclei. The observed abundances of sulphur-bearing molecules, however, greatly exceed those predicted by standard models of cold diffuse molecular clouds, providing further evidence for the enhancement of endothermic reaction rates by elevated temperatures or ion-neutral drift. We have considered the observed abundance ratios in the context of shock and turbulent dissipation region (TDR) models. Using the TDR model, we find that the turbulent energy available at large scale in the diffuse ISM is sufficient to explain the observed column densities of SH and CS. Standard shock and TDR models, however, fail to reproduce the column densities of H$_2$S and SO by a factor of about 10; more elaborate shock models - in which account is taken of the velocity drift, relative to H$_2$, of SH molecules produced by the dissociative recombination of H$_3$S$^+$ - reduce this discrepancy to a factor ~ 3.

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

Velocity-resolved [CII] emission and [CII]/FIR Mapping along Orion with Herschel

We present the first 7.5'x11.5' velocity-resolved map of the [CII]158um line toward the Orion molecular cloud-1 (OMC-1) taken with the Herschel/HIFI instrument. In combination with far-infrared (FIR) photometric images and velocity-resolved maps of the H41alpha hydrogen recombination and CO J=2-1 lines, this data set provides an unprecedented view of the intricate small-scale kinematics of the ionized/PDR/molecular gas interfaces and of the radiative feedback from massive stars. The main contribution to the [CII] luminosity (~85%) is from the extended, FUV-illuminated face of the cloud G_0>500, n_H>5x10^3 cm^-3) and from dense PDRs (G_0~10^4, n_H~10^5 cm^-3) at the interface between OMC-1 and the HII region surrounding the Trapezium cluster. Around 15% of the [CII] emission arises from a different gas component without CO counterpart. The [CII] excitation, PDR gas turbulence, line opacity (from [13CII]) and role of the geometry of the illuminating stars with respect to the cloud are investigated. We construct maps of the [CII]/FIR and FIR/M_Gas ratios and show that [CII]/FIR decreases from the extended cloud component (10^-2-10^-3) to the more opaque star-forming cores (10^-3-10^-4). The lowest values are reminiscent of the "[CII] deficit" seen in local ultra-luminous IR galaxies hosting vigorous star formation. Spatial correlation analysis shows that the decreasing [CII]/FIR ratio correlates better with the column density of dust through the molecular cloud than with FIR/M_Gas. We conclude that the [CII] emitting column relative to the total dust column along each line of sight is responsible for the observed [CII]/FIR variations through the cloud.

preprint2014arXiv

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

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

preprint2014arXiv

Herschel Survey of Galactic OH+, H2O+, and H3O+: Probing the Molecular Hydrogen Fraction and Cosmic-Ray Ionization Rate

In diffuse interstellar clouds the chemistry that leads to the formation of the oxygen bearing ions OH+, H2O+, and H3O+ begins with the ionization of atomic hydrogen by cosmic rays, and continues through subsequent hydrogen abstraction reactions involving H2. Given these reaction pathways, the observed abundances of these molecules are useful in constraining both the total cosmic-ray ionization rate of atomic hydrogen (zeta_H) and molecular hydrogen fraction, f(H2). We present observations targeting transitions of OH+, H2O+, and H3O+ made with the Herschel Space Observatory along 20 Galactic sight lines toward bright submillimeter continuum sources. Both OH+ and H2O+ are detected in absorption in multiple velocity components along every sight line, but H3O+ is only detected along 7 sight lines. From the molecular abundances we compute f(H2) in multiple distinct components along each line of sight, and find a Gaussian distribution with mean and standard deviation 0.042+-0.018. This confirms previous findings that OH+ and H2O+ primarily reside in gas with low H2 fractions. We also infer zeta_H throughout our sample, and find a log-normal distribution with mean log(zeta_H)=-15.75, (zeta_H=1.78x10^-16 s^-1), and standard deviation 0.29 for gas within the Galactic disk, but outside of the Galactic center. This is in good agreement with the mean and distribution of cosmic-ray ionization rates previously inferred from H3+ observations. Ionization rates in the Galactic center tend to be 10--100 times larger than found in the Galactic disk, also in accord with prior studies.

preprint2014arXiv

Ionization toward the high-mass star-forming region NGC 6334 I

Context. Ionization plays a central role in the gas-phase chemistry of molecular clouds. Since ions are coupled with magnetic fields, which can in turn counteract gravitational collapse, it is of paramount importance to measure their abundance in star-forming regions. Aims. We use spectral line observations of the high-mass star-forming region NGC 6334 I to derive the abundance of two of the most abundant molecular ions, HCO+ and N2H+, and consequently, the cosmic ray ionization rate. In addition, the line profiles provide information about the kinematics of this region. Methods. We present high-resolution spectral line observations conducted with the HIFI instrument on board the Herschel Space Observatory of the rotational transitions with Jup > 5 of the molecular species C17O, C18O, HCO+, H13CO+, and N2H+. Results. The HCO+ and N2H+ line profiles display a redshifted asymmetry consistent with a region of expanding gas. We identify two emission components in the spectra, each with a different excitation, associated with the envelope of NGC 6334 I. The physical parameters obtained for the envelope are in agreement with previous models of the radial structure of NGC 6334 I based on submillimeter continuum observations. Based on our new Herschel/HIFI observations, combined with the predictions from a chemical model, we derive a cosmic ray ionization rate that is an order of magnitude higher than the canonical value of 10^(-17) s-1. Conclusions. We find evidence of an expansion of the envelope surrounding the hot core of NGC 6334 I, which is mainly driven by thermal pressure from the hot ionized gas in the region. The ionization rate seems to be dominated by cosmic rays originating from outside the source, although X-ray emission from the NGC 6334 I core could contribute to the ionization in the inner part of the envelope.

preprint2014arXiv

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

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

preprint2014arXiv

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

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

preprint2014arXiv

Widespread Rotationally-Hot Hydronium Ion in the Galactic Interstellar Medium

We present new observations of the (6,6) and (9,9) inversion transitions of the hydronium ion toward Sagittarius B2(N) and W31C. Sensitive observations toward Sagittarius B2(N) show that the high, ~ 500 K, rotational temperatures characterizing the population of the highly-excited metastable H3O+ rotational levels are present over a wide range of velocities corresponding to the Sagittarius B2 envelope, as well as the foreground gas clouds between the Sun and the source. Observations of the same lines toward W31C, a line of sight that does not intersect the Central Molecular Zone, but instead traces quiescent gas in the Galactic disk, also imply a high rotational temperature of ~ 380 K, well in excess of the kinetic temperature of the diffuse Galactic interstellar medium. While it is plausible that some fraction of the molecular gas may be heated to such high temperatures in the active environment of the Galactic center, characterized by high X-ray and cosmic ray fluxes, shocks and high degree of turbulence, this is unlikely in the largely quiescent environment of the Galactic disk clouds. We suggest instead that the highly-excited states of the hydronium ion are populated mainly by exoergic chemical formation processes and temperature describing the rotational level population does not represent the physical temperature of the medium. The same arguments may be applicable to other symmetric top rotors, such as ammonia. This offers a simple explanation to the long-standing puzzle of the presence of a pervasive, hot molecular gas component in the central region of the Milky Way. Moreover, our observations suggest that this is a universal process, not limited to the active environments associated with galactic nuclei.

preprint2013arXiv

Interstellar detection of c-C3D2

We report the first interstellar detection of c-C3D2. The doubly deuterated cyclopropenylidene, a carbene, has been detected toward the starless cores TMC- 1C and L1544 using the IRAM 30m telescope. The J(Ka,Kc) = 3(0,3)-2(1,2), 3(1,3)-2(0,2), and 2(2,1)-1(1,0) transitions of this species have been observed at 3 mm in both sources. The expected 1:2 intensity ratio has been found in the 3(0,3)-2(1,2) and 3(1,3)-2(0,2) lines, belonging to the para and ortho species respectively. We also observed lines of the main species, c-C3H2, the singly deuterated c-C3HD, and the species with one 13C off of the principal axis of the molecule, c-H13CC2H. The lines of c-C3D2 have been observed with high signal to noise ratio, better than 7.5 sigma in TMC-1C and 9 sigma in L1544. The abundance of doubly deuterated cyclopropenylidene with respect to the normal species is found to be (0.4 - 0.8)% in TMC-1C and (1.2 - 2.1)% in L1544. The deuteration of this small hydrocarbon ring is analysed with a comprehensive gas-grain model, the first including doubly deuterated species. The observed abundances of c-C3D2 can be explained solely by gas-phase processes, supporting the idea that c-C3H2 is a good indicator of gas-phase deuteration.

preprint2013arXiv

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

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

preprint2013arXiv

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

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

preprint2013arXiv

MALT90: The Millimetre Astronomy Legacy Team 90 GHz Survey

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

preprint2012arXiv

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

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

preprint2012arXiv

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

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

preprint2012arXiv

Modeling and Analysis Generic Interface for eXternal numerical codes (MAGIX)

The modeling and analysis generic interface for external numerical codes (MAGIX) is a model optimizer developed under the framework of the coherent set of astrophysical tools for spectroscopy (CATS) project. The MAGIX package provides a framework of an easy interface between existing codes and an iterating engine that attempts to minimize deviations of the model results from available observational data, constraining the values of the model parameters and providing corresponding error estimates. Many models (and, in principle, not only astrophysical models) can be plugged into MAGIX to explore their parameter space and find the set of parameter values that best fits observational/experimental data. MAGIX complies with the data structures and reduction tools of ALMA (Atacama Large Millimeter Array), but can be used with other astronomical and with non-astronomical data.

preprint2012arXiv

Molecular line survey of the high-mass star-forming region NGC 6334I with Herschel/HIFI and the SMA

We aim at deriving the molecular abundances and temperatures of the hot molecular cores in the high-mass star-forming region NGC 6334I and consequently deriving their physical and astrochemical conditions. In the framework of the Herschel guaranteed time key program CHESS, NGC 6334I is investigated by using HIFI aboard the Herschel Space Observatory. A spectral line survey is carried out in the frequency range 480-1907 GHz, and auxiliary interferometric data from the SMA in the 230 GHz band provide spatial information for disentangling the different physical components contributing to the HIFI spectrum. The spectral lines are identified with the aid of former surveys and spectral line catalogs. The observed spectrum is then compared to a simulated synthetic spectrum with XCLASS, assuming local thermal equilibrium, and best fit parameters are derived using the model optimization package MAGIX. A total of 46 molecules are identified, with 31 isotopologues, resulting in about 4300 emission and absorption lines. High- energy levels of the dominant emitter methanol and vibrationally excited HCN are detected. The number of unidentified lines remains low with 75, or less than 2 percent of the lines detected. The modeling suggests that several spectral features need two or more components to be fitted properly. Other components could be assigned to cold foreground clouds or to outflows, most visible in the SiO emission. A chemical variation between the two embedded hot cores is found, with more N-bearing molecules identified in SMA1 and O-bearing molecules in SMA2. Spectral line surveys give powerful insights into the study of the interstellar medium. Different molecules trace different physical conditions like the inner hot core, the envelope, the outflows or the cold foreground clouds. The derived molecular abundances provide further constraints for astrochemical models.

preprint2012arXiv

The Abundance, Ortho/Para Ratio, and Deuteration of Water in the High-Mass Star Forming Region NGC 6334 I

We present Herschel/HIFI observations of 30 transitions of water isotopologues toward the high-mass star forming region NGC 6334 I. The line profiles of H_2^{16}O, H_2^{17}O, H_2^{18}O, and HDO show a complex pattern of emission and absorption components associated with the embedded hot cores, a lower-density envelope, two outflow components, and several foreground clouds, some associated with the NGC 6334 complex, others seen in projection against the strong continuum background of the source. Our analysis reveals an H2O ortho/para ratio of 3 +/- 0.5 in the foreground clouds, as well as the outflow. The water abundance varies from ~10^{-8} in the foreground clouds and the outer envelope to ~10^{-6} in the hot core. The hot core abundance is two orders of magnitude below the chemical model predictions for dense, warm gas, but within the range of values found in other Herschel/HIFI studies of hot cores and hot corinos. This may be related to the relatively low gas and dust temperature (~100 K), or time dependent effects, resulting in a significant fraction of water molecules still locked up in dust grain mantles. The HDO/H_2O ratio in NGC 6334 I, ~2 10^{-4}, is also relatively low, but within the range found in other high-mass star forming regions.

preprint2011arXiv

Galactic structure based on the ATLASGAL 870mum survey

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

preprint2011arXiv

Herschel observations of Extra-Ordinary Sources: Methanol as a probe of physical conditions in Orion KL

We have examined methanol emission from Orion KL with of the {\em Herschel}/HIFI instrument, and detected two methanol bands centered at 524 GHz and 1061 GHz. The 524 GHz methanol band (observed in HIFI band 1a) is dominated by the isolated $Δ$J$=$0, K$=-4\rightarrow$-3, v$_t$$=$0 Q branch, and includes 25 E-type and 2 A-type transitions. The 1061 GHz methanol band (observed in HIFI band 4b) is dominated by the $Δ$J$=$0, K$=7\rightarrow$6, v$_t$$=$0 Q branch transitions which are mostly blended. We have used the isolated E-type v$_t$$=$0 methanol transitions to explore the physical conditions in the molecular gas. With HIFI's high velocity resolution, the methanol emission contributed by different spatial components along the line of sight toward Orion KL (hot core, low velocity flow, and compact ridge) can be distinguished and studied separately. The isolated transitions detected in these bands cover a broad energy range (upper state energy ranging from 80 K to 900 K), which provides a unique probe of the thermal structure in each spatial component. The observations further show that the compact ridge is externally heated. These observations demonstrate the power of methanol lines as probes of the physical conditions in warm regions in close proximity to young stars.

preprint2011arXiv

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

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

preprint2011arXiv

On the kinematics of massive star forming regions: the case of IRAS 17233-3606

Direct observations of accretion disks around high-mass young stellar objects would help to discriminate between different models of formation of massive stars. However, given the complexity of massive star forming regions, such studies are still limited in number. Additionally, there is still no general consensus on the molecular tracers to be used for such investigations. Because of its close distance and high luminosity, IRAS 17233-3606 is a potential good laboratory to search for traces of rotation in the inner gas around the protostar(s). Therefore, we selected the source for a detailed analysis of its molecular emission at 230 GHz with the SMA. We systematically investigated the velocity fields of transitions in the SMA spectra which are not affected by overlap with other transitions, and searched for coherent velocity gradients to compare them to the distribution of outflows in the region. Beside CO emission we also used high-angular H2 images to trace the outflow motions driven by the IRAS 17233-3606 cluster. We find linear velocity gradients in many transitions of the same molecular species and in several molecules. We report the first detection of HNCO in molecular outflows from massive YSOs. We discuss the CH3CN velocity gradient taking into account various scenarios: rotation, presence of multiple unresolved sources with different velocities, and outflow(s). Although other interpretations cannot be ruled out, we propose that the CH3CN emission might be affected by the outflows of the region. Higher angular observations are needed to discriminate between the different scenarios. The present observations, with the possible association of CH3CN with outflows in a few thousands AU around the YSOs' cluster, (i) question the choice of the tracer to probe rotating structures, and (ii) show the importance of the use of H2 images for detailed studies of kinematics.

preprint2011arXiv

Spectral imaging of the Central Molecular Zone in multiple 3-mm molecular lines

We have mapped 20 molecular lines in the Central Molecular Zone (CMZ) around the Galactic Centre, emitting from 85.3 to 93.3 GHz. This work used the 22-m Mopra radio telescope in Australia, equipped with the 8-GHz bandwidth UNSW-MOPS digital filter bank, obtaining \sim 2 km/s spectral and \sim 40 arcsec spatial resolution. The lines measured include emission from the c-C3H2, CH3CCH, HOCO+, SO, H13CN, H13CO+, SO, H13NC, C2H, HNCO, HCN, HCO+, HNC, HC3N, 13CS and N2H+ molecules. The area covered is Galactic longitude -0.7 to 1.8 deg. and latitude -0.3 to 0.2 deg., including the bright dust cores around Sgr A, Sgr B2, Sgr C and G1.6-0.025. We present images from this study and conduct a principal component analysis on the integrated emission from the brightest 8 lines. This is dominated by the first component, showing that the large-scale distribution of all molecules are very similar. We examine the line ratios and optical depths in selected apertures around the bright dust cores, as well as for the complete mapped region of the CMZ. We highlight the behaviour of the bright HCN, HNC and HCO+ line emission, together with that from the 13C isotopologues of these species, and compare the behaviour with that found in extra-galactic sources where the emission is unresolved spatially. We also find that the isotopologue line ratios (e.g. HCO+/H13CO+) rise significantly with increasing red-shifted velocity in some locations. Line luminosities are also calculated and compared to that of CO, as well as to line luminosities determined for external galaxies.

preprint2011arXiv

Submillimeter continuum observations of Sagittarius B2 at subarcsecond spatial resolution

We report the first high spatial resolution submillimeter continuum observations of the Sagittarius B2 cloud complex using the Submillimeter Array (SMA). With the subarcsecond resolution provided by the SMA, the two massive star-forming clumps Sgr B2(N) and Sgr B2(M) are resolved into multiple compact sources. In total, twelve submillimeter cores are identified in the Sgr B2(M) region, while only two components are observed in the Sgr B2(N) clump. The gas mass and column density are estimated from the dust continuum emission. We find that most of the cores have gas masses in excess of 100 M$_{\odot}$ and column densities above 10$^{25}$ cm$^{-2}$. The very fragmented appearance of Sgr B2(M), in contrast to the monolithic structure of Sgr B2 (N), suggests that the former is more evolved. The density profile of the Sgr B2(N)-SMA1 core is well fitted by a Plummer density distribution. This would lead one to believe that in the evolutionary sequence of the Sgr B2 cloud complex, a massive star forms first in an homogeneous core, and the rest of the cluster forms subsequently in the then fragmenting structure.

preprint2011arXiv

TIMASSS: The IRAS16293-2422 Millimeter And Submillimeter Spectral Survey. I. Observations, calibration and analysis of the line kinematics

While unbiased surveys observable from ground-based telescopes have previously been obtained towards several high mass protostars, very little exists on low mass protostars. To fill up this gap, we carried out a complete spectral survey of the bands at 3, 2, 1 and 0.8 mm towards the solar type protostar IRAS16293-2422. The observations covered about 200\,GHz and were obtained with the IRAM-30m and JCMT-15m telescopes. Particular attention was devoted to the inter-calibration of the obtained spectra with previous observations. All the lines detected with more than 3 sigma and free from obvious blending effects were fitted with Gaussians to estimate their basic kinematic properties. More than 4000 lines were detected (with sigma \geq 3) and identified, yielding a line density of approximatively 20 lines per GHz, comparable to previous surveys in massive hot cores. The vast majority (~2/3) of the lines are weak and due to complex organic molecules. The analysis of the profiles of more than 1000 lines belonging 70 species firmly establishes the presence of two distinct velocity components, associated with the two objects, A and B, forming the IRAS16293-2422 binary system. In the source A, the line widths of several species increase with the upper level energy of the transition, a behavior compatible with gas infalling towards a ~1 Mo object. The source B, which does not show this effect, might have a much lower central mass of ~0.1 Mo. The difference in the rest velocities of both objects is consistent with the hypothesis that the source B rotates around the source A. This spectral survey, although obtained with single-dish telescope with a low spatial resolution, allows to separate the emission from 2 different components, thanks to the large number of lines detected. The data of the survey are public and can be retrieved on the web site http://www-laog.obs.ujf-grenoble.fr/heberges/timasss.

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 study of three southern high-mass star-forming regions

Based on color-selected IRAS point sources, we have started to conduct a survey of 47 high-mass star-forming regions in the southern hemisphere in 870um dust continuum and molecular line emission in several frequency ranges between 290 GHz and 806 GHz. This paper describes the pilot study of the three sources IRAS12326-6245, IRAS16060-5146, and IRAS16065-5158. To characterize the physical and chemical properties of southern massive star-forming regions, the three high-luminosity southern hemisphere hot cores were observed with APEX in five frequency setups aimed at groups of lines from the following molecules: CH3OH, H2CO, and CH3CN. Using the LTE approximation, temperatures, source sizes, and column densities were determined through modeling of synthetic spectra with the XCLASS program. Dust continuum observations were done with the Large APEX BOlometer CAmera (LABOCA) at 870um and the 3mm continuum was imaged with the ATCA. Based on the detection of high-excitation CH3CN lines and lines from complex organic species, the three sources are classified as line rich, hot core type sources. For all three, the modeling indicates that the line emission emerges from a combination of an extended, cooler envelope, and a hot compact component. All three sources show an overabundance of oxygen-bearing species compared to nitrogen-bearing species. Based on the results obtained in the three sources, which served as templates for the survey, the most promising (and feasible) frequency setups for the remaining 44 sources were decided upon.

preprint2010arXiv

CHESS, Chemical Herschel surveys of star forming regions:Peering into the protostellar shock L1157-B1

The outflow driven by the low-mass class 0 protostar L1157 is the prototype of the so-called chemically active outflows. The bright bowshock B1 in the southern outflow lobe is a privileged testbed of magneto-hydrodynamical (MHD) shock models, for which dynamical and chemical processes are strongly interdependent. We present the first results of the unbiased spectral survey of the L1157-B1 bowshock, obtained in the framework of the key program "Chemical Herschel Surveys of Star Forming Regions" (CHESS). The main aim is to trace the warm and chemically enriched gas and to infer the excitation conditions in the shock region. The CO 5-4 and H2O lines have been detected at high-spectral resolution in the unbiased spectral survey of the HIFI-Band 1b spectral window (555-636 GHz), presented by Codella et al. in this volume. Complementary ground-based observations in the submm window help establish the origin of the emission detected in the main-beam of HIFI, and the physical conditions in the shock.}{Both lines exhibit broad wings, which extend to velocities much higher than reported up to now. We find that the molecular emission arises from two regions with distinct physical conditions: an extended, warm (100K), dense (3e5 cm-3) component at low-velocity, which dominates the water line flux in Band~1; a secondary component in a small region of B1 (a few arcsec) associated with high-velocity, hot (> 400 K) gas of moderate density ((1.0-3.0)e4 cm-3), which appears to dominate the flux of the water line at 179mu observed with PACS. The water abundance is enhanced by two orders of magnitude between the low- and the high-velocity component, from 8e-7 up to 8e-5. The properties of the high-velocity component agree well with the predictions of steady-state C-shock models.

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

Detection of interstellar oxidaniumyl: abundant H2O+ towards the star-forming regions DR21, Sgr B2, and NGC6334

We identify a prominent absorption feature at 1115 GHz, detected in first HIFI spectra towards high-mass star-forming regions, and interpret its astrophysical origin. The characteristic hyperfine pattern of the H2O+ ground-state rotational transition, and the lack of other known low-energy transitions in this frequency range, identifies the feature as H2O+ absorption against the dust continuum background and allows us to derive the velocity profile of the absorbing gas. By comparing this velocity profile with velocity profiles of other tracers in the DR21 star-forming region, we constrain the frequency of the transition and the conditions for its formation. In DR21, the velocity distribution of H2O+ matches that of the [CII] line at 158μ\m and of OH cm-wave absorption, both stemming from the hot and dense clump surfaces facing the HII-region and dynamically affected by the blister outflow. Diffuse foreground gas dominates the absorption towards Sgr B2. The integrated intensity of the absorption line allows us to derive lower limits to the H2O+ column density of 7.2e12 cm^-2 in NGC 6334, 2.3e13 cm^-2 in DR21, and 1.1e15 cm^-2 in Sgr B2.

preprint2010arXiv

Detection of OH$^+$ and H$_2$O$^+$ towards Orion~KL

We report observations of the reactive molecular ions OH$^+$, H$_2$O$^+$, and H$_3$O$^+$ towards Orion KL with Herschel/HIFI. All three $N=1-0$ fine-structure transitions of OH$^+$ at 909, 971, and 1033GHz and both fine-structure components of the doublet {\it ortho}-H$_2$O$^+$ $1_{11}-0_{00}$ transition at 1115 and 1139GHz were detected; an upper limit was obtained for H$_3$O$^+$. OH$^+$ and H$_2$O$^+$ are observed purely in absorption, showing a narrow component at the source velocity of 9 kms$^{-1}$, and a broad blueshifted absorption similar to that reported recently for HF and {\it para}-H$_{2}^{18}$O, and attributed to the low velocity outflow of Orion KL. We estimate column densities of OH$^+$ and H$_2$O$^+$ for the 9 km s$^{-1}$ component of $9 \pm 3 \times 10^{12}$cm$^{-2}$ and $7 \pm 2 \times 10^{12}$cm$^{-2}$, and those in the outflow of $1.9 \pm 0.7 \times 10^{13}$cm$^{-2}$ and $1.0 \pm 0.3 \times 10^{13}$cm$^{-2}$. Upper limits of $2.4\times 10^{12}$cm$^{-2}$ and $8.7\times 10^{12}$cm$^{-2}$ were derived for the column densities of {\it ortho} and {\it para}-H$_3$O$^+$ from transitions near 985 and 1657GHz. The column densities of the three ions are up to an order of magnitude lower than those obtained from recent observations of W31C and W49N. The comparatively low column densities may be explained by a higher gas density despite the assumption of a very high ionization rate.

preprint2010arXiv

First detection of ND in the solar-mass protostar IRAS16293-2422

In the past decade, much progress has been made in characterising the processes leading to the enhanced deuterium fractionation observed in the ISM and in particular in the cold, dense parts of star forming regions such as protostellar envelopes. Very high molecular D/H ratios have been found for saturated molecules and ions. However, little is known about the deuterium fractionation in radicals, even though simple radicals often represent an intermediate stage in the formation of more complex, saturated molecules. The imidogen radical NH is such an intermediate species for the ammonia synthesis in the gas phase. Herschel/HIFI represents a unique opportunity to study the deuteration and formation mechanisms of such species, which are not observable from the ground. We searched here for the deuterated radical ND in order to determine the deuterium fractionation of imidogen and constrain the deuteration mechanism of this species. We observed the solar-mass Class 0 protostar IRAS16293-2422 with the heterodyne instrument HIFI as part of the Herschel key programme CHESS (Chemical HErschel Surveys of Star forming regions). The deuterated form of the imidogen radical ND was detected and securely identified with 2 hyperfine component groups of its fundamental transition in absorption against the continuum background emitted from the nascent protostar. The 3 groups of hyperfine components of its hydrogenated counterpart NH were also detected in absorption. We derive a very high deuterium fractionation with an [ND]/[NH] ratio of between 30 and 100%. The deuterium fractionation of imidogen is of the same order of magnitude as that in other molecules, which suggests that an efficient deuterium fractionation mechanism is at play. We discuss two possible formation pathways for ND, by means of either the reaction of N+ with HD, or deuteron/proton exchange with NH.

preprint2010arXiv

H2CO and CH3OH maps of the Orion Bar photodissociation region

A previous analysis of methanol and formaldehyde towards the Orion Bar concluded that the two molecular species may trace different physical components, methanol the clumpy material, and formaldehyde the interclump medium. To verify this hypothesis, we performed multi-line mapping observations of the two molecules to study their spatial distributions. The observations were performed with the IRAM-30m telescope at 218 and 241 GHz, with an angular resolution of ~11''. Additional data for H2CO from the Plateau de Bure array are also discussed. The data were analysed using an LVG approach. Both molecules are detected in our single-dish data. Our data show that CH3OH peaks towards the clumps of the Bar, but its intensity decreases below the detection threshold in the interclump material. When averaging over a large region of the interclump medium, the strongest CH3OH line is detected with a peak intensity of ~0.06K. Formaldehyde also peaks on the clumps, but it is also detected in the interclump gas. We verified that the weak intensity of CH3OH in the interclump medium is not caused by the different excitation conditions of the interclump material, but reflects a decrease in the column density of methanol. The abundance of CH3OH relative to H2CO decreases by at least one order of magnitude from the dense clumps to the interclump medium.

preprint2010arXiv

Herschel observations of EXtra-Ordinary Sources (HEXOS): Observations of H2O and its isotopologues towards Orion KL

We report the detection of more than 48 velocity-resolved ground rotational state transitions of H2(16)O, H2(18)O, and H2(17)O - most for the first time - in both emission and absorption toward Orion KL using Herschel/HIFI. We show that a simple fit, constrained to match the known emission and absorption components along the line of sight, is in excellent agreement with the spectral profiles of all the water lines. Using the measured H2(18)O line fluxes, which are less affected by line opacity than their H2(16)O counterparts, and an escape probability method, the column densities of H2(18)O associated with each emission component are derived. We infer total water abundances of 7.4E-5, 1.0E-5, and 1.6E-5 for the plateau, hot core, and extended warm gas, respectively. In the case of the plateau, this value is consistent with previous measures of the Orion-KL water abundance as well as those of other molecular outflows. In the case of the hot core and extended warm gas, these values are somewhat higher than water abundances derived for other quiescent clouds, suggesting that these regions are likely experiencing enhanced water-ice sublimation from (and reduced freeze-out onto) grain surfaces due to the warmer dust in these sources.

preprint2010arXiv

Herschel observations of EXtra-Ordinary Sources: The Terahertz spectrum of Orion KL seen at high spectral resolution

We present the first high spectral resolution observations of Orion KL in the frequency ranges 1573.4 - 1702.8 GHz (band 6b) and 1788.4 - 1906.8 GHz (band 7b) obtained using the HIFI instrument on board the Herschel Space Observatory. We characterize the main emission lines found in the spectrum, which primarily arise from a range of components associated with Orion KL including the hot core, but also see widespread emission from components associated with molecular outflows traced by H2O, SO2, and OH. We find that the density of observed emission lines is significantly diminished in these bands compared to lower frequency Herschel/HIFI bands.

preprint2010arXiv

Herschel observations of ortho- and para-oxidaniumyl (H2O+) in spiral arm clouds toward Sgr B2(M)

H2O+ has been observed in its ortho- and para- states toward the massive star forming core Sgr B2(M), located close to the Galactic center. The observations show absorption in all spiral arm clouds between the Sun and Sgr B2. The average o/p ratio of H2O+ in most velocity intervals is 4.8, which corresponds to a nuclear spin temperature of 21 K. The relationship of this spin temperature to the formation temperature and current physical temperature of the gas hosting H2O+ is discussed, but no firm conclusion is reached. In the velocity interval 0 to 60 km/s, an ortho/para ratio of below unity is found, but if this is due to an artifact of contamination by other species or real is not clear.

preprint2010arXiv

Herschel Observations of the W43 "mini-starburst"

Aims: To explore the infrared and radio properties of one of the closest Galactic starburst regions. Methods: Images obtained with the Herschel Space Observatory at wavelengths of 70, 160, 250, 350, and 500 microns using the PACS and SPIRE arrays are analyzed and compared with radio continuum VLA data and 8 micron images from the Spitzer Space Telescope. The morphology of the far-infrared emission is combined with radial velocity measurements of millimeter and centimeter wavelength transitions to identify features likely to be associated with the W43 complex. Results: The W43 star-forming complex is resolved into a dense cluster of protostars, infrared dark clouds, and ridges of warm dust heated by massive stars. The 4 brightest compact sources with L > 1.5 x 10^4 Lsun embedded within the Z-shaped ridge of bright dust emission in W43 remain single at 4" (0.1 pc) resolution. These objects, likely to be massive protostars or compact clusters in early stages of evolution are embedded in clumps with masses of 10^3 to 10^4 Msun, but contribute only 2% to the 3.6 x 10^6 Lsun far-IR luminosity of W43 measured in a 16 by 16 pc box. The total mass of gas derived from the far-IR dust emission inside this region is ~10^6 Msun. Cometary dust clouds, compact 6 cm radio sources, and warm dust mark the locations of older populations of massive stars. Energy release has created a cavity blowing-out below the Galactic plane. Compression of molecular gas in the plane by the older HII region near G30.684-0.260 and the bipolar structure of the resulting younger W43 HII region may have triggered the current mini-star burst.

preprint2010arXiv

Herschel/HIFI discovery of interstellar chloronium (H$_2$Cl$^+$)

We report the first detection of chloronium, H$_2$Cl$^+$, in the interstellar medium, using the HIFI instrument aboard the \emph{Herschel} Space Observatory. The $2_{12}-1_{01}$ lines of ortho-H$_2^{35}$Cl$^+$ and ortho-H$_2^{37}$Cl$^+$ are detected in absorption towards NGC~6334I, and the $1_{11}-0_{00}$ transition of para-H$_2^{35}$Cl$^+$ is detected in absorption towards NGC~6334I and Sgr~B2(S). The H$_2$Cl$^+$ column densities are compared to those of the chemically-related species HCl. The derived HCl/H$_2$Cl$^+$ column density ratios, $\sim$1--10, are within the range predicted by models of diffuse and dense Photon Dominated Regions (PDRs). However, the observed H$_2$Cl$^+$ column densities, in excess of $10^{13}$~cm$^{-2}$, are significantly higher than the model predictions. Our observations demonstrate the outstanding spectroscopic capabilities of HIFI for detecting new interstellar molecules and providing key constraints for astrochemical models.

preprint2010arXiv

Herschel/HIFI observations of spectrally resolved methylidyne signatures toward the high-mass star-forming core NGC6334I

In contrast to extensively studied dense star-forming cores, little is known about diffuse gas surrounding star-forming regions. We study molecular gas in the high-mass star-forming region NGC6334I, which contains diffuse, quiescent components that are inconspicuous in widely used molecular tracers such as CO. We present Herschel/HIFI observations of CH toward NGC6334I observed as part of the CHESS key program. HIFI resolves the hyperfine components of its J=3/2-1/2 transition, observed in both emission and absorption. The CH emission appears close to the systemic velocity of NGC6334I, while its measured linewidth of 3 km/s is smaller than previously observed in dense gas tracers such as NH3 and SiO. The CH abundance in the hot core is 7 10^-11, two to three orders of magnitude lower than in diffuse clouds. While other studies find distinct outflows in, e.g., CO and H2O toward NGC6334I, we do not detect outflow signatures in CH. To explain the absorption signatures, at least two absorbing components are needed at -3.0 and +6.5 km/s with N(CH)=7 10^13 and 3 10^13 cm^-2. Two additional absorbing clouds are found at +8.0 and 0.0 km/s, both with N(CH)=2 10^13 cm^-2. Turbulent linewidths for the four absorption components vary between 1.5 and 5.0 km/s in FWHM. We constrain physical properties of our CH clouds by matching our CH absorbers with other absorption signatures. In the hot core, molecules such as H2O and CO trace gas that is heated and dynamically influenced by outflow activity, whereas CH traces more quiescent material. The four CH absorbers have column densities and turbulent properties consistent with diffuse clouds: two are located near NGC6334, and two are unrelated foreground clouds. Local density and dynamical effects influence the chemical composition of physical components of NGC6334, causing some components to be seen in CH but not in other tracers, and vice versa.

preprint2010arXiv

Ortho-to-para ratio of interstellar heavy water

Despite the low elemental deuterium abundance in the Galaxy, enhanced molecular D/H ratios have been found in the environments of low-mass star forming regions, and in particular the Class 0 protostar IRAS 16293-2422. The CHESS (Chemical HErschel Surveys of Star forming regions) Key Program aims at studying the molecular complexity of the interstellar medium. The high sensitivity and spectral resolution of the HIFI instrument provide a unique opportunity to observe the fundamental 1,1,1 - 0,0,0 transition of the ortho-D2O molecule, inaccessible from the ground, and to determine the ortho-to-para D2O ratio. We have detected the fundamental transition of the ortho-D2O molecule at 607.35 GHz towards IRAS 16293-2422. The line is seen in absorption with a line opacity of 0.62 +/- 0.11 (1 sigma). From the previous ground-based observations of the fundamental 1,1,0 - 1,0,1 transition of para-D2O seen in absorption at 316.80 GHz we estimate a line opacity of 0.26 +/- 0.05 (1 sigma). We show that the observed absorption is caused by the cold gas in the envelope of the protostar. Using these new observations, we estimate for the first time the ortho to para D2O ratio to be lower than 2.6 at a 3 sigma level of uncertainty, to be compared with the thermal equilibrium value of 2:1.

preprint2010arXiv

The distribution of water in the high-mass star-forming region NGC 6334I

We present observations of twelve rotational transitions of H2O-16, H2O-18, and H2O-17 toward the massive star-forming region NGC 6334 I, carried out with Herschel/HIFI as part of the guaranteed time key program Chemical HErschel Surveys of Star forming regions (CHESS). We analyze these observations to obtain insights into physical processes in this region. We identify three main gas components (hot core, cold foreground, and outflow) in NGC 6334 I and derive the physical conditions in these components. The hot core, identified by the emission in highly excited lines, shows a high excitation temperature of 200 K, whereas water in the foreground component is predominantly in the ortho- and para- ground states. The abundance of water varies between 4 10^-5 (outflow) and 10^-8 (cold foreground gas). This variation is most likely due to the freeze-out of water molecules onto dust grains. The H2O-18/H2O-17 abundance ratio is 3.2, which is consistent with the O-18/O-17 ratio determined from CO isotopologues. The ortho/para ratio in water appears to be relatively low 1.6(1) in the cold, quiescent gas, but close to the equilibrium value of three in the warmer outflow material (2.5(0.8)).

preprint2010arXiv

The methanol lines and hot core of OMC2-FIR4, an intermediate-mass protostar, with Herschel-HIFI

In contrast with numerous studies on the physical and chemical structure of low- and high-mass protostars, much less is known about their intermediate-mass counterparts, a class of objects that could help to elucidate the mechanisms of star formation on both ends of the mass range. We present the first results from a rich HIFI spectral dataset on an intermediate-mass protostar, OMC2-FIR4, obtained in the CHESS (Chemical HErschel SurveyS of star forming regions) key programme. The more than 100 methanol lines detected between 554 and 961 GHz cover a range in upper level energy of 40 to 540 K. Our physical interpretation focusses on the hot core, but likely the cold envelope and shocked regions also play a role in reality, because an analysis of the line profiles suggests the presence of multiple emission components. An upper limit of 10^-6 is placed on the methanol abundance in the hot core, using a population diagram, large-scale source model and other considerations. This value is consistent with abundances previously seen in low-mass hot cores. Furthermore, the highest energy lines at the highest frequencies display asymmetric profiles, which may arise from infall around the hot core.

preprint2009arXiv

APEX-CHAMP+ high-J CO observations of low-mass young stellar objects: I. The HH 46 envelope and outflow

AIMS Our aim is to characterize the size, mass, density and temperature profiles of the protostellar envelope of HH~46 IRS 1 and its surrounding cloud material as well as the effect the outflow has on its environment.METHODS The CHAMP+ and LABOCA arrays on the APEX telescope, combined with lower frequency line receivers, are used to obtain a large continuum map and smaller heterodyne maps in various isotopologues of CO and HCO+. The high-J lines of CO (6--5 and 7--6) and its isotopologues together with [C I] 2--1, observed with CHAMP+, are used to probe the warm molecular gas in the inner few hundred AU and in the outflowing gas. The data are interpreted with continuum and line radiative transfer models. RESULTS Broad outflow wings are seen in CO low- and high-J lines at several positions, constraining the gas temperatures to a constant value of ~100 K along the red outflow axis and to ~60 K for the blue outflow. The derived outflow mass is of order 0.4--0.8 M_sol, significantly higher than previously found. The bulk of the strong high-J CO line emission has a surprisingly narrow width, however, even at outflow positions. These lines cannot be fit by a passively heated model of the HH 46 IRS envelope. We propose that it originates from photon heating of the outflow cavity walls by ultraviolet photons originating in outflow shocks and the accretion disk boundary layers. At the position of the bow shock itself, the UV photons are energetic enough to dissociate CO. The envelope mass of ~5 M_sol is strongly concentrated towards HH 46 IRS with a density power law of -1.8.

preprint2009arXiv

APEX-CHAMP+ high-J CO observations of low-mass young stellar objects: II. Distribution and origin of warm molecular gas

The origin and heating mechanisms of warm (50<T<200 K) molecular gas in low-mass young stellar objects (YSOs) are strongly debated. Both passive heating of the inner collapsing envelope by the protostellar luminosity as well as active heating by shocks and by UV associated with the outflows or accretion have been proposed. We aim to characterize the warm gas within protosteller objects, and disentangle contributions from the (inner) envelope, bipolar outflows and the quiescent cloud. High-J CO maps (12CO J=6--5 and 7--6) of the immediate surroundings (up to 10,000 AU) of eight low-mass YSOs are obtained with the CHAMP+ 650/850 GHz array receiver mounted on the APEX telescope. In addition, isotopologue observations of the 13CO J=6--5 transition and [C I] 3P_2-3P_1 line were taken. Strong quiescent narrow-line 12CO 6--5 and 7--6 emission is seen toward all protostars. In the case of HH~46 and Ced 110 IRS 4, the on-source emission originates in material heated by UV photons scattered in the outflow cavity and not just by passive heating in the inner envelope. Warm quiescent gas is also present along the outflows, heated by UV photons from shocks. Shock-heated warm gas is only detected for Class 0 flows and the more massive Class I sources such as HH~46. Outflow temperatures, estimated from the CO 6--5 and 3--2 line wings, are ~100 K, close to model predictions, with the exception of the L~1551 IRS 5 and IRAS 12496-7650, for which temperatures <50 K are found. APEX-CHAMP+ is uniquely suited to directly probe a protostar's feedback on its accreting envelope gas in terms of heating, photodissociation, and outflow dispersal by mapping 1'x1' regions in high-J CO and [C I] lines.

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

Deuterium chemistry in the Orion Bar PDR - "warm" chemistry starring CH2D+

High levels of deuterium fractionation in gas-phase molecules are usually associated with cold regions, such as prestellar cores. Significant fractionation ratios are also observed in hot environments such as hot cores or hot corinos, where they are believed to be produced by the evaporation of the icy mantles surrounding dust grains, and thus are remnants of a previous cold (either gas-phase or grain surface) chemistry. The recent detection of DCN towards the Orion Bar, in a clump at a characteristic temperature of 70K, has shown that high deuterium fractionation can also be detected in PDRs. The Orion Bar clumps thus appear as a good environment for the observational study of deuterium fractionation in luke-warm gas, allowing to validate chemistry models in a different temperature range, where dominating fractionation processes are predicted to be different than in cold gas (< 20K). We aimed at studying observationally in detail the chemistry at work in the Orion Bar PDR, to understand if DCN is produced by ice mantle evaporation, or is the result of warm gas-phase chemistry, involving the CH2D+ precursor ion (which survives higher temperatures than the usual H2D+ precursor). Using the APEX and the IRAM 30m telescopes, we targetted selected deuterated species towards two clumps in the Orion Bar. We confirmed the detection of DCN and detected two new deuterated molecules (DCO+ and HDCO) towards one clump in the Orion Bar PDR. Significant deuterium fractionations are found for HCN and H2CO, but a low fractionation in HCO+. We also give upper limits for other molecules relevant for the deuterium chemistry. (...) We show evidence that warm deuterium chemistry driven by CH2D+ is at work in the clumps.

preprint2009arXiv

Extremely high velocity gas from the massive YSOs in IRAS 17233-3606

Molecular outflows from high-mass young stellar objects provide an excellent way to study the star formation process, and investigate if they are scaled-up versions of their low-mass counterparts. We selected the nearby massive star forming region IRAS 17233-3606 in order to study the kinematics and physics along the molecular outflow(s) originating from this source. We observed IRAS 17233-3606 in CO, a typical tracer of gas associated with molecular outflow, with the Submillimeter Array in the (2-1) transition, and with the APEX telescope in the higher excitation (6-5) line. Additional infrared H2 observations were performed with the UKIRT telescope. The CO data were analysed using a LVG approach. Our data resolve the previously detected molecular outflow in at least three different components, one of them with a high collimation factor ~4, and characterised by emission at extremely high velocities (|v-v_{LSR}|>120 km s^{-1}). The estimate of the kinematical outflow parameters are typical of massive YSOs, and in agreement with the measured bolometric luminosity of the source. The kinematic ages of the flows are in the range 10^2-10^3 yr, and therefore point to young objects that still did not reach the main sequence.

preprint2009arXiv

Increased complexity in interstellar chemistry: Detection and chemical modeling of ethyl formate and n-propyl cyanide in Sgr B2(N)

In recent years, organic molecules of increasing complexity have been found toward the prolific Galactic center source Sagittarius B2. We wish to explore the degree of complexity that the interstellar chemistry can reach in star-forming regions. We carried out a complete line survey of the hot cores Sgr B2(N) and (M) with the IRAM 30 m telescope in the 3 mm range. We analyzed this spectral survey in the LTE approximation. We modeled the emission of all known molecules simultaneously, which allows us to search for less abundant, more complex molecules. We compared the derived column densities with the predictions of a coupled gas-phase and grain-surface chemical code. We report the first detection in space of ethyl formate (C2H5OCHO) and n-propyl cyanide (C3H7CN) toward Sgr B2(N). The abundances of ethyl formate and n-propyl cyanide relative to H2 are estimated to be 3.6e-9 and 1.0e-9, respectively. Our chemical modeling suggests that the sequential, piecewise construction of ethyl and n-propyl cyanide from their constituent functional groups on the grain surfaces is their most likely formation route. Ethyl formate is primarily formed on the grains by adding CH3 to functional-group radicals derived from methyl formate, although ethanol may also be a precursor. The detection in Sgr B2(N) of the next stage of complexity in two classes of complex molecule, esters and alkyl cyanides, suggests that greater complexity in other classes of molecule may be present in the interstellar medium. {Abridged}

preprint2008arXiv

Detection of amino acetonitrile in Sgr B2(N)

Amino acids are building blocks of proteins and therefore key ingredients for the origin of life. The simplest amino acid, glycine, has long been searched for in the interstellar medium but has not been unambiguously detected so far. Since the search for glycine has turned out to be extremely difficult, we aimed at detecting a chemically related species (possibly a direct precursor), amino acetonitrile. With the IRAM 30m telescope we carried out a complete line survey of the hot core regions Sgr B2(N) and (M) in the 3 mm range, plus partial surveys at 2 and 1.3 mm. We analyzed our 30m line survey in the LTE approximation and modeled the emission of all known molecules simultaneously. We identified spectral features at the frequencies predicted for amino acetonitrile lines having intensities compatible with a unique rotation temperature. We also used the VLA to look for cold, extended emission from amino acetonitrile. We detected amino acetonitrile in Sgr B2(N) in our 30m telescope line survey and conducted confirmatory observations of selected lines with the IRAM PdB and the ATCA interferometers. The emission arises from a known hot core, the Large Molecule Heimat, and is compact with a source diameter of 2 arcsec (0.08 pc). We derived a column density of 2.8e16 cm-2, a temperature of 100 K, and a linewidth of 7 km s-1. Based on the simultaneously observed continuum emission, we calculated a density of 1.7e8 cm-3, a mass of 2340 Msun, and an amino acetonitrile fractional abundance of 2.2e-9. The high abundance and temperature may indicate that amino acetonitrile is formed by grain surface chemistry. We did not detect any hot, compact amino acetonitrile emission toward Sgr B2(M) or any cold, extended emission toward Sgr B2, with column-density upper limits of 6e15 and 3e12-14 cm-2, respectively. (abridged)

preprint2006arXiv

APEX 1 mm line survey of the Orion Bar

Unbiased molecular line surveys are a powerful tool for analyzing the physical and chemical parameters of astronomical objects and are the only means for obtaining a complete view of the molecular inventory for a given source. The present work stands for the first such investigation of a photon-dominated region. The first results of an ongoing millimeter-wave survey obtained towards the Orion Bar are reported. The APEX telescope in combination with the APEX-2A facility receiver was employed in this investigation. We derived the physical parameters of the gas through LVG analyses of the methanol and formaldehyde data. Information on the sulfur and deuterium chemistry of photon-dominated regions is obtained from detections of several sulfur-bearing molecules and DCN.

preprint2006arXiv

Dissipative structures of diffuse molecular gas: I - Broad HCO$^+$(1-0) emission

Results: We report the detection of broad HCO+(1-0) lines (10 mK < T < 0.5 K). The interpretation of 10 of the HCO+ velocity components is conducted in conjunction with that of the associated optically thin 13CO emission. The derived HCO+ column densities span a broad range, $10^{11}< N(HCO+)/Δv <4 \times 10^{12} \rm cm^2/(km/s^{-1}$, and the inferred HCO+ abundances, $2 \times 10^{-10}<X(HCO+) < 10^{-8}$, are more than one order of magnitude above those produced by steady-state chemistry in gas weakly shielded from UV photons, even at large densities. We compare our results with the predictions of non-equilibrium chemistry, swiftly triggered in bursts of turbulence dissipation and followed by a slow thermal and chemical relaxation phase, assumed isobaric. The set of values derived from the observations, i.e. large HCO+ abundances, temperatures in the range of 100--200 K and densities in the range 100--1000 cm3, unambiguously belongs to the relaxation phase. The kinematic properties of the gas suggest in turn that the observed HCO+ line emission results from a space-time average in the beam of the whole cycle followed by the gas and that the chemical enrichment is made at the expense of the non-thermal energy. Last, we show that the "warm chemistry" signature (i.e large abundances of HCO+, CH+, H20 and OH) acquired by the gas within a few hundred years, the duration of the impulsive chemical enrichment, is kept over more than thousand years. During the relaxation phase, the \wat/OH abundance ratio stays close to the value measured in diffuse gas by the SWAS satellite, while the OH/HCO+ ratio increases by more than one order of magnitude.

preprint2005arXiv

First astronomical detection of the CF+ ion

We report the first astronomical detection of the CF+ (fluoromethylidynium) ion obtained by recent observations of its J = 1 - 0 (102.6 GHz), J = 2 - 1 (205.2 GHz), and J = 3 - 2 (307.7 GHz) pure rotational emissions toward the Orion Bar. Our search for CF+, carried out using the IRAM 30m and APEX 12m telescopes, was motivated by recent theoretical models that predict CF+ abundances of a few x E-10 in UV-irradiated molecular regions where C+ is present. The measurements confirm the predictions. They provide support for our current theories of interstellar fluorine chemistry, which suggest that hydrogen fluoride should be ubiquitous in interstellar gas clouds.

preprint2003arXiv

The line-of-sight distribution of water in the SgrB2 complex

We report the detection, with the Caltech Submillimeter Observatory, of the 894-GHz HDO(111-000) transition, observed in absorption against the background continuum emission of the SgrB2 cores M and N. Radiative transfer modeling of this feature, together with the published data set of mm and submm HDO and H2O-18 transitions, suggests that ground-state absorption features from deuterated and non-deuterated water trace different gas components along the line of sight. In particular, while the HDO line seems to be produced by the large column densities of gas located in the SgrB2 warm envelope, the H2O-18 ground-state transition detected by SWAS and KAO at 548 GHz (Neufeld et al. 2000; Zmuidzinas et al. 1995a) is instead a product of the hot, diffuse, thin gas layer lying in the foreground of the SgrB2 complex.

preprint2002arXiv

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

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