Issue 
A&A
Volume 551, March 2013



Article Number  A143  
Number of page(s)  7  
Section  Extragalactic astronomy  
DOI  https://doi.org/10.1051/00046361/201220071  
Published online  11 March 2013 
Distribution and evolution of galaxy groups in the Ursa Major supercluster
^{1} Laboratório de Astrofísica Teórica e Observacional Departamento de Ciências Exatas e Tecnológicas Universidade Estadual de Santa Cruz, 45650000 IlhéusBA, Brazil
email: krausefisico@ig.com.br; albr@uesc.br
^{2} Observatório do Valongo, Universidade Federal do Rio de Janeiro, 20080090 Rio de JaneiroRJ, Brazil
email: paal05@gmail.com
Received: 20 July 2012
Accepted: 26 January 2013
Context. We study an SDSS sample of galaxies within ~50 Mpc of the nominal center of the Ursa Major supercluster.
Aims. Our aim is to study galaxy distribution around groups in the supercluster and the link between the distribution of relaxed and nonrelaxed galaxy systems with respect to the supercluster environment.
Methods. Using the FoF algorithm, 40 galaxy groups were identified in this region. Velocity distributions for these groups were studied after applying a 3σclipping routine for outlier removal. We classified the systems according to the results of normality and substructure tests applied to member galaxies. Then, we studied the relative distribution of relaxed and nonrelaxed systems across the supercluster.
Results. We find that 68% of galaxy groups are Gaussian and that all the nonGaussian systems have substructures and probably correspond to multimodal systems in redshift space. We also find that the Gaussian systems inhabit the denser regions of the supercluster, with higher densities of both red and blue galaxies within 2.5 h^{1} Mpc, and have smaller groupgroup pairwise separations.
Conclusions. Our results suggest a spatial segregation of dynamical states, where relaxed systems may have formed and evolved earlier and faster around highdensity peaks, while nonrelaxed systems may be growing slower on the peripheries of lower density peaks. In this picture, galaxy clustering seems to be prompting a continuous internal evolution in the supercluster, with several groups collapsing into the more evolved and contracted regions.
Key words: galaxies: clusters: general / galaxies: groups: general
© ESO, 2013
1. Introduction
Groups of galaxies contain most of galaxies in the Universe and are the link between individual galaxies and largescale structures (e.g., Huchra & Geller 1982; Geller & Huchra 1983; Nolthenius & White 1987; Ramella et al. 1989). The dissipationless evolution of these systems is dominated by gravity. Interactions over a relaxation time tend to distribute the velocities of the galaxy members into a Gaussian distribution (e.g. Bird & Beers 1993). Although the theoretical lineofsight (los) velocity distribution expected in galaxy systems is not exactly Gaussian (e.g. Merritt 1987; Kazantzidis et al. 2004), phenomenological evidence has been suggesting for a long time that normality can be assumed for systems in dynamical equilibrium (e.g. Yahil & Vidal 1977).
Several works have shown that the analysis of velocity distributions can separate galaxy groups into relaxed or nonrelaxed systems (e.g. Hou et al. 2009; Ribeiro et al. 2010, 2011). Actually, several important properties of groups can be studied from this perspective. For instance, Hou et al. (2009, 2012) find rising velocity dispersion profiles for nonGaussian (NG) groups. They have also shown that the majority of NG groups have substructures. Ribeiro et al. (2010) found that galaxies are redder and brighter in Gaussian (G) groups out to 4 R_{200}. Martinez & Zandivarez (2012) have studied the luminosity functions of G and NG groups, concluding that G groups have a brighter characteristic magnitude. Einasto et al. (2012a,b) studied velocity distributions to find signatures of multimodality in SDSS clusters.
The NG nature of the velocity distribution of NG groups must be caused by recent groupgroup mergers. At the same time, NG groups could be contaminated by chance projections. It is therefore important to establish the real nature of NG groups and how they are linked to the environment. We study an Sloan Digital Sky Survey (SDSS) sample of galaxies within ~50 Mpc of the nominal center of the Ursa Major supercluster (UMaS) aiming to establish a relation between largescale environment and the evolution of galaxy groups. After identifying groups in the supercluster using the friendsoffriends (FoF) algorithm, we split up the sample into relaxed and nonrelaxed systems according to the velocity distribution of member galaxies. Then, we investigate the distribution of galaxy groups in the supercluster with respect to their dynamical state. The paper is organized as follows. Section 2 presents the data and methods we have used, Sect. 3 presents our analysis, and Sect. 4 briefly summarizes the main conclusions of the work.
2. Data
In this work, we have used galaxies from the SDSSDR6 (AdelmanMcCarthy et al. 2008) within ~50 Mpc around the coordinates of the cluster Abell 1383 (=177.46°, Dec = 54.62°, and z = 0.059) taken here as the nominal (and arbitrary) “center” of the UMaS. The redshifts are in the range 0.045 < z < 0.075. We have found 18 204 galaxies in this region with apparent magnitude limit m_{r} = 17.77. Clusters were identified using the FoF algorithm (Huchra & Geller 1982), following the choice of linking lengths by Berlind et al. (2006) (b_{⊥} = 0.14 and b_{∥} = 0.75)^{1}, which were optimized for SDSS. The FoF algorithm with these linking lengths finds galaxy groups with N_{gal} ≥ 10 that have an unbiased multiplicity function (see Berlind et al. 2006). These settings led us to find 40 galaxy groups in the UMaS, with numbers of galaxies in the range 15 ≤ N_{gal} ≤ 94. Our catalog contains all but one of the 23 galaxy systems known in the region^{2} (see Kopylova & Kopylov 2009).
Velocity distributions for the identified groups were studied after applying a 3σ clipping routine for removing outliers. First, a virial analysis^{3} is done following Carlberg et al. (1997), where and M_{200} = 3 R_{200}σ^{2}/G. Then, we applied three well known normality tests to the cleaned velocity distributions – the AndersonDarling (AD), ShapiroWilk (SW) and JarqueBera (JB) tests – to decide if a galaxy system is Gaussian or NG. Our criterium is that at least two tests should reject normality at the 95% confidence level to be considered NG; if not, it will be considered Gaussian. Finally, we applied two substructure tests – the β (West et al. 1988) and Δ tests (Dressler & Schectman 1988) – to decide if a galaxy system has substructures. Our criterion is that a galaxy system is considered to have substructures if at least one of the two tests rejects regularity with 95% confidence level. Results are presented in Table 1, where columns are (1) group identification; (2) and (3) equatorial coordinates in degrees; (4) average redshift of members after removing outliers; (5) velocity dispersion in km s^{1}; (6) R_{200} radius in h^{1} Mpc; (7) logarithm of M_{200} in h^{1}M_{⊙}; (8) total number of galaxies after outlier removal; (9) results for the NG nature of the los velocity distribution: 0 for Gaussian and 1 for NG distributions; (10) results for substructure tests: 0 for regular and 1 for irregular systems; (11) Abell cluster identification.
Properties of 40 galaxy groups identified in the UMaS.
Fig. 1 Four examples of multimodality diagnostics for groups in the UMaS region. Vertical bars depict individual galaxies in each group. The solid lines indicate the best Gaussian mixture for for radial velocities distributions in km s^{1}. The dashed lines show the individual components of the multimodal systems. The dotted lines indicate the average Gaussian curve for 1000 resamplings of groups 03 and 34, unimodal but nonGaussian systems. 
3. Analysis
3.1. Multimodality
We find that 27 out of 40, i.e. 68%, of galaxy groups are Gaussian and that all the NG systems have substructures. Also, we find that four groups have substructures (as indicated by the β test) but have Gaussian velocity distributions according to our criteria. If we include them as nonrelaxed systems, the total fraction of virialized groups would fall to 58%. We take a conservative view to define nonrelaxed groups and assume that substructures are dynamically significant only if they are also associated to nonGaussianities in the velocity distribution of a galaxy system.
Fig. 2 Distribution of galaxies and galaxy systems in the UMaS in supergalactic coordinates. Galaxies are in black; Gaussian groups and nonGaussian are in red and blue, respectively. 
Following Ribeiro et al. (2011), we applied the KolmogorovSmirnov (KS) and the Cramérvon Mises (CvM) twosample tests on the G and NG subsamples. We find that NG velocity groups are signifcantly larger and more massive at the 95% confidence level for both tests, with best fit relations:
with [(p = 0.0065 – KS), (p = 0.0019 – CvM)] for the radius distributions, and [(p = 0.0029 – KS) (p = 0.0017 – CvM)] for the mass distributions. This result is similar to the one found by Ribeiro et al. (2011), a possible consequence of applying the virial analysis to nonvirialized objects. Actually, we can assume the velocity distributions of NG groups as Gaussian mixtures with unknown number of components, and use the Dirichlet process mixture (DPM) model to study the velocity distributions (see Ribeiro et al. 2011 for details). The number of components in the mixture is computed using the R language and environment (R Develoment Core Team) under the dpmixsim library (da Silva 2009). We find the following results:

