format: text

text:
# ISOCAM view of the starburst galaxies M 82, NGC 253, and NGC 1808 Based on observations with ISO, an ESA project with instruments funded by ESA member states (especially the PI countries: France, Germany, the Netherlands, and the United Kingdom), and with participation of ISAS and NASA.

## 1 Introduction

As hosts of vigorous episodes of star formation activity, starburst galaxies in the local Universe offer a unique opportunity to study the star formation process in extreme environments reminiscent of the conditions presumably prevailing in primordial galaxies. Studies in recent years have underscored the prominent role of starbursts in galaxy formation and evolution as well as their significant contribution to the extragalactic background, and to the chemical enrichment and heating of the intergalactic medium up to the highest redshifts (e.g. Steidel et al. Ste96 (1996); Puget et al. Pug96 (1996); Smail et al. Sma97 (1997); Franx et al. Fra97 (1997); Pettini et al. Pet98 (1998); Elbaz et al. Elb02 (2002)). Understanding the starburst phenomenon is thus a key issue of local and cosmological relevance.

Nuclear starburst galaxies generally emit the bulk of their luminosity in the infrared, primarily as reprocessed radiation by interstellar dust grains heated by important populations of hot massive stars. In addition, a substantial fraction – if not most – of the star formation in starburst systems is heavily obscured at optical and ultraviolet wavelengths by large amounts of dust. This has been particularly dramatically illustrated by mid-infrared imaging of the colliding pair NGC 4038/4039 (Mirabel et al. Mir98 (1998)). Studies at infrared wavelengths are thus crucial in defining the properties of starburst galaxies and the interplay between the primary energy sources and the reprocessing material. The growing evidence for the presence and significant role of dust in distant galaxies further emphasizes the need for infrared investigations of nearby dust-rich templates (e.g. Rowan-Robinson Row01 (2001); Blain et al. Bla02 (2002); and references therein). However, except for a limited sample, the detailed spatial and spectral energy distributions of starburst galaxies in this regime are still poorly known. Progress has been hindered by instrumental limitations, and by the difficulties inherent to ground-based and air-borne observations at these wavelengths.

The Infrared Space Observatory (ISO; Kessler et al. Kes96 (1996)) has revolutionized the field of infrared astronomy (see the reviews by Cesarsky & Sauvage Ces99 (1999); Genzel & Cesarsky Gen00 (2000)), with unprecedented capabilities compared to its predecessor the Infrared Astronomical Satellite (IRAS). In the mid-infrared (MIR, $`\lambda =5.016.5\mu \mathrm{m}`$), the camera ISOCAM (Cesarsky et al. Ces96 (1996)) provided an improvement in sensitivity and spatial resolution by factors of $`1000`$ and $`60`$, respectively. In addition to broad- and narrow-band filters, ISOCAM featured continuously variable filters (CVFs) allowing arcsecond-scale spectrophotometric imaging with full wavelength coverage at a resolution of $`R\lambda /\mathrm{\Delta }\lambda 40`$. The MIR regime contains a variety of key features tracing different components of the interstellar medium (ISM) and is subject to relatively little extinction by dust ($`A_{517\mu \mathrm{m}}<0.1A_V`$). Spectrophotometric observations permit studies of the most prominent features to characterize the spatial distribution and infer the properties of the ISM and exciting sources, deep into obscured star-forming regions.

In this paper, we present results of ISOCAM CVF observations of M 82, NGC 253, and NGC 1808. M 82 and NGC 253 are two of the nearest and brightest starburst galaxies (at distances of $`D=3.3`$ and 2.5 Mpc, respectively, $`1^{\prime \prime }=16`$ and 12 pc; Freedman & Madore Fre88 (1988); Davidge & Pritchet Dav90 (1990)). Both have been extensively studied in all accessible spectral regions and have often been considered as prototypical objects for the starburst phenomenon (see Telesco Tel88 (1988), Rieke et al. Rie93 (1993), and Engelbracht et al. Eng98 (1998) for reviews). NGC 1808 is a more distant system ($`D=10.9\mathrm{Mpc}`$ for $`H_0=75\mathrm{km}\mathrm{s}^1\mathrm{Mpc}^1`$, $`1^{\prime \prime }=53\mathrm{pc}`$; Sandage & Tammann San87 (1987)) which closely resembles M 82 and NGC 253 in its global properties (e.g. Dahlem et al. Dah90 (1990); Junkes et al. Jun95 (1995)). In all three galaxies, starburst activity takes place within the central $`0.11\mathrm{kpc}`$ regions which are severely obscured at optical and ultraviolet wavelengths. The large extinction levels likely result in part from high disk inclination angles of $`80^{}`$, $`78^{}`$, and $`57^{}`$ for M 82, NGC 253, and NGC 1808 (Götz et al. Got90 (1990); Pence Pen81 (1981); Reif et al. Rei82 (1982)).

At comparable infrared luminosities111as computed from the IRAS fluxes following the prescription of Sanders & Mirabel (San96 (1996)). $`L_{\mathrm{IR}}L_{81000\mu \mathrm{m}}`$ of $`3.7\times 10^{10}`$, $`1.4\times 10^{10}`$, and $`3.4\times 10^{10}\mathrm{L}_{}`$ for M 82, NGC 253, and NGC 1808, these galaxies also share key morphological features including a kiloparsec-scale stellar bar, large amounts of molecular gas in the central $`1\mathrm{kpc}`$ predominantly concentrated in circumnuclear ring-, spiral-, or bar-like structures, populations of compact and luminous “super star clusters” as well as of compact non-thermal radio sources associated with supernova remnants, and a large-scale outflowing starburst wind (Förster Schreiber et al. FS01 (2001); Engelbracht et al. Eng98 (1998); Tacconi-Garman et al. Tac96 (1996); and references therein). The triggering of starburst activity in M 82 is generally attributed to gravitational interaction as evidenced by the extended filamentary tidal features in H i gas distribution threading M 82 and its neighbours M 81 and NGC 3077 or, alternatively, to the stellar bar which may itself have been induced by the interaction (e.g. Yun et al. Yun93 (1993, 1994)). In NGC 253, the bar itself appears to be the primary mechanism responsible for the onset of starburst activity as there are no obvious signatures of interaction with a nearby companion, although a past merger or accretion event involving a small galaxy has been suggested based on kinematical evidence (e.g. Anantharamaiah & Goss Ana96 (1996); Böker et al. Bok98 (1998)). Tidal interaction with NGC 1792 has been proposed for NGC 1808 based on circumstantial evidence (e.g. Dahlem et al. Dah90 (1990); Koribalski et al. Kor93 (1993)). The recent high sensitivity H i observations by Dahlem et al. (Dah01 (2001)) do not indicate the presence of any tidal feature, dwarf galaxy or other debris that could support a close passage in the past, leaving the stellar bar as more probable trigger or a possible merger/accretion event.

In view of their similar nature, we focus on the comparison of the spatial and spectral properties of M 82, NGC 253, and NGC 1808 at MIR wavelengths. Though small, this sample also allows us to probe the transition regime between normal spiral and irregular galaxies and the more extreme luminous and ultraluminous infrared galaxies (LIRGs and ULIRGs for which $`L_{\mathrm{IR}}>10^{11}`$ and $`>10^{12}\mathrm{L}_{}`$, respectively). Section 2 describes the observations and data reduction procedure, and Sect. 3 presents the results. Section 4 discusses the origin and spatial distibution of the continuum and emission features. Section 5 addresses the issues of extinction effects, variations in spectral properties, and indicators of star formation activity. Section 6 summarizes the paper.

## 2 Observations and data reduction

We observed M 82, NGC 253, and NGC 1808 with ISOCAM (Cesarsky et al. Ces96 (1996)) on board ISO (Kessler et al. Kes96 (1996)) as part of the guaranteed time program CAMACTIV (P.I. I. F. Mirabel). We used the CVFs to cover the entire $`\lambda =5.016.5\mu \mathrm{m}`$ range with a resolution of $`R\lambda /\mathrm{\Delta }\lambda 35`$ to 45 and increments $`\delta \lambda =0.06`$ to 0.1$`\mu `$m from short to long wavelengths. For all three sources, twelve single-frame exposures of 2.1 s were recorded per wavelength channel, for total integration times of about $`1^\mathrm{h}05^\mathrm{m}`$. We selected the $`3^{\prime \prime }\mathrm{pixel}^1`$ scale for M 82 and NGC 1808, and the $`1.5^{\prime \prime }\mathrm{pixel}^1`$ scale for NGC 253; total fields of view of $`96^{\prime \prime }\times 96^{\prime \prime }`$ and $`48^{\prime \prime }\times 48^{\prime \prime }`$ are thus covered by the $`32\times 32`$ pixels detector array. Table 1 presents the details of the observations. The areas observed with ISOCAM cover almost entirely the MIR sources in each galaxy (Sect. 4) and are indicated on optical images in Fig. 1. This emphasizes the small MIR source sizes compared to the optical extent of the galaxies and that most of the MIR emission originates in their nuclear regions.

We reduced the data using the CAM Interactive Analysis package (CIA)222CIA is a joint development by the ESA Astrophysics Division and the ISOCAM consortium led by the ISOCAM P.I. C. Cesarsky, Direction des Sciences de la Matière, C.E.A. Saclay, France. following the procedures described in the ISOCAM User’s Manual (Delaney Del97 (1997)). We first subtracted the dark current using a model of the secular evolution of ISOCAM’s dark current (Biviano et al. Biv97 (1997)), removed cosmic ray hits by applying a multiresolution median filtering technique (Starck et al. Sta99 (1999)), and corrected for transient detector memory effects using the “Fouks-Schubert” algorithm (Coulais & Abergel Cou00 (2000)). We then combined the individual exposures for each wavelength channel with the help of a shift-and-add algorithm accounting for the satellite jitter motions (typical amplitude of $`1^{\prime \prime }`$) and performed absolute flux calibration based on the sensitivity calibration files provided in CIA version 3.0.

For proper comparison of fluxes and maps obtained at different wavelengths, we smoothed the CVF channels to a common angular resolution and registered them appropriately. The point spread function (PSF) varies throuhgout the CVF scans, with theoretical full-width at half-maximum $`\mathrm{FWHM}=3.1^{\prime \prime }5.6^{\prime \prime }`$ for the $`3^{\prime \prime }\mathrm{pixel}^1`$ scale and $`2.0^{\prime \prime }5.2^{\prime \prime }`$ for the $`1.5^{\prime \prime }\mathrm{pixel}^1`$ scale, over the range $`5.516.5\mu \mathrm{m}`$. We convolved the data to a final resolution of 5.6″ for M 82 and NGC 1808, and 5.2″ for NGC 253, corresponding to a linear resolution of 90, 297, and 62 pc, respectively. We computed the transfer functions (close to Gaussians in profile) using the library of theoretical PSFs for ISOCAM and a CLEAN algorithm. The small-scale structure is thus smeared out in the higher resolution channels but no important spatial feature is lost.

The source position varies slightly with wavelength due to telescope motions and changes in optical path, with the largest offset occurring near 9.2$`\mu `$m between the short and long wavelength segments of the CVFs. Since no unresolved point source is detected in the data, we determined the position offsets from cross-correlation over the entire field of view, taking as references images integrated over various wavelength intervals. Effects due to the different spatial distributions of emission features were easily identified and ignored in the determination of the instrumental shifts. We fitted straight lines to the cross-correlation offsets for each CVF segment. For the $`3^{\prime \prime }\mathrm{pixel}^1`$ scale, the horizontal, or $`x`$-axis offsets relative to 9.2$`\mu `$m in the short and long wavelength segments did not exceed 0.3 pixel and 0.5 pixel, respectively, with a jump of 0.7 pixel at the transition wavelength. Shifts along the $`y`$-axis were negligible throughout the entire wavelength range ($`<0.2\mathrm{pixel}`$). For the $`1.5^{\prime \prime }\mathrm{pixel}^1`$ scale, both $`x`$\- and $`y`$-axis offsets were $`0.3\mathrm{pixel}`$ within each CVF segment, with jumps at 9.2$`\mu `$m of 0.2 pixel only.

