Free Access
Issue
A&A
Volume 541, May 2012
Article Number L5
Number of page(s) 4
Section Letters
DOI https://doi.org/10.1051/0004-6361/201118615
Published online 26 April 2012

© ESO, 2012

1. Introduction

The edge-on perspective of a spiral galaxy’s disk reveals both the vertical structure of the disk and its faint outskirts, thanks to line-of-sight integration. One typical feature is the mid-plane dust absorption lane. Dust is linked to the cold molecular ISM, mechanically (Allen et al. 1986; Weingartner & Draine 2001), through mutual shielding, and as a catalyst for molecular gas formation. Dalcanton et al. (2004) used the appearance of dust lanes as a probe of vertical stability of spiral disks. They found that in massive spirals, the ISM collapses into a thin dust-lane, while in low-mass disks (vrot < 120 km   s-1) the dust morphology is flocculant.

thumbnail Fig. 1

A color composite of the SPIRE maps of NGC 4244, 250 (red), 350 (green), and 500 μm (blue) with the Hi contours from Zschaechner et al. (2011) (contours at 6.4 × 1019 cm-2 increasing by factors of 2. The SPIRE flux stays within the 8.1 × 1019 cm-2 contour (second to highest).

Open with DEXTER

thumbnail Fig. 2

Major (left) and minor (right) axis profiles, averaged over a strip of 30′′ (solid lines), in the sdss-r (green) and IRAC 3.6 μm (pink, offset by a factor 2 for clarity) and the 250 (blue) and 500 (red) μm SPIRE data, all at native resolution. The major axis profile spans (RC3). We also show the smooth (dot-dash) and clumpy (dashed) TRADING models for comparison (convolved with the SPIRE PSFs at 250 and 500 μm). We assumed PA = 47.3° (de Jong et al. 2007).

Open with DEXTER

One comprehensive approach to investigate this phenomenon is to model the multi-wavelength information on a range of spiral disks with a spectral energy distribution (SED) model. Dust absorption, especially within the thin dark dust lane, is strongest at optical and ultra-violet wavelengths and the challenge for the models is to balance the observed absorption with the dust emission at far-infrared (FIR) and sub-millimeter (mm) wavelengths (e.g., Popescu et al. 2000; Misselt et al. 2001; Gordon et al. 2001; Alton et al. 2004; Bianchi 2008; Baes et al. 2010; Bianchi & Xilouris 2011). Until recently, these models were degenerate in the vertical distributions of stellar light and ISM due to a lack of resolution and wavelength coverage in the FIR and sub-mm, but with the advent of the Herschel Space Observatory, the vertical structure of nearby edge-on disks can now be fully resolved.

Massive edge-on spiral galaxies have been observed by several Herschel surveys (see Boselli et al. 2010; Eales et al. 2010b; Davies et al. 2010) and a dedicated program1 (Verstappen et al. in prep.) which targets seven massive spirals. Our NHEMESES2 program aims to observe a complementary sample of low-mass nearby edge-ons. The combined observations will serve as a test for the different radiative transfer models in the literature; those from Baes et al. (2003, 2011), Bianchi (2007, 2008), Gordon et al. (2001), Misselt et al. (2001) and Pierini et al. (2004), and Popescu et al. (2000, 2011). These models will further test the suggested dichotomy in the vertical dust structure.

Here, we explore the first NHEMESES results on NGC 4244, the prototypical low-mass (vrot = 95 km   s-1), late-type (Sc), edge-on spiral at a distance3 of 4.4 Mpc to compare to the first Herschel results on the massive edge-on spiral NGC 891 (Bianchi & Xilouris 2011). Originally, van der Kruit & Searle (1981) found a single edge-on exponential disk truncated at 9 kpc, later confirmed by Fry et al. (1999). The truncation was also observed with the HubbleSpace Telescope (de Jong et al. 2007). Yoachim & Dalcanton (2006) note a second, thick stellar disk as doubtful, but Comerón et al. (2011) report both a thin and thick disk of equal mass. Its center hosts a rotating nuclear star-cluster (Seth et al. 2008), and the disk is classified as “corrugated” (Florido et al. 1991).

NGC 4244’s ISM is mostly in a thin disk, both the ionized gas (Hoopes et al. 1999) and warm dust, implying a low specific star-formation rate (Kodaira & Yamashita 1996). The Hi disk shows a clear warp, kinematic evidence for a compact dark matter halo, and lagging extra-planar gas (Olling 1996; Zschaechner et al. 2011; Heald et al. 2011a). The H2 mass, based on CO observations, is 1.4 × 107 M (Matthews & Wood 2001).

A recent SED model by MacLachlan et al. (2011) finds that most (80%) of the FIR emission could arise from an optically thin disk with a similar scale-height for both dust and stars (as did Seth et al. 2005, based on stellar populations), a dust scale-length 1.8 times the stellar one and a dust mass of 2.38 × 106 M. They note their model diverges from the observed disk in the FIR and that the dust model is still uncertain without sub-mm observations.

In this Letter, we present new SPIRE sub-mm observations of NGC 4244, adding spatial information on the cold dust, and present two TRADING models to illustrate the role of ISM geometry on the observed overall SED and the sub-mm images.

2. Observations

SPIRE (Spectral and Photometric Imaging REceiver, Griffin et al. 2010) instrument at 250, 350 and 500 μm onboard the ESA Herschel Space Observatory (Pilbratt et al. 2010) observed NGC 4244 at 250, 350 and 500 μm in large map mode, over 25′ × 25′ centered on the object with two cross-scans, using a 30′′ s-1 scan rate. Data-reduction is as in Bianchi & Xilouris (2011); we generated maps with the naïve map-making procedure within HIPE (Ott 2010), with pixel sizes of 6′′, 8′′, 12′′ at 250, 350 and 500 μm, respectively (about 1/3 of the measured FWHM:  ≈ 18′′, 25′′ and 36′′), resulting in uniform backgrounds and sky rms noise of 1, 0.6 and 0.3 MJy sr-1 respectively, Fig. 1). A color correction half-way between that for a point- and extended-source was implemented with a color-correction in between the point- and extended-source (marginally resolved vertically), we derive integrated fluxes of 15.9, 9.8, and 5.1 Jy (250, 350 and 500 μm), and adopt a conservative 15% calibration error.

