Free Access
Issue
A&A
Volume 622, February 2019
Article Number A89
Number of page(s) 13
Section Interstellar and circumstellar matter
DOI https://doi.org/10.1051/0004-6361/201834465
Published online 01 February 2019

© ESO 2019

1 Introduction

Galaxies are inherently three-dimensional (3D) and much of their internal structure is hidden in their projection on the sky. The presence of dust and gas in galaxies of any type further complicates the measurement of their intrinsic shape, and the properties of their stellar populations such as mass, age, metallicity, and initial mass function (IMF). Such information is already blended along the line of sight and in addition effectively obscured by dust. Interstellar dust grains both absorb and scatter light from the ultraviolet (UV) to the near-infrared (NIR) with higher efficiency for shorter wavelengths.

At the same time however, dust is a vital ingredient in galaxy evolution. As dust grains are intermixed with the gas, collisions induce thermal equilibrium between these components (Croxall et al. 2012; Hughes et al. 2015). Gas can cool by transferring kinetic energy to the dust grains, which can more easily convert this into radiative energy at far-infrared (FIR) to millimetre wavelengths. In addition, it is on the dust grains that H2 and other molecules form most efficiently, speeding up molecular cloud formation (Hollenbach & Salpeter 1971; Vidali et al. 2004; Bron et al. 2014). Finally, the already formed molecules are then shielded by the dust as harmful radiation from stars causes photodissociation (see e.g. Draine & Bertoldi 1996; Walch et al. 2015). Taken together, dust thus plays a vital role in the star formation process, and by tracking dust in different environments, it is possible to investigate the efficiency of these processes.

On average, one-third of starlight in galaxies is absorbed and re-emitted at FIR and submillimetre (submm) wavelengths (Popescu et al. 2002; Skibba et al. 2011; Viaene et al. 2016). For early-type galaxies (ETGs), this is usually less (see e.g. Bianchi et al. 2018). About 4% of ETGs show large-scale dust absorption features in optical Sloan Digital Sky Survey (SDSS) images (Kaviraj et al. 2012, 2013). A higher percentage (22%) contain molecular gas (Young et al. 2011), which is usually associated with dust absorption features in high-resolution optical imaging (e.g. Ferrarese et al. 2006). Such dust-lane ETGs occupy the phase-space between actively star-forming galaxies, and passive, massive systems. The minor-merger scenario has long been proposed as the most likely mechanism to explain the presence of dust (and gas) in their interstellar medium (ISM; see e.g. Bertola et al. 1992) In this scenario, small, gas-rich galaxies merge with passive ETGs. These events are thought to spark star formation leading to a complex ISM on short timescales (see e.g. Kaviraj et al. 2007, 2009, 2013; Davis et al. 2015). It is unclear how long the effects of a minor merger remain visible. On the one hand, the presence of dust lanes suggest thatthe ISM has had the time to settle in discs or rings. On the other hand, dust that is not shielded should not survive long in the harsh radiation field of ETGs, many of which are embedded in X-ray haloes (see e.g. Fabbiano 1989; Anderson et al. 2015) or display central far-UV emission (e.g. Burstein et al. 1988). To address these questions, itis important to study individual objects in detail, in complement with the statistical studies mentioned above.

In this paper, we follow the conventional distinction between extinction and attenuation. The combination of absorption and scattering out of the line of sight is defined as dust extinction and effectively reddens the spectrum of a galaxy. Dust attenuation is then defined as the combination of extinction and scattering back into the line of sight. Attenuation therefore depends on the global distribution of stars and dust in a galaxy. It can flatten the intrinsic reddening law as blue light is more efficiently scattered. A3D reconstruction of a galaxy and its intrinsic stellar properties thus requires an adequate correction of dust attenuation.

At the same time, the wavelength-dependent attenuation (the attenuation curve) holds information on the dust mixture itself (see Galliano et al. 2018, for a review). A steep attenuation curve can, for example, indicate an important population of small dust grains (Draine 2003; Jones et al. 2013). However, this information is encrypted by geometric effects. For example, it can also indicate that (along the line of sight) a larger portion of the stars lie behind the dust (see e.g. Viaene et al. 2017a). It is important to break this degeneracy and learn more about the intrinsic dust mix in each galaxy, which is often assumed to be the same for all galaxies.

There are three main modes of dust production known as the winds of asymptotic giant branch (AGB) stars, supernova ejecta, and ISM grain growth; these can all vary in importance for different galaxy types (see Tielens et al. 2005; Finkelman et al. 2012; Rowlands et al. 2014; Hirashita et al. 2015; De Vis et al. 2017; Zhukovska et al. 2018, and references therein). There is also a plethora of dust destruction processes (Jones 2004) depending on the local environment. It is unlikely that dust production and destruction is self-regulating for the whole galaxy population and thus differences in the intrinsic dust mix are expected from galaxy to galaxy, and even inside a single galaxy.

To approach this complex problem, solid measurements of the attenuation levels inside a galaxy and the construction of an absolute attenuation curve are required first. This has been difficult in the past, as classic extinction studies (applied in the Milky Way and Magellanic clouds by e.g. Cardelli et al. 1989; Gordon et al. 2003; Fitzpatrick & Massa 2007) are not possible when individual stars cannot be resolved. As such, measurements of the attenuation curve of galaxies is often done in a statistical way, based on broadband colours and most notably the UV slope, or through the Balmer decrement (see e.g. Calzetti et al. 2000; Wild et al. 2011; Battisti et al. 2016, and references therein). Only the dust-lane ETGs allow an easy empirical measurement of their attenuation profile due to their smooth stellar body and strong, localised absorption features. Historically, their attenuation levels can be determined by constructing a dust-free model from broadband images and comparing that to the observed, reddened galaxy (Goudfrooij et al. 1994; Patil et al. 2007; Finkelman et al. 2008, 2010; Kaviraj et al. 2012; Kulkarni et al. 2014; Viaene et al. 2015).

In addition, state-of-the-art spectral fitting can provide an independent measure of the attenuation levels when properly modelled. This requires a two-component reddening law, which differentiates between line attenuation (by dust in star-forming regions) and continuum attenuation (by diffuse dust). This method was first introduced for broadband modelling of spectral energy distributions (SEDs; Silva et al. 1998; Charlot & Fall 2000) and only later introduced in spectral modelling (see e.g. Oh et al. 2011). However, the useful comparison of spectral-fitting attenuation with empirical measurements (as outlined above) has not been made to our knowledge

Most recently, we applied this empirical method to NGC 5626, a dust-lane ETG (Viaene et al. 2017a). Instead of broadband imaging, however, we made use of integral-field spectroscopic data from MUSE (Bacon et al. 2010). This essentially provides a high-quality image of the galaxy every few Ångstrom, spanning a wide range in optical wavelength. As such, the first absolute attenuation curve at high spectral resolution could be measured for an external galaxy. In this paper, we complement this study by targeting a dust-lane ETG in the Fornax Cluster FCC 167 (NGC 1380). It is one of the few dust-lane ETGs with MUSE data available and contains only a circum-nuclear ISM, making it more comparable to the majority of the gas-rich ETG population (Young et al. 2011; Davis et al. 2013). The dust in FCC 167 could therefore be very different from the large-scale dust lane observed in NGC 5626.

The additional advantage of IFU observations is that the morphology and kinematics of the ionised gas can be studied simultaneously. The total flux of the ionised gas highlights its projected morphology, while the kinematic moment maps reveal its dynamics. This offers an insight into the 3D distribution of the ISM. The ionised gas may be diffuse or located in photodissociation regions (PDRs) around new stars. In both cases the gas can be associated with dust and thus set a prior on thegeometry of the dust component. This can in turn be used to select a limited set of possible 3D dust distributions and feed these distributions into radiative transfer simulations. Only by solving the radiative transfer equations in3D is it possible to take non-local absorption and scattering (in and out of the line of sight) into account. In this way, dust properties and dust distribution can be disentangled from each other. This was successfully applied by Viaene et al. (2015) to NGC 4370, a dust-lane ETG in the Virgo Cluster. They found that a large-scale ring of dust was able to explain the optical attenuation and was consistent with dust mass estimates from FIR emission.

The goal of this paper is to explore whether it is possible to infer qualitative information about the kind of dust in a galaxy by combining constraints on the geometry of the ISM and advanced radiative transfer modelling. This in turn is informed by an accurate and spatially resolved measurement for the attenuation curve and by the ionised-gas kinematics. While answering this question, we also hope to better understand the origin of the ISM within FCC 167. The observations and data analysed in this study are discussed in Sect. 2. In Sect. 3, we outline our measurements of the attenuation in the dust lane and we discuss other ISM tracers in Sect. 4. Radiative transfer modelling of the central regions of FCC 167 is presented in Sect. 5 before a more general discussion in Sect. 6. Finally, we formulate our conclusions in Sect. 7.

thumbnail Fig. 1

Different stages of the MGE modelling to produce attenuation maps for every wavelength slice of the MUSE datacube. Left panel: reconstructed image at 5509 Å (the V -band central wavelength). The green ellipses indicate the inner and outer edge of the dust lane (red contour). Middle panel: corresponding dust-free MGE model fit. Right panel: attenuation map at 5509 Å obtained from Eq. (1).

2 Observations and data

