Methanol masers : Reliable tracers of the early stages of high-mass star formation

The GLIMPSE and MSX surveys have been used to examine the mid-infrared properties of a statistically complete sample of 6.7 GHz methanol masers. The GLIMPSE point sources associated with methanol masers are clearly distinguished from the majority, typically having extremely red mid-infrared colors, similar to those expected of low-mass class 0 young stellar objects. The intensity of the GLIMPSE sources associated with methanol masers is typically 4 magnitudes brighter at 8.0 micron than at 3.6 micron. Targeted searches towards GLIMPSE point sources with [3.6]-[4.5] > 1.3 and an 8.0 micron magnitude less than 10 will detect more than 80% of class II methanol masers. Many of the methanol masers are associated with sources within infrared dark clouds (IRDC) which are believed to mark regions where high-mass star formation is in its very early stages. The presence of class II methanol masers in a significant fraction of IRDC suggests that high-mass star formation is common in these regions. Different maser species are thought to trace different evolutionary phases of the high-mass star formation process. Comparison of the properties of the GLIMPSE sources associated with class II methanol masers and other maser species shows interesting trends, consistent with class I methanol masers tracing a generally earlier evolutionary phase and OH masers tracing a later evolutionary phase.Comment: 45 pages, 19 figures, accepted for publication in Ap

ISOCAM observations of the L1551 star formation region

The results of a deep mid-IR ISOCAM survey of the L1551 dark molecular cloud are presented. The aim of this survey is a search for new YSO (Young Stellar Object) candidates, using two broad-band filters centred at 6.7 and 14.3 micron. Although two regions close to the centre of L1551 had to be avoided due to saturation problems, 96 sources were detected in total (76 sources at 6.7 micron and 44 sources at 14.3 micron). Using the 24 sources detected in both filters, 14 were found to have intrinsic mid-IR excess at 14.3 micron and were therefore classified as YSO candidates. Using additional observations in B, V, I, J, H and K obtained from the ground, most candidates detected at these wavelengths were confirmed to have mid-IR excess at 6.7 micron as well, and three additional YSO candidates were found. Prior to this survey only three YSOs were known in the observed region (avoiding L1551 IRS5/NE and HL/XZ Tau). This survey reveals 15 new YSO candidates, although several of these are uncertain due to their extended nature either in the mid-IR or in the optical/near-IR observations. Two of the sources with mid-IR excess are previously known YSOs, one is a brown dwarf MHO 5 and the other is the well known T Tauri star HH30, consisting of an outflow and an optically thick disk seen edge on.Comment: 14 Pages, 8 Figure

ISOCAM observations of the rho Ophiuchi cloud: Luminosity and mass functions of the pre-main sequence embedded cluster

We present the results of the first extensive mid-infrared (IR) imaging survey of the rho Ophiuchi embedded cluster, performed with the ISOCAM camera on board the ISO satellite. The main molecular cloud L1688, as well as L1689N and L1689S, have been completely surveyed for point sources at 6.7 and 14.3 micron. A total of 425 sources are detected including 16 Class I, 123 Class II, and 77 Class III young stellar objects (YSOs). Essentially all of the mid-IR sources coincide with near-IR sources, but a large proportion of them are recognized for the first time as YSOs. Our dual-wavelength survey allows us to identify essentially all the YSOs with IR excess in the embedded cluster down to Fnu ~ 10 - 15 mJy. It more than doubles the known population of Class II YSOs and represents the most complete census to date of newly formed stars in the rho Ophiuchi central region. The stellar luminosity function of the complete sample of Class II YSOs is derived with a good accuracy down to L= 0.03 Lsun. A modeling of this lumino- sity function, using available pre-main sequence tracks and plausible star for- mation histories, allows us to derive the mass distribution of the Class II YSOs which arguably reflects the IMF of the embedded cluster. We estimate that the IMF in rho Ophiuchi is well described by a two-component power law with a low- mass index of -0.35+/-0.25, a high-mass index of -1.7 (to be compared with the Salpeter value of -1.35), and a break occurring at M = 0.55+/-0.25 Msun. This IMF is flat with no evidence for a low-mass cutoff down to at least 0.06 Msun.Comment: A&A Document Class -- version 5.01, 27 pages, 10 figures v2: typos added including few changes in source numberin

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

Mapping the column density and dust temperature structure of IRDCs with Herschel

Infrared dark clouds (IRDCs) are cold and dense reservoirs of gas potentially available to form stars. Many of these clouds are likely to be pristine structures representing the initial conditions for star formation. The study presented here aims to construct and analyze accurate column density and dust temperature maps of IRDCs by using the first Herschel data from the Hi-GAL galactic plane survey. These fundamental quantities, are essential for understanding processes such as fragmentation in the early stages of the formation of stars in molecular clouds. We have developed a simple pixel-by-pixel SED fitting method, which accounts for the background emission. By fitting a grey-body function at each position, we recover the spatial variations in both the dust column density and temperature within the IRDCs. This method is applied to a sample of 22 IRDCs exhibiting a range of angular sizes and peak column densities. Our analysis shows that the dust temperature decreases significantly within IRDCs, from background temperatures of 20-30 K to minimum temperatures of 8-15 K within the clouds, showing that dense molecular clouds are not isothermal. Temperature gradients have most likely an important impact on the fragmentation of IRDCs. Local temperature minima are strongly correlated with column density peaks, which in a few cases reach NH2 = 1 x 10^{23} cm^{-2}, identifying these clouds as candidate massive prestellar cores. Applying this technique to the full Hi-GAL data set will provide important constraints on the fragmentation and thermal properties of IRDCs, and help identify hundreds of massive prestellar core candidates.Comment: Accepted for publication in A&A Herschel special issu

