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

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

Astronomy & Astrophysics
Special feature
Letter to the Editor

Herschel: the first science highlights

FIR colours and SEDs of nearby galaxies observed with Herschel
A. Boselli1 ,L.Ciesla1 ,V. Buat1 ,L.Cortese2 , R. Auld2 , M. Baes3 ,G.J. Bendo4 ,S.Bianchi5 ,J. Bock6 , D.J.Bomans7 , M. Bradford6 , N. Castro-Rodriguez8 ,P.Chanial4 , S. Charlot9 ,M.Clemens10 , D. Clements4 , E. Corbelli5 , A. Cooray11 , D. Cormier12 , A. Dariush2 ,J. Davies2 , I. De Looze3 ,S.di SeregoAlighieri5 ,E.Dwek13 , S. Eales2 ,D.Elbaz12 , D. Fadda14 , J. Fritz3 , M. Galametz12 , F. Galliano12 ,D.A.Garcia-Appadoo15 , G. Gavazzi16 ,W.Gear2 ,C.Giovanardi5 , J. Glenn17 ,H.Gomez2 , M. Griffin2 , M. Grossi18 ,S.Hony12 , T. M. Hughes2 ,L.Hunt5 , K. Isaak2,19 , A. Jones20 , L. Levenson6 ,N.Lu6 ,S.C.Madden12 , B. O'Halloran4 , K. Okumura12 , S. Oliver21 ,M.Page22 , P. Panuzzo12 , A. Papageorgiou2 ,T. Parkin23 , I. Perez-Fournon8 , D. Pierini24 , M. Pohlen2 , N. Rangwala17 ,E.Rigby25 , H. Roussel9 , A. Rykala2 ,S.Sabatini26 , N. Sacchi26 , M. Sauvage12 , B. Schulz27 , M. Schirm23 ,M. W. L.Smith2 , L. Spinoglio26 , J. Stevens28 , S. Sundar9 , M. Symeonidis22 ,M.Trichas4 , M. Vaccari29 , J. Verstappen3 ,L.Vigroux9 , C. Vlahakis30 , C. Wilson23 , H. Wozniak31 , G. Wright27 , E. M. Xilouris32 , W. Zeilinger33 , and S. Zibetti34
(Affiliations are available in the online edition) Received 29 March 2010 / Accepted 9 April 2010
ABSTRACT

We present infrared colours (in the 25-500 m spectral range) and UV to radio continuum spectral energy distributions of a sample of 51 nearby galaxies observed with SPIRE on Herschel. The observed sample includes all morphological classes, from quiescent ellipticals to active starbursts. Active galaxies have warmer colour temperatures than normal spirals. In ellipticals hosting a radio galaxy, the far-infrared (FIR) emission is dominated by the synchrotron nuclear emission. The colour temperature of the cold dust is higher in quiescent E-S0a than in star-forming systems probably because of the different nature of their dust heating sources (evolved stellar populations, X-ray, fast electrons) and dust grain properties. In contrast to the colour temperature of the warm dust, the f 350/ f 500 index sensitive to the cold dust decreases with star formation and increases with metallicity, suggesting an overabundance of cold dust or an emissivity parameter <2 in low metallicity, active systems.
Key words. galaxies: ISM ­ galaxies: spiral ­ galaxies: elliptical and lenticular, cD ­ infrared: galaxies

