Äîêóìåíò âçÿò èç êýøà ïîèñêîâîé ìàøèíû. Àäðåñ îðèãèíàëüíîãî äîêóìåíòà : http://www.arcetri.astro.it/~sperello/AA518_L50.pdf Äàòà èçìåíåíèÿ: Mon Jul 19 18:08:28 2010 Äàòà èíäåêñèðîâàíèÿ: Tue Oct 2 02:17:51 2012 Êîäèðîâêà: Ïîèñêîâûå ñëîâà: m 103

A&A 518, L50 (2010) DOI: 10.1051/0004-6361/201014533
c ESO 2010

Astronomy & Astrophysics
Special feature
Letter to the Editor

Herschel: the first science highlights

The Herschel Virgo Cluster Survey
III. A constraint on dust grain lifetime in early-type galaxies
M. S. Clemens1 , A.P.Jones2 , A. Bressan1 , 3 , M. Baes4 ,G.J. Bendo5 ,S.Bianchi6 ,D.J. Bomans7 , A. Boselli8 , E. Corbelli6 ,L.Cortese9 , A. Dariush9 , J. I. Davies9 , I. De Looze4 , S. di Serego Alighieri6 ,D.Fadda10 , J. Fritz4 , D. A. Garcia-Appadoo11 , G. Gavazzi12 , C. Giovanardi6 , M. Grossi13 , T. M. Hughes9 , L.K.Hunt6 ,S.Madden14 , D. Pierini15 , M. Pohlen9 ,S.Sabatini16 ,M.W.L.Smith9 , J. Verstappen4 , C. Vlahakis17 , E. M. Xilouris18 , and S. Zibetti19
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19

INAF-Osservatorio Astronomico di Padova, Vicolo dell'Osservatorio 5, 35122 Padova, Italy e-mail: marcel.clemens@oapd.inaf.it Institut d'Astrophysique Spatiale (IAS), BÁtiment 121, UniversitÈ Paris-Sud 11 and CNRS, 91405 Orsay, France e-mail: Anthony.Jones@ias.u-psud.fr SISSA-ISAS, International School for Advanced Studies, via Beirut 4, 34014 Trieste, Italy Sterrenkundig Observatorium, Universiteit Gent, Krijgslaan 281 S9, 9000 Gent, Belgium Astrophysics Group, Imperial College London, Blackett Laboratory, Prince Consort Road, London SW7 2AZ, UK INAF-Osservatorio Astrofisico di Arcetri, Largo Enrico Fermi 5, 50125 Firenze, Italy Astronomical Institute, Ruhr-University Bochum, Universitaetsstr. 150, 44780 Bochum, Germany Laboratoire d'Astrophysique de Marseille, UMR 6110 CNRS, 38 rue F. Joliot-Curie, 13388 Marseille, France Department of Physics and Astronomy, Cardiff University, The Parade, Cardiff, CF24 3AA, UK NASA Herschel Science Center, California Institute of Technology, MS 100-22, Pasadena, CA 91125, USA ESO, Alonso de Cordova 3107, Vitacura, Santiago, Chile UniversitÞ di Milano-Bicocca, piazza della Scienza 3, 20100, Milano, Italy CAAUL, ObservatÑrio AstronÑmico de Lisboa, Universidade de Lisboa, Tapada da Ajuda, 1349-018, Lisboa, Portugal Laboratoire AIM, CEA/DSM ­ CNRS ­ UniversitÈ Paris Diderot, Irfu/Service d'Astrophysique, 91191 Gif sur Yvette, France Max-Planck-Institut fuer extraterrestrische Physik, Giessenbachstrasse, Postfach 1312, 85741, Garching, Germany INAF-Istituto di Astrofisica Spaziale e Fisica Cosmica, via Fosso del Cavaliere 100, 00133, Roma, Italy Leiden Observatory, Leiden University, PO Box 9513, 2300 RA Leiden, The Netherlands National Observatory of Athens, I. Metaxa and Vas. Pavlou, P. Penteli, 15236 Athens, Greece Max-Planck-Institut fuer Astronomie, Koenigstuhl 17, 69117 Heidelberg, Germany

Received 29 March 2010 / Accepted 30 April 2010
ABSTRACT

