# The Herschel Virgo Cluster Survey V. Star-forming dwarf galaxies - dust in metal-poor environments

M. Grossi, L. K. Hunt, S. Madden, C. Vlahakis, D. J. Bomans, M. Baes, G. J.

Bendo, S. Bianchi, A. Boselli, M. Clemens, et al.

► To cite this version:

M. Grossi, L. K. Hunt, S. Madden, C. Vlahakis, D. J. Bomans, et al. The Herschel Virgo Cluster Survey V. Star-forming dwarf galaxies - dust in metal-poor environments. Astronomy and Astrophysics - A&A, 2010, 518, pp.L52. 10.1051/0004-6361/201014653. hal-01439178

HAL Id: hal-01439178

https://hal.science/hal-01439178

Submitted on 23 Sep 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
LETTER TO THE EDITOR

# The Herschel Virgo Cluster Survey

## V. Star-forming dwarf galaxies – dust in metal-poor environments*

M. Grossi<fnref n="1" />, L. K. Hunt², S. Madden³, C. Vlahakis⁴, D. J. Bomans⁵, M. Baes⁶, G. J. Bendo⁷, S. Bianchi², A. Boselli⁸,
M. Clemens⁹, E. Corbelli², L. Cortese¹⁰, A. Dariush¹⁰, J. I. Davies¹⁰, I. De Looze⁶, S. di Serego Alighieri², D. Fadda¹¹,
J. Fritz⁶, D. A. Garcia-Appadoo¹², G. Gavazzi¹³, C. Giovanardi², T. M. Hughes¹⁰, A. P. Jones¹⁴, D. Pierini¹⁵,
M. Pohlen¹⁰, S. Sabatini¹⁶, M. W. L. Smith¹⁰, J. Verstappen⁶, E. M. Xilouris¹⁷, and S. Zibetti¹⁸

*(Affiliations are available in the online edition)*

Received 31 March 2010 / Accepted 2 May 2010

ABSTRACT

We present the dust properties of a small sample of Virgo cluster dwarf galaxies drawn from the science demonstration phase data set of the Herschel Virgo Cluster Survey. These galaxies have low metallicities ($7.8 < 12 + \log(O/H) < 8.3$) and star-formation rates $\le 10^{-1} M_{\odot} yr^{-1}$. We measure the spectral energy distribution (SED) from 100 to 500 $\mu$m and derive dust temperatures and dust masses. The SEDs are fitted by a cool component of temperature $T \lesssim 20$ K, implying dust masses around $10^5 M_{\odot}$ and dust-to-gas ratios $D$ within the range $10^{-3}–10^{-2}$. The completion of the full survey will yield a larger set of galaxies, which will provide more stringent constraints on the dust content of star-forming dwarf galaxies.

**Key words.** galaxies: dwarf – galaxies: ISM – dust, extinction – infrared: ISM

1. Introduction

Late-type dwarf galaxies are metal-poor systems, and they represent unique environments to investigate the properties of dust in the low-metallicity regime. Mid/far-infrared (MIR/FIR) and submillimetre (submm) observations have provided information about the different dust components in dwarfs, showing that even metal-poor galaxies may host a significant amount of dust (Thuan et al. 1999; Houck et al. 2004; Engelbracht et al. 2008). While the majority of the studies to date have focused mainly on bright and isolated dwarfs (Galliano et al. 2003; Hunt et al. 2005; Galliano et al. 2005; Galametz et al. 2009), little is known about the interplay between the environment and the dust properties of low-mass and low-metallicity systems.