FCC 167 (αJ2000 = 03h36m27.55s, δJ2000 = − 34°58′33.88′′, D = 18.8 Mpc) is part of the Fornax Cluster and was observed in the context of the Fornax 3D project (F3D, PI: M. Sarzi). This is a Director’s Discretionary Time programme on the Multi Unit Spectroscopic Explorer (MUSE; Bacon et al. 2010), an integral-field spectrograph mounted on the Unit Telescope 4 of the Very Large Telescope (VLT). It is a magnitude-limited survey (mr > 17) of all galaxies within the virial radius of the Fornax Cluster. The central pointing of FCC 167, used in this paper, was observed on December 31, 2016 in five separate exposures of 720 s. Each subsequent exposure was rotated by 90° to mitigate instrumental effects. The data reduction was performed using the MUSE pipeline (version 1.6.2; Weilbacher et al. 2016). Additional cleaning of contaminating sky lines was performed using the Zurich Atmospheric Purge (ZAP; Soto et al. 2016) algorithm. For more information on the observing strategy and data reduction process, we refer to the F3D survey paper (Sarzi et al. 2018). We measure a full width half maximum (FWHM) of four MUSE pixels or 0.8 arcsec for the instrument point spread function (PSF), corresponding to a physical resolution of 70 pc.

During this investigation, we also made use of several ancillary observations of FCC 167 (see Sects. 4 and 5). To estimate the structure and content of molecular gas we relied on ALMA observations from the ALMA Fornax Cluster Survey (PI: T. Davis). These observations map the CO(1–0) transition in the galaxy. For details of the survey including moment maps we refer to Zabel et al. (2018). To constrain the dust heating by hot gas electrons and the associated dust survival timescale we investigated archival X-ray data from the Chandra X-ray Observatory.

Finally, our panchromatic radiative transfer models predict the emission in the FIR and submm wavelength regime. To estimate the energy balance in the galaxy and the dust heating mechanisms, we compared the model SED to fluxes from the Herschel space telescope (Pilbratt et al. 2010). These fluxes were measured from the PACS (Poglitsch et al. 2010) and SPIRE (Griffin et al. 2010) images provided by the DustPedia (Davies et al. 2017) database. We performed our own aperture photometry using the Python photutils package (Bradley et al. 2018) to match the MUSE field of view and adequately compare themodel and observed fluxes.

3 Spectrally resolved optical attenuation

The attenuation along a line of sight is defined as the logarithmic difference of observed ( Iλobs$I_{\lambda}^{\text{obs}}$) and intrinsic (Iλintr$I_{\lambda}^{\text{intr}}$) intensity of light Aλ=2.5log(Iλobs/Iλintr).\begin{equation*}A_{\lambda} \;{=}\; -2.5\log(I_{\lambda}^{\text{obs}}/I_{\lambda}^{\text{intr}}).\vspace*{-2pt} \end{equation*}(1)

Measuring this value is particularly difficult because the intrinsic flux is not known. However, it is possible to estimate this value by interpolating the surface brightness of dust-free regions in the galaxy, or by finding a model for the unattenuated spectrum. The former case requires a smooth stellar body, which is why this method is primarily used for ETGs (Goudfrooij et al. 1994; Patil et al. 2007; Finkelman et al. 2008, 2010; Kaviraj et al. 2012; Kulkarni et al. 2014; Viaene et al. 2015). The latter is more broadly applicable, but works best when spectra (rather than broadband fluxes) are fitted (Sarzi et al. 2006).

A zoom of the central pointing of FCC 167 is presented in Fig. 1 (left panel) and shows the symmetric and smooth stellar surface brightness distribution. The dark dust lane disturbs this distribution but is in itself also symmetric and regular, making this an ideal test case to study circumnuclear dust. We combine two ellipses as projected inner and outer boundaries of the dust lane to construct a mask. The dust lane is not visible on the other side of the major axis because of its high inclination. We do not mask these corresponding pixels as the expected attenuation on the stellar flux is minimal. We use the remaining, unmasked pixels to construct a dust-free model of the galaxy. The same technique as in Viaene et al. (2017a) is used, which we briefly summarise.

The MUSE datacube is binned by a factor of five in wavelength space to boost signal-to-noise ratio (S/N), creating independent slices of Δλ = 6.25 Å. Each slice is then treated as a separate image, with the dust lane masked (Fig. 1, left panel). A multi-Gaussian expansion (MGE; Emsellem et al. 1994; Cappellari 2002) fit is then performed to create a dust-free image from the unmasked pixels (Fig. 1 middle panel). The global attenuation in this wavelength bin is then computed by measuring the flux in the dust lane area for both the observed and dust-free image and applying Eq. (1). At the same time, attenuation maps are created for every wavelength bin by subtracting both images in log space (following Eq. (1) per pixel). An example attenuation map is shown in the right panel of Fig. 1, highlighting the enhanced attenuation in the dust lane area.

The combined global attenuation measurements (in the dust-lane area) are shown in Fig. 2 (top panel, red line). For reference, we also plot the total flux in the dust lane in the bottom panel of Fig. 2. Although our method provides an absolute measure of the attenuation in every wavelength bin, we normalise this by the 5509 Å bin, to replicate the often shown AλAV ratio. This allows us to compare the observed attenuation curve in FCC 167 with literature curves. The most direct comparison to make is the only other spectrally resolved attenuation curve of a galaxy: the Viaene et al. (2017a) curve for NGC 5626, which was derived with the same method (Fig. 2, green line). It is immediately evident that FCC 167 has a flatter curve than NGC 5626. This could point to different dust mixes or, more prosaically, to a different dust geometry. While both dust lanes can be traced back to a ring of dust (see also Sect. 5), in NGC 5626 the dust lane has a lower inclination (~ 45°) than FCC 167 (~77°). Steeper curves can be found when relatively more stars are behind the dust lane (Witt & Gordon 2000; Viaene et al. 2017a). We could thus naively expect that, at a given dust mix, FCC 167 produces a steeper attenuation curve, which is not the case. Alternatively, flatter attenuation curves could point to, on average, larger grains or a lack of smaller grains (Jones et al. 2013). These considerations thus suggest that the small grains were destroyed in FCC 167. Indeed, it is possible that the survival time of small grains is shorter in the centre of a galaxy than in large-scale dust lanes at larger galactocentric radii such as in NGC 5626.

In Fig. 2, we also highlight the wavelength of several important emission lines (shaded green) and of the NaD absorption line (shaded purple). The emission lines all show a dip in the attenuation curve. This was also noted by Viaene et al. (2017a) for NGC 5626 and attributed to spatial asymmetries in our nearly monochromatic images owing to the emission-line flux. At these wavelengths, the assumption of a symmetric and predictable brightness profile reaches its limit and the intrinsic flux behind the dust lane is underestimated. On the other hand, the NaD lines can also signal additional interstellar absorption due to cold-gas material (in which neutral Sodium can exist, see e.g. Davis et al. 2012; Sarzi et al. 2017; Nedelchev et al. 2017). This may well be the reason for the enhanced attenuation at these wavelengths in the attenuation curve of NGC 5626. Conversely, the absence of such a NaD excess in FCC 167 means that the NaD absorption lines in the MUSE spectra are purely photospheric in origin.

We can also compare these two empirically measured, spectrally resolved attenuation curves with the classic average extragalactic attenuation law from Calzetti et al. (2000). Their method is based on the modelling global measurements of the UV slope and Balmer decrement in star-forming galaxies, and was recently applied on a large set of SDSS galaxies by Battisti et al. (2016, 2017). Both attenuation curves are similar (Fig. 2, blue lines). In addition, we plot the extinction law of The Heterogeneous dust Evolution Model for Interstellar Solids (THEMIS) dust model for Milky Way (MW) dust (Ysard et al. 2015; Jones et al. 2017). It is important to note that this is in fact an extinction law and not an attenuation curve. Nevertheless, normalised to the V -band central wavelength, an extinction law mimics the effect of a thin foreground screen of dust, which is also the implicit assumption in the Calzetti method. The THEMIS curve also lies close to the Calzetti and Battisti curves in the optical (although they diverge in the UV spectrum). These literature curves seem to find the middle ground between the steeper NGC 5626 curve and the flatter FCC 167 curve. From this sparse comparison sample, it is possible to see the validity of Calzetti-like attenuation curves as an average attenuation law. However, users should be aware that large differences can occur on the level of individual galaxies.

The maps of Aλ for every wavelength bin allow us to investigate the attenuation curve on a pixel-by-pixel basis in the dust lane of FCC 167. We focus on two main properties: the slope and the strength of the attenuation. To measure this, we fit a generalised Calzetti law to the attenuation curve of each pixel as follows: Aλ=AV(k(λ)kV)(λ0.55μm)δ,\begin{equation*}A_{\lambda} \;{=}\; A_V \left(\frac{k(\lambda)}{k_V}\right) \left( \frac{\lambda}{0.55\, \mu\mathrm{m}}\right)^{\delta}, \end{equation*}(2)

where k(λ) is the Calzetti et al. (2000) curve, and AV and δ can change the shape of k(λ) by varying the strength and slope, respectively. The parameter maps for AV and δ are shown in Fig. 3.

The model AV map shows a regular morphology and closely resembles the attenuation map at 5509 Å in Fig. 1. This is expected and validates the fitting of Eq. (2) to each pixel independently. The morphology of AV also suggests a ring (rather than e.g. an exponential disc) of dust as the 3D structure behind the dust lane. Interestingly the slope map (δ, right panel in Fig. 3) has high values (flatter slopes) along the east-side edge and the lowest values (steepest slopes) towards the south and south-west edges. The middle regions of the dust lane have Calzetti-like attenuation curves (δ ≈ 0). Naively, this could be interpreted as variations in the extinction curve of the dust, and thus as variations in the dust grain size distributions. However, the change in δ is particularly smooth while dust mix variations are expected to be more localised and associated with clumpy dust structures. In addition, this assumes that the star-dust distribution along the line of sight is the same in each pixel, which is unlikely to be the case.