Passive early-type galaxies (ETGs) provide an ideal laboratory for studying the interplay between dust formation around evolved stars and its subsequent destruction in a hot gas. Using Spitzer-IRS and Herschel data we compare the dust production rate in the envelopes of evolved AGB stars with a constraint on the total dust mass. Early-type galaxies which appear to be truly passively evolving are not detected by Herschel. We thus derive a distance independent upper limit to the dust grain survival time in the hostile environment of ETGs of <46 ± 25 Myr for amorphous silicate grains. This implies that ETGs which are detected at far-infrared wavelengths have acquired a cool dusty medium via interaction. Given likely time-scales for ram-pressure stripping, this also implies that only galaxies with dust in a cool (atomic) medium can release dust into the intra-cluster medium.
Key words. galaxies: elliptical and lenticular, cD ­ galaxies: ISM ­ infrared: galaxies

1. Introduction
In this paper we compare the dust outflow rate from the evolved AGB stars in passive ETGs, with the total dust mass, to estimate dust grain lifetimes. Both quantities are derived from observations; the outflow rate from Spitzer-IRS spectra, and the total dust mass from Herschel (Pilbratt et al. 2010) maps of the
Herschel is an ESA space observatory with science instruments provided by European-led Principal Investigator consortia and with important participation from NASA.

Herschel Virgo Cluster Survey, HeViCS (Davies et al. 2010, www.hevics.org). We assume that evolved stars are the only source of dust in these passive systems as there is very little evidence that type Ia supernovae produce dust. IRAS data have been used in the past in a similar way. Soifer et al. (1986) estimated the mass-loss from evolved stars in the bulge of M 31 from 12 and 25 m fluxes and the total dust mass from the 60 and 100 m fluxes. They concluded that given that the observed dust mass could be produced by stellar outflows in only 107 yr, there must be some mechanism that depletes the bulge of inter-stellar matter. Jura et al. (1987) made similar Page 1 of 4

Article published by EDP Sciences


A&A 518, L50 (2010)

calculations for a sample of elliptical galaxies. They estimated dust outflow rates by assuming that half of the 12 m flux came from dusty stellar envelopes. They also found that the stellar mass-loss rate would be sufficient to produce the observed cool ISM in much less than a Hubble time. Dust mass production rate and total dust mass can only reasonably be compared in objects that are truly passively evolving. Temi et al. (2007) find, that for a given blue luminosity, the 70 and 160 m luminosities vary by 2 orders of magnitude, indicating that many ETGs have a significant dust component not directly attributable to mass-loss from evolved stars. For our purposes, these ETGs are not "passive" because other processes (mergers?) have probably contributed to the dust mass. Recent studies of ETGs with Spitzer (Bressan et al. 2006; Clemens et al. 2009) have shown that even in samples of ETGs selected to be the most passive objects (i.e. lying on the colourmagnitude relation) a significant fraction show evidence of either on-going or recent past star formation. These results indicate that ETGs in which dust features are seen in the optical (e.g. Sadler & Gerhard 1985) and many that have been detected by IRAS at 60 and 100 m(Knapp et al. 1989)and Spitzer at 70 and 160 m (Kaneda et al. 2007) may also host low levels of star formation. Intriguigly, of the 7 ellipticals detected at 70 and 160 m by Kaneda et al. (2007) the only object not showing PAH emission features is actually a radio source, so that the far-infrared emission is probably synchrotron! Tsai & Mathews (1995) studied dust destruction in ETGs via thermal sputtering in the hot (106 -107 K), low density (nH 0.1 cm-3 ) gas, and found that the destruction time-scale is short compared to any cooling flow or dust transport time-scale. They therefore concluded that dust is destroyed "on the spot" before it has time to migrate within a galaxy. In a following work, Tsai & Mathews (1996) found dust-to-gas mass ratios that are orders of magnitude less than that in typical stellar ejecta or in the ISM of the Milky Way. Their interpretation of the IRAS data was, at that time, hampered by the resolution of the existing observations; they were therefore not able to tie down the origin and distribution of the observed dust emission. As part of the Herschel science demonstration phase (SDP), a 4 â 4 deg2 field, centred approximately on M 87, has been observed as part of the HeViCS with both the PACS (Poglitsch et al. 2010) and SPIRE (Griffin et al. 2010) instruments. See Davies et al. (2010) for details.

