F. P. Israel

Black hole accretion and star formation as drivers of gas excitation and chemistry in Markarian 231

We present a full high resolution SPIRE FTS spectrum of the nearby ultraluminous infrared galaxy Mrk 231. In total 25 lines are detected, including CO J = 5-4 through J = 13-12, 7 rotational lines of H2O, 3 of OH+ and one line each of H2O+, CH+, and HF. We find that the excitation of the CO rotational levels up to J = 8 can be accounted for by UV radiation from star formation. However, the approximately flat luminosity distribution of the CO lines over the rotational ladder above J = 8 requires the presence of a separate source of excitation for the highest CO lines. We explore X-ray heating by the accreting supermassive black hole in Mrk 231 as a source of excitation for these lines, and find that it can reproduce the observed luminosities. We also consider a model with dense gas in a strong UV radiation field to produce the highest CO lines, but find that this model strongly overpredicts the hot dust mass in Mrk 231. Our favoured model consists of a star forming disk of radius 560 pc, containing clumps of dense gas exposed to strong UV radiation, dominating the emission of CO lines up to J = 8. X-rays from the accreting supermassive black hole in Mrk 231 dominate the excitation and chemistry of the inner disk out to a radius of 160 pc, consistent with the X-ray power of the AGN in Mrk 231. The extraordinary luminosity of the OH+ and H2O+ lines reveals the signature of X-ray driven excitation and chemistry in this region.

The Spitzer Survey of the Small Magellanic Cloud: S3MC Imaging and Photometry in the Mid- and Far-Infrared Wave Bands

We present the initial results from the Spitzer Survey of the Small Magellanic Cloud (S^3MC), which imaged the star-forming body of the SMC in all seven MIPS and IRAC wave bands. We find that the F_8/F_(24) ratio (an estimate of PAH abundance) has large spatial variations and takes a wide range of values that are unrelated to metallicity but anticorrelated with 24 μm brightness and F_(24)/F_(70) ratio. This suggests that photodestruction is primarily responsible for the low abundance of PAHs observed in star-forming low-metallicity galaxies. We use the S3MC images to compile a photometric catalog of ~400,000 mid- and far-infrared point sources in the SMC. The sources detected at the longest wavelengths fall into four main categories: (1) bright 5.8 μm sources with very faint optical counterparts and very red mid-infrared colors ([5.8] - [8.0] > 1.2), which we identify as YSOs; (2) bright mid-infrared sources with mildly red colors (0.16 ≾ [5.8] - [8.0] < 0.6), identified as carbon stars; (3) bright mid-infrared sources with neutral colors and bright optical counterparts, corresponding to oxygen-rich evolved stars; and (4) unreddened early B stars (B3-O9) with a large 24 μm excess. This excess is reminiscent of debris disks and is detected in only a small fraction of these stars (≾5%). The majority of the brightest infrared point sources in the SMC fall into groups 1-3. We use this photometric information to produce a catalog of 282 bright YSOs in the SMC with a very low level of contamination (~7%).

THE SPITZER SURVEY OF THE SMALL MAGELLANIC CLOUD: FAR-INFRARED EMISSION AND COLD GAS IN THE SMALL MAGELLANIC CLOUD