We did not attempt to correct for instrumental flat-field, straylight, and ghost images generated by multiple reflections between the detector and filters since appropriate calibration files for the CVF mode were not available. Based on theoretical models for point-like and extended sources as well as observations of stars and of the zodiacal light (Biviano et al. 1998a ; 1998b ; Okumura Oku00 (2000)), the combined flat-field, straylight, and ghost effects integrated over the entire detector array vary between $`35\%`$ and 10% of the source flux from 5 to 16$`\mu `$m, and exhibit qualitatively similar spatial structures at different wavelengths. The r.m.s. noise estimated from existing flat-field images obtained with ISOCAM’s broad-band filters is $`10\%`$ for both pixel scales (Biviano et al. 1998b ).

Based on available calibration accuracy reports (Blommaert & Cesarsky Blo98 (1998) and references therein), the systematic errors of the data processing and flux calibration described above sum up to total uncertainties on the absolute photometry of $`50\%`$. For our analysis, we will consider formal effective uncertainties on the relative fluxes estimated as follows. We measured the dispersion around the mean flux among the individual reduced exposures for each wavelength channel at each detector pixel, included an additional 10% r.m.s. noise for intrinsic pixel-to-pixel variations in sensitivity, and accounted for uncertainties in relative spectral calibration. For M 82 and NGC 253, the latter were evaluated from the comparison with data from the ISO Short Wavelength Spectrometer (SWS; de Graauw et al. deG96 (1996)) in Sect. 3, amounting to 5% and 11%, respectively. For NGC 1808, very low flux levels and high noise around 10$`\mu `$m, below 6$`\mu `$m, and outside of the main emission region led to poor correction for transient memory effects because of the difficulty in determining the initial and stabilized fluxes and of the bad behaviour of the algorithm wherever fluxes become negative. In this case, we estimated the errors in relative spectral calibration by comparing results when applying the transient correction to the data as they were and to the data where all values lower than 0.01 were set to 0.01, adding a conservative 10%. The median and average differences are of 16% and a factor of 2. Experimentation showed that the choice of the threshold is of little consequence as long as it is sufficiently smaller than the signal where the transient correction is satisfactory.

Overall, the formal effective uncertainties per pixel and wavelength channel within the brighter emission regions of interest (defined by the “starburst cores” in Sect. 3) have median values of 18%, 29%, and 45% for M 82, NGC 253, and NGC 1808. Computed over the entire ISOCAM field of view, the medians increase to 32%, 88%, and a factor of 2.6 as a result of the higher noise towards the edges of the array where little flux is detected.

## 3 Results

### 3.1 Mid-infrared spectra and feature identification

We extracted from the reduced ISOCAM data cubes spectra integrated over various regions: a “starburst core” centered on each galaxy’s nucleus333as defined by the peak of the near-infrared 2$`\mu `$m emission; see table 2 for coordinates and references. and covering most of the emission source, and a smaller region at the position of the peak observed for the integrated MIR emission. For M 82 and NGC 253, we further included a “disk” corresponding to an annulus outside of the starburst core and the ISOCAM field of view. Due to the higher noise level in the NGC 1808 data, only the starburst core and MIR peak spectra are of sufficient quality for analysis. The disk, core, and MIR peak regions were chosen to sample relatively quiescent to intense star-forming activity as traced by e.g. the 15$`\mu `$m continuum emission (see Sects. 3.2 and 5.2). The starburst cores correspond to the apertures generally used in previous studies when referring to the global properties of the starburst in each galaxy.

Figure 2 indicates the selected regions on maps of the $`14.815.2\mu \mathrm{m}`$ emission. Figure 3 shows the spectra, all normalized to a total $`6.06.6\mu \mathrm{m}`$ flux density of unity to facilitate comparison of the relative strength of the emission features at $`\lambda 11\mu \mathrm{m}`$ and of the continuum intensity at longer wavelengths. The formal effective uncertainties are plotted along with the spectra. Table 2 gives the parameters for the synthetic apertures, the normalizing fluxes, and the average and median uncertainties for each spectrum. Figure 3 also shows for M 82 and NGC 253 the average spectrum of all individual resolution elements within the ISOCAM field of view, normalized as described above, together with the dispersion around the mean, the full range observed, and the typical (median) uncertainties. The resolution elements for M 82 and NGC 253 correspond to rebinned pixels of size $`2\times 2`$ and $`4\times 4`$ in original detector pixels, respectively.

All spectra look very similar and exhibit the “classical” characteristics observed towards star-forming regions and galaxies: conspicuous broad emission features (the “unidentified infrared bands” or UIBs), a featureless continuum rising importantly at $`\lambda 11\mu \mathrm{m}`$, and an apparent dip near 10$`\mu `$m. The UIBs are attributed to a family of particles, the nature of which still is debated, stochastically heated by single ultraviolet photons while the long-wavelength continuum is ascribed to very small dust grains between the transient heating and thermal equilibrium regimes depending on grain properties and radiation field intensity (e.g. Léger et al. Leg89 (1989); Allamandola et al. All89 (1989); Désert et al. Des90 (1990); Duley & Williams Dul91 (1991); Tielens et al. Tie99 (1999); see also reviews by Puget & Léger Pug89 (1989); Cesarsky & Sauvage Ces99 (1999); Genzel & Cesarsky Gen00 (2000)). We will hereafter refer to these ISM components as “PAHs,” adopting the currently popular model in which the UIB carriers consist of polycyclic aromatic hydrocarbon molecules, and “VSGs.”

A number of weaker emission features are also detected in the spectra of Fig. 3. However, their identification is problematic at such low spectral resolution due to possible feature blends. To emphasize this point and secure the identifications, Fig. 4 shows the spectra of M 82 and NGC 253 obtained at $`R5001000`$ with the ISO-SWS (from Sturm et al. Stu00 (2000) and Förster Schreiber et al. FS01 (2001)). The SWS data are also compared with ISOCAM spectra taken in the same apertures, after convolution to the same spectral resolution. The ISOCAM and SWS data agree very well, confirming the accuracy of the absolute and relative flux calibration for both instruments. For M 82, the differences are 5% on average (24% at most). For NGC 253, they are of 11% on average (30% at most) at $`\lambda 5.4\mu \mathrm{m}`$ while they reach a factor of two at shorter wavelengths, probably due to residual transient memory effects.

From this comparison (see also Sturm et al. Stu00 (2000)), the features at 6.2, 7.7, 8.6, and 11.3$`\mu `$m are unambiguously identified with PAH emission. The broad feature near 12.7$`\mu `$m clearly results from the blending of the PAH 12.7$`\mu `$m band and of the \[Ne ii\] 12.81$`\mu `$m fine-structure line. The \[Ne iii\] 15.56$`\mu `$m line is also blended with the nearby PAH 15.7$`\mu `$m feature. From the SWS data, and with the continuum and integration bandpasses given below, the PAH 12.7$`\mu `$m accounts for about 50% of the flux in the 12.7$`\mu `$m blend for both galaxies. The PAH 15.7$`\mu `$m contains about 30% of the flux in the 15.6$`\mu `$m blend for M 82 (the noisier SWS spectrum of NGC 253 makes an estimate difficult). The least contaminated fine-structure line detected with ISOCAM is \[Ar ii\] 6.99$`\mu `$m, with the underlying PAH at 7.0$`\mu `$m and the $`\mathrm{H}_2`$ $`00S(5)`$ rotational line at 6.91$`\mu `$m contributing $`20\%`$ and 5%, respectively, to the blend flux in both galaxies. The weak features at 5.65, 13.55, and 14.25$`\mu `$m are identified with PAH bands; the latter is definitely not due to the high excitation \[Ne v\] 14.32$`\mu `$m line which is undetected in the SWS data. The final identifications in the ISOCAM spectra are given in Fig. 3. Because of the excellent correspondence between features seen in all three galaxies, the identifications for M 82 and NGC 253 are adopted for NGC 1808 as well.

We performed various continuum and feature measurements on the spectra. Table 3 gives the intervals used in the computations and tables 4 and 5 report the results. The flux densities in the $`6.09.0\mu \mathrm{m}`$ and $`13.515.0\mu \mathrm{m}`$ bands are dominated by emission from PAHs and VSGs, respectively. We chose the narrow continuum bands centered at 5.5 and 15.0$`\mu `$m so as to minimize the contribution from PAHs and other emission lines. In particular, we obtain the same 5.5$`\mu `$m flux densities within $`5\%`$ from the ISOCAM data and from the higher resolution SWS spectra for the SWS field of view in both M 82 and NGC 253, indicating negligible contribution from the adjacent PAH 5.65$`\mu `$m feature and \[Fe ii\] 5.34$`\mu `$m line in the lower resolution data. We measured the fluxes in the PAH features at 6.2, 7.7, 8.6, and 11.3$`\mu `$m, in the \[Ar ii\] 6.99$`\mu `$m line, and in the PAH 12.7$`\mu `$m $`+`$ \[Ne ii\] 12.81$`\mu `$m and PAH 15.7$`\mu `$m $`+`$ \[Ne iii\] 15.56$`\mu `$m blends by integrating the flux under the feature profiles after subtracting a continuum baseline obtained by linear interpolation between adjacent spectral intervals. More sophisticated methods such as profile fitting (e.g. Uchida et al. Uch00 (2000)) are not necessary for our purposes and are difficult to apply to the lower signal-to-noise (S/N) ratio data of individual pixels in generating linemaps (Sect. 3.2).

We did not compute fluxes for the individual \[Ne ii\] 12.81$`\mu `$m and \[Ne iii\] 15.56$`\mu `$m lines. The wavelength sampling is too coarse for reliable profile decomposition, with the feature peaks of each blend sampled by adjacent wavelength channels. In addition, our attempts to subtract the PAH contribution by attributing an excess in the blend 12.7$`\mu `$m/PAH 11.3$`\mu `$m ratio to the \[Ne ii\] 12.81$`\mu `$m line proved too sensitive to the definition of “pure” PAH 12.7$`\mu `$m/11.3$`\mu `$m ratio (e.g. as measured outside of the starburst cores where comparatively little fine-structure line emission from H ii regions is expected). Complications further arise from possible intrinsic variations in the PAH ratios, extinction effects ($`A_{12.7\mu \mathrm{m}}/A_{11.3\mu \mathrm{m}}0.55`$), and unconstrained fine-structure line emission from disk H ii regions. The fluxes for \[Ar ii\] 6.99$`\mu `$m are much more reliable because possible contributions by other features in our sample galaxies are substantially smaller, as mentioned above.

### 3.2 Mid-infrared images: mapping the features

We obtained broad- and narrow-band images as well as maps of the PAHs and \[Ar ii\] 6.99$`\mu `$m line emission from the ISOCAM data cubes by applying to each pixel the procedures described above for the spectra. Figures 5, 6, and 7 present selected images and ratio maps for M 82, NGC 253, and NGC 1808. Contours corresponding to the same levels relative to the peak intensity are plotted for all continuum and emission feature maps for ease of comparison by visual inspection. A $`3\sigma `$ contour is also shown to delineate regions where observed small-scale structures are reliable. In the following, we describe the various maps source by source; their interpretation will be discussed in subsequent sections.

#### 3.2.1 M 82

The continuum and emission feature maps in M 82 show a globally smooth spatial distribution, centered and peaking roughly 5″ southwest of the nucleus. The PAH and 5.5$`\mu `$m continuum emission are the most extended and symmetric, with disk-like isophotes elongated along the galactic plane ($`\mathrm{P}.\mathrm{A}.70^{}`$). In contrast, the 15$`\mu `$m continuum and \[Ar ii\] 6.99$`\mu `$m line emission have more compact distributions which are more asymmetric relative to the nucleus. The 15$`\mu `$m/5.5$`\mu `$m ratio map outlines well the difference in peak morphology between the short and long wavelength continuum. The \[Ar ii\] 6.99$`\mu `$m distribution is the most compact, with centroid (determined from the emission out to a radius of 25″) displaced 3.5″ southwest of that of the PAH and continuum emission, and showing only a slight extension towards the east.

