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

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

Determining the Parameters of Massive Protostellar Clouds via Radiative Transfer Modeling

A one-dimensional method for reconstructing the structure of prestellar and protostellar clouds is presented. The method is based on radiative transfer computations and a comparison of theoretical and observed intensity distributions at both millimeter and infrared wavelengths. The radiative transfer of dust emission is modeled for specified parameters of the density distribution, central star, and external background, and the theoretical distribution of the dust temperature inside the cloud is determined. The intensity distributions at millimeter and IR wavelengths are computed and quantitatively compared with observational data. The best-fit model parameters are determined using a genetic minimization algorithm, which makes it possible to reveal the ranges of parameter degeneracy as well. The method is illustrated by modeling the structure of the two infrared dark clouds IRDC-320.27+029 (P2) and IRDC-321.73+005 (P2). The derived density and temperature distributions can be used to model the chemical structure and spectral maps in molecular lines.Comment: Accepted for publication in Astronomy Report

How do methanol masers manage to appear in the youngest star vicinities and isolated molecular clumps?

General characteristics of methanol (CH3OH) maser emission are summarized. It is shown that methanol maser sources are concentrated in the spiral arms. Most of the methanol maser sources from the Perseus arm are associated with embedded stellar clusters and a considerable portion is situated close to compact HII regions. Almost 1/3 of the Perseus Arm sources lie at the edges of optically identified HII regions which means that massive star formation in the Perseus Arm is to a great extent triggered by local phenomena. A multiline analysis of the methanol masers allows us to determine the physical parameters in the regions of maser formation. Maser modelling shows that class II methanol masers can be pumped by the radiation of the warm dust as well as by free-free emission of a hypercompact region hcHII with a turnover frequency exceeding 100 GHz. Methanol masers of both classes can reside in the vicinity of hcHIIs. Modelling shows that periodic changes of maser fluxes can be reproduced by variations of the dust temperature by a few percent which may be caused by variations in the brightness of the central young stellar object reflecting the character of the accretion process. Sensitive observations have shown that the masers with low flux densities can still have considerable amplification factors. The analysis of class I maser surveys allows us to identify four distinct regimes that differ by the series of their brightest lines.Comment: 8 pages, 4 figures, invited presentation at IAU242 "Astrophysical Masers and their environments

Physical properties of Southern infrared dark clouds

It is commonly assumed that cold and dense Infrared Dark Clouds (IRDCs) likely represent the birth sites massive stars. Therefore, this class of objects gets increasing attention. To enlarge the sample of well-characterised IRDCs in the southern hemisphere, we have set up a program to study the gas and dust of southern IRDCs. The present paper aims at characterizing the continuuum properties of this sample of objects. We cross-correlated 1.2 mm continuum data from SIMBA@SEST with Spitzer/GLIMPSE images to establish the connection between emission sources at millimeter wavelengths and the IRDCs we see at 8 $\mu$m in absorption against the bright PAH background. Analysing the dust emission and extinction leads to a determination of masses and column densities, which are important quantities in characterizing the initial conditions of massive star formation. The total masses of the IRDCs were found to range from 150 to 1150 $\rm M_\odot$ (emission data) and from 300 to 1750 $\rm M_\odot$ (extinction data). We derived peak column densities between 0.9 and 4.6 $\times 10^{22}$ cm$^{-2}$ (emission data) and 2.1 and 5.4 $\times 10^{22}$ cm$^{-2}$ (extinction data). We demonstrate that the extinction method fails for very high extinction values (and column densities) beyond A$_{\rm V}$ values of roughly 75 mag according to the Weingartner & Draine (2001) extinction relation $R_{\rm V} = 5.5$ model B. The derived column densities, taking into account the spatial resolution effects, are beyond the column density threshold of 3.0 $\times 10^{23}$ cm$^{-2}$ required by theoretical considerations for massive star formation. We conclude that the values for column densities derived for the selected IRDC sample make these objects excellent candidates for objects in the earliest stages of massive star formation.Comment: Accepted for publication in Astronomy & Astrophysic