It is still unclear whether most of the dust mass in these systems is at very low temperature ($T \lesssim 10$ K), as the excess emission in the submm observed in some dwarfs might suggest (Madden 2002; Galliano et al. 2005; Galametz et al. 2009). This cold dust component may either be associated with the star-forming regions, residing in clumpy molecular complexes (Galliano et al. 2003), or extend beyond the optical disc following the distribution of the neutral hydrogen (HI, Popescu et al. 2002). However, the main objection to this interpretation is the high dust mass associated with the cold dust, and alternative explanations have been proposed (see Bendo et al. 2010, and references therein). For example, Lisenfeld et al. (2002) used a mixed dust model (Desert et al. 1990) of large and small grains, finding that a high abundance of small grains (compared to the Milky Way) with a shallow emissivity ($\beta < 2$) may account for the

excess submm emission in NGC 1569. Zhu et al. (2009) applied a similar approach to the starburst spiral galaxy NGC 3310 without introducing a cold dust component. Only additional measurements in the FIR/submm regime and the analysis of a larger sample of galaxies can help us properly assess the issue of dust at low temperatures.

The Herschel Space Observatory (Pilbratt et al. 2010) allows us for the first time to study the wavelength range between 200 and 500 $\mu$m, where the emission from cold dust is expected to be the most prominent. As part of the Herschel open time key program, the Herschel Virgo Cluster Survey<fnref n="1" /> (HeViCS, Davies et al. 2010) will map an area of 64 square degrees of the Virgo cluster with PACS (Poglitsch et al. 2010) and SPIRE (Griffin et al. 2010) to investigate the dust content of the different morphological types within the cluster. The dwarf galaxy population of the Virgo cluster is dominated by dwarf ellipticals, but also contains a non-negligible fraction (~10%) of late-type dwarfs with signs of current star formation activity (Binggeli et al. 1987). In this Letter, we present the dust properties of three of these galaxies detected with PACS and SPIRE in the HeViCS science demonstration phase (SDP) data set, a $4^\circ \times 4^\circ$ field covering the central region of the cluster.

2. Sample selection

The HeViCS SDP field contains 28 late-type low-luminosity galaxies classified in the Virgo cluster catalog (VCC, Binggeli et al. 1985) as Sm (7), Im (12), blue compact dwarfs (BCDs) (6), and dwarf irregulars (3). Three of the BCDs have been detected with PACS and SPIRE and are analysed in this work: VCC 562,

<footnote>* Herschel is an ESA space observatory with science instruments provided by European-led Principal Investigator consortia and with important participation from NASA.</footnote>

<footnote n="1">¹ http://www.hevics.org</footnote>
<figure><img src="image_2.png" /></figure>

**Fig. 1.** SDSS and Herschel images of VCC 562, VCC 1179, and VCC 1356 at 100, 160, 250, 350, and 500 µm. The field size is 150". The FWHM of the beam at each wavelength is displayed at the bottom-right corner. The different shape of the PACS beams depend on the single-scan (ellipse) or crossed-link scan (circle) coverage of the area.

VCC 1179, and VCC 1356. They have *B* magnitudes between 15.5 and 16.7 ($M_B$ between −14.5 and −15.6), and HI masses of $10^7-10^8 M_\odot$ assuming a distance to Virgo of 16.5 Mpc (Mei et al. 2007). We inferred star-formation rates (SFRs) of between $10^{-1}$ and $10^{-2} M_\odot$ yr$^{-1}$ from extinction-corrected $H\alpha$ and UV fluxes (Gavazzi et al. 2004; Boselli et al. 2009). Although these dwarfs are compact, thus classified as BCDs by Binggeli et al. (1985), their SFRs are lower than the typical values for BCDs (Hunter & Elmegreen 2004; Zitrin et al. 2009). Metallicities were estimated using the Sloan Digitised Sky Survey (SDSS) spectroscopic database (York et al. 2000). Since the [OII] λ3727 nebular line is outside the observed wavelength range of the SDSS spectra, we estimated the oxygen abundances following the method of Pilyugin & Thuan (2007). Our estimates agree, within the errors, with the measurements available in the literature for VCC 562² and VCC 1179 (Vílchez & Iglesias-Páramo 2003). The basic parameters of the three galaxies are given in Table 1.