Our ISOCAM maps provide an important complementary dataset to existing MIR images in the literature, which were mostly obtained in different bands or over limited regions (although with higher angular resolution up to $`1^{\prime \prime }`$). The distributions observed in our PAH maps and for the PAH 3.29$`\mu `$m feature by Normand et al. (Nor95 (1995)) and Satyapal et al. (Sat95 (1995)) at $`1^{\prime \prime }`$ resolution are consistent with each other. Our 15$`\mu `$m continuum map and the 19.2$`\mu `$m image of Telesco et al. (Tel91 (1991)) are similar. Maps of the $`N`$-band (10.8$`\mu `$m) and $`11.813.0\mu \mathrm{m}`$ emission generated from the ISOCAM data cubes agree well with those of Telesco et al. (Tel91 (1991)) and Telesco & Gezari (Tel92 (1992)) within the regions covered by the latter two images.

Our \[Ar ii\] 6.99$`\mu `$m linemap globally resembles those of other tracers of ionized gas from H ii regions at mid- and near-infrared wavelengths such as \[Ne ii\] 12.81$`\mu `$m, \[Ar iii\] 8.99$`\mu `$m, \[S iv\] 10.51$`\mu `$m, Br$`\alpha `$, Br$`\gamma `$, and Pa$`\beta `$ (Achtermann & Lacy Ach95 (1995); Satyapal et al. Sat95 (1995)). At a resolution of $`1^{\prime \prime }2^{\prime \prime }`$, these maps reveal a rich sub-structure dominated by prominent sources $`6^{\prime \prime }`$ and 12″ southwest of the nucleus and $`6^{\prime \prime }`$ to the northeast (labeled W1, W2, and E1 by Achtermann & Lacy Ach95 (1995)). While their intensity ratio depends somewhat on the emission line considered, W1 and W2 are together about three times brighter than E1. Radial velocity data of the \[Ne ii\] 12.81$`\mu `$m and Br$`\gamma `$ emission (Larkin et al. Lar94 (1994); Achtermann & Lacy Ach95 (1995)) suggest that most sources reside in a nearly edge-on rotating ring at radius coinciding in projection with W1 and E1 as well as along the stellar bar at larger radii, where the most recent starburst episode took place about 5 Myr ago (e.g. Förster Schreiber et al. FS02 (2002)). Within the positional uncertainties and resolution limitations, the spatial distribution of our \[Ar ii\] 6.99$`\mu `$m map peaks between W1 and W2 and encompasses E1, and thus traces well the youngest starburst regions.

Little differences are seen between the PAH maps but the PAH 6.2$`\mu `$m/7.7$`\mu `$m and PAH 8.6$`\mu `$m/11.3$`\mu `$m ratio maps reveal spatial variations at the 20% and 60% level, respectively. The variations are significant within the brighter emission regions along the disk, where the relative uncertainties are $`<15\%`$. Structures seen towards the map edges are much less reliable as they become comparable in amplitude to the uncertainties. The PAH 6.2$`\mu `$m/7.7$`\mu `$m ratio is lower along the disk and reaches minima on each side of the nucleus. The overall morphology appears to curve northwards away from the nucleus and extensions are hinted at above and below the galactic plane. Similar shape and spurs are observed notably in the large-scale distribution of the molecular and ionized gas line emission and of the radio continuum emission (e.g. Shen & Lo She95 (1995); Achtermann & Lacy Ach95 (1995); Wills et al. Wil99 (1999)). The spatial variations in the PAH 8.6$`\mu `$m/11.3$`\mu `$m ratio match roughly those of the PAH 6.2$`\mu `$m/7.7$`\mu `$m ratio, with higher values along the disk and maxima flanking the nucleus. Noticeably, the western peak lies closer to the nucleus than the western PAH 6.2$`\mu `$m/7.7$`\mu `$m minimum, and the apparent curving and large-scale extensions have no counterpart in the PAH 8.6$`\mu `$m/11.3$`\mu `$m map.

Figure 8 compares the PAH ratio maps with the CO $`J:10`$ millimetric emission of Shen & Lo (She95 (1995)) convolved at the resolution of the ISOCAM maps. The overall correspondence of the PAH 6.2$`\mu `$m/7.7$`\mu `$m minima and PAH 8.6$`\mu `$m/11.3$`\mu `$m maxima with the peaks in CO emission, as well as the curved shape and northeastern extension for the PAH 6.2$`\mu `$m/7.7$`\mu `$m ratio, is quite striking. We believe that the observed variations in PAH ratios are mostly real. Characteristic patterns expected for ghosts are not seen (unresolved ring- or arc-like features most prominent in the presence of point-like sources). Artifacts due to the uncorrected flat field and straylight could produce extended lobes on each side of an axis at $`\mathrm{P}.\mathrm{A}.150^{}`$ for the M 82 data, i.e. roughly the minor axis (Biviano et al. 1998a ; 1998b ; Okumura Oku00 (2000)). However, such lobes would have a much larger extent than the features seen in our PAH ratio maps and the chromatic dependence between 6.2 and 7.7$`\mu `$m, and 8.6 and 11.3$`\mu `$m is predicted and observed to be smaller than the measured variations at 20% and 60% levels, respectively.

#### 3.2.2 NGC 253

In NGC 253, the emission in the continuum bands, PAH features, and \[Ar ii\] 6.99$`\mu `$m line is characterized by a strong peak within $`1^{\prime \prime }2^{\prime \prime }`$ of the nucleus, embedded in a diffuse envelope elongated along the galactic plane ($`\mathrm{P}.\mathrm{A}.50^{}`$). The 15$`\mu `$m continuum and \[Ar ii\] 6.99$`\mu `$m line distributions do not seem to extend as far as the PAH emission in the outer parts of the source. The noisy $`\lambda <6\mu \mathrm{m}`$ channels prevent reliable assessment of the lower level, large-scale 5.5$`\mu `$m continuum emission. The centroids in the various images are essentially indistinguishable (differences $`<0.5\mathrm{pixel}`$).

Though the images reveal little spatial structure because of the intrinsically small source size and the limited angular resolution, they are consistent with expectations at this resolution based on previously published MIR maps (mostly obtained at $`1^{\prime \prime }`$ resolution). These include broad-band images at 10.8 ($`N`$ band) and 19.2$`\mu `$m, narrow-band images at 8.5, 10.0, 12.5, and 20.2$`\mu `$m tracing PAH and/or continuum emission, maps of the 3.29 and 11.3$`\mu `$m PAH emission, and images of the \[Ne ii\] 12.81$`\mu `$m line and of the underlying continuum $`+`$ PAH emission (Piña et al. Pin92 (1992); Telesco et al. Tel93 (1993); Keto et al. Ket93 (1993, 1999); Kalas & Wynn-Williams Kal94 (1994); Böker et al. Bok98 (1998)). We note that the \[Ne ii\] linemaps of Böker et al. (Bok98 (1998)) and Keto et al. (Ket99 (1999)) likely contain a contribution from PAH 12.7$`\mu `$m emission because the observations were made at low spectral resolution ($`\mathrm{FWHM}0.2\mu \mathrm{m}`$).

The spatial distribution seen in all our maps is evidently dominated by a very prominent compact MIR source, better outlined in arcsecond resolution images. Keto et al. (Ket99 (1999)) associated this source with a bright off-nucleus super star cluster resolved by optical Hubble Space Telescope observations. Based on their data, the source has a size of $`20\mathrm{pc}`$ (marginally resolved) and accounts for $`20\%`$ of the total continuum emission at 12 and 20$`\mu `$m, and 12% of the total \[Ne ii\] 12.81$`\mu `$m flux. From our own measurements, nearly 25% of the total continuum emission detected with ISOCAM between 13.5 and 15$`\mu `$m originates in the 70 pc-diameter MIR peak; this fraction varies between 10% and 30% for our various broad- and narrow-band and emission feature measurements (tables 4 and 5). NGC 253 is quite remarkable in the compactness of its main MIR emitting region, with a $`\mathrm{FWHM}150\mathrm{pc}`$ while the optical disk extends over $`20\mathrm{kpc}`$. Although of a different nature, the interacting system NGC 4038/4039 offers a similar example, with 15% of its $`12.518\mu \mathrm{m}`$ luminosity being produced within a 100 pc-size star-forming knot in the overlapping region between the galaxies (Mirabel et al. Mir98 (1998)).

Additional small-scale structure of NGC 253 at MIR wavelengths includes a secondary much fainter peak in the 12 and 20$`\mu `$m continuum emission nearly coinciding with the nucleus, about 2″ northeast of the prominent source discussed above. The \[Ne ii\] 12.81$`\mu `$m line emission differs somewhat from the continuum, except possibly for the brightest peak, showing a bilobal or arc-like structure also hinted at in Br$`\gamma `$ images and suggestive of a circumnuclear star-forming ring (Böker et al. Bok98 (1998); Engelbracht et al. Eng98 (1998); Keto et al. Ket99 (1999)). None of these features, however, is resolved with ISOCAM.

The PAH 6.2$`\mu `$m/7.7$`\mu `$m and PAH 8.6$`\mu `$m/11.3$`\mu `$m ratio maps show variations of about 40% and nearly a factor of 2, respectively. These are statistically significant in view of the corresponding formal uncertainties of $`20\%`$ and $`<30\%`$. The ratio maps differ markedly, with the PAH 6.2$`\mu `$m/7.7$`\mu `$m image indicating a general increase from south to north of the nucleus whereas the PAH 8.6$`\mu `$m/11.3$`\mu `$m ratio appears more centrally concentrated, peaking near the nucleus and slightly more extended northeast. The reality of the arc-like feature $`10^{\prime \prime }`$ northwest of the nucleus in the PAH 8.6$`\mu `$m/11.3$`\mu `$m map is dubious because it is barely resolved and could be due to ghost effects given the strong unresolved peak of the emission. No corresponding structure that could perhaps support this feature is seen in maps of the molecular gas emission, of the radio continuum emission, of tracers of ionized gas or even of extinction (e.g. Peng et al. Pen96 (1996); Ulvestad & Antonucci Ulv97 (1997); Engelbracht et al. Eng98 (1998)). The limited region for which $`\mathrm{S}/\mathrm{N}3`$ prevents us from examining the larger-scale distribution in PAH ratios.

#### 3.2.3 NGC 1808

The starburst in NGC 1808 covers a region of comparable physical size to that of M 82 but being three times more distant, less structural details are resolved by ISOCAM. The PAH emission appears the most extended and oriented parallel to the major axis of the galaxy ($`\mathrm{P}.\mathrm{A}.140^{}`$), peaking at the nucleus. The 15$`\mu `$m emission follows closely the PAH emission. The \[Ar ii\] 6.99$`\mu `$m line emission is the most distinct in that the peak clearly is off-nucleus, about 5″ to the southeast, and the centroid of the emission region lies 2.5″ southeast of that for the 15$`\mu `$m continuum and PAH emission. The \[Ar ii\] 6.99$`\mu `$m map agrees well in peak position and extent with the global distribution of the most intense star-forming regions, or hot spots, as traced by H recombination lines and radio continuum emission (e.g. Saikia et al. Sai90 (1990); Krabbe et al. Kra94 (1994); Kotilainen et al. Kot96 (1996)). Due to the limited number of pixels with $`\mathrm{S}/\mathrm{N}3`$, no useful ratio maps could be made for NGC 1808.

Previously published $`N`$ band images of NGC 1808 (Telesco et al. Tel93 (1993); Krabbe et al. Kra01 (2001)) show an overall similar morphology as our 15$`\mu `$m continuum and PAH feature maps with, at $`1^{\prime \prime }`$ resolution, a strong point-like source at the nucleus and a southeastern extension covering the starburst regions. ISOCAM broad-band LW4 (6$`\mu `$m) polarisation observations were presented by Siebenmorgen et al. (Sie01 (2001)) along with a CVF spectrum of the central 25″ of NGC 1808 which is essentially identical in shape to ours of the starburst core (30″ aperture; Fig. 3). Siebenmorgen et al. (Sie01 (2001)) assigned all detected features to PAHs and successfully reproduced them with PAH emission alone in their radiative transfer models. As also noted by these authors, the contribution of nebular gas emission lines blended with PAH features is uncertain for NGC 1808, and we cannot constrain it using, e.g., SWS data as for M 82 and NGC 253. However, our \[Ar ii\] 6.99$`\mu `$m map clearly differs from those of the PAH emission and is consistent with the spatial distribution of the bulk of H ii regions, supporting the idea that the 7$`\mu `$m feature is indeed dominated by the \[Ar ii\] line instead of PAHs.

