Conference Report CSIRO PUBLISHING Publications of the Astronomical Society of Australia, 2005, 22, 326–334 www.publish.csiro.au/journals/pasa

Galaxy Groups: Proceedings from a Swinburne University Workshop

Virginia A. KilbornA,F, Kenji BekkiB, Sarah BroughA, Marianne T. DoyleC, Ekaterina A. EvstigneevaC, Duncan A. ForbesA, Bärbel S. KoribalskiD, Matthew S. OwersB, Chris PowerA, Michael J. DrinkwaterC, David J. RohdeC, Christopher A. BlakeE, Warrick J. CouchB, Michael B. PracyB, and Brad K. GibsonA

A Centre for Astrophysics and Supercomputing, Swinburne University, Melbourne VIC 3122, Australia B School of Physics, University of New South Wales, Sydney NSW 2052, Australia C Department of Physics, University of Queensland, Brisbane QLD 4072, Australia D Australia Telescope National Facility, CSIRO, Epping NSW 1710, Australia E Department of Physics and Astronomy, University of British Columbia, Vancouver V6T 1Z1, Canada F Corresponding author. Email: [email protected]

Received 2005 September 9, accepted 2005 October 24

Abstract: We present the proceedings from a two-day workshop held at Swinburne University on 2005 May 24–25. The workshop participants highlighted current Australian research on both theoretical and observa- tional aspects of groups. These proceedings include short one-page summaries of a number of the talks presented at the workshop. The talks presented ranged from reconciling N-body simulations with observa- tions, to the Hi content of in groups and the existence of ‘dark galaxies’. The formation and existence of ultra-compact dwarfs in groups, and a new supergroup in were also discussed.

Keywords: galaxies: clusters: general — methods: N-body simulations — cosmology: theory — dark matter — gravitation: catalogues — surveys — galaxies: photometry — radio lines: galaxies

The majority of galaxies in the Universe lie in groups in groups might be one of the most important factors in (Tully 1987; Eke et al. 2004a). However, the physical their evolution. processes operating in groups are poorly understood. Observations of neutral hydrogen (Hi) in galaxy groups For example, to what extent do gravitational interactions can provide information about the processes that have and the intra-group medium determine the morphology occurred in the groups. In particular, evidence of previous and formation properties of galaxies residing in and galaxy interactions is often present in the neutral hydrogen around groups? Simulations of groups of galaxies have not where none is obvious in the optical (e.g. M81;Yun,Ho, & been reconciled with the observations: In general, many Lo 1994). In the dense cluster environment, ram pressure more satellite galaxies are predicted than are seen in galaxy stripping of Hi from gas-rich spiral galaxies is observed groups (e.g. Klypin et al. 1999; Moore et al. 1999). (e.g. Vollmer et al. 2001; Oosterloo & van Gorkum 2005), Environment appears to play a key role in the evolu- corresponding with the observation that spiral galaxies tion of galaxies. The star formation rate of galaxies in and near the centre of clusters tend to be Hi-deficient (Solanes around (out to two to three virial radii) clusters is lower et al. 2001). There is little corresponding data on the Hi than that of the field (Lewis et al. 2002; Gomez et al. 2003). content of spiral galaxies in groups. The Hi content of The galaxy density at these radii are comparable to that of galaxies in compact groups tends to be deficient accord- groups. There are several possible mechanisms for the star ing to Verdes-Montenegro et al. (2001), whereas Stevens formation to be quenched such as removal of the gas via et al. (2004) find no evidence of Hi deficiencies in this ram pressure stripping (Abadi, Moore, & Bower 1999; environment. The few studies of loose groups tend to Vollmer et al. 2001; Kenney, van Gorkum, & Vollmer support moderate Hi deficiencies, caused by gravitational 2004), gravitational interactions (Toomre & Toomre 1972; interactions (Kilborn et al. 2005; Omar & Dwarakanath Vollmer 2003; Bekki et al. 2005a, 2005b), and galaxy 2005). mergers (Zabludoff & Mulchaey 1998). The latter two Gravitational interactions may also play a part in the are more likely to occur in galaxy groups where the rel- formation of a new class of galaxies, the ultra-compact ative velocities of galaxies are lower. The observed drop dwarfs (UCD). UCDs were first discovered in the For- in star formation rate at large cluster radii suggests that nax cluster by Hilker et al. (1999) and Drinkwater et al. ram pressure stripping cannot be solely responsible, as (2000), and while they are unresolved in ground-based this only occurs in the dense intra-cluster medium near optical images, their spectra show they are located at the the centre of clusters. Thus the pre-processing of galaxies cluster velocity. UCDs have luminosities between that of