1. Introduction
The energetic output of any extragalactic source can be determined by constructing its spectral energy distribution (SED). The stellar component emits in the ultraviolet (UV) to nearinfrared (NIR) domain, young and massive stars dominating the ultraviolet UV1 and old stars the NIR. Dust, produced by the aggregation of metals injected into the interstellar medium (ISM) by massive stars through stellar winds and supernovae, efficiently absorbs the stellar light, in particular that at short wavelengths, and re-emits it in the infrared domain (5 m-1 mm). At longer wavelengths, the emission of normal galaxies is generally dominated by the loss of energy of relativistic electrons accelerated in supernovae remnants (Lequeux 1971; Kennicutt 1983) (synchrotron emission). Reconstructing SEDs is thus of fundamental importance for quantifying the relative contribution of the different emitting sources to the bolometric emission of galaxies and studying the physical relations between the various galaxy components (e.g., interstellar radiation field, metallicity, dust and gas content, magnetic fields). In particular, the
Herschel is an ESA space observatory with science instruments provided by Principal Investigator consortia. It is open for proposals for observing time from the worldwide astronomical community. 1 The UV emission of early-type galaxies is however due to an old stellar population (UV upturn; e.g. O'Connell 1999; Boselli et al. 2005).

importance of the infrared domain explored by Herschel resides in the dust that, by means of the absorption and scattering of UV, optical and NIR photons, modifies the stellar spectra of galaxies. SEDs are thus crucial for quantifying dust extinction and reconstructing the intrinsic distribution of the different stellar populations within galaxies. Furthermore, fitting infrared SEDs is necessary for measuring the dust properties such as mass, temperature, fraction of PAHs, and hardness of the interstellar radiation field (ISRF), all crucial ingredients in the study of the physical processes within the ISM (e.g., Draine et al. 2007). The interpretation of the infrared SEDs of normal galaxies has already been the subject of several studies (Dale & Helou 2002; Dale et al. 2001, 2005, 2007; Chary & Elbaz 2001) even within the Virgo cluster region (Boselli et al. 1998, 2003). These however were generally limited in the infrared to 170 m, the spectral domain covered by ISO or Spitzer. These new Herschel data allow us to extend, for the first time for such a large sample of normal galaxies, to the sub-mm domain ( 500 m) where the emission is dominated by the coldest dust component. This unexplored spectral range is crucial for determining the total mass of dust and for an accurate determination of the total infrared luminosity. This paper presents the first, Herschel-based statistical study of the properties of the SED of a sample of nearby, normal galaxies spanning a large range of morphological type and luminosity. Page 1 of 5

Article published by EDP Sciences


A&A 518, L61 (2010)

2. The data
Galaxies analysed in this work were observed during the Herschel (Pilbratt et al. 2010) SPIRE (Griffinet al. 2010) science demonstration phase as part of the Herschel Reference Survey (HRS), a guaranteed time key project designed to observe with SPIRE a volume-limited, K -band-selected, complete sample of nearby galaxies (Boselli et al. 2010), and the Herschel Virgo Cluster Survey (HeViCS), an open time key project focused on covering 60 sq. deg of the Virgo cluster with PACS and SPIRE (Davies et al. 2010). To these, we added two galaxies of the Very Nearby Galaxy Sample, M 81 (Bendo et al. 2010) and M 82 (Roussel et al. 2010). The present sample is thus composed of 51 objects with photometric data in the three SPIRE bands out of which 33 are Virgo members, 13 background, 3 isolated and 2 (M 81 and M 82) nearby galaxies. The Herschel data used in this work were processed using the Level 1 procedures described in Pohlen et al. (2010), fluxes being multiplied by a factor of 1.02, 1.05, and 0.94 at 250, 350, and 500 m respectively, to take into account the updated flux calibrations. Integrated flux densities were extracted using the QPHOT task of IRAF. We assume a conservative uncertainty in the flux density of 30%. This includes uncertainties on the absolute flux calibration (15%; Swinyard et al. 2010), uncertainties introduced by the map-making technique and the possible contamination of background objects (on and off source), which might affect the flux extraction procedures. Independent observations of the galaxies NGC 4438 and NGC 4435 (the first one being an interacting system with extended morphology, the second one a point-like early-type source) performed during both the HRS and the HeViCS surveys, provided consistent results to within 10%. This should thus be considered as an upper limit to the uncertainty introduced by map making and flux extraction. Absolute flux calibration uncertainties, being systematic in the three bands, do not affect the observed trends in the SPIRE colour-colour diagrams. The dataset analysed here includes SPIRE data and measurements available at other frequencies, from UV to radio centimetre. Most of these data were taken from the GOLDMine database (Gavazzi et al. 2003). Despite this sample not being complete in any sense, and being dominated by Virgo cluster galaxies for which perturbations induced by the cluster environment may lead to systematic differences in the emission properties relative to similar isolated objects even in the FIR spectral domain (Boselli & Gavazzi 2006; Cortese et al. 2010), this is the first sample observed with Herschel that is suitable for a statistical analysis since it consists of well-known nearby galaxies spanning a wide range of both morphological type and luminosity.

Fig. 1. The infrared colours of our targets. Galaxies are coded according to their morphological type: magenta empty circles for E-S0a, red filled circles for Sa-Sb, green triangles for Sbc- Scd, blue squares for Sd, Im, BCD, and Irr galaxies. The black dotted line indicates the colour expected from the Dale & Helou (2002) empirical SED, the red longdashed line those from Chary & Elbaz (2001), the blue-short dashed, and the green dashed-dotted line the colours of the morphology- and luminosity-dependent templates of Boselli et al. (2003). Table 1. Median colours with 1 standard deviation in the colour distribution for different morphological classes.
Type E-S0a1 Sa-Sb2 Sbc-Scd Sd-Irr f 0. 0. 0. 0. 60/ f 37 ± 34 ± 43 ± 54 ± 100 0.06 0.12 0.10 0.23 f 250/ 2.60 ± 2.50 ± 2.29 ± 2.54 ± f 350 0.38 0.24 0.40 0.35 f 350/ 3.13 ± 2.86 ± 2.70 ± 2.97 ± f 500 0.39 0.22 0.44 0.35 f 0. 0. 0. 0. 25/ f 09 ± 08 ± 08 ± 24 ± 250 0.01 0.09 0.04 0.24 f 100/ 1.82± 1.10 ± 1.21 ± 1.84 ± f 250 0.73 0.42 0.45 0.63 f 60/ f 500 5.69 ± 2.96 3.14 ± 2.66 3.66 ± 2.60 10.30 ± 9.07
(2)

3. Far infrared colours
A phenomenological, model-independent technique for quantifying the spectral properties of galaxies is that of determining their colours. To do this, we combine SPIRE and IRAS flux densities, the latter being sensitive to the warm dust component. Figure 1 shows the IR colours of the sample galaxies. A first analysis of Fig. 1 indicates that in star-forming galaxies the flux density ratios f 60/ f 500, f 25/ f 250, or f 100/ f 250 are strongly correlated with the generally used IRAS colour index f 60/ f 100 (panels a-c). The dynamic range covered by f 60/ f 500, however, is a factor of about 30 larger than that covered by the f 60/ f 100 flux density ratio, and is thus a much clearer tracer of the average temperature of the bulk of the dust component. Starburst galaxies, generally defined as those
Page 2 of 5

Notes. (1) Excluding the synchrotron-dominated M 84 and M 87; cluding the starburst NGC 4491.

ex-

objects with f 60/ f 100 > 0.5 (Rowan-Robinson & Crawford 1989) have f 60/ f 500 spanning from 3 to 30. The prototype starburst galaxy in the local universe M 82 has a f 60/ f 500 of 26, significantly larger than NGC 4491, a starburst in the Virgo cluster, and VIIIZw182, a background merging system at z = 0.07. Early spirals (Sa-Sb, red filled dots, see Table 1) have f 60/ f 500 colours generally colder than Sbc-Scd (green triangles) and Sd, Im, BCD, and Irr (blue squares). Early-types with a synchrotron-dominated IR emission (M 87, M 84) are well separated in all IRAS-SPIRE or SPIRE colour diagrams with respect to the other dust-dominated E-S0a. The remaining early-types


A. Boselli et al.: FIR colours and SEDs of nearby galaxies observed with Herschel

M87 M84 M86

NGC4435 NGC4476

M58 M90 NGC4491

NGC4388 NGC3982

IC3476 IC3311 NGC4294

M100 M99

NGC4532 NGC4438 IC3583

NGC3448 VIIIZw182

Fig. 2. The UV to radio centimetre observed SEDs of a subsample of galaxies, in 4 different panels according to their morphological type. The interacting systems NGC 3448 (Arp 205) and NGC 4438 (Arp 120) are included in the Sd-BCD-Pec sample. SPIRE data are marked in red. The continuum line in the optical domain is the medium resolution integrated spectra with emission lines in star-forming systems taken from Gavazzi et al. (2004). The dashed line shows the combination of two modified black bodies ( B(, T ) with = 2) of temperature T = 20 and 40 K respectively normalized to the 100 and 60 m flux densities. The dotted line in the radio domain indicates the best-fit synchrotron emission.

have colour indices indicating that the cold dust has a higher temperature than in star-forming systems. As for the interpretation of the FIR properties of the early-type galaxies in the SINGS galaxy sample (Draine et al. 2007), higher dust-weighted mean starlight intensities can explain the high FIR colour temperatures of E/S0 galaxies. However, the relative importance of X-ray heating (Wolfire et al. 1995), stochastic heating, heating from fast electrons in the hot gas, as in supernovae, and the (unknown) size-distribution of dust grains in these environments with low-density ISM needs further exploration. The empirical SEDs of Dale & Helou (2002), Chary & Elbaz (2001), and Boselli et al. (2003), despite possible uncertainties in the absolute flux calibration (15%), only qualitatively cover the wide range of infrared colours observed in our sample even excluding the radio galaxies M 87 and M 84, and underpredict the f 250/ f 350 ratio for a given f 100/ f 250 ratio (panel d). Furthermore, these models do not reproduce the coldest colour temperatures observed in the SPIRE colour diagram f 350/ f 500 versus f 250/ f 350 (panel f). It is also interesting that even the most active galaxies such as M 82 and NGC 4491, which are expected to be dominated by warm dust heated by the dominant starburst, host a cold dust component as traced by the 500 m emission which is underestimated by models (see panel a).

4. Spectral energy distributions
Combining integrated flux densities from UV to radio centimetre we constructed the observed SED of the target galaxies. Figure 2

shows some examples of UV to radio centimetre SEDs of galaxies according to their morphological type. Figure 2 shows that in the elliptical galaxies M 87 and M 84, the sub-mm domain is dominated by synchrotron emission. M 87 is a powerful radio galaxy (Virgo A), where synchrotron dominates down to 10 m (Baes et al. 2010). M 84 is a moderately active radio galaxy with a luminosity at 20 cm of 2 â 1023 WHz-1 .Inspirals,the SPIRE data closely follow a modified black body ( = 2) of temperature T 20 K (dashed line) (e.g., Bendo et al. 2003, Bend10). This however can be taken just as a first order approximation since quantitative data in relation to dust can be determined only after an accurate SED fitting. To identify the heating sources of the emitting dust, we can use any tracer of the hardness of the interstellar radiation field. Here we adopt the birthrate parameter b, defined as the ratio of the present star formation rate (SFR) to the SFR averaged along the life t0 of the galaxy2, hence b SFRt0 /M . Following Boselli et al. (2001, 2009), SFR is proportional to the extinction-corrected UV or H flux, and M to the NIR flux. Therefore b is tightly related to the hardness of the UV radiation field. Figure 3 shows the relationship between the two colour indices f 60/ f 100 and f 350/ f 500 and the birthrate parameter, this last determined for late-type galaxies only. Figure 3 shows that the colour index f 60/ f 100, sensitive to the presence of warm dust, increases with b, indicating that galaxies with the warmest dust temperature are those at present most active in star formation (b 1). In contrast, the temperature of the cold dust appears to be anticorrelated with b, indicating
2

The birthrate parameter is also called the specific star formation rate. Page 3 of 5


A&A 518, L61 (2010)

5. Conclusions
Our analysis has enabled us to reach the following conclusions: a) the infrared colour index f 60/ f 500 is more capable of detecting a starburst than f 60/ f 100. b) Normal galaxies show a gradual increase in their dust temperature along the Hubble sequence, from Sa to Sc-Im-BCD with the exception of E-S0a, where the dust temperature is higher than in star-forming systems probably because of the different nature of their dust heating sources. c) SPIRE colours can be used to discriminate thermal from synchrotron emission in radio galaxies. d) In contrast to the warm dust, the colour temperature f 350/ f 500 index decreases with star formation activity and increases with metallicity. This admittedly weak evidence might be indicative of an overabundance of cold dust or an emissivity parameter < 2 in low metallicity, active systems.
Acknowledgements. SPIRE has been developed by a consortium of institutes led by Cardiff University (UK) and including Univ. Lethbridge (Canada); NAOC (China); CEA, LAM (France); IFSI, Univ. Padua (Italy); IAC (Spain); Stockholm Observatory (Sweden); Imperial College London, RAL, UCLMSSL, UKATC, Univ. Sussex (UK); and Caltech/JPL, IPAC, Univ. Colorado (USA). This development has been supported by national funding agencies: CSA (Canada); NAOC (China); CEA, CNES, CNRS (France); ASI (Italy); MCINN (Spain); Stockholm Observatory (Sweden); STFC (UK); and NASA (USA).