We present new FIR maps of the SMC at 24, 70, and 160 μm obtained as part of the Spitzer Survey of the Small Magellanic Cloud (S^(3)MC). These maps cover most of the star formation in the SMC bar and wing. We combine our maps with literature data to derive the dust mass surface density across the SMC. We find a total dust mass of M_(dust) = 3 × 10^5 M_☉, implying a dust-to-hydrogen ratio over the region studied of log_(10)(D/H) = -2.86, or 1 : 700, which includes H_2. Assuming the dust to trace the total gas column, we derive H_2 surface densities across the SMC. We find a total H_2 mass M_(H_2) = 3.2 × 10^7 M_☉ in a distribution similar to that of the CO, but more extended. We compare profiles of CO and H_2 around six molecular peaks; on average H_2 is more extended than CO by a factor of ~1.3. The implied CO-to-H_2 conversion factor over the whole SMC is X_(CO) = (13 ± 1) × 10^(21) cm^(-2) (K km s^(-1))^(-1). Over the volume occupied by CO the conversion factor is lower, X_(CO) = (6 ± 1) × 10^(21) cm^(-2) (K km s^(-1))^(-1), but still a few times larger than that found using virial mass methods. The molecular peaks have H_2 surface densities Σ_(H_2) ≈ 180 ± 30 M pc^(-2), similar to those in Milky Way GMCs, and correspondingly low extinctions, A_V ~ 1-2 mag. The theory of photoionization-regulated star formation predicts A_V ~ 6, which would require the GMCs to be ~3 times smaller than our 46 pc resolution element. For a given hydrostatic gas pressure, the SMC has a 2-3 times lower ratio of molecular to atomic gas than spiral galaxies. Combined with lower mean densities, this results in this galaxy having only 10% of its gas in the molecular phase.

Surveying the Agents of Galaxy Evolution in the Tidally Stripped, Low Metallicity Small Magellanic Cloud (SAGE-SMC). I. Overview

The Small Magellanic Cloud (SMC) provides a unique laboratory for the study of the lifecycle of dust given its low metallicity (~1/5 solar) and relative proximity (~60 kpc). This motivated the SAGE-SMC (Surveying the Agents of Galaxy Evolution in the Tidally Stripped, Low Metallicity Small Magellanic Cloud) Spitzer Legacy program with the specific goals of studying the amount and type of dust in the present interstellar medium, the sources of dust in the winds of evolved stars, and how much dust is consumed in star formation. This program mapped the full SMC (30 deg^2) including the body, wing, and tail in seven bands from 3.6 to 160 μm using IRAC and MIPS on the Spitzer Space Telescope. The data were reduced and mosaicked, and the point sources were measured using customized routines specific for large surveys. We have made the resulting mosaics and point-source catalogs available to the community. The infrared colors of the SMC are compared to those of other nearby galaxies and the 8 μm/24 μm ratio is somewhat lower than the average and the 70 μm/160 μm ratio is somewhat higher than the average. The global infrared spectral energy distribution (SED) shows that the SMC has approximately 1/3 the aromatic emission/polycyclic aromatic hydrocarbon abundance of most nearby galaxies. Infrared color-magnitude diagrams are given illustrating the distribution of different asymptotic giant branch stars and the locations of young stellar objects. Finally, the average SED of H II/star formation regions is compared to the equivalent Large Magellanic Cloud average H II/star formation region SED. These preliminary results will be expanded in detail in subsequent papers.

HERschel Inventory of The Agents of Galaxy Evolution (HERITAGE): The Large Magellanic Cloud dust

The HERschel Inventory of The Agents of Galaxy Evolution (HERITAGE) of the Magellanic Clouds will use dust emission to investigate the life cycle of matter in both the Large and Small Magellanic Clouds (LMC and SMC). Using the Herschel Space Observatory’s PACS and SPIRE photometry cameras, we imaged a 2° × 8° strip through the LMC, at a position angle of ~22.5° as part of the science demonstration phase of the Herschel mission. We present the data in all 5 Herschel bands: PACS 100 and 160 μm and SPIRE 250, 350 and 500 μm. We present two dust models that both adequately fit the spectral energy distribution for the entire strip and both reveal that the SPIRE 500 μm emission is in excess of the models by ~6 to 17%. The SPIRE emission follows the distribution of the dust mass, which is derived from the model. The PAH-to-dust mass (f_(PAH)) image of the strip reveals a possible enhancement in the LMC bar in agreement with previous work. We compare the gas mass distribution derived from the HI 21 cm and CO J = 1−0 line emission maps to the dust mass map from the models and derive gas-to-dust mass ratios (GDRs). The dust model, which uses the standard graphite and silicate optical properties for Galactic dust, has a very low GDR = 65^(+15) _(−18) making it an unrealistic dust model for the LMC. Our second dust model, which uses amorphous carbon instead of graphite, has a flatter emissivity index in the submillimeter and results in a GDR = 287^_(+25)_(−42) that is more consistent with a GDR inferred from extinction.