© Astronomical Society of Australia 2005 10.1071/AS05030 1323-3580/05/04326 Galaxy Groups: Proceedings from a Swinburne University Workshop 327

globular clusters and dwarf galaxies (Bekki et al. 2003). 6 × 1010 M. The centre of the frame coincides with the Bekki et al. (2003) proposed that UCDs were formed centre of the larger galaxy, and the companion galaxy has from galaxy ‘threshing’of dwarf galaxies, where the outer gone away from the larger one so that it can not be seen in parts of the galaxy were tidally removed, leaving just the this frame. It is clear from this Figure that new can be nucleus. formed in the denser part of the intragroup Hi ring owing This workshop addressed some of these current issues to tidal compression of the Hi gas. The gaseous region in galaxy groups, and summaries of nine of the presenta- around these new stars can be identified as intragroup Hii tions follow. regions as observed in NGC 1533 by Ryan-Weber et al. (2003). These results imply that galaxy–galaxy interaction 1 Formation of Intergroup and Intragroup and the resultant tidal stripping can be responsible for the Stellar and Gaseous Objects formation of intragroup Hi gas clouds and intragroup Hii Kenji Bekki regions. We investigate the formation of apparently isolated Hi 2 A Supergroup in 6dFGS gas and Hii regions in group of galaxies based on gas Sarah Brough dynamical simulations of disk galaxies in groups. Details on the disk galaxy models and treesph code adopted in Hierarchical structure formation leads us to presume that the present study have already been described in Bekki et clusters of galaxies are built up from the accretion and al. (2002) and Bekki et al. (2005a, 2005b), so we give only merger of smaller structures like galaxy groups (e.g. Blu- a brief review here. We investigate the dynamical evolu- menthal et al. 1984). Although we observe clusters of tion of stellar and gaseous components (including globular galaxies forming along filaments and accreting galaxy clusters) in an interacting pair of late-type disk galaxies group-like structures (e.g. Kodama et al. 2001) we lack with the mass ratio smaller than 0.1 (with the larger disk clear examples of groups merging to form cluster: ‘super- similar to the Galaxy) in groups of galaxies. Star forma- groups’. tion is modelled as the Schmidt law (Schmidt 1959) with The Eridanus cloud lies at a distance of ∼21 Mpc and an exponent of 1.5 and a threshold gas density for star for- includes two optically classified groups of galaxies, NGC mation (Kennicutt 1998) is also included. It is found that 1407 and NGC 1332. These groups are part of the Group (1) massive isolated Hi clouds with the masses of 109 M Evolution Multiwavelength Study (Osmond & Ponman can be formed as a result of tidal stripping during galaxy– 2004) in which their X-ray properties were analyzed. The galaxy interaction in groups (Bekki et al. 2005a, 2005b), NGC 1407 group shows X-ray emission from intra-group (2) star formation can occur in some of these isolated Hi gas, indicating the presence of a massive structure. In con- gas clouds if the companion galaxies are more massive and trast, the X-ray emission from the NGC 1332 group is have a larger amount (5 × 109 M) of gas initially before associated with NGC 1332 itself, not with intra-group gas. tidal interaction (Bekki et al. 2005a), and (3) intragroup Omar & Dwarakanath (2005) suggest that there is intra- GCs can be formed by tidal stripping (Yahagi & Bekki group gas associated with NGC 1395, however, no group 2005). has previously been associated with this galaxy. Figure 1 shows the final distributions of gas and new At present there is some debate as to the nature of the stars after tidal stripping of Hi gas from the companion (the Eridanus cloud, with Wilmer et al. (1989) describing it as smaller galaxy) in the model with the larger disk mass of a cluster made up of three-four subclumps and Omar & Dwarakanath (2005) describing it as a loose group at an evolutionary stage intermediate to that of Ursa Major and Fornax. The X-ray information suggests that these are distinct systems and a possible candidate for a supergroup. To determine what this structure is, it was impor- tant to define which galaxies are associated with which structure. We obtained galaxy positions, velocities, and K-band magnitudes from the 6dF Galaxy Survey (6dFGS) DR2 (Jones et al. 2005a) and, as 6dFGS is not yet complete, NED1, catalogues for a circle of radius 15◦ (6.6 × 6.6 Mpc), centred on the position of NGC 1332, in the velocity range 500–2500 km s−1, obtaining 513 galaxies. We then used the ‘friends-of-friends’ method (FOF; Huchra & Geller 1982) to determine which galaxy was associated with which structure.

