Issue 
A&A
Volume 668, December 2022



Article Number  A37  
Number of page(s)  14  
Section  Cosmology (including clusters of galaxies)  
DOI  https://doi.org/10.1051/00046361/202244239  
Published online  01 December 2022 
Quasispherical superclusters
^{1}
Tuorla Observatory, Department of Physics and Astronomy, University of Turku, Vesilinnantie 5, 20014 Turku, Finland
email: pekheina@utu.fi
^{2}
Institut d’Astrophysique Spatiale, CNRS (UMR8617) Université ParisSud, Bâtiment 121, 91405 Orsay, France
^{3}
Biodiversity Unit, University of Turku, 20014 Turku, Finland
^{4}
Tartu Observatory, Tartu University, 61602 Tõravere, Estonia
^{5}
Estonian Academy of Sciences, Kohtu 8, 10130 Tallinn, Estonia
Received:
10
June
2022
Accepted:
29
September
2022
Context. Superclusters are systems with varied properties and varied fractional overdensities. Their dynamical state evolves under the influence of two components: dark energy and gravitational force. The dominant component at any spatial location and cosmic epoch is determined by the total mass and the local overdensity of the system. However, generally the dynamical state of superclusters is poorly known.
Aims. We study properties of superclusters and select a sample of quasispherical superclusters, the dynamics of which can be studied using the Λ significance diagram.
Methods. We extracted our supercluster sample with an adaptive local threshold density method from the Sloan Digital Sky Survey Data Release 7 (SDSS DR7) data and estimated their masses using the dynamical masses for member galaxies and groups. We used topological analysis based on Minkowski functionals and the positions of galaxies and galaxy groups in superclusters. Finally, we highlight the dynamical state of a few exceptional types of superclusters found in this study using the Λ significance diagram.
Results. Our final sample contains 65 superclusters in the distance range of 130−450 Mpc. Supercluster masses range between 1.1 × 10^{15} M_{⊙} and 1.4 × 10^{16} M_{⊙} and sizes between 25 Mpc and 87 Mpc. We find that pancaketype superclusters form the lowluminosity, small, poor and lowmass end of superclusters. We find four superclusters of unusual types, exhibiting exceptionally spherical shapes. These socalled quasispherical systems contain a highdensity core surrounded by a relatively spherical density and galaxy distribution. The masstolight ratio of these quasisphericals is higher than those of the other superclusters, suggesting a relatively high dark matter content. Using the Λ significance diagram for oblate and prolate spheroids, we find that three quasispherical superclusters are gravitationally bound at the present epoch.
Conclusions. Quasispherical superclusters are among the largest gravitationally bound systems found to date, and form a special class of giant systems that, dynamically, are in between large gravitationally unbound superclusters and clusters of galaxies in an equilibrium configuration.
Key words: largescale structure of Universe / dark energy / dark matter
© P. Heinämäki et al. 2022
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
This article is published in open access under the SubscribetoOpen model. Subscribe to A&A to support open access publication.
1. Introduction
According to the current cosmological paradigm, the great majority of the contents of the Universe is governed by the dark sector. In this picture, formation and evolution of the cosmic structure is driven by the competition between the gravitational attraction of the dark matter (DM) and the expansion of the space itself. In the standard ΛCDM cosmology, the latter is assumed to be currently accelerated due to the constant energy density of the vacuum, which appears as a cosmological constant, Λ, in Einstein’s field equations. During the matterdominated era of cosmic history, the gravitational attraction of the initial primordial density fluctuations of the DM leads to the formation of increasingly large concentrations of matter.
While galaxies and clusters are highly concentrated virialized systems with substantial fractional overdensities δ_{m} ∼ 10^{6} − 10^{2}, superclusters as a whole are perceived as large and relatively loose, unrelaxed agglomerations of galaxies and clusters, with overdensities of just a few times the mean density of the Universe. Superclusters are generally surrounded by voids, that is, underdense regions with δ_{m} < 1. Due to their small fractional overdensity, superclusters as a whole remained in the linear regime until recent times and finally around z ∼ 0.5 (Frieman et al. 2008), when the dark energy (DE) density started to dominate over the matter density, and the size of the evolving concentrations (the growth of linear perturbations) was halted at supercluster scales. On the other hand, because of the nonuniform distribution of gravitating matter, the DE (cosmic ‘antigravity’) can be stronger than gravity also locally on scales of ∼1 − 10 Mpc (Chernin 2001; Byrd et al. 2012). Similarly, we know that superclusters may have collapsing cores and one or several merging clusters in their highdensity cores (Einasto et al. 2021 and references therein). Therefore, the presence of DE along with gravitating matter influences the formation of the largescale structure on all scales from groups of galaxies to superclusters (Teerikorpi et al. 2015).
Superclusters form a heterogeneous sample of structures where the observed properties vary considerably in size, richness, luminosity, and morphology (Einasto et al. 1980, 2011a; Davis et al. 1982; Liivamägi et al. 2012). Due to the small fractional overdensities of superclusters, their extraction from the observed galaxy surveys requires careful analysis of the galaxy density field. Several methods have been introduced for reconstruction of the suitable density field from large galaxy surveys: gridbased density with appropriate kernel smoothing (e.g. Basilakos et al. 2001), luminosityweighted density field (e.g. Einasto et al. 2003), the Wiener filtering technique (e.g. Erdogdu et al. 2004), and the Voronoi tesselationbased methods (e.g. Nadathur & Crittenden 2016). A dynamical definition of superclusters was introduced by Tully et al. (2014). Using a combination of the density and peculiar velocity fields of galaxies (Cosmicflows2 catalogue) in our local Universe, these latter authors defined the Laniakea supercluster as the volume (with a characteristic diameter of about 160 Mpc) within which the galaxies have an inward peculiar velocity (v_{pec}) component. Outside the outer surface of the Laniakea supercluster, the matter is expanding along with the Hubble flow and velocity (u = Hr) or it has a peculiar velocity in the direction of some other mass concentrations. Within this concept, the regions of coherent inward motions are called ‘basins of attraction’ while the regions of coherent outflow are called ‘basins of repulsion’ (e.g. Dupuy et al. 2019). Detailed peculiar velocity information for galaxies is available only for our local Universe.
To probe the dynamical state of superclusters in much larger cosmological volumes, we must use some other methods. For instance, if one considers only relatively spherical systems, one may use the luminosity density field data together with analytical approaches, such as the spherical collapse assumption (e.g. Einasto et al. 2020) and the concept of the Λ significance diagram (Teerikorpi et al. 2015) to characterise the dynamical state of a supercluster or its parts (e.g. Gramann et al. 2015; Einasto et al. 2015, 2021; Chon et al. 2015; Teerikorpi et al. 2015). The Λ significance diagram presents a useful way to distinguish dynamically different (expanding or contracting) regions from each other on different spatial and mass scales around superclusters. In the dynamical picture of structures, at one end we have large gravitationally unbound superclusters defined as systems with inward peculiar velocity components relative to the Hubble flow, and at the other end there are clusters of galaxies in an equilibrium configuration. It is expected that between these extremes there is a continuum of structures with different dynamical states.
In this paper, our aim is to use topological tools and the Λ significance diagram in order to study the properties and dynamical state of the superclusters from the Sloan Digital Sky Survey Data Release 7 (SDSS DR7) sample. To achieve the latter objective, we are especially interested in the most spherical superclusters of the sample. We use the Hubble constant H_{0} = 67 km^{−1} Mpc^{−1}, and for the global density of the dark energy, ρ_{Λ} = 6 × 10^{−30} g cm^{−3} (Planck Collaboration XVI 2014). The dimensionless Hubble parameter h = 0.67 is used in all calculations made in this paper.
2. The sample of SDSS superclusters
2.1. Data
We use the fluxlimited supercluster sample drawn from the SDSS DR7 main survey by Liivamägi et al. (2012). This was constructed using the 3D luminosity density field of galaxies smoothed with a 8 h^{−1} Mpc kernel. Before compiling the supercluster catalogue, galaxy data were reduced to minimise various selection effects. First, to calculate the comoving distances of the galaxies, their redshifts were corrected into the CMB frame. Secondly, the SDSS sample is very narrow at small distances, and becomes very diluted at the far end of the sample. Therefore, we selected the lower and upper distance limits of the sample to be z = 0.009 and z = 0.2. Finally, the absolute magnitudes of galaxies were calculated taking into account galactic extinction and k+ecorrections.
As a next step, groups and clusters of galaxies were determined using a modified friendsoffriends (FoF) algorithm described in detail in Tago et al. (2010). In the redshift space, both radial and transverse linking lengths were determined separately. Due to the use of a fluxlimited sample, the linking length was increased slightly as a function of distance (for more details, see Tago et al. 2010). Distances based on spectroscopic redshifts are affected by peculiar velocities of galaxies in groups and clusters. This effect of the radial redshift distortions (the fingerofgod effect) on the distances was corrected using the standard deviation σ_{r} of the projected galaxy distances from the group centre and the standard deviation σ_{v} of the radial velocities of galaxies (both in physical coordinates at the group location). The corrected radial distance of each galaxy of the group was calculated using the ratio of the mentioned standard deviations multiplied by the distance between the group centre and the galaxy. This procedure efficiently suppresses the apparent elongation of groups in the redshift space (for more details, see Liivamägi et al. 2012; Tempel et al. 2012). The luminosity density field calculations are carried out in the Cartesian grid. For that purpose, the corrected galaxy distances were also transferred from the SDSS angular coordinates to the Cartesian coordinates.
Luminosities of the fluxlimited sample of galaxies were corrected for the distance selection effect using the procedure by Tempel et al. (2011). Here the luminosity of each galaxy was corrected by multiplying it by a distancedependent weighting factor. This factor is defined as the ratio of the expected total luminosity to the luminosity within the magnitude limits at the distance of the galaxy. The final luminosity density field was calculated using a B_{3} spline kernel of the scale of 8 h^{−1} Mpc and the corrected luminosities (see more details in Liivamägi et al. 2012).
The final sample of superclusters (a total of 1313) was determined from the luminosity density field using an adaptive density threshold method which gives each supercluster its own local threshold limit. It is important to note that unlike in the ordinary method based on fixed density levels, in the adaptive method every supercluster has its own limiting density level, similarly to other astronomical objects (see Liivamägi et al. 2012 for details). As a result, many superclusters defined using the fixed density level method are split into several superclusters in the adaptive case. Thus, the number of superclusters is higher (∼10%) in the adaptive method catalogue.
In addition to total luminosities (the total weighted luminosity of the galaxies and galaxy groups) of the superclusters, this catalogue gives for example the maximum sizes and volumes of superclusters. The volume of a supercluster is calculated from the luminosity density field as the summed volumes of the connected 1 (h^{−1} Mpc)^{3} grid cells (the resolution of the density field calculation) within this supercluster. This method does not assume any specific shape for the structure, making morphological analysis possible. The catalogue gives morphology information for the superclusters based on Minkowski functionals and shapefinders derived from them.
2.2. Masses of superclusters and selection biases
The original catalogue by Liivamägi et al. (2012) did not contain mass estimates for superclusters. In this study, we approximate their masses using dynamical masses. The dynamical mass is obtained by summing together the dynamical masses of all the galaxy groups belonging to a supercluster. For that purpose, we first identified galaxy groups of Tago et al. (2010) in our superclusters. We estimate the dynamical masses of the groups within the radius where the average density is 200 times the critical density (ρ_{c}) using Eq. (8) in Tempel et al. (2014):
where σ_{v} is the 1D velocity dispersion and the gravitational radius R_{g} = 4.582σ_{sky}. According to Tempel et al. (2014), the velocity dispersion is estimated using the lineofsight velocities of all detected galaxies of a galaxy group (or cluster) following the standard formula , where v_{mean} and z_{m} are the mean group velocity and redshift, respectively, v_{i} is the velocity of an individual group member, and n is the number of galaxies with observed velocities within the group. Assuming dynamical symmetry, the real (3D) velocity dispersion in groups would be . In order to estimate R_{g} of the galaxy groups, the rms deviation σ_{sky} of the projected distances of group galaxies are calculated using the formula , where r_{i} is the projected distance in the sky from the group centre (in comoving coordinates). The equation of the gravitational radius R_{g} comes from mass estimations based on the Hernquist profile, where σ_{sky} is the rms deviation of the projected distances of group galaxies in the sky from the group centre. Within the mass range 10^{9} h^{−1} M_{⊙} − 10^{15} h^{−1} M_{⊙}, the Hernquist masses are 1.55−1.75 times higher than the NFW masses. We use the average value 1.67 to convert the masses to approximate the NFW masses, which are commonly used estimates of the dark matter halo masses. For the group dynamical masses ≥10^{14} M_{⊙}, the errors are estimated using the results based on the cluster mass reconstruction analysis by Old et al. (2014, 2015), where over 20 different clustermassestimation techniques were analysed (including the method used here) aiming to study statistical uncertainties in clusters masses. For the method used here, they found an intrinsic scatter of 0.3 dex (∼a factor of 2) for masses ≥10^{14} M_{⊙}. For the less massive < 10^{14} M_{⊙} groups, the errors are larger and we used onesigma confidence level errors for all groups with masses < 10^{14} M_{⊙} and richer than five members. Old et al. (2014, 2015) found in their analysis that, in general, the method used here has a tendency to underestimate dynamical masses.
For small groups (N < 5), the observables σ_{v} and R_{vir} are not well defined, which leads to large scatter in the small group masses. Following Gramann et al. (2015) and Einasto et al. (2015), we assigned to each group with less than five members the median mass of such groups, which in our case is 3.3 × 10^{12} ± 5.6 × 10^{10} M_{⊙}.
Each supercluster also contains single galaxies. We assume that, due to the survey magnitude limit, every single galaxy is a member of an unresolved small group and therefore we also assigned to each such galaxy the median mass of the groups that have less than five members. In addition, as in Gramann et al. (2015), Einasto et al. (2015, 2021), we assume that the intracluster gas increases the total mass of the supercluster by 10%. Although the total luminosities of the groups are corrected via a weighting procedure (the ratio of the expected total luminosity to the expected luminosity in the visibility window) because of the apparent magnitude limit of the sample, the group richness decreases inevitably as a function of distance. This means that the virial mass estimation is not equally reliable in the whole range of distances in the SDSS survey. Tago et al. (2010) showed that the estimated group richness remains reliable up to a distance of about 450 Mpc. After applying this limit, the number of superclusters in the original sample was reduced to 163 superclusters.
Einasto et al. (2011a) estimated that the minimum richness (number of galaxies in the supercluster) to resolve the morphological properties of superclusters is near the value 300. We also used MonteCarlo simulations to test the reliability of the correlation of the Minkowski shape parameters and the galaxy number density in superclusters and found an agreement with the requirement of the minimum richness (∼300). Based on this, and in order to maximise the number of the superclusters in our analysis, we required that a supercluster richness be more than 280 galaxies.
Finally, we removed all the superclusters that have a contact with the edge of the survey mask. The original sample contains 1206 superclusters, and after applying the selection criteria the final sample for our analysis contains 65 superclusters. Figure 1 shows the mass distribution of the final supercluster sample. The supercluster masses range between 1.1 × 10^{15} M_{⊙} and 1.4 × 10^{16} M_{⊙}. The median mass is 2.6 × 10^{15}. The lower panel in Fig. 1 shows the supercluster size (diameter) distribution as a function of the comoving distance. Size is defined as a maximum distance between the galaxies inside a supercluster. Sizes range between 25 Mpc and 87 Mpc, while the median size is about 49 Mpc. The figure shows a large amount of scatter for supercluster size as a function of distance. Due to the survey geometry, the number of superclusters increases with distance as the survey volume increases. Our result agrees with those of Liivamägi et al. (2012), who showed that the average supercluster size remains approximately the same across the survey volume. The gap in the distribution around 150−250 Mpc is the known void region between the nearby superclusters (Einasto et al. 2011a).
Fig. 1. Mass distribution and size–distance scatter plot of the superclusters. The upper figure shows the mass histogram of the final sample of the 65 superclusters and the lower figure shows supercluster diameter as a function of distance. 
The masses of the large wellknown superclusters have been studied by several authors. The mass of the Shapley supercluster was estimated by Sheth & Diaferio (2011) and Ragone et al. (2006) for example, who found values of 1 − 2 × 10^{16} M_{⊙}. Einasto et al. (2021, 2022) found that the total mass of the Corona Borealis supercluster and the superclusters in the Sloan Great Wall and in the BOSS Great Wall lie in the range of 0.47 − 1 × 10^{16} M_{⊙}. Recently, Böhringer & Chon (2021) used Xrayluminous galaxy groups and clusters to estimate the masses of eight superclusters. The four most prominent of them (the PerseusPisces, the Centaurus, the Coma, and the Hercules supercluster) have masses of between 5 × 10^{15} M_{⊙} and 2.2 × 10^{16} M_{⊙} with sizes from about 40 to more than 100 Mpc. For the four smaller ones, the estimated masses were in the range of 2 × 10^{15} M_{⊙} − 7 × 10^{15} M_{⊙} with sizes of 20 to 46 Mpc. The Laniakea supercluster (Tully et al. 2014) has an estimated mass of ∼10^{17} M_{⊙} and a diameter of about 160 Mpc. We note that in that study, mass and size estimates cover the region of the ‘basins of attraction’, while superclusters in general are extracted using the criterion of fractional overdensity.
3. Morphology
While clusters and groups are, in general, relatively spherical systems, the spherical shape obviously is not the usual case for superclusters. Instead, superclusters most commonly show a filamentary (or elongated prolate) morphology (Jaaniste et al. 1998; Basilakos 2003; Einasto et al. 2011a). All adaptively selected superclusters used in our analysis are also elongated at some level.
A useful tool to study supercluster morphology is the shapefinder based on a set of Minkowski functionals (more details are provided in Sahni et al. 1998; Saar et al. 2007; Einasto et al. 2007; CostaDuarte et al. 2011). Briefly, the final shapefinder is a vector K = (K_{1}, K_{2}) whose components K1 and K2 are based on the shapefinders H_{i}, (i = 1, 2, 3) as and . The shapefinders H_{i} are determined by geometrical quantities (V = volume, S = surface area, C = integrated mean curvature) as H_{1} = V/S, H_{2} = S/C, and H_{3} = C. One can use combinations of the shapefinder parameters K_{1} and K_{2} to measure the shape (planarity, filamentary, sphericity) of the 3D surface of the structure.
The ratio of the vector components K_{1}/K_{2} > 1 corresponds to a ‘pancake’type structure (for an ideal pancake K ≈ (1, 0)) while 0 < K_{1}/K_{2} < 1 corresponds to a filamenttype structure (for an ideal filament K ≈ (0, 1)). When both K_{1} and K_{2} are close to zero, the shape is close to a sphere. A perfect sphere has the parameters K_{1} = 0 and K_{2} = 0. If we consider a triaxial ellipsoid with axes a, b, and c, a sphere K ≈ (0, 0) is described as (a, b = a, c = a) (Sahni et al. 1998; Shandarin et al. 2004). At the other end of the diagonal K1 = K2, the combination K_{1} = 1 and K_{2} = 1 corresponds to the ribbontype morphology. Moving along the diagonal towards K_{1} = 0 and K_{2} = 0, the ribbons are filled and grow to spheres. Figure 2 shows schematically the correspondence between the vector components K_{1} and K_{2} and different general shapes.
Fig. 2. Correspondence between the vectors K_{1} and K_{2} and different shapes. 
Though there is no unambiguous correspondence between the values of the shapefinder vectors and certain shapes, the vectors constitute a continuum of values which indicate different shapes. The fourth Minkowski functional V_{3} describes the sum of the clumps, voids, and tunnels in the system (e.g. Saar et al. 2007). Thus a high value of V_{3} indicates a complicated clumpy morphology, while V_{3} = 1 means a smooth structure with a single centre.
4. Properties of the superclusters
Figure 3 shows the luminosity, mass, number of groups, volume, diameter, and masstolight values of the final sample of 65 superclusters smoothed onto the grid on the K1 − K2 plane. The distribution of the supercluster properties is illustrated with coloured contours. White crosses show the K1 − K2 distribution of the final sample of our superclusters. Above the black inclined line, the superclusters are of filamentary type and below the line they are of pancake type (Fig. 2).
Fig. 3. Coloured contours show the distribution of the values of supercluster properties in the K1 − K2 plane. Quasispherical superclusters K_{1}, K_{2} < (0.05, 0.05) are shown with red points, white points are all the superclusters. Luminosities are in units of 10^{10} L_{⊙}, masses in 10^{15} M_{⊙}, volumes in Mpc^{3}, and diameter in megaparsecs (Mpc). The dark line indicates K_{1} = K_{2} and divides the morphologies between pancakes and filaments. 
A majority of the superclusters are of filamentary type (84%). Among the superclusters, those of filamentary type are more luminous, richer, and larger, in agreement with the results of CostaDuarte et al. (2011) and Einasto et al. (2011b) who used fixed supercluster threshold limits. Our analysis shows that the dynamical mass follows the same trend. However, for the masstolight ratio (M/L), the division between the filamentary and pancake structures is not as clear.
The red points represent superclusters which have both K1 and K2 values of less than 0.05. The choice of the limit 0.05 is based on our finding that only among superclusters with K1 and K2 of less than this value can one find systems with a single dominant centre, which is natural for spherical or pancaketype systems. For such systems, the fourth Minkowski functional V_{3} = 1. We also required that these systems do not split further at higher density thresholds. To distinguish these superclusters from the other ones, we term them quasispherical (QS) superclusters (also referred to hereafter as QSs). Figure 3 shows that these QS (a total of four large red dots) superclusters are smaller, poorer, less luminous, and less massive than the filamentarytype superclusters, on average.
In Fig. 3 from left to right the red dots represent: QS 524, QS 1832, QS 849, and QS 550. For QS 524 K_{1}/K_{2} < 1, which indicates mild filamentarity. In Table 1 we present data for these systems. Their diameters are between 24 and 33 Mpc and the masses range from 1.2 × 10^{15} to 2.7 × 10^{15} M_{⊙}.
Quasispherical superclusters: Coordinates and sizes.
The result of the positionbased crossmatching between QSs and known Abell clusters is shown in Table 1. Figure 4 shows the sky distribution of the galaxies in the region of the QS superclusters viewed from three different orthogonal directions in Cartesian coordinates. Before the transformation of the angular coordinates to the Cartesian one, redshift distortions in the groups and clusters were corrected as described in Sect. 2.
Fig. 4. Distribution of galaxies in the regions of superclusters QS 550, QS 849, QS 524, and QS 1832 viewed from three different orientations (x − y, x − z, y − z) in Cartesian coordinates. Distribution of galaxies in the plane of the sky are shown in Appendix B. 
All QS superclusters contain a known Abell cluster in their central region. The crossmatching with other supercluster catalogues based on a fixed density level showed that QSs can also be identified as a part of the previously known superclusters. QS 550 is in the Hercules supercluster or in the SCl 154 based on the catalogue of superclusters of Abell clusters by Einasto et al. (2001). QS 849 is in the zone of the supercluster SCl 541, extracted by a fixed density threshold limit by Liivamägi et al. (2012). QS 524 is located and centred on the richest cluster A1436 in the Ursa Major supercluster (Einasto et al. 2001) and lastly, QS 1832 is in the Bootes supercluster.
Clearly, the adaptive threshold density method extracts more concentrated superclusters, which represent highdensity regions (or cores) of the superclusters defined by a fixed density level. The QSs have their masses in the range of M ≈ 1 − 3 × 10^{15} M_{⊙}. This is the same mass range as Einasto et al. (2016) found for the components (highdensity cores) in the four superclusters of the Sloan Great Wall (SGW). These latter authors concluded that the core regions of some of these elongated components with radii smaller than ≈12 Mpc may already be collapsing.
An interesting point is that, contrary to the general trend of the properties, the M/L ratios of the QS superclusters are not the smallest ones among the superclusters. Their mean M/L ratio is about 364 M_{⊙}/L_{⊙}, while for the rest of the superclusters of our sample (white crosses in Fig. 3) ⟨M/L⟩ = 319 M_{⊙}/L_{⊙}. The M/L ratios of the QSs are given in Table 2. We see that QS 849 has the lowest M/L ratio. In Sect. 5, we show that QS 849 may also be in a different dynamical phase from the other three QSs. If we leave out QS 849 and consider only QS 550, QS 524, and QS 1832, their mean ⟨M/L⟩ is equal to 406 M_{⊙}/L_{⊙}. This relatively high ratio may indicate an exceptional amount of dark matter in these systems.
Quasispherical superclusters: Properties.
In Table 2, we also present the number of galaxy groups and individual galaxies in each QS. The third column shows the estimated total mass of the QS. The mass errors are calculated as explained in Sect. 2.1. We note that the most massive groups (M_{tot} ≥ 10^{14} M_{⊙}) dominate the total mass budget of the QSs, except in QS 849. The fraction of groups more massive than 10^{14} M_{⊙} is 77% in QS 550, 72% in QS 524, and 55% in QS 1832. In QS 849, the majority (63%) of its total mass is in groups with masses of less than 10^{14} M_{⊙}.
5. The Λ significance diagram
Due to the large variation in supercluster morphology, a complex ensemble of substructures, and our almost total lack of information on the velocity fields outside of the local Universe, the study of the dynamical state of superclusters as a whole is difficult. An exception to this, are superclusters that are relatively spherical and do not contain several massive substructures, similar to the four quasispherical superclusters found above. In those cases, an analytical approach like the spherical collapse model (e.g. Gramann et al. 2015; Einasto et al. 2020) or the Λ significance diagram (Teerikorpi et al. 2015) can be applied. Below, we use the Λ significance diagram to estimate their dynamical state.
The dynamical state of the gravitating system in the expanding Universe can be analysed studying the motion of a test particle in the gravity field of the system treated as a pointlike mass (or a spherical mass distribution) in the background of the antigravitating force due to DE. The main parameter in this approach is the zeroacceleration radius R_{ZG} (e.g. Chernin 2001; Chernin et al. 2006):
which refers to the radius (in real space in 3D) where the radial peculiar component of the test particle velocity is equal to the Hubble expansion velocity in the gravitation field of mass M. A useful parameter, which measures the influence of DE, is the energy density ratio ⟨ρ_{M}⟩/ρ_{Λ}, as calculated for the system under inspection. For that purpose, Teerikorpi et al. (2015) introduced the Λ significance diagram (a log(⟨ρ_{M}⟩/ρ_{Λ}) versus log R graph) where R is the radius of a system, ⟨ρ_{M}⟩ is its average mass density, and ρ_{Λ} is the DE density equal to the global value ρ_{Λ} ≈ 6 × 10^{−30} g cm^{−3}.
The location of a galaxy system (or an inner part of it) in the Λ significance diagram indicates whether its overall dynamics is dominated by gravity or by the outward ‘antigravity’ expulsion of DE. The diagram also displays a few relevant scales, which appear for each fixed mass of a spherical or a slightly flattened system in the ΛCDM cosmology. In Teerikorpi et al. (2015), a few example galaxy systems were also shown in a similar diagram. In the present paper, we also use this graph as a tool to characterise the significance of the DM gravity–DE expulsion interplay for the dynamics of quasispherical superclusters.
In logarithms, the ratio ⟨ρ_{M}⟩/ρ_{Λ} becomes
where M is the mass within the radius R of a system. In the log(⟨ρ_{M}⟩/ρ_{Λ}) versus log R diagram, this relation forms a family of inclined straight lines for different values of the mass M. The intersections of the lines with horizontal lines of constant log(⟨ρ_{M}⟩/ρ_{Λ}) give the following scales or radii for the mass M:

The zerovelocity radius R_{ZV} (at this distance, a test particle systemcentric radial velocity u = Hr − v_{pec} is v_{pec} = Hr and u = 0). This corresponds to the turnaround radius in the spherical collapse model.

The zerogravity radius R_{ZG} (gravity force equal to Einstein’s antigravity force). Thus, acceleration is zero, du/dt = 0, which indicates the maximum radius of a gravitationally bound system at the present epoch.

The EinsteinStraus radius R_{ES} (the radial velocity reaches the Hubble velocity, u = Hr). This last, longest distance corresponds to the spherical volume where the mass M produces an average density that is equal to the cosmic global density.
In Fig. 5, these three different radii are shown for the systems of mass 10^{14} M_{⊙}. For spherical and homogenous systems, the horizontal lines ZV, ZG, and ES correspond to the constant ratios of ⟨ρ_{M}⟩/ρ_{Λ} = 6, 2, and 3/7.
Fig. 5. Λ significance diagram for SDSS clusters of galaxies (purple) with QS superclusters marked in green. The Laniakea supercluster, the central regions of the supercluster A2142 and the Coma cluster with three different mass and size estimates are shown. Inclined lines show the effect of different example mass values. Green vertical lines show the values of the R_{ZV} (zero velocity radius), R_{ZG} (zero gravity radius), and the R_{ES} (the EinsteinStraus radius for the mass 10^{14} M_{⊙}) and the horizontal lines show logarithmic values for the constant ratios of the ⟨ρ_{M}⟩/ρ_{Λ} = 6 (ZV), 2 (ZG), and 3/7 (ES). 
Figure 5 shows different types of structures in the Λ significance diagram. Clusters of galaxies (from the SDSS survey, closer than 450 Mpc) are shown with purple circles. These are well inside the R_{ZV}, the turnaround radius in the spherical collapse model. Gramann et al. (2015) and Einasto et al. (2015) showed in their dynamical analysis that the supercluster A2142 has a relatively spherical core region (inside the radius of about 8 Mpc) and that it is in a collapsing state. We show the core region of the supercluster A2142 as a blue triangle. It is, as expected, close to but above the ⟨ρ_{M}⟩/ρ_{Λ} = 6 (ZV) line and thus inside its zerovelocity radius R_{ZV}. According to Tully et al. (2014), after removal of the mean cosmic expansion and longrange flows, the Laniakea supercluster has a coherent, approximately spherical peculiar inward velocity flow field of matter (defining the basin of attractor) within a radius of about 80 Mpc. Outside this region, matter moves along the Hubble flow or is once again gravitationally attracted by some mass concentrations. Consistently, in Fig. 5, the Laniakea system is located near the ⟨ρ_{M}⟩/ρ_{Λ} = 3/7 (ES) line. The Coma cluster is known for its spherical symmetry. Following Teerikorpi et al. (2015) and Chernin et al. (2013), three different mass estimates at different scales around Coma are plotted in Fig. 5. The Coma cluster seems to be inside the turnaround radius up to nearly 14 Mpc.
The quasispherical superclusters found in this paper are marked with green dots and ID numbers. Only QS 524 is located inside the R_{ZG} radius and perhaps also QS 550, if the error bars are considered. The error bars for the log(⟨ρ_{M}⟩/ρ_{Λ}) values are obtained using the estimated errors for masses in Sect. 2.1. The FingerofGod effect is during the group reconstruction procedure. We take into account the possible Kaiser effect including a 10% error for the QS radii. We also note that the adopted radii are defined as half of the diameter (defined as the maximum distance between galaxies in a QS). Thus, any deviation from the spherical shape of the QS will decrease the actual radii and therefore the value of log(⟨ρ_{M}⟩/ρ_{Λ}). Table 2 shows the estimated log(⟨ρ_{M}⟩/ρ_{Λ}) for all QSs.
In view of the observations that suggest that the more local Hubble constant H_{0} is somewhat larger than the Planck result (the ‘Hubble tension’), we note that the parameter ⟨ρ_{M}⟩/ρ_{Λ} (Eq. (3)) has the useful property that it does not directly depend on H_{0}. This may be seen as follows. The dynamical (virial) mass for each subsystem of a supercluster is proportional to a size parameter, and therefore to , and then the density of the whole system is . On the other hand, the critical density is . Hence, in the ratio ⟨ρ_{M}⟩/ρ_{Λ} the strong quadratic dependence on cancels out and the only cosmological parameter that remains is Ω_{Λ} in the nominator of the ratio. We assume that Ω_{Λ} does not differ very much from its standard value. If it did, for example by 10%, then the shift in the yaxis of the Λ significance diagram (Fig. 5) would be only 0.04. A different H_{0} would make a shift along the xaxis, but being constant this does not affect our discussion.
5.1. Shapes
In the adaptive density threshold method, the superclusters were identified using a 8 h^{−1} Mpc kernel. Once the superclusters were extracted, the Minkowski functionals were applied to the smoothed (with a 4 h^{−1} Mpc kernel) luminosity density field of the galaxies in each supercluster. This approach efficiently sorts the general morphologies from the large variety of supercluster shapes. Smoothing also averages the shapes. Although we used a tight criterion for the shapefinder vectors K_{1, 2} to select the most spherical superclusters, all four QSs were found to be flattened at least in one axial direction when only galaxy positions were used to determine the ellipticity. Therefore, we calculated the covariance matrix and eigenvalues for the galaxy distribution in each supercluster in order to define the principal semiaxes of the ellipsoid.
The superclusters SC 550 and SC 849 are found to be considerably elongated systems (namely c > a, c > b and their ), and may be interpreted as prolate spheroids, although the β value indicates a mild triaxiality. The superclusters SC 524 and SC 1832 more closely resemble triaxialtype ellipsoids as they have and 0.53, respectively.
We note that 3D ellipse fitting to the galaxy positions results in more precise shapes than the method based on the combination of the luminosity density and the Minkowski functionals. In the former method, 3D ellipse fitting, galaxies are considered as points, while the latter method instead follows the shape of the galaxy luminosity field, which represents the mass and therefore the gravitational potential field of the supercluster. Presumably, the real shape of the gravitational potential field of a QS supercluster is somewhere between the results given by these two different methods. We also tested the convex hull method; it gives more spherical shapes for all QSs.
5.2. Zerogravity radius for nonspherical spheroids
In Fig. 5, we assume that the QS superclusters are spherical. However, our analysis of the shapes (Sect. 5.1) shows that the QSs are relatively prolate or triaxial spheroids. We know that for a flattened (oblate) system, the zerogravity balance line log(⟨ρ_{M}⟩/ρ_{Λ}) = 0.3 (ZG) will be shifted down. In the limiting case of a homogenous infinitesimally flat disc, the balance line ZG should be shifted down by 0.37 so that log(⟨ρ_{M}⟩/ρ_{Λ}) = − 0.07 (Teerikorpi et al. 2015). In Fig. 6, this line is marked with the symbol ZG_{e}. Not only the amount of flatness, but also the matter distribution inside a nonspherical system affects the zerogravity balance line. In practise, this means that all systems have their own individual zerogravity balance line and thus the zerogravity radius. A more detailed mathematical description of how these properties shift the position of the balance line ZG in the Λ significance diagram is given in the appendix. Here we give a brief summary.
Fig. 6. Λ significance diagram for four QS superclusters. The purple horizontal line shows the shapecorrected ZG balance line for QS 550 and QS 1832 and the green line shows the ZG value for QS 849 and QS 524. Inclined lines show the QS masses (see Table 2). Arrows show estimated upper limits for the QSs in the graph. The log(⟨ρ_{M}⟩/ρ_{Λ}) ZV line is shown as in Fig. 5. The ZGe shows the theoretical limit of the zerogravity balance line for completely flat systems. 
At first, we calculate the force factor f, which denotes the force ratio between a spheroid (oblate or prolate) and an ideal spherical system with the same mass and radius (along the longest axis of the system). This force factor is calculated for different combinations of axial ratios and density profiles. For oblate systems, the force factor f(e, n) can be described as a function of the flatness (defined as e^{2} = 1 − (b/a)^{2}, a = c) and the matter density profile ρ = w^{n} of the system. This results in the balance line Δ⟨ρ_{M}⟩/ρ_{Λ} which shifts smoothly as a function of the flatness. The shift is larger for systems with a flatter shape and flatter density profile (highlighted in Fig. A.2).
Next we estimate how the balance line ZG will be shifted for prolate systems. First, we consider the force factor f_{α} for prolate homogenous spheroids with different axial ratios (α = 0 denotes homogenous mass distribution). We also consider different cases where a fraction of the total mass of a system is located at its centre and the rest is distributed uniformly. Presumably, the real distribution is somewhere between the cases α = 0 and α = 0.5 (half of the mass at the centre; see Appendix A.2 and Table A.1).
In Fig. 6 we show how the balance line ZG of each QS is shifted when α = 0.5 is assumed. In such a case, only QS 849 is clearly below its balance line ZG (green horizontal line) and therefore it is also outside of its R_{ZG}. Three other QSs are above their individual ZG line. Our analysis of the axial ratios in Sect. 5 showed that QS 524 and QS 1832 are relatively triaxial while QS 550 is quite a prolate system. Thus, QS 524 and QS 1832 are not pure prolate or oblate systems, but between these two types of spheroids. Figure A.2 shows the general trend that the shift Δ⟨ρ_{M}⟩/ρ_{Λ} increases when moving from spherical to oblate systems and further to prolate ones. In Fig. 5, the QSs are considered as purely spherical systems and those in Fig. 6 as purely prolate systems. This indicates that the actual balance lines ZG for QS 524 and QS 1832 are situated between these two extremes. In that case, QS 1832 might not be inside its zerogravity radius R_{ZG}.
It is also important to note that the virial masses of the component galaxy groups (and hence the supercluster masses) might be systematically underestimated. This is because the size and virial mass of a system will be undervalued if the dark halo is not fully probed by the observed galaxies (Chernin et al. 2012). The possible amount of mass missed in this way was studied by Teerikorpi et al. (2018) who concluded that, in comparison with the mass up to R_{200}, the maximum mass could be roughly a factor of two larger. Kashibadze & Karachentsev (2018) considered the missing halo mass around galaxy groups, including also the possible retarded outflow component up to the EinsteinStraus distance. These authors ‘calibrated’ this fraction using the mass and average matter density determinations in the local Universe, as compared with the average global dark matter density. The conclusion was that the total mass of a typical group, including its extended halo, could be a factor of three larger than the usually determined virial mass. This is shown in Fig. 6 with arrows pointing to the positions in the Λ significance diagram when the supercluster masses are multiplied by a factor of three, corresponding to a plausible upper limit. As mentioned before, previous studies have shown that the group mass estimation method adopted by us also systematically underestimates masses. If the masses of the QS superclusters were for example two times larger than adopted here, then they would be found close to or above the ZV limit line for spherical systems. This corresponds to the turnaround radius, the boundary where the system completely decouples from the darkenergydriven expansion and virialized. All this means that log(⟨ρ_{M}⟩/ρ_{Λ}) values in Fig. 6 can be considered as lower limits and it is probable that at least QS 524 and QS 550 are gravitationally bound objects. From a dynamical point of view, these objects represent intermediate systems between clusters and superclusters in general.
While zero gravity puts a minimum limit on the overdensity of a gravitationally bound system with mass and size in the current epoch, it does not necessarily prevent the system from expanding in the future. The futures of collapsing systems in the simulations and observations have been studied previously, by for example Dünner et al. (2006), Luparello et al. (2011), Chon et al. (2015), and Gramann et al. (2015). In these studies, the authors showed that, to ensure the collapse of the system in the future, the characteristic overdensity of the systems at the present epoch (Ω_{m} = 0.3 and Ω_{Λ} = 0.7) should be at least ⟨ρ⟩/ρ_{m} ≈ 7.8 or ⟨ρ⟩/ρ_{c} ≈ 2.4 or ⟨ρ⟩/ρ_{Λ} ≈ 3.4. Simulations show that in the future, such systems will get more and more spherical and finally collapse. The value of the density contrast for the future collapsing system at the present epoch log(⟨ρ⟩/ρ_{Λ})≈0.5 is somewhat higher than the values found for QS 524 log(⟨ρ⟩/ρ_{Λ}) = 0.43 (see Table 2 and Fig. 6). This means that we cannot predict whether or not QS 524 will collapse in the future.
6. Discussion
6.1. Gas in superclusters
The warmhot intergalactic medium (WHIM) is expected to be distributed in relatively lowdensity nonvirialized stuctures within superclusters and filaments of galaxies (Cen & Ostriker 1999). In principle, the allsky ymaps of the diffuse gas can be used to trace WHIM. However, yvalues of at least two orders lower as compared with those of clusters decrease the signal to the noise level, hampering the analysis.
To avoid this problem, Tanimura et al (2019) resorted to statistical analysis and stacked the Planckymaps of the regions of 580 superclusters. After careful masking of the clusters and groups they reported the first detections of the intercluster gas in superclusters at the 2.1σ level. In general, superclusters have very multimodal shapes, which unfortunately increases the noise. We realised that this effect is minimised for QS superclusters. Secondly, we found above that the QS superclusters are in or near to the ZG radius and thus at the gravitationally bound stage. For such systems, the density contrast (the mean matter density of the system over the mean matter density of the Universe) ⟨ρ⟩/ρ_{m} ≈ 5 or ⟨ρ⟩/ρ_{c} ≈ 1.5 is required (Gramann et al. 2015). According to simulations, this value approximately corresponds to the baryon overdensity ⟨ρ/ρ_{b}⟩≈10 (Cen & Ostriker 2006) or ⟨ρ/ρ_{c}⟩≈0.2 (Haider et al. 2016) of the system and this in turn corresponds to the lower limit of the baryon density where most of the WHIM mass between temperatures 10^{5} − 10^{7} K should be located. Therefore, it would be interesting to see whether or not the QSs give a positive SunyaevZeldovich signal.
With this aim in mind, we follow the procedure used by Tanimura et al (2019): First, we stack the four supercluster positions in the Planckymap (Planck Collaboration XXII 2016) in taking into account the cluster mask, the point source mask, and the ICB mask of Planck. We then compute the aperture photometry of the inner part. Finally, 2000 random stacks are made to build the probability distribution of the photometry.
Figure 7 shows the result for the four QSs using three different centre definitions compared with random positions. The horizontal error bars are obtained by bootstrap. The photometry is taken at half of the diameter (see Table 1); this was found to give the most prominent signal. According to our analysis (Fig. 6), the radius R of the QSs approximately corresponds to the zerogravity radius of these systems. Their zerovelocity radius is about half of the ZG radius. Figure 7 shows that the errors are large and the signal is marginal (about 1.5σ detection). However, even this result is promising, because we have used only four objects for stacking.
Fig. 7. ysignal distribution of the random sample superimposed by the four stacked QS superclusters. G: location of the galaxy near the highest density peak in the supercluster, D: position of the density peak of the supercluster, L: luminosity centre of the supercluster. 
6.2. Dark matter in superclusters
The masstolight ratios determined for clusters and superclusters show large variations, from about 150 to several hundred (e.g. Bahcall et al. 1995; Carlberg et al. 1996; Gavazzi et al. 2004; Schirmer et al. 2011; Bahcall & Kulier 2014; Einasto et al. 2015, 2016). Despite the scatter, the results imply that the mean masstolight ratio ⟨M/L⟩_{r} increases as a function of scale, and even at cluster scales begins to approach the Universal value (Bahcall & Kulier 2014). The fact that ⟨M/L⟩ remains about the same from clusters to superclusters indicates that the superclusters do not contain a significant amount of additional dark matter between the member clusters (we considered this point when estimating supercluster masses from the summed dynamical masses of the member systems). According to Blanton et al. (2003), the SDSS DR7 mean luminosity density in the rband is Mpc^{−3}. Using this value and with Ω_{m} = 0.3, Ω_{Λ} = 0.68, the mean rband masstoluminosity ratio on cosmologically representative scales will be ⟨M/L⟩_{r} ≃ 330 ± 40 (M_{⊙}/L_{⊙}). This value is close to the average supercluster masstolight ratio ⟨M/L⟩_{r} ≃ 319 (M_{⊙}/L_{⊙}) found in Sect. 4. A considerably higher ⟨M/L⟩_{r} value was found for the smallest and the most spherical superclusters in the sample. This may indicate that the QS superclusters contain more dark matter than superclusters in general, in analogy with dwarf galaxies on much smaller scales.
7. Summary
The elongated distribution of data points in the diagram of Fig. 5 has, not unexpectedly, a similar slope to that of the de Vaucouleurs (de Vaucouleurs 1970) universal density–radius relation (−1.7), which he showed to extend from compact galaxies up to the Local supercluster and the local Universe as known at the time. A rather similar result was obtained by Karachentsev (1968) from pairs of galaxies up to superclusters. This relation, showing that smaller systems have a higher mean density, evidently reflects the hierarchical structure of the galaxy Universe, suggesting that the groups appearing in the diagram, for example, are actually parts of larger systems and are therefore naturally denser than their host systems, on the average.
The SDSS groups as well as prominent clusters in general are usually located deep within the gravitydominated radius (R ≪ R_{ZV}), which is also a requirement for them to be virialized. The virialized region starts when log⟨ρ_{M}⟩/ρ_{Λ} > 0.78. The systems above the zerovelocity radius line are totally within the collapsing region. The outer parts of the systems below the line have not yet been retarded down to zero velocity. For example, (Einasto et al. 2020) estimated that the R_{ZG} scale for the supercluster A2142 is about 11.5 Mpc.
In reality, there could be gravitationally bound superclusters that are not spherical, but have some other type of morphology. In any case, it is interesting that the superclusters with almost spherical morphology appear to be gravity dominated. This is an extension of the general trend: the gravitationally bound celestial bodies, from planets and stars to galaxies and galaxy clusters, tend to exhibit spherical or nearspherical symmetries or axialsymmetries, which may be traced backed to the isotropic effect of gravity around a mass concentration. This is also the case for systems still collapsing under gravity.
Dynamically, the largescale structures of the Universe evolve under an interplay between two factors, dark energy and gravitational force. We used a sample of 1313 superclusters constructed using the 3D luminosity density field and an adaptive threshold method (Liivamägi et al. 2012). This method makes it possible to determine for each supercluster its own density threshold limit. As the cosmological density field and peculiar velocity fields are closely related, we assume that the adaptive density threshold method may provide a way to extract superclusters that are physically more interconnected than systems found using a fixed density threshold method. After consideration of various selection effects, we selected a subsample of 65 superclusters for our final analysis. Using the Minkowski functionals, we studied supercluster morphologies as a function of their physical properties, such as total mass, luminosity, and volume. We extracted a small sample of the most spherical (quasispherical) superclusters. The Λ significance diagram was originally intended for spherical systems (similarly to the spherical collapse model). However, we also found simple analytical solutions for nonspherical spheroids in order to also characterise their dynamical state using this diagram.
The results of our analysis of superclusters extracted with an adaptive local threshold density method can be summarised as follows: (1) Superclusters with roughly spherical or pancaketype shapes are smaller, less massive, less luminous, and contain less galaxies and groups than other superclusters. (2) In our sample of 65 superclusters, we found four potential QS superclusters indicating that such systems are relatively rare, representing just a few percent of all superclusters. These systems have relatively high masstolight ratios, which may indicate that these small and exceptionally spherical superclusters contain relatively large amounts of dark matter. The QSs found in this work are identified as highdensity regions of the known extended superclusters defined with a fixed density level. (3) The thermal SunyaevZeldovich effect signal found using the four stacked QS superclusters is not formally significant, but promising enough to encourage a followup study of selected superclusters and stacking of their ysignals. For example, the number of systems can be increased by allowing a wider range for the K1 and K2 parameters and the clumpiness of the system. (4) Using the Λ significance diagram and lower limits for the masses of the QSs, we found that two of them, QS 524 and QS 550, and perhaps also QS 1832, are dynamically bound systems that form their own entity. If so, then these objects are among the largest gravitationally bound systems found to date.
QSs are large, symmetrically bound systems with a smooth matter distribution. This indicates special smooth initial conditions and a ‘quiet’ formation history. QSs represent a special class of giant systems that are dynamically between large gravitationally unbound superclusters and clusters of galaxies in an equilibrium configuration. Comparing our observational results with data from large cosmological simulations could be an interesting cosmological test.
Acknowledgments
We thank the referee for a report which helped us to improve the paper. We thank Juhan Liivamägi for helpful discussions. This work was supported by grants of the Vilho, Yrjö and Kalle Väisälä Foundation of the Finnish Academy of Science and Letters. The present study was also supported by the ETAG project PRG1006, and by the European Structural Funds grant 79 for the Centre of Excellence “The Dark Side of the Universe” (TK133). We applied in this study the R statistical environment (Ihaka & Gentleman 1996).
References
 Bahcall, N. A., & Kulier, A. 2014, MNRAS, 439, 2505 [Google Scholar]
 Bahcall, N. A., Lubin, L. M., & Dorman, V. 1995, ApJ, 447, 81 [NASA ADS] [CrossRef] [Google Scholar]
 Basilakos, S. 2003, MNRAS, 344, 602 [NASA ADS] [CrossRef] [Google Scholar]
 Basilakos, S., Plionis, M., & RowanRobinson, M. 2001, MNRAS, 323, 47 [Google Scholar]
 Blanton, M. R., Hogg, D. W., Bahcall, N. A., et al. 2003, ApJ, 592, 819 [Google Scholar]
 Byrd, G. G., Chernin, A. D., Teerikorpi, P., & Valtonen, M. J. 2012, Paths to Dark Energy: Theory and Observation (Berlin: de Gruyter) [CrossRef] [Google Scholar]
 Böhringer, H., & Chon, G. 2021, A&A, 656, A144 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Carlberg, R. G., Yee, H. K. C., Ellingson, E., et al. 1996, ApJ, 462, 32 [NASA ADS] [CrossRef] [Google Scholar]
 Cen, R., & Ostriker, J. P. 1999, ApJ, 514, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Cen, R., & Ostriker, J. P. 2006, ApJ, 650, 560 [Google Scholar]
 Chernin, A. D. 2001, PhysicsUspekhi, 44, 1099 [NASA ADS] [CrossRef] [Google Scholar]
 Chernin, A. D., Teerikorpi, P., & Baryshev, Yu. V. 2006, A&A, 456, 13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Chernin, A. D., Teerikorpi, P., Valtonen, M. J., et al. 2012, A&A, 539, A4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Chernin, A. D., BisnovatyiKogan, G. S., Teerikorpi, P., et al. 2013, A&A, 553, A101 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Chon, G., Böhringer, H., & Zaroubi, S. 2015, A&A, 575, L14 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 CostaDuarte, M. V., Sodré, L., Jr., & Durret, F. 2011, MNRAS, 411, 1716 [NASA ADS] [CrossRef] [Google Scholar]
 Davis, M., Huchra, J., Lathman, D. W., & Tonry, J. 1982, ApJ, 253, 423 [NASA ADS] [CrossRef] [Google Scholar]
 Dupuy, A., Courtois, H. M., Florent, D., et al. 2019, MNRAS, 489, L1 [CrossRef] [Google Scholar]
 Dünner, R., Araya, P. A., Meza, A., & Reisenegger, A. 2006, MNRAS, 366, 803 [NASA ADS] [CrossRef] [Google Scholar]
 Einasto, J., Jõeveer, M., & Saar, E. 1980, MNRAS, 193, 353 [NASA ADS] [CrossRef] [Google Scholar]
 Einasto, M., Einasto, J., Tago, E., Müller, V., & Andernach, H. 2001, AJ, 122, 2222 [NASA ADS] [CrossRef] [Google Scholar]
 Einasto, J., Hütsi, G., Einasto, M., et al. 2003, A&A, 405, 425 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Einasto, M., Saar, E., & Liivamägi, L. J. 2007, A&A, 476, 697 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Einasto, M., Liivamägi, L. J., Tago, E., et al. 2011a, A&A, 532, A5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Einasto, M., Liivamägi, L. J., Saar, E., et al. 2011b, A&A, 535, A36 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Einasto, M., Gramann, M., Liivamägi, L. J., et al. 2015, A&A, 580, A69 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Einasto, M., Lietzen, H., Gramann, M., et al. 2016, A&A, 595, A70 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Einasto, M., Deshev, B., Tenjes, P., et al. 2020, A&A, 641, A172 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Einasto, M., Kipper, R., Tenjes, P., et al. 2021, A&A, 649, A51 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Einasto, M., Tenjes, P., Gramann, M., et al. 2022, A&A, 666, A52 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Erdogdu, P., Lahav, O., Zaroubi, S., et al. 2004, MNRAS, 352, 939 [NASA ADS] [CrossRef] [Google Scholar]
 Frieman, J. A., Turner, M. S., & Huterer, D. 2008, ARA&A, 46, 385 [Google Scholar]
 Gavazzi, R., Mellier, Y., Fort, B., Cuillandre, J.C., & DantelFort, M. 2004, A&A, 422, 407 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Gramann, M., Einasto, M., Heinämäki, P., et al. 2015, A&A, 581, A135 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Haider, M., Steinhauser, D., Vogelsberger, M., et al. 2016, MNRAS, 457, 3024 [Google Scholar]
 Hofmeister, A. M., Criss, R. E., & Criss, E. M. 2018, Planet. Space Sci., 152, 68 [NASA ADS] [CrossRef] [Google Scholar]
 Ihaka, R., & Gentleman, R. 1996, J. Comput. Graph. Stat., 5, 299 [Google Scholar]
 Jaaniste, J., Tago, E., Einasto, M., et al. 1998, A&A, 336, 35 [NASA ADS] [Google Scholar]
 Karachentsev, I. D. 1968, Publ. Byurakan Obs., 39, 76 [NASA ADS] [Google Scholar]
 Kashibadze, O. G., & Karachentsev, I. D. 2018, A&A, 609, A11 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Liivamägi, L. J., Tempel, E., & Saar, E. 2012, A&A, 539, A80 [Google Scholar]
 Luparello, H., Lares, M., Lambas, D. G., & Padilla, N. 2011, MNRAS, 415, 964 [Google Scholar]
 MacMillan, W. D. 1930, The Theory of the Potential (New York: McGrawHill) [Google Scholar]
 Nadathur, S., & Crittenden, R. 2016, ApJ, 830, L10 [Google Scholar]
 Old, L., Skibba, R. A., Pearce, F. R., et al. 2014, MNRAS, 441, 1513 [Google Scholar]
 Old, L., Wojtak, R., Mamon, G. A., et al. 2015, MNRAS, 449, 1897 [Google Scholar]
 Planck Collaboration XVI. 2014, A&A, 571, A16 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XXII. 2016, A&A, 594, A22 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Ragone, C. J., Muriel, H., Proust, D., et al. 2006, A&A, 445, 819 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Saar, E., Martínez, V. J., Starck, J.L., & Donoho, D. L. 2007, MNRAS, 374, 1030 [NASA ADS] [CrossRef] [Google Scholar]
 Shandarin, S. F., Sheth, J. V., & Sahni, V. 2004, MNRAS, 353, 162 [NASA ADS] [CrossRef] [Google Scholar]
 Sahni, V., Sathyaprakash, B. S., & Shandarin, S. F. 1998, ApJ, 495, L5 [NASA ADS] [CrossRef] [Google Scholar]
 Schirmer, M., Hildebrandt, H., Kuijken, K., & Erben, T. 2011, A&A, 532, A57 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Sheth, R., & Diaferio, A. 2011, MNRAS, 417, 4 [Google Scholar]
 Tago, E., Saar, E., Tempel, E., et al. 2010, A&A, 514, A102 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Tanimura, H., Aghanim, N., Douspis, M., et al. 2019, A&A, 625, A67 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Teerikorpi, P., Heinämäki, P., Nurmi, P., et al. 2015, A&A, 577, A144 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Teerikorpi, P., Heinämäki, P., & Chernin, A. D. 2018, Astron. Nachr., 339, 705 [NASA ADS] [CrossRef] [Google Scholar]
 Tempel, E., Saar, E., Liivamägi, L. J., et al. 2011, A&A, 529, A53 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Tempel, E., Tago, E., Liivamägi, L. J., et al. 2012, A&A, 550, A106 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Tempel, E., Tamm, A., Gramann, M., et al. 2014, A&A, 566, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Tully, R. B., Courtois, H., Hoffman, Y., & Pomarède, D. 2014, Nature, 513, 71 [Google Scholar]
 de Vaucouleurs, G. 1970, Science, 167, 1203 [NASA ADS] [CrossRef] [Google Scholar]
 Woltjer, L. 1967, Lect. Appl. Math. (Am. Math. Soc.), 9, 1 [NASA ADS] [Google Scholar]
