A&A 477, 105-116 (2008)
DOI: 10.1051/0004-6361:20077925
C. del Burgo1 - D. Carter2 - G. Sikkema3
1 - School of Cosmic Physics, Dublin Institute for Advanced
Studies, Dublin 2, Ireland
2 -
Astrophysics Research Institute, Liverpool John Moores University,
12 Quays House, Egerton Wharf, Birkenhead, CH41 1LD, UK
3 -
Kapteyn Astronomical Institute, University of Groningen, PO Box 800,
9700 AV Groningen, The Netherlands
Received 22 May 2007 / Accepted 8 October 2007
Abstract
Aims. Shells in Ellipticals are peculiar faint sharp edged features that are thought to be formed by galaxy mergers. We determine the shell and dust distributions, and colours of a well-resolved shell and the underlying galaxy in NGC 5982, and compare the spatial distributions of the dust and gas phases, in order to investigate the origin of shells and dust in NGC 5982.
Methods. We use Spitzer data in the wavelength range from 3.6 to 160
m and HST/ACS optical data. Surface photometry, galaxy models and residual images are determined using IRAF task ELLIPSE and GALPHOT task ELLIPFIT. After subtracting the galaxy models, residual images are used to identify the shells. Excess emission is identified by subtracting a scaled stellar photosphere emission template. Colours
[3.6]-[4.5],
[3.6]-[8.0] and
[5.8]-[8.0] of the underlying galaxy and
[5.8]-[8.0] of the excess emission are obtained using techniques similar to those of Pahre et al. The
[3.6]-[4.5] and V-I colours of a well-resolved shell are determined using techniques of Sikkema et al.
Results. Shells are visible in the 3.6
m image, fainter at 4.5
m and undetected at 5.8
m. We find two new shells, the outermost identified so far in NGC 5982. All shells apart from shell 24 are on the photometric major axis. Shell 24 is the only one for which reliable colours are obtained. It has colours
[3.6]-[4.5] and V-I bluer than the underlying galaxy. Excess emission at 4.5, 5.8 and 8
m is widely extended. The ratio of excess to total emission decreases towards the centre, which could be related to a lower mass loss rate from AGB stars. The
[3.6]-[4.5] colour indicates a major contribution from late-type (K2-M0III) stars towards the centre. Dust traced by 24
m emission follows the stellar spatial distribution, supporting a circumnuclear origin. We find extended 160
m emission from cold dust, possibly forming a disk inclined to the principal axes, with a mass of a few
.
The warmer big grains and/or smaller grains traced by the 70
m emission are confined to a small region elongated along the minor axis and shifted a few kpc from the centre.
Conclusions. We detect for the first time shells from mid-infrared data. The very different distributions of dust, warm gas and HI gas together with the presence of shells and a kinematically decoupled core suggest a minor merger in NGC 5982.
Key words: galaxies: elliptical and lenticular, cD - galaxies: photometry - galaxies: ISM - galaxies: formation - galaxies: interactions
Shell galaxies (Malin & Carter 1980, 1983) are elliptical galaxies with faint, often sharp edged features in their envelopes, which are generally interpreted as the remnants of mergers of low mass, low velocity dispersion galaxies with the elliptical (Quinn 1984; Hernquist & Quinn 1988). The study of the structure and morphology of these features is useful in interpreting the parameters of minor mergers, and the importance of mergers in determining the present day structure of ellipticals (Dupraz & Combes 1986; Nulsen 1989).
The interstellar medium (ISM) in ellipticals mostly consists of X-ray emitting
hot gas (
K) but it may also contain significant amounts of
warm and cold gas and dust. Recent studies show that approximately 70% of
ellipticals harbour HI gas (Morganti et al. 2006); a similar proportion
contains warm gas (Sarzi et al. 2006). For a similar galaxy sample,
Combes et al. (2007) find a CO detection rate of 28% in the central regions. The far-infrared emission in ellipticals
gauges the bulk of dust, which has masses from 10
to 10
(Temi et al. 2004, 2007).
Our knowledge of the far-infrared emission from early-type galaxies
is largely based on Infrared Astronomical Satellite (IRAS) data at 12, 25, 60 and 100
m. IRAS was successful in detecting emission
from dust with a wide range of temperatures in many bright elliptical
galaxies (Knapp et al. 1989). These observations were limited, however,
by sensitivity, angular resolution (
4
on the sky) and the lack of
filterbands to gauge and resolve the emission from cold dust.
Infrared Space Observatory ISO
improved this situation, providing
1.5
angular resolution
observations at 150 and 200
m which are more sensitive to the peak emission
from ellipticals. It also included filterbands comparable to
IRAS and mid-infrared imaging at an angular resolution of a few arcseconds,
to study the spatial distribution of the mid-infrared carriers
(see Ferrari et al. 2002; Xilouris et al. 2004).
ISO could not always resolve the structure in elliptical galaxies,
which may have emission from Active Galactic Nuclei (AGNs), dense clouds, circumstellar
dust and diffuse interstellar dust (Temi et al. 2004). The modern Spitzer instruments
IRAC (filterbands at 3.6, 4.5, 5.8, 8,0
m)
and MIPS (broad filterbands at 24, 70, 160
m) offer
unprecedent sensitivity and angular resolution to much better differentiate the
different dust components in a sample of ellipticals with a wide luminosity range
(see, e.g., Pahre et al. 2004, hereinafter PAFW04). Spitzer IRS is used to
observe the spectral energy distribution (SED) of the mid-infrared carriers in
early-type galaxies (Kaneda et al. 2005; Bressan et al. 2006).
It is thought that dust is produced by evolving red giant stars in the old stellar population
typical of ellipticals, and it is dispersed into the interstellar medium (ISM).
Conversely, dust is destroyed by sputtering in the host interstellar hot gas of
ellipticals. During their sputtering lifetimes big dust grains in the ISM
are heated by the interstellar radiation field and thermal electron
impacts to temperatures that radiate in the far-infrared (see Temi et al. 2004,
and references therein). Small optically
thick clouds have been also observed in
50% of bright ellipticals
(van Dokkum & Franx 1995; Lauer et al. 2005). These clouds are likely formed
from dusty gas expelled from stars in the central Kiloparsec of ellipticals
(Mathews & Brighenti 2003).
Dust can be also produced in galaxy mergers and it is largely accepted that
ellipticals are formed by mergers. However, it is not straightforward to
explain all kinematical and structural properties in ellipticals by the
current merger models.
NGC 5982 is a shell galaxy with regular major axis shells (Sikkema et al. 2007, hereinafter SCP07), which identifies it as type I in the classification scheme of Prieur (1988). The redshift is cz=3017 km s-1 (Smith et al. 2000), which corresponds to a distance D=43 Mpc assuming H0=70 km s-1 Mpc-1. NGC 5982 is the brightest galaxy of LGG 402 with 4 members (García 1993). There is a possible very low activity level LINER type nucleus (Ho et al. 1997). This galaxy harbours a Kinematically Decoupled Core (KDC: Wagner et al. 1988), and probably a rotating central disk (Emsellem et al. 2004). Outside the core, the galaxy is slowly rotating (Emsellem et al. 2004) and presents central boxy isophotes (SCP07).
This paper presents a photometric analysis of the shells and the underlying galaxy in the elliptical (E3) galaxy NGC 5982 from Spitzer infrared data, and compares this with our earlier analysis of HST/ACS V- and I-band photometric data (SCP07). The observations and data processing are presented in Sect. 2. Section 3 is devoted to the analysis of the data. The results and their interpretation are presented in Sect. 4. Section 5 contains the discussion of our results, and finally, Sect. 6 presents our conclusions.
We use public Spitzer and Hubble Space Telescope (HST) data. The first were retrieved from the Spitzer Data Archive using the Archive Tool Leopard. The HST observations belong to our programme to study shell ellipticals (GO9399, PI: D. Carter) and are described in two earlier publications (Sikkema et al. 2006; SCP07).
NGC 5982 has been observed using the Infrared Array Camera (IRAC;
Fazio et al. 2004) on the Spitzer Space Telescope (Werner et al. 2004)
during the campaign IRAC005700 of GO programme number 3403 (PI: J. Surace).
NGC 5982 was imaged using 14 exposures of 26.8 s (effective exposure
time per pixel) centred on different positions with the four IRAC filterbands
at 3.6, 4.5, 5.8 and 8
m.
We use the standard pipeline
(version 14.0) post Basic Calibrated Data (post-BCD) products.
The mosaics resulting from
the addition of the 14 exposures have effective areas between 103 and 104 arcmin2 and they are less noisy in an oversampled
7
square region centred on NGC 5982.
The image mosaics at 3.6 and 5.8
m roughly sample
that square region and a similar-sized region to the ENE side;
those at 4.5 and 8.0
m cover the square region and
a similar sized region to the WSW side.
These image mosaics are large enough to derive
reliable values of the sky background to subtract from the images.
The angular resolution of the IRAC images is
.
The existing calibration of IRAC data is based on point-source
observations (Fazio et al. 2004). Extended sources require a different
calibration. We use a similar approach to PAFW04
in order to determine instrumental magnitudes relative to Vega.
The aperture corrections of PAFW04 are based on large
aperture photometry of NGC 777 assuming that it has no dust emission
and that the 2MASS near-infrared and IRAC colours match those
corresponding to a M0III star. This is a suitable stellar type to match
the SED from an early-type galaxy in the K-band
(Frogel et al. 1978). We also use the aperture corrections
of Reach et al. (2005, hereinafter RMC05) and compare the results
for both aperture corrections.
The most significant differences are for the aperture
correction at 5.8
m. See Sect. 3.2 for a
detailed discussion.
![]() |
Figure 1:
Surface brightness at 3.6 |
| Open with DEXTER | |
The angular resolution of the MIPS images is
,
18
and
40
at 24, 70 and 160
m, respectively.
The same observations have already been used
by Temi et al. (2007) to derive luminosities of NGC 5982
at 24, 70 and 160
m of
L24=1040.86 erg s-1,
L70=1040.92 erg s-1 and
L160=1040.66 erg s-1.
![]() |
Figure 2:
From top to bottom excess emission at 4.5, 5.8 and 8 |
| Open with DEXTER | |
HST ACS_WFC camera data (Pavlovsky et al. 2004) obtained
with the filters F606W (V band) and F814W (I band) are used here. For a detailed
description of the observations and the data processing see Sikkema et al. (2006)
and SCP07, here we just provide a summary.
The field of view was of
,
with a
pixel size of 0.049 arcsec. The exposure times were of 1314 and 1020 s
in the V- and I-bands, respectively. The 80% completeness levels were
of 26.08 and 24.80 mag in V- and I-bands, respectively. In this paper we
use the V- and I-band images to extract shell profiles in these filterbands
(see Sect. 4.4) that were not presented in SCP07. We use the technique
introduced by SCP07 to determine the surface brightness locally in wedges for
the underlying galaxy and the shells.
Surface brightness images of NGC 5982 of
in size on the sky are
presented in Figs. 1a-d, f-h corresponding to the emission
at 3.6, 4.5, 5.8, 8.0, 24, 70 and 160
m, respectively.
At first view, the surface brightness distribution
at 3.6, 4.5, 5.8 and 8.0
m is very similar.
The 24
m surface brightness distribution is also very similar to them,
but very different from those at 70 and 160
m. The latter two are elongated along the
photometric minor axis of the galaxy, with the 70
m emission not resolved
along the perpendicular axis.
There is a double peak emission at 70
m that is shifted
with respect to the photometric centre, with the bulk of the emission to the North.
The 160
m emission is more extended than that at 70
m, with
emission to both North and South of the centre.
The sky background is subtracted from each image. Hereafter, we refer to each image corrected for the sky background as galaxy emission, observed emission, or simply observed image.
A scaled stellar photosphere emission template is subtracted from the 4.5, 5.8 and 8.0
m images to obtain any excess emission
at these wavelengths.
We follow a similar approach to PAFW04 but using
the image at 3.6
m only instead of a combination of the images at
3.6 and 4.5
m as our stellar template.
We also use the GALPHOT galaxy model at 3.6
m (see next section) as a stellar template.
For each template, we use the aperture corrections provided by PAFW04
as well as those of RMC05 and compare the results.
The stellar template is scaled to match the corresponding theoretical colour
of M0III stars, before being subtracted from the observed emission at 4.5, 5.8 and 8.0
m (PAFW04).
Figure 2 shows the excess emission at 4.5, 5.8 and 8.0
m in NGC 5982
for three of the four possible combinations of templates and aperture corrections.
Figure 1e shows the contours of the
excess emission of NGC 5982 at 8.0
m obtained with the image at 3.6
m as
stellar template and the RMC05 aperture corrections for comparison with
the observed images shown in Fig. 1.
We highlight some remarks resulting from our comparison of the excess emission obtained with the two stellar templates combined with each of the aperture corrections:
![]() |
Figure 3:
Global parameters for NGC 5982:
a) surface brightness at 3.6 |
| Open with DEXTER | |
![]() |
Figure 4:
Residual images of NGC 5982 at 3.6 |
| Open with DEXTER | |
Surface photometry at wavelengths 3.6, 4.5, 5.8 and 8.0
m
is obtained using an ellipse fit with the task ELLIPFIT of
GALPHOT (Jørgensen et al. 1992) and the task ELLIPSE from the IRAF STSDAS package.
The latter is used to measure the flux using the same isophotal
shapes defined from a reference image (i.e., the 3.6
m image).
ELLIPFIT provides us a galaxy model and a table with the fitted parameters, including
ellipticity (
), position angle (PA), surface brightness (
)
and the coefficients C3, C4, S3 and S4 (Carter 1978) for different distances
from the centre of the galaxy.
Point-like sources, background galaxies, dust features,
bright shells and artifacts are identified and masked by hand in an
iterative way. Remaining faint shell structures have a
negligible effect on the results. The best fits are obtained by allowing
centre, PA and ellipticity to vary with radius.
There is an excellent agreement between the values
(e.g. of
and PA) obtained using ELLIPSE and ELLIPFIT.
The inner region (
)
is poorly
fitted compared with the rest of the galaxy, especially for the 3.6
m image.
The differences between the observed image and galaxy model are up to
3%.
This fitting problem was also noted at optical wavelengths by SCP07,
who suggested that it could be due to the KDC, which is
10
in size
(Mc Dermid et al. 2006). Alternatively it may be related to the presence of dust
or the high density of point-like objects in the innermost region of NGC 5982.
The ellipticity gets significantly lower in the inner 4
of the
galaxy and it was not possible to fit smoothly the central region. We obtain the best
result using 6 harmonics. The values of the C4 coefficient for
are very similar to those determined from the
HST/ACS I-band image (SCP07), and indicate a boxy structure (Carter 1978).
Global surface brightness profiles for NGC 5982 are shown by plotting
for
each fitted ellipse as function of semi-major axis a. Figures 3a, c
and d respectively shows
,
and PA versus a at 3.6
m.
To investigate the effect of the IRAC Point Source Response Function (PSF)
(
)
upon the surface brightness, ellipticity and position angle
profiles, we measure
,
,
PA, S4 and C4 in the HST/ACS V-band image
and also in the image resulting from convolving it with the IRAC PSF at 3.6
m.
We show only the results for
and PA (see Figs. 3c and d,
respectively). The ellipticity drops toward smaller radii partly
as result of degradation of the angular resolution. The effect on PA is
much smaller.
We use a Sérsic (1968) profile to model the luminosity at 3.6
m.
This profile is defined by the equation:
![]() |
(1) |
Table 1 lists the least-squares fitted values for
,
n
and
for
.
The difference
between the luminosity profile at 3.6
m and the fitted model is plotted in
Fig. 3b.
We compute residual images from the difference between the observed
image at 3.6, 4.5 and 5.8
m and the corresponding galaxy models
at these wavelengths obtained from ELLIPFIT (see previous section).
Figure 4 (left) shows the residual
image at 3.6
m where faint shell-like structures are observed.
The innermost shells are much less resolved in the Spitzer images than in
the R-band (see Fig. 4, right) and the HST/ACS images.
All the shells are much fainter at longer wavelengths. The residual
image at 5.8
m does not show any shells.
We note that the innermost shells are in a region where some 4.5
m excess emission is
observed (see Figs. 2 (top) and 4).
Residual images are only used to identify and locate the shells. For the determination of the shell profiles and colours we work within wedges as in SCP07 to avoid the effect of large scale fluctuations in the sky background that affect the extraction of faint shell properties.
Colour profiles
[3.6]-[4.5],
[3.6]-[8.0],
[5.8]-[8.0] for the underlying stellar emission
and the
[5.8]-[8.0] colour for the excess emission (see Sect. 3.2)
are determined using the isophotal shapes obtained from the 3.6
m
image. The colours are, therefore, measured for the same regions of the galaxy.
Figures 3e-g show the colours
[3.6]-[4.5],
[3.6]-[8.0]and
[5.8]-[8.0], respectively, using the RMC05 aperture corrections.
Figure 3g also shows the
[5.8]-[8.0] colour of the excess emission obtained
from the PAFW04 aperture corrections. The colour profiles are very sensitive
to the sky background subtracted. Colours
[3.6]-[8.0] and
[5.8]-[8.0]
are especially sensitive to the sky background emission at 8.0
m that may
modify the colour profile from
.
Colours are also
sensitive to the aperture corrections (see Sect. 4.1),
most significantly in the case of
the
[5.8]-[8.0] colour of the excess emission (see Fig. 3g).
However, the differences in this colour using the two different
stellar templates (see Sect. 3.2) are not significant.
Azimuthally averaged surface brightness profiles practically wash away
the shells visible in Fig. 4.
The shells are too faint in NGC 5982 to affect the colour
profiles (see Figs. 3e-g). As mentioned above,
we use the technique of SCP07 based on wedges to extract faint shell colours
(see Sect. 4.4).
Table 1:
Surface photometry at 3.6
m and 4.5
m.
The isophotes of the excess emission at 5.8 and 8.0
m show a slightly different orientation
of the major axis in comparison with the total emission (i.e., stellar
plus excess emission) at these wavelengths (see Figs. 1d and e).
The emission profiles at 5.8 and 8
m are significantly different from
those at 3.6
m. This is due to the observed excess emission at 5.8 and 8
m.
For each wavelength, the surface brightness ratio between the excess emission and the
total emission varies with the radius. The change is different for the two aperture corrections
used here. If applying those of PAFW04, the excess emission at 8
m is
17% and
37% of the total emission in the centre and for radii between
4
and
25
on the sky, respectively.
These values are 3% higher when using the RMC05 aperture corrections.
At 5.8
m, the values in the centre and outer region are respectively of
5% and
20% when using the PAFW04 aperture corrections. The values are 10% higher
when using the RMC05 aperture corrections.
The excess emission at 4.5
m is more elongated that that at 5.8 and 8
m
(see Fig. 2), and represents
1% and
5% of the total
emission at this wavelength when using the PAFW04 and RMC05 aperture corrections,
respectively. The excess emission in the centre is significantly fainter, and is
undetectable when using the RMC05 aperture corrections.
We discuss in Sect. 5.1 the possible cause of the observed decrease
in the excess emission towards smaller radii.
The
[5.8]-[8.0] colour of NGC 5982 for the total emission
is moderately red. Once the dominant stellar emission is subtracted,
the
[5.8]-[8.0] colour of the excess component is redder at smaller radii.
If using the PAFW04 aperture corrections, the
[5.8]-[8.0] colour
of the excess emission approaches the values found in the
star-forming galaxy M 81 (Willner et al. 2004) for
,
and it reaches a value as high as
mag at
.
A value of
2.06 mag is predicted by
Li & Draine (2001) for the Polycyclic Aromatic Hydrocarbon (PAH) emission, for
a range of 5 orders of magnitude for UV radiation field.
However, the
[5.8]-[8.0] colour in the inner region is significantly less red
when using the aperture corrections of RMC05 (see Fig. 3g).
The
[5.8]-[8.0] colour profile of NGC 4203 (PAFW04)
is also sentitive to the aperture correction and shows a similar trend to NGC 5982
when using the same aperture correction.
The mid-infrared colours provide information about the mean spectral type of the composite stellar population in the galaxy. The [3.6]-[4.5] colour (see Fig. 3e) is redder towards larger radii in NGC 5982. The colour change is comparable to a small change in the stellar population among K giant stars, typical for early-type systems. The [3.6]-[4.5] colour for the outer regions is similar to that of KOIII stars. Towards inner regions the [3.6]-[4.5] colour approaches values corresponding to K2-M0III stars.
The
[3.6]-[8.0] colour is nearly constant in the inner bulge and
similar to that found by PAFW04 in NGC 4203.
However, NGC 5982 does not present the same strong gradient from
.
The colour gradient from
is very sensitive to the subtracted sky background
at 8
m.
![]() |
Figure 5:
Emission profile at 160 |
| Open with DEXTER | |
The dust carriers responsible for the emission at 24
m follow the same distribution as
the stellar component emitting at 3.6
m, which suggests that the dust is circumstellar.
Emission at 24
m in the diffuse ISM of the Milky Way is due to very small grains
(Léger & Puget 1984; see also review of Draine 2003).
At 18
m there is a broad silicate feature which is produced by Asymptotic Giant Branch
(AGB) stars in ellipticals (Bregman et al. 2006).
The emission images of NGC 5982 at 70 and 160
m reveal a different distribution for the
corresponding carriers at these wavelengths from those responsible for the emission
at shorter wavelengths. Figure 5 shows the 160
m emission profile along axis
together with the 3.6
m emission profile smoothed to the same angular
resolution.
The lack of correspondence between the 70 and 160
m images hinders the
estimation of dust temperature. The bulk of the dust is traced by the 160
m emission.
The carriers at this wavelength are likely to be large cold silicate grains which may have ice mantles
(see, e.g., Draine 2003). The peak emission of warmer silicate grains will move towards
shorter wavelengths.
In the Milky Way emission at 70
m mostly traces very small grains (Draine 2003).
If the dust emitting at 160
m has a temperature of
20 K
and properties similar to those found in most elliptical galaxies from ISO data
(Temi et al. 2004) we estimate a dust mass of a few times
.
This is in the mass range derived from Spitzer data for a sample of ellipticals
(Temi et al. 2007).
At 160
m we find a bright extended source with a peak emission of 4.8 MJy sr-1, spatially co-incident within the errors with 2MASX J15383273+5924022
from the 2MASS catalog of extended sources (Cutri & Skrutskie 1998).
This has counterparts at 3.6, 4,5, 5.8 and 24
m according to the images analyzed here.
No redshift is available for this source, but it may well be in the background.
Shell structures are visible in the 3.6
m residual image
(see Fig. 4, left). These correlate well with the shells found in
SCP07, in particular in the ground-based R band image shown in their
Fig. 5 and our Fig. 4 (right).
Below 40 arcsec radius, due to the lower angular resolution of the IRAC images,
the shells are not as clear as in the HST/ACS V- and I-band data.
In general, the shells are much fainter in the 4.5
m residual image than in the
3.6
m residual image. We can not distinguish any shells at 5.8
m. To the best
of our knowledge this is the first time that shells are observed in the mid-infrared.
We have discovered two probable outer shells from the mid-infrared emission,
labelled 25 and 26 (see Fig. 4 left) which are not visible in
the HST/ACS images (shell 26 is mostly outside the field of view),
although they do appear to be visible in the ground-based R-band image
(see Fig. 4 right).
These shells are circular arcs, located at distances of about 126
(shell 25) and 150
(shell 26) on opposite sides along the major axis. Shell 24,
the outermost shell detected in the HST/ACS images, is
clearer in both the 3.6 image and the R-band image.
This is the only shell to lie on the minor axis of this galaxy.
![]() |
Figure 6:
Left: mean surface luminosity in circular wedges with an opening angle
of |
| Open with DEXTER | |
![]() |
Figure 7:
Shell profiles in the V-band ( left) and at 3.6 |
| Open with DEXTER | |
Using the methodology outlined in SCP07, we determine the luminosity profile of
shell 24 by performing surface photometry in circular and elliptical wedges,
of a constant opening angle chosen to fully include the azimuthal extent of the
shell, centred on the photometric centre of the galaxy, and extending to beyond
the shell radius. Figure 6 (left) shows the resulting surface brightness
profiles at 3.6 and 4.5
m of shell 24 using circular wedges.
Figure 6 (right) shows the residuals (i.e., difference between the
observed surface brightness and the galaxy model extracted within wedges)
versus radius (i.e., distance from the photometric centre)
at 3.6 and 4.5
m as well as in the R-band for shell 24.
The shell profiles are very similar at these wavelengths and are in good agreement
with those derived from HST/ACS V- and I-band data.
An analysis of the shell profiles to learn about their 3D density
distribution is carried out using the methodology described in SCP07.
In the first model, it is assumed that the shell has
an intrinsic
(r0-r)-0.5 density distribution, where r0 is the radius of maximum
stellar density. This model is predicted by theory to describe the inner
parts of phase wrapped shells (Dupraz & Combes 1986; Prieur 1988). In the second model,
it is assumed that shells have an intrinsic Gaussian shape
with thickness
.
Although this model lacks a physical basis, it has
been used to describe spatially wrapped shells (Turnbull & Collet 1999).
Figure 7 shows the model fits for the V-band and at 3.6
m.
The Gaussian model is a better fit to the observations at all the wavelengths.
The first model suffers from "noise'' and a sudden drop in flux at r0 due to its
discontinuous nature, and to the noise in the data (see SCP07).
The mean colours of shell 24 are obtained. The
[3.6]-[4.5] colour of shell 24
is determined using both circular and elliptical wedges.
In each case we find a consistent mean colour difference between
the shell and the underlying galaxy of
mag, in the sense that the shell
is bluer than the underlying galaxy.
The mean colour of the underlying galaxy is
[3.6]-[4.5]=-0.08 mag at the shell position.
The V-I colour of this shell and the underlying galaxy at this position
are also determined. We find that that V-I is
mag and
1.25 mag
for the shell and the underlying galaxy respectively,
indicating that the shell is also bluer than the galaxy at optical wavelengths.
The stellar population in NGC 5982 is dominated by old stars. The galaxy
is classified as type YP, i.e., containing a young stellar
population mixed with the old one, by Michard & Prugniel (2004).
Mc Dermid et al. (2006) found that the stellar population in the inner 2 kpc, corresponding to
the size of the KDC, has a luminosity weighted mean age of
8 Gyr. Mc Dermid et al. also
claim that the KDC is unlikely to be the result of a
recent merger, although formation from a recent "dry'' merger is not ruled out.
Sikkema et al. (2006) observed that the globular cluster populations in NGC 5982 are
consistent with an old (>5 Gyr) and evolved system.
We find that the surface brightness profile at 3.6
m resembles that
corresponding to the HST/ACS I-band image. This indicates that the emission
from old stellar photospheres dominates the bulk emission at these wavelengths.
This further supports our assumption that the 3.6
m emission is entirely
photospheric and unaffected by dust emission, and therefore it can be
used as a stellar template to determine any excess emission at 4.5, 5.8
and 8.0
m (see Sect. 3.2).
We find that the emission at 4.5
m is mostly photospheric, with a faint excess emission.
The
[3.6]-[4.5] colour is bluer towards the inner regions, with
a profile that matches the colours of late K or M giants (see Fig. 3e).
Inner regions could have a major contribution from low-mass AGB stars with lower
mass-loss rates, which would explain their low excess to total emission
ratios (see Sect. 4.1).
Colour
[3.6]-[8.0] is too red to be explained by any stellar population.
This may be a signature of the extended excess emission detected in NGC 5982.
A similar value for
[3.6]-[8.0] was found in NGC 4203 by PAFW04 who interpreted
it as due to warm dust emission at 8.0
m.
Systematic uncertainties in the aperture corrections,
especially at 5.8
m, hinder the analysis and interpretation of
the excess emission colours (see Sects. 3.2
and 4.1).
The excess emission at 8
m obtained using the two aperture corrections
is very similar. However, the ratio of excess to
total emission at 5.8
m is considerably higher when using
the aperture corrections of RMC05 and it is more difficult
to explain its origin since the emission below 8
m is expected
to be dominated by stellar photospheres (Bressan et al. 2006).
The mid-infrared SEDs of passive early-type galaxies contain emission
from mass-losing oxygen-rich AGB stars, with a silicate emission feature
at 10
m (Bressan et al. 1998, 2006; Bregman et al. 2006).
Some excess emission at 5.8 and 8
m may originate in the moderately
thick circumstellar envelopes of the AGB stars,
where silicate dust grains are formed and heated to
300 K
(Bregman et al. 2006).
The IRAC filterband at 8
m is sensitive to the 9-12
m emission from
silicates. At shorter wavelengths the continuum is produced near the stellar
photosphere (Bregman et al. 2006). The lower excess emission observed in
the centre would be due to a lower mass-loss rate from AGBs (see previous
section).
Excess emission at 5.8 and 8
m may be due to
PAH emission. There is a series of emission bands, emitting
at 3.3, 6.2, 7.7, 8.6, 11.3, 12.7 and 16.4
m, observed in
interacting early type galaxies (NGC 4435: Panuzzo et al. 2006), starbursts
and AGNs (M87: Bressan et al. 2006). Bands at 6.2, 7.7
m,
8.2 and 8.6
m lie within the IRAC filterbands at 5.8 and 8
m.
The central
[5.8]-[8.0] colour in NGC 5982 obtained using the aperture
corrections of PAFW04 resembles that from PAH emission. This is not confirmed,
however, when using the RMC05 aperture corrections.
The excess emission could be due to very warm dust heated by a central AGN. It is not clear, however, whether there is a central LINER in NGC 5982 (Ho et al. 1997). A low level AGN may influence its surroundings. UV radiation from an active nucleus could be hard enough to destroy PAHs and even very small grains. Only big grains would survive. The dust destruction by the AGN may be restricted, however, to a few hundred parsecs, i.e., a region well within the IRAC PSF. Sputtering by hot gas can affect regions much farther out in the galaxy. On the other hand, a depletion of small particles could be also produced by dust coagulation, with the smaller particles being adsorbed onto the surface of big grains, as is observed in dense regions of the Milky Way (see, e.g., Stepnik et al. 2003; del Burgo et al. 2003, and references therein).
The lack of definitive observations to determine the level of AGN activity
makes our discussion more difficult. It is plausible
that AGB mass-loss could largely be responsible from the extended excess
emission in NGC 5982, as is observed in many other early-type galaxies
(Bressan et al. 2006). This is also in agreement with 6.7 and 15
m
extended emission observations in early-type galaxies (Xilouris et al. 2004).
We find similar surface brightness distributions at 3.6
m and 24
m.
This suggests that the 24
m emission is mostly due to circumstellar dust
associated with the stellar component traced at 3.6
m. In the diffuse ISM
of our Galaxy emission at 24
m is interpreted as due to very small
(a few nanometer in size) grains that are transiently heated by the stellar radiation
field and have a peak emission at
60
m. Such very small grains are likely
to be the dust population traced by the 24
m emission in NGC 5982.
Big silicate grains can be traced by
the Spitzer/MIPS filterbands at 70 and 160
m. Temi et al. (2004) used a
two dust temperature model to fit the far-infrared emission in a rather inhomogeneous
ISO sample of 53 early-type galaxies. They found temperatures of around 40 K and 20 K,
the latter being the dominant component by mass. Generally, big grains would present
a continuous temperature range.
In dense/cold regions of the Milky Way, the dust absorption/coagulation of small grains
onto the surfaces of big grains with ice mantles is favoured. This results in grains
with enhanced far-infrared emissivities with respect to the grains
in the diffuse interstellar medium (Stepnik et al. 2003; del Burgo et al. 2003).
The different surface brightness distributions observed in NGC 5982
at 70 and 160
m may be partly due to changes
in the abundances and properties (temperature, emissivity) of dust.
This far-infrared emitting dust, consisting of grains of various
size and temperature, appears to be distributed
in a disk inclined to the principal axes.
If the average emissivity of the cold big grains were enhanced, as has
been proposed to explain the SEDs in spirals (Dasyra et al. 2005), less
dust would be necessary to account for the far-infrared emission.
The surface brightness distributions at 70 and 160
m suggest a dust
source additional to mass-loss from the underlying stellar population.
In Sect. 5.4 we discuss the origin of this dust.
Shells in NGC 5982 are expected to be the result of phase wrapping. The well defined type I shells in NGC 5982 are thought to be created by a minor merger with a small, non-rotating galaxy.
The large number of shells found in NGC 5982 indicates an old age.
We have revised the dynamical age of the shells in this galaxy
found by SCP07 taken into account our discovery of the outermost shells.
The dynamical age, estimated from Eq. (88) of
Nulsen (1989), using the outer shell radii found in this paper,
and assuming a mass for NGC 5982 of
from
the fundamental plane method (Eq. (7) of van Dokkum & Stanford 2003)
is
years.
We analyze the colours [3.6]-[4.5] and V-I in the shell 24, the only one found on the photometric minor axis. Both colours of the shell are bluer with respect to the underlying galaxy. Blue shells may indicate the presence of a younger stellar population, not necessarily confined to the shells themselves. It may also be related to the presence of a significant metal-poor component within the shell. This hypothesis is supported by the positive correlation between metallicity and galaxy mass (Sandage & Visvanathan 1978) and the expected much smaller mass of the galaxy intruder with respect to the main galaxy. Finally, it could be that the shells contain less dust per unit stellar mass than the underlying galaxy or that the dust in the shells has different properties. It is suggestive that the [3.6]-[4.5] in the shell 24 resembles that in the central region of NGC 5982 where a faint excess emission is observed. Blue V-I colours have been found in other shell ellipticals, for example, in the South side of NGC 1344, which may be evidence for a relatively recent merger event (SCP07).
Stickel et al. (2004) claimed to find a large amount (
)
of
cold (
K) dust associated with atomic and molecular gas in the northern
shell of NGC 5128. We observe only some extended and extremely faint emission
at 160
m (surface brightness of
0.5 MJy sr-1 above the sky
background) at the position of shell 24. We do not find any excess
emission at 4.5, 5.8 or 8.0
m associated with this shell. We do find, however, a
counterpart in the mid-infrared for the HI gas cloud found by Morganti et al. (2006),
which is close to the shell 13 found by SCP07 (see Sect. 5.5).
The far-infrared emission in ellipticals is not entirely explained by steady state models of dust production by mass loss from the old population of red giant stars and dust destruction by grain sputtering in the hot interstellar gas (Temi et al. 2007). The excess of cold dust could have an external origin, from a recent merger with another galaxy. However, Temi et al. (2007) argue that the most plausible origin of the far IR emission is an internal mechanism, whereby cold dust originating from mass-loss from stars in the core is accumulated, and then buoyantly transported to distances of a few kpc by intermittent energy outbursts from AGNs.
More precise spectroscopic observations are required to verify or exclude the existence of a low-level AGN in the nucleus of NGC 5982, and to quantify the impact on the origin of the excess emission.
NGC 5982 contains warm gas (traced by H
,
[OIII]
5007) in
the central 4 kpc that lies in a rapidly rotating central component and in a
filament extending from the core to the Southern side (Sarzi et al. 2006).
The ionized gas in the centre rotates in the same sense as the stars and there is
a velocity gradient along the filament, with positive velocities in the Southern edge
and negative velocities in the innermost part.
Temi et al. (2007) found that optical emission line fluxes from warm
gas correlate weakly with L70 and L160 in ellipticals, suggesting that the dust
may be responsible for cooling the warm gas. They also suggest that the warm gas has
an internal origin, being expelled by orbiting red giant stars. Both results are,
however, not definitive due to the poor quality of the dataset. In NGC 5982
we find that the spatial distributions of the warm gas and the dust are different.
There is no dust emission along the filament of warm gas. In the central 2.5 kpc region,
the warm gas lies on some faint 70
m emission, and it encloses the irregularly
distributed patches traced by optical absorption (SCP07). There is also a small patch
of warm gas which lies on the southern peak of the 70
m emission. Dusty warm gas
forms possibly optically thick central disks in ellipticals (Lauer et al. 2005;
Sarzi et al. 2004).
The observations of SCP07 in a sample of shell ellipticals indicate that the dust
seems to be out of dynamical equilibrium, forming irregular patches and lanes. SCP07
conclude that shell ellipticals have a higher dust prevalence than normal
early type galaxies and that external influences are responsible for the ubiquitous
presence of dust.
No HI gas is found in the central region of NGC 5982. In the next section we describe in more detail the HI cloud found by (Morganti et al. 2006).
In the internal dust origin scenario of Temi et al. (2007),
cooling to below 104 K and consequent molecular formation are expected.
Conversely, Combes et al. (2006) suggest that the CO emission found in some ellipticals
has been accreted in many of them, and that it is the material fueling star formation.
No CO is detected in the central region of NGC 5982, with beam
resolutions of
and 12 arcsec (Combes et al. 2006).
We note, however, that the upper limit for the H2 gas of
(Combes et al. 2006) is approximately two orders of magnitude above
our estimated cold dust mass.
On the other hand, the CO/H2 ratio is very uncertain and the H2 content may be
underestimated by enhanced photodissociation of the CO molecule.
The central cold dust emission traced at 70 and 160
m
represents the bulk of dust in NGC 5982.
If the dust seen in Figs. 1g and h forms a disk inclined to the
principal axes of NGC 5982 then it would be expected to precess to the major
axis on a timescale of between 0.5 and 2 times the precession timescale
(Steiman-Cameron & Durisen 1990). The simulations of Varnas
(1990) suggest that this will happen in about
years. If the dynamical age (see Sect. 5.3) is greater than the
precession time then we would not expect to see the dust morphology found,
however both of these estimated ages have uncertainties of at least a factor two,
and we cannot conclude that they are inconsistent. The dust survival time
(i.e., sputtering lifetime) should be also consistent. It is expected to be low in
the centres of early type galaxies (
108 Gyr) although the self-shielding
in dust clouds may increase the estimates considerably (SCP07). If the sputtering
lifetime is much shorter than the merger age and the dust has a external
origin, then it was accreted in the late stages of the merger or in a more recent merger.
There is a small (
4 kpc in size) HI cloud located
6 kpc to the ENE
of the optical photometric centre of the galaxy, with a mass of
and at a velocity of 2830 km s-1, only
200 km s-1
lower than the systemic velocity (Morganti et al. 2006). The HI cloud is close to the
shell 13 found in SCP07, located 27.3
(i.e., 5.7 kpc) E.
We find emission at 3.6, 4.5, 5.8, 24 and 70
m which is very close to
the centre of the HI beam (see Fig. 8, top and bottom).
There is no ionized gas detected on the HI cloud although the ionized gas velocities
are consistent with the HI velocity (Morganti et al. 2006).
If the HI gas and cold dust observed in NGC 5982 have the same origin,
their very different spatial distributions observed
today suggest a severe disruption of these components.
We note that the dust mass estimated from the HI gas cloud is,
dependent upon the assumed gas-to-dust ratio, of a few times
,
a value of the same order as our estimated dust mass (see Sect. 4.3).
A search for CO emission at the position of the HI cloud would be valuable.
The very different distributions of the ISM components, the presence of a KDC (Wagner et al. 1988) and shells, and the proximity of other galaxies around NGC 5982 (García 1993) suggest a merger scenario to explain the shell formation and the origin of dust in this galaxy.
![]() |
Figure 8:
HST/ACS V-band image ( top). The centre of the image corresponds to the
position of the HI gas cloud detected with a beam of
|
| Open with DEXTER | |
Acknowledgements
We acknowledge the anonymous referee for his/her useful comments. This work is based on observations made with the Spitzer Space Telescope, which is operated by the Jet Propulsion Laboratory, CaliforniaInstitute of Technology, under a contract with NASA. Also it is based on observations made with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope Institute, which is operated by the Association of Universities for research in Astronomy, Inc., under NASA contract NAS 5-26555. These observations are associated with program GO9399.