arXiv:astro-ph/0511050v1 1 Nov 2005 o.Nt .Ato.Soc. Astron. R. Not. Mon. eevdJl 2005 July Received Raychaudhury Somak and Porter C. Scott Sloan and Survey surveys Redshift Sky Galaxy Digital 2dF of the filament remarkable in a galaxies Supercluster: Pisces-Cetus The hsipista htteUies snthmgnoson homogeneous not 1998). is al. Universe et below the Jenkins found scales 1994; that be al. that et to implies Einasto seem This 1993; fil- galaxies al. of most et network which (Zucca a in of voids presence and the aments revealed structure have large-scale Universe the the of simulations and Observations INTRODUCTION 1 eg,Bgh ta.20;Bade l 03 ilre al. et Miller ( 2003; clusters the rich than individual al. larger of et much radii Brand are virial that 2002; structures Large-scale claimed al. 2004). been et have dis- Bagchi these the than (e.g., though larger far galaxies, structures of of of scale covery structures the is common which largest 2004), Sahni the & Sheth, Shandarin, 2004;
hspoet xed.Te r nesnilto ostudy to ⋆ scales tool what essential to an clear are not They is extends. it property though this Universe, the in literature. tems the paper. in this entities in various terminology term denote this this adopt to though We used even galaxies, been of has superclusters called ten c colo hsc n srnm,Uiest fBirmingham, of University Astronomy, and Physics of School -al [email protected] somak@star.sr.bham.ac [email protected]; E-mail: 05RAS 2005 uecutr r h ags rvttoal on sys- bound gravitationally largest the are Superclusters ∼ 100 h 70 − 1 p Baawj hva,&Sheth & Bhavsar, (Bharadwaj, Mpc 000 –1(05 rne oebr21 M L (MN 2018 November 7 Printed (2005) 1–11 , h icsCtssprlse (redshift supercluster Pisces-Cetus The ABSTRACT iprin n pia ufc rgtespols h oe ii t constitu limit the of lower 5 mass A be The virial total including to profiles. the clusters, to cluster brightness amounting the filaments, Abell surface of these 5 masses optical virial m and and th the two 11 estimate dispersions in identify and of suppressed We respectively, poor 85 consisting is studies. the Abell previous supercluster, formation among with this star higher consistent in contrast, be is to In which appears ones. clusters, effect richer super This the the We field. in groups. in the groups and and to in sky belonging clusters those the galaxies constituent to in on its formation supercluster of star distribution the enhanced of the of filament map ong extent remarkable the and the and a explore survey redshift We as galaxy Survey. two-degree-field emerges the and of edges galaxies, of superclusters tutr fteUies;Xry:glxe:clusters. galaxies: X-rays: Universe; the of structure n 19600 and e words: Key δ M > ∼ > 4 10 . . 2 and 7 × h 70 − h 10 1 70 − p)aeof- are Mpc) 3 15 δ aais lses eea;Glxe:eouin omlg:large- Cosmology: evolution; Galaxies: general; clusters: Galaxies: Mpc M h 70 − > 1 3 .uk M 1 epciey hscrepnst asoedniiso Ω of overdensities mass to corresponds This respectively. . ⊙ o h w lmnsmkn ptesupercluster. the up making filaments two the for 3 of n 6 and imnhmB52T UK 2TT, B15 Birmingham . ⋆ oihswihd o aeti emtia etr into feature geometrical this take 1994; not al. consideration. friends-of-friends do et or Einasto which percolation using gorithms (e.g., 1993) al. constructed et cat Zucca usually supercluster are Nevertheless, pre-eminence logues the structures. 2002 show al. filamentary 2004) et of Donovan Bagchi in- & (e.g., (e.g., of Barrett, superclusters 2dfGRS Ebeling, studies of or the systems 2004) the like dividual Hawkrigg as surveys & well Drinkwater, large-volume Pimbblet, as Pogosyan 2005), of & Connolly Kofman, analyses & Bond, Krughoff, (e.g., Colberg, 1996; structure evolutio large-scale the of of simulations Numerical structures. mentary (e.g., 1988). models Bahcall be- 2002; formation Plionis discriminating structure & in Basilakos, and Kolokotronis, useful on matter light being dark to thus mass tween scales, of mass ratio largest high-end and the the systems the collapsing Super- quantifying of tractable. in function analytically useful in be be evolved can can clusters and always regime are have These linear they 2000). the since Bahcall (e.g., interesting Universe ri particularly the given have in structure that to perturbations density largest-scale the 0 × 10 uecutr r xetdt epeoiatyfila- predominantly be to expected are Superclusters 15 h 70 − z 1 ∼ M A 0 T . E ⊙ 6 soeo h ihs ery( nearby richest the of one is 06) tl l v2.2) file style X vrvlmso lot6750 almost of volumes over igSonDgtlSky Digital Sloan oing n lses un out turns clusters, ent rmtervelocity their from s aaisi rich in galaxies e nrdhf space, redshift in lse compared cluster n vdneof evidence find h assof masses the o aaisa the at galaxies rgop than groups er jrfilaments ajor el13and 133 bell h 70 − / < z Ω 3 crit Mpc scale 0 . al- 1) se a- ≡ n 3 ; 2 Porter and Raychaudhury
Figure 1. (a, Left) Clusters of galaxies belonging to Pisces-Cetus supercluster. The minimal spanning tree used to identify the supercluster (Raychaudhury et al. in preparation) is shown connecting the clusters. Edges in the tree with length 20 h−1
Recently, Raychaudhury et al. (in preparation) have In this paper, we seek to characterise the extent of the compiled a catalogue of supercluster-size structures from a supercluster in 2 from the 2dFGRS and SDSS surveys, and § redshift survey of Abell clusters (z < 0.1) using the min- from the SuperCosmos Sky Survey (Hambly et al. 2001). In imal spanning tree statistic, which contains all the infor- 3 we calculate mean velocities and velocity dispersions for § mation used in a percolation analysis (Bhavsar & Splinter the constituent clusters, while in 4 we look at evidence § 1996), and is particularly sensitive to filaments. In this of enhanced star formation in galaxies belonging to groups catalogue, only one of the long filaments (containing 8 or along the filaments of the supercluster, as an indicator of more Abell clusters) is in the 2dF galaxy redshift survey galactic interaction in superclusters. In 5 we calculate virial § (Colless et al. 2003) [2dFGRS] or Sloan Digital Sky Survey masses of the constituent clusters based on their optical lu- (Abazajian et al. 2005) [SDSS] regions. In fact, it straddles minosity profiles and velocity dispersion, to estimate a lower both of them. limit to the mass overdensity in the supercluster. We sum- More than a decade ago, in a series of papers, Tully marise our major conclusions in 6. Detailed notes and X- § (1986, 1987, 1988) pointed out the Pisces-Cetus superclus- ray images of the constituent clusters can be found in the −1 ter to be one of the five richest systems of clusters of Appendix. Throughout this paper we use H0 = 70 km s −1 galaxies in the z < 0.1 Universe, and speculated that it Mpc . −1 could be as large as 300 h70 Mpc across, and might include the Local Group. Subsequent work (e.g., Zucca et al. 1993; Einasto et al. 1994) has pointed out interesting features of 2 CLUSTERS OF GALAXIES BELONGING TO significantly smaller substructures of this remarkable web THE PISCES-CETUS SUPERCLUSTER of galaxies, groups and clusters, and at least one unusually large void (Burns et al. 1988). The clusters comprising the Pisces-Cetus supercluster from In the supercluster catalogue of Raychaudhury et al. the MST analysis of Raychaudhury et al. (in preparation) (in preparation), the Pisces-Cetus supercluster occupies a form the core of our sample, and are listed in Table 1. Unlike region of about 700 sq. deg. of the (equatorial) southern sky. MST analyses of galaxy catalogues, where completeness is While using the Minimal Spanning Tree (MST) in construct- an important consideration, the supercluster list obtained −1 ing this catalogue, if a maximum edge-length of 20 h100 Mpc in this catalogue is considered reasonably complete since it is used, the system breaks up into two filaments (A, B) of 11 arises from an analysis of rich Abell clusters in the volume and 5 clusters respectively, at mean redshifts of 0.0625 and closer than z < 0.1 (Galactic latitude b > 20◦). | | 0.0545 respectively. If, however, adjacent clusters on the tree To the member list of the Pisces-Cetus Superclus- −1 are allowed to be separated by a distance of up to 25 h100 ter from the above list, We added a few clusters, in the Mpc, the supercluster becomes a long filament of 19 clus- same volume as the filament, from the Supplementary ters at a mean redshift of 0.0591, making the supercluster Abell cluster catalogue (Abell, Corwin, & Olowin 1989) and the third richest in the catalogue (following the Shapley and the Edinburgh-Durham Southern Galaxy Cluster Catalogue Horologium-Reticulum superclusters). This is illustrated in (Lumsden et al. 1997). These extra clusters were taken as Fig. 1a, where the clusters excluded by the former consider- potential members of the supercluster if their redshift as de- ation are shown to be connected to the main tree by dashed rived in 3 was found to be consistent with the other mem- § lines. bers.
c 2005 RAS, MNRAS 000, 1–11 The Pisces-Cetus Supercluster 3
ulations where groups are seen to preferentially form on fil- aments and move along filaments to fall into clusters. For comparison, Fig. 1b shows all galaxies with spectroscopic redshifts from the SDSS and 2dFGRS catalogues, where only galaxies within 1000 km/s of the upper and lower redshift bounds of the supercluster have been plotted. Overdensities of galaxies are clearly visible in the regions around the richer clusters and the filamentary structure of the supercluster is also visible. Figs. 2 show the same galaxies, but using the redshift information described below. A combined look at Fig. 1 and Fig. 2 reveals that the supercluster is made of two main fil- aments. From Figs. 2a, which plots redshift against Right Ascension for the groups and clusters in the region, it is evident that the filament containing A85 is at a lower red- shift than the longer filament, which has A133 at the near end and A2716 at the more distant end. The poorer clusters E0332 and S1155 appear to be in a clump separated from both of the main filaments, with a filament of groups joining these three clusters to the principal filament.
3 SPECTROSCOPIC DATA Positions from the NASA/IPAC Extragalactic Database (NED) were taken as a starting point for the cluster cen- tres, and then they were located on UK Schmidt Telescope (UKST) BJ plates. The NED positions, particularly those from Abell, Corwin, & Olowin (1989), were largely inade- quate. Therefore, centres were determined from available X- ray data whenever possible (see Appendix), and from a vi- sual inspection of the appropriate UKST plates otherwise. Our adopted coordinates can be found in Table 1. Redshifts were taken from all avail- able sources including NED, ZCAT (http://cfa-www.harvard.edu/ huchra/zcat/), the first ∼ data release of the 6dF Galaxy Survey Database (Jones et al. 2004) [6dFGS DR1], the 2dF galaxy red- shift survey (Colless et al. 2003) and the second release (DR2) of the Sloan Digital Sky Survey (Abazajian et al. 2005). The galaxies were initially extracted from these −1 catalogues within an Abell radius (2.14 h70 Mpc) of the new cluster centres and filtered for repeated galaxies. Figure 2. (a,top) Clusters of galaxies belonging to Pisces-Cetus supercluster, covering the declination range −32◦ <δ < −9◦. The 2PIGG groups of galaxies are shown as little circles, but only in 3.1 Mean redshift and velocity dispersion the declination range covered by the 2dFGRS (−32◦ <δ < −25◦). (b, bottom) Galaxies with redshifts (mostly from 2dFGRS, see §3) Mean recessional velocities for clusters that have redshifts ◦ −1 in the Pisces-Cetus supercluster region (Declination range −32 < available for > 4 galaxies, within the Abell radius of 2.14 h70 δ < −20◦). The mean positions of the clusters in the supercluster Mpc, were calculated from the compilation described above, are marked. The 2dFGRS only goes up to δ < −25◦, which is using an iterative process clipping 3σ about the median, and why some of the clusters don’t seem to be associated with obvious finally calculating the mean. If > 12 galaxies with known concentrations. redshift remained, the velocity dispersion of the cluster and its error were calculated, using a method that takes into ac- count observational errors on individual redshift measures A plot of these clusters within the supercluster region (Danese, de Zotti, & di Tullio 1980). The resulting mean ve- and redshift range, including the connecting MST, can be locities, velocity dispersions and the number of cluster mem- seen in Fig. 1a. Also plotted are the groups of galaxies bers used to calculate them can be found in Table 1 for each in the region of the supercluster covered by the 2dFGRS cluster belonging to the Pisces-Cetus supercluster. [δ < 25◦] from the 2dFGRS Percolation-inferred Galaxy Velocity histograms of each cluster can be seen in Fig. 3. − Group (2PIGG) catalogue (Eke et al. 2004), which are seen Where more than a peak is obviously present, the velocity to follow a filamentary structure (though not necessarily the dispersion of the peak closest to the mean recessional veloc- tree defined by the MST), as expected from numerical sim- ity of the supercluster has been adopted. However, in some
c 2005 RAS, MNRAS 000, 1–11 4 Porter and Raychaudhury cases, membership remains ambiguous. For example, even 4.1 Rate of star formation: the η parameter though a large number of galaxies were found with mea- To investigate the variation in the rate of star forma- sured redshifts within the Abell radius of A87, it is likely tion rate with local environment, we used the spectral pa- that a large fraction of them are in fact members of A85, −1 rameter η, which is found from a principal components since the centres of these two are only 2 h Mpc apart. ∼ 70 analysis (Madgwick et al. 2002) of the 2dFGRS spectra. Fig. 2b shows all the galaxies with redshifts from 2dF- ◦ ◦ Madgwick et al. (2003) found a tight correlation between GRS and ZCAT in the declination range 32 <δ < 20 . − − η and the equivalent width of the Hα emission line, sug- The diagram shows clearly the filamentary nature of the gesting that η is a measure of relative star formation. Here supercluster and the positions of the Abell clusters on the we will look at galaxies with η > 0, which corresponds to filaments. There are some clusters not obviously associated W (Hα)>4A,˚ so that we can assume that these galaxies have with overdensities– these are in the declination range not a significantly∼ enhanced rate of star formation compared to covered by the 2dFGRS. the mass of the older stellar population.
4 THE INFLUENCE OF THE SUPERCLUSTER 4.2 Samples of galaxies in various environments ENVIRONMENT ON STAR FORMATION IN Seven samples of galaxies were complied from the 2dFGRS GALAXIES (Colless et al. 2003) catalogue, for which the values of η pa- Studies of the environmental effects on galaxy evolution in rameter were extracted. These were: the past have shown that star formation is suppressed in (i) SC-PG: Galaxies (n = 573) belonging to poor groups, the cores of rich clusters (Dressler 1980; Couch & Sharples extracted from the 2PIGG catalogue of groups found in the 1987; Balogh et al. 1998). More recently, from the high- 2dFGRS (Eke et al. 2004), which have 4 6 N 6 10 mem- quality spectra of the 2dFGRS and SDSS archives, the dif- bers, and lying within the Pisces-Cetus supercluster region ference between the properties of star-forming and quiescent (Fig. 1); galaxies have been further characterised, and shown to de- (ii) NSC-PG: Same as above, but for poor 2PIGG groups pend on the local density of galaxies. For instance, studies that do not lie in superclusters (n = 3822). This was (Lewis et al. 2002; G´omez et al. 2003; Balogh et al. 2004) achieved by omitting all galaxies within a liberal radius of all identifying star formation activity with the equivalent width supercluster centres in the catalogue of Raychaudhury et al. of the Hα emission line W (Hα), show that the fraction of (in preparation); ˚ galaxies with W (Hα)>4A steadily declines with increase in (iii) SC-RG: Galaxies (n = 6130) belonging to rich the three-dimensional density of galaxies. Kauffmann et al. 2PIGG groups, which have N > 10 members, lying in the (2004) show similar effects, with the star formation history Pisces-Cetus supercluster region; being quantified from line indices measured elsewhere in the (iv) NSC-RG: Same as above, but for groups not in su- spectrum as well, including the H absorption line and the δ perclusters (defined as in (ii)); strength of the 4000A˚ break. (v) FIELD: All galaxies (n = 44363) in the 2PIGG cata- It is apparent from various other studies (e.g. logue not flagged as being a member of any group; Miles et al. 2004) that there is evidence for an enhanced rate (vi) SC-C: Galaxies (n = 369) from Abell clusters that of galaxy mergers in the not-so-high-density environments of are part of the Pisces-Cetus supercluster, falling within the groups, particularly in the poor ones, where the low velocity 2dFGRS region. dispersion of galaxies favours tidal interaction in a dynami- (vii) NSC-C: Galaxies (n = 552) belonging to Abell clus- cally sluggish system. This leads to, for instance, a depletion ters in the 2dFGRS regions, with similar redshifts, but not in the intermediate-mass galaxies as evident from bimodal belonging to any supercluster in the Raychaudhury et al. luminosity functions. In such systems, we would expect the (in preparation) catalogue. rate of merger-induced star formation to be enhanced with respect to the field, even though star formation is found to Since the rich “groups” and “clusters” are taken from be quenched in high density environments typical of rich different lists, some “cluster” galaxies would be included clusters. within the “rich group” samples, but this small number will Indeed, in the hierarchically merging scenario supported be swamped by the much larger number of group galaxies. by large-volume numerical simulations such as Jenkins et al. Assuming that η > 0 for a galaxy signifies evidence (1998), much of this merger-driven evolution in the history of enhanced star formation, we investigate the fraction of of a galaxy is expected to occur in small groups forming in galaxies that have positive η in the various group and clus- and falling along a supercluster filament towards regions of ter samples, and compare them with the galaxies in the field higher density. In contrast with the behaviour of galaxies in sample. In Fig. 4, we plot cumulative histograms of the frac- the field and among the highest densities in clusters, where tion of galaxies with a value of η between 0 <η<η0 in each one expects to find an abundance of quiescent galaxies, we sample. The histograms are normalised by the total number would expect to find an increase in the fraction of galaxies of galaxies in each sample. with active star formation, induced by interaction, in the To find a formal measure of the difference between the groups that belong to a supercluster filament. Here we in- histograms, we used both the familiar Kolmogorov-Smirnov vestigate this effect by selecting a sample of galaxies that (K-S) statistic and the Kuiper statistic (Press et al. 1992). belong to groups belonging to the part of the Pisces-Cetus Since we are looking for differences in our distributions in supercluster in the 2dFGRS volume, and comparing their the high-η tail, the Kuiper statistic would in principle fare properties with those of galaxies elsewhere. better than the K-S statistic in quantifying differences.
c 2005 RAS, MNRAS 000, 1–11 The Pisces-Cetus Supercluster 5
−1 Figure 3. Radial velocities of galaxies lying within the Abell radius of 2.14 h70 Mpc of each galaxy cluster, belonging to the Pisces-Cetus supercluster, compiled from NED, 6dFGS, 2dFGRS and SDSS. For each cluster, the mean recessional velocity is shown as a dotted line, and for those with >12 redshifts, the velocity dispersion is shown on the upper right.
c 2005 RAS, MNRAS 000, 1–11 6 Porter and Raychaudhury
Figure 4. (Left, a) Cumulative distribution of the spectral parameter η for galaxies with 0 <η<η0, normalised by the total number of galaxies in each sample. A high value of the parameter η indicates a higher width of the Hα emission line, and thus a higher rate of star formation. The thickest solid line refers to the FIELD sample (see section 4.2). The thinnest and medium thick solid histograms refer to poor group galaxies within the supercluster, (SC-PG) and rich group galaxies within the supercluster (SC-RG) respctively. The dashed histograms refer to poor group galaxies not within the supercluster (NSC-PG, thinner line), and to rich group galaxies not within the supercluster (NSC-RG, thicker line). (Right, b) Cumulative distribution of η, as in (a), with the thick histogram being the field sample as in the left panel, the dash-dotted line representing cluster galaxies within the supercluster region (SC-C) and the dotted histogram standing for cluster galaxies not within the supercluster region (NSC-C).
It can be seen in Fig 4a that there is a marked differ- and 0.1 according to the K-S and Kuiper statistics respec- ence, in the distribution of the η parameter, between the tively. group galaxies within and outside of the supercluster. For It can be seen in Fig 4b that the cluster galaxies both poor and rich groups, galaxies in those that belong to have a smaller percentage of their galaxies in the higher the supercluster show evidence of enhanced star formation. η values than the field sample, implying star formation is For the SC-PG and FIELD samples, the K-S and Kuiper suppressed in the cluster environment in agreement with statistics indicate that the probability of the samples to be Dressler (1980); Couch & Sharples (1987); Balogh et al. drawn from the same parent distribution is 0.01 and 10−8 (1998). However, the cluster galaxies within the superclus- respectively. This would be consistent with the filamentary ter seem to have an even lower star formation rate than structure of the supercluster leading to an enhanced merger the cluster galaxies outside of the supercluster. This is the rate among galaxies belonging to the supercluster filament. opposite trend to that found for the group galaxies. For both SC and NSC groups, galaxies in poorer groups seem to have the higher star formation rates. This would support evidence of more rapid merging within the poorer 5 RADIAL PROFILES OF LIGHT groups, which would have generally lower velocity disper- For all clusters belonging to the Pisces-Cetus supercluster, sions, which would in turn encourage tidal interactions. This listed in Table 1, all galaxies within a radius of 2.14 h−1 is in agreement with the difference found between the lumi- 70 Mpc (equal to the Abell radius, defined as 1.5 h−1 Mpc) nosity functions of galaxies by Miles et al. (2004) in X-ray 100 of our adopted centres were extracted from the Reliable bright and X-ray faint groups. Galaxies catalogue of the SuperCosmos Science archive Even smaller probabilities are obtained when the SC- (Hambly et al. 2001). Galaxies were radially binned, with PG and SC-C samples, i.e. galaxies belonging to poor groups variable radius such that there was an equal number of and to clusters in the supercluster, are compared. On the galaxies (N 10) within each annulus. The luminosity per ∼ other hand, the difference between the SC-PG and SC-RG unit area in each annulus was then calculated to plot the galaxies, however, is marginal, but possibly significant, the surface brightness profiles, an example of which is shown in probability of being drawn from the same distribution is 0.53 Fig. 5. A 10% error is shown associated with each point.
c 2005 RAS, MNRAS 000, 1–11 The Pisces-Cetus Supercluster 7
5.2 Virial mass estimates The mass within a radius r from the centre of the cluster can be derived from the modelled velocity dispersion and spatial density profiles λ (r) (Koranyi & Geller 2000). On assumption of a constant mass-to-light ratio Υ M/L, the ≡ mass density ρ (r) Υ λ (r), following (1), results in the ≡ enclosed mass profile,
3 r r/rc M( The core radius rc has already been found in the fitting of the surface brightness profile above. The quantity δc and c are linked by the requirement that the mean density within r200 should be 200 ρc, i.e., × 200 c3 δc = , (5) 3 [ln(1 + c) c/(1 + c)] − Figure 5. The radial surface brightness profile for galaxies in where c = r200/rc. the Abell 2734, with the best-fit projected NFW model (Eq. 2) Assuming that the virial theorem holds, an estimate of superposed. Member galaxies are binned in annuli around the r200 can be gained from cluster centre, such that each annulus contains an equal number 3 of galaxies. M200rv 2 (1 + z) 3 = 100H0 2 , (6) Mν r200 3σr 2 where σr is the radial velocity dispersion of the cluster and 5.1 NFW Surface brightness profiles rv is the approximate virial radius determined by The radial surface brightness profiles were fit, using χ2 min- π N(N 1) rv − , (7) imisation, with the Navarro, Frenk & White (1996) (NFW) 2 −1 ∼ i i