Free Access
Issue
A&A
Volume 541, May 2012
Article Number A10
Number of page(s) 5
Section Extragalactic astronomy
DOI https://doi.org/10.1051/0004-6361/201218904
Published online 18 April 2012

© ESO, 2012

1. Introduction

M 82 is a nearby edge-on starburst galaxy in the M 81–M 82 group. Because of its proximity (3.53 Mpc, 1″ = 17 pc; Karachentsev et al. 2002), M 82 has been observed over a wide range of wavelengths (e.g., X-ray: Strickland et al. 2004; UV: Hoopes et al. 2005; optical: Ohyama et al. 2002; infrared: Engelbracht et al. 2006; Kaneda et al. 2010; Roussel et al. 2010; radio: Yun et al. 1994; Walter et al. 2002). The starburst activity in M 82 was presumably triggered by close encounter with M 81 about 100 Myr ago (Yun et al. 1993). As a result of the starburst activity, M 82 shows prominent galactic superwind, pushing out a significant amount of the interstellar medium from the galactic disk. The outflows are observed in ionized gas (e.g., Hα: Ohyama et al. 2002; X-ray: Strickland et al. 2004) with filamentary structures, which entrains neutral gas (e.g., CO: Walter et al. 2002; H2: Veilleux et al. 2009) and dust (e.g., Alton et al. 1999; Engelbracht et al. 2006; Leeuw & Robson 2009; Kaneda et al. 2010; Roussel et al. 2010; Yoshida et al. 2011).

Polycyclic aromatic hydrocarbons (PAHs) show prominent emission features in the mid-infrared (mid-IR), which are important to characterize the interstellar medium especially in star-forming galaxies (e.g., Lu et al. 2003; Smith et al. 2007). With Spitzer, Engelbracht et al. (2006) showed that the emission of PAH extends to 6 kpc from the galactic plane. With AKARI, Kaneda et al. (2010) found a tight correlation between the PAH and the Hα emission in the halo of the galaxy. Since the plasma responsible for the Hα emission shows outflows with a radial velocities of  ~600 km   s-1 (Shopbell & Bland-Hawthorn 1998), the PAHs are also very likely to be flowing out of the galaxy. The PAHs are, however, expected to be destroyed quite easily in X-ray hot plasma with temperature of  ≳ 106 K, since PAHs are small (~6 Å; Tielens 2008) particles with two-dimensional structures. For example, Micelotta et al. (2010) calculated that the timescale of the PAH destruction by collisions with electrons in the X-ray superwind of M 82 is less than a thousand years. Therefore it is unlikely that the PAHs coexist with the X-ray hot plasma in the halo. The PAHs are probably within dense molecular clouds, protected against the destruction (Veilleux et al. 2009), or in a much cooler (Hα) plasma phase even if they coexist with the plasma (Kaneda et al. 2010).

In this paper, we present the near-IR (2.5–4.5 μm) spectra of M 82 obtained with the Infrared Camera (IRC; Onaka et al. 2007) on board AKARI (Murakami et al. 2007). The spectral range includes the PAH 3.3 μm emission feature and the 3.4–3.6 μm features. We clearly detect these features from the halo of M 82. Both features are attributed to C-H vibration modes of carbonaceous grains; the former is caused by aromatic hydrocarbons, the latter by aliphatic hydrocarbons (Duley & Williams 1981). Based upon these features, we discuss the properties of carbonaceous grains including PAHs in the galactic superwind of M 82.

2. Observations and data reductions