1 The NASA/IPAC Extragalactic Database (NED) is operated by the Jet Figure 1 The final distributions of gas (cyan) and new stars Propulsion Laboratory, California Institute ofTechnology, under contract (magenta) after tidal stripping of Hi gas. with NASA. 328 V. A. Kilborn et al.

Table 1. HOPCAT optical galaxy matching results Ϫ18 NGC 1407 group

Optical galaxy match category Percentage Ϫ 20 Optically matched with velocity NGC 1332 group Single match 42 Compact group member 16 Ϫ22 Optically matched with no velocity

Dec (degrees) Single match 20 Compact group member 6 No matches Ϫ 24 No match; multiple galaxies present 11 Blank field; no visible galaxy 5 58 56 54 52 50 RA (degrees)

Figure 2 6dFGS & NED galaxies with v<2500 km s−1 and the survey covers the whole of the southern sky up to ◦ FOF group finder output. Dec =+2 . The HIPASS Optical Catalogue (HOPCAT) is based on the HIPASS catalogue (HICAT; Meyer et al. 2004; Zwaan et al. 2004), which represents the largest Hi Figure 2 shows that the FOF algorithm finds three dis- selected catalogue at this time. tinct groups. The NGC 1407 and NGC 1332 groups are One motivation for the HIPASS survey is to investi- centred on the large ellipticals of the same name and their gate the existence of dark galaxies. For the purposes of X-ray centroids. The Eridanus group is not centred on any this paper we define a dark galaxy as any Hi source that large elliptical or the X-rays associated with NGC 1395. contains gas (and dark matter) but no detectable stars, and It is possible to determine the dynamical parameters is sufficiently far away from other galaxies, groups or clus- of these structures. Consistent with its X-ray informa- ters such that a tidal origin can be excluded, i.e. isolated tion, NGC 1407 is a massive group of 19 galaxies with a dark galaxies. For further details on HOPCAT see Doyle − high velocity dispersion (σ = 372 km s 1) and a low cross- et al. (2005). = −1 ing time (tc 0.03 H0 ), consistent with it being virial- To identify the optical galaxy match for each HICAT ized. The NGC 1332 group has fewer galaxies (N = 10) detection we use 15 × 15 arcmin SuperCOSMOS images −1 but forms a compact structure (σ = 163 km s , tc = (Hambly et al. 2001a; Hambly, Irwin, & MacGillivray −1 0.04 H0 ). The Eridanus group is made up of more 2001b; Hambly et al. 2001c) to allow for uncertainity galaxies (N = 31) but is a much looser, irregular struc- in the original HIPASS position. However this creates a ture. This is echoed in its low velocity dispersion (σ = problem with matching the correct optical galaxy with −1 = −1 156 km s ) and high crossing time (tc 0.06 H0 ). its corresponding original HIPASS detection due to mul- The NGC 1407 group is a massive group at a late stage tiple galaxies present in the majority of the images. To in group evolution, whilst Eridanus and NGC 1332 appear overcome this problem we cross-check HICAT veloc- to each be at an early stage of their evolution, suggesting ity measurements with optical and high-resolution radio that this is likely to be a supergroup. velocities from NED and the 6dFGS (Wakamatsu et al. We also examined the colours of the galaxies with 2003) to validate the galaxy-match choice. respect to their density. The density was measured as the An automated visual interactive program (adric), projected surface density of the five nearest neighbours where images centred on each HIPASS source position − to each galaxy within ±1000 km s 1. B-band magnitudes are viewed by several people to minimize galaxy selection were obtained from HyperLEDA2 in order to calculate bias, has been developed. adric utilizes the SuperCOS B − K colours. The galaxies are redder in denser environ- images, SExtractor image analysis (Bertin & Arnouts ments. The density at which this occurs is of the order 1996), and NED and 6dFGS velocities for cross-checking, − ∼2 Mpc 2, equivalent to that on the outskirts of this to reliably match HICAT Hi sources with their optical structure. This indicates that changes in galaxy properties counterparts. The results for the matching process are are occuring at densities equivalent to those of galaxies shown in Table 1. HOPCATcontains optical galaxy choice currently infalling. categories that not only describe the type of optical galaxy match but the quality of the resulting magnitudes. 3 HOPCAT and Isolated Dark Galaxies The selection criteria used to search for isolated dark galaxy candidates are extinction cut-off and the blank field Marianne T. Doyle, Michael J. Drinkwater, category. We use an ABj extinction cut at 1 mag because, and David J. Rohde beyond this extinction, optically faint galaxies will be dust The blind Hi Parkes All-Sky Survey (HIPASS), on the obscured.