Unimodal groups: 3, 9, 34

Bimodal groups: 4, 6, 18, 19, 22, 25, 26, 28, 30

Trimodal group: 1.
In Fig. 1 we plot the velocity distributions for typical cases of NG groups. Multimodality seems to be the prominent cause of nonGaussianity in galaxy group velocity distributions, corresponding to 10 out of 13, i.e. 77% of the NG cases in the UMaS. There are three unimodal systems (#03, #09, and #34) all having negative kurtosis excess (k = −1.26, − 0.87, and k = −1.05, respectively) – the mean values for Gaussian and NG groups are k_{G} = −0.34 and k_{NG} = −1.16 A negative kurtosis corresponds to a flattened and usually multimodal distribution. However, if the separation between the modes is too small, the DPM model is unable to recognize distinct components. In Fig. 1, we illustrate two cases of NG unimodal systems, also plotting the Gaussian curve resulting from 1000 resamplings keeping 90% of the data^{4}. For each group we apply the normality tests mentioned in Sect. 2 to the new velocity distribution. The reduced samples classified as normal distributions are used to compute the average Gaussian we plot in Fig. 1.
Fig. 3 Galaxy number density profiles in redshift space. Filled and open circles denote G and NG systems, respectively. Error bars are obtained from 1000 bootstrap resamplings of galaxy groups. Top, middle and bottom plots refer to all, red, and blue galaxies, respectively. 
3.2. Supercluster environment
Now, we try to link the dynamical state of galaxy systems and the supercluster environment. The spatial distribution of galaxy groups in the UMaS region is presented in Fig. 2 in supergalactic Cartesian coordinates^{5}, where we see galaxies and groups. From this figure, it seems that G groups inhabit the denser regions of the supercluster, approximately following the central archshaped structure in the plot. That could be hinting at some environmental influence in the galaxy groups spatial distribution.
We have used two simple indicators to probe the possible relation between galaxy groups and the supercluster environment.
3.2.1. Distribution of galaxies around groups
First, we built the stacked G and NG galaxy number density (in redshift space) up to 2.5 h^{1} Mpc, around the groups’ centers, a distance corresponding to ≃ for our sample. Results are presented in Fig. 3, where we see similar shapes for the profiles around G and NG groups, but with densities systematically higher in the case of G groups. Error bars in this plot are obtained from a bootstrap technique with 1000 resamplings applied to the galaxy groups. This result indicates that G systems present significantly richer surroundings than NG sytems, with ratio of number densities of galaxies around G and NG groups following log [ν_{NG}(r)/ν_{G}(r)] ≃ 0.35 ± 0.10 (ν is the galaxy number density).
We also classified galaxies into blue or red objects following the criterion of Strateva et al. (2001) – a red galaxy should have u − r ≥ 2.22. This enables us to rebuild the galaxy density profiles for these two populations. The results for red galaxies are also presented in Fig. 3 (the middle plot), where we see that G groups have more red galaxies within 2.5 h^{1} Mpc than NG groups ± 0.11, but note that the difference are smaller at larger radii, reflecting (i) the well known correlation between colors and environment – the fraction of red galaxies decreases toward the outer part of clusters (e.g. Balogh et al. 2000; De Propis et al. 2004); and (ii) the quicker decline for G systems. At the same time, for blue galaxies the difference between G and NG systems is smaller, but still statistically significant: , see Fig. 3 (on the bottom plot). These results show that G groups are denser and redder than NG groups up to 2.5 h^{1} Mpc.
3.2.2. Distribution of groups through the UMaS
Also, we studied the relative separation between groups in both G and NG subsets. We computed the distance from each group to the first (k = 1), second (k = 2), and third (k = 3) nearest groups in the spatial distribution. In Fig. 4, we illustrate our results for the supergalactic YZ projection^{6}. The mean pairwise distances of G groups show just a small variation with k: 11 h^{1} Mpc (k = 1), 12 h^{1} Mpc (k = 2), and 14 h^{1} Mpc (k = 3); while we see a wider variation for NG groups: 10 h^{1} Mpc (k = 1), 19 h^{1} Mpc (k = 2), and 26 h^{1} Mpc (k = 3). This suggests that G groups are closer to each other than NG groups.
We estimated the significance of the pairwise separation test by scrambling the identities of the G and NG groups 10 000 times, preserving the overall numbers of 27 Gaussian and 13 nonGaussian. We calculated the number of replicas with more extreme statistics than the input data (i.e. larger differences between G and NG mean pairwise separations). For k = 1, we found 4356 cases, indicating that G and NG systems have the same distributions with respect to the first neighbor. However, for k = 2 and k = 3 we found 279 and 232 cases of more extreme statistics. We can therefore reject the null hypothesis of the same spatial distribution for G and NG systems at the ~2.8% and ~2.3% signficance levels, respectively. The difference between the spatial distributions of G and NG groups is reinforced by the result of the twodimensional KS test (e.g. Press et al. 1992; Fasano & Franceschini 1987), which evaluates the null hypothesis that a pair of data sets are drawn from the same distribution. We found that this is rejected at the 95% confidence level for all 2D data projections.
Combined, these results indicate that G groups probably inhabit denser regions, with higher densities of red galaxies, and are closer to each other (at least for k = 2 and k = 3) than NG groups – a trend that suggests an important environmental effect on structure formation through the UMaS volume.
Fig. 4 Pairwise distances from each group to the first (k = 1), second (k = 2) and third (k = 3) nearest neighbors in the supergalactic YZ projection. The arrows point to the kth nearest neighbor among groups of the same class. The red and blue arrows indicate distances for G and NG systems, respectively. The mean distance between groups appears on the right side of each box. The size of all boxes is 100 × 100 h^{1} Mpc. 
4. Discussion
Superclusters are the largest structures identified in the universe. They correspond to complex systems comprising a few up to several galaxy clusters and groups experiencing longterm gravitational evolution. In this work, we studied a galaxy sample ~50 Mpc around an arbitrary center (the Abell cluster 1383) of the UMaS. We have addressed two questions: galaxy distribution around groups in the supercluster and the link between the distribution of relaxed and nonrelaxed galaxy systems with respect to the largescale environment in which they are embedded. We find that 68% of galaxy groups are Gaussian and that all the NG systems have substructures. This corresponds to a lower fraction of relaxed groups in comparison to our previous studies (Ribeiro et al. 2010, 2011)^{7}. This result cannot be due to the different criteria employed to split up G and NG groups, since we are being more conservative to define nonrelaxed systems in the present work. Alternatively, this may result from the supercluster environment, which could be favoring interactions and mergers between the groups. Actually, finding multimodality as the main cause of nonGaussinity seems consistent with the idea of ongoing processes in nonrelaxed systems (e.g. Ribeiro et al. 2011; Hou et al. 2012; Einasto et al. 2012a).
In this respect, an important issue relates to the possibility of also having NG groups as the result of spurious members in the velocity distributions. Wojtak et al. (2007) show that most statistical methods (including the 3σclipping procedure used in the present work) exclude 60%–70% of unbound galaxies on average. The probability of chance alignments with the remaining interlopers was estimated with the following procedure:

1.
We assume that all groups are 100% reliable in the membership.

2.
We place all galaxies within them at the center, making their velocity dispersions nil.

3.
And we rotate the cube around a random axis.

4.
Then, we disperse the los velocities according to a Gaussian distribution for each group with the original velocity dispersion for the group.

5.
Finally, we apply an FoF code with the same parameters mentioned in Sect. 2, and analyze the Gaussianity of the resulting groups.
After 1000 repetitions of this procedure, we found that ~15% of the mock groups, which are Gaussian by construction, become NG, once additional members are included. This means that about two NG groups in our sample would correspond to Gaussian systems with spurious members in their velocity distributions. This fraction is too low to signicantly modify the results we found in the previous sections, so most of NG groups are probably portraits of multimodal systems.
Even presuming this is so, Gaussian groups are still dominant across the supercluster, which suggests the system is not young enough. Actually, Einasto et al. (2012b) show that superclusters have two main morphological types, spiders and filaments. Clusters in superclusters of spider morphology have higher probabilities to have substructure and higher peculiar velocities of their main galaxies than clusters in superclusters of filament morphology. The UMaS is classified as filamenttype in that work, which suggests it might be dynamically evolved. In the same direction, we found that G groups inhabit denser regions, with higher densities of both red and blue galaxies within 2.5 h^{1} Mpc and are closer to each other than NG groups, suggesting that relaxed systems may have formed and evolved earlier around highdensity peaks, while nonrelaxed systems may be growing more slowly in the peripheries of lower density peaks. This shows the complex nature of superclusters so that further studies are needed to understand the role of the largescale environment in galaxy group distribution and evolution.
But consistent with Hou et al. (2009) for CNOC2 data.
Acknowledgments
We thank the referee Gary Mamon for very useful suggestions. We are thankful to Jason Pinkney for making the substructure codes available. A.L.B.R. acknowledges the support of the CNPq, grants 306870/20100 and 478753/20101. P.A.A.L. acknowledges the support of the FAPERJ, grant E26/110.237/2010, and CNPQ, grant 304692/20115.
References
 AdelmanMcCarthy, J. K., AdelmanMcCarthy, J. K., Agüeros, M. A., et al. 2008, ApJS, 175, 297 [NASA ADS] [CrossRef] [Google Scholar]
 ArayaMelo, P. A., Reisenegger, A., Meza, A., et al. 2009, MNRAS, 399, 97 [Google Scholar]
 Balogh, M., Navarro, J., & Morris, S. L. 2000, ApJ, 540, 113 [NASA ADS] [CrossRef] [Google Scholar]
 Berlind, A. A., Frieman, J., Weinberg, D. H., et al. 2006, ApJS, 167, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Bird, C., & Beers, T. 1993, AJ, 105, 1596 [NASA ADS] [CrossRef] [Google Scholar]
 Carlberg, R. G., Yee, H. K. C., & Ellingson, E. 1997, ApJ, 478, 462 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 da Silva, A. F. 2009, Comput. Methods Programs Biomed., 94, 1 [CrossRef] [Google Scholar]
 De Propis, R., De Propris, R., Colless, M., et al. 2004, MNRAS, 351, 125 [NASA ADS] [CrossRef] [Google Scholar]
 Dressler, A., & Shectman, S. A. 1988, AJ, 95, 985 [NASA ADS] [CrossRef] [Google Scholar]
 Dünner, R., Araya, P. A., Meza, A., & Reisenegger, A. 2006, MNRAS, 366, 803 [NASA ADS] [CrossRef] [Google Scholar]
 Einasto, M., Vennik, J., Nurmi, P., et al. 2012a, A&A, 540, A123 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Einasto, M., Liivamägi, L. J., Tempel, E., et al. 2012b, A&A, 542, A36 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Fasano, G., & Franceschini, A. 1987, MNRAS, 225, 155 [NASA ADS] [CrossRef] [Google Scholar]
 Geller, M. J., & Huchra, J. P. 1983, ApJS, 52, 61 [NASA ADS] [CrossRef] [Google Scholar]
 Hou, A., Parker, L., Harris, W., & Wilman, D. J. 2009, ApJ, 702, 1199 [NASA ADS] [CrossRef] [Google Scholar]
 Hou, A., Parker, L. C., Wilman, D. J., et al. 2012, MNRAS, 421, 3594 [NASA ADS] [CrossRef] [Google Scholar]
 Huchra, J. P., & Geller, M. J. 1982, ApJ, 257, 423 [NASA ADS] [CrossRef] [Google Scholar]
 Kazantzidis, S., Magorrian, J., & Moore, B. 2004, ApJ, 601, 37 [NASA ADS] [CrossRef] [Google Scholar]
 Merritt, D. 1987, ApJ, 313, 121 [NASA ADS] [CrossRef] [Google Scholar]
 Kopylova, F. G., & Kopylov, A. I. 2009, Astrophysical Bulletin, 64, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Martínez, H. J., & Zandivarez, A. 2012, MNRAS, 419, L24 [Google Scholar]
 Nolthenius, R., & White, S. D. M. 1987, MNRAS, 225, 505 [NASA ADS] [CrossRef] [Google Scholar]
 Press, W. H., Teukolsky, S. A., Vetterling, W. T., & Flannery, B. P. 1992, Numerical Recipes in C, 2nd edn. (Cambridge: Cambridge University Press) [Google Scholar]
 Ramella, M., Geller, M., & Huchra, J. P. 1989, ApJ, 344, 57 [NASA ADS] [CrossRef] [Google Scholar]
 Ribeiro, A. L. B., Lopes, P. A. A., & Trevisan, M. 2010, MNRAS, 409, L124 [NASA ADS] [Google Scholar]
 Ribeiro, A. L. B., Lopes, P. A. A., & Trevisan, M. 2011, MNRAS, 413, L81 [NASA ADS] [Google Scholar]
 Strateva, I., Ivezić, Ž., Knapp, G. R., et al. 2001, AJ, 122, 1861 [Google Scholar]
 West, M., Oemler, A., & Dekel, A. 1988, ApJ, 327, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Wojtak, R., Lokas, E. L., Mamon, G. A., et al. 2007, A&A, 466, 437 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Yahil, A., & Vidal, N. V. 1977, ApJ, 214, 347 [NASA ADS] [CrossRef] [Google Scholar]
All Tables
All Figures
Fig. 1 Four examples of multimodality diagnostics for groups in the UMaS region. Vertical bars depict individual galaxies in each group. The solid lines indicate the best Gaussian mixture for for radial velocities distributions in km s^{1}. The dashed lines show the individual components of the multimodal systems. The dotted lines indicate the average Gaussian curve for 1000 resamplings of groups 03 and 34, unimodal but nonGaussian systems. 

In the text 
Fig. 2 Distribution of galaxies and galaxy systems in the UMaS in supergalactic coordinates. Galaxies are in black; Gaussian groups and nonGaussian are in red and blue, respectively. 

In the text 
Fig. 3 Galaxy number density profiles in redshift space. Filled and open circles denote G and NG systems, respectively. Error bars are obtained from 1000 bootstrap resamplings of galaxy groups. Top, middle and bottom plots refer to all, red, and blue galaxies, respectively. 

In the text 
Fig. 4 Pairwise distances from each group to the first (k = 1), second (k = 2) and third (k = 3) nearest neighbors in the supergalactic YZ projection. The arrows point to the kth nearest neighbor among groups of the same class. The red and blue arrows indicate distances for G and NG systems, respectively. The mean distance between groups appears on the right side of each box. The size of all boxes is 100 × 100 h^{1} Mpc. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.