Among the Sm and Im types, only one bright IBm galaxy, IC 3583 (VCC 1686), was detected in the SDP field, but given its different properties from the three dwarfs, it will be analysed in more detail in a separate work.

### 3. Photometry

The SDP data were reduced using the Level I procedures described in Pohlen et al. (2010). More details on the observations can be found in Davies et al. (2010). The angular resolution for PACS in fast scan parallel mode is 7" × 12"/7 and 11"/6 × 15"/7, at 100, and 160 µm, respectively. For SPIRE, the PSF FWHM is 18".1, 25".2, and 36".9 at 250, 350, and 500 µm respectively. Pixel sizes are 3".2, 6".4, 6", 10", and 14" for the five bands. The rms levels of empty sky regions at 100, 160, 250, 350,

and 500 µm are ~2 mJy pix$^{-1}$, ~5 mJy pix$^{-1}$, ~8 mJy beam$^{-1}$, ~7 mJy beam$^{-1}$, and ~9 mJy beam$^{-1}$, respectively.

Figure 1 compares PACS and SPIRE data of the three dwarfs to the SDSS images. VCC 1179 is not detected at 100 µm, nor VCC 562 at 500 µm. The other galaxies are marginally detected at 500 µm with a low signal-to-noise ratio (S/N) (4σ and 5σ for VCC 1179 and VCC 1356, respectively). The VCC 1179 images show an additional feature to the north extending beyond the edge of the optical disc. It is not clear whether this feature is associated with the galaxy or not, but it also does not seem to be related to any background source in the SDSS images.

Photometry from 100 to 500 µm was derived by means of the standard growth curve analysis. The aperture is centred on the brightness peak, and the sky level is determined by considering the value that minimizes the radial variation in successively larger apertures. At large radii, the photometry should be relatively constant, and the asymptotic value corresponds to the total flux. The value of the sky thus derived is always within 1.5 standard deviations of the sky background measured in empty regions around the galaxies. The uncertainty in the flux density is assumed to be 30% (Boselli et al. 2010; Swinyard et al. 2010).

### 4. Dust temperatures and masses

To derive the temperatures of the dust, we fitted the PACS+SPIRE SEDs with a single modified Planck function and an emissivity law $k_\nu \propto \nu^2$ (Fig. 2, upper panels). We constrained the dust masses using the 250 µm flux densities, and (Li & Draine 2001). The resulting dust masses are approximately $10^5 M_\odot$, and the dust temperatures are around 20 K (see Table 1).

However, Fig. 2 shows that 500 µm fluxes of VCC 1179 and VCC 1356 tend to be underestimated by the single-temperature fits. This difference could be due to either thermal or non-thermal radio emission (Condon 1992; Hunt et al. 2005), an additional cold (~10 K) dust component (Galliano et al. 2003; Galametz et al. 2009), an enhanced abundance of small grains (Lisenfeld et al. 2002), or the different optical properties of the amorphous dust grains (Meny et al. 2007).

<footnote n="2">² Note that according to Vílchez & Iglesias-Páramo (2003) the uncertainties in the measurement of VCC 562 are also compatible with a higher metal abundance, $12 + \log(O/H) = 8.4-8.6$.</footnote>
<table><caption>**Table 1.** The properties of the three star-forming dwarf galaxies in the Virgo cluster under study.</caption><thead><tr><th>ID</th><th>m<sub>B</sub><br>mag</th><th>12 + log(O/H)</th><th>SFR<br>M<sub>☉</sub> yr<sup>-1</sup></th><th>log(M<sub>HII</sub>/M<sub>☉</sub>)</th><th>T<sub>d</sub><sup>1BB</sup><br>K</th><th>log(M<sub>d</sub>/M<sub>☉</sub>)</th><th>log(D)</th></tr></thead><tbody><tr><td>VCC 562</td><td>16.74</td><td>7.8 ± 0.2</td><td>4.6 × 10<sup>-2</sup></td><td>7.67 ± 0.03</td><td>20.2 ± 0.2</td><td>5.18 ± 0.13</td><td>-2.49 ± 0.13</td></tr><tr><td>VCC 1179</td><td>15.46</td><td>8.3 ± 0.2</td><td>9.4 × 10<sup>-2</sup></td><td>7.39 ± 0.06</td><td>16.2 ± 0.6</td><td>5.48 ± 0.14</td><td>-1.91 ± 0.15</td></tr><tr><td>VCC 1356</td><td>16.15</td><td>8.0 ± 0.3</td><td>4.6 × 10<sup>-2</sup></td><td>8.35 ± 0.01</td><td>17.9 ± 1.1</td><td>5.53 ± 0.16</td><td>-2.82 ± 0.16</td></tr></tbody></table>