### 3.3 Salient features

We summarize here the most important aspects of the results presented above:

* in terms of features present, the MIR spectral energy distribution is essentially identical among the three galaxies as well as within each of them,
* the largest spectral variations are observed in the $`\lambda 11\mu \mathrm{m}`$ continuum intensity relative to the shorter wavelength emission, while the shape of the $`\lambda =511\mu \mathrm{m}`$ region is nearly invariant,
* the PAH emission is spatially more extended than the 15$`\mu `$m continuum and \[Ar ii\] 6.99$`\mu `$m line emission,
* the \[Ar ii\] 6.99$`\mu `$m line emission traces well the young starburst regions, and
* the relative intensities of the PAH features can exhibit complex spatial variations.

## 4 Origin and spatial distribution of the PAH and continuum emission

The ubiquity of PAH features – seen towards sources as diverse as circumstellar regions, diffuse atomic clouds, H ii regions, molecular clouds, normal galaxies, starbursts, and ULIRGs – and the close association of the MIR VSG continuum emission with Galactic and extragalactic actively star-forming sites now seem well established; the ISO mission played a major role in this recognition (see Geballe Geb97 (1997); Tokunaga Tok97 (1997); Tielens et al. Tie99 (1999); Cesarsky & Sauvage Ces99 (1999); Genzel & Cesarsky Gen00 (2000) for reviews). Spatial mapping of Galactic sources has revealed in more detail the origin of the main MIR emission components in star-forming regions: PAH features arise predominantly in the photodissociation regions (PDRs) at the interface between H ii regions and molecular clouds while the VSG $`\lambda 11\mu \mathrm{m}`$ continuum is more intense and steeper closer to the exciting source (e.g. Verstraete et al. Ver96 (1996); Cesarsky et al. 1996a ; 1996b ; Crété et al. Cre99 (1999)). By extension, in star-forming galaxies, the emission components are believed to trace PDRs and H ii regions, respectively. This is supported by spectral decomposition of starburst galaxies (Tran Tra98 (1998); Laurent et al. Lau00 (2000)), and by photometric and spectrophotometric imaging of spiral galaxies and starburst systems (e.g. Mirabel et al. Mir98 (1998); Roussel et al. 2001a ). Low excitation fine-structure line emission such as \[Ar ii\] 6.99$`\mu `$m, \[Ar iii\] 8.99$`\mu `$m, \[Ne ii\] 12.81$`\mu `$m, and \[Ne iii\] 15.56$`\mu `$m usually accompanies the long-wavelength VSG continuum and originates primarily in H ii regions (e.g. Sturm et al. Stu00 (2000)).

The interpretation of the ISOCAM maps of M 82, NGC 253, and NGC 1808 within the above framework is not straightforward. Figure 9 complements the various maps with light profiles of the $`69\mu \mathrm{m}`$, $`13.515\mu \mathrm{m}`$, and \[Ar ii\] 6.99$`\mu `$m emission taken along the major and minor axes of the galaxies. The profiles clearly support the extension along the major axis indicated by the maps, and show that the emission in these tracers along the minor axis is but marginally resolved and essentially indistinguishable within each galaxy.

Somewhat surprisingly given their common origin in H ii regions, a close spatial correlation between the VSG continuum and the \[Ar ii\] 6.99$`\mu `$m line emission is not observed in our sample galaxies. The exception may be NGC 253 but this probably results from the intrinsically small MIR source size and resolution limitations (Sect. 3.2.2). Although the extension to the northeast is real, the differences between the various profiles even along the major axis are smaller than the PSF FWHM. We conclude for NGC 253 that on scales comparable to and larger than the resolution of our ISOCAM maps (5.2″, or $`60\mathrm{pc}`$), no significant differences are observed in the spatial distribution of the PAH, VSG, and \[Ar ii\] 6.99$`\mu `$m emission within the central 400 pc.

For M 82, the \[Ar ii\] 6.99$`\mu `$m emission along the galactic plane peaks at a different position and is more compact than the VSG and PAH emission. The different distributions of the VSG and \[Ar ii\] 6.99$`\mu `$m emission could be interpreted in terms of differences in average energy of the exciting photons: the \[Ar ii\] 6.99$`\mu `$m line traces the $`16\mathrm{eV}`$ radiation field while the VSGs may be heated by UV photons at lower energies as well (Désert et al. Des90 (1990); Siebenmorgen & Krügel Sie92 (1992); Boulanger et al. Bou94 (1994); Dwek et al. Dwe97 (1997); Jones & d’Hendecourt Jon00 (2000)). The brighter PAH emission around the nucleus could be originating predominantly in the PDRs associated with the H ii regions within the most active starburst sites. Beyond a radius of $`300\mathrm{pc}`$, the PAH profile is resolved out from the VSG profile and more extended. These differences could be attributed to excitation of aromatic band carriers by yet softer photons or to an increased filling factor for the PDRs compared to the H ii regions in less intense starburst regions, and perhaps to a contribution from diffuse cirrus clouds as observed in the Milky Way and in some spiral galaxies (e.g. Ristorcelli et al. Ris94 (1994); Onaka et al. Ona96 (1996); Mattila et al. Mat99 (1999); Roussel et al. 2001a ; Li & Draine Li02 (2002)). Alternative possibilities could include destruction by photodissociation or depletion of the ISM in PAHs and VSGs in the most intense starburst regions, flattening the spatial distribution around the peak (e.g. Boulanger et al. Bou88 (1988); Carral et al. Car94 (1994); Normand et al. Nor95 (1995)).

In NGC 1808, the bulk of the \[Ar ii\] 6.99$`\mu `$m emission is clearly displaced to the southeast compared to the emission in the PAH and VSG bands defined above, which have virtually identical distributions. As discussed in Sect. 3.2.3, the \[Ar ii\] 6.99$`\mu `$m distribution agrees well with the off-nucleus location of the main star-forming sites. The spectrum integrated over the starburst core and MIR peak in NGC 1808 is also much flatter than observed towards the corresponding regions in M 82 and NGC 253 and resembles more the spectrum of their outer disks (Fig. 3). The long-wavelength continuum in NGC 1808 may in fact not be dominated by VSG emission from pure H ii regions but rather produced primarily in PDRs by small particles akin to the carriers of the main UIB bands. This is reminiscent of the situation in disks of spiral galaxies. There, the 15$`\mu `$m emission correlates with the 7$`\mu `$m PAH-dominated emission and lies in a distinct regime compared to more active regions including circumnuclear regions in spiral (especially barred) galaxies and starbursts, characterized by an excess in the 15$`\mu `$m/7$`\mu `$m colour (e.g. Roussel et al. 2001a ; Dale et al. Dal01 (2001); Förster Schreiber, Roussel, & Sauvage, in prep.). From radiative transfer modeling of NGC 1808’s $`3200\mu \mathrm{m}`$ spectrum, Siebenmorgen et al. (Sie01 (2001)) concluded that the MIR range is dominated by PAHs and large dust grains (radii between 100 and 2400 Å) with negligible contribution by very small graphites and silicates (radii of $`580`$ Å). The key feature of their models is the inclusion of hot spots where large grains are heated locally by massive stars to higher temperatures giving rise to the 25$`\mu `$m emission and gradually contributing less towards shorter wavelengths where PAHs take over as main emitters.

We note that the morphology of the \[Ar ii\] 6.99$`\mu `$m emission for M 82 and NGC 1808, especially the off-nucleus peak position, reflects the spatial distribution of the H ii regions and not variations in the excitation state of the ionized gas whereby, e.g., \[Ar iii\] emission would become important at the expense of \[Ar ii\]. As discussed in Sect. 3.2.1, the H ii regions in M 82 are concentrated in a circumnuclear ring-like structure at radius of $`85\mathrm{pc}`$. Direct tracers such as near-infrared to radio H recombination lines, MIR Ne, Ar, and S fine-structure lines, the He i 2.06$`\mu `$m recombination line, and radio thermal continuum all show a similar distribution characterized by peaks flanking the nucleus, with the southwestern complex being the most intense, and little emission at the nucleus (e.g. Satyapal et al. Sat95 (1995); Achtermann & Lacy Ach95 (1995); Seaquist et al. Sea96 (1996); Förster Schreiber et al. FS01 (2001)). In particular, the \[Ar iii\] 8.99$`\mu `$m map of Achtermann & Lacy (Ach95 (1995)) tracing the $`28\mathrm{eV}`$ radiation field clearly shows this bilobal distribution and relative lack of emission at the nucleus. Less information is available for NGC 1808 but Br$`\gamma `$ imaging clearly reveals the H ii region complexes to lie predominantly to the east and southeast of the nucleus, thus indicating a circumnuclear location for the starburst (Krabbe et al. Kra94 (1994); Kotilainen et al. Kot96 (1996)).

Extinction effects could partly account for differences in the PAH, VSG, and \[Ar ii\] 6.99$`\mu `$m emission on small spatial scales (comparable to the resolution elements in our ISOCAM data) but are not likely to affect the relative concentrations. The highest levels of obscuration are measured in the near surroundings of the MIR peaks in all three galaxies and generally decrease at larger distances (e.g. Larkin et al. Lar94 (1994); Engelbracht et al. Eng98 (1998); Krabbe et al. Kra94 (1994)). Depending on the extinction law applicable in the $`310\mu \mathrm{m}`$ region (Sect. 5.1 below), such a radial gradient in extinction would either affect very little the relative global distributions of the tracers discussed here, or suppress more severely the shorter wavelength emission in the central regions, which is inconsistent with the overall trends for the 15$`\mu `$m continuum and \[Ar ii\] 6.99$`\mu `$m line.

The nature of the short-wavelength continuum in normal and pure starburst galaxies is still debated. It is much weaker than in systems hosting an active galactic nucleus (AGN) where it is attributed to hot dust ($`1501700\mathrm{K}`$) close to the AGN (Genzel & Cesarsky Gen00 (2000) and references therein). In non-AGN galaxies, the existence of $`35\mu \mathrm{m}`$ featureless continuum from small transiently heated dust grains underlying the PAH features has been proposed (e.g. Helou et al. Hel00 (2000)). On the other hand, the PAH features are best fitted with Lorentzian profiles whose wide wings can account for the apparent continuum pedestal (e.g. Boulanger et al. 1998b ). Our narrow-band 5.5$`\mu `$m continuum map of M 82 provides only tentative evidence of a more extended distribution than for the 15$`\mu `$m continuum emission which would suggest a closer association of the emitting particles with PAHs than VSGs (see also Mattila et al. Mat99 (1999)).

## 5 PAH and continuum emission: relations to the ISM conditions and star-forming activity

The two most striking aspects seen in the spectra presented in Fig. 3 are their remarkable similarity in the $`\lambda =511\mu \mathrm{m}`$ regime dominated by PAH emission, and the spread in relative intensity of the longer wavelength emission attributed to VSGs. In this section, we focus on these two spectral components after discussing extinction effects at MIR wavelengths.

### 5.1 The effects of extinction at mid-infrared wavelengths

Interstellar extinction can significantly affect the MIR spectral energy distribution (SED) in obscured sources. In this regard, Rigopoulou et al. (Rig99 (1999)) noted the gradual suppression of the PAH 8.6$`\mu `$m feature in progressively obscured starbursts and ULIRGs. They further proposed that variations in the PAH 6.2$`\mu `$m/7.7$`\mu `$m ratio are dominated by extinction effects on the basis of the trend of decreasing ratio with increasing extinction observed in a subset of their sample for which reliable extinction determinations were available.