The IRC spectroscopic observations of M 82 were carried out in October 2007 and August 2008. A summary of the observations is listed in Table 1. All data were taken from the AKARI archives; the first three observational data are publicly available on the official site1, while the last two are not because they were taken during the instrumental calibration time. To obtain the near-IR (2.5–4.5 μm) spectra, the grism spectroscopic mode was used in the observations. Although the data covered the wavelength range of 2.5–5.0 μm, we did not use the range of 4.5–5.0 μm, since the signal-to-noise ratios were relatively low at these wavelengths. We extracted spectral data from the two slits Ns and Nh with the sizes of 5″  ×  48″ and 3″  ×  60″, respectively (Ohyama et al. 2007). They have different spectral resolutions for diffuse sources (Ns: R ~ 100, Nh: R ~ 150) due to the difference in the slit width. Depending on the signal intensity, we divided the slit aperture area into two, three, or six and integrated the sub-slit data to derive a spectrum. As a result, we obtained 34 spectra from the whole galaxy. Figure 1 shows the positions of the slit sub-apertures on the AKARI 7 μm band map, from which we created spectra. The AKARI 7 μm band covers the wavelength range of 5.9–8.4 μm (Onaka et al. 2007), which includes the PAH 6.2 and 7.7 μm features. Hereafter, we define the center, disk, and halo regions by the areas demarcated by the 7 μm surface brightness levels of 40% and 0.32% of the peak, which correspond to the third and tenth contour from the top in Fig. 1. According to this definition, sub-aperture regions A2-4 belong to the center, regions A1, A5-6, B1-6, C1-6, D1-3, and E1-2 to the disk, and regions E3, F1-2, G1-2, H1-2, I1-2, and J1-2 to the halo. The basic spectral analyses were performed with the official pipeline prepared for reducing phase 3 data (IRC Spectroscopy Toolkit for Phase 3 data Version 20110301). In addition to the basic pipeline process, we applied the custom procedures to improve signal-to-noise ratios for each spectrum as described in Yamagishi et al. (2011). Finally, we applied smoothing with a boxcar kernel of  ~0.05 μm (R ~ 60) to the Ns and Nh spectra to mitigate the effect of the difference in the slit width (i.e. resolution) between the Ns (R ~ 100) and Nh (R ~ 150) slits.

Table 1

Summary of the observations and sub-slit apertures used in the present study.

thumbnail Fig. 1

Positions of the sub-slit apertures used to derive the 34 spectra in Fig. 2, overlaid on the AKARI 7 μm band contour map (Kaneda et al. 2010). North is up and west to the right. The size of each box represents the area of each sub-slit aperture, which is listed in Table 1. The name of each region is indicated together, where the incremental order of the number is from the north to the south.

3. Results

Figure 2 shows the spectra obtained from various regions of M 82. Many of the spectra show the aromatic hydrocarbon feature at 3.3 μm and the aliphatic features at 3.4–3.6 μm. In particular, the spectrum in the halo of region F1 clearly shows the aromatic emission although the region is located at the distance of about 2 kpc away from the galactic center. To quantify the 3.3–3.6 μm features, we approximated a linear baseline that is determined at the wavelength ranges of 3.0–3.2 and 3.7–3.9 μm on both sides of the features, shown by the blue lines in Fig. 2. Then, for the halo spectra, we calculated a standard deviation from a linear baseline at 3.8–4.5 μm and adopted the 5σ error in evaluating the significance of the 3.3–3.6 μm excess emission. For example, the spectrum in region F1 shows significance of 12σ for the 3.3–3.6 μm excess emission. As a result, the excess emission is significantly detected from the spectra of all regions except F2, H2, I1, and I2.

thumbnail Fig. 2

AKARI 2.5–4.5 μm spectra taken from the 34 regions in M 82, the positions of which are shown in Fig. 1. The blue linear baselines are determined at the wavelength ranges of 3.0–3.2 μm and 3.7–3.9 μm. For the spectra with significant detection of the excess emission, the red curves indicate the best-fit Drude profiles representing the contribution of the aromatic 3.3 μm emission, which are determined at the wavelength range of 3.20–3.35 μm.

The PAH emission extended in the halo was already reported from the mid-IR imaging observations (Engelbracht et al. 2006; Kaneda et al. 2010). Spectroscopically, however, Engelbracht et al. (2006) confirmed only the PAH 17 μm features in the halo region. Therefore we show for the first time the aromatic 3.3 μm emission in the halo. Near the center region, hydrogen recombination lines (Brα at 4.05 μm and Brβ at 2.62 μm) and absorption features of H2O ice at 3.05 μm and CO2 ice at 4.27 μm are clearly detected in the spectra. The strong recombination lines are consistent with the well-known intense starburst activity in the center of M 82, while spatially resolved H2O and CO2 ices are another new finding, although the H2O ice itself was already reported from ISO data (Sturm et al. 2000). These will be summarized in a separate paper by using all the AKARI spectra of the ices in nearby galaxies.

