No Arabic abstract
Determining the efficiency with which gas is converted into stars in galaxies requires an accurate determination of the total reservoir of molecular gas mass. However, despite being the most abundant molecule in the Universe, H$_2$ is challenging to detect through direct observations and indirect methods have to be used to estimate the total molecular gas reservoir. These are often based on scaling relations from tracers such as CO or dust, and are generally calibrated in the Milky Way. Yet, evidence that these scaling relations are environmentally dependent is growing. In particular, the commonly used CO-to-H$_2$ conversion factor (X$_{rm CO}$) is expected to be higher in metal-poor and/or strongly UV-irradiated environments. We use new SOFIA/FIFI-LS observations of far-infrared fine structure lines from the ionised and neutral gas and the Meudon photodissociation region model to constrain the physical properties and the structure of the gas in the massive star-forming region of 30 Doradus in the Large Magellanic Cloud, and determine the spatially resolved distribution of the total reservoir of molecular gas in the proximity of the young massive cluster R136. We compare this value with the molecular gas mass inferred from ground-based CO observations and dust-based estimates to quantify the impact of this extreme environment on commonly used tracers of the molecular gas. We find that the strong radiation field combined with the half-solar metallicity of the surrounding gas are responsible for a large reservoir of CO-dark molecular gas, leaving a large fraction of the total H$_2$ gas (> 75%) undetected when adopting a standard X$_{rm CO}$ factor in this massive star-forming region.
With an aim of probing the physical conditions and excitation mechanisms of warm molecular gas in individual star-forming regions, we performed Herschel SPIRE FTS observations of 30 Doradus in the LMC. In our FTS observations, important FIR cooling lines in the ISM, including CO J=4-3 to 13-12, [CI] 370 micron, and [NII] 205 micron, were clearly detected. In combination with ground-based CO data, we then constructed CO spectral line energy distributions (SLEDs) on 10 pc scales over a 60 pc x 60 pc area and found that the shape of the observed CO SLEDs considerably changes across 30 Doradus, e.g., the peak transition varies from J=6-5 to 10-9, while the slope characterized by the high-to-intermediate J ratio ranges from 0.4 to 1.8. To examine the source(s) of these variations in CO transitions, we analyzed the CO observations, along with [CII] 158 micron, [CI] 370 micron, [OI] 145 micron, H2 0-0 S(3), and FIR luminosity data, using state-of-the-art models of PDRs and shocks. Our detailed modeling showed that the observed CO emission likely originates from highly-compressed (thermal pressure ~ 1e7-1e9 K cm-3) clumps on 0.7-2 pc scales, which could be produced by either UV photons (UV radiation field ~ 1e3-1e5 Mathis fields) or low-velocity C-type shocks (pre-shock medium density ~ 1e4-1e6 cm-3 and shock velocity ~ 5-10 km s-1). Considering the stellar content in 30 Doradus, however, we tentatively excluded the stellar origin of CO excitation and concluded that low-velocity shocks driven by kpc scale processes (e.g., interaction between the Milky Way and the Magellanic Clouds) are likely the dominant source of heating for CO. The shocked CO-bright medium was then found to be warm (temperature ~ 100-500 K) and surrounded by a UV-regulated low pressure component (a few (1e4-1e5) K cm-3) that is bright in [CII] 158 micron, [CI] 370 micron, [OI] 145 micron, and FIR dust continuum emission.
We present deep Hubble Space Telescope (HST) NICMOS 2 F160W band observations of the central 56*57 (14pc*14.25pc) region around R136 in the starburst cluster 30 Dor (NGC 2070) located in the Large Magellanic Cloud. Our aim is to derive the stellar Initial Mass Function (IMF) down to ~1 Msun in order to test whether the IMF in a massive metal-poor cluster is similar to that observed in nearby young clusters and the field in our Galaxy. We estimate the mean age of the cluster to be 3 Myr by combining our F160W photometry with previously obtained HST WFPC2 optical F555W and F814W band photometry and comparing the stellar locus in the color-magnitude diagram with main sequence and pre-main sequence isochrones. The color-magnitude diagrams show the presence of differential extinction and possibly an age spread of a few megayears. We convert the magnitudes into masses adopting both a single mean age of 3 Myr isochrone and a constant star formation history from 2 to 4 Myr. We derive the IMF after correcting for incompleteness due to crowding. The faintest stars detected have a mass of 0.5 Msun and the data are more than 50% complete outside a radius of 5 pc down to a mass limit of 1.1 Msun for 3 Myr old objects. We find an IMF of dN/dlog(M) M^(-1.20+-0.2) over the mass range 1.1--20 Msun only slightly shallower than a Salpeter IMF. In particular, we find no strong evidence for a flattening of the IMF down to 1.1 Msun at a distance of 5 pc from the center, in contrast to a flattening at 2 Msun at a radius of 2 pc, reported in a previous optical HST study. We examine several possible reasons for the different results. If the IMF determined here applies to the whole cluster, the cluster would be massive enough to remain bound and evolve into a relatively low-mass globular cluster.
We present the first fully calibrated H$_2$, 1-0 S(1) image of the entire 30 Doradus nebula. The observations were conducted using the NOAO Extremely Wide-Field Infrared Imager on the CTIO 4-meter Blanco Telescope. Together with a NEWFIRM Br$gamma$ image of 30 Doradus, our data reveal the morphologies of the warm molecular gas and ionized gas in 30 Doradus. The brightest H$_2$-emitting area, which extends from the northeast to the southwest of R136, is a photodissociation region viewed face-on, while many clumps and pillar features located at the outer shells of 30 Doradus are photodissociation regions viewed edge-on. Based on the morphologies of H$_2$, Br$gamma$, $^{12}$CO, and 8$mu$m emission, the H$_2$ to Br$gamma$ line ratio and Cloudy models, we find that the H$_2$ emission is formed inside the photodissociation regions of 30 Doradus, 2 - 3 pc to the ionization front of the HII region, in a relatively low-density environment $<$ 10$^4$ cm$^{-3}$. Comparisons with Br$gamma$, 8$mu$m, and CO emission indicate that H$_2$ emission is due to fluorescence, and provide no evidence for shock excited emission of this line.
While the CO(1-0) transition is often used to deduce the total molecular hydrogen in galaxies, it is challenging to detect in low metallicity galaxies, in spite of the star formation taking place. In contrast, the [CII] 158 micron line is relatively bright, highlighting a potentially important reservoir of H2 that is not traced by CO(1-0), but residing in the C+ - emitting regions. We explore a method to quantify the total H2 mass (MH2) in galaxies and learn what parameters control the CO-dark gas reservoir. We present Cloudy grids of density, radiation field and metallicity in terms of observed quantities, such as [OI], [CI], CO(1-0), [CII], total infrared luminosity and the total MH2 and provide recipes based on these models to derive total MH2 mass estimates from observations. The models are applied to the Herschel Dwarf Galaxy Survey, extracting the total MH2 for each galaxy which is compared to the H2 determined from the observed CO(1-0) line. While the H2 traced by CO(1-0) can be negligible, the [CII] 158 micron line can trace the total H2. 70% to 100% of the total H2 mass is not traced by CO(1-0) in the dwarf galaxies, but is well-traced by [CII] 158 micron line. The CO-dark gas mass fraction correlates with the observed L[CII]/LCO(1-0) ratio. A conversion factor for [CII] luminosity to total H2 and a new CO-to-total-MH2 conversion factor, as a function of metallicity, is presented. A recipe is provided to quantify the total mass of H2 in galaxies, taking into account the CO and [CII] observations. Accounting for this CO-dark H2 gas, we find that the star forming dwarf galaxies now fall on the Schmidt-Kennicutt relation. Their star-forming efficiency is rather normal, since the reservoir from which they form stars is now more massive when introducing the [CII] measures of the total H2, compared to the little amount of H2 in the CO-emitting region.
The most luminous Spitzer point sources in the 30 Doradus triggered second generation are investigated coherently in the 3-8 micron region. Remarkable diversity and complexity in their natures are revealed. Some are also among the brightest JHK sources, while others are not. Several of them are multiple when examined at higher angular resolutions with HST NICMOS and WFPC2/WFC3 as available, or with VISTA/VMC otherwise. One is a dusty compact H II region near the far northwestern edge of the complex, containing a half dozen bright I-band sources. Three others appear closely associated with luminous WN stars and causal connections are suggested. Some are in the heads of dust pillars oriented toward R136, as previously discussed from the NICMOS data. One resides in a compact cluster of much fainter sources, while another appears monolithic at the highest resolutions. Surprisingly, one is the brighter of the two extended mystery spots associated with Knot 2 of Walborn et al. Masses are derived from YSO models for unresolved sources and lie in the 10-30 M_{sun} range. Further analysis of the IR sources in this unique region will advance understanding of triggered massive star formation, perhaps in some unexpected and unprecedented ways.