Fig. 3. The relationship between the infrared colours f 60/ f 100 and f 350/ f 500 and the birthrate parameter b (a), b)) and the metallicity (c), d)). Symbols are coded as in Fig. 1. The dashed lines give the linear best fits to the data: f 60/ f 100 = 0.18 ± 0.04 log b + 0.41 ± 0.02 (R = 0.69; = 99.99%); f 350/ f 500 = -0.23 ± 0.10 log b + 2.72 ± 0.06 (R = 0.39; = 95.16%); f 60/ f 100 = -0.44 ± 0.18 [12 + log(O/H)] + 4.22 ± 1.55 (R = 0.55; = 96.59%); f 350/ f 500 = 1.41 ± 0.55 [12 + log(O/H)] - 9.56 ± 4.70 (R = 0.55; = 97.40%), where R is the correlation coefficient and is the Spearman's probability that the two variables are correlated. Best fits with metallicity were performed excluding the 2 outliers M 82 and the perturbed system NGC 3448 (Arp 205) since their uncertain metallicity is probably not representative of normal galaxies.

References
Bendo, G., Joseph, R., Wells, M., et al. 2003, AJ, 125, 2361 Bendo, G., Dale, D., Draine, B., et al. 2006, ApJ, 652, 283 Bendo, G. J., et al. 2010, A&A, 518, L65 Boselli, A., & Gavazzi, G. 2006, PASP, 118, 517 Boselli, A., Lequeux, J., Sauvage, M., et al. 1998, A&A, 335, 53 Boselli, A., Gavazzi, G., Donas, J., & Scodeggio, M. 2001, AJ, 121, 753 Boselli, A., Gavazzi, G., & Sanvito, G. 2003, A&A, 402, 37 Boselli, A., Cortese, L., Deharveng, J. M., et al. 2005, ApJ, 629, L29 Boselli, A., Boissier, S., Cortese, L., et al. 2009, ApJ, 706, 1527 Boselli, A., Eales, S., Cortese, L., et al. 2010, PASP, 122, 261 Chary, R., & Elbaz, D. 2001, ApJ, 556, 562 Cortese, L., Davies, J. I., Pohlen, M., et al. 2010, A&A, 518, L49 Dale, D., & Helou, G. 2002, ApJ, 576, 159 Dale, D., Helou, G., Contursi, A., Silbermann, N., & Kolhatkar, S. 2001, ApJ, 549, 215 Dale, D., Bendo, G., Engelbracht, C., et al. 2005, ApJ, 633, 857 Dale, D., Gil de Paz, A., Gordon, K., et al. 2007, ApJ, 655, 863 Davies, J. I., Baes, M., Bendo, G. J., et al. 2010, A&A, 518, L48 Draine, B., & Li, A. 2007, ApJ, 657, 810 Draine, B., Dale, D., Bendo, G., et al. 2007, ApJ, 663, 866 Galametz, M., et al. 2010, A&A, 518, L55 Galliano, F., Madden, S., Jones, A., Wilson, C., & Bernard, J. 2005, A&A, 434, 867 Gavazzi, G., Boselli, A., Donati, A., Franzetti, P., & Scodeggio, M. 2003, A&A, 400, 451 Gavazzi, G., Zaccardo, A., Sanvito, G., Boselli, A., & Bonfanti, C. 2004, A&A, 417, 499 Griffin, M. J., et al. 2010, A&A, 518, L3 Kennicutt, R. 1983, A&A, 120, 219 Kewley, L., & Ellison, S. 2008, ApJ, 681, 1183 Lequeux, J. 1971, A&A, 15, L42 O'Connell, R. 1999, ARA&A, 37, 603 O'Halloran, B., et al. 2010, A&A, 518, L58 Pettini, M., & Pagel, B. 2004, MNRAS, 348, L59 Pilbratt, G. L., et al. 2010, A&A, 518, L1 Pohlen, M., et al. 2010, A&A, 518, L72 Roussel, H., et al. 2010, A&A, 518, L66 Rowan-Robinson, M., & Crawford, J. 1989, MNRAS, 238, 523 Swinyard, B. M., Ade, P., Baluteau, J. P., et al. 2010, A&A, 518, L4 Wolfire, M., Hollenbach, D., McKee, C., Tielens, A., & Bakes, E. 1995, ApJ, 443, 152