Another important consideration concerns the trough centered near 9.7$`\mu `$m traditionally attributed to absorption by interstellar silicate dust grains. Recent studies have cast doubt on the reliability of extinction estimates based on the observed feature depth. Sturm et al. (Stu00 (2000)) discussed this issue in detail based on SWS data of M 82 and NGC 253, from a comparison with the spectrum of the Galactic reflection nebula NGC 7023 and from considerations of the relative optical depths expected for silicate absorption near 9.7 and 18$`\mu `$m. In particular, they showed that the M 82 SWS spectrum, including the dip at 9.7$`\mu `$m, can be well reproduced by a combination of the NGC 7023 spectrum and a power law rising from 8.5$`\mu `$m longwards without invoking large extinctions. More generally, the overall spectral invariance in the $`511\mu \mathrm{m}`$ region among a variety of Galactic and extragalactic sources which are known or expected to cover a large range in extinction seems to argue in favour of the 9.7$`\mu `$m “absorption feature” being predominantly due to the gap between the main PAH emission complexes (e.g. Helou et al. Hel00 (2000)).

These suggestions raise the important question of how such an interpretation can be reconciled with the large and variable extinction in M 82, NGC 253, and NGC 1808 determined by alternative methods based, for instance, on the relative intensities of H recombination lines. To address this issue, we explored quantitatively the effects of extinction at MIR wavelengths by applying a range of extinction to a representative template spectrum constructed in a similar way as Sturm et al. (Stu00 (2000)). We added to a PDR component (from the ISOCAM data of NGC 7023) a power law $`f_\nu (\lambda 8.5)^\alpha `$ which represents well the VSG emission seen in our data. We adopted an index $`\alpha =1.5`$ and scaled the two components so as to obtain resulting spectra similar to that of M 82. For the purpose of illustrating the effects of different levels of extinction, the choice of these parameters is irrelevant.

We considered two representative geometries: a uniform foreground dust screen and a homogeneous mixture of dust and sources. For these models, the observed intensity $`I(\lambda )`$ of the emerging radiation is proportional to the total intrinsic intensity of the sources $`I_0(\lambda )`$ by the factors $`e^{\tau _\lambda }`$ and $`(1e^{\tau _\lambda })/\tau _\lambda `$, respectively, where $`\tau _\lambda `$ is the optical depth related to the extinction in magnitude via $`A_\lambda =1.086\tau _\lambda `$. We adopted the extinction law proposed by Draine (Dra89 (1989)) and investigated the effects of deviations from this law in the $`310\mu \mathrm{m}`$ region as found towards the Galactic Center (Lutz Lut99 (1999); hereafter simply “GC law”). Such deviations are consistent with the H recombination line spectrum observed with the SWS in M 82 (Förster Schreiber et al. FS01 (2001)). We varied the level of extinction specified in visual magnitudes $`A_V`$ for the mixed model in the range $`A_V^{\mathrm{MIX}}=0700\mathrm{mag}`$. For a meaningful comparison, we varied $`A_V`$ for uniform foreground obscuration in the range $`A_V^{\mathrm{UFS}}=050\mathrm{mag}`$, which results in the same attenuation at 10$`\mu `$m as the range considered for the mixed model. We normalized each extincted model spectrum to the integrated flux between 6.0 and 6.6$`\mu `$m, as for the spectra of Fig. 3.

Figure 10 presents our simulations. The plots outline the distinct behaviours for the two geometries considered with increasing $`A_V`$. They especially emphasize the fact that for the mixed model, differential extinction effects over the spectrum rapidly reach an asymptotic regime where large variations in $`A_V`$ become less and less perceptible as $`A_V`$ increases. In addition, except for the most extreme cases ($`A_V>20\mathrm{mag}`$), the choice of dust and source geometry has little impact. This is because the absolute levels of extinction are fairly small and the wavelength dependence of extinction is relatively weak; the wavelength range covered is thus too limited to probe appreciably different optical depths which could allow the discrimination between different dust and sources distributions.

For the Draine extinction law, the effects are largest between 8 and 13$`\mu `$m. Our models illustrate well how the PAH 8.6 and 11.3$`\mu `$m features are substantially suppressed as extinction increases, indicating that their ratio with the PAH features at 6.2 and 7.7$`\mu `$m can be much affected by obscuration. Moreover, for a continuum level defined over the limited $`813\mu \mathrm{m}`$ interval, relatively small optical depths would be inferred even if the extinction is in fact large. Assuming the GC law at $`310\mu \mathrm{m}`$ results in a remarkably different behaviour. The extinction effects are significantly smaller in the $`813\mu \mathrm{m}`$ region and the suppression of the PAH 8.6 and 11.3$`\mu `$m features is much less important. The largest effects are observed in the relative level between the short and long wavelength regions. This is a consequence of the much flatter GC law between 3 and 10$`\mu `$m.

Our simulations indicate that dust obscuration adds a significant degree of degeneracy in the interpretation of the observed MIR emission, complicating the determination of intrinsic properties. Depending on the extinction law used, the shape of the PAH 8.6 and 11.3$`\mu `$m features as well as their flux ratio with the PAH 6.2 and 7.7$`\mu `$m features vary considerably. Therefore, extinction estimates based on diagnostics involving the PAH 8.6 or 11.3$`\mu `$m feature may be strongly biased by the choice of extinction law. Furthermore, part of the spread in intensity observed in the long-wavelength continuum relative to the shorter wavelength emission could be attributed to varying extinction levels (see Fig. 3 and Sect. 5.2 below).

The interpretation of the 9.7$`\mu `$m trough as largely due to the gap between the flanking PAH complexes appears to hold over a wide range of extinction, especially for the GC extinction law. The near-invariance of the $`511\mu \mathrm{m}`$ spectrum is therefore not inconsistent with large and/or variable obscuration among and within galaxies. We illustrate this with the case of M 82, for which a mixed model with $`A_V=52\mathrm{mag}`$ best reproduces the observed H recombination lines from radio to optical wavelengths while the best-fit uniform foreground screen extinction of $`A_V=5\mathrm{mag}`$ provides a much poorer fit (Förster Schreiber et al. FS01 (2001)). We applied these two extinction models to the same template as for the simulations described above. The resulting spectra are plotted along with the observed spectrum of M 82 in the bottom panels of Fig. 10. We did not attempt to fine-tune the models by formal fitting. Given the uncertainties on the exact nature of the emitting particles, specific assumptions on model parameters are not well constrained and a simple empirical approach is sufficient for our purposes. The comparison shows that strictly from the point of view of the $`516\mu \mathrm{m}`$ range, both extinction models reproduce equally well the observed SED of M 82, thus demonstrating that a high extinction cannot be excluded from the overall shape of the MIR spectrum alone.

We wish to emphasize that we do not dismiss the possibility of silicate absorption around 9.7$`\mu `$m in general but merely want to point out the difficulties involved in the interpretation of the observed feature. The SEDs of moderately to highly obscured sources including ULIRGs often exhibit a strong dip near 9.7$`\mu `$m together with suppressed PAH 8.6 and 11.3$`\mu `$m features relative to those at 6.2 and 7.7$`\mu `$m, consistent with the presence of silicate grains (e.g. Dudley & Wynn-Williams Dud97 (1997); Dudley Dud99 (1999); Laurent et al. Lau00 (2000); Le Floc’h et al. LeF02 (2002)). Given the near invariance of the PAH spectrum in a wide range of environments, it should be possible to define an indicator measuring the differential extinction between the 9.7$`\mu `$m region and adjacent less affected intervals that quantifies the absolute extinction. One must however be aware of the importance of sufficient wavelength coverage and resolution to assess properly the impact of the PAH complexes which can make the silicate absorption look artificially deep, of the dependence of the inferred $`A_V`$ on the assumed extinction law and spectral intervals used to measure the feature depth, and of the limited sensitivity of this diagnostic leading to rather large uncertainties in the derived $`A_V`$. The upcoming launch of SIRTF in 2003 will provide $`540\mu \mathrm{m}`$ spectroscopy with an increase in sensitivity by two orders of magnitude and will help us address this issue by enabling a better sampling of the continuum emission as well as both the 9.7 and 18$`\mu `$m silicate bands.

### 5.2 The $`\lambda \mathit{11}\mu m`$ continuum emission

The $`\lambda 11\mu \mathrm{m}`$ continuum emission in our data exhibits a large spread in intensity relative to the shorter wavelength emission. Similar variations have been seen in normal spirals and starburst galaxies observed with ISO instruments and, based on evidence provided by Galactic sources, are generally interpreted in terms of the relative contribution of H ii regions to the MIR emission (e.g. Laurent et al. Lau00 (2000); Roussel et al. 2001a ; Dale et al. Dal01 (2001)). A close link between the 15$`\mu `$m/7$`\mu `$m colour as measured through the ISOCAM broad-band LW3 and LW2 filters ($`\lambda =1215\mu \mathrm{m}`$ and $`58.5\mu \mathrm{m}`$, respectively) and the IRAS 25$`\mu `$m/12$`\mu `$m colour has been emphasized by Dale et al. (Dal01 (2001)) and Roussel et al. (2001a ).

In order to assess quantitatively the relationship between the long-wavelength continuum properties and the H ii regions, we investigated the relationship between the 15$`\mu `$m narrow-band continuum and the \[Ar ii\] 6.99$`\mu `$m line emission. The \[Ar ii\] 6.99$`\mu `$m is the most direct probe of H ii regions and the least contaminated by PAH emission available from our ISOCAM data sets. The $`{}_{}{}^{}\mathrm{Ar}_{}^{0}`$ ionization potential of 15.8 eV is close to that of H (13.6 eV). Argon being a noble element, it is, like neon, not expected to be significantly depleted onto dust grains; while the argon gas-phase abundance may increase with time as a result of star formation activity and differ among galaxies, it is not likely to vary significantly over scales of $`1001000\mathrm{pc}`$ as covered by the starburst regions in our sample. For the relatively low nebular excitation in all three sources (Fig. 4; see also Sturm et al. Stu00 (2000); Förster Schreiber et al. FS01 (2001)), argon is mostly singly-ionized so that the \[Ar ii\] 6.99$`\mu `$m emission should trace the H ii regions in much the same way as H recombination lines and its luminosity, scale with the radiation field intensity of the young stars. In addition, since the gas-phase abundances of the galaxies are solar within a factor of a few (e.g. Webster & Smith Web83 (1983); Forbes et al. For92 (1992); Förster Schreiber et al. FS01 (2001)), the proportionality factor should be roughly similar.

Figure 11 presents the 15$`\mu `$m continuum versus \[Ar ii\] 6.99$`\mu `$m luminosities for our galaxies, normalized to unit projected surface area ($`\mathrm{\Sigma }_{15\mu \mathrm{m}}`$ and $`\mathrm{\Sigma }_{[\mathrm{Ar}\mathrm{II}]}`$). The data are shown for the selected regions as well as for the individual resolution elements enclosed within the radii defining the outer limit of the disk regions in M 82 and NGC 253. Resolution elements with measurements at $`<3\sigma `$ were excluded. The data points follow a well defined distribution, with a remarkable overlap for the different galaxies. Least-squares fitting to the resolution elements’ data accounting for the individual formal uncertainties yields the relation

$$\mathrm{log}(\mathrm{\Sigma }_{15\mu \mathrm{m}})=1.01\mathrm{log}(\mathrm{\Sigma }_{[\mathrm{Ar}\mathrm{II}]})+0.807,$$
(1)
with dispersion around the fit of 0.203 dex. The extinction paths in the diagram of Fig. 11 are nearly parallel to the locus formed by the data points; however, implausibly large differential extinctions would be required in order to account entirely for the observed distribution.

The essentially linear proportionality we find between $`\mathrm{\Sigma }_{15\mu \mathrm{m}}`$ and $`\mathrm{\Sigma }_{[\mathrm{Ar}\mathrm{II}]}`$ indicates that the 15$`\mu `$m emission provides a good quantitative indicator of the star-forming activity in starburst environments, to within uncertainties determined by the dispersion of the data and by extinction (a factor of 2.5 for an $`A_V=50\mathrm{mag}`$ assuming purely foreground extinction). Roussel et al. (2001b ) reached a similar conclusion for more quiescent spiral disks where the $`1215\mu \mathrm{m}`$ emission scales linearly with the H$`\alpha `$ line flux. Part of the scatter in Fig. 11 may result from variations in the physical conditions and exact composition of the gas and dust within and between the galaxies. More important factors, however, are variations in relative spatial distribution of the emission and possibly of the excitation state of the gas.