Viaene et al. (2017a) investigated the geometric effects of a thick dust slab embedded in a stellar body. They concluded that steep slopes (lower/negative δ) are caused by relatively more stars lying behind the dust slab along the line of sight. Vice versa, flatter slopes (higher/positive δ) can be explained by a dust structure, which lies further away and is deeper embedded in the stellar body. In this interpretation the east side of the dust lane lies further away than the west side. This is consistent with an inclined dust ring as there is also no visible dust lane on the east side of the galaxy.

In addition, the difference in δ would also imply that the north side of the dust lane must lie further away than the south side. For an axisymmetric ring of dust of uniform density, the measured AV should then be lower on the north side. This is however not what we observe. Unfortunately, within this simple framework we can not make the distinction between a possible dust density gradient from north to south, or changing dust grain properties.

thumbnail Fig. 2

Spectrally resolved attenuation curve for FCC 167 (red line, this work) and for NGC 5626 (green line, Viaene et al. 2017a). Three reference curves are plotted: The Calzetti et al. (2000) and Battisti et al. (2017) average attenuation curves for different samples of galaxies (blue solid and dashed line, respectively) and the THEMIS extinction curve (black dotted line). All curves are normalised to the V -band central wavelength; AV = 0.35 for FCC 167. For reference, the emission spectrum of the dust lane area is shown in the bottom panel. Several important emission lines (shaded green) and the NaD absorption line (shaded purple) are highlighted as well.

thumbnail Fig. 3

Parameter maps for the fitted pixel-by-pixel attenuation curves. Left panel: V -band attenuation AV. Right panel: slope parameter δ, where steeper slopes are associated with more negative δ and a Calzetti-like attenuation curve has δ = 0.

4 Interstellar medium gas properties

4.1 Ionised gas

Dust and gas are usually associated with each other in galaxies. The MUSE spectra allow us to study the ionised gas through emission lines associated with this component. The F3D analysis pipeline produces a pPXF (Cappellari & Emsellem 2004; Cappellari 2017) and GandALF (Sarzi et al. 2006) best-fit spectral model for each spatial bin. The spatial bins are determined by Voronoi binning to achieve similar signal to noise in the stellar continuum. This already provides resolved maps of the ionised gas properties.

However, to further push in spatial resolution we proceed to fit each individual MUSE spectrum with GandALF while adopting the stellar kinematics and the best-fitting stellar population model returned by pPXF and GandALF, respectively, for the Voronoi bins that contain the MUSE spaxel that is being re-fitted (see also Sarzi et al. 2018). The results are high-quality measurements of the line fluxes and velocity moments in each spaxel. The only limiting factor is the instrumental spatial resolution of 70 pc.

Our main ionised-gas measurements and results are shown in Fig. 4. Starting with the ionised-gas surface brightness (top left panel), the map for the observed flux of [NII] λ6583 line traces the smooth distribution of an inclined disc of ionised gas with a radially declining flux. The disc structure extends beyond the contours of the dust lane and a remarkable tail extends from the disc in the south-west direction.

The ionised gas reveals clear and regular rotation (Fig. 4, top second panel) and has speeds of around 250 km s−1 in the main body. A small twist is visible in the velocity minor axis towards the outer regions of the disc. This is consistent with a spiral-like shape in the velocity dispersion map (top last panel). The velocity dispersion is further enhanced in the centre of the galaxy, due mostly to beam smearing but possibly also owing to shocks or kinetic energy input associated with a central active galactic nucleus (AGN) activity.

Figure 4 also shows a preliminary axisymmetric model for the ionised-gas kinematics (top third panel). This model was constructed using the KINematic Molecular Simulation (KinMS; Davis et al. 2013) package, which can be used to simulate arbitrary gas distributions. The stellar gravitational potential is based on our best MGE model and a mass-follow-light assumption (see also Sect. 5). Whereas the model is not perfect, we can already use this to highlight the higher velocity values in the south compared to the observations. From this we conclude that the tail of ionised gas is blue-shifted with respect to the main ionised gas disc.

One of the benefits of the GandALF fitting procedure is a simultaneous fit of dust reddening in the lines, and reddening of the stellar continuum emission. The bottom left panel of Fig. 4 shows the reddening by diffuse dust (i.e. excluding Balmer decrement line reddening) as quantified by the E(BV) metric. This map bears clear morphological resemblance with the AV map in Fig. 3, and the dust lane in Fig. 1. For diffuse dust, stronger reddening coincides with higher attenuation suggesting – to first order – that the line-of-sight geometry and dust grain properties are similar throughout the dust lane.

In contrast, the distribution of the line reddening (bottom second panel in Fig. 4) does not follow a ring shape. It is more uniformly distributed across the dust lane region. Reddening in both components in the tail is marginally detected indicating only low levels of dust there. Interestingly, the line reddening is roughly six times higher than the reddening by diffuse dust across the ISM. The sources of the ionised gas are thus much more embedded than the main stellar body. We further discuss the implications of diffuse versus clumpy ISM on the observations in Sect. 6.

The last two bottom panels of Fig. 4, showing maps for the [NII]λ6583/Hα ratio and the equivalent width (EW) of Hα, convey information on the sources of ionisation at work in FCC 167. A standard BPT classification based on the [NII]/Hα versus [OIII]λ5007/Hβ diagram (Baldwin et al. 1981; Veilleux & Osterbrock 1987) and the division lines drawn and adopted by Schawinski et al. (2007) reveals that the ionised-gas emission everywhere falls in the category of low-ionization nuclear emission-line regions (LINER). The exception are the regions 3′′ away from the centre that corresponds to the local minima in the [NII]/Hα map, just inside the dust-lane regions. The nebular emission is classified as composite and is likely powered by some combination of low-level star formation and the radiation from either old stars or a central AGN emission. The EW of Ha map confirms the presence of all three such ionising sources. This indeed shows both rather constant EW values as expected in the case of LINER-like emission powered by post-AGB stars (Binette et al. 1994; Stasińska et al. 2008; Sarzi et al. 2010) and three distinct peaks that point to the additional ionising sources. The central peak is due to a true LINER AGN and the other two owing to ongoing star formation as their location corresponds to the minima in the [NII]/Hα map.

The dominant morphology in the E(BV) map corresponds remarkably well with the derived AV maps in Figs. 1 and 3. This again suggests a ring-like structure for the dust, even though this is not obvious from the ionised gas surface density. Indeed, we note that the observed nebular emission only traces gas that is ionised and not the total hydrogen reservoir. The fraction of gas that is ionised obviously depends on the density of ionising photons. In this case of FCC 167 this correlates with the main stellar body, which is centrally peaked, where ionisation is due to old stars. In addition there is the source of ionising photons from the low-luminosity AGN in the centre and from any OB stars being formed. The result is a stronger ionisation field in the centre, which explains why the [NII] line (Fig. 4, top left panel) looks more like a centrally peaked disc of material. Despite these caveats, it is still possible that the ionised gas is tracing a centrally peaked gaseous disc rather than a ring, and that the dust distribution looks more like a ring due to AGN ionisation keeping the disc in the centre too hot for dust or molecules to survive. We therefore explore the molecular gas in FCC 167 as well.

thumbnail Fig. 4

Ionised gas in FCC 167. Top left panel: surface brightness (in log scale) of the [NII] λ6583 line already highlighting a tail of ionised gas connecting to the main reservoir. Top second panel: velocity map indicating organised motion along the major axis; there are several additional kinematic features towards the outskirts. Top third panel: an axisymmetric kinematic model of the measured velocity map. Top right panel: velocity dispersion of the [NII] line showing a central peak and hints of a spiral-like shape. Bottom left panel: reddening, E(BV), by diffuse dust (as opposed to emission-line reddening) derived from the spectral fitting. Blue contours correspond to molecular gas emission measured by ALMA (see Fig. 5). Bottom second panel: reddening, E(BV), of the emission lines which is much stronger than the diffuse reddening component, and without any distinct features. Bottom third and right panels show the [NII]/Hα line ratio and Hα EW, respectively. Together these two maps indicate low level star formation inside the dust ring and nuclear ionisation due to an AGN.

thumbnail Fig. 5

Moment maps of the CO(1-0) line observations with ALMA. Left panel: total line intensity with the region of the dust lane overlaid. Middle panel: velocity map indicating a regular rotation. Right panel: linewidth map (which peaks in the centre owing to beam smearing). The size of the ALMA beam is shown in each frame as a black ellipse.

4.2 Molecular gas

The ALMA Fornax Cluster Survey (PI: T. A. Davis) mapped the molecular gas in FCC 167 with ALMA (Zabel et al. 2018). The survey targets the CO(1–0) line and the continuum emission around 3 mm for all galaxies in the Fornax Cluster with a stellar mass >3 × 108 M that contain dust (Fuller et al. 2014) or HI (Serra et al. in prep.). Figure 5 shows the CO(1–0) line observations for FCC 167.

There is strong, extended CO emission from the centre of FCC 167 which overlaps with the dust lane area (Fig. 5, left panel). It also extends towards the east, where the other side of the dust ring resides, but does not cause significant optical attenuation. The ALMA morphology is broader than the dust or ionised gas maps because of the ~ 4′′ beam. The velocity and width of the CO line also can be seen in Fig. 5 (middle and right panels). The molecular gas exhibits regular rotation, with line-of-sight rotational velocities up to 250 km s−1. Typical line widths are between 40−160 km s−1 with the highest values in the centre due to beam smearing.