Because of the sensitivity of the emission at the shorter PACS wavelengths to small masses of relatively warm dust, we consider only the SPIRE images at 250, 350 and 500 m to determine the dust mass detection limit. Our mass detection limit depends on whether we consider the detection of an extended or point source. For passive galaxies, the dust should follow the stellar distribution, and thus be strongly centrally peaked. As the core radius of an elliptical galaxy at the distance of the Virgo cluster is typically less than the resolution of the SPIRE images (respectively 18 , 25 and 37 ) the most likely detection would take the form of a point source. Unless the dust is very cold ( 17 K) the tightest limit on the dust emission comes from the 250 m map. This has a point source detection limit of S rms = 9.6 mJy. We calculate the dust mass detection limit as (Hildebrand 1983): Md S rms ()D2 d () B(, T d) (1)

where D is the distance to the Virgo cluster (16.5 Mpc, Mei et al. 2007) and d () is the dust mass opacity coefficient. The value of d () depends on the grain composition but the mid-infrared spectra (Sect. 2.2) are consistent with pure silicate grains. For dust composed only of amorphous silicate grains d (250 m) = 0.6517 m2 kg-1 (Draine & Lee 1984). Adopting this value, and T d = 30 K, we arrive at a 2 mass detection limit of Md 8.7 â 103 M 2 .
2.2. Dust production rate from AGB stars

Assuming a spherically symmetric, stationary stellar wind of ve CE locity v , and dust mass-loss rate Md the circumstellar dust density is, CE = d CE Md 4r 2 v


·

(2)

The optical depth of the circumstellar envelope, , is obtained by multiplying Eq. (2) by the dust opacity coefficient and integrating from rin to rout , where the former is the dust sublimation radius and the latter is a convenient outer radius rin . We thus have, 1 = CE Md k1 4vrin (3)

2. Results Bressan et al. (2006) observed a sample of ETGs in the Virgo cluster with Spitzer-IRS. The objects were selected to lie on the optical colour-magnitude relation and were expected to be passively evolving objects. The mid-infrared spectra for this sample showed that only 76% were actually passive objects, with the remainder showing signs of either on-going or recent past star formation (PAH features) or AGN activity (atomic lines). Of this sample of 17 galaxies, 9 lie in the HeViCS SDP field. 6 of these show IRS spectra consistent with passively evolving stellar populations (with the possible exception of NGC 4371), and it is these 6 objects, listed in Table 1, on which we base our study.
2.1. Dust mass detection limit

where 1 is the dust opacity at 1m, following the notation of Bressan et al. (1998), and we have neglected the very small term 1/rout . The dust sublimation radius is rin L1/2 (Granato & Danese 1994) and the proportionality constant depends mainly on the dust composition (through its sublimation temperature), the size distribution and only weakly on the shape of the stellar spectrum. Inserting the value of rin into Eq. (3) we obtain, CE Md = A 1 v


L4

(4) L . Th e the mixBressan (Fig. 1),

Of the 6 passive ETGs in the HeViCS SDP field, none are detected at any wavelength (100, 160, 250, 350 and 500 m)1 .
We note that several ETGs are detected in the SDP field, including M 87, M 84, M 86, NGC 4459, NGC 4435 and NGC 4550, but none of these are truly passively evolving. Page 2 of 4
1

where L4 is the bolometric luminosity in units of 104 factor A is mainly determined by the dust mixture. For ture of silicate grains of M-type AGB stars used in et al. (1998), suitable for the spectra of passive ETGs A 6.5 â 10-10 M yr-1 /(km s-1 ).
2

Were the dust composed of amorphous carbon (Serra Diaz-Caßo & Jones 2008)then d (250 m) = 1.932 (Rouleau & Martin 1991)and the mass detection limit Md 2.9 â 103 M .


M. S. Clemens et al.: The Herschel Virgo Cluster Survey. III. Table 2. Parameter values for calculation of grain lifetime and the power of their dependence on the derived lifetime. S rms mJy 9.6 1 d (250) m2 kg-1 0.6517 ­1 Td K 30 ­1 A 10-10 6.5 1 1 fit 1 v km s- 10 1 L4 104 L 0.3 0.5