The observed relationship holds remarkably well in view of the differences in morphology of the 15$`\mu `$m continuum and \[Ar ii\] 6.99$`\mu `$m line emission as seen in Figs. 5, 6, 7, and 9. The correlation extends however over two to three orders of magnitude, much larger than the dispersion of a factor of 1.6. Undoubtedly, the relative spatial variations between both tracers contribute significantly to the scatter. Within the starburst cores of M 82 and NGC 253, where the steeply rising SEDs at $`\lambda 11\mu \mathrm{m}`$ are dominated by VSG emission, the differences may be attributed to different ranges in exciting photon energies for VSGs and \[Ar ii\] 6.99$`\mu `$m line as well as to extinction effects (Sect. 4). For NGC 1808, the spatial distributions differ probably because a different dust/particles population produces the flat long-wavelength continuum (Sect. 4). Nevertheless, the data for NGC 1808 and for the disk regions of M 82 and NGC 253, characterized by flat $`11\mu \mathrm{m}`$ SEDs, are well described by Eq. 1.

Variations of the excitation state of the photoionized nebulae may influence the $`\mathrm{\Sigma }_{15\mu \mathrm{m}}`$ versus $`\mathrm{\Sigma }_{[\mathrm{Ar}\mathrm{II}]}`$ relationship as well. In particular, the \[Ar iii\] 8.99$`\mu `$m line is fairly strong in the ISO-SWS spectrum of M 82 (Fig. 4) but the \[Ar iii\] 8.99$`\mu `$m/\[Ar ii\] 6.99$`\mu `$m ratio of 0.18 (0.26 after extinction correction) is low and abundance estimates indicate that $`25\%`$ only of the argon is doubly ionized (Förster Schreiber et al. FS01 (2001)). As emphasized in Sect. 4, the spatial distributions of various MIR fine-structure lines and Br$`\alpha `$ are similar and suggest a roughly constant excitation state of the H ii regions across M 82 (Achtermann & Lacy Ach95 (1995)). This is further confirmed by the nearly uniform He i 2.06$`\mu `$m/Br$`\gamma `$ ratio in the central $`260\times 160\mathrm{pc}`$ which corresponds closely to the excitation derived from the SWS \[Ne iii\] 15.56$`\mu `$m/\[Ne ii\] 12.81$`\mu `$m, \[Ar iii\] 8.99$`\mu `$m/\[Ar ii\] 6.99$`\mu `$m, and \[S iv\] 10.5$`\mu `$m/\[S iii\] 18.7$`\mu `$m ratios within larger apertures up to $`430\times 225\mathrm{pc}`$ (Förster Schreiber et al. FS01 (2001)). Similar arguments are more difficult for NGC 253 and NGC 1808 because of the lack of relevant data. However, the \[Ar iii\] 8.99$`\mu `$m/\[Ar ii\] 6.99$`\mu `$m $`0.03`$ measured from the SWS spectrum of NGC 253 is substantially lower than for M 82, as is the \[Ne iii\] 15.56$`\mu `$m/\[Ne ii\] 12.81$`\mu `$m ratio (0.06 compared to 0.18; Thornley et al. Tho00 (2000); see also Giveon et al. Giv02 (2002) for the relationship between these two MIR line ratios). For NGC 1808, we can only note from the ISOCAM spectra (Fig. 3) that the weakness of the 15.7$`\mu `$m feature and the non-detection of \[Ar iii\] 8.99$`\mu `$m indicate low excitation supporting that argon mostly is singly ionized. Summarizing, the excitation state of the gas might cause for M 82 a general but small offset to the left in the relationship of Fig. 11 compared to the other galaxies because of the larger fraction of $`\mathrm{Ar}^{++}`$ while spatial variations are not likely to introduce scatter larger than this offset. These effects are expected to be smaller for both NGC 253 and NGC 1808.

The relationship between the 15$`\mu `$m continuum and star formation intensity is supported by the analysis of a larger sample including spiral and starburst galaxies that will be presented in a forthcoming paper (Förster Schreiber, Roussel, & Sauvage, in prep.). In this paper, based on the previous work of Roussel et al. (2001b ) for spiral disks, we find that the 15$`\mu `$m continuum as well as the 7$`\mu `$m PAH-dominated emission correlate well with the star formation rate over nearly 7 orders of magnitude, with the interesting distinction of two regimes that correspond to quiescent star formation in disks and more intense activity in circumnuclear regions and starbursts.

### 5.3 The $`\lambda =\mathit{5}\mathit{11}\mu m`$ emission

The near invariance of the $`511\mu \mathrm{m}`$ spectrum seen in our data of M 82, NGC 253, and NGC 1808 has been noted in a number of other studies of different types of galaxies powered mainly by star formation, as well as in a variety of Galactic sources (see Tielens et al. Tie99 (1999) for a review; see also e.g. Boulanger et al. 1998a ; Helou et al. Hel00 (2000); Uchida et al. Uch00 (2000)). This indicates that the PAHs are very stable under a wide range of physical conditions despite their small sizes (typically $`100`$ atoms) and is taken as direct observational evidence for the stochastic nature of the emission processes involved. Furthermore, if different types of PAHs coexist, the near constancy in the $`511\mu \mathrm{m}`$ region suggests that their relative abundances vary little. Finally, it has consequences on the interpretation of and extinction measurements from the 9.7$`\mu `$m silicate absorption feature because of the intrinsic gap between the main $`69\mu \mathrm{m}`$ and $`1113\mu \mathrm{m}`$ PAH complexes (Sect. 5.1).

Variations of the relative intensities of the PAH features do exist, however. We focus on the PAH 6.2$`\mu `$m/7.7$`\mu `$m and PAH 8.6$`\mu `$m/11.3$`\mu `$m ratios, maps of which were presented for M 82 and NGC 253 in Sect. 3.2. Of particular interest is also the enhancement of the PAH 11.3$`\mu `$m relative to the other features in NGC 1808 compared to M 82 and NGC 253, evident in Fig. 3. The variations of PAH ratios in our sample are comparable to those observed within and between Galactic sources of similar types (e.g. Cesarsky et al. 1996a ; Boulanger et al. 1998a ; Lu Lu98 (1998); Crété et al. Cre99 (1999); Uchida et al. Uch00 (2000)), along the disk of the spiral galaxy NGC 891 (Mattila et al. Mat99 (1999)), and among a sample of 15 starbursts and ULIRGs (Rigopoulou et al. Rig99 (1999)). A detailed discussion of PAH ratios is beyond the scope of this paper, but we briefly mention possible interpretations of the variations seen in our data in the light of some recent theoretical and empirical work.

Since the emission mechanism is stochastic in nature, PAH ratios are not directly related to the SED of the incident radiation (e.g. Boulanger et al. 1998a ; Uchida et al. Uch00 (2000)). On the other hand, they may carry indirect information since PAHs exposed to hard and intense radiation fields can be ionized, lose H atoms, or be photodissociated (e.g. Léger et al. Leg89 (1989); Allamandola et al. All89 (1989); Schutte et al. Sch93 (1993); Allamandola et al. All99 (1999); Hudgins & Allamandola Hud99 (1999)). In particular, the PAH 8.6$`\mu `$m/11.3$`\mu `$m is believed to trace the fraction of singly ionized to neutral PAHs, presumably driven by the strength of the radiation field from the OB stars (see references above, and Joblin et al. Job96 (1996)). Draine & Li (Dra01 (2001)) presented a thorough calculation of the expected PAH spectrum as a function of various parameters including PAH size, charging conditions, and starlight intensity. Their results indicate that the PAH 6.2$`\mu `$m/7.7$`\mu `$m depends primarily on PAH size while the PAH 11.3$`\mu `$m/7.7$`\mu `$m is mainly sensitive to the fraction of ionized versus neutral PAHs and only modestly to PAH size. The effects of radiation field become noticeable only at high intensities ($`10^5`$ the average local Galactic far-UV flux in Habing units) and for large PAHs with $`10^2`$ carbon atoms. Mattila et al. (Mat99 (1999)) suggested that PAH 6.2$`\mu `$m/7.7$`\mu `$m variations may be due to differences in average temperature of the PAHs during their temperature spikes, related to PAH size or mean exciting UV photon energies. Alternatively, they could be attributed to broadening of the 7.7$`\mu `$m feature at low radiation field energy densities (Uchida et al. Uch00 (2000)); this would reduce the PAH 7.7$`\mu `$m flux and increase the continuum level within our fixed bandpasses. In addition to intrinsic variations, extinction can significantly alter the shape and relative intensities of the PAH features as shown in Sect. 5.1.

Of our sample galaxies, M 82 offers the most interesting case for interpretation. The overall bilobal structure in the maps of PAH ratios and CO gas distribution (Fig. 8) might result from the different composition and physical processes that the emitting PAHs undergo when exposed to the varying conditions across M 82, from the starburst core to the more quiescent disk via the transition regions marked by the molecular gas ring. The maxima in PAH 8.6$`\mu `$m/11.3$`\mu `$m ratio lie at smaller radii than the minima in PAH 6.2$`\mu `$m/7.7$`\mu `$m ratio, possibly indicating a higher degree of PAH ionization within the most intense starburst sites at the inner edge of the molecular ring; this is particularly striking southwest of the nucleus where the peak PAH 8.6$`\mu `$m/11.3$`\mu `$m coincides very well with the location of the most prominent H ii region complexes (Sect. 3.2.1). The variations of PAH 6.2$`\mu `$m/7.7$`\mu `$m ratio could reflect differences in the PAH size distribution, combined with extinction effects, where larger PAHs can better form and survive in denser, more shielded environments associated with molecular gas concentrations. In NGC 253, the peak PAH 8.6$`\mu `$m/11.3$`\mu `$m ratio at the nucleus where the starburst is mainly occurring could also be due to a larger fraction of ionized PAHs. Finally, the enhanced PAH 11.3$`\mu `$m feature in NGC 1808 compared to M 82 and NGC 253 could be explained by a more neutral mixture of PAHs and a more diffuse radiation field.

The relative intensity of PAH and continuum emission can also hold information on the physical environment within astronomical sources. In particular, the ratio of the peak intensity of the PAH 7.7$`\mu `$m to the underlying continuum, hereafter PAH 7.7$`\mu `$m L/C ratio, has been shown to constitute a powerful discriminator between star formation activity or an AGN as the main source of the bulk of infrared luminosity (Rigopoulou et al. Rig99 (1999); see also Genzel et al. Gen98 (1998); Laurent et al. Lau00 (2000); Tran et al. Tra01 (2001)). Starburst-dominated objects are characterized by ratios $`1`$ while AGN-dominated ones have ratios $`1`$. We compared measurements of the PAH 7.7$`\mu `$m L/C ratio from our data of M 82, NGC 253, and NGC 1808 with those obtained by Rigopoulou et al. (Rig99 (1999)) for a large sample including starbursts, AGNs, and ULIRGs. Our motivation was to assess how much this indicator depends on the source luminosity at the faint end and compare the spread in global ratios among different galaxies with the spatial variations within individual galaxies.