We constructed a dynamical model for the molecular gas again using KinMS (Davis et al. 2013). The output of KinMS is a 3D RA-Dec-velocity data cube that can be compared directly with the ALMA data cube. We use Markov chain Monte Carlo (MCMC) methods, in particular the adamet Python package of Cappellari et al. (2013), which implements the Adaptive Metropolis algorithm of Haario et al. (2001), to constrain the morphology of the molecular gas in the centre of the galaxy. Both exponential discs and ring morphologies were fitted to the data and only a ring-shaped model results in a good fit. The best-fit results for the various morphological parameters are listed in Table 1. The inclination was fixed at 76.3°, based on our models of the dust lane (Sect. 5).

Boizelle et al. (2017) mapped the CO(2–1) transition in the centre of FCC 167, also using ALMA. Their high-resolution observations suggest a gravitationally stable and slightly warped disc of molecular gas; a brighter ring/tightly wrapped spiral structure is superimposed in the inner region. The ring morphology fitted to the CO(1–0) emission is fully consistent with that seen in CO(2–1), although the latter appears to be somewhat more centrally concentrated. This may suggest a central depression in the low density CO(1–0) material, that is perhaps carved by the central AGN (see e.g. Davis et al. 2018; Trani et al. 2018).

With increasing distance from the centre of the galaxy, it seems that the different ISM components are distributed in a nested structure of increasing spatial extent. Most of the molecular gas resides within the dust lane radius. Beyond that, the molecular gas surface density drops, suggesting that the ISM coincident with the more extended gas disc is primarily atomic and around this is a skin of ionised material.

Table 1

Best-fit model parameters for the ISM in FCC 167.

5 Three-dimensional radiative transfer models

In this section, we set out to model the 3D distribution of stars and dust in FCC 167 with the aim to reproduce the attenuation curve in the dust lane. Dust continuum radiative transfer simulations allow a proper treatment of scattering and non-local absorption. These simulations have been widely used to model broadband images of galaxies (see Steinacker et al. 2013, for a review) and study attenuation and dust heating sources in detail (De Looze et al. 2014; Viaene et al. 2017a). Solving the inverse radiative transfer problem is of particular interest for our study, as it finds the best parametric model given a set of observations. This technique was pioneered by De Geyter et al. (2013, 2014) with their code FitSKIRT, an optimiser for the radiative transfer code SKIRT (Baes et al. 2011; Camps & Baes 2015). The FitSKIRT code simultaneously fits radiative transfer models to a series of images (i.e. oligochromatic fitting), and thereby breaks the degeneracy between dust reddening and colour gradients in the underlying stellar populations. Its optimisation scheme is based on genetic algorithms, which are well suited to compare noisy models (from the Monte Carlo radiative transfer) to noisy data. This method has been successfully applied to edge-on galaxies of various types (De Geyter et al. 2014; Viaene et al. 2015; Mosenkov et al. 2018).

5.1 FitSKIRT ring model

We use FitSKIRT to find the best-fitting dust ring model as the evidence in the above sections support this geometry. We explored the possibility of an exponential dust disc, but this provided no good fits to the data. To constrain the model, we select five image slices from the MUSE data cube, at 4602, 5509, 6509, 8059, and 9340 Å. The choice of the slices is meant to span the widest range in wavelength, avoid strong emission or absorption lines, and provide a roughly equal spacing in the wavelength domain. The number of slices is a trade-off between sampling the wavelength range and computational efficiency. We found that with these five images the models converge well towards an optimal solution. Each fit is run ten times with a different random seed to assess the uncertainty in the best-fit parameters.

To obtain a parametrisation for the stellar body, we first fit a MGE model to the integrated light image of FCC 167, masking out the dust lane. This combination of 2D Gaussians can be easily deprojected to a 3D profile assuming an inclination angle, which is a free parameter. The MGE model describes the overall shape of the stellar body and the relative brightness ratios of the different Gaussian functions. The total luminosities (at each wavelength) of the stellar component are left as free parameters. This approach makes no assumptions about the underlying stellar populations and their associated colours, and limits the number of free parameters in the stellar component.

The ring of dust is described by an axisymmetric density profile in cylindrical coordinates r and z. It is defined by a Gaussian with radius R and FWHM W and a vertical scale height H, all of which are free parameters, i.e. ρd(r,z)=ρ0exp[4ln2(rR)2W2]exp(|z|H).\begin{equation*}\rho_d(r,z) \;{=}\; \rho_0\,\exp\left[ -\frac{4\ln{2}(r-R)^2}{W^2} \right] \exp\left( -\frac{|z|}{H} \right). \end{equation*}(3)

The FWHMW is trivially related to the Gaussian standard deviation σ = FWHM/(22ln2)$/(2\sqrt{2\ln{2}})$. The value ρ0 is the central density and also left as a free parameter linked to the total dust mass. Thus, there are five free geometric parameters and five free luminosity parameters. This is a modest task for FitSKIRT, which was tested on fitting more complex edge-on spiral galaxies (De Geyter et al. 2014) with up to 11 geometric parameters and only 4 luminosity constraints; the SDSS griz bands. This is comparable to the wavelength range provided by MUSE. We should therefore be able to find the best configuration, within the limits ofthe chosen geometry (a ring).

In Fig. 6 we present the model for an example image slice at the V -band central wavelength (5509 Å). The results and model quality are very similar across the spectrum. In this first model, we use the THEMIS mix for diffuse MW dust (Ysard et al. 2015; Jones et al. 2017). Visually, the model images resemble the observed images remarkably well. Monte Carlo noise and PSF convolution of the model images leads to a basic level of noise in the model maps. Adding to that the inevitable observational noise, it becomes speculative to interpret structures below the 10% residual level. Still, small deviations can be spotted in the residual maps, which are mostly associated with the densest parts of the dust ring. The model underestimates the attenuation in these regions by 10− 20%. In addition there is a discrepancy in the central zone, right above the dust lane: the model pixel values are 15− 25% lower than the observed values. This could be an overestimate of the attenuation in the inner parts of the dust ring and could be a misfit in the stellar body. In the latter case, the MGE model is not completely capable to match the brightest pixels in the core of FCC 167. This may be related to the presence of a low-luminosity AGN. Nevertheless, we consider this to be a very good radiative transfer fit as model deviations lie well below 25%.

Our second model starts from the exact same geometry, but relies on a stripped version of the dust mix. To assert maximum contrast, weremoved all carbonaceous material from the THEMIS MW model. There are no amorphous hydrocarbons as small grains (the most significant alteration), and the large silicate grains are stripped from their carbon-rich mantles, leaving only the bare silicate grains. The average grain size of this dust mix is larger, which further affects the attenuation and emission properties. The best-fit model for this dust mix is shown in the bottom row of Fig. 6. Visually, this model looks almost identical to the first and matches the observations remarkably well. The residual maps are slightly better with most pixels below 20% and only a few in the 20−25% range. The largest deviations can again be associated with the nuclear regions, where the MGE component seems to fall short to explain the sharp rise in brightness towards the central pixels.

The best-fit parameters and their uncertainties are summarised in Table 1. The uncertainties were derived from running the fit ten times and computing the standard deviation on the spread in each parameter. The two models produce marginally consistent dust masses of (5.01 ± 0.34) × 104M vs. (6.30 ± 0.60) × 104M; the silicate model has a higher mass. Both models also retrieve the same inclination angle (77.6° ± 1.8° and 76.3°± 1.2° for the MW and silicate models, respectively). The inclination angle is perfectly consistent with the value of 77.6° obtained from the stellar dynamical modelling presented in Sarzi et al. (2018).

The structural parameters for the dust ring are also consistent between the two dust mix models, butthe uncertainties are larger. On the one hand, the MW dust model finds a radius of 321 ± 40 pc, a FWHM of 276 ± 99 pc, and a scale height of 16.6 ± 6.0 pc. On the other hand, a smaller (R = 285 ± 53 pc) and more narrow (W = 217 ± 151 pc) ring is found for the silicate model and has a slightly lower scale height (H = 13.6 ± 8.2 pc). While the relative uncertainties on the radius are reasonable, they are much higher for the ring width and scale height, especially for the silicate model (up to 70%). The uncertain values for the scale height may be linked to the resolution limit of the data. The best-fitting scale heights are comparable with the pixel scale of 17.8 pc and well below the PSF FWHM of 70 pc. The dust scale heights are an order of magnitude smaller than the best-fit model for the molecular gas. This can mostly be attributed to the difference in beam size. The molecular scale height corresponds to 0.6 times the 4′′ ALMA beam.

There may also be some weak degeneracies in the fitting; a dust lane may be represented by a wide ring and smaller radius or by a larger radius and more narrow ring. A key diagnostic to break this degeneracy is an accurate representation of the dust attenuation in the central regions and of the peak attenuation in the ring. Both aspects are influenced by the underlying stellar distribution. As we noted before, the MGE model is unable to perfectly retrieve the central peak and the global brightness level at the same time. As such, FitSKIRT can only retrieve the structural parameters of the dust ring up to a certain level of accuracy. Still, the level of accuracy seems reasonable given the complexity of the problem.

An important comparison to differentiate the two FitSKIRT models is their optical attenuation curve. One of the output products of the fitting is a dust-free image at every wavelength bin. We measured the flux in the dust lane area for the model image and the dust-free counterpart and applied Eq. (1) to obtain the attenuation curve. This can directly be compared to the empirically derived global attenuation curve from Sect. 3. The attenuation curves for both dust mixes are shown in Fig. 7, along with their model uncertainties. The model with the MW mix overestimates the observed attenuation level by 0.05−0.08 mag in the blue part of the spectrum, but only by 0−0.02 mag in the red part. The general slope of this model curve is significantly steeper than the observed slope. The comparison is much better for the bare silicate model; the offsets are consistently between 0 and 0.03 mag across the wavelength range. The best-fit attenuation curve for this model perfectly matches the observed slope. This is perhaps the clearest difference between the two radiative transfer models for FCC 167. While MW dust results in an attenuation curve that is clearly too steep, the bare silicate dust mix naturally produces the observed slope.