1

Notes. E.g., tg 1/T d . () The units of A are M yr-1 /(km s-1 ). () In the Rayleigh-Jeans limit; which is a bad approximation for T d 20 K.

"average" AGB star, LCE (10 m). Thus, the total dust mass-loss d sampled by the slit is, s Mdlit = A 1 v


Fig. 1. Spitzer-IRS spectrum of NGC 4473. The fitted model components are the stellar photosphere, and dusty circumstellar envelope (CE). The spectrum of NGC 4371 is shifted by 0.2 dex.

L4

s 4 D2 fd lit (10 m)

LCE (10 m) d

·

(5)

Table 1. Dust mass loss rates and grain lifetimes, tg . NGC 4371 4442 4473 4474 4551 4564 Md /L (K) 3. 3. 3. 3. 2. 4. 80 48 57 71 61 28


M yr-1 1026 WHz-1

Md /M 0. 0. 0. 0. 0. 0.



M yr-1 1012 M

00855 00712 00742 00806 00468 00959

Md (beam) 10-4 M yr- 1.0 1.7 1.9 0.59 0.36 1.3

1

tg 106 yr 86 51 46 150 240 65

s Summarizing, the spectral fit provides 1 and fd lit (10 m). From 1 , which characterizes the dusty envelope, and the assumed bolometric stellar luminosity, we obtain LCE (10 m) and finally d s Mdlit , for a specified value of v . We note that, since LCE (10 m) d scales approximately with the bolometric luminosity of the star, the total mass-loss has only a weak dependence on the assumed luminosity of the typical AGB star. We finally note that recent work (Boyer et al. 2010) has shown dust to be released into the ISM only by AGB stars, contrary to some claims that RGB stars also contribute.

2.3. Dust grain lifetimes

Notes. () For an assumed dust temperature of 30 K. () The specific dust mass-loss rate per unit K -band luminosity. () The mass loss per 1012 M of galaxy mass. A dust mass opacity coefficient suitable for amorphous silicate grains has been used.

To determine the dust mass-loss rate from the MIR spectrum we fit a combination of a pure photospheric atmosphere, a MARCS model (Gustafsson et al. 2008) of 4000 K, plus the emission from a dusty envelope selected from the library of Bressan et al. (1998), as described in Bressan et al. (2007). We s thus obtain the fraction of flux, say at 10 m, fd lit (10 m), due to dusty circumstellar envelopes, and, from the distance of the galaxy, the corresponding "dust" luminosity sampled by the slit, s 4 D2 fd lit (10 m). Surprisingly, though ETGs contain a mixture of evolved stars of different metallicity and age, a single dusty envelope is typically enough for a good fit to the MIR spectrum. The shape of the broad feature near 10 m in general requires 1 2­4, and so the envelopes are optically thin at 10 m. From the shape of the 10 m feature we thus derive the optical depth of the "average" AGB star. To compute its dust mass-loss rate through Eq. (4) we need to assume a bolometric luminosity and wind velocity. Typical values for evolved, old stars are L4 0.2­ 0.4 and v 5-15 km s-1 . The total dust mass-loss sampled by the slit is finally obtained by multiplying the value for the "average" AGB star by the number of dusty AGB stars sampled by the slit. This is given by the ratio of the dust luminosity of the galaxy to that of the