Appendix A: Zerogravity radius for oblate and prolate matter spheroids
The Λ significance diagram is originally intended for spherically symmetric systems, but it also has limited use for flattened or elongated spheroids as we show here. A system cannot be bound as a whole if the antigravity force between its centre and a part of the edge is greater than the gravity pulling the edge inwards. Here we consider the gravity force at the rim of a disc (an oblate spheroid) and at the end of a ‘cigar’ (a prolate spheroid), in comparison with the corresponding antigravity force.
It is expected that, at the rim of the disc and at the end of the cigar, the gravity force is larger than at the surface of a sphere of the same mass and the radius equal to the corresponding (half) major axis of the spheroid. This may be understood by noting that when a sphere is transformed into a spheroid, the average distance of the mass points from the considered fixed point on the surface of the sphere is decreased. Also, the projections of the forces along the considered axis are larger. The change is bigger for a prolate spheroid.
A.1. Oblate spheroid
A disc exerts a higher gravity force on a particle at its rim than a sphere of the same mass and radius. If the ratio of these forces is f(> 1), then the balance condition at the equatorial rim R = a is
or
Here, ⟨ρ_{M}⟩ is the mean density caused by the mass M within the sphere of radius a (< the mean density of the disc). Formally, it may be used in the Λ significance diagram, provided that the ‘gravity = antigravity’ line is correspondingly lowered for a flattened system, which here we take to be a spheroid with the axes a = c and b < a. We consider the case of homogenous mass distribution as well as the case of the density decreasing from the centre.
We calculate the factor α = f(e, n) for spheroids characterised by the eccentricity e (defined as e^{2} = 1 − (b/a)^{2}) and the density decreasing from the centre as the power law w^{n}, where −3 < n ≤ 0 and w parametrizes the spheroidal shells inside the outer surface. The value n = 0 corresponds to a homogenous spheroid. It is known (see e.g. Woltjer 1967; Hofmeister et al. 2018) that then the gravity force at the equator, at the edge of the spheroid with the equatorial radius a is
hence
For different values of n (here 0, −2), this allows us to calculate the factor f(e, n) and hence to see how sensitive the location of the ‘gravity = antigravity’ line is when the eccentricity changes from the value = 0 (a sphere) towards more flattened spheroids. These are:
For a flat disc, with the density increasing inwards, the ratio f(1, n) < 3π/4, and the case of constant density gives f(1, 0) = 3π/4 = 2.35. In this limiting case, the ‘gravity = antigravity’ balance line in the Λ significance diagram would be shifted down to the value ⟨ρ_{M}⟩/ρ_{Λ} = −0.07 in order to roughly take very flattened systems into account. Figures A.1 and A.2 show results for less flattened spheroids as a function of the axis ratio b/a = (1 − e^{2})^{1/2}.
Fig. A.1. Force excess factor f vs. the axis ratios b/a (oblate) and a/c (prolate) of a spheroid for different density laws (n = 0 and α = 0 mean homogeneity). The corresponding values of log⟨ρ_{M}⟩/ρ_{Λ}, when there is the ‘gravity = antigravity’ balance at the maximum extension of the spheroid, are shown for three values of f. We note that for the spherical case (b/a = a/c = 1) f = 1, and log⟨ρ_{M}⟩/ρ_{Λ} = 0.3, as generally used in the Λ significance diagram. 
We note that when the zerogravity condition is fulfilled at the equatorial rim of the spheroid, this guarantees that the whole spheroid lies in the gravitydominated region, because at the polar point the gravity force is greater than at the rim and the distance to the centre is smaller.
A.2. Prolate spheroid
Analogously to our analysis of oblate spheroids, we want to calculate the ratio f of the force at the end of an elongated spheroid (having the axes a = b and c > a) and the force at the surface of a sphere of the radius c, both having the same mass M. The density of the spheroid is then M/((4/3)πa^{2}c).
MacMillan (1930) discussed a general case of the potential field around a prolate homogenous spheroid. Here we are satisfied with calculating the gravitation force at the end of such a spheroid (density = ρ), which also serves to check the result by MacMillan in this special case. Using the general potential function as derived by MacMillan (1930) in his equation 32.13, by derivation along the zaxis at z = c we obtain the following formula for the present case:
Here we have in addition used the expression M = (4/3)πa^{2}c)ρ for the mass of the homogenous spheroid. The eccentricity is e^{2} = 1 − (a/c)^{2}
Using this expression, the excess factor f for a homogenous prolate spheroid becomes simply
It may be confirmed that f → 1, when e → 0 (a sphere). Also, we checked the formula for F by direct numerical integration, summing the forces from thin circular discs, making up the spheroid. The resulting total force is found to be the same.
If a more realistic mass distribution is such that a fraction α of the total mass M is situated as a point mass at the centre and the remaining mass fills the spheroid uniformly, then the ratio f_{α} is equal to
where f_{0} comes from the previous equation. Table A.1 shows the ratios f_{0} and f_{α} (α = 0.5) for different axis ratios a/c. As expected, the excess factor f is larger at the end point of an elongated spheroid than at the rim of a flattened spheroid of the same maximum radius. This is because in the former case the mass is closer, on average, to the surface point considered. Figures A.1 and A.2 show the results of our calculations as graphs.
Fig. A.2. Shift of the ‘gravity = antigravity’ balance line in the Δlog⟨ρ⟩/ρ_{Λ} vs. log R diagram as a function of the axis ratio b/a (oblate) and a/c (prolate). For example, when a/c is greater than 0.5, the shift is less than about 0.3 for all the cases considered. 
Excess force factors f_{0} and f_{0.5} for different a/c.
Appendix B: Quasisphericals in the plane of the sky
Fig. B.1. Distribution of galaxies in the regions of superclusters QS 550, QS 849, QS 524, and QS 1832 in the plane of the sky. Orange dots correspond to supercluster member galaxies and blue dots show the galaxies outside the superclusters. A red spherical circle illustrates the spherical diameter of a supercluster (see table 1). Abell clusters are labeled. 
All Tables
All Figures
Fig. 1. Mass distribution and size–distance scatter plot of the superclusters. The upper figure shows the mass histogram of the final sample of the 65 superclusters and the lower figure shows supercluster diameter as a function of distance. 