an excess of the cold dust emission in the most active galaxies. These trends with b may be non-universal since they might be related to the presence of cluster galaxies, which are characterised by a reduced star formation activity (b 0.1) because of their interaction with the cluster environment (e.g., Boselli & Gavazzi 2006). An opposite behaviour of the f 60/ f 100 and f 350/ f 500 colour indices is also present when plotted versus the gas metallicity index 12 + log(O/H) (determined using the prescriptions of Kewley & Ellison 2008; based on the Pettini & Pagel 2004; calibration and using mainly the Gavazzi et al. 2004, spectroscopic data), i.e., while f 60/ f 100 decreases with metallicity (with the exception of the starburst M 82), f 350/ f 500 seems to increase with 12 + log(O/H), with a possible exception for the interacting system NGC 3448 (Arp 205). A similar radial trend with metallicity is also observed for both M 99 and M 100 (Pohlen et al. 2010). This result is consistent with the presence of emission at > 850 m that could originate in <10 K dust (Galliano et al. 2005; Galametz et al. 2010; O'Halloran et al. 2010) or dust with < 2 (e.g. Bendo et al. 2006), which may be more prominent in low metallicity galaxies. A value of < 2 implies an enhanced contribution from carbonaceous dust because amorphous hydrocarbons have values of in the range 1-1.5. Before attempting to determine the origin of this cold dust component, this interesting result should be confirmed on a more robust statistical basis.

Page 5 is available in the electronic edition of the journal at http://www.aanda.org
Page 4 of 5


A. Boselli et al.: FIR colours and SEDs of nearby galaxies observed with Herschel
19

1

2 3 4 5 6

7 8 9 10 11 12 13 14 15 16 17 18

Laboratoire d'Astrophysique de Marseille, UMR6110 CNRS, 38 rue F. Joliot-Curie, 13388 Marseille, France e-mail: Alessandro.Boselli@oamp.fr School of Physics and Astronomy, Cardiff University, Queens Buildings The Parade, Cardiff CF24 3AA, UK Sterrenkundig Observatorium, Universiteit Gent, Krijgslaan 281 S9, 9000 Gent, Belgium Astrophysics Group, Imperial College, Blackett Laboratory, Prince Consort Road, London SW7 2AZ, UK INAF - Osservatorio Astrofisico di Arcetri, Largo Enrico Fermi 5, 50125 Firenze, Italy Jet Propulsion Laboratory, Pasadena, CA 91109, USA; Department of Astronomy, California Institute of Technology, Pasadena, CA 91125, USA Astronomical Institute, Ruhr-University Bochum, Universitaetsstr. 150, 44780 Bochum, Germany Instituto de AstrofÌsica de Canarias, C/VÌa LÀctea s/n, 38200 La Laguna, Spain Institut d'Astrophysique de Paris, UMR7095 CNRS, UniversitÈ Pierre & Marie Curie, 98bis boulevard Arago, 75014 Paris, France INAF - Osservatorio Astronomico di Padova, Vicolo dell'Osservatorio 5, 35122 Padova, Italy Department of Physics & Astronomy, University of California, Irvine, CA 92697, USA Laboratoire AIM, CEA/DSM ­ CNRS ­ UniversitÈ Paris Diderot, Irfu/Service d'Astrophysique, 91191 Gif-sur-Yvette, France Observational Cosmology Lab, Code 665, NASA Goddard Space Flight Center Greenbelt, MD 20771, USA NASA Herschel Science Center, California Institute of Technology, MS 100-22, Pasadena, CA 91125, USA ESO, Alonso de Cordova 3107, Vitacura, Santiago, Chile Universita' di Milano-Bicocca, piazza della Scienza 3, 20100 Milano, Italy Department of Astrophysical and Planetary Sciences, CASA CB-389, University of Colorado, Boulder, CO 80309, USA Centro de Astronomia e AstrofÌsica da Universidade de Lisboa, ObservatÑrio AstronÑmico de Lisboa, Tapada da Ajuda, 1349-018 Lisboa, Portugal

20 21 22 23 24 25 26 27

28

29 30 31

32 33 34

ESA Astrophysics Missions Division, ESTEC, PO Box 299, 2200 AG Noordwijk, The Netherlands Institut d'Astrophysique Spatiale (IAS), Batiment 121, Universite Paris-Sud 11 and CNRS, 91405 Orsay, France Astronomy Centre, Department of Physics and Astronomy, University of Sussex, UK Mullard Space Science Laboratory, University College London, Holmbury St Mary, Dorking, Surrey RH5 6NT, UK Dept. of Physics & Astronomy, McMaster University, Hamilton, Ontario, L8S 4M1, Canada Max-Planck-Institut fuer extraterrestrische Physik, Giessenbachstrasse, Postfach 1312, 85741 Garching, Germany School of Physics & Astronomy, University of Nottingham, University Park, Nottingham NG7 2RD, UK INAF-Istituto di Astrofisica Spaziale e Fisica Cosmica, via Fosso del Cavaliere 100, 00133 Roma, Italy Infrared Processing and Analysis Center, California Institute of Technology, Mail Code 100-22, 770 South Wilson Av, Pasadena, CA 91125, USA Centre for Astrophysics Research, Science and Technology Research Centre, University of Hertfordshire, College Lane, Herts AL10 9AB, UK University of Padova, Department of Astronomy, Vicolo Osservatorio 3, 35122 Padova, Italy Leiden Observatory, Leiden University, PO Box 9513, 2300 RA Leiden, The Netherlands Observatoire Astronomique de Strasbourg, UMR 7550 UniversitÈ de Strasbourg ­ CNRS, 11 rue de l'UniversitÈ, 67000 Strasbourg, France Institute of Astronomy and Astrophysics, National Observatory of Athens, I. Metaxa and Vas. Pavlou, P. Penteli, 15236 Athens, Greece Institut fÝr Astronomie, UniversitÄt Wien, TÝrkenschanzstr. 17, 1180 Wien, Austria Max-Planck-Institut fuer Astronomie, Koenigstuhl 17, 69117 Heidelberg, Germany

Page 5 of 5