**Notes.** Dust temperature, masses, and dust-to-gas ratios were obtained by fitting a single modified black-body function with $\beta = 2$. The distance to Virgo is assumed to be 16.5 Mpc.

We first considered the possibility of radio emission. Following Condon (1992), we estimated the non-thermal and thermal radio flux at 500 µm that would be expected given the SFR of each galaxy. Because of the steep frequency fall-off of the non-thermal radio component toward higher frequencies (~$v^{-0.8}$), the estimated thermal flux at 500 µm is higher than the non-thermal one by a factor of ~9. Nevertheless, in both cases, given the low SFR of these BCDs (≤10⁻¹ M☉ yr⁻¹), the expected radio flux at 500 µm is of the order of 200 µJy, and thus cannot be responsible for the 500 µm excess.

The emission from an additional cold dust component, ~10 K, would peak between 200 and 300 µm and could cause the excess we observe at 500 µm. To explore this possibility, we refined the fiducial SED model by adding a second modified black body at a lower temperature. For this two-component grey-body fit, we also used $\beta = 2$. Since we have four points in the SED of VCC 1179 and the uncertainties in the flux densities are large, the fit could not be tightly constrained, and we discuss only the results obtained for VCC 1356. More than one set of parameters were able to provide a reasonable fit to the SED of this galaxy. We inferred mean temperatures of $T_c = 8$ K and $T_w = 19.1$ K for the two dust components and the fit for these values is displayed in Fig. 2 (lower panel). The corresponding total dust mass is $M_d = 1.6 \times 10^7$ M☉. This high mass is inconsistent with current chemical evolution models (see next section). The least extreme value of cold dust temperatures providing a reasonable fit ($T_c = 11$ K, with $T_w = 19.6$ K) still yields a total dust mass $M_d = 2 \times 10^6$ M☉, which is six times higher than the value for single-temperature fit. Fitting the data with a more sophisticated dust model to test alternative explanations is beyond the scope of this Letter; this analysis will be performed in a future work when the observations of the field are completed. At this stage, given the low S/N of the detections at 500 µm and the large error bars, it is difficult to discriminate between the different scenarios.

## 5. Dust-to-gas mass ratios

The dust-to-gas-mass ratio $D(M_{dust}/M_{gas})$ not only provides information about the amount of metals that are locked in dust grains, but also provides an indication of the star-formation history of a galaxy, reflecting the net balance between the formation and destruction of dust (Hirashita et al. 2002). The ratio $D$ is known to correlate with the oxygen abundance (Lisenfeld & Ferrara 1998; Edmunds 2001; Hirashita et al. 2002; James et al. 2002). Here we check whether the gas and dust masses we obtained for these galaxies in the Virgo cluster are compatible with the predictions of dust formation models.

To derive $D$, we used neutral hydrogen masses available from the ALFALFA catalog (Giovanelli et al. 2007; Kent et al. 2008, see Table 1). The angular resolution of the 21-cm data (3.3 × 3.8) includes all dust emission detected in these objects. We assume the total gas mass is given by the atomic component

