A&A 444, 777-790 (2005)
DOI: 10.1051/0004-6361:20042193
L. Metcalfe1,2 - B. O'Halloran3,4 - B. McBreen3 - M. Delaney3 - M. Burgdorf5 - K. Leech2 - P. Barr6 - J. Clavel6 - D. Coia1,3 - L. Hanlon3 - P. Gallais7 - R. Laureijs6 - N. Smith8
1 - XMM-Newton Science Operations Centre, European Space Agency, Villafranca del Castillo, PO Box 50727, 28080 Madrid, Spain
2 -
ISO Data Centre, European Space Agency, Villafranca del Castillo, PO Box 50727, 28080 Madrid, Spain
3 -
Physics Department, University College Dublin, Dublin 4, Ireland
4 -
Dunsink Observatory, Castleknock, Dublin 15, Ireland
5 -
SIRTF Science Center, California Institute of Technology, 220-6, Pasadena, CA 91125, USA
6 -
Research and Scientific Support Department, European Space Agency, ESTEC, Keplerlaan 1, 2200 AG Noordwijk ZH, The Netherlands
7 -
CEA Saclay/Service d'Astrophysique, Orme des Merisiers, 91191 Gif-sur-Yvette Cedex, France
8 -
Cork Institute of Technology, Cork, Ireland
Received 17 October 2004 / Accepted 9 August 2005
Abstract
Markarian (Mkn) 297 is a complex system comprised of two
interacting galaxies that has been modelled with a variety of
scenarios. Observations of this system were made with the Infrared
Space Observatory (ISO) using the ISOCAM, ISOPHOT and LWS
instruments. ISOCAM maps at 6.7 m, 7.7
m,
12
m and 14.3
m are presented which, together with
PHT-S spectrometry of the central interacting region, probe the
dust obscured star formation and the properties of the organic
dust. The ISOCAM observations reveal that the strongest emission in
the four bands is at a location completely unremarkable at visible
and near-IR (e.g. 2MASS) wavelengths, and does not coincide with the
nuclear region of either colliding galaxy. This striking characteristic
has also been observed in the overlap
region of the colliding galaxies in the Antennae (NGC 4038/4039),
the intragroup region of Stephan's Quintet, and in IC 694 in the
interacting system Arp 299, and again underlines the importance
of infrared observations in understanding star formation in
colliding/merging systems. At 15
m, the hidden source in
Mkn 297 is, respectively, 14.6 and 3.8 times more luminous than the hidden
sources in the Antennae (NGC 4038/4039) and Stephan's Quintet.
Numerical simulations of the Mkn 297 system indicate that a co-planar
radial penetration between two disk galaxies yielded the observed
wing formation in the system about 1.5
108 years
after the collision. A complex emission pattern with knots and
ridges of emission was detected with ISOCAM. The 7.7
m map
predominantly shows the galaxy in emission from the 7.7
m
feature attributed to PAHs (Polycyclic Aromatic Hydrocarbons). The
14.3/7.7
m ratio is greater than unity over most of the
galaxy, implying widespread strong star formation. Strong
emission features were detected in the ISOPHOT spectrum, while
[O I], [O III] and [C II] emission lines were seen with LWS.
Using data from the three instruments, luminosities and masses for two
dust components were determined. The total infrared luminosity is
approximately 1011
,
which (marginally) classifies
the system as a luminous infrared galaxy (LIRG). A supernova that
exploded in 1979 (SN 1982aa) gave rise
to one of the most powerful known radio remnants which falls close to
the strongest mid-infrared source and is identified with star forming
region 14 in the optical. This supernova explosion may have been accompanied
by a gamma-ray burst (GRB), consistent with the idea that GRBs are
associated with supernovae in star forming regions, and a search
for a GRB consistent with the direction to Mkn 297, in satellite
data from July to December 1979, is recommended.
Key words: galaxies: general - galaxies: interactions - galaxies: starburst - ISM: dust, extinction - infrared: galaxies
Mkn 297 (also known as NGC 6052, Arp 209 and UGC 10182) has been the subject of many investigations because of its peculiar morphology (Weedman 1972). It has been referred to as a typical clumpy irregular galaxy (Taniguchi & Tamura 1987) and a giant HII region galaxy (Thuan & Martin 1981; Thuan 1983). It has a complex optical structure consisting of forty three identified bright optical knots within a common envelope (Hecquet et al. 1987). The galaxy has no obvious companion and appears to be isolated. However, it has two distinct components revealed through optical spectroscopy (Duflot 1976; Burenkov 1988), CO (Sage et al. 1993) and HI line emission (Sofue et al. 1990), and consequently the system has been variously modelled as a colliding system consisting of two late-type spiral galaxies (Alloin & Duflot 1979); the collision of a spiral galaxy with an irregular galaxy (Burenkov 1988); or more recently, using numerical N-body simulations, as the collision of two disk galaxies (Taniguchi & Noguchi 1991). Throughout this paper we will refer to the component identified as a spiral by Burenkov (1988) as galaxy A and the other component as galaxy B.
Table 1:
The log of the ISO observations of Mkn 297. The
nine columns list the TDT number (a unique identifier for an ISO
observation), the AOT number (which identifies the observing mode
used), the filter label, the wavelength range (
),
the reference wavelength of the filter where appropriate, the duration of
the measurement (including both the on- and off-source measurements), the
field-of-view of the instrument in the configuration employed, the positional
offset of the background measurement(s) with respect to the target reference
position, and any additional notes.
The two components have a velocity separation of about
200 km s-1 and modelling reveals that their closest
approach occurred about 1.5 108 years ago
(Taniguchi & Noguchi 1991). Mkn 297 has similar infrared properties to
the Antennae Galaxies (NGC 4038/4039) and, as in the
Antennae, the strongest infrared source does not coincide with the
nuclei of the interacting galaxies
(Vigroux et al. 1996; Mirabel et al. 1998; Kunze et al. 1996) nor, indeed, with any other feature
notable at visual or near-IR wavelengths (see Figs. 1-3 below).
Heeschen et al. (1983)
have found a peculiar compact and variable radio source, Mkn 297A,
lying several arcseconds north of the nuclei of the
interacting galaxies (see Fig. 8). They interpreted the compact source as a
complex of supernova remnants associated with a region of
vigorous star formation. Yin & Heeschen (1991) and Deeg et al. (1993) interpret Mkn 297A as a single luminous
radio supernova which exploded in July 1979. Mkn 297 was observed
by ROSAT, and has an X-ray luminosity
erg s-1 (Hensler et al. 1997) with ROSAT PSPC
around 1keV, a value that is comparable to the cumulative flux
of the Antennae observed with Chandra at similar energy (Zezas et al. 2002).
The capabilities of the Infrared Space Observatory (ISO) (Kessler et al. 1996) permitted observations of Mkn 297 across a wide range of IR wavelengths with three complementary instruments. The observations and data reduction techniques are presented in Sect. 2. The results are presented in Sect. 3 and discussed in Sect. 4. A value of H0 = 72 km s -1 Mpc-1 has been adopted and yields a distance of 66 Mpc to Mkn 297.
The ISO observations were obtained using: (a) the mid-infrared camera ISOCAM (Cesarsky et al. 1996), mainly in beam-switching mode, though one larger-area raster observation has been used; (b) the spectrometric channel (PHT-S) of the ISO photopolarimeter (ISOPHOT) (Lemke et al. 1996) in triangular chopping mode and (c) the medium-resolution grating mode of the long wavelength spectrometer LWS (Clegg et al. 1996) employing dedicated on- and off-source measurements. The parameters of these observations are listed in Table 1 and discussed in the following sections.
A detailed technical description of the ISOCAM instrument and its observing
modes can be found in (Blommaert et al. 2003). For the present
work the 1.5'' per pixel plate scale of the instrument has been
used yielding a field of view of 48''
48''. As
explained in Sect. 2.4 an available lower-resolution
6
/pixel map has been used only to rule out the
presence of very extended mid-IR emission from the system.
Observations were obtained in four filters, with reference
wavelengths (Moneti et al. 1997) at 6.7
m (LW2), 7.7
m
(LW6), 12
m (LW10) and 14.3
m (LW3). Each filter
observation consisted of two measurements, one centred 1.5
away from the target, and the other centred on the target
coordinates. Each of these beams was observed for 140 on-chip
integrations of 2.1 s each, for a total dwell-time per
filter of 588 s. The diameter of the ISO point spread function
(PSF) central maximum, in arcseconds, at the first Airy minimum is
0.84
m). The FWHM is about half that, and
Okumura (1998) obtained a value of 2.3'' at 6.7
m and
4.6'' at 14.3
m for the PSF FWHM in the 1.5'' per pixel
configuration.
All data processing used the CAM Interactive Analysis (CIA) software (Ott et al. 1997,2001; Delaney et al. 2002), as follows: (i) dark subtraction was performed using a dark model with correction for slow drift of the dark-current throughout the mission; (ii) glitch effects due to cosmic rays were removed following the method described in Starck et al. (1998a); (iii) transient correction for signal attenuation due to the lag in the detector response was performed by the method described by Abergel et al. (2000); (iv) the maps were flat-fielded using library flat-fields; (v) pixels affected by glitch residuals and other persistent effects were manually suppressed; (vi) the on-target and off-target beams were subtracted and (vii) the images were deconvolved using a multi-resolution transform method (Starck et al. 1998b).
Considerable care was needed in the preparation of
the 14.3/6.7 m ratio map. Apart from deconvolution
and background subtraction (O'Halloran et al. 2002), it was
necessary to eliminate any misalignment between the 14.3 and
6.7
m images due mainly to some jitter
in the
positioning of the CAM wheels, leading to small offsets in
astrometry. This was done by choosing a suitable reference source in both
images and locating its position using a Gaussian fit.
The coordinates of the reference object in each image were established
at sub-pixel accuracy by the fit and the vector distance between the
reference source positions in successive images was determined, establishing
the extent of misalignment. A PSF was
chosen, from a library of model PSFs, offset from the centre by
the negative of the distance vector. This PSF was used to
deconvolve and align the misaligned images. The reference image was then
deconvolved with a centred PSF so the resultant images were
aligned and had the same beam. Noise spikes in the ratio map, due to the
division of background pixels having noise-related low apparent signal
values were suppressed by applying the same positive threshold to the background
pixels in both images before division. Note that photometry of sources
was performed on the un-deconvolved maps using the CIA routine XPHOT, employing
aperture photometry scaled for the PSF.
The durations of the observations in Table 1 include
instrumental, but not spacecraft,
overheads. The flux densities quoted in Table 2 for four source regions
assume a spectrum having
,
and have a
photometric accuracy of about
15%.
Table 2:
The ISOCAM flux densities for MIR clumps in Mkn 297.
The photometry was performed before deconvolution and is
accurate to about 15%.
![]() |
Figure 1:
A deconvolved ISOCAM 14.3 ![]() ![]() ![]() |
Open with DEXTER |
PHT-S consists of a dual grating spectrometer with a resolving
power of 90 (Laureijs et al. 2003). Band SS covers the
range 2.5-4.9 m, while band SL covers
the range 5.8-11.6
m. The PHT-S spectrum of
the core region of the Mkn 297 system was obtained
by operating the
aperture of PHT-S in
triangular chopping mode using the ISOPHOT focal-plane chopper.
In this mode, the satellite pointed to the given target coordinates
centred between two off-source-positions, and the chopper moved
alternately from the target to the first off-source position, back
to the target and
then to the second off-source position. Consequently, 50% of the exposure
time is on target and 50% is distributed equally over the two
off-source positions.
In the on-target phase of the chop the PHT-S aperture was centred on the
centre of the ISOCAM pointing, a position about 3 arcsec
East of the centre of the 15
m image presented in Fig. 1, as indicated
in the figure. The
PHT-S aperture
covered completely ISOCAM source regions 1, 2 and 4, and vignetted
slightly the other regions. (See Sects. 2.4 and 3.3 below.)
PHT data processing was performed using the ISOPHOT Interactive Analysis (PIA) system, version 8.1 (Gabriel 2002). Data reduction consisted primarily of the removal of instrumental effects such as ionising particle impacts which result in a spurious increase in two or more consecutive read-out voltage values. The disturbance is usually very short and the slope of the integration ramp after a glitch is similar to the slope before it. Further reduction of the spectrum was performed by manual deglitching, where the slopes of the ramps were considered for a given PHT-S channel, and two possible corrections were applicable. First a strong startup transient was searched for (Acosta-Pulido et al. 2000) and, if present, suppressed by deleting either the initial 64 or 128 s of integration. The remaining signals were inspected for large deviations from the mean, indicating strong cosmic ray hits which could not be removed by the automatic deglitching. These anomalous signals were suppressed manually by masking data affected by the glitch. The mean differential signal (i.e. the difference of the on-source and off-source chopper beams) for each of the 128 pixels of PHT-S was derived. To absolutely calibrate the resulting spectrum it was scaled to match the signal level of the ISOPHOT-S Highly Processed Data Products for this target which can be found in the ISO archive (Richards & Klaas 2002) associated with the Mkn 297 archive entry for this observation. The HPDP are somewhat noisier than our PIA produced product, but benefit from being part of a systematic reprocessing of a large sample of PHT-S spectra, and we consider that they are likely to have a better overall calibration. Relative (from point to point in the spectrum) and absolute (overall normalisation of the spectrum to flux standards) spectrometric accuracies of better than 10% are indicated. The results were plotted to obtain the spectrum for the region of Mkn 297 in the PHT-S aperture. Flux values derived from this spectrum have to be further scaled, as described in Sect. 2.4 below, to yield values more appropriate to the galaxy as a whole.
A detailed technical description of the LWS instrument and its observing
modes can be found in (Gry et al. 2003). A LWS spectrum
of Mkn 297 spanning the range from 45 to 180 m was obtained
with the LWS aperture centred on the ISOCAM map, at a position about
3'' East of the nucleus of galaxy A, as indicated by
the filled square in Fig. 1. The beam of LWS was
slightly elliptical and its FWHM varied between 65'' and 85'',
depending on wavelength and direction (Swinyard et al. 1998). It was
assumed that the source was completely included in the beam of
the LWS instrument, so that no extended-source correction was necessary.
Because of an expected significant background signal in the LWS aperture due to combined Zodiacal and Cirrus emission in this sky region, a dedicated background measurement was made with the LWS aperture centred about 12 arcmin away from the on-source measurement (at a position determined to be of similar background based on the IRAS maps), and the resulting background was subtracted from the on-target case.
The grating was scanned 6 times over the entire wavelength range. The spectral sampling was set to give 4 samples per resolution element in each of the scans.
We made use of the LWS Highly Processed Data Products (HPDP) for this target, available from the ISO archive (Lloyd et al. 2003). Because of the small flux from the source and the resulting poor signal-to-noise ratio of the LWS continuum spectrum, the data were rebinned to one point for each of the 10 LWS detectors, employing a scan-averaging method described by Sidher et al. (2000), yielding 10 samples of the continuum spanning the LWS spectral range and plotted in Fig. 7.
The rather different apertures, or fields-of-view, of the
instruments warrant comment. Observatory pipeline products from the
low-resolution ISOCAM 6
/pixel raster observation
of ISO TDT 09101068 have been used to establish that substantially less
than 10% of the total galaxy MIR light falls outside of the
1.5
/pixel, 48
field-of-view
upon which the results reported here
rely. The much larger LWS aperture (see Sect. 2.3 above) is
assumed to have covered the full extent of the galaxy, even at the
longer wavelengths addressed by LWS. The PHT-S aperture is a square
24
on a side, and was placed as shown in
Fig. 1. It certainly missed some of the integrated galaxy signal in the
PHT-S wavelength range, although it included most of the knots of
mid-infrared emission found in the ISOCAM map. But the PHT-S aperture
has a vignetted profile (Laureijs et al. 2003). So some scaling has
to be applied to PHT-S signals to adapt them to the signals that
would have been recorded in corresponding spectral ranges over the full
area of the ISOCAM footprint. To determine the scaling factor the
PHT-S spectrum was used to
synthesise a 7.7
m (LW6) ISOCAM bandpass measurement by
integrating the CAM filter profile over the PHT-S spectrum. The
result was compared with the CAM 7.7 micron global galaxy
photometry and found to differ by a factor of 1.46. PAH feature
strengths derived from the PHT-S spectrum were scaled by this factor
(Sect. 3.3).
The deconvolved 14.3 m map of Mkn 297 overlaid on a HST
image (filter F555W at 5407 Å) is presented in Fig. 1, while
the deconvolved maps at 6.7, 7.7, 12 and 14.3
m are
presented in Figs. 2i-2iv. The figures reveal a complex
emission pattern of overlapping sources and seven source regions
are numbered on the maps. The integrated CAM flux densities for the system
are presented in Table 2 along with estimated flux densities from
the four brightest MIR clumps. Knots that appear well defined in one bandpass
may not be so sharp in another (see Fig. 2). Representative flux-density values have
been derived by performing aperture photometry about the locations
numbered in 1 to 4 in Figs. 1 and 2 and scaling each result to the PSF that would
deliver the same signal in the aperture used. Where no well-defined
point-source appears in some bandpass these "point-source'' flux
densities nevertheless serve to characterise relative brightness
over the map. Source 4
coincides approximately with the nuclear region of galaxy A and
source 2 with the nuclear region of galaxy B. The strongest
source, labeled 1 on all the CAM maps, does not coincide with
the nucleus of either galaxy but is displaced by about 2 kpc
from the nucleus of galaxy A and peaks over a region devoid of
optical emission knots. (We stress, by the way, that this source is clearly
distinguished also in the un-deconvolved maps.) Hecquet et al. (1987) reported numerous optical
emission knots in Mkn 297, and knots 7 and 8 in their numbering
scheme have the second and third highest values of B-R (1.8 and 2.3 respectively)
for the system and are within a few
arcseconds of the mid-infrared peak. Knot 14 is the site of the variable radio
source reported by Heeschen et al. (1983). These optical knots may represent the visible
portions of a much vaster obscured starburst. A similar
result has been obtained
(Vigroux et al. 1996; Mirabel et al. 1998) for the interacting galaxies
NGC 4038/4039 (The Antennae) where the strongest mid-infrared source
does not coincide with either of the nuclei of the
constituent galaxies.
![]() |
Figure 2:
Deconvolved ISOCAM contour maps overlaid on the HST image:
i) 6.7 ![]() ![]() ![]() ![]() |
Open with DEXTER |
In order to compare the ISOCAM map with previously published
results at other wavelengths
(Lonsdale et al. 1992; Heeschen et al. 1983; Deeg et al. 1993,1997; Hecquet et al. 1987) it is
necessary to consider some inconsistencies in the published
astrometry of features in the Mkn 297 system. The first three of the
above-referenced papers use astrometry reported by Heeschen et al. (1983),
derived ultimately from glass
copies of original photographic plates and attributed in
Hecquet et al. (1987)
to a private communication from Casini in 1980. However, a careful comparison
of the Hecquet et al. (1987) astrometry (their Table 1 and Fig. 2b) with
2MASS and HST images reveals a roughly 6
offset of their coordinates to the North West.
The resulting positions for optical knots in Mkn 297 differ
systematically from positions derived by
Deeg et al. (1993) based on independent CCD photometry. Similarly,
the radio supernova reported by Heeschen et al. (1983)
and Yin (1994), which from Figs. 1, 2 and 8 of this paper
clearly corresponds to H
source 14 of Hecquet et al. (1987), only falls at that location subject to the
above adjustment of the astrometry of Hecquet et al. (1987). The astrometry of Casini, Heeschen et al. (1983), Yin (1994) and Deeg et al. (1993)
appears to be mutually consistent with all recent imaging data.
In fact, Lonsdale et al. (1992) noted the same discrepancy and
refer to a private communication from Hecquet & Coupinot (1991)
in which they concur.
![]() |
Figure 3:
An overlay of the ISO 15 ![]() |
Open with DEXTER |
Figure 1 therefore adopts
the 2MASS astrometry referenced to ISOCAM source 4, assumed to
correspond to the brightest knot in the 2MASS K-band data. This
places ISO source 2 over the nucleus of galaxy B. Figure 3 presents an
overlay of the ISO 15 m (white contours) and 2MASS K-band (black contours)
on the HST image of Mkn 297 used in Fig. 1. It is clearly seen that
ISO source 1 is not evident in the 2MASS map. (See also the J, H and K
maps published by Smith et al. (1996).)
The K-band data of Cairós et al. (2003) also show a source
roughly 4 arcsec North of our source 1 and consistent with the
location of the 2MASS source and the optical knots.
Within the accuracy of the various sources of positional information used to construct Fig. 1 these several detections are consistent with a single extended emitting region suffering, at different positions, different degrees of extinction.
Figure 2 presents the 6.7 m, 7.7
m, 12
m and 14.3
m
ISOCAM maps. The 7.7
m map essentially shows the galaxy in the
emission of a spectral feature usually attributed to polycyclic
aromatic hydrocarbon (PAH) emission. The 14.3
m map includes contributions
from warm dust, nebular line emission such as [NeII] and PAH
feature emission. The galaxy has the
approximate shape of an inverted tuning fork with regions 1, 2, 3
and 4 dominating the emission to different degrees at different
wavelengths. IRAS flux densities for the system are given in
Table 3 for reference.
Table 3: IRAS flux densities for Mkn 297 (from Sanders et al. 2003). The four columns list the filter used, the wavelength range, the reference wavelength for the filter and the flux density measured.
![]() |
Figure 4:
The 14.3/6.7 um ratio map overlaid on the HST image of Mkn 297. Deconvolved
6.7 and 14.3 ![]() ![]() |
Open with DEXTER |
The 14.3/6.7 m ratio map presented in Fig. 4 shows values
greater than unity across much of the map, which is indicative of
active star formation over a wide area. The 14.3/6.7
m ratio
generally decreases as interactions develop and starbursts age because PAHs,
which dominate the shorter-wavelength band,
are no longer destroyed by the highly ionizing O-stars, and dust
emission declines (Vigroux et al. 1999; Charmandaris et al. 1999; Helou 1999; Cesarsky & Sauvage 1999).
The ratio is highest for source regions 1, 2, 4 and 3 respectively,
with a peak value of 3.4 at source 1 and a global value of 1.5 to 2
over much of the galaxy.
The PHT-S spectrum is presented in Fig. 5, and line-fluxes derived for the features in the spectrum are given in Table 4. To allow for the limited (vignetted) aperture size of PHT-S with respect to the extent of the galaxy (see Fig. 1), the PHT-S signal was normalised to the ISOCAM 7.7 micron flux density as described in Sect. 2.4.
![]() |
Figure 5: PHT-SL spectrum of Mkn 297, with the PHT-S aperture centred 3'' East of the nucleus of galaxy A. The PAH features are indicated. The spectrum has been Hanning smoothed. The flux levels are those recorded in the PHT-S aperture, before scaling to the larger ISOCAM field-of-view (see Sect. 2.4). |
Open with DEXTER |
The most marked characteristic of the PHT-S spectrum is the
presence of the family of infrared bands at 6.2, 7.7, 8.6 and
11.3 m that are generally attributed to emission from
PAHs. There is a possible feature, seen at about the 3-sigma level, corresponding to
the nebular line emission from [Ar III] at a rest-frame wavelength 8.99
m. There is
no convincing evidence of [S IV] at 10.51
m. These features have been
observed in other starburst galaxies (Metcalfe et al. 1996; O'Halloran et al. 2000; Laureijs et al. 2000; O'Halloran et al. 2002; Helou 1999; O'Halloran et al. 2005).
Table 4: PHT-SL fluxes for Mkn 297. The columns give, respectively, the line identification, the wavelength range, and the integrated line-flux and the integrated line-flux scaled to the full ISOCAM footprint, as described in Sect. 2.4.
The LWS line fluxes for the far-infrared lines at 157.7 m [CII],
88.8
m [O II] and for [OI] at 63.1
m and
146.7
m are listed in Table 5 and corresponding spectra are
presented in Fig. 6.
Table 6 lists the continuum source strength recorded in the ten LWS detectors by binning data across the spectral range of each detector. The spectral energy distribution of Mkn 297 using ISOCAM, PHT-SL, LWS and IRAS fluxes is presented in Fig. 7.
Table 5: LWS line fluxes for Mkn 297. The three columns give the line identification, the wavelength and the line-flux respectively.
![]() |
Figure 6:
LWS spectra of a) the [OI] line at 63.1 ![]() ![]() ![]() ![]() |
Open with DEXTER |
Taniguchi & Noguchi (1991), using numerical N body
simulations, were able to account for the morphology and velocity
structure of Mkn 297. The simulations suggested that a coplanar
radial penetration between two disk galaxies yielded the wing
formation which is seen now, about 1.5 108 years after
the collision, when the disk of an edge-on (from the Earth's perspective)
galaxy has been deformed into a wing structure. ISOCAM source 2 (Fig. 1)
corresponds to the nucleus of the edge-on galaxy, and the North-South elongated
structure corresponds to poorly
formed tails emerging from the nucleus of this galaxy (Taniguchi & Noguchi 1991).
One of the best examples of a galaxy exhibiting such a tail morphology
is the Antennae (NGC 4038/4039), where the two well developed
tails are considered tidal debris from each component galaxy. The tail of
Mkn 297 differs from the Antennae because it is not so well developed
and is more luminous in the infrared, with a high 14.3/6.7
m
ratio indicative of massive star formation.
Table 6: Flux densities measured with LWS, binned to yield a single number per LWS detector.
![]() |
Figure 7:
Spectral energy distribution of Mkn 297 using
ISO and IRAS flux data, including the key for the different
symbols. The PHT-S spectrum has been normalised to the CAM
7.7 ![]() |
Open with DEXTER |
The future evolution of the two galaxies has been predicted based on the simulations. After the wing phase, the stellar disk is deformed into a bar-like structure, gas clouds tend to sink into the nuclear regions of the edge-on galaxy and the nuclear starburst becomes dominant. It is expected that Mkn 297 will evolve into two galaxies, one a starburst nucleus galaxy and the other a ring galaxy. (See also Efstathiou et al. 2000.)
The simulations also suggest a separation between the two galaxies of about four galactic radii, and hence the orbital plane of the collision must be significantly inclined from the plane of the sky (by about 70 degrees, to account for the apparent closeness of the two galactic nuclei).
Mkn 297 has an optical spectrum typical of a HII region, with
line intensities consistent with excitation by massive young
stars. The nuclear region of galaxy A (source 4 in Fig. 1) has
the highest optical surface brightness, with a strong stellar
continuum and an abundance of [O/H] 8.7, which is
near the solar value and typical of late-type spiral galaxies
(Lonsdale et al. 1992). The nuclear region of galaxy B (source 2 in
Fig. 1) has a much bluer spectrum, with an excitation
comparable to that of the HII regions in the Large
Magellanic Cloud (Burenkov 1988).
Compact CO line emission from Mkn 297 has been mapped and, at
the resolution of the systems used (of the order of 20
HPBW) has
been associated with the optical knots (Sage et al. 1993; Sofue et al. 1990). The
estimated mass of H2 is
and is approximately 75% of the HI mass of
.
Atomic oxygen and ionized carbon are the principal coolants of the
gaseous interstellar medium via their fine structure lines in the far
infrared (FIR). In particular, [OI] at 63 m and [CII] at
158
m dominate the cooling in the photodissociation regions
associated with massive stars such as Wolf Rayets, along with [OI] at
146
m and [OIII] at 88
m (Sanders et al. 2003). The [OI]
and [CII] features are also produced in the warm atomic gas behind
dissociative shocks, in HII regions or in photodissociation regions
(PDRs), while [OIII] is more associated with denser environments
within HII regions (Braine & Hughes 1999; Malhotra et al. 1997).
The [OI] 63 m and [CII] 158
m lines were well detected in Mkn 297
and were used to determine the luminosity ratios
/
and (
)/
,
which probe the nature of the environment within
the galaxy. Values of
/
and (
)/
were determined using the IRAS flux-densities (Malhotra et al. 1997).
The values found are
consistent with those from other starburst galaxies
(Braine & Hughes 1999; Malhotra et al. 1997), though for higher dust
temperatures and star formation rates these ratios decrease (Malhotra et al. 1997).
Mkn 297 falls at the higher star
formation/dust temperature end of this
correlation, which provides further evidence of strong ongoing star
formation within this galaxy.
The [CII] line emission is the dominant gas coolant in most regions
of atomic interstellar gas and therefore reflects the heating input
to the gas (Malhotra et al. 2001). The ratio of [CII]/far-infrared
(FIR), as a function of the ratio of the flux density at 60 mm to
100
m, R(60/100), has been measured for a large sample of galaxies
and used to study the radiation field (Helou et al. 2001).
The
[CII]/FIR ratio for Mkn 297 is
and R(60/100) is 0.64 and
these values are typical of a star forming galaxy as shown in
Fig. 1 in Helou et al. (2001). A strong correlation was also
found between [CII] and the integrated mid-infrared flux in the range 5-10
m
in a large sample of star forming galaxies. The mid-infrared flux
is dominated by aromatic features (AFEs), that are generally
associated with the smallest interstellar grains. The ratio of the
two quantities [CII]/AFE, where AFE is
over the wavelength
range 5-10
m, is nearly constant at 1.5% over a wide range in
values of R (60/100) (Helou et al. 2001; Dale et al. 2000).
In Mkn 297, the value of
was
synthesised from the PHT-SL spectrum (Fig. 5), extrapolating the
continuum to 5 microns in the integral. The resulting value was scaled to the
area of the ISOCAM footprint, as described in Sect. 2.4, resulting in a high value
of [CII]/AFE, 4%.
Although the value of 4% might be regarded as an upper limit because the PHT-SL
measurement covered a much smaller field of view than the LWS spectrometer
(Table 1). However, as argued in Sect. 2.4, it is unlikely that much flux in
the PHT-SL range falls outside of the ISOCAM footprint to which it has been scaled.
The value of 4% is an additional argument in favour of star formation
dominating the mid-infrared emission.
The unidentified infrared bands (UIBs) dominate the mid-infrared
emission from Mkn 297 (Fig. 5) with little or no emission in the
wavelength region between 9 and 10 m from very small grains
(VSGs). This result is also observed in IC 694 (Source A in Arp
299) but not in knot A in the Antennae. In cases where the
mid-infrared emission is dominated by UIBs, the SFR can be obtained
from the following relationship (Roussel et al. 2001):
The SFR can also be obtained from the H
luminosity
using (Lee et al. 2002):
A dust model for Mkn 297 is shown in Fig. 7, denoted by the solid curve. It includes two separate dust populations: a warm dust component at 130 K and a cooler dust component at 38 K (Krugel & Siebenmorgen 1994; Siebenmorgen et al. 1999).
There is growing evidence for the existence of several components in
the dust distribution of galaxies (Popescu et al. 2002; Klaas et al. 2001) that can be broken
down into warm dust components with K associated with
star formation regions, and a spatially extended distribution of cold
dust with
K.
Very small dust grains and PAHs are transiently heated by the
single-photon emission process to temperatures much higher than
40 K, up to several hundred degrees Kelvin (Désert et al. 1990), out of
thermal equilibrium with their environment (Calzetti et al. 2000). These
account for the bulk of the emission in the PHT-S waveband and for
emission at wavelengths shorter than 40 m. The "warmer'' or
"hot dust'' component, the 130 K component in Fig. 7, is due to small grains.
Dust emission from HII regions is dominated by large dust grains heated by the intense radiation field. These grains can reach temperatures in excess of 20 K up to more than 40 K. Due to the nature of the blackbody, this emission will outshine the cold dust in these regions. The 38 K dust component in Fig. 7 comes from large grains in HII regions.
The coldest dust emission is associated with large
grains emitting at wavelengths in excess of 80 m
and in thermal equilibrium with their environment
(Calzetti et al. 2000). A small contribution at 100
m is also made
by very small dust grains. The large dust grains account
for practically all the emission longward of 80
m in
galaxies, with the maximum of the SED occurring close to
200
m. We have not detected this very cold component in Mkn 297
within the sensitivity of the current measurements (and given
the occurrence of cirrus and zodiacal light confusion at the
longest wavelengths).
The luminosity of each component was found by integrating over a
model greybody, fitting temperatures and calculating a
dust mass for the component (Klaas & Elsaesser 1991). Using
this method, the IR luminosities and dust component masses in
similar starburst galaxies such as the Antennae and
NGC 6240 (Klaas et al. 1997) were calculated in order to check
the validity of the derived values. The values obtained were
generally in good agreement
The dust model yields for Mkn 297 a total
of
.
For each component, the luminosities and dust masses are given in
Table 7. The warmer dust emission is concentrated within
the disk and associated with the strongly ionizing sources found
in star forming regions.
Table 7:
Total infrared luminosity and dust component masses for Mkn 297. The total FIR
luminosity of
places this galaxy right on the
threshold of LIRG status.
Deeg et al. (1997) report radio data for a source within 3
of
Source 1 having flux densities of 244
20 mJy and
104
7 mJy at 0.325 GHz and 1.489 GHz respectively.
This is consistent with the variable radio source reported by
Heeschen et al. (1983), which varied at 4.885 GHz (6 cm) from 4.2 mJy on January 12 1980
to 12.2 mJy on April 23 of that year. See also Yin et al. (1994), and Fig. 8.
Pierini et al. (2003) give expressions to relate warm and cold FIR emission to radio flux, and the values for FIR emission for Mkn 297 given in Table 7 fall on the correlations in each case. From Condon et al. (1990), log (L(1.4 Ghz)) = 22.74, while Table 7 yields log (L(FIR)) = 37.85. For the warm and cold dust respectively: log (L(FIR)) = 36.74 and log (L(FIR)) = 37.82. All lie close to the relevant correlations, though it must be remarked that the total FIR luminosity for Mkn 297, and particularly the cold component, fall towards the high FIR extreme of the Pierini et al. correlations.
![]() |
Figure 8: An overlay of the Yin (1994) 6 cm radio contours on the HST image of Mkn 297. |
Open with DEXTER |
The JHK 2MASS images of Mkn 297 reveal a source that is North of
galaxy A and coincides with the star forming regions 7 and 8 of
Hecquet et al. (1987). An overlay of the ISO 15 m map and
the 2MASS K-band map on the HST image is shown in Fig. 3. The
K-band peak clearly corresponds to the optically evident HII region,
and not at all to the hidden MIR source - Source 1. The star
forming region that is nearest to the hidden source is region 16
of Hecquet et al. (1987) and has the largest extinction of all
the star forming regions in Mkn 297. The hidden source in Mkn 297
is inconspicuous in the optical and also in the JHK 2MASS images.
There are three other sources where strong mid-infrared emission has been detected that are unremarkable in the optical. The first source is knot A in the overlap region of the colliding galaxies in the Antennae (Vigroux et al. 1996; Mirabel et al. 1998), the second is in the intragroup region of Stephan's Quintet that could be an intruding galaxy (Xu et al. 1999) and the third is source A or IC694 in the interacting galaxies in Arp 299 (Charmandaris et al. 2002; Gallais et al. 2004). The hidden source in Mkn 297 is important in this context and is very similar to those in the Antennae and Stephan's Quintet because the hidden sources do not coincide with the nuclei of any of the interacting galaxies.
At 15 m (LW3), the hidden source in Mkn 297 is 14.6 times more
luminous than knot A in the Antennae, and 3.8 times more luminous than
the hidden source in Stephan's Quintet. However, its luminosity is
only 0.09 times that of the hidden source IC 694 in Arp 299. However
the ratio of the LW3 flux from the hidden source to the total flux is
similar for all sources with values of 0.23, 0.19, 0.26 and 0.10
for Mkn 297, the Antennae, Arp 299 and Stephan's Quintet
respectively.
The hidden source in Mkn 297 appears to be more similar in properties to knot A in the Antennae, than it is to the nuclear source IC 694 (Charmandaris et al. 2002; Gallais et al. 2004). The source in IC 694 becomes visible as a point source in the JHK 2MASS images and gets brighter from J to H to K and also contributes significantly to the near infrared radiation from the galaxy. However this is not the case in Mkn 297 or the Antennae, where knot A is inconspicuous in the near infrared and becomes dominant in the mid-infrared. Table 8 lists comparative infrared properties for Mkn 297 and three sources with somewhat similar morphological characteristics.
Table 8:
A comparison of the properties the hidden sources (h.s.) in four systems.
Column 3 lists the total 15 m emission from the host galaxy in each case.
The differences between the hidden sources suggest that dust
extinction is lower in the hidden source in IC 694 than in Mkn 297 or
knot A in the Antennae. An alternative interpretation is that the
thermal emission from an enshrouded AGN becomes visible in the near
infrared. This phenomenon was observed in the nucleus of NGC 1068
(Alonso-Herrero et al. 2003) and will be discussed further in
Sect. 4.6. The variable and compact radio source Mkn 297A (Fig. 8) that
falls on the edge of the hidden source in Fig. 1 has been studied
extensively (Hummel et al. 1987; Condon et al. 1991; Yin & Heeschen 1991; Yin 1994). This source
was unresolved with the resolutions of the VLA and the flux density
shows a rapid rise followed by a slower decline with maximum flux
density occurring first at the shortest wavelength (Yin 1994). The
flux density, S, was modelled as a power law function of both
frequency and time of the form
The VLA map of Mkn 297 (Fig. 8) shows that the variable radio source coincides with star forming region 14 of Hecquet et al. (1987).
Several diagnostic tools have been proposed to probe the nature of the
activity within a central starburst region. Lu et al. (1999)
cite the ratio of the integrated PAH luminosity and the 40 to 120 m
IR luminosity as providing a tool to discriminate between
starbursts, AGN and normal galaxies. The lower the ratio the
more likely the galaxy harbours an AGN, due to the dominance of very
small dust grain emission powered by the AGN in the 40 to 120
m
region (Vigroux et al. 1999). For Mkn 297, the ratio
is 0.11, which is typical of a strong starburst. Similarly,
Lutz et al. (1998), Genzel et al. (1998),
Laureijs et al. (2000)
state that the ratio of the 7.7
m PAH flux to the nearby continuum can
provide a measure of the level of star formation activity within
the nucleus, given the ubiquity of strong PAH features in regions of
high star formation (Clavel et al. 2000).
The ratio for Mkn 297 is 3.7. Plotting this value, and the
5.9
m/60
m flux densities
(0.061 Jy in the PHT-S aperture/6.79Jy from IRAS) in Fig. 3 of
Laureijs et al. (2000) indicates strong star formation,
suggesting that an AGN within Mkn 297 is unlikely (O'Halloran et al. 2000;
Laureijs et al. 2000).
Another set of diagnostics uses empirical criteria based on the fact that mid-infrared emission from star forming or active galaxies arises mostly from HII regions, photo-dissociation regions (PDRs) and AGNs (Laurent et al. 2000). The diagnostic diagrams provide quantitative estimates of the contributions of AGN, PDR and HII regions. These diagrams are referred to as the Laurent diagnostic diagrams (Peeters et al. 2004). In Mkn 297 the ratio of LW3/LW2 is 2 and that of LW2/LW4 is roughly 2, where the LW4 value was synthesised from the PHT-SL spectrum, making approximate allowance for the fact that the LW4 bandpass extends to shorter wavelength than PHT-S. These values show that the dominant contribution to the mid-infrared emission is from PDR and HII regions in Mkn 297 and not from an AGN.
However, mid-IR diagnostics cannot probe beyond values of about 10 magnitudes of extinction, and so miss any AGN embedded deep into a star forming region with, say, 100 magnitudes of extinction. The absence of an AGN within Mkn 297 can therefore only be tentatively suggested until it is confirmed by hard X-ray observations.
Heeschen et al. (1983) found a peculiar, compact and variable radio source which falls about 2 kpc north of the nucleus of galaxy A, on the edge of the strongest infrared source in Fig. 1 (source 1). They suggested, among other possibilities, that it might be a complex of supernova remnants. In subsequent radio VLA observations (Yin 1994; Yin & Heeschen 1991) a decay in flux density from 14 mJy to 4.7 mJy at 20 cm was recorded between 1983 and 1990, and the source (Fig. 8) was attributed to a single, very luminous, radio supernova that exploded in mid-1979. The radio SN was later cataloged as SN 1982aa, of uncertain type because of the paucity of radio results within its first 8 years and the absence of optical observations (Green 1994). VLBI observations (Lonsdale et al. 1992) are consistent with the single supernova occurring in a molecular cloud.
The most surprising aspect of the supernova interpretation was the very high
radio luminosity of
.
As
a radio supernova it is one of the most energetic events ever
observed, alongside SN 1986 J, 41. 9 + 58 in M82, SN 1998bw and
SN 2003dh (Wilkinson & de Bruyn 1990; Berger et al. 2003; Galama et al. 1998). These powerful supernovae
are more energetic, by a factor of 10, than canonical supernovae which
release
1052 erg during the explosion. Indeed, SN 1982aa
has a total released energy of 10
erg, as derived
from the radio flux reported by (Yin 1994).
So far, many of the most powerful supernovae have been of type Ib/c, the best known example of which has been SN 1998bw, one of the most luminous radio SN ever observed (Weiler et al. 2002,2001). Interestingly, SN 1998bw has been strongly associated with a weak gamma-ray burst (GRB), GRB 980425 (Soderberg et al. 2004; Kouveliotou et al. 2004; Galama et al. 1998). It has only been exceeded in radio luminosity by SN 1982aa (Yin 1994) and SN 2003dh, which is identified with the cosmological burst GRB 030329 (Berger et al. 2003).
A well known model for the progenitors of GRBs, linked to powerful supernovae, and by extension to regions of massive star formation, is the "collapsar'' model. In this model the core of a massive star collapses to a black hole, which then powers the GRB by the accretion of an additional solar mass (McFadyean & Woosley 1999), and signatures consistent with spin-up and spin-down of black holes have been found in the light profiles of GRBs (McBreen et al. 2002). For SN 1982aa, given the amount of energy released it is quite plausible that a GRB may have been associated with this candidate obscured supernova. Furthermore, Mkn 297 is morphologically similar to host galaxies of GRBs seen in HST images, which show most GRB host galaxies to be either interacting or disturbed (Wainwright 2005). In 1979, many spacecraft had GRB detectors and a search is recommended in data from the period July through December 1979, for a GRB consistent with the direction to Mkn 297. It is also important to continue radio observations of Mkn 297A so that a proper comparison can be made of its long term radio profile with that of SN 1998bw, SN 2003dh and other supernovae (Frail et al. 2000) with associated GRBs that have yet to be discovered.
Acknowledgements
We would like to thank the referee, Dr. Vassilis Charmandaris, for detailed comments which greatly improved the paper. We gratefully thank Dr. Qi Feng Yin for providing us with his VLA radio maps for study and for the radio data reproduced in Fig. 8. We thank Dr. Avishay Gal-Yam for helpful comments ans suggestions. The ISOCAM data presented in this paper were analysed using CIA, a joint development by the ESA Astrophysics Division and the ISOCAM Consortium. The ISOCAM Consortium is led by the ISOCAM PI, C. Cesarsky. The ISOPHOT data presented in this paper was reduced using PIA, which is a joint development by the ESA Astrophysics Division and the ISOPHOT consortium. LIA is a joint development of the ISO-LWS Instrument Team at Rutherford Appleton Laboratories and the Infrared Processing and Analysis Center (IPAC). The research of M. Burgdorf was carried out at the Jet Propulsion Laboratory, California Institute of Technology, and sponsored by the National Aeronautics and Space Administration. This publication makes use of data products from the Two Micron All Sky Survey, which is a joint project of the University of Massachusetts and the Infrared Processing and Analysis Center/California Institute of Technology, funded by the National Aeronautics and Space Administration and the National Science Foundation.