For completeness, we also plot the corresponding intrinsic extinction curves κ for each dust mix in Fig. 7. We normalise these geometry-free curves to the model attenuation at 5509 Å for ease of comparison. This shows clearly the effect of the ring geometry for both models, converting a steeper extinction curve onto a flatter attenuation curve. Interestingly, the slope difference in extinction is significant, but not that large. This is then amplified by the geometry into a large slope difference in attenuation.

thumbnail Fig. 6

Best-fit models for a ring of dust with a THEMIS MW dust mix (top row) and for a THEMIS bare sillicate dust mix (bottom row). We only show a representative image slice at the V -band central wavelength (5509 Å) as the results are very similar across the spectrum. Both observed images in the first column are the same. The second column shows the corresponding model images and the third column shows the residual maps following IobsImodIobs$\frac{I_{\text{obs}}-I_{\text{mod}}}{I_{\text{obs}}}$ in pixel flux.

thumbnail Fig. 7

Attenuation curves of the best-fit radiative transfer models for a dust ring geometry. The THEMIS MW dust mix produces a steeper curve (blue line) than a dust mix with only bare silicate grains (orange line). The dashed lines correspond to the corresponding intrinsic extinction curves of the dust mixes, normalised to the model attenuation at 5509 Å. The observed attenuation curve of FCC 167 (from Sect. 3) is shown in black.

5.2 Panchromatic model and energy balance

The centre of a massive galaxy like FCC 167 could be a hostile environment for dust grains. The high stellar density results in a strong interstellar radiation field and the Hα EW map shows further ionising radiation from low-level star formation and a LINER AGN. Both stellar sources contribute to the heating of the dust grains. In addition, there may be other heating mechanisms at work, such as through cosmic rays or collisions with hot gas atoms. To assess the importance of these contributions, we measured the energy deficit in FCC 167. In a scenario where the dust is only heated by the evolved stellar body, there should be a perfect match between the observed FIR/submm flux and the predicted flux from the FitSKIRT model. To simulate the panchromatic emission for each model, we adopted the fitted parameters for the dust component and assign an emission template to the stellar component. To estimate the age and metallicity of the stellar component we used pPXF to fit the MUSE global spectrum extracted in the same central region used for our radiative-transfer modelling. Using the Bruzual & Charlot (2003) SED templates, this fit converged to 10-Gyr-old, solar metallicity average stellar population in the central region of FCC 167. This template was normalised to the total 2MASS K-band flux, assuming this band only contains flux from the evolved stellar populations with no significant dust attenuation. The transfer of radiation is then simulated in 3D with SKIRT.

Figure 8 compares the panchromatic SED models with the observed fluxes. The fluxes were measured from the DustPedia maps (see Sect. 2), and the optical spectrum was extracted from the MUSE data cube. Both models perform well – by design – in the optical and NIR. They also perform reasonably well in the UV with a slight overprediction of the far-UV radiation. The discrepancies are much larger in the infrared, where neither of the model SEDs reach the observed fluxes. In the mid-infrared (MIR), the bare silicate model lacks any features of transiently heated grains because there are no small grains in the dust mix. The MIR is therefore significantly lower than the MW dust model. In the FIR, both models peak at the same wavelengths, and the silicate model exhibits a slightly steeper Rayleigh-Jeans tail. In fact, the Rayleigh-Jeans slope of the silicate model lies closer to the observed slope than the MW model. But in both cases the dust emission is up to an order-of-magnitude lower than observed, and the dust peaks at shorter wavelengths, indicating a higher temperature.

We calculate the energy deficit by integrating the emission from 12 to 500 μm equivalent to the WISE 3 – SPIRE 500 μm bands. The MW model only produces 35% of the observed flux, and the silicate model produces 33%. Another way to quantify the energy deficit is to infer the total dust mass from the observed SED. We fitted a modified black-body function to the Herschel fluxes; a variable β matches the steep Rayleigh-Jeans slope. We find that the opacity to convert flux to dust mass is κ350 = 0.348 m2 kg−1 for both dust models. Both dust mixes have a similar value because at 350 μm (the commonly used wavelength to convert flux to dust mass) the emission is dominated by large silicate grains. Within the uncertainties, both dust models thus give the same dust mass from the observed fluxes, which is 8.31.8+2.2×105M$8.3^{+2.2}_{-1.8} \times 10^5 M_{\odot}$. This is already more than an order of magnitude higher than the dust masses derived from visual attenuation (Table 1). We discuss this further in the next section.

thumbnail Fig. 8

Panchromatic radiative transfer model predictions based on the best-fit FitSKIRT models for THEMIS MW dust and the bare silicate dust mix. A single (10 Gyr) stellar population was used for the intrinsic SED of the stellar component. No additional dust heating sources are adopted in the model. The green line represents a modified black-body fit to the Herschel fluxes from 100 to 500 μm; a free β matches the steep Rayleigh-Jeans slope.

6 Discussion

6.1 Dust distribution and survival

In Sects. 3 and 5 we measured the global optical attenuation curve in FCC 167 and found a radiative transfer model that replicates this attenuation curve. The best match is obtained with a dust mix that contains no small grains. However, the models predict dust emission spectra which lie an order of magnitude below the observed fluxes at long wavelengths. The energy balance problem is widely known in radiative transfer models (see Mosenkov et al. 2018, for the most recent discussion on this topic) and varies between a factor of 1.5 and 5 for spiral galaxies. Little research has been carried out to assess the energy balance in ETGs, but the discrepancy in FCC 167 (a factor of ~ 3) is comparable with these studies for spiral galaxies. Several factors seem to be at play.

First, there could be a diffuse dust component that is not associated with the dust lane, but with the general stellar body. This was originally postulated by Goudfrooij & de Jong (1995) to explain similar order-of-magnitude differences in attenuation-based and emission-based dust masses for ETGs. This theory is supported by Kaviraj et al. (2012) as they find similar massdeficits for a large sample of dust-lane ETGs. However, the most reasonable origin for such a diffuse dust component is mass loss from AGB stars, but mass-loss yields for these stars fall drastically short of explaining the observed dust masses (Finkelman et al. 2012; Rowlands et al. 2014). Herschel observations of FCC 167, while coarse in resolution, show a clear, unresolved FIR/submm peak at the centre of the galaxy, and no extended component. Both the ionised gas and molecular gas are confined to the centre of the galaxy as well. In the optical spectra of FCC 167 we find only weak reddening for the stellar continuum, and the vast majority of it happens inside the dust lane.

However, evolved stars could also host warm dust in their atmospheres as opposed to dust expelled by winds. Such a dust component would be brightest in the MIR (see e.g. Temi et al. 2009). Using Eq. (1) from Davis et al. (2014) we compute a 3− 27% contribution of this dust to the WISE 22 μm flux. The radiative transfer models do not include such warm dust around evolved stars. This could in part explain the 45% deficit between the observed WISE 22 μm flux and the MW dust model flux, and the 78% deficit for the bare silicate model. Still, the discrepancy is too large to be explained solely by this component.

Second, the dust may be clumpy, which is an explanation that has also been put forward to rectify the energy balance problem in spiral galaxies (Bianchi 2008; Baes et al. 2010; De Looze et al. 2012). Cold and dense clumps would contribute little to the attenuation of starlight since they are compact. These clumps have an equilibrium temperature and thus emit FIR/submm light as a modified black body. The difference in dust mass suggests that roughly 7.7 × 105M or 93% of the observed dust mass should be locked up in clumps. This percentage seems to be particularly high, but it could be a result of intense dust destruction (see below), where only dense clumps can survive. A clumpy medium is also favoured by a much higher line reddening than continuum reddening (Fig. 4, bottom left panels) and by evidence of opaque molecular clouds in this galaxy (see discussion in Sect. 4.1 of Boizelle et al. 2017). Higher line-to-diffuse reddening ratios are seen in other galaxies, especially for those with higher mass (e.g. Koyama et al. 2019). If line emission originates from compact ionised-gas clouds, thenthese seem to be much dustier than the diffuse ISM of FCC 167. An additional effect of clumpy dust is a shift of the FIR peak to longer wavelengths as the dust is generally colder. This would bring the current panchromatic models closer to the observed SED, however, we lack the spatial resolution to test this properly in our simulations.

Third, it is likely that the dust is not only heated by the evolved stellar populations. The BPT line ratios reveals low level star formation just inside the dust lane, although the UV and absorption-line signature of young stellar populations may have been lost in the global spectra and broadband measurements that we extracted in the central regions of FCC 167. Additionally, the large nebular-line reddening suggest that such star-forming regions may be heavily obscured. Ongoing star formation can heat the dust in two ways: through conventional photon heating as young stellar populations emit most of their light as energetic UV and optical radiation and through generation of cosmic rays, which can get caught by dust grains, dissipating energy as heat.

Finally, it is also not uncommon for high-mass galaxies to contain a hot gas reservoir. When a hot gas electron hits a dust grain, it dissipates away its energy onto the grain, thus heating the dust. An important additional effect is that these electrons can damage the surface of the dust grains. Such “sputtering” by hot gas is an effective dust destruction mechanism and is particularly effective in destroying the small grains (Jones 2004). FCC 167 was detected in one of the fields observed by the Chandra X-ray Observatory. The X-ray emission is clearly separated from the main intracluster X-ray emission, which is centred at NGC 1399 (Scharf et al. 2005). In order to assess and separate the relative importance of X-ray binaries and diffuse hot gas to such X-ray emission, we extracted the X-ray spectrum in the nuclear region of the galaxy as shown in Fig. 9.