Parkes Radio Telescope, was completed in 2000. This A total of 3692 galaxies have an ABj extinction less than 1 mag, with only 13 galaxies also in the blank field 2 http://leda.univ-lyon1.fr/ category. From these, 11 are found to be in over-crowded Galaxy Groups: Proceedings from a Swinburne University Workshop 329

fields. One object on close inspection does have a faint for us to explore not only higher density early-type envi- optical counterpart. The final dark galaxy candidate has ronments such as Dorado and NGC 1407, but also a low recently been confirmed as a false Hi detection. Our con- density spiral environment such as NGC 0681. One of our clusion is that from the 4315 Hi radio detections in HICAT selected groups, NGC 4038, includes the famous Anten- no isolated dark galaxies have been found. nae system of colliding galaxies, and is very important To conclude, we present the largest catalogue to date in understanding UCDs as supermassive star clusters in of optical counterparts for Hi radio-selected galaxies, interacting systems. HOPCAT. Of the 4315 Hi radio-detected sources from the We have carried out a search for UCDs using the 2dF HIPASS catalogue, we find optical counterparts for 3618 spectrograph on the Anglo-Australian Telescope (AAT) (84 per cent) galaxies. Isolated ‘dark galaxy’ candidates in each group in a single two-degree field centred on the are investigated using an extinction cut of ABj < 1 mag group centre of mass. To improve our chances of detecting and the ‘blank field’ optical galaxy match category. We UCDs in the limited observing time, we defined a subsam- conclude that there are no isolated optically dark galaxies ple of objects looking similar to bright Fornax and Virgo present within the limits of the HIPASS survey. UCDs: They are unresolved in Schmidt photographic plates and have approximately the same luminosity range − − 4 Ultra-Compact Dwarf Galaxies in as bright Fornax andVirgo UCDs ( 14

The motivation for this project comes from the fact (kmVelocity s that groups are poorly studied relative to clusters and yet 1000 contain most galaxies in the Universe. Recent large sur- veys like 2dF and SDSS have shown that star formation 0 suppression occurs at group-like densities. 0 0.5 1 1.5 2 2.5 3 But many questions remain: What is physical mecha- Distance from group centre (Mpc) nism? What is the timescale? A key property of galaxy groups is the virial radius. Figure 3 Velocity versus distance plot for NGC 5044 group galax- This can be calculated from the velocity dispersion of the ies. The solid line shows the mean group velocity, and the dashed ∼ galaxies or from the group X-ray temperature (Osmond & lines the velocity dispersion. The virial radius occurs at 0.75 Mpc. Ponman 2004). The combination of Hi mapping (e.g. McKay et al. 2004) and 6dFGS spectra have increased dramatically the M =−22.0, L /L = 32 and over 2B mag brighter than number of confirmed group members. We are now in a B X B the second ranked galaxy. It has an old stellar popula- good position to compare the virial radii derived from tion and is featureless, suggesting any merger/collapse velocity dispersion with those from X-ray temperatures. happened along time ago. Figure 3 shows an example for the NGC 5044 group (a massive group). The virial radius occurs at ∼0.75 Mpc for both methods. Note the number of group galaxies beyond 6 Where is the Hi in Galaxy Groups? one virial radius. To better understand processes in groups we are also Virginia A. Kilborn investigating isolated elliptical galaxies, to act as a control Deep wide-field Hi observations of galaxy groups can sample. We selected galaxies with T<−3, D<130 Mpc, help us understand the evolution and formation of galax- B<14, and no neighbours within 700 km s−1, 0.67 Mpc ies within groups. The existence of Hi tails and bridges projected radius, and 2B mag (see Reda et al. 2004 for provide evidence of previous interactions, and the Hi con- details). The aim is to understand their formation, namely tent of galaxies can tell us about their history. However are they old ‘primordial’ systems, recent mergers, or it is only with the advent of multibeam receivers in the collapsed groups? last ten years that such surveys are now possible. While Our optical imaging study indicated that some iso- HIPASS surveyed the whole southern sky for Hi (Meyer lated galaxies appear relatively undisturbed while others et al. 2004), its sensitivity for detecting low mass galaxies have much ‘fine structure’ such as plumes and shells indi- was limited to the local volume. Thus to more thoroughly cating a recent merger. Most isolated galaxies obey the explore the Hi content of groups, we have conducted fur- cluster elliptical colour–magnitude and fundamental plane ther wide-field Hi observations of 16 groups in the Group relations however again there were a few notable excep- Evolution Multiwavelength Study (GEMS; see Forbes, tions. These deviant galaxies tended to be ones with fine this meeting, for a fuller description). The Hi observa- structure and evidence for young central stellar popula- tions were made at the Parkes radiotelescope, using the tions (Reda et al. 2005). A future work will examine the multibeam receiver. An area of ∼5.5 × 5.5 degrees was radial kinematics for isolated ellipticals (Hau & Forbes, observed around each group centre. The observations are in preparation). about twice as deep as HIPASS, and the velocity resolu- The crossing times for some groups are much shorter tion is ten times that of HIPASS. The sensitivity of the Hi than the age of the Universe. So if some groups have survey is ∼1– 5 × 108 M, depending on the distance to collapsed already, how would they appear? Simulations the group. suggest that multiple mergers will result in a large, We detected 210 Hi sources in the 16 datacubes. The relatively isolated with a group-like X- positions of the Hi sources were cross correlated with both ray halo (the hot gas cooling time is generally longer than the NED database, and the 6dFGRS (Jones et al. 2005a). the age of the Universe). Do any of our isolated galaxies We found the majority of the Hi detections corresponded look like collapsed groups? with one or more previously optically catalogued galaxies, Our typical isolated galaxy has MB =−20.5 (and often with known .Wealso found 14 Hi detections LX/LB = 30), which is too low luminosity to be an entire that do not match with previously catalogued galaxies. group. One potential collapsed group is NGC 1132 with A further eight Hi detections provided the redshift for Galaxy Groups: Proceedings from a Swinburne University Workshop 331 previously catalogued galaxies, thus across the 16 groups 7 Neutral Hydrogen in the M83 the total number of new group members was 22. In some Bärbel S. Koribalski groups, there has been an addition of up to 50% of known group members, when adding in the new Hi and 6dFGS ATCA and Parkes Hi line observations of the grand redshifts. design spiral galaxy M83 (NGC 5236, HIPASS J1337-29) Optical images of the 14 new Hi sources were obtained reveal a very extended Hi distribution with a diameter of ∼ from the digital sky survey (DSS). The majority of the 100 kpc, several times larger than the optical Holmberg new detections could be tentatively matched with a visible radius. While the inner disk of M83 rotates remarkably optical galaxy, however there were four cases where there regularly, the Hi gas dynamics appear increasingly pecu- was no obvious optical counterpart on the DSS images.We liar towards the outer regions which show clear signs of have obtained high resolution Hi images from the ATCA tidal disruption. The most prominent tidal features are the for one of these Hi sources, GEMS_NGC 3783_12, in the asymmetric outer Hi arm bending towards the east of M83 NGC 3783 galaxy group. The structure of this object is and a spectacular stellar stream, consisting of mainly old reminiscent of a tidal Hi cloud, as it is irregular and has a stars, to the north. M83 is surrounded by numerous dwarf tail. The formation of this object is uncertain, and we find galaxies and, given its large dynamical mass, is likely to no Hi bridge to any optical galaxy to the sensitivity of the attract and accrete them in regular intervals. datacube (∼5 × 108 M). The cloud is ∼ 50 kpc in size, M83 is a late-type spiral galaxy with an unusually large ∼ = and has an Hi mass of ∼109 M. There is one further Hi Hi envelope of 100 kpc (for D 4.5 Mpc), at least five cloud candidate in the NGC 3783 group, and we have two times larger than its optical Holmberg diameter. It is a more candidates in the NGC 5044 galaxy group. ATCA member of the nearby CentaurusA group which appears observations of these Hi cloud candidates are scheduled to consist of two subgroupings, one around M83 and the for 2006 January. other around CenA (NGC 5128). Figure 4 shows the Wehave begun to investigate the Hi content of the group deep Parkes multibeam Hi data of the M83 subgrouping; galaxies, and whether Hi deficiencies are seen similar to no low-surface brightness Hi extensions were detected that in clusters (Solanes et al. 2001) and compact groups between M83 and its neighbouring galaxies down to an (Verdes-Montenegro et al. 2001). We have looked at one group in depth, NGC 1566. We find the total Hi content of the group is consistent with the optical members of the 30Ј group. However two galaxies are about ten times more Hi ESO444–G084 deficient than would be expected from their optical type Ϫ28°00Ј and size (Kilborn et al. 2005). There is no diffuse X-ray emission in this galaxy group, and the Hi deficient galaxies 30Ј do not lie near the centre of the group, thus we expect the IC4316 Hi deficiency is caused by tidal stripping of the gas by Ϫ ° Ј either the group potential or other galaxies. 29 00 UGCA365 We can place some limits on the existence of intra- group neutral hydrogen, down to the sensitivity of our 30Ј NGC 5264 survey. Such intra-group Hi is rare, and makes up less than 2% of the Hi detections we found in our survey. However, Ϫ30°00Ј our survey cannot place limits on the number of low-mass Dec (J2000) Hi galaxies in groups. 30Ј So to answer the question ‘Where is the Hi in galaxy M83 groups’ we can say that to the limit of ∼1– 5 × 108 M,it Ϫ31°00Ј is mostly contained within galaxies, and Hi clouds, tidal NGC 5253 tails, and Hi bridges add only a small part to the total Hi 30Ј content of galaxy groups. I acknowledge the GEMS team for the general survey, Ϫ32°00Ј and thank Heath Jones for supplying the 6dF data.

13h45m 42m 39m 36m 33m 30m RA (J2000) Table 2. GEMS Hi survey results Figure 4 Hi distribution of M83 and its surroundings based on Number of groups observed 16 data from the Parkes HIDEEP survey. The contour levels on the −1 −1 Total galaxies detected 210 left are 0.5, 1, 2, 4, 8, 16, 32, 64, 128, and 256 Jy beam km s Previously uncatalogued 14 (the first contour corresponds to an Hi column density of 18 −2 New redshifts 8 ∼0.7 × 10 atoms cm ). We measure a peak Hi column density of 5 × 1020 atoms cm−2 occurs at the centre of M83. The gridded Confirmed Hi clouds 1  Total new group members 22 Parkes beam of 15.5 is indicated at the bottom left. Galaxies detected in Hi are marked with a cross. 332 V. A. Kilborn et al.

Hi column density limit of ∼1018 atoms cm−2 (assuming near-neighbours, projected surface densities, group mem- the Hi gas fills the beam). bership, correlation with the remaining 2dFGRS galaxies, The Hi distribution of M83, as revealed with the Aus- and the luminosity function. The analyses are also per- tralia Telescope Compact Array (Koribalski et al., in formed on a random sample of 2dFGRS galaxies which preparation) is most remarkable. No longer does this act as a benchmark for comparison to our starburst sample. grand-design spiral look regular and undisturbed. The Hi The near neighbour (NN) analysis involves finding the maps show streamers, irregular enhancements, an asym- transverse separation of the starburst and its nearest faint metric tidal arm, diffuse emission, and a thoroughly and bright neighbours. A faint NN is defined as a galaxy ∗ ∗ twisted velocity field, much in contrast to its regular with b (z) + 1

9 The Subhalo Mass Function in Galaxy Clusters: ( 20

N M : 15.00 M : 14.79 M : 14.89 Another Success for the Cold Dark 15 vir vir vir Matter Model? 10 Chris Power, Brad K. Gibson 5 20 One of the defining characteristics of the Cold Dark Matter M : 14.97 M : 14.81 Average 15 vir vir (CDM) model is the hierarchical manner in which massive LCDM 10 systems are assembled through the merging and accretion NS04 5 of less massive progenitors. However, this merging pro- 0 cess is incomplete and the remnant subhaloes constitute 11 12 13 11 12 13 11 12 13 ∼ Ϫ1 10% of the virial mass in a typical CDM halo, following Log10 M [h M ] a power-law mass function, N(M) ∝ M−0.8 with M the subhalo mass (Gao et al. 2004a). Figure 5 Comparison of simulated mass functions with averaged In a recent study, Natarajan & Springel (2004) have observational mass functions from Natarajan & Springel (2004). combined strong and weak lensing observations of a sample of five rich clusters taken from Natarajan et al. be directly compared (Gao et al. 2004b). Therefore, it is (2004) to construct an averaged mass function of the unclear that the discrepancy is significant at all without a host dark matter halos of the most massive cluster galax- proper treatment of the galaxy formation process. This is ies. They compared this averaged mass function with very much work in progress! the subhalo mass function derived from the cluster sim- We warmly thank the Virgo Consortium, Alexander ulations of Springel et al. (2001), who studied a single Knebe, and Stuart Gill for the use of their simulated rich CDM cluster simulated at progressively higher mass clusters. and force resolution; the best resolved contained ∼20 million particles within the virial radius. The authors References claimed good agreement between the amplitudes and Abadi, M. G., Moore, B., & Bower, R. G. 1999, MNRAS, 308, 947 slopes of the respective mass functions over the mass range Baldwin, J., Phillips, M., & Terlevich, R. 1981, PASP, 93, 5 Bertin, E., & Arnouts, S. 1996, A&AS, 117, 393 1011 ≤ M/M ≤ 1012.5, and concluded that such concor- Bekki, K., Forbes, D. A., Beasley, M. A., & Couch, W. J. 2002, dance provided further observational evidence in favour MNRAS, 335, 1176 of the CDM model. Bekki, K., Couch, W. J., Drinkwater, M. J., & Shioya, Y. 2003, We have investigated the significance of this claim MNRAS, 344, 399 by analyzing a sample of eleven rich cluster mass CDM Bekki, K., Koribalski, B. S., Ryder, S. D., & Couch, W. J. 2005a, haloes, each resolved with ∼1 million particles. In contrast MNRAS, 357, L21 Bekki, K., Koribalski, B. S., & Kilborn, V. A. 2005b, MNRAS, in to Natarajan et al. (2004), we find that subhalo mass func- press (astro-ph/0505580) tions derived from simulated cluster haloes disagree with Blake, C. A., et al. 2004, MNRAS, 355, 713 the mass functions derived from lensing observations in a Blumenthal, G. R., Faber, S. M., Primack, J. R., & Rees, M. J. 1984, systematic manner, as can be seen in Figure 5. If we con- Natur, 311, 517 sider subhaloes within a projected radius of 500h−1 kpc Colless, M., et al. 2001, MNRAS, 328, 1039 3 Colless, M., et al. 2003, astro-ph/0306581 from the centre of the simulated cluster , we find that the Doyle, M. T., et al. 2005, MNRAS, 361, 34 number of a given mass (solid histogram) underestimates Drinkwater, M. J., Jones, J. B., Gregg, M. D., & Phillipps, S. 2000, the number of cluster galaxy haloes (hatched histogram) PASA, 17, 227 by a factor of about three over the range of overlap.Agree- Drinkwater, M. J., Gregg, M. D., Couch, W. J., et al. 2004, PASA, ment between observationally derived and simulated mass 21, 375 Eke, V. R., et al. 2004a, MNRAS, 348, 866 functions is possible only if we consider all subhaloes Eke, V. R., et al. 2004b, MNRAS, 355, 769 within the virial radius of the host (dotted histogram). Gao, L., White, S. D. M., Jenkins, A., Stoehr, F., & Springel, V. How significant is this discrepancy for the CDM 2004a, MNRAS, 355, 819 model? Recent studies have demonstrated that galax- Gao, L., De Lucia, G., White, S. D. M., & Jenkins, A. 2004b, ies and subhaloes represent different populations with MNRAS, 352, L1 Gomez, P. L., et al. 2003, ApJ, 584, 210 distinct spatial and kinematic distributions, and should not Huchra, J. P., & Geller, M. J. 1982, ApJ, 257, 423 Hambly, N. C., Davenhall, A. C., Irwin, M. J., & MacGillivray, H. T. 3 Comparable to Natarajan et al. (2004). 2001a, MNRAS, 326, 1315 334 V. A. Kilborn et al.

Hambly, N. C., Irwin, M. J., & MacGillivray, H. T. 2001b, MNRAS, Reda, F., Forbes, D., Beasley, M., O’Sullivan, E., & Goudfrooij, P. 326, 1295 2004, MNRAS, 354, 851 Hambly, N. C., et al. 2001c, MNRAS, 326, 1279 Reda, F., Forbes, D., & Hau, G. 2005, MNRAS, in press Hilker, M., Infante, L., Viera, G., et al. 1999, A&AS, 134, 75 Ryan-Weber, E.,Webster, R., & Bekki, K. 2003, in ‘The IGM/Galaxy Jones, D. H., Saunders, W., Read, M., & Colless, M. 2005a, PASA, Connection: The Distribution of Baryons at z = 0’, ASSL Conf. 22, 277 Proc.Vol.281, eds. J. L. Rosenberg, & M. E. Putman (Dordrecht: Jones, J. B., Drinkwater, M. J., Jurek, R., et al. 2005b, MNRAS, in Kluwer), p. 223 press Schmidt, M. 1959, ApJ, 129, 243 Kauffmann, G., et al. 2003, MNRAS, 346, 1055 Solanes, J. M., Manrique, A., García-Gómez, C., González- Kenney, J. D. P., van Gorkum, J. H., & Vollmer, B. 2004, AJ, 127, Casado, G., Giovanelli, R., & Haynes, M. P. 2001, ApJ, 548, 97 3361 Springel, V., White, S. D. M., Tormen, G., & Kauffmann, G. 2001, Kennicutt, R. C., Jr. 1998, ARA&A, 36, 189 MNRAS, 328, 726 Kilborn, V. A., Koribalski, B. S., Forbes, D. A., Barnes, D. G., & Stevens, J. B., Webster, R. L., Barnes, D. G., Pisano, D. J., & Musgrave, R. C. 2005, MNRAS, 356, 77 Drinkwater, M. J. 2004, PASA, 21, 318 Klypin, A., Kravtsov, A. V., Valenzuela, O., & Prada, F. 1999, ApJ, Toomre, A., & Toomre, J. 1972, ApJ, 178, 623 522, 82 Tully, R. B. 1987, ApJ, 321, 280 Kodama, T., Smail, I., Nakata, F., Okamura, S., & Bower, R. G. Verdes-Montenegro, L., Yun, M. S., Williams, B. A., 2001, ApJL, 562, 9 Huchtmeier, W. K., Del Olmo, A., & Perea, J. 2001, A&A, Lewis, I., et al. 2002, MNRAS, 334, 673 377, 812 McKay, N., et al. 2004, MNRAS, 352, 1121 Vollmer, B., Cayatte, V., Balkowski, C., & Duschl, W. J. 2001, ApJ, Meyer, M. J., et al. 2004, MNRAS, 350, 1195 561, 708 Mieske, S., Hilker, M., & Infante, L. 2004, A&A, 418, 445 Vollmer, B. 2003, A&A, 398, 525 Moore, B., Ghigna, S., Governato, F., Lake, G., Quinn, T., Wakamatsu, K., Colless, M., Jarrett, T., Parker, Q., Saunders, W., & Stadel, J., & Tozzi, P. 1999, ApJ, 524, 19 Watson, F. 2003, ASPC, 289, 97 Natarajan, P., Kneib, J., Smail, I., & Ellis, R. 2004, astro-ph/0411426 Willmer, C. N. A., Focardi, P., da Costa, L. N., & Pellegrini, P. S. Natarajan, P., & Springel, V. 2004, ApJL, 617, 13 1989, AJ, 98, 1531 Norberg, P., et al. 2002, MNRAS, 336, 907 Yahagi, H., & Bekki, K. 2005, ApJL, submitted Omar, A., & Dwarakanath, K. S. 2005, JApA, 26, 1 Yun, M. S., Ho, P. T. P., & Lo, K. Y. 1994, Natur, 372, 530 Oosterloo, T., & van Gorkum, J. 2005, A&A, 437, 190 Zabludoff, A. I., & Mulchaey, J. S. 1998, ApJ, 496, 39 Osmond, J. P. F., & Ponman, T. J. 2004, MNRAS, 350, 1511 Zwaan, M. A., et al. 2004, MNRAS, 350, 1210