Figure 12 shows the PAH 7.7$`\mu `$m L/C ratio versus peak PAH 7.7$`\mu `$m luminosity for selected regions of our galaxies, for individual resolution elements of M 82 and NGC 253 (within the radii defining the outer limit of the disk regions in table 2), and for the Rigopoulou et al. (Rig99 (1999)) sample. We measured the PAH 7.7$`\mu `$m peak intensity and underlying continuum according to the definitions of Rigopoulou et al. We note that with these definitions, the continuum might be underestimated in highly obscured sources as a result of extinction affecting notably the 11$`\mu `$m region (as also discussed by Laurent et al. Lau00 (2000)). We chose the PAH 7.7$`\mu `$m peak luminosity as other characteristic property, taken as an approximate indicator of the total infrared luminosity. Although the PAH to infrared luminosity ratio can vary by up to factors of several in different environments (Rigopoulou et al. Rig99 (1999)), our assumption has little consequences on the interpretation of Fig. 12 since the data span a range in luminosity extending over more than six orders of magnitude444Additional though indirect support for our assumption is provided by the relationship between the 7$`\mu `$m flux measured through the ISOCAM LW2 filter ($`58.5\mu \mathrm{m}`$) and the far-infrared (FIR, $`43123\mu \mathrm{m}`$) flux obtained by Roussel et al. (2001b ) for galactic disks, with $`F_{7\mu \mathrm{m}}(F_{\mathrm{FIR}})^{1.2}`$. Furthermore, Dale et al. (Dal01 (2001)) derived a typical value of $`15\%`$ for the $`520\mu \mathrm{m}`$ to FIR luminosity ratio in normal galaxies while Charmandaris et al. (Cha02 (2002)) found this ratio to be $`5\%`$ for three ULIRGs..

Our data of M 82, NGC 253, and NGC 1808 extend very well the trend defined by the global properties of pure starburst galaxies and ULIRGs, populating the PAH 7.7$`\mu `$m L/C $`1`$ region down to luminosities about an order of magnitude lower. The ratios in all three galaxies lie well above the starburst–AGN separation at a ratio of unity and form a tighter distribution. The average and $`1\sigma `$ dispersion for the individual resolution elements in M 82 and NGC 253 are $`5.0\pm 0.9`$ compared to $`3.7\pm 2.5`$ for the pure starbursts of Rigopoulou et al. and $`2.8\pm 1.0`$ for their ULIRGs (excluding those with ratios $`1`$ or with only limits on the measurements). It is also interesting to note that the starburst trend holds for regions on spatial scales ranging from $`60\mathrm{pc}`$ for the smallest individual regions in NGC 253 up to several kiloparsecs for the largest starbursts and ULIRGs (see the near-infrared images of Rigopoulou et al. Rig99 (1999)).

## 6 Summary and conclusions

We have presented ISOCAM $`\lambda =5.016.5\mu \mathrm{m}`$ spectrophotometric imaging of the starburst galaxies M 82, NGC 253, and NGC 1808. The spectrum of all three objects, down to the smallest scales of $`100\mathrm{pc}`$ accessible from our data, exhibit similar characteristics including prominent PAH bands, a featureless continuum most obvious at $`\lambda 11\mu \mathrm{m}`$, and a trough in the 10$`\mu `$m region. We securely identified the main emission features detected in the $`R40`$ ISOCAM data of M 82 and NGC 253 based on their $`R5001000`$ SWS spectra. The comparison emphasizes the caution that should be exercised when interpreting low resolution data because of potential blends or misidentifications between emission lines from ionized gas and PAH features originating in PDRs, a notable example being \[Ne ii\] 12.81$`\mu `$m and PAH 12.7$`\mu `$m. Using a simple model combining a template PDR spectrum and a power-law $`f_\nu (\lambda 8.5)^{1.5}`$, we constructed a representative starburst SED and explored the effects of extinction at MIR wavelengths. Our simulations illustrate the importance of the assumed extinction law (e.g. the widely used Draine Dra89 (1989) law versus the Galactic Center law of Lutz Lut99 (1999)) and of the intrinsic PAH spectrum (especially the gap between the main $`69\mu \mathrm{m}`$ and $`1113\mu \mathrm{m}`$ complexes) in shaping the SED of astronomical sources. This complicates the interpretation of PAH ratios as well as extinction measurements relying on the silicate dust absorption at 9.7$`\mu `$m (see also Sturm et al. Stu00 (2000)).

As observed previously in a wide range of Galactic and extragalactic sources, the $`511\mu \mathrm{m}`$ spectrum in our galaxies is nearly invariant. The relative PAH intensities exhibit nevertheless measureable and significant variations of $`20\%100\%`$ which may be attributed to various, possibly interrelated effects including the intensity of the incident radiation field and the PAH size distribution, ionization, and dehydrogenation. In our sample, M 82 probably best illustrates variations of PAH ratios due to an increased fraction of ionized PAHs within the most intense starburst sites and, admittedly speculatively, perhaps also to differences in typical PAH sizes depending on the molecular gas concentrations. The PAH 7.7$`\mu `$m L/C ratio in all three galaxies clearly lies in the range observed for pure starburst systems and extends the trend reported previously by Rigopoulou et al. (Rig99 (1999)) to lower luminosities.

In contrast, the $`\lambda 11\mu \mathrm{m}`$ region varies most among our sample galaxies and the ISOCAM maps show a comparatively more compact 15$`\mu `$m continuum distribution relative to the PAH emission. Strong 15$`\mu `$m continuum in M 82 and NGC 253 indicates an important contribution by VSGs contrary to the case of NGC 1808 where the long-wavelength emission is much flatter presumably because of negligible contribution by VSGs. We find, however, that in all three galaxies the 15$`\mu `$m continuum and \[Ar ii\] 6.99$`\mu `$m line fluxes satisfy a linear relationship. We infer from this that the 15$`\mu `$m continuum provides a good indicator of star formation activity in starbursts, complementing the similar results of Roussel et al. (2001b ) for galactic disks. In a broader perspective, our galaxies fit well in the trend of increasing ISOCAM 15$`\mu `$m/7$`\mu `$m ratios with higher levels of star formation activity found among normal disk galaxies and starburst-powered LIRGs/ULIRGs (e.g. Laurent et al. Lau00 (2000); Roussel et al. 2001b ; Dale et al. Dal01 (2001)).

The value of the ISOCAM spectrophotometric imaging of M 82, NGC 253, and NGC 1808 presented in this paper also lies in that it complements existing MIR data with maps of PAH features, fine-structure lines, and continuum components not previously imaged. The poorer angular resolution of ISOCAM compared to that achieved with large ground-based telescopes is compensated by the larger field of view revealing more of the large scale emission. In that respect, we stress the small size of the MIR source relative to the optical extent of all three galaxies: the ISOCAM maps cover the central $`1.5\mathrm{kpc}`$ for M 82, 0.6 kpc for NGC 253, and 5 kpc for NGC 1808 while the optical diameters are about 10 kpc for M 82 and 20 kpc for the other two. By measuring the flux density within the ISOCAM LW10 filter bandpass, equivalent to the IRAS 12$`\mu `$m band, we recover all of the IRAS 12$`\mu `$m emission in the entire field of view for M 82 and NGC 1808, and about 70% for NGC 253. This suggests that the total MIR emission is strongly dominated by the starburst sites in the nuclear regions while the more quiescent star formation taking place in the disk at larger radii does not contribute much. It will be interesting to see whether the MIR cameras on board SIRTF, to be launched in 2003, will confirm this result or discover faint diffuse emission, especially at shorter ($`35\mu \mathrm{m}`$) wavelengths.

###### Acknowledgements.

We are grateful to E. Sturm for making the SWS spectrum of NGC 253 available to us in electronic form. We warmly thank H. Roussel, A. Vogler, S. Madden, and especially D. Tran for many interesting discussions on various aspects of this work. We also wish to thank the referee, Dr. N. Bergvall, for his valuable comments that improved the quality of the paper. VC would like to acknowledge the partial support of JPL contract 960803.


text:
# The End of Unified Dark Matter?

## I Introduction

Despite the broad interest in dark matter and dark energy, their physical properties are still poorly understood. Indeed, it has never even been shown that the two are in fact two separate substances. The goal of this paper is to show that they are.

There is strong evidence from a multitude of observations that there is about six times more cold dark matter (CDM) than baryons in the cosmic matter budget, making up of order 30% of critical density . In addition to this clustering dark component, observations of supernovae, the cosmic microwave background fluctuations and galaxy clustering provide mounting evidence of a uniformly distributed dark energy with negative pressure which has come to dominate the universe recently (at redshifts $`z1`$) and caused its expansion to accelerate. It currently constitutes about two thirds of the critical density .

Although the dark energy can be explained by introducing the cosmological constant ($`\mathrm{\Lambda }`$) into general relativity (lending the standard model the name $`\mathrm{\Lambda }`$CDM), this “solution” has two severe problems, frequently triggering anthropic explanations and general unhappiness. The first problem is explaining its magnitude, since theoretical predictions for $`\mathrm{\Lambda }`$ lie many orders of magnitude above the observed value. The second problem is the so-called cosmic coincidence problem: explaining why the three components of the universe (matter, radiation and $`\mathrm{\Lambda }`$) presently are of similar magnitudes although they all scale differently with the Universe’s expansion.

As a response to these problems, much interest has been devoted to models with dynamical vacuum energy, so-called quintessence. These models typically involve scalar fields with a particular class of potentials, allowing the vacuum energy to become dominant only recently. Although quintessence is the most studied candidate for the dark energy, it generally does not avoid fine tuning in explaining the cosmic coincidence problem. Recently several alternative models have been proposed, such as the condensate models of.

An alternative to quintessence which has attracted great interest lately is the so-called generalized Chaplygin gas (hereafter GCG)  (see also the related earlier work of). Rather than fine tuning some potential, the model explains the acceleration of the Universe via an exotic equation of state causing it to act like dark matter at high density and like dark energy at low density. The model is interesting for phenomenological reasons but can be motivated by a brane-world interpretation. An attractive feature of the model is that it can explain both dark energy and dark matter in terms of a single component, and has therefore been referred to as unified dark matter (UDM) or “quartessence” . (See also .)