Probing the Early Evolution of Young High-Mass Stars

Near-infrared imaging surveys of high-mass star-forming regions reveal an amazingly complex interplay between star formation and the environment (Churchwell et al. 2006; Alvarez et al. 2004). By means of near-IR spectroscopy the embedded massive young stars can be characterized and placed in the context of their birth site. However, so far spectroscopic surveys have been hopelessly incomplete, hampering any systematic study of these very young massive stars. New integral field instrumentation available at ESO has opened the possibility to take a huge step forward by obtaining a full spectral inventory of the youngest massive stellar populations in star-forming regions currently accessible. Simultaneously, the analysis of the extended emission allows the characterization of the environmental conditions. The Formation and Early Evolution of Massive Stars (FEMS) collaboration aims at setting up a large observing campaign to obtain a full census of the stellar content, ionized material, outflows and PDR's over a sample of regions that covers a large parameter space. Complementary radio, mm and infrared observations will be used for the characterization of the deeply embedded population. For the first eight regions we have obtained 40 hours of SINFONI observations. In this contribution, we present the first results on three regions that illustrate the potential of this strategy.Comment: To appear in ASP Conf. Proceedings of "Massive Star Formation: Observations confront Theory", H. Beuther et al. (eds.), held in Heidelberg, September 200

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\,$\mu$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.Comment: 13 pages, 11 figures, accepted for Astronomy and Astrophysics, high-resolution version of paper at http://www.mpia.de/homes/beuther/papers.htm

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

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

ATLASGAL - Ammonia observations towards the southern Galactic Plane

Context: The initial conditions of molecular clumps in which high-mass stars form are poorly understood. In particular, a more detailed study of the earliest evolutionary phases is needed. The APEX Telescope Large Area Survey of the whole inner Galactic disk at 870 μm, ATLASGAL, has therefore been conducted to discover high-mass star-forming regions at different evolutionary phases. Aims: We derive properties such as velocities, rotational temperatures, column densities, and abundances of a large sample of southern ATLASGAL clumps in the fourth quadrant. Methods: Using the Parkes telescope, we observed the NH3 (1, 1) to (3, 3) inversion transitions towards 354 dust clumps detected by ATLASGAL within a Galactic longitude range between 300° and 359° and a latitude within ± 1.5°. For a subsample of 289 sources, the N2H+ (1–0) line was measured with the Mopra telescope. Results: We measured a median NH3 (1, 1) line width of ~ 2 km s-1, rotational temperatures from 12 to 28 K with a mean of 18 K, and source-averaged NH3 abundances from 1.6 × 10-6 to 10-8. For a subsample with detected NH3 (2, 2) hyperfine components, we found that the commonly used method to compute the (2, 2) optical depth from the (1, 1) optical depth and the (2, 2) to (1, 1) main beam brightness temperature ratio leads to an underestimation of the rotational temperature and column density. A larger median virial parameter of ~ 1 is determined using the broader N2H+ line width than is estimated from the NH3 line width of ~ 0.5 with a general trend of a decreasing virial parameter with increasing gas mass. We obtain a rising NH3 (1, 1)/N2H+ line-width ratio with increasing rotational temperature. Conclusions: A comparison of NH3 line parameters of ATLASGAL clumps to cores in nearby molecular clouds reveals smaller velocity dispersions in low-mass than high-mass star-forming regions and a warmer surrounding of ATLASGAL clumps than the surrounding of low-mass cores. The NH3 (1, 1) inversion transition of 49% of the sources shows hyperfine structure anomalies. The intensity ratio of the outer hyperfine structure lines with a median of 1.27 ± 0.03 and a standard deviation of 0.45 is significantly higher than 1, while the intensity ratios of the inner satellites with a median of 0.9 ± 0.02 and standard deviation of 0.3 and the sum of the inner and outer hyperfine components with a median of 1.06 ± 0.02 and standard deviation of 0.37 are closer to 1