In the text 
Fig. 2. Correspondence between the vectors K_{1} and K_{2} and different shapes. 

In the text 
Fig. 3. Coloured contours show the distribution of the values of supercluster properties in the K1 − K2 plane. Quasispherical superclusters K_{1}, K_{2} < (0.05, 0.05) are shown with red points, white points are all the superclusters. Luminosities are in units of 10^{10} L_{⊙}, masses in 10^{15} M_{⊙}, volumes in Mpc^{3}, and diameter in megaparsecs (Mpc). The dark line indicates K_{1} = K_{2} and divides the morphologies between pancakes and filaments. 

In the text 
Fig. 4. Distribution of galaxies in the regions of superclusters QS 550, QS 849, QS 524, and QS 1832 viewed from three different orientations (x − y, x − z, y − z) in Cartesian coordinates. Distribution of galaxies in the plane of the sky are shown in Appendix B. 

In the text 
Fig. 5. Λ significance diagram for SDSS clusters of galaxies (purple) with QS superclusters marked in green. The Laniakea supercluster, the central regions of the supercluster A2142 and the Coma cluster with three different mass and size estimates are shown. Inclined lines show the effect of different example mass values. Green vertical lines show the values of the R_{ZV} (zero velocity radius), R_{ZG} (zero gravity radius), and the R_{ES} (the EinsteinStraus radius for the mass 10^{14} M_{⊙}) and the horizontal lines show logarithmic values for the constant ratios of the ⟨ρ_{M}⟩/ρ_{Λ} = 6 (ZV), 2 (ZG), and 3/7 (ES). 