It appears in Fig. 2 that compared with the aromatic 3.3 μm emission, the aliphatic 3.4–3.6 μm emission is stronger in the halo than in the center and disk. Hence we calculated the fractions of the aliphatic features to the total excess (=aromatic + aliphatic features). In the calculation, the contribution of the aliphatic features was determined by subtracting the Drude profile, which represents the contribution of the aromatic 3.3 μm emission (Smith et al. 2007), from the total excess flux above the linear continuum between 3.3–3.6 μm. We used the wavelength range of 3.20–3.35 μm to fit the Drude profile; the central wavelength of the Drude profile was fixed at 3.295 μm. We first fitted the profile to the spectrum of region A3 at the galactic center and then adopted the best-fit width to the other spectra. The best-fit Drude profile to each spectrum is shown in the red curves in Fig. 2. We considered the errors given by the official pipeline plus the systematic errors originating in the baseline continuum; the latter errors were evaluated by changing the continuum fitting regions with shifts of  ± 0.05 μm.

Figure 3a shows the ratio of the aliphatic features to the total excess flux plotted against the distance from the galactic center. We measured the distance from the peak in the 7 μm emission to the center of each sub-aperture. The figure shows a significant correlation between the ratio and distance with a linear correlation coefficient of R =  + 0.75; the ratio increases with the distance from the galactic center. Hence Fig. 3a reveals that the strength of the aliphatic features relative to the aromatic feature is increasing in the halo of M 82. The colors denote north or south positions with respect to the galactic disk, i.e. the major axis of M 82 with the position angle of 64° (Mayya et al. 2005). In Fig. 3a, the ratios in the northern halo are higher than those in the southern halo, while there is no such systematic difference between the north and south near the center.

thumbnail Fig. 3

a) Ratio of the aliphatic features to the total excess (=aromatic + aliphatic features) flux plotted against the distance from the center of the galaxy. The data points in the halo are labeled with the region names. The diamonds, squares, and asterisks represent the regions in the center, disk, and halo, respectively. The colors of the symbols are changed depending on the location; the black for the on-plane (A3, B3, C3-4, and E1), blue for the north side (A1-2, B1-2, C1-2, and J1-2), and red for the south side from the galactic plane (A4-6, B4-6, C5-6, D1-3, E2-3, F1-2, G1-2, H1-2, and I1-2). b) Ratio of the aromatic 3.3 μm feature intensity to the AKARI 7 μm surface brightness (Kaneda et al. 2010) plotted against the distance from the center of the galaxy. The ratios are normalized to unity in the galactic center of region A3. The definition of the symbols and colors is the same as Fig. 3a. The data points with non-detection of the PAH 3.3 μm emission are shown with the 3σ upper limits.

4. Discussion

In general, the 3.3 μm feature is sensitive to PAHs of small sizes compared with the features at longer wavelengths (Schutte et al. 1993). Our result on the detection of the 3.3 μm feature from the halo spectra unambiguously confirms the presence of very small PAHs in the galactic superwind of M 82. With AKARI near- to mid-IR spectroscopy, Onaka et al. (2010) also detected the PAH features including the 3.3 μm feature in the filamentary outflow of the starburst dwarf galaxy NGC 1569. Thus, to our knowledge, this is the second clear detection of the PAH 3.3 μm feature in galactic superwind far outside a galactic disk. Onaka et al. (2010) suggested that the PAHs in the filament of NGC 1569 might be produced from the fragmentation of large carbonaceous grains in shocks. Our halo spectra in Fig. 2 exhibit unusually high ratios of the aliphatic 3.4–3.6 μm features to the aromatic 3.3 μm feature, supporting the picture that hydrocarbon grains are highly processed in the halo.