The XSPEC (Arnaud 1996) package was used to fit a two-component (wa*po+mekal) model to the observed spectrum. We used an absorbed power-law (wa*po) to account for X-ray photons from unresolved binaries and from the (potential) AGN. To describe the hot gas we use a single-temperature mekal model, which includes line emissions from several elements. The model in Fig. 9 shows that a hot gas component is a necessary addition to the power-law continuum. Our fit yields a hot gas density of 7 × 10−3 cm−3 with a temperature of 3.2 × 107 K. Dust heating by hot gas could thus be an important factor in FCC 167.

These hot gas properties can also be used to estimate the dust lifetime following Jones (2004). Dust grains, fully exposed to the hot gas in this galaxy would be destroyed on a timescale of 0.4 Myr (very small grains) to 14 Myr (large grains). At a characteristic galactocentric radius of 285 pc and rotational velocity of 250 km s−1, the orbital timescale is ~7 Myr, which is comparable to the dust destruction timescale. The efficient dust destruction in the centre of FCC 167 raises the important question of why there is a dust lane in the first place. If the dust lane is stable and long-lived, then it must be self-shielding or partly shielded by the gas. This supports the suggestion that a good fraction of the dust is organised in dense clumps. Ongoing star formation inside the dust lane also produces some dust and replenishes the current reservoir.

thumbnail Fig. 9

Fit tothe Chandra X-ray spectrum of FCC 167. Top panel: observed photon counts per channel (crosses). The black line shows the total model flux. Bottom panel: unfolded model (solid line) and its separate components (dotted line). A wa*po+mekal model was used to represent a power law and a hot gas line spectrum. The unfolded model is adjusted for spectral resolution and filter response function before comparison to the data in the top panel.

6.2 Origin of the ISM

The results acquired in this study allow us to speculate about the origin of the gas and dust in FCC 167. Statistically, dust lanes in ETGs areattributed to gas-rich minor mergers (e.g. Kaviraj et al. 2012; van de Voort et al. 2018). A relic of such a scenario would be the blue-shifted tail of ionised gas falling onto the main reservoir from behind the galaxy. This means we are observing a peculiar moment in the merger history of FCC 167, right after a minor merger. If gas and dust were externally accreted, it would be heavily processed by the hot gas on its way inward. This could explain the lack of small grains we derive from the radiative transfer fit. Once the material settles in a dust ring, it can cool down and start forming stars in the shielded inner regions. In the very centre however, fuelling of the supermassive black hole induces feedback and again dissociates the molecular gas. However, such a minor merger is disfavoured by the fact the galaxy lives in a cluster. Mergers and accretion events are relatively rare in such an environment compared to the field.

Alternatively, we must consider that the gas and dust was already present in the galaxy, and is now being stripped. The gas-to-dust ratio of ~80 in FCC 167 is more in line with discs of spiral galaxies (see Fig. 2 in Davis et al. 2015, for a comparison) and ram-pressure stripping is a common process in a clusters like Fornax. It is however unlikely that FCC 167 is radially falling into the cluster centre as the tail of ionised gas is pointing towards the central galaxy NGC 1399. To better assess this, we must look at the problem in 3D. Distances based on surface brightness fluctuations measured by Ferrarese et al. (2000) indicate that FCC 167 (18.8 Mpc) is closer to us than NGC 1399 (20.5 Mpc). However, FCC 167 is moving away from us much faster (1877 km s−1) than NGC 1399 (1425 km−1). Thus, based on the radial velocity component, FCC 167 seems to move towards NGC 1399. We can infer the tangential velocity component of the galaxy from the tail of ionised gas, which points south towards NGC 1399 and is blue-shifted. The tangential component must point north, so gas is pushed out on the south side and towards us. This would mean that FCC 167 is in fact moving through the cluster around NGC 1399 and is on its way to apocentre. As a consequence, the gas and dust in FCC 167 may be a relic of a much larger reservoir, pre-dating its interaction with the cluster. In the above stripping scenario, the dust gradually gets exposed to the harsh radiation field and can transform into a clumpy, self-shielding medium with little small grains.

7 Summary and conclusions

The goal of this paper was to address the origin and fate of the nuclear dust lane in FCC 167. To this end we have studied several aspects of the ISM in this ETG. Our investigation was centred around new MUSE observations from the Fornax 3D project. We first measured the attenuation curve in the dust lane area of the galaxy by treating each slice in the MUSE cube as a separate image. A dust-free surface brightness model was compared to the observations to determine the dust attenuation.

– We find that the attenuation curve in FCC 167 is quite flat compared to literature curves. Given the high inclination of the galaxy, this is already hard to explain with a canonical MW dust grain mix and points towards a larger average grain size.

Pixel-to-pixel variations in the attenuation curve strength and slope seem consistent with an inclined ring of dust. To better assess the ISM distribution we also map the ionised and molecular gas in FCC 167. The molecular gas was detected by ALMA in the CO(1-0) line and is more compact than the dust lane. It shows regular rotation and a centrally peaked velocity dispersion. A dynamical model suggests again a ring-like morphology, potentially because nuclear feedback is preventing molecules to form in the very centre of FCC 167. The ionised gas is organised in a regularly rotating disc that extends slightly beyond the dust lane.

– The line reddening is roughly six times stronger than reddening of the stellar continuum. Reddening is limited to the area of the dust lane and again consistent with a ring-like, perhaps clumpy dust distribution. We detect low-level star formation in the inner parts of the gas disc and find traces of AGN activity through increased ionisation rates.

The ISM components are distributed in a nested structure of increasing spatial extent. The molecular gas resides closest to the centre of the galaxy, within the radius of the dust lane. Beyond that we find the ionised gas, with an interesting tail towards the south side of the galaxy and slightly blue-shifted with respect to the main ionised disc. The position and relative velocity of FCC 167 to the cluster centre indicate that the galaxy is moving through the cluster on its way to apocentre.

– Two very different scenarios emerge to explain this nested ISM structure. On the one hand, gas could be falling into the central potential of FCC 167. The ionised gas gradually settles in a ring where it can cool down and form molecules and stars in the shielded inner regions However, this scenario is unlikely to happen in a cluster environment. On the other hand, the current ISM could be a relic of a far greater reservoir that predates the interaction of the galaxy with the cluster environment. The gas and dust can gradually be stripped by ram pressure, exposing the dust to sputtering and ionising the gas.

We performed 3D radiative transfer simulations to fit the intrinsic structure of the dust and break the geometry-dust mix degeneracy that shapes the attenuation curve.

– A model without small carbonaceous grains fits the optical attenuation much better than the Milky Way dust mix. Still, both models underestimate the FIR dust emission by ~ 67%. We attribute the energy deficit to a clumpy dust distribution after ruling out a second, diffuse, dust component.

Our interpretation is supported by the high line-to-continuum reddening ratio and by efficient dust destruction timescales. Indeed, hot gas is present in the central regions of this galaxy. A fit to the X-ray spectrum yield gas properties that would destroy the dust grains by sputtering of the order of an orbital timescale (~7 Myr).

– Efficient dust destruction explains why almost no small grains are present in the dust lane and requires the remaining dust to be self-shielding if the ISM is to survive.

Our study has revealed that we may be witnessing a very particular moment in the history of FCC 167, where dust and gas are quickly processed as they are exposed to a hostile environment. Our findings show how detailed analyses of individual systems can complement statistical studies of dust-lane ETGs. Together, these studies can pinpoint the origin of the ISM and place dust-lane ETGs on the general timeline of galaxy formation.

Acknowledgements

S.V. is supported by the UGent-BOF fund and the UHerts visitor fund. E.M.C. and L.M. acknowledge financial support from Padua University through grants DOR1699945/16, DOR1715817/17, DOR1885254/18, and BIRD164402/16. GvdV acknowledges funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme under grant agreement No 724857 (Consolidator Grant ArcheoDyn). We thank Anthony Jones for several stimulating discussions. Based on observations collected at the European Organisation for Astronomical Research in the Southern Hemisphere under ESO programme 296.B-5054(A). This research has made use of data obtained from the Chandra Data Archive. This paper makes use of the following ALMA data: ADS/JAO.ALMA#2015. 1.00497.S. ALMA is a partnership of ESO (representing its member states), NSF (USA) and NINS (Japan), together with NRC (Canada) and NSC and ASIAA (Taiwan) and KASI (Republic of Korea), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO, and NAOJ.

