A&A 431, 433-449 (2005)
DOI: 10.1051/0004-6361:20041782
D. Coia 1 - B. McBreen1 - L. Metcalfe2, 3 - A. Biviano4 - B. Altieri 2 - S. Ott5 - B. Fort6 - J.-P. Kneib7,8 - Y. Mellier6,9 - M.-A. Miville-Deschênes10 - B. O'Halloran1,11 - C. Sanchez-Fernandez2
1 - Department of Experimental Physics, University College, Belfield, Dublin 4, Ireland
2 - XMM-Newton Science Operations Centre, European Space Agency, Villafranca del Castillo, PO Box 50727,
28080 Madrid, Spain
3 - ISO Data Centre, European Space Agency, Villafranca del Castillo, PO Box 50727, 28080 Madrid, Spain
4 - INAF/Osservatorio Astronomico di Trieste, via G.B. Tiepolo 11, 34131, Trieste, Italy
5 - Science Operations and Data Systems Division of ESA, ESTEC, Keplerlaan 1, 2200 AG Noordwijk,
The Netherlands
6 - Institut d'Astrophysique de Paris, 98 bis boulevard Arago, 75014 Paris, France
7 - Observatoire Midi-Pyrénées, 14 avenue Edouard Belin, 31400 Toulouse, France
8 - California Institute of Technology, Pasadena, CA 91125, USA
9 - Observatoire de Paris, 61 avenue de l'Observatoire, 75014 Paris, France
10 - Canadian Institute for Theoretical Astrophysics, 60 St-George Street, Toronto, Ontario M5S 3H8, Canada
11 - Dunsink Observatory, Castleknock, Dublin 15, Ireland
Received 3 August 2004 / Accepted 26 October 2004
Abstract
Observations of the core of the massive cluster Cl 0024+1654, at a redshift
,
were obtained with the Infrared Space Observatory using ISOCAM at 6.7
m (hereafter 7
m) and
14.3
m (hereafter 15
m). Thirty five sources were detected at 15
m and thirteen of them are
spectroscopically identified with cluster galaxies. The remaining sources consist of four stars, one quasar, one
foreground galaxy, three background galaxies and thirteen sources with unknown redshift. The sources with
unknown redshift are all likely to be background sources that are gravitationally lensed by the cluster.
The spectral energy distributions (SEDs) of twelve cluster galaxies were fit from a selection of
20 models using the program GRASIL. The ISOCAM sources have best-fit SEDs typical of spiral or starburst models
observed 1 Gyr after the main starburst event. The star formation rates were obtained for cluster members. The
median infrared luminosity of the twelve cluster galaxies is
,
with 10 having
infrared luminosity above
,
and so lying near or above the
threshold for identification as a luminous infrared galaxy (LIRG). The [O II] star formation rates obtained for 3 cluster galaxies are one to two orders of magnitude lower than the infrared values, implying that most of the
star formation is missed in the optical because it is enshrouded by dust in the starburst galaxy.
The cluster galaxies in general are spatially more concentrated than those
detected at 15 m. However the velocity distributions of the two categories are comparable. The
colour-magnitude diagramme is given for the galaxies within the ISOCAM map. Only 20% of the galaxies that are
significantly bluer than the cluster main sequence were detected at 15
m, to the limiting sensitivity
recorded. The counterparts of about half of the 15
m cluster sources are blue, luminous, star-forming
systems and the type of galaxy that is usually associated with the Butcher-Oemler effect. HST images of these
galaxies reveal a disturbed morphology with a tendency for an absence of nearby companions. Surprisingly the
counterparts of the remaining 15
m cluster galaxies lie on the main sequence of the colour-magnitude
diagramme. However in HST images they all have nearby companions and appear to be involved in interactions and
mergers. Dust obscuration may be a major cause of the 15
m sources appearing on the cluster main sequence.
The majority of the ISOCAM sources in the Butcher-Oemler region of the colour-magnitude diagram are best fit by
spiral-type SEDs whereas post-starburst models are preferred on the main sequence, with the
starburst event probably triggered by interaction with one or more galaxies.
Finally, the mid-infrared results on Cl 0024+1654 are compared with four other clusters observed with ISOCAM. Scaling the LIRG count in Cl 0024+1654 to the clusters Abell 370, Abell 1689, Abell 2218 and Abell 2390 with reference to their virial radii, masses, distances, and the sky area scanned in each case, we compared the number of LIRGs observed in each cluster. The number in Abell 370 is smaller than expected by about an order of magnitude, even though the two clusters are very similar in mass, redshift and optical richness. The number of LIRGs detected in each of Abell 1689, Abell 2218 and Abell 2390 is 0, whereas 3 were expected from the comparison with Cl 0024+1654. A comparison of the mid-infrared sources in Abell 1689 and Abell 2218 shows that the sources in Abell 1689 are more luminous and follow the same trend identified in the comparison between Cl 0024+1654 and Abell 370. These trends seem to be related to the dynamical status and history of the clusters.
Key words: galaxies: clusters: general - galaxies: clusters: individual: Cl 0024+1654 - infrared: galaxies
Clusters of galaxies contain thousands of members within a region a few Mpc in diameter, and are the largest
known gravitationally bound systems of galaxies, having masses up to 1015
for the richest systems.
In hierarchical models clusters of galaxies grow by accreting less massive groups falling along filaments at a
rate governed by the initial density fluctuation spectrum, the cosmological parameters and the nature and amount
of dark matter. In the cluster environment, newly added galaxies are transformed from blue, active star forming
systems, to red, passive ellipticals, undergoing a morphological evolution stronger than that of field galaxies
at a similar redshift (Gavazzi & Jaffe 1987; Byrd & Valtonen 1990;
Abraham et al. 1996b). The cluster galaxy population is also characterized by a lower star
formation rate (SFR) than field galaxies of similar physical size and redshift (Couch et al.
2001; Lewis et al. 2002).
Butcher & Oemler (1978) showed that clusters of galaxies generally have a fraction fB of blue galaxies that increases with cluster redshift,
ranging from a value near 0 at z = 0, to 20% at z = 0.4 and to 80% at z = 0.9, suggesting a
strong evolution in clusters (Rakos & Schombert 1995). The galaxies responsible for the
Butcher-Oemler effect (hereafter BO effect) are generally luminous, spirals, and emission-line systems, with
disturbed morphologies. The high resolution imaging achieved by the Hubble Space Telescope (HST) has greatly
improved morphological studies of the cluster galaxy population over a wide range in redshift. The mixture of
Hubble types in distant clusters is significantly different from that seen in nearby systems. The population of
star-forming and post-starburst galaxies are disk dominated systems, some of which are involved in interactions
and mergers. (Abraham et al 1996a; Stanford et al. 1998; Couch et al.
1998; Van Dokkum et al. 1998; Morris et al.
1998; Poggianti et al. 1999; Best 2000).
Many mechanisms have been proposed to explain the complicated processes that occur in clusters, including ram
pressure stripping of gas (Gunn & Gott 1972), galaxy harassment (Moore et al.
1996; Moss & Whittle 1997), galaxy infall (Ellingson et al.
2001), cluster tidal forces (Byrd & Valtonen 1990; Fujita
1998) and interactions with other cluster galaxies (Icke 1985; Moss
& Whittle 1997). The main processes responsible for the morphological and spectral
evolution of cluster galaxies have yet to be determined. Ram pressure and tidal effects can quench the star
formation activity gradually because they operate over a period longer than 1 Gyr (e.g. Ghigna et al.
1998; Ramirez & de Souza 1998). Galaxy-galaxy interactions, galaxy
harassment and cluster mergers can enhance it and produce changes in galaxy properties over timescales of
100 Myr (e.g. Lavery & Henry 1986; Moore et al. 1996). Recent
changes in the properties of the galaxies may be detectable in the mid-infrared if associated with a burst of
star formation.
In the context of our ongoing exploitation of mid-infrared cluster data obtained with ESA's Infrared Space
Observatory (ISO, Kessler et al. 1996), we have analysed ISO observations of Abell 370,
Abell 2218 and Abell 2390 (Metcalfe et al. 2003; Altieri et al. 1999;
Biviano et al. 2004), Abell 2219 (Coia et al. 2005), and Cl 0024+1654 (this paper). In
common with other surveys we have found that mid-infrared observations of field galaxies from deep surveys
reveal a population of starburst galaxies that evolve significantly with redshift (Aussel et al.
1999; Oliver et al. 2000; Serjeant et al.
2000; Lari et al. 2001; Gruppioni et al. 2002;
Elbaz & Cesarsky 2003; Metcalfe et al. 2003; Sato et al.
2003). This class of sources are Luminous and Ultraluminous infrared galaxies (LIRGs and
ULIRGs, Sanders & Mirabel 1996; Genzel & Cesarsky 2000), have SFRs
of 100
(Oliver et al. 2000; Mann et al. 2002) and seem
to be almost always the result of galaxy-galaxy interactions at least in the local Universe (Veilleux et al.
2002). The ISOCAM sources account for most of the contribution of the mid-infrared to the
Cosmic Infrared Background (CIRB, Altieri et al. 1999; Franceschini et al.
2001; Metcalfe et al. 2001,2003; Elbaz et al.
2002). Studies of the global SFR show a decline by a factor of 3-10 since the peak of
star formation at z= 1-2 (Madau et al. 1996; Steidel et al. 1999;
Elbaz & Cesarsky 2003). The downturn in the global SFR and population of LIRGs and ULIRGs
may be caused by galaxies running out of gas available for star formation and the buildup of large scale
structure in the Universe that changed the environment of galaxies.
A study of the impact of the environment on galaxies in clusters could help in understanding the global SFR.
![]() |
Figure 1:
The 15 ![]() ![]() ![]() ![]() |
Open with DEXTER |
In this work we focus on the mid-infrared properties of the galaxy cluster Cl 0024+1654. The paper is organized as follows: Sect. 2 contains a description of the cluster. Section 3 describes the infrared observations and outlines the data reduction, source extraction and calibration processes. Section 4 presents the results, the model spectral energy distributions (SEDs) and star formation rates for cluster galaxies. Section 5 describes the spatial, redshift and colour properties of cluster galaxies and contains a description of Hubble Space Telescope images of some galaxies detected by ISOCAM. Section 6 makes a comparison between Cl 0024+1654 and other clusters studied, including Abell 1689, Abell 370, Abell 2390 and Abell 2218. The conclusions are in Sect. 7. The Appendix contains additional comments on some of the ISOCAM sources.
We adopt H0=70 km s-1 Mpc-1,
and
.
With this cosmology, the
luminosity
distance to the cluster is
DL = 2140 Mpc and 1
corresponds to 5.3 kpc at the
cluster redshift. The age of the Universe at the cluster redshift of 0.39 is 9.3 Gyr.
Cl 0024+1654 is a rich cluster of galaxies at redshift
(Humason & Sandage 1957; Gunn
& Oke 1975; Smail et al. 1993). It has a spectacular system of
gravitationally lensed arcs (Fig. 1) that were first observed by Koo (1988)
and subsequently studied by Mellier et al. (1991), Kassiola et al.
(1994), Wallington et al. (1995) and Smail et al.
(1997). The main arc is split into segments (Colley et al. 1996;
Tyson et al. 1998) and is the lensed image of a background blue galaxy at redshift z =
1.675 (Broadhurst et al. 2000). Unlike other clusters of galaxies, such as Abell 2218 or
Abell 2390 that have arc-like features, Cl 0024+1654 does not have a central dominant cD galaxy.
![]() |
Figure 2:
The fractional distribution of velocities of
cluster galaxies in Cl 0024+1654. The counterparts of the cluster members detected at 15 ![]() |
Open with DEXTER |
Cl 0024+1654 was one of the two clusters examined by Butcher & Oemler (1978) in their first published work on galaxy colours. The fraction of blue galaxies in Cl 0024+1654 is fB = 0.16which is much larger than the values of fB = 0.04 and fB = 0.03 for the Virgo (z = 0.003) and Coma (z = 0.02) clusters respectively (Butcher & Oemler 1984; Dressler et al. 1985; Schneider et al. 1986).
The determination of the spectral types of galaxies has been made for the cluster and the field (Czoske et al.
2002; Balogh et al. 2002). Using a wide field HST survey of Cl 0024+1654, Treu et al. (2003) found that the fraction of early-type galaxies (E+S0) is
highest ()
in the cluster core, declines rapidly to about 50% at
1 Mpc, and reaches the
background value of
at the periphery of the cluster (at a radius of
5 Mpc).
The mass profile of the cluster has been inferred from gravitational lensing analyses, kinematical analyses of
the redshifts of the cluster galaxies, and X-ray observations. The low value of the X-ray luminosity implies a
mass 2 to 3 times smaller than that predicted by lensing models (Soucail et al. 2000; Ota
et al. 2004). There are also differences among the lensing models. Broadhurst et al.
(2000) favour a cuspy NFW (Navarro et al. 1997) profile for
the mass distribution of the cluster. Tyson et al. (1998) argue that this result requires
an average cluster velocity distribution much higher than the measured value of
km s-1 (Dressler et al. 1999) and favour a uniform-density core for the halo
mass profile. The weak lensing analysis made by Kneib et al. (2003), based on an extensive HST survey of Cl 0024+1654 (Treu et al. 2003), indicates that the region within 5 Mpc from
the centre of the cluster is well fit by a steep NFW-like profile, and the isothermal fit is strongly rejected.
Ota et al. (2004) analysed the X-ray emission from the cluster and found that an
isothermal model with a mean temperature of 4.4 keV is a good fit to the temperature distribution. They also
estimated the mass of the cluster on the assumption that the intracluster medium is in hydrostatic equilibrium,
and confirmed that the mass inferred from X-ray observations is smaller than the mass predicted by lensing
models by a factor of about 3. The complex structure of the central region of the cluster might explain the
discrepancy (Zhang et al. 2005).
During an extensive spectroscopic survey of Cl 0024+1654, Czoske et al. (2001) detected a significant number of galaxies at a redshift slightly lower than that of the cluster (Fig. 2). This concentration (component B) is interpreted as a small, less massive cluster superimposed on the main cluster (component A) and lying at a mean redshift of z = 0.381. Numerical simulations indicate that the two clusters in Cl 0024+1654 are involved in a high-speed collision along the line-of-sight to the cluster (Czoske et al. 2002). The relative velocity of the two components, derived from the difference in redshift of the two main galaxy concentrations, is about 3000 km s-1. The collision scenario is not unusual because 30-40% of galaxy clusters have substructures that are detected in optical, X-ray and radio observations (e.g. Dressler & Shectman 1988; Schuecker et al. 2001; Mercurio et al. 2003). These substructures indicate that the hosting clusters are not fully dynamically relaxed and might be currently undergoing, or have recently undergone, mergers. The important results obtained by Czoske et al. (2001,2002) may give a new insight into the BO effect in medium and high redshift clusters. The impact of the two clusters could have triggered star formation in cluster galaxies, especially those on the leading edge of the smaller cluster, thus generating the large number of blue galaxies observed in Cl 0024+1654. The differences between the results from the methods used to determine the mass of Cl 0024+1654 are naturally explained by the cluster collision (Czoske et al. 2002).
Table 1:
Observational parameters used. The observations were made with the ISOCAM LW2
(5-8.5 m) and
LW3 (12-18
m) filters at their respective reference wavelengths of 6.7
m and 14.3
m. On-chip
integration time was always 5.04 seconds and the 3
per-pixel-field-of-view was used. M
and N are the
number of steps along each dimension of the raster, while dm and dn are the increments for a raster step. The
sensitivity in
Jy is given at the 5
level. The table also includes the total area covered and the
number of readouts per raster step. The rasters were repeated k times. Tot. t is the total time dedicated to
each filter observation.
The core of the cluster Cl 0024+1654 was observed at 7 m and 15
m using the LW2 and LW3 filters of
the camera ISOCAM (Cesarsky et al. 1996) on board ISO. The ISOCAM long wavelength (LW)
detector consisted of a
SiGa array. The observations were made on June 9 and June 15, 1997, in
raster mode, with an on-chip integration time of 5.04 s in the 3
per pixel field of view, and cover
an area of approximately 38 square arcminutes. A discussion of the relevant observing strategy for ISOCAM can be
found in Metcalfe et al. (2003) (with the exception that the Cl 0024+1654 rasters were not
"micro-scanned'' with sub-pixel finesse. I.e. the raster step sizes were multiples of the array pixel size,
whereas Metcalfe et al. describe rasters with step sizes involving fractions of a pixel's dimensions). For each
pointing of the raster, 14 readouts were performed at 7
m and 10 at 15
m. Fifty readouts were taken
at the beginning of each raster to allow for detector stabilization. The parameters used for the observations
are given in Table 1, which also lists the 5
sensitivity reached after data processing.
The diameter of the point spread function (PSF) central maximum at the first Airy minimum is
m) arcseconds. The FWHM is about half that amount and Okumura (1998) obtained
values of 3.3
at 7
m and 5
at 15
m for the PSF FWHM in the 3
per pixel
field-of-view. The data were reduced and, to take advantage of a slight deliberate offset between the two
15
m maps, were rebinned so that the final map has a pixel size of 1
,
with potential for slightly
improved spatial resolution. The maximum depth is reached toward the centre of the rasters, where the dwell time
per position on the sky is greatest. As recorded in Table 1, two rasters were made at 15
m
and a single raster at 7
m that had a larger step size and less observation time (by a factor of 7).
The data were reduced using the ISOCAM Interactive Analysis System or CIA (Delaney &
Ott 2002; Ott et al. 1997) in conjunction with dedicated routines, following the
method described in Metcalfe et al. (2003). The two 15
m raster maps were
normalized to their respective redundancy maps
and merged
into a single raster, thus increasing the sensitivity of the map to faint sources.
The 15 m sources were found by visually inspecting the merged and the two individual maps. Detections
common to the three maps are given in the primary source list (Table 2). The selection
criterion adopted is very conservative because it requires the source to be detected in both individual rasters,
resulting in the rejection of faint sources that can be detected only in the merged map. Therefore we include a
supplementary list containing sources found in the merged map and in at least one of the individual
rasters (Table 3). The fluxes of the sources were computed by aperture photometry using a circular
aperture of 9
diameter, centred on the infrared source.
Since there was only one raster at 7 m, an independent comparison of source detections was not possible
and real sources were considered to be only those having a 15
m counterpart.
The optical data were obtained on the ESO/VLT during October 2001 with the VLT/FORS instrument in service
mode. The data have been processed at the TERAPIX data
center
. Pre-calibrations, astrometric and photometric calibrations as well
as image stacking and catalog production were done using standard CCD image processing algorithms and tools
available at the TERAPIX center (see McCracken et al. 2003 for details).
Table 2:
The primary list of sources detected at 15 m. The columns are (from left to right):
source identification number in order of increasing Right Ascension; source flux and precision before
calibration (in ADUgs); source flux and precision after calibration (in mJy); Right Ascension and Declination of
the infrared source (J2000); redshift and name of the optical counterpart, if known. The redshift of source
ISO_Cl0024_02 was taken from Schmidt et al. (1986). All other redshifts are from Czoske
et al. (2001) or provided by T. Treu (private communication). The names of the stars are
from HEASARC.
Table 3:
The supplementary list of 15 m sources. The meaning of the columns is the same as in
Table 2.
Table 4:
The list of 7 m sources with source flux exceeding 5
of the local noise, and having
15
m counterparts. From left to right: 15
m source identification number as listed in
Table 2; signal and precision of flux in ADU; flux and precision in mJy; Right Ascension, Declination;
nature of the source; and finally the [7
m]/[15
m] colour ratio.
Monte Carlo simulations were performed to calibrate the complex data reduction process by characterizing how it
affects model point sources with known properties inserted into the raw data, and to characterize the precision
with which source signals could be recovered from the data. The procedure adopted is fully described in Metcalfe
et al. (2003), the only exception being that all photometry performed for the present work
employed the XPHOT routine in CIA, whereas Metcalfe et al. used SExtractor (Bertin &
Arnouts 1996). The fluxes of the inserted fake sources were measured in exactly the same
way as the fluxes of the real sources. The simulations were performed independently for all rasters. The fluxes
of the inserted fake sources were measured by manual photometry adopting the same parameters, i.e. same
aperture, used for the real sources. The simulations establish the relationship between source signal detected
in the photometric aperture, and total source signal actually collected in the detector. Signal in the
photometric aperture can therefore be scaled to total signal, but at this point the units are detector units, or
ADUgs. The relationship between source signal deposited in the detector and actual source flux-density in mJy is
determined by two further scaling factors. (a) A filter specific ISOCAM calibration factor relating stabilised
source signal, in detector units, to mJy, i.e. 1 ADU per gain per second
mJy for LW2 and
1 ADUgs
0.51 mJy for LW3 (Delaney & Ott 2002); and (b) a correction for detector responsive
transient effects which cause the signal from faint sources to fall below that expected on the basis of bright
reference source measurements. The derivation of that scaling factor, in the context of the reduction algorithm
applied here, is described in Metcalfe et al. (2003) and references therein.
The merged 15 m map (Fig. 1) is overlaid on a V-band image of the cluster taken with the Very Large
Telescope (VLT).
The list of sources detected in the two individual 15 m maps is given in Table 2 and a
supplementary list is given in Table 3 (see Sect. 3). The name of the ISOCAM source is
composed of the satellite acronym (ISO), the partial name of the cluster (Cl0024), and an identification number
assigned to each source.
The 7 m fluxes for the stars, and their [7
m]/[15
m] flux ratios, are given in
Table 4. The labels of the 7
m sources are taken to be the same as their 15
m counterparts.
The fourth star detected at 15
m (ISO_Cl0024_01) is outside the boundary of the 7
m raster. No
other sources were detected above the 5
limit of
400
Jy.
The column named "Precision'' in Tables 2 and 3 reflects the repeatability of the
photometric results in the fake source simulations. For each source brightness the precision is
the 1scatter found in the recovered fluxes of fake sources of similar brightness inserted into the raw data. It
should be noted that the ratio signal/precision is not, however, the signal-to-noise ratio; nor is it the
significance of a source detection. The precision, as defined here, includes, for example, signal deviations
caused by residual glitches around the map and by variations in residual background level across the map, which
are of course factors that can influence the accuracy of source photometric measurements. However, the
significance of an individual source detection must be determined through an examination of the pixel-to-pixel
signal variations in the map local to the source. All of the sources listed in Tables 2
and 3 have significance values greater than 5
.
The quoted precision is a measure of the relative
calibration accuracy within the sample recorded. The absolute flux calibration is, in addition, affected by
factors such as imperfections in the knowledge of the detector response and responsive transient correction, and
differences between the spectral shapes of measured sources with respect to the canonical calibration source for
ISOCAM, which has a stellar spectrum. Experience over many surveys suggests an absolute calibration accuracy of
15% for ISOCAM faint-source photometry due to accumulated systematic effects. The measured colour ratios
for the stars listed in Table 4 are consistent with this, and give some confidence in the absolute
calibration at the
15% level.
The number of sources in the primary and supplementary lists are 30 and 5, respectively, yielding a total of 35 sources. A search was performed for counterparts of the ISOCAM sources at other wavelengths, using a search
radius of 6
in HEASARC
and NASA/IPAC
Extragalactic Database
. Of the 30 sources in the primary list (Table 2), 19 are identified with catalogued
sources. Of these nineteen, four are stars, one is a quasar (PC 0023+1653), ten are cluster members, one is a
foreground galaxy and three are background galaxies. As for the 5 supplementary sources (Table 3),
three are identified with cluster galaxies. There is a total of 13 sources with unknown redshift and all have
optical counterparts (Fig. 1).
The median flux-density of the unknown-redshift sources is about 250 Jy. The expected number of 15
m sources lensed by the cluster down to that flux limit is about
,
based on the log N - log Sdistribution of Metcalfe et al. (2003). So the expected number of background sources is
sufficient to account for the observed number of unknown-redshift sources.
Table 5: References to ISOCAM sources in existing archives. The first column lists the source number from Tables 2 and 3. The key to the references in Cols. 2 to 8 is Dressler & Gunn (1992, DG), Czoske et al. (2001, CKS), Schneider et al. (1986, SDG), Smail et al. (1997, SDC), McLean & Teplitz (1996, MT), Butcher & Oemler (1978, BO), Soucail et al. (2000, S) and Pickles & van der Kruit (1991, P). The table does not include the four stars or the quasar PC 0023+1653.
![]() |
Figure 3:
SEDs for cluster members in Cl 0024+1654 detected at 15 ![]() |
Open with DEXTER |
Spectral energy distributions were computed for the cluster galaxies using the 15 m fluxes in
Tables 2 and 3, and a 5
upper limit of 400
Jy for the 7
m fluxes. The
spectral range of the ISOCAM sources was extended by including measurements in the optical photometric bands and
the near-infrared. The values were retrieved from the NASA/IPAC Extragalactic Database. The references to
extensive observations of the ISOCAM sources from a wide range of catalogues are listed in Table 5.
The SEDs were modelled using the program GRASIL (Silva et al. 1998). These models have already been used by Mann et al. (2002) to fit the SEDs of galaxies detected by ISOCAM in the Hubble Deep Field South. Given the limited photometric accuracy and, for some galaxies, the limited amount of available photometric data, it is not possible to obtain an accurate model of the SED of each individual galaxy. Instead, observed SEDs are compared with models representative of broad classes of spectral types. The number of models to be considered needs to be adequate to reflect the quality of the data and the questions to be answered (see, e.g., Rowan-Robinson 2003). The redshifts of the galaxies are known and hence it was possible to use more than the five models of Mann et al. (2002) without the risk of multiple solutions.
In this work, 20 models were considered. They were taken from either the public GRASIL library, or built by
running the publicly available GRASIL code, or kindly provided by L. Silva (private comm.). The
20 models
reproduce the SEDs of several kinds of galaxies, from early-type, passively evolving ellipticals (labelled "E'',
in the following), to spiral galaxies (labelled "S''), and starburst galaxies (similar to the local examples
Arp 220 and M 82), as seen either at the epoch of the starburst event, or 1 Gyr after (labelled "SB'' or "SB+1'',
respectively). At variance with Mann et al. (2002) we avoid considering models older than
the age of the Universe at the cluster redshift. Specifically, we consider models for galaxies with a formation
redshift z=1 (this is the redshift when stars start forming in the GRASIL models), corresponding to an age at
the cluster redshift of 3 Gyr, and models with a formation redshift z=4, corresponding to an age at the
cluster redshift of
8 Gyr.
The best-fitting SED model was determined by
minimization, leaving the normalization of the model free.
In principle, a given SED model cannot be freely rescaled, since the parameters in the GRASIL code depend on the
mass of the galaxy. However, a fine tuning of the values of the GRASIL parameters only makes sense if the
photometric uncertainties of the observed galaxy are small enough, and the SED is very well covered by
observations, which is not the case here. As a matter of fact, the model SEDs of two galaxies of very different
mass, such as M 82 and NGC 6090, could not be distinguished after rescaling, in the wavelength range covered by
our observations, and with the accuracy of our photometric data points.
The SEDs are plotted in Fig. 3 and summarized in Table 6 for all cluster galaxies detected with ISOCAM, except for ISO_Cl0024_12, for which there are not enough photometric data. All the fits are acceptable at the 1% confidence level, and all, except two, at the 5% confidence level. S and SB+1 models provide the best-fit to most SEDs, but the fits are, in general, not unique. Among the acceptable fits (at the 5% confidence level) there are also SB and E models, but the latter always underestimate the mid-infrared fluxes. We are unable to distinguish between models with different formation redshifts.
Table 6:
Results from model SEDs for known cluster galaxies in Cl 0024+1654. The columns list: LW3
identification number; best-fit SED model; number of data points in the observed SED;
of the fit;
rejection probability of the fit; other acceptable models (rejection probability 95% or less); 15
m
k-correction; luminosity at 15
m (h=0.7,
,
)
using the transformation from
Boselli et al. (1998), and k-corrected using best-fit SED; total infrared luminosity
including the k-correction using the best-fit SED and computed using the transformation from Elbaz et al.
(2002); SFR from total IR luminosity using the transformation of Kennicutt
(1998); position in the V-I vs. I colour-magnitude diagram (Fig. 7) where
"MS'' stands for "Main Sequence'', "BO'' for "Butcher-Oemler'', "U'' for "Unknown''; EW [O II] line from
Czoske et al. (2001); SFR computed from the [O II] emission line using the
transformation to the optical SFR; ratio infrared to optical SFR.
Overall, the results show that most ISOCAM cluster sources have SEDs typical of actively star-forming galaxies, although not necessarily observed in an exceptional starbursting phase.
The SEDs of the nine cluster sources from the primary list were combined assuming the mean redshift of the
ISOCAM cluster sources, and normalising the individual SEDs with the flux density in the rest-frame H band,
either observed or as predicted by the best-fit SED models (see also Biviano et al. 2004). We did
not consider the three supplementary sources; since their 15 m flux values are rather uncertain, including
these sources would increase the scatter of the average 15
m flux density value. The resulting average
observed SED is shown in Fig. 4, along with a SB+1 model fit to guide the eye
(S models fit the
average SED equally well). The average SED provides a synthetic view of the spectral shape of the 15
m
cluster galaxies, but we refrain from a physical interpretation of it, which would require a detailed analysis
of the systematics inherent to its construction.
![]() |
Figure 4:
The average SED of 9 ISOCAM cluster
sources (the sources from the supplementary list have not been considered). Error bars give the rms values for
the flux density. An upper limit of 400 ![]() ![]() |
Open with DEXTER |
The mid-infrared emission, free from dust extinction, is a reliable tracer of star formation (Genzel & Cesarsky 2000). The infrared emission from a galaxy is the sum of various contributions including:
![]() |
Figure 5:
The total infrared luminosity distribution for
cluster members. The median value for
![]() ![]() ![]() ![]() |
Open with DEXTER |
The luminosities at 15 m, in units of solar luminosity, were obtained using the relationship (Boselli et al. 1998):
Half of the 15 m cluster galaxies are LIRGs and four more are within 1
of the luminosity of a
LIRG. Star forming episodes that are enshrouded by dust are commonly found in these galaxies in the field. The
cause of the starburst phase of LIRGs in the local universe is still debated, but in the range of infrared
luminosities measured with ISOCAM the starburst seems to be linked to mergers or interactions between pairs of
spiral galaxies of unequal size and single galaxies in about one-quarter of the sources (Hwang et al.
1999; Ishida & Sanders 2001). The cause of the starburst in the
15
m sources of Cl 0024+1654 seems to be consistent with what is observed in the local universe because of
the evidence for interactions and mergers in the HST maps (see Sect. 5.2).
![]() |
Figure 6:
a) Adaptive Kernel map for cluster members. Black boxes represent the 15 ![]() ![]() |
Open with DEXTER |
The Equivalent Width (EW) of the [O II] line was available from Czoske et al. (2001) for
only three of the mid-infrared cluster members. The [O II] line is redshifted into the V band and can be used to
compute the optical SFR provided that the V-band magnitude is known. The luminosity
was derived
from the V-band magnitude in units of erg s-1
using:
![]() |
(4) |
In our sample, AGN activity does not seem to be a major contributor to the energy output of the systems.
Measurements of the H
and [N II] emission lines are required to further determine the AGN contribution to
the mid-infrared emission (Balogh et al. 2002) and future radio observations could also
help in determining the AGN contribution (Rengarajan et al. 1997; Dwarakanath & Owen
1999).
The spatial and radial distributions of the cluster galaxies, with and without detectable 15 m emission,
are given in Fig. 6. The isodensity contours and dots in Fig. 6a refer to
galaxies belonging to component A, while the black boxes represent the cluster members with 15
m emission.
The figure shows that the cluster members detected at 15
m are less concentrated than the cluster galaxies
not detected in the mid-infrared. The Kolmogorov-Smirnov test reveals that there is only a 4% probability that
the two distributions are drawn from the same parent one.
The velocity distributions for cluster galaxies with and without detectable 15 m emission are given in
Fig. 2. The counterparts of the 15
m sources have redshifts which place them in Component
A - the larger of the two interacting components. One source, ISO_Cl0024_12, with z = 0.386 is close to the
boundary. The velocity dispersion computed in the cluster rest-frame (Fig 2), using the
biweight estimator (Beers et al. 1990) in the central region of the ISO map, is
km s-1 for a sample of galaxies, while it is
km s-1 for the 15
m galaxies. The 15
m sources are
associated mainly with spiral and emission-line galaxies, and studies of nearby clusters indicate a larger
velocity distribution for this population relative to the general cluster population
(Biviano et al. 1997).
![]() |
Figure 7:
V - I colour-magnitude diagramme for component-A galaxies (empty circles) within the region mapped by
ISOCAM. The eleven 15 ![]() ![]() |
Open with DEXTER |
![]() |
Figure 8:
Miniature maps obtained from HST. The maps are
centred on the 15 ![]() ![]() |
Open with DEXTER |
The colour-magnitude diagramme for cluster galaxies is given in Fig. 7. It provides an important
link between the optical and mid-infrared properties of the ISOCAM galaxies. Apparent magnitudes in the V and I bands are available from Czoske et al. (2001) for 11 of the
13 ISOCAM cluster galaxies.
The circles represent component-A galaxies within the region mapped by ISOCAM. The spectroscopically confirmed
mid-infrared cluster galaxies are indicated by filled circles and are labelled with the ISO source number.
Section M of Fig. 7 contains the galaxies that satisfy the BO definition, given in Sect. 1. Six of the 15 m galaxies fall in section M and satisfy the BO requirement. Only 20% of the
cluster galaxies in section M are detected at 15
m. However, this fraction increases to
38% when
only galaxies brighter than
mI= 19.8 are considered. The non-detections of fainter I-band sources
are probably due to the sources falling below the sensitivity limits of ISOCAM.
All 15 m cluster sources have in common an excess of mid-infrared emission that sets them apart
from the rest of the cluster members. In at least some of these systems, one would expect that the star forming
activity is primarily triggered by interactions and mergers. Recent changes in the galaxy population reveal
themselves in the mid-infrared and help identify the processes that cause the burst of obscured star formation.
The infall of smaller groups into the cluster environment provides a way of promoting slow encounters and
mergers within clusters (Mihos 2004). Slow encounters between galaxies are better able to
drive instabilities than fast encounters. Complex interactions of the kind that may occur between galaxies in
this cluster have been studied in the local Universe in Hickson Compact Groups (see, e.g., López-Sánchez et al. 2004, and references therein).
Table 7:
Summary of ISOCAM observations and results at 15 m for five clusters of galaxies. The data was
obtained from Metcalfe et al. (2003) for Abell 370, Abell 2218 and Abell 2390, Fadda et al. (2000) and Duc et al. (2002) for Abell 1689, and this paper for
Cl 0024+1654. The content of the columns in the table are as follows: name and redshift of the cluster, total
area scanned, sensitivity reported at the
level, flux of the weakest reported source
in
Jy, total
observation time, total number of sources detected including sources without redshifts and stars. Then number of
cluster galaxies, virial radius, virial mass, number of sources with
detected and expected. The expected number of sources was obtained by comparison with Cl 0024+1654 as described
in the text. Virial radii and masses are from Girardi & Mezzetti (2001) and King et al.
(2002).
We looked for indications of ongoing and past interaction events among the 15 m cluster sources, either in
their morphologies, or in their spectra. A HST image was available for part of Cl 0024+1654 from the MORPHS
collaboration (Smail et al. 1997). This image is limited to the South-Eastern region of
the cluster core and contains seven of the 15
m cluster sources. For these galaxies, miniature maps
centred on the ISOCAM source coordinates and having an area of
arcsec2 were obtained. The
miniature maps are presented in Fig. 8.
Miniature maps are available for 3 of the BO galaxies in section M of Fig. 7. The HST maps show
that ISO_Cl0024_18 and ISO_Cl0024_22 do not have companions within a radius of 5
while
ISO_Cl0024_08 appears to be interacting with at least another galaxy.
HST miniature maps for four of the non-BO 15 m emitters in section N are also given in Fig. 8 and
each has at least one nearby companion. ISO_Cl0024_4s appears to be merging because it has multiple nuclei and
has at least two other companions. ISO_Cl0024_24 and ISO_Cl0024_3s show evidence of ongoing interactions
(e.g. tidal bridges) with at least one nearby galaxy.
Even though the number of infrared sources with HST maps is small, there seems to be a difference between the BO
and non-BO galaxies in the sense that the former tend not to have nearby companions whereas the latter do. The
optical counterparts of the 15 m sources on the main colour-magnitude sequence are not passive early type
galaxies as one would expect from their position in the colour-magnitude diagramme. Rather, the optical
counterparts appear to be galaxies that are involved in interactions and mergers in the HST images.
There is indirect support for this apparent difference between the BO and non-BO galaxies from the results of
the SED fitting process. Five of the six 15 m galaxies in the BO region of the diagram have SEDs that are
best-fit by S-type models. On the other hand, two out of four 15
m galaxies on the main sequence of the
colour-magnitude diagram have SEDs that are best-fit by SB+1-type. If the starburst events are triggered by
interactions, there is marginal evidence from the SED analysis that the 15
m galaxies on the main
colour-magnitude sequence have suffered more interaction events in the recent past than the BO 15
m
galaxies.
Given the scarce morphological data, and the fact that the SED fitting results are not very well
constrained (SB+1-type models also provide acceptable fits to the SEDs of BO galaxies, see Table 6),
we consider the evidence for a difference in the interaction properties of BO and non-BO 15 m galaxies only
tentative. We conclude from our data that interaction events do trigger the star formation activity (and
stimulate the IR emission) in at least some 15
m cluster galaxies.
Obscuration by dust could explain why some 15 m sources lie on the cluster main sequence. We tested this
possibility by computing the de-reddened V and I magnitudes of the 15
m galaxies, using the dust-free
best-fit model SEDs, as given by the GRASIL code. As expected, once de-reddened the 15
m galaxies on the
colour-magnitude sequence move in the colour-magnitude diagramme toward the BO region, and in one case
(ISO_Cl0024_3s), it even becomes part of the BO category (Fig. 7). Dust obscuration at least
partially explains the unusual colours of the 15
m sources on the cluster main sequence.
Observations of clusters of galaxies at high redshift (z>0.4) with SCUBA have yielded an excess of
sub-millimeter sources with high values of the SFR (Best 2002). Here we compare the
results of the observations of Cl 0024+1654 with four clusters of galaxies that were observed with ISOCAM at
15 m (Fadda et al. 2000; Metcalfe et al. 2003) to determine if
Cl0024+1654 has more LIRGs than the other clusters. We have not included the cluster J 1888.16 Cl in Table 7 because the relevant data have not been published (Duc et al. 2004). However this
cluster seems to be comparable to Cl0024+1654 because it has many infrared sources and at least 6 of them are
confirmed LIRGs. The characteristics of the observations, and results for the five clusters, are summarized in
Table 7. The observed areas and sensitivities are comparable for the first three clusters in
Table 7, but the area is somewhat smaller for Abell 2218, and much smaller for the ultra-deep
observations of Abell 2390.
The number of 15 m sources that are identified with cluster galaxies ranges from 13 in Cl 0024+1654 to 1 in
Abell 370. The number of 15
m sources that have fluxes consistent with LIRGs, within the precision of the
measurements is listed in Table 7. In Cl 0024+1654 the number of LIRGs is 10 (Table 7) and
includes 6 sources with
and 4 more with
between
and
(Table 6). The number of LIRGs detected in the other clusters is
either 0 or 1 (Table 7). The virial mass and radius of Cl 0024+1654, Abell 370 and Abell 1689 are the
same to within 20% and are considerably smaller than the values for the more massive clusters Abell 2218 and
Abell 2390.
However the ratio of the virial mass to area within the virial radius varies by less than 50% for the 5
clusters. We now compare the number of LIRGs in Cl 0024+1654 with those in the other clusters. The comparison
is limited to LIRGs because it eliminates selection effects for the 15 m sources. In any case most of the
15
m sources in the nearer clusters are too faint to be detected in the more distant clusters Cl 0024+1654
and Abell 370. In the comparison, the number of LIRGs in Cl 0024+1654 was multiplied by the ratios of:
Large numbers of cluster sources are detected by ISOCAM at 7 m in both Abell 1689 and Abell 2218. In Abell
2218 the SEDs of most of these 7
m sources are well fit by models of quiescent ellipticals with negligible
SFRs and a median luminosity of
(Biviano et al. 2004).
Cl 0024+1654 is involved in a recent merger, as shown by the two velocity components in Fig. 2.
Abell 1689 also has a complex and broad velocity structure, with evidence for three distinct groups that overlap
spatially and are well separated in velocity space (Girardi et al. 1997). The two clusters
Abell 370 and Abell 2218 show no evidence for major merger activity. The number of luminous 15 m sources
can be, in part, explained by considering the dynamical status of the cluster. Simulations made by Bekki
(1999) show that the time-dependent tidal gravitational field existing in cluster-group
mergers induces secondary starbursts by efficiently transferring large amounts of gas from the disk to the
nucleus. The model establishes a link between the population of starburst or post-starburst galaxies and the
presence of substructures in clusters. The existence of starburst and post-starburst galaxies in a merging
cluster, and spread over a wide area, is an important prediction of this model which seems to be in agreement
with the results from Cl 0024+1654 and Abell 1689.
Furthermore it can be expected that merging clusters at much larger redshifts than Cl 0024+1654 will contain LIRGs and ULIRGs because the interacting and merging galaxies will be more gas rich (see, e.g. Balogh et al. 2000; Kauffmann & Haehnelt 2000; Francis et al. 2001).
In this context it is interesting to note that Dwarakanath & Owen (1999) found different radio source populations in two very similar clusters. The number density of radio sources in Abell 2125 (z = 0.246) exceeded that in Abell 2645 (z = 0.25) by almost an order of magnitude. The cluster Abell 2125, with the larger number of radio sources, also shows evidence for a merger.
Many infrared sources were found in the cluster Abell 1689 (Fadda et al. 2000) that is at
a much smaller redshift than Cl 0024+1654. No galaxy in Abell 1689 has total infrared luminosity above
(we use our cosmology and Eq. (2) to convert the 15
m fluxes from
Fadda et al. (2000) to total infrared luminosities), while half of the 15
m cluster
sources in Cl 0024+1654 are above this luminosity. All of the 15
m sources detected in Cl 0024+1654 would
have been easily detected in Abell 1689, if present, whereas most of the sources in Abell 1689 would not have
been detected in the observations of Cl 0024+1654 when allowance is made for different sensitivities and
distances. The two clusters have very similar virial radii and masses (Table 7). However a much larger
part of the periphery of Cl 0024+1654 was observed because it is at a higher redshift. It is possible that some
of the differences between the 15
m sources arises from a change in the source population between the core
and periphery of the two clusters.
The population of 15 m sources may also be influenced by the history of the cluster. Various processes
such as ram pressure stripping of gas from cluster galaxies, tidal effects and previous starbursts will
inevitably leave less gas in the interacting and merging galaxies to fuel the starburst (e.g. Gunn & Gott
1972; Byrd & Valtonen 1990; Fujita 1998).
These processes may vary from one cluster to another depending on its history and provide dispersion in the
luminosity of the starbursts. The large number of LIRGs in Cl 0024+1654 are fuelled by the gas rich progenitor
galaxies. It is conceivable that in Abell 1689 the starbursts are less luminous either because there is a
smaller supply of gas to fuel the outburst, or the luminosity has decreased because a longer time has elapsed
since the last merger or there is a significant difference between the sources in the core and periphery.
The average value of the ratio
is reasonably comparable
for both clusters (it is
15 for Cl 0024+1654 and
11 for Abell 1689). The internal properties of
the 15
m sources in Cl 0024+1654 and Abell 1689 are broadly similar even though their distributions of the
luminosities of 15
m sources are different. The colour-magnitude diagrammes of the 15
m sources in
the two clusters are also quite similar.
The cluster Cl 0024+1654 was observed with ISO. A total of 35 sources were detected at 15 m and all have
optical counterparts. Sources with known redshift include four stars, one quasar, three background galaxies, one
foreground galaxy and thirteen cluster galaxies. The remaining 13 sources are likely to be background sources
lensed by the cluster.
The spatial, radial and velocity distributions were obtained for the cluster galaxies. The ISOCAM cluster
galaxies appear to be less centrally grouped (in the cluster) than those not detected at 15 m and the
Kolmogorov-Smirnov test reveals that there is only a 4% probability that the two distributions are drawn from
the same parent population. No statistically significant differences were found between the velocity
distributions of the 15
m sources and other cluster galaxies in the region mapped by ISOCAM.
Spectral energy distributions were obtained for cluster members and used as indicators of both morphological type and star forming activity. The ISOCAM sources have as best-fit SEDs predominantly those of spiral or starburst models observed 1 Gyr after the main starburst. Star formation rates were computed from the infrared and the optical data. The SFRs inferred from the infrared are one to two orders of magnitude higher than those computed from the [O II] line emission, suggesting that most of the star forming activity is hidden by dust.
A colour-magnitude diagramme is given for cluster sources falling within the region mapped by ISOCAM. V and
I-band magnitudes are available for 11 of the cluster sources, and 6 of these have colour properties that are
consistent with Butcher-Oemler galaxies and best-fit SEDs that are typical of spiral models. The remaining
15 m cluster galaxies have colours that are not compatible with Butcher-Oemler galaxies and have best-fit
SEDs that are typical of starburst galaxies 1 Gyr after the main burst. HST images are available for these
latter systems and all have nearby companion galaxies. These results suggest that interactions and mergers are
responsible for some of the luminous infrared sources in the cluster.
The 15 m sources in Cl 0024+1654 were compared with four other clusters observed with ISOCAM. The results
show that the number of LIRGs in Abell 370 is smaller than expected by about one order of magnitude, if Abell
370 were to be comparable with Cl 0024+1654. Furthermore no LIRGs were detected in Abell 1689, Abell 2218 and
Abell 2390 when a total of 3 was expected, based on the results from Cl 0024+1654. A comparison of the
mid-infrared sources in Abell 1689 and Abell 2218 shows that the sources in Abell 1689 are more luminous than in
Abell 2218 and follow the same trend identified in the comparison between Cl 0024+1654 and Abell 370. There is
clear evidence for an ongoing merger in Cl 0024+1654 and Abell 1689. The number and luminosities of the
mid-infrared cluster sources seem to be related to the dynamical status and history of the clusters.
Acknowledgements
D.C. and B.McB. thank Enterprise Ireland for support. We thank Laura Silva for valuable discussions regarding the spectral energy distributions. We also thank T. Treu and S. Moran for redshift informations. D.C. gratefully acknowledges the hospitality of ESA's ISO Data Centre (IDC) at Vilspa, Spain, where part of this work was completed. JPK acknowledges support from CNRS as well as from Caltech. The ISOCAM data presented in this paper was 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. We thank the anonymous referee for the many comments that improved the content of the paper.