<figure><img src="image_3.png" /><figcaption>**Fig. 2.** SEDs of the VCC dwarfs compared to a single modified black body with emissivity index $\beta = 2$ (top), and a two-component model fit with $\beta = 2$, $T_w = 19.1$ K, $T_c = 8$ K for VCC 1356 (bottom).</figcaption></figure>

only: independent of the lack of CO detections, this assumption is justified by the low H₂-to-HI mass fractions expected in metal-poor environments (e.g., Robertson & Kravtsov 2008; Krumholz et al. 2009a,b).

The resulting values of $D$ are reported in Table 1. Figure 3 shows $D$ for the Virgo BCDs (filled circles) as a function of nebular oxygen abundance. Also illustrated in Fig. 3 are data from the literature (James et al. 2002; Hunt et al. 2005; Walter et al. 2007). Selected model predictions for the $D$-metallicity correlation are also shown (Edmunds 2001; Hirashita et al. 2002). The Hirashita et al. (2002) models, shown as long-dashed and dot-dashed curves, depend on the dust destruction efficiency, which begins to be effective when the oxygen abundance 12 + log (O/H) ~ 8. Were the interstellar dust mass an approximately constant fraction of the ISM metal abundance, as proposed by Edmunds (2001) and James et al. (2002), the relation between $D$ and O/H would be linear. A solid line indicates their prediction for dust production by SNe, and the short-dashed line represents that by only evolved low- and intermediate-mass stars. Alternatively, a non-linear trend is predicted by Lisenfeld & Ferrara (1998) and Hirashita et al. (2002) because of the additional effects of outflows and dust destruction efficiency. The two most metal-poor objects plotted, I Zw 18 and SBS 0335–052 (at 12 + log (O/H) ~ 7.2), have a very low dust-to-gas ratio for their metallicity (Hunt et al. 2005), but are approximately consistent with linearity.

Our estimates of $D$ for the three dwarfs, based on the single-temperature fit of the data, are indeed consistent with what would be expected by the linear models. If these galaxies hosted
<figure><img src="image_4.png" /><figcaption>**Fig. 3.** Dust-to-gas ratio *D* versus oxygen abundance; HeViCS BCDs are shown as filled circles. The vertical dotted line connects different dust-mass estimates for VCC 1356 as described in the text; the top and middle point correspond to dust masses calculated for a two-component modified black body with *T*<sub>c</sub> = 8 K and *T*<sub>c</sub> = 11 K, respectively, and the bottom point to a single-temperature fit. Data from the literature are also plotted: squares are from James et al. (2002), stars from Hunt et al. (2005), and plus signs from Walter et al. (2007). The solid and short-dashed lines correspond to the models by James et al. (2002), and the long-dashed and dot-dashed lines to models by Hirashita et al. (2002).</figcaption></figure>

a cold dust component with $T_c < 11$ K, as discussed in the previous section (and shown in Fig. 3), the corresponding *D* would be much larger than the model predictions. In this case, a significant missing gaseous component (H₂, cold HI), a few times more massive than the amount of detected HI, should be introduced to explain this discrepancy.

## 6. Comparison with far-IR properties of known star-forming dwarfs

The brightest and most metal-poor BCDs, such as IZw18 and SBS 0335-052E, have been studied with *Spitzer*, and dust has been detected even in these low abundance environments (Houck et al. 2004; Wu et al. 2007). While the SED peak of SBS 0335-052E is shifted to shorter wavelengths (λ ~ 28 µm, Houck et al. 2004), I Zw 18 has also been detected at 70 µm (Wu et al. 2007; Engelbracht et al. 2008), suggesting that a cool dust component might be present in this galaxy. Hunt et al. (2005) studied the global SEDs of seven BCDs showing that far-IR emission in these galaxies peaks at or shortward of ~60 µm.