3. Results and analysis

Figure 1 shows a SPIRE color-composite, together with the Hi contours from Zschaechner et al. (2011). The dust emission falls entirely within the second highest Hi contour, concentrated in a flat, single disk with no emission in the extended warped Hi disk and envelope, in contrast to the more massive NGC 891, where Popescu & Tuffs (2003) and Bianchi & Xilouris (2011) find evidence for dust throughout the Hi disk and the extended envelope (Oosterloo et al. 2007). Dust is concentrated in the thin Hi disk, but we note that dust and Hi clumps do not coincide.

Figure 2 shows the radial and vertical profiles of the sub-mm data, sdss-r data (Aihara et al. 2011), and the Spitzer 3.6 μm image from a related Spitzer program4. The sub-mm flux is contained within the D25 radius, and the profile shows clumps of emission at both ends of the disk. These clumps of sub-mm emission are the last reliable detection of the disk before the truncation radius (420′′, 9 kpc van der Kruit & Searle 1981). The vertical sub-mm profiles are similar in width to those in Spitzer 3.6 μm emission, a good tracer of the stellar mass (e.g. Meidt et al. 2012). Figure 2 suggest that dust is confined to the thin disk, in contrast to NGC 891, where a second vertical component is observed (Kamphuis et al. 2007; Bianchi & Xilouris 2011), as well as a thin dust lane (Xilouris et al. 1999; Schechtman-Rook et al. 2012a,b).

Figure 3 shows the optical and infrared SED of the galaxy. The SPIRE flux densities show a spectral index Fν ∝ ν1.6, flatter than the more massive NGC 891 (Fν ∝ ν2.7), and thus indicative of a colder dust temperature (or of different dust properties). When compared to the radiative transfer model of MacLachlan et al. (2011), the observed fluxes at 350 and 500 μm are higher by a factor two (Fig. 3, Holwerda et al. 2011).

We model the surface brightness distribution and SED of NGC 4244 using the radiative transfer code TRADING (Bianchi 2008). We initialize the stellar distribution with the scale-length and height (h,z) as derived by Fry et al. (1999) (scaled to 4.4 Mpc), and adopted an intrinsic SED that, after being processed by dust, reproduces the observed stellar fluxes in the optical and NIR. We truncate the model, for both stars and dust, at 420′′ (9 kpc) since it is to this radius we have reliable fluxes for both.