References

  1. Anderson, M. E., Gaspari, M., White, S. D. M., Wang, W., & Dai, X. 2015, MNRAS, 449, 3806 [NASA ADS] [CrossRef] [Google Scholar]
  2. Arnaud, K. A. 1996, in Astronomical Data Analysis Software and Systems V, ASP Conf. Ser., 101, 17 [NASA ADS] [Google Scholar]
  3. Bacon, R., Accardo, M., Adjali, L., et al. 2010, in Ground-Based and Airborne Instrumentation for Astronomy III, Proc. SPIE, 7735, 773508 [Google Scholar]
  4. Baes, M., Fritz, J., Gadotti, D. A., et al. 2010, A&A, 518, L39 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  5. Baes, M., Verstappen, J., De Looze, I., et al. 2011, ApJS, 196, 22 [NASA ADS] [CrossRef] [Google Scholar]
  6. Baldwin, J. A., Phillips, M. M., & Terlevich, R. 1981, PASP, 93, 5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  7. Battisti, A. J., Calzetti, D., & Chary, R.-R. 2016, ApJ, 818, 13 [NASA ADS] [CrossRef] [Google Scholar]
  8. Battisti, A. J., Calzetti, D., & Chary, R. R. 2017, ApJ, 840, 109 [NASA ADS] [CrossRef] [Google Scholar]
  9. Bertola, F., Buson, L. M., & Zeilinger, W. W. 1992, ApJ, 401, L79 [NASA ADS] [CrossRef] [Google Scholar]
  10. Bianchi, S. 2008, A&A, 490, 461 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  11. Bianchi, S., De Vis, P., Viaene, S., et al. 2018, A&A, 620, A112 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  12. Binette, L., Magris, C. G., Stasińska, G., & Bruzual, A. G. 1994, A&A, 292, 13 [NASA ADS] [Google Scholar]
  13. Boizelle, B. D., Barth, A. J., Darling, J., et al. 2017, ApJ, 845, 170 [NASA ADS] [CrossRef] [Google Scholar]
  14. Bradley, L., Sipocz, B., Robitaille, T., et al. 2018, astropy/photutils: v0.5 [Google Scholar]
  15. Bron, E., Le Bourlot, J., & Le Petit F. 2014, A&A, 569, A100 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Bruzual, G., & Charlot, S. 2003, MNRAS, 344, 1000 [NASA ADS] [CrossRef] [Google Scholar]
  17. Burstein, D., Bertola, F., Buson, L. M., Faber, S. M., & Lauer, T. R. 1988, ApJ, 328, 440 [NASA ADS] [CrossRef] [Google Scholar]
  18. Calzetti, D., Armus, L., Bohlin, R. C., et al. 2000, ApJ, 533, 682 [NASA ADS] [CrossRef] [Google Scholar]
  19. Camps, P., & Baes, M. 2015, Astron. Comput., 9, 20 [NASA ADS] [CrossRef] [Google Scholar]
  20. Cappellari, M. 2002, MNRAS, 333, 400 [NASA ADS] [CrossRef] [Google Scholar]
  21. Cappellari, M. 2017, MNRAS, 466, 798 [NASA ADS] [CrossRef] [Google Scholar]
  22. Cappellari, M., & Emsellem, E. 2004, PASP, 116, 138 [NASA ADS] [CrossRef] [Google Scholar]
  23. Cappellari, M., Scott, N., Alatalo, K., et al. 2013, MNRAS, 432, 1709 [NASA ADS] [CrossRef] [Google Scholar]
  24. Cardelli, J. A., Clayton, G. C., & Mathis, J. S. 1989, ApJ, 345, 245 [NASA ADS] [CrossRef] [Google Scholar]
  25. Charlot, S., & Fall, S. M. 2000, ApJ, 539, 718 [NASA ADS] [CrossRef] [Google Scholar]
  26. Croxall, K. V., Smith, J. D., Wolfire, M. G., et al. 2012, ApJ, 747, 81 [NASA ADS] [CrossRef] [Google Scholar]
  27. Davies, J. I., Baes, M., Bianchi, S., et al. 2017, PASP, 129, 044102 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  28. Davis, T. A., Krajnović, D., McDermid, R. M., et al. 2012, MNRAS, 426, 1574 [NASA ADS] [CrossRef] [Google Scholar]
  29. Davis, T. A., Alatalo, K., Bureau, M., et al. 2013, MNRAS, 429, 534 [NASA ADS] [CrossRef] [Google Scholar]
  30. Davis, T. A., Young, L. M., Crocker, A. F., et al. 2014, MNRAS, 444, 3427 [NASA ADS] [CrossRef] [Google Scholar]
  31. Davis, T. A., Rowlands, K., Allison, J. R., et al. 2015, MNRAS, 449, 3503 [NASA ADS] [CrossRef] [Google Scholar]
  32. Davis, T. A., Bureau, M., Onishi, K., et al. 2018, MNRAS, 473, 3818 [NASA ADS] [CrossRef] [Google Scholar]
  33. De Geyter, G., Baes, M., Fritz, J., & Camps, P. 2013, A&A, 550, A74 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  34. De Geyter, G., Baes, M., Camps, P., et al. 2014, MNRAS, 441, 869 [NASA ADS] [CrossRef] [Google Scholar]
  35. De Looze, I., Baes, M., Bendo, G. J., et al. 2012, MNRAS, 427, 2797 [NASA ADS] [CrossRef] [Google Scholar]
  36. De Looze, I., Fritz, J., Baes, M., et al. 2014, A&A, 571, A69 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  37. De Vis, P., Gomez, H. L., Schofield, S. P., et al. 2017, MNRAS, 471, 1743 [NASA ADS] [CrossRef] [Google Scholar]
  38. Draine, B. T. 2003, ARA&A, 41, 241 [NASA ADS] [CrossRef] [Google Scholar]
  39. Draine, B. T., & Bertoldi, F. 1996, ApJ, 468, 269 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  40. Emsellem, E., Monnet, G., & Bacon, R. 1994, A&A, 285, 723 [NASA ADS] [Google Scholar]
  41. Fabbiano, G. 1989, ARA&A, 27, 87 [NASA ADS] [CrossRef] [Google Scholar]
  42. Ferrarese, L., Mould, J. R., Kennicutt, Robert C., Jr., et al. 2000, ApJ, 529, 745 [NASA ADS] [CrossRef] [Google Scholar]
  43. Ferrarese, L., Côté, P., Jordán, A., et al. 2006, ApJS, 164, 334 [NASA ADS] [CrossRef] [Google Scholar]
  44. Finkelman, I., Brosch, N., Kniazev, A. Y., et al. 2008, MNRAS, 390, 969 [NASA ADS] [CrossRef] [Google Scholar]
  45. Finkelman, I., Brosch, N., Kniazev, A. Y., et al. 2010, MNRAS, 409, 727 [NASA ADS] [CrossRef] [Google Scholar]
  46. Finkelman, I., Brosch, N., Funes, J. G., et al. 2012, MNRAS, 422, 1384 [NASA ADS] [CrossRef] [Google Scholar]
  47. Fitzpatrick, E. L., & Massa, D. 2007, ApJ, 663, 320 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  48. Fuller, C., Davies, J. I., Auld, R., et al. 2014, MNRAS, 440, 1571 [NASA ADS] [CrossRef] [Google Scholar]
  49. Galliano, F., Galametz, M., & Jones, A. P. 2018, ARA&A, 56, 673 [NASA ADS] [CrossRef] [Google Scholar]
  50. Gordon, K. D., Clayton, G. C., Misselt, K. A., Landolt, A. U., & Wolff, M. J. 2003, ApJ, 594, 279 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  51. Goudfrooij, P., & de Jong T. 1995, A&A, 298, 784 [NASA ADS] [Google Scholar]
  52. Goudfrooij, P., de Jong, T., Hansen, L., & Norgaard-Nielsen, H. U. 1994, MNRAS, 271, 833 [NASA ADS] [CrossRef] [Google Scholar]
  53. Griffin, M. J., Abergel, A., Abreu, A., et al. 2010, A&A, 518, L3 [Google Scholar]
  54. Haario, H., Saksman, E., & Tamminen, J. 2001, Bernoulli, 7, 223 [CrossRef] [MathSciNet] [Google Scholar]
  55. Hirashita, H., Nozawa, T., Villaume, A., & Srinivasan, S. 2015, MNRAS, 454, 1620 [NASA ADS] [CrossRef] [Google Scholar]
  56. Hollenbach, D., & Salpeter, E. E. 1971, ApJ, 163, 155 [NASA ADS] [CrossRef] [Google Scholar]
  57. Hughes, T. M., Foyle, K., Schirm, M. R. P., et al. 2015, A&A, 575, A17 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  58. Jones, A. P. 2004, Astrophysics of Dust, ASP Conf. Ser., 309, 347 [Google Scholar]
  59. Jones, A. P., Fanciullo, L., Köhler, M., et al. 2013, A&A, 558, A62 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  60. Jones, A. P., Köhler, M., Ysard, N., Bocchio, M., & Verstraete, L. 2017, A&A, 602, A46 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  61. Kaviraj, S., Schawinski, K., Devriendt, J. E. G., et al. 2007, ApJS, 173, 619 [NASA ADS] [CrossRef] [Google Scholar]
  62. Kaviraj, S., Peirani, S., Khochfar, S., Silk, J., & Kay, S. 2009, MNRAS, 394, 1713 [NASA ADS] [CrossRef] [Google Scholar]
  63. Kaviraj, S., Ting, Y.-S., Bureau, M., et al. 2012, MNRAS, 423, 49 [NASA ADS] [CrossRef] [Google Scholar]
  64. Kaviraj, S., Rowlands, K., Alpaslan, M., et al. 2013, MNRAS, 435, 1463 [NASA ADS] [CrossRef] [Google Scholar]
  65. Koyama, Y., Shimakawa, R., Yamamura, I., Kodama, T., & Hayashi, M. 2019, PASJ, in press [Google Scholar]
  66. Kulkarni, S., Sahu, D. K., Chaware, L., Chakradhari, N. K., & Pandey, S. K. 2014, New Astron., 30, 51 [NASA ADS] [CrossRef] [Google Scholar]
  67. Mosenkov, A. V., Allaert, F., Baes, M., et al. 2018, A&A, 616, A120 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  68. Nedelchev, B., Sarzi, M., & Kaviraj, S. 2017, MNRAS, submitted [arXiv:1705.07994] [Google Scholar]
  69. Oh, K., Sarzi, M., Schawinski, K., & Yi, S. K. 2011, ApJS, 195, 13 [NASA ADS] [CrossRef] [Google Scholar]
  70. Patil, M. K., Pandey, S. K., Sahu, D. K., & Kembhavi, A. 2007, A&A, 461, 103 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  71. Pilbratt, G. L., Riedinger, J. R., Passvogel, T., et al. 2010, A&A, 518, L1 [CrossRef] [EDP Sciences] [Google Scholar]
  72. Poglitsch, A., Waelkens, C., Geis, N., et al. 2010, A&A, 518, L2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  73. Popescu, C. C., Tuffs, R. J., Völk, H. J., Pierini, D., & Madore, B. F. 2002, ApJ, 567, 221 [NASA ADS] [CrossRef] [Google Scholar]
  74. Rowlands, K., Dunne, L., Dye, S., et al. 2014, MNRAS, 441, 1017 [NASA ADS] [CrossRef] [Google Scholar]
  75. Sarzi, M., Falcón-Barroso, J., Davies, R. L., et al. 2006, MNRAS, 366, 1151 [NASA ADS] [CrossRef] [Google Scholar]
  76. Sarzi, M., Shields, J. C., Schawinski, K., et al. 2010, MNRAS, 402, 2187 [NASA ADS] [CrossRef] [Google Scholar]
  77. Sarzi, M., Falcón-Barroso, J., Davies, R. L., et al. 2017, GANDALF: Gas AND Absorption Line Fitting, Astrophysics Source Code Library [Google Scholar]
  78. Sarzi, M., Iodice, E., Coccato, L., et al. 2018, A&A, 616, A121 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  79. Scharf, C. A., Zurek, D. R., & Bureau, M. 2005, ApJ, 633, 154 [NASA ADS] [CrossRef] [Google Scholar]
  80. Schawinski, K., Thomas, D., Sarzi, M., et al. 2007, MNRAS, 382, 1415 [NASA ADS] [CrossRef] [Google Scholar]
  81. Silva, L., Granato, G. L., Bressan, A., & Danese, L. 1998, ApJ, 509, 103 [NASA ADS] [CrossRef] [Google Scholar]
  82. Skibba, R. A., Engelbracht, C. W., Dale, D., et al. 2011, ApJ, 738, 89 [NASA ADS] [CrossRef] [Google Scholar]
  83. Soto, K. T., Lilly, S. J., Bacon, R., Richard, J., & Conseil, S. 2016, MNRAS, 458, 3210 [NASA ADS] [CrossRef] [Google Scholar]
  84. Stasińska, G., Vale Asari, N., Cid Fernandes, R., et al. 2008, MNRAS, 391, L29 [NASA ADS] [Google Scholar]
  85. Steinacker, J., Baes, M., & Gordon, K. D. 2013, ARA&A, 51, 63 [NASA ADS] [CrossRef] [Google Scholar]
  86. Temi, P., Brighenti, F., & Mathews, W. G. 2009, ApJ, 707, 890 [NASA ADS] [CrossRef] [Google Scholar]
  87. Tielens, A. G. G. M., Waters, L. B. F. M., & Bernatowicz, T. J. 2005, in Chondrites and the Protoplanetary Disk, eds. A. N. Krot, E. R. D. Scott, & B. Reipurth, ASP Conf. Ser., 341, 605 [NASA ADS] [Google Scholar]
  88. Trani, A. A., Mapelli, M., & Ballone, A. 2018, ApJ, 864, 17 [NASA ADS] [CrossRef] [Google Scholar]
  89. van de Voort, F., Davis, T. A., Matsushita, S., et al. 2018, MNRAS, 476, 122 [NASA ADS] [CrossRef] [Google Scholar]
  90. Veilleux, S., & Osterbrock, D. E. 1987, ApJS, 63, 295 [NASA ADS] [CrossRef] [Google Scholar]
  91. Viaene, S., De Geyter, G., Baes, M., et al. 2015, A&A, 579, A103 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  92. Viaene, S., Baes, M., Bendo, G., et al. 2016, A&A, 586, A13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  93. Viaene, S., Sarzi, M., Baes, M., Fritz, J., & Puerari, I. 2017a, MNRAS, 472, 1286 [NASA ADS] [CrossRef] [Google Scholar]
  94. Viaene, S., Baes, M., Tamm, A., et al. 2017b, A&A, 599, A64 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  95. Vidali, G., Roser, J. E., Manicò, G., & Pirronello, V. 2004, Adv. Space Res., 33, 6 [NASA ADS] [CrossRef] [Google Scholar]
  96. Walch, S., Girichidis, P., Naab, T., et al. 2015, MNRAS, 454, 238 [NASA ADS] [CrossRef] [Google Scholar]
  97. Weilbacher, P. M., Streicher, O., & Palsa, R. 2016, MUSE-DRP: MUSE Data Reduction Pipeline, Astrophysics Source Code Library [Google Scholar]
  98. Wild, V., Charlot, S., Brinchmann, J., et al. 2011, MNRAS, 417, 1760 [NASA ADS] [CrossRef] [Google Scholar]
  99. Witt, A. N., & Gordon, K. D. 2000, ApJ, 528, 799 [Google Scholar]
  100. Young, L. M., Bureau, M., Davis, T. A., et al. 2011, MNRAS, 414, 940 [NASA ADS] [CrossRef] [Google Scholar]
  101. Ysard, N., Köhler, M., Jones, A., et al. 2015, A&A, 577, A110 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  102. Zabel, N., Davis, T. A., Smith, M. W. L., et al. 2018, MNRAS, 3074 [Google Scholar]
  103. Zhukovska, S., Henning, T., & Dobbs, C. 2018, ApJ, 857, 94 [NASA ADS] [CrossRef] [Google Scholar]