However, the dwarfs detected with *Herschel* are more “quiescent” than the typical BCDs analysed in these studies. Walter et al. (2007) observed a sample of dwarf irregular galaxies in the M 81 group with SFRs and metallicities more similar to the Virgo dwarfs. They derived dust masses and temperatures of the M 81 dwarfs using only the 70 and 160 µm *Spitzer/MIPS* bands. Their dust-to-gas mass ratios, taking into account the total HI mass of the galaxies, are displayed in Fig. 3 (plus signs), and are about one order of magnitude smaller than those of the Virgo dwarfs with a similar metal abundance. Finally, six Im-BCD galaxies in Virgo analysed with *ISO* data by Popescu et al. (2002) show similar properties to the dwarfs in the current work.

In particular, their 60–170 µm SEDs indicate the presence of warm dust and a cooler component with a median temperature ~18 K, in agreement with our results. The resulting dust-to-gas mass ratios in a few cases are rather large, ≳0.1.

## 7. Conclusions

We have presented PACS and SPIRE observations of three star-forming dwarf galaxies in Virgo detected in the far-IR/submm regime with the science demonstration phase data set for the HeViCS survey. The data indicate the presence of cool dust with a temperature ≤20 K and dust masses around $10^5 M_\odot$. We have discussed the possibility that these galaxies host an additional cold ($T \le 10$ K) component to explain the excess at 500 µm in two of the dwarfs. However, the low *S/N* of the 500 µm detections precludes us from drawing firm conclusions. The completion of the full area of the survey will enable us to place more stringent constraints on the dust content of star-forming dwarf galaxies in a dense cluster environment.

## References

Bendo, G. J., Wilson, C. D., Warren, B. E., et al. 2010, MNRAS, 402, 1409

Binggeli, B., Sandage, A., & Tammann, G. A. 1985, AJ, 90, 1681

Binggeli, B., Tammann, G. A., & Sandage, A. 1987, AJ, 94, 251

Boselli, A., Boissier, S., Cortese, L., et al. 2009, ApJ, 706, 1527

Boselli, A., et al. 2010, A&A, 518, L61

Condon, J. J. 1992, ARA&A, 30, 575

Davies, J. I., et al. 2010, A&A, 518, L48

Desert, F., Boulanger, F., & Puget, J. L. 1990, A&A, 237, 215

Edmunds, M. G. 2001, MNRAS, 328, 223

Engelbracht, C. W., Rieke, G. H., Gordon, K. D., et al. 2008, ApJ, 678, 804

Galametz, M., Madden, S., Galliano, F., et al. 2009, A&A, 508, 645

Galliano, F., Madden, S. C., Jones, A. P., et al. 2003, A&A, 407, 159

Galliano, F., Madden, S. C., Jones, A. P., Wilson, C. D., & Bernard, J. 2005, A&A, 434, 867

Gavazzi, G., Zaccardo, A., Sanvito, G., Boselli, A., & Bonfanti, C. 2004, A&A, 417, 499

Giovanelli, R., Haynes, M. P., Kent, B. R., et al. 2007, AJ, 133, 2569

Griffin, M. J., et al. 2010, A&A, 518, L3

Hirashita, H., Tajiri, Y. Y., & Kamaya, H. 2002, A&A, 388, 439

Houck, J. R., Charmandaris, V., Brandl, B. R., et al. 2004, ApJS, 154, 211

Hunt, L., Bianchi, S., & Maiolino, R. 2005, A&A, 434, 849

Hunter, D. A., & Elmegreen, B. G. 2004, AJ, 128, 2170

James, A., Dunne, L., Eales, S., & Edmunds, M. G. 2002, MNRAS, 335, 753

Kent, B. R., Giovanelli, R., Haynes, M. P., et al. 2008, AJ, 136, 713

Krumholz, M. R., McKee, C. F., & Tumlinson, J. 2009a, ApJ, 693, 216

Krumholz, M. R., McKee, C. F., & Tumlinson, J. 2009b, ApJ, 699, 850