Dust mass-loss rates were derived from Spitzer-IRS spectra that were taken through a slit of dimensions 18 â 3. 6 centred on the galaxy. In order to compare these values with the dust masses within one 250 m beam we assume that the light distribution in the mid and far-infrared are similar, and that both are similar to that at K -band. For each galaxy, we use the 2MASS K -band image, smoothed to the Spitzer-IRS spatial resolution, to determine the flux within the IRS slit. We then smooth the K -band image to the Herschel 250 m resolution and calculate the flux within one beam. We use the ratio of these two fluxes to derive the mass-loss rate from the area of one Herschel beam. We derive dust grain lifetimes simply as the ratio of dust mass detection limit to dust mass-loss rate, Md / Md . Our dust 3 mass detection limit of Md 8.7 â 10 M for silicate grains results in grain lifetimes as given in Table 1. The tightest constraint, of <46 ± 25 Myr, is provided by NGC 4473. The error reflects the uncertainties in v and L4 .In Table 2 we provide the values used to arrive at this grain lifetime and their dependencies. For dust temperatures in the range 20-40 K, for example, the grain lifetime lies in the range 69-31 Myr. If dust grains are destroyed via sputtering in the 106 -107 K gas in ETGs then the theoretically, expected dust grain lifetime in years is t = ki rg /nH (e.g., Tielens et al. 1994), where rg is the grain radius in nm, nH is the gas density in cm-3 and ki = 310 for silicate dust (210 for amorphous carbon, 1500 for graphite) so the grain lifetime should be longer in galaxies with a lower density hot ISM. Of our sample, NGC 4371 and NGC 4474 both lie in the least X-ray luminous category as defined by Irwin & Sarazin (1998). NGC 4564 has been observed by Chandra, and the central hot gas density estimated to be 0.011 ± 0.006 cm-3 by Soria et al. (2006). These authors analyzed a sample of "quiescent" galaxies and found central gas densities in the range
Page 3 of 4


A&A 518, L50 (2010)

0.011­0.03 cm-3 . These values overlap with the low end of estimates from Chandra data of more general ETG samples (Pellegrini 2005). We thus assume a hot gas density of 0.02 cm-3 for all our ETGs. The expected dust grain lifetime, for "typical" 100 nm silicate grains is 1.55 Myr (1.05 for amorphous hydrocarbon, 7.5 for graphite) and thus consistent with our nondetections.

3. Discussion and conclusions
We interpret the ratio of dust mass to dust production rate as a measure of grain lifetime. In principle, however, dust could be removed from the galaxy rather than destroyed. A wind of hot gas produced by type Ia supernovae probably cannot drive dust (or gas) out of the potential well of massive ETGs (Mathews & Loewenstein 1986). However, good evidence of ram-pressure stripping of the hot ISM in elliptical galaxies has been found in the form of X-ray tails (Randall et al. 2008; Machacek et al. 2006). If dust is found in this hot gas then it too may be removed. However, the timescale for such stripping is estimated to be of the order a few 108 yr (Takeda et al. 1984; Murakami & Babel 1999, Eq. (12)) which is an order of magnitude longer than our upper limit to the grain lifetime. The observation that the X-ray emission of most ETGs is not displaced from the stellar distribution, in fact, suggestes that either the stripping time-scale is longer than the hot gas production time-scale or that the ram-pressure is typically insufficient to remove the hot medium. It is now clear that if dust is isolated from the hot medium in denser, cooler clouds, than ram-pressure does have time to strip dust (Cortese et al. 2010) but the passive ETGs in our study do not show any evidence of accumulating dust in such clouds. Passive ETGs do not directly pollute the intra-cluster medium with dust. Grain lifetimes of <46 Myr have immediate implications for ETGs that are detected in the FIR. Either they produce dust at rates much higher than those given in Table 1 or their dust is shielded from the hot gas in cool clouds. For example, the dust mass of NGC 4435 of 1.2 â 106 M would require a total dust production rate of 0.03 M yr-1 , more than an order of magnitude greater than any value in Table 1. This object is known to host star formation (Panuzzo et al. 2007) and so dust is likely to be in cool clouds more typical of late-type galaxies; grain lifetimes in this case may be much longer. Dust may have been acquired from an interaction, but in any case it probably arrived as part of a cooler medium. In a passive ETG, the ISM, and the dust within it, come mainly from mass-loss from evolved stars. Unless the dust-togas ratio in this released material is very variable, one would expect dust and gas production to follow one another. If the X-ray luminosity, Lx n2 one would naively expect a relae 2 tion Lx Md . However, if the grain lifetime, tg 1/nH , then we would expect only a linear correlation between Lx and Md . Kaneda et al. (2007) find evidence of an anti-correlation for a small sample of elliptical galaxies detected by Spitzer-MIPS, but this may be expected if the dust were of external origin.