First, we try a smooth exponential dust disk with the same scale-height (zd) as for stars, and modify the scale-length (hd) and the central V-band face-on optical depth though the whole disk () to approximate both the sub-mm surface brightness profiles (Fig. 2) and the SED (Fig. 3). As MacLachlan et al. (2011), we require a dust scale-length larger than the stellar one to match the observed major axis profiles. However, contrary to MacLachlan et al. (2011), we need a higher opacity (and hence dust mass) to match the peak of the SED5. Ours and MacLachlan et al. (2011) models are compared in Table 1. While the smooth model can broadly match the peak of the SED and the profiles, it lacks hot dust to explain emission at 70 μm. And it under-predicts the 500 μm flux, implying a lower dust temperature6.

thumbnail Fig. 3

The spectral energy distribution of NGC 4244. The Herschel/SPIRE fluxes from this work (black points) are show together with literature data (red symbols, Dale et al. 2009; MacLachlan et al. 2011). We show the smooth model (blue line); the clumpy model with and without embedded stellar emission (solid and dashed green line, respectively), and the model from MacLachlan et al. (black dotted line).

Open with DEXTER

Secondly, we consider an exclusively clumpy disk model; an exponential distribution of Giant Molecular Clouds (GMC), similar to the clumpy part of the smooth+clumpy model in Bianchi (2008) with MGMC = 106.5 M, gas-to-dust mass ratio = 100, RGMC ≃ 100 pc. Bianchi (2008) find the SED does not depend critically on the GMC size distribution. This model matches the sub-mm SED well but fails in the FIR (Fig. 3, dashed line). To match the 70 μm emission, we heat 60% of the clouds internally with a hot young starcluster (T = 35   000 K, 1.56 × 106 L, Fig. 3 solid green line), similar to the second dense dust disk with young stars introduced by Popescu et al. (2000) to reconcile FIR and sub-mm fluxes. We note that PAH fluxes (e.g., Spitzer μm) are not yet reproduced by this model but this may be solved with a range of cloud sizes and heating source types in a future iteration. The resulting heated clumpy model requires a large dust scale-length, hd ~ rtrunc, an increase in total dust mass (three times the smooth model mass) and some 60% of the GMCs to be internally heated (equivalent to ~12% Lbol of the stellar disk).

The smooth and partially heated clumpy SED models illustrate how, in the case of a small disk as NGC 4244, high-resolution sub-mm data reveals a dust disk with more cold dust than a FIR-only model, which is distributed more in clumps and throughout the stellar disk (similar scale-length and height as the stellar disk). A complete fit of the SED (including future Herschel/PACS observations) will be presented in the context of our complete low-mass edge-on sample.

Table 1

Results from the SED models.

4. Discussion

Stacked observations in the Herschel H-ATLAS survey (Eales et al. 2010a) point to specific characteristics (Td ~ 15 K, Md ~ 107 M) for blue, low-mass spirals (~109  M stellar mass Bourne et al. 2012). NGC 4244’s stellar and Hi mass are 6 and 2.5 × 109 M (Strickland et al. 2004; Zschaechner et al. 2011). Compared to the H-ATLAS results, our SED dust mass is reasonable and the gas-to-dust ratio (535–170 for the smooth and clumped model respectively) is in line with those found for smaller spiral disks (Leroy et al. 2009).

The central optical depth of is similar to the face-on optical depths in occulting galaxies (White et al. 2000; Holwerda et al. 2007b) or distant galaxy counts (Holwerda et al. 2005, 2007a,c). We therefore feel confident that our SED model accurately reflects the true dust mass of NGC 4244.

The scale-length of the dust in the SED models is consistently much larger to the stellar one; a factor 1.8 to 4.9, depending on the model. For NGC 891, a ratio of 1.5–2 is typically found (Xilouris et al. 1999; Popescu et al. 2000, 2011; Driver et al. 2007; Bianchi 2008; Bianchi & Xilouris 2011). Similarly, Holwerda et al. (2009) and Holwerda et al. (2012) report a clumpy, extended and flat dust profile in a similar mass face-on galaxy occulting a distant bright bulge (2MASX J00482185-2507365, near NGC 253). These flat radial profiles of dust suggest radial mixing of the ISM to a similar level found in the stellar disks (de Jong et al. 2007; Roškar et al. 2008, Radburn-Smith et al. 2012).