Star Formation in the Milky Way and Nearby Galaxies

We review progress over the past decade in observations of large-scale star formation, with a focus on the interface between extragalactic and Galactic studies. Methods of measuring gas contents and star formation rates are discussed, and updated prescriptions for calculating star formation rates are provided. We review relations between star formation and gas on scales ranging from entire galaxies to individual molecular clouds.Comment: 55 pages, 15 figures, in press for Annual Reviews of Astronomy and Astrophysics; Updated with corrected equation 5, improved references, and other minor change

Strong absorption by interstellar hydrogen fluoride: Herschel/HIFI observations of the sight-line to G10.6-0.4 (W31C)

We report the detection of strong absorption by interstellar hydrogen fluoride along the sight-line to the submillimeter continuum source G10.6-0.4 (W31C). We have used Herschel's HIFI instrument, in dual beam switch mode, to observe the 1232.4763 GHz J=1-0 HF transition in the upper sideband of the Band 5a receiver. The resultant spectrum shows weak HF emission from G10.6-0.4 at LSR velocities in the range -10 to -3 km/s, accompanied by strong absorption by foreground material at LSR velocities in the range 15 to 50 km/s. The spectrum is similar to that of the 1113.3430 GHz 1(11)-0(00) transition of para-water, although at some frequencies the HF (hydrogen fluoride) optical depth clearly exceeds that of para-H2O. The optically-thick HF absorption that we have observed places a conservative lower limit of 1.6E+14 cm-2 on the HF column density along the sight-line to G10.6-0.4. Our lower limit on the HF abundance, 6E-9 relative to hydrogen nuclei, implies that hydrogen fluoride accounts for between ~ 30 and 100% of the fluorine nuclei in the gas phase along this sight-line. This observation corroborates theoretical predictions that - because the unique thermochemistry of fluorine permits the exothermic reaction of F atoms with molecular hydrogen - HF will be the dominant reservoir of interstellar fluorine under a wide range of conditions.Comment: Accepted for publication in Astronomy and Astrophysics (Herschel special issue). This revised version corrects a typographic error in the HTML abstract, in which the lower limit on the HF abundance (should be 6E-9) was previously misstated. The abstract in the PDF version is correct and the latter has not been modifie

Chemistry in Infrared Dark Cloud Clumps: a Molecular Line Survey at 3 mm

We have observed 37 Infrared Dark Clouds (IRDCs), containing a total of 159 clumps, in high-density molecular tracers at 3 mm using the 22-meter ATNF Mopra Telescope located in Australia. After determining kinematic distances, we eliminated clumps that are not located in IRDCs and clumps with a separation between them of less than one Mopra beam. Our final sample consists of 92 IRDC clumps. The most commonly detected molecular lines are (detection rates higher than 8%): N2H+, HNC, HN13C, HCO+, H13CO+, HCN, C2H, HC3N, HNCO, and SiO. We investigate the behavior of the different molecular tracers and look for chemical variations as a function of an evolutionary sequence based on Spitzer IRAC and MIPS emission. We find that the molecular tracers behave differently through the evolutionary sequence and some of them can be used to yield useful relative age information. The presence of HNC and N2H+ lines do not depend on the star formation activity. On the other hand, HC3N, HNCO, and SiO are predominantly detected in later stages of evolution. Optical depth calculations show that in IRDC clumps the N2H+ line is optically thin, the C2H line is moderately optically thick, and HNC and HCO+ are optically thick. The HCN hyperfine transitions are blended, and, in addition, show self-absorbed line profiles and extended wing emission. These factors combined prevent the use of HCN hyperfine transitions for the calculation of physical parameters. Total column densities of the different molecules, except C2H, increase with the evolutionary stage of the clumps. Molecular abundances increase with the evolutionary stage for N2H+ and HCO+. The N2H+/HCO+ and N2H+/HNC abudance ratios act as chemical clocks, increasing with the evolution of the clumps.Comment: Accepted to ApJ. 29 page

Mass-Losing Semiregular Variable Stars in Baade's Windows

By cross-correlating the results of two recent large-scale surveys, the general properties of a well defined sample of semi-regular variable stars have been determined. ISOGAL mid-infrared photometry and MACHO lightcurves are assembled for approximately 300 stars in the Baade's Windows of low extinction towards the Galactic bulge. These stars are mainly giants of late M spectral type, evolving along the asymptotic giant branch (AGB). They are found to possess a wide and continuous distribution of pulsation periods and to obey an approximate log~period -- bolometric magnitude relation or set of such relations. Approximate mass-loss rates in the range of 1e-8 to 5e-7 M_sun per year are derived from ISOGAL mid-infrared photometry and models of stellar spectra adjusted for the presence of optically-thin circumstellar silicate dust. Mass-loss rates depend on luminosity and pulsation period. Some stars lose mass as rapidly as short-period Miras but do not show Mira-like amplitudes. A period of 70 days or longer is a necessary but not a sufficient condition for mass loss to occur. For AGB stars in the mass-loss ranges that we observe, the functional dependence of mass-loss rate on temperature and luminosity is found to be in agreement with recent theoretical predictions. If we include our mass-loss rates with a sample of extreme mass-losing AGB stars in the Large Magellanic Cloud, we get the general result for AGB stars that mass-loss rate is proportional to luminosity^{2.7}, valid for AGB stars with 10^{-8} to 10^{-4} M_sun per year (Abridged).Comment: to appear in The Astrophysical Journal, 51 pages, 9 figures, 3 tables; table 1 will be available in machine-readable format at the electronic Ap