All Tables

Table 1

Best-fit model parameters for the ISM in FCC 167.

All Figures

thumbnail Fig. 1

Different stages of the MGE modelling to produce attenuation maps for every wavelength slice of the MUSE datacube. Left panel: reconstructed image at 5509 Å (the V -band central wavelength). The green ellipses indicate the inner and outer edge of the dust lane (red contour). Middle panel: corresponding dust-free MGE model fit. Right panel: attenuation map at 5509 Å obtained from Eq. (1).

In the text
thumbnail Fig. 2

Spectrally resolved attenuation curve for FCC 167 (red line, this work) and for NGC 5626 (green line, Viaene et al. 2017a). Three reference curves are plotted: The Calzetti et al. (2000) and Battisti et al. (2017) average attenuation curves for different samples of galaxies (blue solid and dashed line, respectively) and the THEMIS extinction curve (black dotted line). All curves are normalised to the V -band central wavelength; AV = 0.35 for FCC 167. For reference, the emission spectrum of the dust lane area is shown in the bottom panel. Several important emission lines (shaded green) and the NaD absorption line (shaded purple) are highlighted as well.

In the text
thumbnail Fig. 3

Parameter maps for the fitted pixel-by-pixel attenuation curves. Left panel: V -band attenuation AV. Right panel: slope parameter δ, where steeper slopes are associated with more negative δ and a Calzetti-like attenuation curve has δ = 0.

In the text
thumbnail Fig. 4

Ionised gas in FCC 167. Top left panel: surface brightness (in log scale) of the [NII] λ6583 line already highlighting a tail of ionised gas connecting to the main reservoir. Top second panel: velocity map indicating organised motion along the major axis; there are several additional kinematic features towards the outskirts. Top third panel: an axisymmetric kinematic model of the measured velocity map. Top right panel: velocity dispersion of the [NII] line showing a central peak and hints of a spiral-like shape. Bottom left panel: reddening, E(BV), by diffuse dust (as opposed to emission-line reddening) derived from the spectral fitting. Blue contours correspond to molecular gas emission measured by ALMA (see Fig. 5). Bottom second panel: reddening, E(BV), of the emission lines which is much stronger than the diffuse reddening component, and without any distinct features. Bottom third and right panels show the [NII]/Hα line ratio and Hα EW, respectively. Together these two maps indicate low level star formation inside the dust ring and nuclear ionisation due to an AGN.

In the text
thumbnail Fig. 5

Moment maps of the CO(1-0) line observations with ALMA. Left panel: total line intensity with the region of the dust lane overlaid. Middle panel: velocity map indicating a regular rotation. Right panel: linewidth map (which peaks in the centre owing to beam smearing). The size of the ALMA beam is shown in each frame as a black ellipse.

In the text
thumbnail Fig. 6

Best-fit models for a ring of dust with a THEMIS MW dust mix (top row) and for a THEMIS bare sillicate dust mix (bottom row). We only show a representative image slice at the V -band central wavelength (5509 Å) as the results are very similar across the spectrum. Both observed images in the first column are the same. The second column shows the corresponding model images and the third column shows the residual maps following IobsImodIobs$\frac{I_{\text{obs}}-I_{\text{mod}}}{I_{\text{obs}}}$ in pixel flux.

In the text
thumbnail Fig. 7

Attenuation curves of the best-fit radiative transfer models for a dust ring geometry. The THEMIS MW dust mix produces a steeper curve (blue line) than a dust mix with only bare silicate grains (orange line). The dashed lines correspond to the corresponding intrinsic extinction curves of the dust mixes, normalised to the model attenuation at 5509 Å. The observed attenuation curve of FCC 167 (from Sect. 3) is shown in black.

In the text
thumbnail Fig. 8

Panchromatic radiative transfer model predictions based on the best-fit FitSKIRT models for THEMIS MW dust and the bare silicate dust mix. A single (10 Gyr) stellar population was used for the intrinsic SED of the stellar component. No additional dust heating sources are adopted in the model. The green line represents a modified black-body fit to the Herschel fluxes from 100 to 500 μm; a free β matches the steep Rayleigh-Jeans slope.

In the text
thumbnail Fig. 9

Fit tothe Chandra X-ray spectrum of FCC 167. Top panel: observed photon counts per channel (crosses). The black line shows the total model flux. Bottom panel: unfolded model (solid line) and its separate components (dotted line). A wa*po+mekal model was used to represent a power law and a hot gas line spectrum. The unfolded model is adjusted for spectral resolution and filter response function before comparison to the data in the top panel.

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.