Only NGC 4371 is detected by MIPS at 70 m (marginally at 160 m). The detection is consistent with our Herschel detection limits. This object actually shows extremely weak PAH features in its Spitzer-IRS spectrum, and is therefore probably not totally passive. The object is included for comparison, and its IRS spectrum included in Fig. 1. We have assumed that AGB stars are the only source of dust. Were type Ia supernovae able to produce dust from the metals they produce, we estimate that the contribution to the dust production rate would be, at most, similar to that of AGB stars. Our limit on grain lifetime would be a factor of 2 lower in this case. Although dust grains are rather short-lived in ETGs, the prospects for survival in the intra-cluster medium are rather better because the gas densities are typically at least 2 orders of magnitude lower than in ETGs (BÆhringer et al. 1994). Therefore, if dust can be removed from cluster galaxies, or form in the intra-cluster medium, there is a good chance that such dust will be detected by the HeViCS, especially at larger cluster radii where gas densities are lower.
Acknowledgements. M.C. and A.B. acknowledge support from contract ASI/INAF I/016/07/0.

References
BÆhringer, H., Briel, U. G., Schwarz, R. A., et al. 1994, Nature, 368, 828 Boyer, M. L., van Loon, J. Th., McDonald, I., et al. 2010, ApJ, 711, L99 Bressan, A., Granato, G. L., & Silva, L. 1998, A&A, 332, 135 Bressan, A., Panuzzo, P., Buson, L., et al. 2006, ApJ, 639, L55 Bressan, A., Panuzzo, P., Silva, L., et al., 2007, ASPC, 374, 333 Clemens, M. S., Bressan, A., Panuzzo, P., et al. 2009, MNRAS, 392, 982 Cortese, L., et al. 2010, A&A, 518, L49 Davies, J. I., et al. 2010, A&A, 518, L48 Draine, B. T., & Lee, H. M. 1984, ApJ, 285, 89 Granato, G. L., & Danese, L. 1994, MNRAS, 268, 235 Griffin, M. J., et al. 2010, A&A, 518, L3 Gustafsson, B., Edvardsson, B., Eriksson, K., et al. 2008, A&A, 486, 951 Hildebrand, R. H. 1983, QJRAS, 24, 267 Irwin, J. A., & Sarazin, C. L. 1998, ApJ, 499, 650 Jura, M., Kim, D. W., Knapp, G. R., & Guhathakurta, P. 1987, ApJ, 312, L11 Kaneda, H., Onaka, T., Kitayama, T., Okada, Y., & Sakon, I. 2007, PASJ, 59, 107 Knapp, G. R., Guhathakurta, P., Kim, D.-W., & Jura, M. A. 1989, ApJS, 70, 329 Machacek, M., Jones, C., Forman, W. R., & Nulsen, P. 2006, ApJ, 644, 155 Mathews, W. G., & Loewenstein, M. 1986, ApJ, 306, L7 Mei, S., Blakeslee, J. P., CÒtÈ, P., et al. 2007, ApJ, 655, 144 Murakami, I., & Babul, A. 1999, MNRAS, 309, 161 Panuzzo, P., Vega, O., Bressan, A., et al. 2007, ApJ, 656, 206 Pellegrini, S. 2005, ApJ, 624, 155 Pilbratt, G. L., et al. 2010, A&A, 518, L1 Poglitsch, A., et al. 2010, A&A, 518, L2 Randall, S., Nulsen, P., Forman, W. R., et al. 2008, ApJ, 688, 208 Rouleau, F., & Martin, P. G. 1991, ApJ, 377, 526 Sadler, E. M., & Gerhard, O. E. 1985, MNRAS, 214, 177 Serra DÌaz-Cano, L., & Jones, A. P. 2008, A&A, 492, 127 Soifer, B. T., Rice, W. L., Mould, J. R., et al. 1986, ApJ, 304, 651 Soria, R., Fabbiano, G., Graham, A. W., et al. 2006, ApJ, 640, 126 Takeda, H., Nulsen, P. E. J., & Fabian, A. C. 1984, MNRAS, 208, 261 Temi, P., Brighenti, F., & Mathews, W. G. 2007, ApJ, 660, 1215 Tielens, A. G. G. M., McKee, C. F., Seab, C. G., & Hollenbach, D. J. 1994, ApJ, 431, 321 Tsai, J. C., & Mathews, W. G. 1995, ApJ, 448, 84 Tsai, J. C., & Mathews, W. G. 1996, ApJ, 468, 571

Page 4 of 4