DUST HEATING SOURCES IN GALAXIES: THE CASE OF M33 (HERM33ES)

Dust emission is one of the main windows to the physics of galaxies and to star formation as the radiation from young, hot stars is absorbed by the dust and reemitted at longer wavelengths. The recently launched Herschel satellite now provides a view of dust emission in the far-infrared at an unequaled resolution and quality up to 500 mu m. In the context of the Herschel HERM33ES open time key project, we are studying the moderately inclined Scd local group galaxy M33 which is located only 840 kpc away. In this article, using Spitzer and Herschel data ranging from 3.6 mu m to 500 mu m, along with H I, H alpha maps, and Galaxy Evolution Explorer ultraviolet data, we have studied the emission of the dust at the high spatial resolution of 150 pc. Combining Spitzer and Herschel bands, we have provided new, inclination-corrected, resolved estimators of the total infrared brightness and of the star formation rate from any combination of these bands. The study of the colors of the warm and cold dust populations shows that the temperature of the former is, at high brightness, dictated by young massive stars but, at lower brightness, heating is taken over by the evolved populations. Conversely, the temperature of the cold dust is tightly driven by the evolved stellar populations.

The James Clerk Maxwell Telescope Nearby Galaxies Legacy Survey. I. Star-Forming Molecular Gas in Virgo Cluster Spiral Galaxies

We present large-area maps of the CO J = 3-2 emission obtained at the James Clerk Maxwell Telescope for four spiral galaxies in the Virgo Cluster. We combine these data with published CO J = 1-0, 24 μm, and Hα images to measure the CO line ratios, molecular gas masses, and instantaneous gas depletion times. For three galaxies in our sample (NGC 4254, NGC 4321, and NGC 4569), we obtain molecular gas masses of 7 × 108 – 3 × 109 M ☉ and disk-averaged instantaneous gas depletion times of 1.1-1.7 Gyr. We argue that the CO J = 3-2 line is a better tracer of the dense star-forming molecular gas than the CO J = 1-0 line, as it shows a better correlation with the star formation rate surface density both within and between galaxies. NGC 4254 appears to have a larger star formation efficiency (smaller gas depletion time), perhaps because it is on its first passage through the Virgo Cluster. NGC 4569 shows a large-scale gradient in the gas properties traced by the CO J = 3-2/J = 1-0 line ratio, which suggests that its interaction with the intracluster medium is affecting the dense star-forming portion of the interstellar medium directly. The fourth galaxy in our sample, NGC 4579, has weak CO J = 3-2 emission despite having bright 24 μm emission; however, much of the central luminosity in this galaxy may be due to the presence of a central active galactic nucleus.

Evidence for CO Shock Excitation in NGC 6240 from Herschel SPIRE Spectroscopy

We present Herschel SPIRE FTS spectroscopy of the nearby luminous infrared galaxy NGC 6240. In total 20 lines are detected, including CO J = 4-3 through J = 13-12, 6 H2O rotational lines, and [C I] and [N II] fine-structure lines. The CO to continuum luminosity ratio is 10 times higher in NGC 6240 than Mrk 231. Although the CO ladders of NGC 6240 and Mrk 231 are very similar, UV and/or X-ray irradiation are unlikely to be responsible for the excitation of the gas in NGC 6240. We applied both C and J shock models to the H-2 v = 1-0 S(1) and v = 2-1 S(1) lines and the CO rotational ladder. The CO ladder is best reproduced by a model with shock velocity v(s) = 10 km s(-1) and a pre-shock density n(H) = 5 x 10(4) cm(-3). We find that the solution best fitting the H-2 lines is degenerate. The shock velocities and number densities range between v(s) = 17-47 km s(-1) and n(H) = 10(7)-5x10(4) cm(-3), respectively. The H-2 lines thus need a much more powerful shock than the CO lines. We deduce that most of the gas is currently moderately stirred up by slow (10 km s(-1)) shocks while only a small fraction (less than or similar to 1%) of the interstellar medium is exposed to the high-velocity shocks. This implies that the gas is rapidly losing its highly turbulent motions. We argue that a high CO line-to-continuum ratio is a key diagnostic for the presence of shocks.