We compared the distribution of the PAH 3.3 μm feature with that of the PAHs features in the AKARI 7 μm band. As described above, the AKARI 7 μm band includes the PAH 6.2 and 7.7 μm features. Figure 3b displays the ratio of the 3.3 μm feature to the AKARI 7 μm band brightness plotted against the distance from the galactic center. Here we considered that the 3.3 μm feature intensities from the spectra taken in the instrumental calibration time could be systematically overestimated by 20% (Mori et al. 2011), and therefore we added the 20% errors to the lower sides of the corresponding errorbars. Considering that the surface brightness in the 7 μm band decreases by three orders of magnitude from the center to the halo, the ratios are strikingly constant over such a wide brightness range. The result therefore indicates that the unusual hydrocarbon spectra at 3.3–3.6 μm are caused by unusually strong aliphatic 3.4–3.6 μm features, but not by an unusually faint aromatic 3.3 μm feature. The systematic difference between the northern and southern disk regions at  <0.5 kpc from the center is most probably caused by an extinction effect; it is known that M 82 is highly inclined at  ~80° and the north side of the galaxy is close to us (Lynds & Sandage 1963; Shopbell & Bland-Hawthorn 1998). Thus the southern disk regions are more affected by the interstellar extinction than the northern disk regions. Considering the extinction effect, the ratios tend to be more constant throughout the galaxy.

The halo spectra in Fig. 2, especially those in regions J1 and J2, show remarkable resemblance to the spectra obtained at the molecular loop near the center of our Galaxy where the CO radio observations showed that the violent motion and the shock heating of gas took place (Fukui et al. 2006). Kaneda et al. (2012) found that the loop spectra have unusual properties at 3.3–3.6 μm similarly to our halo spectra, i.e., faint 3.3 μm emission and broad excess above the continuum, from which they suggested that the destruction of small PAHs and the shattering of carbonaceous grains at the foot point of the molecular loop. In M 82, Kaneda et al. (2010) and Roussel et al. (2010) already reported the far-IR emission widely extended in the halo, showing large cold grains entrained in the galactic superwind. Hence it is possible that the unusually strong aliphatic 3.4–3.6 μm features in the halo are caused by fragmentation of larger grains into small carbonaceous ones. At the same time, since we do detect the PAH features in the halo, a significant amount of PAHs including very small ones must be protected in dense molecular clouds against shocks in the halo. Alternatively, Kwok & Zhang (2011) pointed out that mixed aromatic-aliphatic organic nanoparticles can be carriers of both 3.3 and 3.4–3.6 μm features. In this case, the observed unusual features may simply reflect the processing of these nanoparticles.

As can be seen in Fig. 3a, the spectra in the northern halo (J1 and J2) show higher ratios of the aliphatic to the aromatic emission than those in the southern halo (E3, F1, G1-2, and H1). This result may suggest that the hydrocarbon grains are exposed to different interstellar conditions between the northern and southern halos. Comparing the spatial distribution of the H2v = 1−0 S(1) emission at 2.12 μm with that of the PAH emission, Veilleux et al. (2009) pointed out that the dominant excitation mechanism of molecular hydrogen may be different between the northern and southern halos. They suggested that molecular hydrogen may be excited mostly by UV radiation in the southern halo because they have similar distributions there; it is also expected that photoionization by OB stars is dominant in the southern halo region (Shopbell & Bland-Hawthorn 1998). In contrast, their distributions are significantly different in the northern halo (Veilleux et al. 2009). Our result in Fig. 3a favors the scenario that shocks are more important for the heating of molecular hydrogen in the northern halo where hydrocarbon grains are more efficiently processed and destroyed.

5. Conclusions