Li, A., & Draine, B. T. 2001, ApJ, 554, 778

Lisenfeld, U., & Ferrara, A. 1998, ApJ, 496, 145

Lisenfeld, U., Israel, F. P., Stil, J. M., & Sievers, A. 2002, A&A, 382, 860

Madden, S. C. 2002, Ap&SS, 281, 247

Mei, S., Blakeslee, J. P., Côté, P., et al. 2007, ApJ, 655, 144

Meny, C., Gromov, V., Boudet, N., et al. 2007, A&A, 468, 171

Pilbratt, G. L., et al. 2010, A&A, 518, L1

Pilyugin, L. S., & Thuan, T. X. 2007, ApJ, 669, 299

Poglitsch, A., et al. 2010, A&A, 518, L2

Pohlen, M., et al. 2010, A&A, 518, L72

Popescu, C. C., Tuffs, R. J., Völk, H. J., Pierini, D., & Madore, B. F. 2002, ApJ, 567, 221

Robertson, B. E., & Kravtsov, A. V. 2008, ApJ, 680, 1083

Swinyard, B. M., Ade, P., Baluteau, J.-P., et al. 2010, A&A, 518, L4

Thuan, T. X., Sauvage, M., & Madden, S. 1999, ApJ, 516, 783

Vilchez, J. M., & Iglesias-Páramo, J. 2003, ApJS, 145, 225

Walter, F., Cannon, J. M., Roussel, H., et al. 2007, ApJ, 661, 102

Wu, Y., Charmandaris, V., Hunt, L. K., et al. 2007, ApJ, 662, 952

York, D. G., Adelman, J., Anderson, Jr., J. E., et al. 2000, AJ, 120, 1579

Zhu, M., Papadopoulos, P. P., Xilouris, E. M., Kuno, N., & Lisenfeld, U. 2009, ApJ, 706, 941

Zitrin, A., Brosch, N., & Bilenko, B. 2009, MNRAS, 399, 924
¹ CAAUL, Observatório Astronómico de Lisboa, Universidade de Lisboa, Tapada da Ajuda, 1349-018, Lisboa, Portugal
e-mail: grossi@oal.ul.pt

² INAF-Osservatorio Astrofisico di Arcetri, Largo Enrico Fermi 5, 50125 Firenze, Italy

³ Laboratoire AIM, CEA/DSM - CNRS - Université Paris Diderot, Irfu/Service d'Astrophysique, 91191 Gif sur Yvette, France

⁴ Leiden Observatory, Leiden University, PO Box 9513, 2300 RA Leiden, The Netherlands

⁵ Astronomical Institute, Ruhr-University Bochum, Universitaetsstr. 150, 44780 Bochum, Germany

⁶ Sterrenkundig Observatorium, Universiteit Gent, Krijgslaan 281 S9, 9000 Gent, Belgium

⁷ Astrophysics Group, Imperial College London, Blackett Laboratory, Prince Consort Road, London SW7 2AZ, UK

⁸ Laboratoire d'Astrophysique de Marseille, UMR 6110 CNRS, 38 rue F. Joliot-Curie, 13388 Marseille, France

⁹ INAF-Osservatorio Astronomico di Padova, Vicolo dell'Osservatorio 5, 35122 Padova, Italy

¹⁰ 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

¹³ Universita' di Milano-Bicocca, piazza della Scienza 3, 20100, Milano, Italy

¹⁴ Institut d'Astrophysique Spatiale (IAS), Batiment 121, Universite Paris-Sud 11 and CNRS, 91405 Orsay, 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

¹⁷ Institute of Astronomy and Astrophysics, National Observatory of Athens, I. Metaxa and Vas. Pavlou, P. Penteli, 15236 Athens, Greece

¹⁸ Max-Planck-Institut fuer Astronomie, Koenigstuhl 17, 69117 Heidelberg, Germany