The Structure of a Low-metallicity Giant Molecular Cloud Complex

To understand the impact of low metallicities on giant molecular cloud (GMC) structure, we compare far-infrared dust emission, CO emission, and dynamics in the star-forming complex N83 in the Wing of the Small Magellanic Cloud (SMC). Dust emission (measured by Spitzer as part of the Spitzer Survey of the SMC and Surveying the Agents of a Galaxys Evolution in the SMC surveys) probes the total gas column independent of molecular line emission and traces shielding from photodissociating radiation. We calibrate a method to estimate the dust column using only the high-resolution Spitzer data and verify that dust traces the interstellar medium in the H I-dominated region around N83. This allows us to resolve the relative structures of H2, dust, and CO within a GMC complex, one of the first times such a measurement has been made in a low-metallicity galaxy. Our results support the hypothesis that CO is photodissociated while H2 self-shields in the outer parts of low-metallicity GMCs, so that dust/self-shielding is the primary factor determining the distribution of CO emission. Four pieces of evidence support this view. First, the CO-to-H2 conversion factor averaged over the whole cloud is very high 4-11 × 1021 cm–2 (K km s–1)–1, or 20-55 times the Galactic value. Second, the CO-to-H2 conversion factor varies across the complex, with its lowest (most nearly Galactic) values near the CO peaks. Third, bright CO emission is largely confined to regions of relatively high line-of-sight extinction, AV 2 mag, in agreement with photodissociation region models and Galactic observations. Fourth, a simple model in which CO emerges from a smaller sphere nested inside a larger cloud can roughly relate the H2 masses measured from CO kinematics and dust.

Determining dust temperatures and masses in the Herschel era: The importance of observations longward of 200 micron

Context. The properties of the dust grains (e.g., temperature and mass) can be derived from fitting far-IR SEDs (≥100 μm). Only with SPIRE on Herschel has it been possible to get high spatial resolution at 200 to 500 μm that is beyond the peak (~160 μm) of dust emission in most galaxies. Aims. We investigate the differences in the fitted dust temperatures and masses determined using only 200 μm data (new SPIRE observations) to determine how important having >200 μm data is for deriving these dust properties. Methods. We fit the 100 to 350 μm observations of the Large Magellanic Cloud (LMC) point-by-point with a model that consists of a single temperature and fixed emissivity law. The data used are existing observations at 100 and 160 μm (from IRAS and Spitzer) and new SPIRE observations of 1/4 of the LMC observed for the HERITAGE key project as part of the Herschel science demonstration phase. Results. The dust temperatures and masses computed using only 100 and 160 μm data can differ by up to 10% and 36%, respectively, from those that also include the SPIRE 250 & 350 μm data. We find that an emissivity law proportional to λ^(−1.5) minimizes the 100–350 μm fractional residuals. We find that the emission at 500 μm is ~10% higher than expected from extrapolating the fits made at shorter wavelengths. We find the fractional 500 μm excess is weakly anti-correlated with MIPS 24 μm flux and the total gas surface density. This argues against a flux calibration error as the origin of the 500 μm excess. Our results do not allow us to distinguish between a systematic variation in the wavelength dependent emissivity law or a population of very cold dust only detectable at λ ≥ 500 μm for the origin of the 500 μm excess.