With AKARI, we have performed near-IR (2.5–4.5 μm) spectroscopic observations of the nearby edge-on starburst galaxy M 82. We obtain 34 spectra in the galaxy including the northern and southern halos to investigate the properties of carbonaceous grains in the galactic superwind. In many spectra, aromatic 3.3 μm emission and aliphatic 3.4–3.6 μm emission are clearly detected. In particular, the aromatic and aliphatic emissions are detected in halo regions 2 kpc far away from the galactic center, which clearly demonstrates that even small PAHs survive in a harsh environment of the galactic superwind. We find that the fraction of the aliphatic to the aromatic feature increases with the distance from the galactic center. At the same time, we show that the strength of the PAH 3.3 μm feature relative to that of the PAH 6.2 and 7.7 μm features does not decrease toward the

halo. Therefore we conclude that the unusual hydrocarbon spectra at 3.3–3.6 μm are caused by unusually strong aliphatic 3.4–3.6 μm features, but not by an unusually faint aromatic 3.3 μm feature. The result suggests that in the galactic superwind, small carbonaceous grains emitting the aliphatic features are produced by fragmentation of larger grains in shocks, while an appreciable amount of small PAHs are protected in dense molecular clouds against shocks.


Acknowledgments

We thank all the members of the AKARI project. We also express many thanks to the anonymous referee for the useful comments. AKARI is a JAXA project with the participation of ESA. This work is supported by Grants-in-Aid for Japan Society for the Promotion of Science Fellows No. 23005457 and Research No. 22340043.

References

  1. Alton, P. B., Davies, J. I., & Bianchi, S. 1999, A&A, 343, 51 [NASA ADS] [Google Scholar]
  2. Duley, W. W., & Williams, D. A. 1981, MNRAS, 196, 269 [NASA ADS] [CrossRef] [Google Scholar]
  3. Engelbracht, C. W., Kundurthy, P., Gordon, K. D., et al. 2006, ApJ, 642, L127 [NASA ADS] [CrossRef] [Google Scholar]
  4. Fukui, Y., Yamamoto, H., Fujishita, M., et al. 2006, Science, 314, 106 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  5. Hoopes, C. G., Heckman, T. M., Strickland, D. K., et al. 2005, ApJ, 619, L99 [NASA ADS] [CrossRef] [Google Scholar]
  6. Kaneda, H., Ishihara, D., Suzuki, T., et al. 2010, A&A, 514, A14 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  7. Kaneda, H., Ishihara, D., Mouri, A., et al. 2012, PASJ, 64, 25 [Google Scholar]
  8. Karachentsev, I. D., Dolphin, A. E., Geisler, D., et al. 2002, A&A, 383, 125 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  9. Kwok, S., & Zhang, Y. 2011, Nature, 479, 80 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  10. Leeuw, L. L., & Robson, E. I. 2009, AJ, 137, 517 [NASA ADS] [CrossRef] [Google Scholar]
  11. Lu, N., Helou, G., Werner, M. W., et al. 2003, ApJ, 588, 199 [NASA ADS] [CrossRef] [Google Scholar]
  12. Lynds, C. R., & Sandage, A. R. 1963, ApJ, 137, 1005 [NASA ADS] [CrossRef] [Google Scholar]
  13. Mayya, Y. D., Carrasco, L., & Luna, A. 2005, ApJ, 628, L33 [NASA ADS] [CrossRef] [Google Scholar]
  14. Micelotta, E. R., Jones, A. P., & Tielens, A. G. G. M. 2010, A&A, 510, A37 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  15. Mori, T., Kaneda, H., Ishihara, D., Oyabu, S., & Wada, T. 2011, PASP, 123, 942 [NASA ADS] [CrossRef] [Google Scholar]
  16. Murakami, H., Baba, H., Barthel, P., et al. 2007, PASJ, 59, 369 [Google Scholar]
  17. Ohyama, Y., Taniguchi, Y., Iye, M., et al. 2002, PASJ, 54, 891 [NASA ADS] [CrossRef] [Google Scholar]
  18. Ohyama, Y., Onaka, T., Matsuhara, H., et al. 2007, PASJ, 59, 411 [NASA ADS] [Google Scholar]
  19. Onaka, T., Matsuhara, H., Wada, T., et al. 2007, PASJ, 59, 401 [NASA ADS] [Google Scholar]
  20. Onaka, T., Matsumoto, H., Sakon, I., & Kaneda, H. 2010, A&A, 514, A15 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  21. Roussel, H., Wilson, C. D., Vigroux, L., et al. 2010, A&A, 518, L66 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  22. Schutte, W. A., Tielens, A. G. G. M., & Allamandola, L. J. 1993, ApJ, 415, 397 [NASA ADS] [CrossRef] [Google Scholar]
  23. Shopbell, P. L., & Bland-Hawthorn, J. 1998, ApJ, 493, 129 [NASA ADS] [CrossRef] [Google Scholar]
  24. Smith, J. D. T., Draine, B. T., Dale, D. A., et al. 2007, ApJ, 656, 770 [NASA ADS] [CrossRef] [Google Scholar]
  25. Strickland, D. K., Heckman, T. M., Colbert, E. J. M., Hoopes, C. G., & Weaver, K. A. 2004, ApJS, 151, 193 [NASA ADS] [CrossRef] [Google Scholar]
  26. Sturm, E., Lutz, D., Tran, D., et al. 2000, A&A, 358, 481 [NASA ADS] [Google Scholar]
  27. Tielens, A. G. G. M. 2008, ARA&A, 46, 289 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  28. Veilleux, S., Rupke, D. S. N., & Swaters, R. 2009, ApJ, 700, L149 [NASA ADS] [CrossRef] [Google Scholar]
  29. Walter, F., Weiss, A., & Scoville, N. 2002, ApJ, 580, L21 [NASA ADS] [CrossRef] [Google Scholar]
  30. Yamagishi, M., Kaneda, H., Ishihara, D., et al. 2011, ApJ, 731, L20 [NASA ADS] [CrossRef] [Google Scholar]
  31. Yoshida, M., Kawabata, K. S., & Ohyama, Y. 2011, PASJ, 63, 493 [Google Scholar]
  32. Yun, M. S., Ho, P. T. P., & Lo, K. Y. 1993, ApJ, 411, L17 [NASA ADS] [CrossRef] [Google Scholar]
  33. Yun, M. S., Ho, P. T. P., & Lo, K. Y. 1994, Nature, 372, 530 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]