We find a dust scale-height similar to the stellar one, as do MacLachlan et al. (2011). They attribute this to either the flattening of dust disks (Dalcanton et al. 2004) or a smaller stellar scale-height due to inefficient dynamical heating of the stellar disk in lower mass galaxies (Bizyaev & Kajsin 2004). To distinguish between these two scenarios, we compare the oblateness (h/z) of the stellar and dust disk (Table 1) to those in Fig. 2 in Dalcanton et al. (2004). We also include the oblateness of the 3.6 and 8.0 μm Spitzer observations (Holwerda et al., in prep.). We note that the value for z of 300 pc is the scale-height of the thin disk found by Comerón et al. (2011) and RGB stars (Seth et al. 2005). NGC 4244’s stellar disk oblateness is at the highest, i.e., flattest value for low-mass disks, similar to the oblateness of massive disks, arguing against the less efficient disk heating compared to massive disks. The dusty disk oblateness is extremely high, even compared to the PAH emission (h/z = 10 at 8 μm), but this is predominantly because of the large scale-length of the dust disk. The similar stellar and dust scale-height and the best match by clumpy SED model are both consistent with the scenario from Dalcanton et al. (2004); in a low-mass disk such as this one, the ISM does not collapse into a thin plane but remains fractured and spread throughout the height of the stellar disk.

5. Summary and conclusions

Based on new sub-mm observations with the Herschel Space Observatory of NGC 4244, and some illustrative models with the TRADING code, we find that:

  • i.

    The cold dust is confined to the inner, flat Hi disk(Fig. 1).

  • ii.

    The dusty ISM is confined to the stellar radius (Fig. 2).

  • iii.

    A smooth model under-predicts 70 and 500 μm fluxes (Fig. 3).

  • iv.

    A clumpy dust disk model with a fraction of the clumps heated by embedded sources fits both the general SED (Fig. 3) and the morphology of the disk (Fig. 2) in sub-mm well.

  • vi.

    The inferred scales of stars and dust are consistent with a scenario that includes some radial migration of dust and a vertical ISM scale supported by turbulence (Table 1).

Future improvements of the NGC 4244 SED model includes new FIR (Herschel/PACS) information to constrain the cold dust structure and the number of GMCs and the stellar type of their heating sources. We plan to put these results in the context of the full NHEMESES and HEROES Herschel observations to infer if there is indeed a phase-change with mass in the dust structure and additionally to benchmark current SED models.

The authors thank ESA for the Fellowship Program and operation of the Herschel Space Observatory, John MacLachlan for his SED data and Laura Zschaechner for her Hi map and the anonymous referee for his or her diligent work.


1

HERschel Observations of Edge-on Spirals, HEROES.

2

New HErschel Multi-wavelength Extragalactic Survey of Edge-on Spirals (Holwerda et al. 2011).

3

The mean of the distances in NED from Heald et al. (2011b). Radburn-Smith et al. (2011) derive a distance of 4.7 Mpc.

4

P.I. R. S. de Jong, Holwerda et al., in prep.

5

Although the model of MacLachlan et al. (2011) includes a distribution of compact sources, with a dust mass about ten times than that of their diffuse disk, its contribution to the peak SED, as estimated from their Fig. 12, is smaller than 15%.

6

Our assumption is that the dust grain population of NGC 4244 and Milky Way are identical but differences in sub-mm flux and optical extinction can be due to different grain distributions as well.

References

All Tables

Table 1

Results from the SED models.

All Figures

thumbnail Fig. 1

A color composite of the SPIRE maps of NGC 4244, 250 (red), 350 (green), and 500 μm (blue) with the Hi contours from Zschaechner et al. (2011) (contours at 6.4 × 1019 cm-2 increasing by factors of 2. The SPIRE flux stays within the 8.1 × 1019 cm-2 contour (second to highest).

Open with DEXTER
In the text
thumbnail Fig. 2

Major (left) and minor (right) axis profiles, averaged over a strip of 30′′ (solid lines), in the sdss-r (green) and IRAC 3.6 μm (pink, offset by a factor 2 for clarity) and the 250 (blue) and 500 (red) μm SPIRE data, all at native resolution. The major axis profile spans (RC3). We also show the smooth (dot-dash) and clumpy (dashed) TRADING models for comparison (convolved with the SPIRE PSFs at 250 and 500 μm). We assumed PA = 47.3° (de Jong et al. 2007).

Open with DEXTER
In the text
thumbnail Fig. 3

The spectral energy distribution of NGC 4244. The Herschel/SPIRE fluxes from this work (black points) are show together with literature data (red symbols, Dale et al. 2009; MacLachlan et al. 2011). We show the smooth model (blue line); the clumpy model with and without embedded stellar emission (solid and dashed green line, respectively), and the model from MacLachlan et al. (black dotted line).

Open with DEXTER
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.