Arxiv:1505.03150V2 [Astro-Ph.GA] 1 Jun 2015 Describe These Geometries

Astronomy & Astrophysics manuscript no. NGC4370arXiv c ESO 2018 October 15, 2018 NGC4370: a case study for testing our ability to infer dust distribution and mass in nearby galaxies S. Viaene1, G. De Geyter1, M. Baes1, J. Fritz1; 2, G.J. Bendo3, M. Boquien4, A. Boselli5, S. Bianchi6, L. Cortese7, P. Côté8, J.-C. Cuillandre9, I. De Looze1; 10, S. di Serego Alighieri6, L. Ferrarese8, S. D. J. Gwyn11, T. M. Hughes1, and C. Pappalardo12 1 Sterrenkundig Observatorium, Universiteit Gent, Krijgslaan 281, B-9000 Gent, Belgium e-mail: [email protected] 2 Centro de Radioastronomía y Astrofísica, CRyA, UNAM, Campus Morelia, A.P. 3-72, C.P. 58089, Michoacán, Mexico 3 Jodrell Bank Centre for Astrophysics, School of Physics and Astronomy, University of Manchester, Oxford Road, Manchester M13 9PL, UK 4 Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge, CB3 0HA, UK 5 Laboratoire d’Astrophysique de Marseille, UMR 6110 CNRS, 38 rue F. Joliot-Curie, F-13388 Marseille, France 6 Osservatorio Astrofisico di Arcetri – INAF, Largo E. Fermi 5, 50125 Firenze, Italy 7 Centre for Astrophysics & Supercomputing, Swinburne University of Technology, Mail H30 – PO Box 218, Hawthorn, VIC 3122, Australia 8 National Research Council of Canada, Herzberg Astronomy and Astrophysics, 5071 W. Saanich Road, Victoria, BC V9E2E7, Canada 9 Canada-France-Hawaii Telescope, 65-1238 Mamalahoa Hwy, Kamuela, HI 96743 USA 10 Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge, CB3 0HA, UK 11 Canadian Astronomy Data Centre, 5071 West Saanich Rd, Victoria BC, V9E 2E7, Canada 12 Centro de Astronomia e Astrofísica da Universidade de Lisboa, Observatório Astronómico de Lisboa, Tapada da Ajuda, 1349-018 Lisbon, Portugal October 15, 2018 ABSTRACT Context. A segment of the early-type galaxy population hosts a prominent dust lane, often decoupled from its stellar body. Methods of quantifying the dust content of these systems based on optical imaging data usually yield dust masses that are an order of magnitude lower than dust masses derived from the observed far IR (FIR) emission. The discrepancy is often explained by invoking a diffuse dust component that is hard to trace in the UV or optical. Aims. High-quality optical data from the Next Generation Virgo cluster Survey (NGVS) and FIR/sub-mm observations from the Herschel Virgo Cluster Survey (HeViCS) allow us to revisit previous methods of determining the dust content in galaxies and explore new ones. NGC 4370 is an edge-on, early-type galaxy with a conspicuous dust lane and regular morphology, making it suitable for several (semi-)analytical modelling techniques. We aim to derive the dust mass from both optical and FIR data and to investigate the need to invoke a putative diffuse dust component. Methods. We used different methods to determine the total dust mass in the dust lane. We used our exquisite optical data to create colour and attenuation maps, which are converted to approximate dust mass maps based on simple dust geometries. Dust masses were also derived from SED fits to FIR to sub-mm observations. Finally, inverse radiative transfer fitting was performed to investigate more complex dust geometries, such as an exponential dust disc and a dust ring and to treat the dust-starlight interaction in a self-consistent way. 5 Results. We find that the empirical methods applied to the optical data yield lower limits of 3:4 × 10 M , an order of magnitude below the total dust masses derived from SED fitting. In contrast, radiative transfer models yield dust masses that are slightly lower, but fully consistent with the FIR-derived mass. We find that the effect of a nuclear stellar disc on the derivation of the total dust mass is minor. Conclusions. Dust is more likely to be distributed in a ring around the centre of NGC 4370 as opposed to an exponential disc or a simple foreground screen. Moreover, by using inverse radiative transfer fitting, we are able to constrain most of the parameters that arXiv:1505.03150v2 [astro-ph.GA] 1 Jun 2015 describe these geometries. The resulting dust masses are high enough to account for the dust observed at FIR/sub-mm wavelengths, so that no diffuse dust component needs to be invoked. We furthermore caution against interpreting dust masses and optical depths based on optical data alone, when using overly simplistic star-dust geometries and the neglect of scattering effects. Key words. galaxies: individual: NGC4370 - galaxies: ISM - infrared: ISM - galaxies: fundamental: parameters - dust, extinction - methods: observational 1. Introduction magnetic spectrum have revealed significant amounts of matter in the interstellar medium (ISM) and even star forming activity Early-type galaxies (ETGs) are not red and dead as was once the in these objects (see e.g. Rampazzo et al. 2005; Combes et al. common understanding. Deeper observations across the electro- 2007; Young et al. 2011). The likely mechanism to form early- Article number, page 1 of 16 type galaxies is found in major or minor mergers of spiral galax- as the dust lane (based on IRAS and ISO mass estimates, re- ies. During these events, interstellar matter of the progenitors is spectively). Such a vast reservoir of cold dust should then be partly expelled and partly consumed in waves of fast star forma- observable with the Herschel Space Observatory (Pilbratt et al. tion (see e.g. Barger et al. 1996; De Lucia et al. 2006). After a 2010). However, in most ETGs, the FIR and sub-mm emission few Gyr, the galaxy relaxes to become a stable, non-star-forming is much more compact (see e.g. Smith et al. 2012b). As a result, object with a depleted ISM. Many ETGs, especially the more dust emission is not spatially resolved in the vertical direction massive ones, are embedded in a halo of hot, X-ray-emitting gas in Herschel observations. Therefore, a direct detection of a dif- (Forman et al. 1985; Sarzi et al. 2013), which is the remnant of fuse dust component following the stellar distribution has not yet the violent merger processes. Interstellar dust cannot survive in been reported in ETGs. the harsh environment of the hot gas. Measurements of signifi- We can now revisit these problems using high-quality data cant amounts of dust of various morphologies in over half of the covering a broad range of wavelengths and advanced models. ETGs (Goudfrooij et al. 1994c,a; van Dokkum & Franx 1995; In the UV, optical and NIR, large surveys such as SDSS (York Ferrari et al. 1999; Tran et al. 2001; Patil et al. 2007; Finkelman et al. 2000), 2MASS (Skrutskie et al. 2006), GALEX (Martin et al. 2008; Kaviraj et al. 2012; Kulkarni et al. 2014) then lead et al. 2005), and WISE (Wright et al. 2010) have made multi- to the conclusion that the dust reservoir must be replenished in wavelength datasets available for many objects. In the FIR/sub- some way. Since internal production mechanisms seem to fall mm regime, Herschel proved a tremendous leap forward in terms short, a merger origin of the dust seems the most likely option of sensitivity. Consequently, the dust content in ETGs has been (see e.g. Finkelman et al. 2012, for an elaborate account). studied intensely over the past few years. Dedicated investi- The dust mass determination in ETGs is complicated by gations of individual targets (e.g. M87 Baes et al. 2010, M86 the discrepancy found between the different dust mass measure- Gomez et al. 2010, and NGC 4125 Wilson et al. 2013) were ments. Numerous FIR observations of ETGs have been obtained supplemented with statistical studies. For example, from the using IRAS, ISO, and Spitzer (see e.g. Knapp et al. 1989; Leeuw Herschel Reference Survey (HRS, Boselli et al. 2010): Smith et al. 2004; Temi et al. 2004, 2007; Xilouris et al. 2004). Temi et al.(2012b); Cortese et al.(2012b), from the Herschel Virgo et al.(2004) found that dust masses derived from ISO observa- Cluster Survey (HeViCS, Davies et al. 2010): Clemens et al. tions (2:5 − 240 µm) are more than an order of magnitude higher (2010); di Serego Alighieri et al.(2013), and from the Herschel than IRAS-derived (12 − 100 µm) dust masses. They argue that Astrophysical Terahertz Large Area Survey (H-ATLAS, Eales ISO detects much more cold dust, which could explain the dif- et al. 2010): Agius et al.(2013); Rowlands et al.(2012). As the ference. Indeed, IRAS was only sensitive to wavelengths shorter SPIRE (Griffin et al. 2010) and PACS (Poglitsch et al. 2010) in- than 100 µm, making it insensitive to the emission from cold struments on-board the Herschel spacecraft cover wavelengths (< 25 K) dust (Devereux & Young 1990). Because cold dust from 70-500 µm, the peak of cold dust emission could be suffi- is actually by far the most massive dust component in galaxies, ciently sampled, and more reliable dust masses could be derived IRAS-derived dust masses consequently yield lower limits (Gor- for all of these objects. don et al. 2010). IRAS-derived dust masses are in turn an order of magnitude On the modelling side, significant progress has been made in higher than those computed using extinction models (see e.g. the field of dust radiative transfer (for a review, see Steinacker Patil et al. 2007; Finkelman et al. 2012, and references therein). et al. 2013). Multiple efforts have already succeeded in fitting Deriving dust masses from ultraviolet (UV), optical, and near- radiative transfer models to individual galaxies (Xilouris et al. infrared (NIR) data can be difficult due to the effects of scatter- 1997, 1999; Popescu et al. 2000; Bianchi 2007; De Looze et al. ing and the complicated star-dust distribution (see e.g.

Load more