All Tables

Table 1

Summary of the observations and sub-slit apertures used in the present study.

All Figures

thumbnail Fig. 1

Positions of the sub-slit apertures used to derive the 34 spectra in Fig. 2, overlaid on the AKARI 7 μm band contour map (Kaneda et al. 2010). North is up and west to the right. The size of each box represents the area of each sub-slit aperture, which is listed in Table 1. The name of each region is indicated together, where the incremental order of the number is from the north to the south.

In the text
thumbnail Fig. 2

AKARI 2.5–4.5 μm spectra taken from the 34 regions in M 82, the positions of which are shown in Fig. 1. The blue linear baselines are determined at the wavelength ranges of 3.0–3.2 μm and 3.7–3.9 μm. For the spectra with significant detection of the excess emission, the red curves indicate the best-fit Drude profiles representing the contribution of the aromatic 3.3 μm emission, which are determined at the wavelength range of 3.20–3.35 μm.

In the text
thumbnail Fig. 3

a) Ratio of the aliphatic features to the total excess (=aromatic + aliphatic features) flux plotted against the distance from the center of the galaxy. The data points in the halo are labeled with the region names. The diamonds, squares, and asterisks represent the regions in the center, disk, and halo, respectively. The colors of the symbols are changed depending on the location; the black for the on-plane (A3, B3, C3-4, and E1), blue for the north side (A1-2, B1-2, C1-2, and J1-2), and red for the south side from the galactic plane (A4-6, B4-6, C5-6, D1-3, E2-3, F1-2, G1-2, H1-2, and I1-2). b) Ratio of the aromatic 3.3 μm feature intensity to the AKARI 7 μm surface brightness (Kaneda et al. 2010) plotted against the distance from the center of the galaxy. The ratios are normalized to unity in the galactic center of region A3. The definition of the symbols and colors is the same as Fig. 3a. The data points with non-detection of the PAH 3.3 μm emission are shown with the 3σ upper limits.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.