This approach has been thoroughly investigated for its impact on the 0th order cosmology, i.e., the cosmic expansion history (quantified by the Hubble parameter $`H[z])`$ and corresponding spacetime-geometric observables. An interesting range of models was found to be consistent with SN Ia data  and CMB peak locations .

Some work has also studied constraints from 1st order cosmology (the growth of linear perturbations), finding an interesting range of models to be consistent with cosmic microwave background (CMB) measurements . There is, however, a fatal flaw in UDM models that manifests itself only at recent times and on smaller (Galactic) scales and has therefore not been revealed by these studies. As we will see, this flaw rules out all GCG models except those that are for all practical purposes identical to the usual $`\mathrm{\Lambda }`$CDM model.

The rest of this Letter is organized as follows. In the next section, we review the fundamentals of the GCG model. We then consider in section III the evolution of density inhomogeneities in the model and use the predicted matter power spectrum to constrain it with observational data. Finally, we describe how the basic flaw that rules out the GCG model is indeed a generic feature of a broad class of unified dark matter models.

## II The Chaplygin gas

A standard assumption in cosmology is that the pressure of a single substance is, at least in linear perturbation theory, uniquely determined by its density. A generalized Chaplygin gas  is simply a substance where this relation $`p(\rho )`$ is a power law

$$p=A\rho ^\alpha $$
(1)
with $`A`$ a positive constant. The original Chaplygin gas had $`\alpha =1`$. The standard $`\mathrm{\Lambda }`$CDM model has two separate dark components, both with $`\alpha =1`$, giving a constant equation of state $`wp/\rho `$ that equals $`0`$ for dark matter and $`1`$ for dark energy.

By inserting equation (1) into the energy conservation law, one finds that the GCG density evolves as

$$\rho (t)=\left[A+\frac{B}{a^{3(1+\alpha )}}\right]^{\frac{1}{1+\alpha }},$$
(2)
where $`a(t)`$ is the cosmic scale factor normalized to unity today, i.e., $`a=(1+z)^1`$ where $`z`$ denotes redshift. Here $`B`$ is an integration constant. The striking feature here is that although the GCG has $`\rho a^3`$ when sufficiently compressed, it’s density will never drop below the value $`A^{1/1+\alpha }`$ no matter how much you expand it. Defining

$$\mathrm{\Omega }_m^{}\frac{B}{A+B},\rho _{}(A+B)^{\frac{1}{1+\alpha }},$$
(3)
equation (2) takes the form

$$\rho (a)=\rho _{}\left[(1\mathrm{\Omega }_m^{})+\mathrm{\Omega }_m^{}a^{3(1+\alpha )}\right]^{\frac{1}{1+\alpha }}.$$
(4)
For comparison, a standard flat model with current CDM density parameter $`\mathrm{\Omega }_m`$ as well as dark energy density $`(1\mathrm{\Omega }_m)`$ whose equation of state $`w_{}`$ is constant gives

$$\rho (a)=\rho _{}\left[(1\mathrm{\Omega }_m)a^{3(1+w_{})}+\mathrm{\Omega }_ma^3\right].$$
(5)
We see that the last two equations bear a striking similarity even though the former involves a single substance and the latter involves two. Both have two free parameters. Both have the current density $`\rho (1)=\rho _{}`$. Making the identification $`\mathrm{\Omega }_m^{}=\mathrm{\Omega }_m`$, both have $`\rho (a)\mathrm{\Omega }_m\rho _{}a^3`$ at early times as $`a0`$ (for $`w_{}<0`$), showing that $`\mathrm{\Omega }_m^{}`$ can be interpreted as an effective matter density in the GCG model. Indeed, for the special case $`\alpha =0`$ and $`w_{}=1`$, we see that both models coincide with standard $`\mathrm{\Lambda }`$CDM. For $`\alpha =0`$ the GCG model becomes equivalent to $`\mathrm{\Lambda }`$CDM not only to 0th order in perturbation theory as above but to all orders, even in the nonlinear clustering regime.

The 0th order cosmology determined by equation (4) together with the Friedman equation

$$H\frac{\dot{a}}{a}=\left[\frac{8\pi G}{3}\rho \right]^{1/2}$$
(6)
(which determines $`a(t)`$ and the spacetime metric to 0th order) has been thoroughly in previous work, and by studying constraints from supernovae observations, Makler et al.  have placed interesting constraints on the $`(\alpha ,\mathrm{\Omega }_m^{})`$ parameter space.

## III Growth of inhomogeneities

Let us now consider the evolution of density perturbations in this UDM model. Following the standard calculations of , we obtain for the relativistic analog of the Newtonian 1st order perturbation equation in Fourier space that a density fluctuation $`\delta _k`$ with wave vector $`𝐤`$ evolves as

$`\ddot{\delta }_k+H\dot{\delta }_k[23(2wc_s^2)]`$ (7)
$`{\displaystyle \frac{3}{2}}H^2\delta _k\left[16c_s^23w^2+8w\right]=\left({\displaystyle \frac{kc_s}{a}}\right)^2\delta _k,`$ (8)
where the equation of state $`wp/\rho `$ and the squared sound speed $`c_s^2p/\rho `$ are evaluated to 0th order and hence depend only on time, not on position. (We use units where the speed of light $`c=1`$ throughout.) This equation is valid on subhorizon scales $`|𝐤|H/c`$. In other words, the growth of density fluctuations is completely determined by the two functions $`w(a)`$ and $`c_s^2(a)`$. Combining equation (1) and equation (4), these two functions are

$`w`$ $`=`$ $`\left[1+{\displaystyle \frac{\mathrm{\Omega }_m^{}}{1\mathrm{\Omega }_m^{}}}a^{3(1+\alpha )}\right]^1,`$ (9)
$`c_s^2`$ $`=`$ $`\alpha w=\alpha \left[1+{\displaystyle \frac{\mathrm{\Omega }_m^{}}{1\mathrm{\Omega }_m^{}}}a^{3(1+\alpha )}\right]^1.`$ (10)
This shows a second reason why the GCG has been considered promising for cosmology it starts out behaving like pressureless CDM ($`w0`$, $`c_s0`$) early on (for $`a1`$) and gradually approaches cosmological constant behavior ($`w1`$) at late times. There is also an intermediate state where the effective equation of state is $`p=\alpha \rho `$. (Going beyond 1st order perturbation theory, the GCG that gets gravitationally bound in galactic halos maintains its density high enough to keep acting like CDM forever.)

To solve equation (7) numerically, we change the independent variable from $`t`$ to $`\mathrm{ln}a`$. Using the properties

$$\frac{d}{dt}=H\frac{d}{d\mathrm{ln}a},\ddot{\delta }_k=H^2\delta ^{\prime \prime }+\frac{1}{2}(H^2)^{}\delta ^{},$$
(11)
where $`{}_{}{}^{}d/d\mathrm{ln}a`$ and

$$\xi \frac{(H^2)^{}}{2H^2}=\frac{3}{2}\left(1+(1/\mathrm{\Omega }_m^{}1)a^{3(1+\alpha )}\right)^1,$$
(12)
equation (7) takes the form

$`\delta _k^{\prime \prime }`$ $`+`$ $`[2+\xi 3(2wc_s^2)]\delta _k^{}`$ (13)
$`=`$ $`\left[{\displaystyle \frac{3}{2}}(16c_s^2+8w3w^2)\left({\displaystyle \frac{kc_s}{aH}}\right)^2\right]\delta _k.`$ (14)
Even before solving this, it s obvious that a non-zero sound speed, if present for a sufficiently long time-span, is going to have a dramatic effect on the $`k`$-dependence of the perturbation growth. If $`c_s^2>0`$, then fluctuations with wavelength below the Jeans scale $`\lambda _J=\sqrt{\pi |c_s^2|/G\rho }`$ will be pressure-supported and oscillate rather than grow. This oscillation is confirmed by the numerical solutions, and is analogous to the acoustic oscillations in the photon-baryon fluid in the pre-decoupling epoch. If $`c_s^2<0`$, corresponding to negative $`\alpha `$, fluctuations below this wavelength will be violently unstable and grow exponentially .

A key point which has apparently been overlooked in prior work is that whereas all the other terms in equation (LABEL:firstorder) are of order unity or smaller, the sound speed term $`(kc_s/aH)^2`$ can be much larger even if the sound speed is tiny $`|c_s|1`$. This is because $`c_s`$ is multiplied by the prefactor $`k/aH`$ which can be enormous, since it it the Horizon scale divided by the perturbation scale. Defining a critical wavelength $`\lambda _c`$ by

$$\lambda _c^2\frac{c_s^2}{(aH)^2}=\frac{\alpha w}{(aH)^2},$$
(16)
the pressure term in equation (LABEL:firstorder) becomes simply $`(\lambda _ck)^2`$, so we expect oscillations or exponential blowup in the power spectrum on scales $`k\mathrm{}>\lambda _c^1`$. These are created mainly during the recent transition period when both $`a`$ and $`w`$ are or order unity (growing from 0 to 1), and since neither effect is seen in observed data, we therefore expect to obtain constraints of order $`|\alpha |\mathrm{}<(H/k)^2`$, the squared ratio of the perturbation scale to the horizon scale. This heuristic argument thus suggests that Galaxy clustering constraints on scales down to $`10h^1\mathrm{Mpc}`$ would give the constraint $`|\alpha |\mathrm{}<(10h^1\mathrm{Mpc}/3000h^1\mathrm{Mpc})^210^5`$ — we will see that this approximation is in fact fairly accurate.

For our numerical calculations, we evolved a scale invariant Harrison-Zeldovich spectrum up to redshift $`z=100`$ (before which the GCG is indistinguishable from $`\mathrm{\Lambda }`$CDM) with CMBfast  to correctly include all the relevant effects (early super-horizon evolution, pre-recombination acoustic oscillations, Silk-damping, etc..), with cosmological parameters given by the concordance model of . We then used equation (LABEL:firstorder) to evolve the fluctuations from $`z=100`$ until today. Results for a sample of $`\alpha `$-values are plotted in Figure 1, and show how tiny non-zero values of $`\alpha `$ result in large changes on small scales as expected.

We constrain $`\alpha `$ by making a $`\chi ^2`$ fit of the theoretically predicted power spectrum against that observed with the 2dF 100k Galaxy Redshift Survey as analyzed by . For each $`\alpha `$, we use the best fitting power normalization to ensure that our constraints come only from the shape of the power spectrum, not from the overall amplitude which involves mass-to-light bias. To be conservative and stay firmly in the linear regime, we discard data with $`k>0.3h/`$Mpc. We run our code for a fine grid of models with $`1<\alpha <1`$ to find the corresponding $`\chi ^2`$ values. The likelihood function $`e^{\mathrm{\Delta }\chi ^2/2}`$ is plotted in Figure 2. It predictably peaks around $`\alpha 0`$, and the observed skewness is simply due to the fact that the oscillating solution ($`\alpha >0`$)is easier to fit than the exponentially unstable solution ($`\alpha <0`$). Setting $`\mathrm{\Delta }\chi ^2=1`$ cutoff as in a crude Bayesian analysis gives the constraints $`0.00000081<\alpha <0.0000079`$.

To place this result in context, Figure 3 shows the 0th order constraints from Makler et al.  with our new strict constraints superimposed.

## IV Conclusions

Above we showed that GCG models with $`|\alpha |10^5`$ are ruled out by observation, since they cause fluctuations or blowup in the matter power spectrum that are not observed. Let us now examine the assumptions that went into this and the broader implications.

First of all, our extremely tight constraints imply that that the narrow range of allowed GCG models will be completely indistinguishable from $`\mathrm{\Lambda }`$CDM to both to 0th order and at the early times when primary CMB anisotropies are produced. This means that the corresponding standard constraints on cosmological parameters from CMB, SN Ia etc. apply also to the GCG models making the identification $`\mathrm{\Omega }_m^{}=\mathrm{\Omega }_m`$, so that there are no interesting degeneracies between $`\alpha `$ and other parameters that can significantly widen the allowed $`\alpha `$-range. We have therefore used standard constraints $`0.2<\mathrm{\Omega }_m<0.4`$ for the allowed region in Figure 3.

Second, our limiting our constraints to linear scales $`k<0.3h/\mathrm{Mpc}`$ was probably overly conservative. As reviewed in , there are quite strong constraints on the linear power spectrum on much smaller scales from weak lensing, from the Lyman $`\alpha `$ forest and perhaps even from lensing of halo substructure  which if used would tighten our upper limit on $`\alpha `$ to $`10^6`$, $`10^7`$ and $`10^{10}`$, respectively.

Third, we saw that all that really mattered in equation (LABEL:firstorder) as far as the constraints were concerned was the pressure term $`(\lambda _ck)^2`$. This means that our results apply more generally than merely to the GCG case: any unified dark matter model where $`p`$ is a unique function of $`\rho `$ is ruled out if the sound speed is non-negligible, i.e., if the function $`p(\rho )`$ departs substantially from a constant over the range where pressure has an effect — quantitatively, if $`|d\mathrm{ln}p/d\mathrm{ln}\rho |\mathrm{}>10^5`$, again rendering it indistinguishable from standard $`\mathrm{\Lambda }`$CDM.

In contrast, standard quintessence models have no such problems. Although they typically have high sound speeds causing oscillations as above, this does not prevent the dark matter from clustering since it is a separate dynamic component. Quintessence models would fail as above if there the two components were tightly coupled, and this is effectively what happens with UDM since the two are one and the same substance.

To salvage the UDM idea in some form, its pressure must not be uniquely determined by its density — not even on subhorizon scales. As worked out in detail by Hu , the effective sound speed can under some circumstances differ from the adiabatic sound speed, and it is only the former that must approximately vanish to satisfy our constraints. Although it may be possible to concoct such models, say by introducing another physical field upon which $`p`$ depends and making it fluctuate in such a way as to cancel the problematic pressure gradients, this would be giving up much of the elegance and simplicity that gave the unified dark matter idea its appeal, essentially substituting one extra field for another.

In conclusion, precision data is gradually allowing us to test rather than assume the physics underlying modern cosmology. We have taken a step in this direction by ruling out a broad class of so-called unified dark matter models. Our results indicate that dark energy is either indistinguishable from a pure cosmological constant or a separate component from the dark matter with a life of its own.

The authors wish to thank Raul Abramo for helpful comments. This work was supported by NSF grants AST-0071213, AST-0134999, AST-0098606 and PHY-0116590, NASA grants NAG5-9194 & NAG5-11099, MT and MZ are David and Lucile Packard Foundation fellows and MT is a Cottrell Scholar of Research Corporation.