In the text 
Fig. 6. Λ significance diagram for four QS superclusters. The purple horizontal line shows the shapecorrected ZG balance line for QS 550 and QS 1832 and the green line shows the ZG value for QS 849 and QS 524. Inclined lines show the QS masses (see Table 2). Arrows show estimated upper limits for the QSs in the graph. The log(⟨ρ_{M}⟩/ρ_{Λ}) ZV line is shown as in Fig. 5. The ZGe shows the theoretical limit of the zerogravity balance line for completely flat systems. 

In the text 
Fig. 7. ysignal distribution of the random sample superimposed by the four stacked QS superclusters. G: location of the galaxy near the highest density peak in the supercluster, D: position of the density peak of the supercluster, L: luminosity centre of the supercluster. 

In the text 
Fig. A.1. Force excess factor f vs. the axis ratios b/a (oblate) and a/c (prolate) of a spheroid for different density laws (n = 0 and α = 0 mean homogeneity). The corresponding values of log⟨ρ_{M}⟩/ρ_{Λ}, when there is the ‘gravity = antigravity’ balance at the maximum extension of the spheroid, are shown for three values of f. We note that for the spherical case (b/a = a/c = 1) f = 1, and log⟨ρ_{M}⟩/ρ_{Λ} = 0.3, as generally used in the Λ significance diagram. 

In the text 
Fig. A.2. Shift of the ‘gravity = antigravity’ balance line in the Δlog⟨ρ⟩/ρ_{Λ} vs. log R diagram as a function of the axis ratio b/a (oblate) and a/c (prolate). For example, when a/c is greater than 0.5, the shift is less than about 0.3 for all the cases considered. 

In the text 
Fig. B.1. Distribution of galaxies in the regions of superclusters QS 550, QS 849, QS 524, and QS 1832 in the plane of the sky. Orange dots correspond to supercluster member galaxies and blue dots show the galaxies outside the superclusters. A red spherical circle illustrates the spherical diameter of a supercluster (see table 1). Abell clusters are labeled. 

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.