Abstract
We present a new statistical analysis of the large-scale stellar mass distribution in the Sloan Digital Sky Survey (data release 7). A set of volume-limited samples shows that the stellar mass of galaxies is concentrated in a range of galaxy luminosities that is very different from the range selected by the usual analysis of galaxy positions. Nevertheless, the two-point correlation function is a power-law with the usual exponent , which varies with luminosity. The mass concentration property allows us to make a meaningful analysis of the angular distribution of the full flux-limited sample. With this analysis, after suppressing the shot noise, we extend further the scaling range and thus obtain and a clustering length . Fractional statistical moments of the coarse-grained stellar mass density exhibit multifractal scaling. Our results support a multifractal model with a transition to homogeneity at about .
1. Introduction
It is claimed that we are entering the era of precision cosmology, in which the fundamental parameters of the universe are known within a few percent precision and only remain to progressively refine them [1, 2]. In addition to the global parameters of the FLRW model, we have the parameters that determine the primordial density perturbations and hence the formation of large-scale structure. A particularly interesting parameter is the amplitude of the primordial density perturbations, usually measured in terms of the present linear-theory mass dispersion on a scale of h−1 Mpc, named . The scale of 8 h−1 Mpc comes from Peebles’ observation that galaxy counts on this scale have an rms fluctuation approximately equal to one [3]. In general, the mass dispersion over the length , , defines a scale such that . This scale separates the linear and quasi-homogeneous regime of the evolution of density fluctuations from the nonlinear regime of strong clustering.
A similar scale is the distance such that , where is the reduced galaxy-galaxy correlation function. Actually, this function is well approximated by the power-law:for distances not much larger than [3]. The scale such that the rms dispersion of galaxy counts is equal to one is related to by a -dependent factor [3]. Historically, the power-law form of was deduced from the angular positions of galaxies, yielding and h−1 Mpc [4, 5]. Nowadays, good galaxy redshift surveys are available, but the angular positions are still useful. In particular, the angular correlation function is used, in combination with other data, to determine precision values of [2, 6].
The amplitude of the primordial density perturbations is not theoretically constrained and determines the present scale of transition to homogeneity, which can have any value, once given the global cosmological parameters. The transition to homogeneity has been the subject of numerous studies, some of them motivated by the bold proposal that the scale of homogeneity is too large to be accessible or even that no such transition exists [7–9]. This proposal, namely, that the mass distribution is fully scale invariant, at all scales in a Newtonian cosmology, has been called the “fractal universe” [3]. Actually, equation (1) corresponds to a fractal distribution on scales , for any magnitude of .
The scale and specific form of transition to homogeneity determine the size of the largest structures, although these structures can have a length much larger than [10]. An examination of the literature shows values of the scale of transition to homogeneity that go from the standard value of 5 h−1 Mpc [3] to more than 20 h−1 Mpc [11–14]. The range that these values span is quite long. In contrast, the current precision values of are assigned a relative error of about one percent [2]. This precision is surprisingly good, in comparison, given that determines the scale of transition to homogeneity.
Most of the literature about the large-scale structure of the universe based on the statistical analysis of galaxy catalogs is limited to the distribution of the galaxy positions. We make here, following up on Gaite [15], a new statistical analysis of the large-scale structure, adding an important ingredient: the stellar masses of galaxies, that is to say, we study the large-scale distribution of stellar mass. The distribution of stellar mass in the SDSS has already been studied by Li and White [16], employing the projected correlation function , on scales h−1Mpc ( is the separation perpendicular to the line of sight [17]). Li and White find that is very well represented, over the range h−1 kpc h−1 Mpc, by a power law that corresponds to equation (1) with and h−1 Mpc. These results basically agree with the analysis of the SDSS galaxy positions by Zehavi et al. [18], but Li and White deem the scaling of the stellar mass to be better. At any rate, the stellar mass and galaxy number distributions can be shown to be very different, thus warranting further statistical analysis.
Besides, we must take into account that the “projection” that produces is performed on the correlation function and is not associated to an actual projected distribution but to a set of local orthogonal projections on the tangent planes to the surface of the unit sphere (a rather unwieldy construction from the mathematical standpoint). Arguably, the direct study of the angular projection of the stellar mass distribution is preferable. On the other hand, there are statistical measures not based on the two-point correlation function. Our study combines the analysis of a set of volume-limited samples with an angular analysis and, moreover, involves various statistical measures.
The importance of galaxy masses in the study of galaxy clustering was recognized years ago [8]. Furthermore, Pietronero [8] argued that a full multifractal analysis of galaxy clustering is necessary. A multifractal model was also proposed by Jones et al. [19], without considering the galaxy masses. Pietronero and collaborators initiated the multifractal analysis with masses of galaxies derived from the observed luminosities by assuming a simple mass-luminosity relation [9, 20]. Meanwhile, the large-scale structure has been known to consist not just of clusters but of a web structure [21, 22]. This structure can be described as a multifractal [15, 23]. Multifractality has become essential to the scaling analysis of the distribution of galaxies, although normally without considering the galaxy masses [24, 25].
As we shall see, most of a multifractal structure is preserved in an angular projection. Pietronero [8] already considered some basic properties of the angular projection of a fractal galaxy distribution. This question was further studied by Coleman and Pietronero [9] and Durrer et al. [26], assuming a monofractal model.
Our galaxy data come from the Sloan Digital Sky Survey, data release 7 (SDSS DR7) [27], as provided by the New York University Value-Added Galaxy Catalog (NYU-VAGC) [28]. The NYU-VAGC contains the stellar masses of the galaxies, calculated with the method described by Blanton and Roweis [29] (which gives similar results to the method of [30]). In fact, we use the same data as Li and White [16]; which facilitates the comparison of the respective results. Moreover, the SDSS DR7 is well studied and, in particular, there is a careful analysis of its galaxy-galaxy angular correlation function [31]. We have already employed the same data to make a multifractal analysis of the stellar mass distribution based on the convergence of multifractal spectra for a shrinking scale [15]. The present treatment of scaling laws is much more elaborate, allowing us to compare with the preceding studies of scaling in the SDSS and to study the transition to homogeneity.
Our plan is as follows: we begin in Section 2 with a short discussion of scaling, fractality, and homogeneity in galaxy surveys and of the role played by the scale-dependent mass variance. In Section 3, we apply these concepts to the study of a set of volume-limited samples from the SDSS DR7 and to the analysis of scaling of the respective mass variances. Multifractality is proved by the analysis of fractional statistical moments in Section 4. We next consider the theory of projection of fractal distributions in Section 5, and we connect this theory with the standard theory of the angular two-point correlation function of galaxies (Section 5.1). After this theoretical study, we undertake the analysis of the angular projection of the SDSS DR7 in Section 6, where we deal with the increased shot noise (Section 6.1) and we obtain from the variance of the coarse-grained angular density (Section 6.2). We end with a general discussion (Section 7).
2. Scaling, Fractality, and Homogeneity in Galaxy Surveys
It is pertinent to begin with some general considerations about scaling and homogeneity in galaxy surveys, and, in particular, about how to characterize fractality and how to determine the scale of transition to homogeneity. Let us assume that we can obtain the three-dimensional stellar mass distribution, which requires a redshift survey with galaxy stellar masses and the construction of volume-limited samples.
Taking the two-point correlation function as the basic statistic, the fundamental scaling law is the power-law reduced correlation function in equation (1). This law is supposed to hold when is not small, that is to say, for not much larger than (but it could also be valid for ). The coarse-grained mass fluctuation , namely, the scale-dependent rms fluctuation of the stellar mass in a cell of linear dimension , derives from the reduced two-point correlation function of the mass distribution. If this function follows the power law (1), the mass variance follows a power law with the same exponent, namely,where the form of depends on the shape of the cell [3].
Regardless of any scaling property, must decrease with and the scale of homogeneity can be defined just by the condition . However, this condition is insufficient for homogeneity. Indeed, let us notice that a positive random variable with an rms dispersion equal to its mean value cannot be approximately Gaussian; for example, let us consider the lognormal distribution (see [32]; p. 5). A suitable criterion for homogeneity is the mass variance , as assumed by Gaite [15]; indeed, a lognormal distribution function with this variance is hardly skewed ([32]; p. 5). Whether or not the scale defined by is close to the one defined by depends on how sharp the transition to homogeneity is. To analyze this question, we naturally choose the power-law form in equation (2).
For a spherical cell of radius , in equation (2) is given by Peebles [3]. The scale such that equals a given number is a definite function of . For , is only a little larger than , whereas, for , it is between and (the larger, the smaller is). In conclusion, the scale such that is roughly equivalent to , but the scale where the probability distribution is approximately Gaussian can be several times larger. This observation reveals that the values h−1 Mpc and are compatible with the presence of relatively large structures, for example, cosmic voids of size h−1 Mpc [3]. Actually, even larger structures can be observed, provided that does not fall too rapidly for [10].
After clarifying our concept of homogeneity, let us now recall several aspects of scale invariance and fractal geometry. First of all, let us remark that the fractal geometry of a mass distribution consists in some sort of scale invariance in the strongly inhomogeneous or strong clustering regime. Its most general form is called multifractality and is expressed in terms of the multifractal spectrum or, alternatively, of the Rényi dimension spectrum [33, 34]. The multifractal spectrum is a basic function that represents how mass concentrates, in terms of the number of mass concentrations of a given “strength,” measured by the local dimension, while the Rényi dimensions provide a sort of averaged information. The functions and are related by a Legendre transform in terms of the variables and . The Rényi dimension expresses the scaling of the statistical -moment of the coarse-grained mass distribution and is, in general, a nonincreasing function of . The second-order moment has an intuitive meaning, and the corresponding Rényi dimension is . General -moments and their scaling relations are studied in Section 4.
It is to be remarked that the primary concept of fractal dimension and its various definitions are concerned with the small-scale behavior of sets or mass distributions, and fractality is roughly equivalent to asymptotic scaling in the limit of vanishing scale, as described in standard fractal geometry textbooks [33, 34]. In consequence, multifractal measures are defined for any mass distribution, regardless of its behavior on large scales. An arbitrary mass distribution is not expected to be such that it becomes homogeneous on large scales, tending to a uniformly homogeneous state, as occurs in cosmology because of the cosmological principle. This is why raw statistical -moments are employed in fractal geometry, instead of central moments, which require a globally defined mean density.
To make the cosmological principle compatible with fractal models of strong clustering, Mandelbrot [7] proposed the conditional cosmological principle, in which “every possible observer” is replaced with “every observer located at a material point.” Thus, the natural measure for the fractal analysis of strong clustering is the conditional density, namely, the average density at a distance from an occupied point [7, 9, 20]. It can be expressed as
But, it is conceptually independent of the mean density in spite of the fact that it appears in these expressions.
Motivated by these ideas, some authors look for scaling of the conditional density of the galaxy distribution (or of its integral in the ball of radius ). When calculated with this method, the scale of homogeneity is often large [11, 13]. The conditional density is not used by traditional cosmologists, who are accustomed to the function . If , then the conditional density is conceptually useful but can be expressed in terms of , by equation (3), and it scales as . Naturally, the condition for fractal scaling is , which makes and equivalent. However, the condition is never fulfilled in a sufficiently long range of . Therefore, the values of and obtained from power-law fits of either function differ, namely, is larger when calculated from . Let us illustrate this point with an example, employing a sample from Gaite [15].
Gaite [15] uses a method of multifractal analysis based on coarse-graining cells that are adapted to the SDSS sample geometry. These cells are not spherical nor have a simple shape but have equal volume, which serves to measure their size. Let this volume be . The second moment of the coarse-grained density isand the density variance (mass variance) is
The scaling law (1) is equivalent to the scaling law for the mass variance (2). Likewise, a scaling law for is equivalent to a scaling law for . An example of both the scaling laws for and is displayed in Figure 1, which we now explain.
Figure 1 refers to a volume-limited sample of 1765 galaxies in the redshift range , defined by Gaite [15] and called VLS1. This sample was adequate to compute a quite complete multifractal spectrum and is now useful to test the two options for the scaling law that yields . The solid blue lines in Figure 1 are the graphs of and . The fit to the power-lawin the interval yields and , that is to say, (the power-law fit is a least-squares linear fit of versus ). An analogous fit to the power-law,yields and (). Both fits are represented by dashed red lines in Figure 1. The latter fit, with and Mpc/h, basically agrees with the canonical values of and (from the reduced two-point correlation function of galaxy positions). However, both fits are questionable for calculating , because the fractal regime demands , but the fitted range extends beyond the small scales where this condition holds. Therefore, some intermediate values of and should be more appropriate. These values are uncertain because we have too small a range of scales with both and negligible discreteness effects. These effects make shrink towards , as the steeper left ends of the graphs reflect.
The uncertainties of 20% in and 30% in due to the limited scaling range can be compared with the uncertainties in galaxy positions and masses. Gaite [15] has analyzed how these uncertainties affect the multifractal spectrum of VLS1, with the result that only the right-hand side part of the spectrum is affected (left-hand side of Figure 4 in [15]). This implies that only the (Rényi) dimensions with are affected. Therefore, the evaluation of should not be affected. To confirm it, we calculate now the uncertainty of due to the uncertainty of galaxy positions and masses. To do it, we employ the same procedure as Gaite [15], namely, we generate ten variant samples, with values of redshift and stellar mass given by, respectively, Gaussian and lognormal distributions with variances in accord with the literature. The variant samples yield ten values of for each value of , and we compute for each the standard deviation of those ten values. The result is that the relative error of is always smaller than 4.2% (the corresponding error bars in Figure 1 would hardly be visible).
(a)
(b)
(c)
As regards errors in the linear fits, their magnitude grows as one tries to fit more points, and the natural criterion to stop is to minimize the rms error per degree of freedom (the number of degrees of freedom is the number of points minus two). In most of the fits here, these errors are quite small in comparison with errors from other sources. The relative small uncertainty due to statistical errors in galaxy positions and masses combined with the smallness of errors in the linear fits shows that the major uncertainty in the fractal analysis is due to discreteness effects, that is to say, to limitations imposed by a small fractal scaling range.
Gaite [15] observes that the low-redshift VLS1 obtains a multifractal spectrum that is more complete than the ones obtained from higher-redshift volume-limited samples. We find that higher- samples are also less useful to directly calculate as above. The higher is in a volume-limited sample the more luminous the galaxies in it and the smaller their number density (as we shall discuss when we construct a set of volume-limited samples in Section 3). A smaller number density implies that discreteness effects take over growing ranges of the smaller scales, reducing the interval of scales where the fractal regime can be measured, that is to say, the interval where .
However, we may consider the scaling of just , in the longest possible interval, without demanding that . For example, one can dismiss the upper graph in Figure 1, corresponding to , and think of extending the scaling of below the line . Likewise, equation (1) usually extends to (weak correlations). Besides, we remark that the frequently employed projected correlation function (e.g., the one in [16] or [18]) is a dimensional function that renders obscure the strength of correlations. The range with is the appropriate one for the molecular correlations in statistical physics, and scale invariance takes place in critical phenomena. In this regard, we can expect that equation (1) holds from the small scales where to the large scales where , encompassing the fractal and critical-fluid regimes [10].
At any rate, what we conclude from the analysis of VLS1 is that the statistical analysis of the stellar mass distribution in the SDSS DR7 by means of volume-limited samples is plagued with errors, especially, discreteness errors. An alternative to it is the study of the projected angular distribution, which has two advantages. On the one hand, we avoid the influence of redshift space distortions on the statistics, which is considerable [18]. On the other hand, we have all points in the angular region, whereas volume-limited samples span the same angular region but occupy very different three-dimensional regions (which are awkwardly shaped, in addition). We see next the information obtained from a systematic study of volume-limited samples and leave the study of the angular projection to Section 5.
3. Volume-Limited Samples
To be systematic in the study of the three-dimensional mass distribution, we construct volume-limited samples in successive intervals of absolute magnitude; for example, in the SDSS, one can take intervals in which it changes by one unit. We follow the procedure and conventions of Zehavi et al. [18] (the absolute magnitude is denoted by , not to be confused with our notation of as a mass).
We employ the same SDSS-DR7 data from the NYU-VAGC as in Gaite [15], and we also select the apparent magnitude range . The angular ranges are defined in terms of equal-area coordinates and radians, where the angles and constitute the SDSS coordinate system. The available ranges of our coordinates are and (radians), providing a total solid angle steradians. We further restrict the sample to , obtaining galaxies.
The main characteristics of our volume-limited samples are reported in Table 1. Naturally, the redshift intervals are almost coincident with the ones of Zehavi et al. [18], but the number of galaxies in each sample is larger, because the angular region is larger. The number of galaxies grows with absolute luminosity, up to , with . This growth might suggest that discreteness effects are reduced up to this point, but the number density shrinks and the overall effect is that discreteness progressively hinders the detection of strong clustering, as noticed in Section 2. For instance, the sample with and number density 0.00159 has a volume per galaxy of 629 , that is to say, the equivalent length is 8.6 , as large as the expected value of .
As regards number density, the best sample is the first one, with and number density 0.08550 (this sample is somewhat like the VLS1 studied in Section 2). However, we are studying here the distribution of mass rather than the distribution of individual galaxies and the sample represents a fraction of the total mass density that only amounts to 0.009. Thus, the clustering properties of the full mass distribution are hardly influenced by the galaxies with , in spite of their abundance.
The last variable in Table 1 is the galaxy mass, and we can observe the evident correlation between galaxy mass and absolute luminosity. This correlation is significant in the multifractal analysis, because this analysis distinguishes the strength of mass concentrations, measured by the local dimension . To be precise, the concentration strength is measured in coarse multifractal analysis by the logarithm of the mass in each cell, namely, , ([34], applied to cosmology by [35]). However, galaxies are not equal-size mass concentrations and, actually, the galaxy sizes are not even part of the data. Nevertheless, it is to be expected that the stellar mass in the more luminous samples is more concentrated. Therefore, we also expect that the set of values of and hence the dimension of each sample are smaller the more luminous the galaxies in it.
3.1. Scaling of Cell Mass Variances
Here, we examine if cell mass variances are power-law functions of the cell volume , according to equation (7) and employing the coarse-graining method of Gaite [15] as in Section 2. We focus on the three galaxy samples in Table 1 in the interval , whose contribution to the total mass density is dominant (more than 78%). But, we also examine, to be thorough, the two adjacent samples, to include most of the total mass density as well as most galaxies in the set of samples.
The results of the calculation and plotting of appear in Table 2 and Figure 2. To wit, in Figure 2 are the log-log plots of for the three samples with , and in Table 2 are the numerical results of the log-log linear fits corresponding to the five volume-limited samples of Table 1 with higher mass density. We obtain a nontrivial scaling for intermediate values of , which crosses over, for small , to the trivial scaling corresponding to isolated points (). It is remarkable that the scalings displayed in Figure 2 hold in considerable ranges that go across the point of unit variance, extending from the moderately strong clustering regime to the quasi-homogeneous regime.
Regarding the values of the scaling exponent in Table 2, we can observe a definite growing trend, except in the sample with (an anomaly in this range is also noted by [18]). The length scale , which we can take as a homogeneity scale, also grows with luminosity. The growing trends of both and the homogeneity scale are also obtained in the analysis of galaxy positions (Table 1 of [18]).
In the (moderately) strong clustering regime, the scaling of the mass variances denote fractality, with , the lower values for the larger luminosities. This decreasing trend of fractal dimension as a function of luminosity is expected, as explained at the end of Section 3. Let us notice that the galaxies with , which concentrate most of the mass, correspond to a narrow interval of the scaling exponent, namely, , and hence to a narrow interval of . The concentration of mass in a small range of fractal dimensions is indeed expected in a multifractal [33, 34]. The consequences of the phenomenon of mass concentration for a general statistical analysis are considered in the next section.
4. Fractional Statistical Moments
The analysis in the preceding section is based on the order two statistical moment and, specifically, on the variance . However, the variance is a sufficient statistic only for the normal distribution. Some conclusions about the importance of other moments can be drawn in general, without appeal to scaling arguments.
It has been remarked earlier that a positive random variable with an rms dispersion equal to its mean value is not nearly normal. From the general statistical moment inequality ([36]; p. 55), we deduce a lower bound to the skewness:
Hence, implies that the probability distribution function is very skewed. The lognormal distribution is an example that goes from nearly normal for to very skewed for . Furthermore, this distribution is heavy-tailed, namely, not exponentially bounded, so its moment generating function is ill-defined [32]. For such distributions, high-order moments are not meaningful. Actually, heavy-tailed probability distributions may not have high-order integral moments, but they can be analyzed using fractional low-order moments [37]. This type of analysis connects with the analysis of multifractals, which are mass distributions with very skewed probability functions and are related to the lognormal model [38, 35].
Indeed, a complete multifractal analysis requires the full set of fractional -moments. In the standard method of coarse multifractal analysis [33, 34], a region that contains the mass distribution is coarse-grained with a partition in cells of linear size . Using this partition, fractional statistical moments are defined aswhere the index runs over the set of nonempty cells, is the mass in the cell , and is the total mass. Then, multifractal behavior is given, for , byand the Rényi dimension spectrum is given by(in particular, ).
The fractional moments of the mass distribution, namely,are related to bywhere is the cell volume (e.g., [23]). We restrict ourselves to moments with , which are more interesting and easier to calculate from data. From equations (10)–(13), the scaling exponent of is
Let us remark that this exponent is negative for , so that grows as shrinks, and the probability function becomes more skewed.
Conversely, as . Naturally, cumulants, which express the departure from Gaussianity, are useful in this limit and are standard in cosmology; e.g., the skewness and the excess kurtosis [3]. However, cumulants are not useful to study very skewed probability distributions, such as the mass distribution in the strong clustering regime. If we want to replace with some moment that vanishes at homogeneity, we may consider the central absolute moments, defined asor consider just , because both the quantities generalize the variance to (but are not equal).
Actually, there is more information about the multifractal properties of a mass distribution in its -moments with noninteger and close to one than in its -moments with integer . This is because the mass is concentrated on the singularities that fulfill , where is the local dimension corresponding to [33, 34]. However, (at any scale) and is close to one for close to one. Therefore, we have to take not too close to one to have a measurable change with scale. Let us see what happens in one example.
Taking and employing the sample with absolute magnitudes , for example, we find that both and the corresponding central absolute moment do have scaling ranges that are almost as large as possible, with exponents and , respectively. However, on closer inspection, the scaling with exponent splits into one part in the strong clustering regime, with exponent , and another in the quasi-homogeneous regime, with exponent . Both these exponents can be explained. The first value, using equation (14) for the exponent of , gives , that is to say, it corresponds to the vanishing dimension of a set of isolated points. The exponent should give, according to equation (14), (now hardly compatible with a vanishing dimension), but it is not sensible to calculate a fractal dimension only in the quasi-homogeneous regime. Nevertheless, we can understand the value −0.57 of the exponent as follows.
In the quasi-homogeneous regime, any mass distribution approaches normality. This is realized, in particular, by the lognomal model, whose moments are easily calculated [32], giving
For small , , and we obtain the ratio
Indeed, the exponent of is close to , which is the exponent of (equal to , with from Table 2). Moreover, if we compute the ratio along the quasi-homogeneous interval, then we find that it goes from 0.33 to 0.38, in progressively better accord with equation (17), which gives 0.375.
Regarding the scaling ofwith exponent , we have no interpretation, but it is evident that it cannot correspond to a sensible value of : to see this, let us use expression (14), recalling that is a nonincreasing function of and .
Since or the central absolute moment have no useful scaling in the quasi-homogeneous regime, we may try to find the scaling of just in the strong clustering regime. Here and henceforth, we employ the sample with absolute magnitudes . A fit in the moderately strong clustering regime yields the exponent and hence . This value is not very precise but is in accord with the also quite imprecise value deduced from the multifractal spectrum found by Gaite [15].
The dimension of the mass concentrate, , cannot be calculated as the with , because . However, the limit is easily taken in equation (11), using equations (9) and (10), and it yields the so-called entropy dimension [33]With a linear fit of the numerator of this formula versus , again in the moderately strong clustering regime, we obtain . It is somewhat smaller than the value found by Gaite [15], but it is compatible with it.
Statistical moments with integer give us little information. Of course, low-order cumulants are useful in the quasi-homogeneous regime. Moreover, we observe that the skewness has better scaling behavior than, for example, or the corresponding central absolute moment. The scaling of higher-order cumulants might play a role. Nevertheless, the relevant information about the fractal regime is provided by fractional low-order moments.
Although the above quoted results refer to the sample with , we have conducted an exploration of general scaling properties, for other samples and values of , finding that the slight rise of scaling exponent with luminosity already seen for in Section 3.1 holds in general. This general rule confirms the expected multifractal behavior, namely, the smaller values of for more luminous galaxy subsamples (Section 3). At any rate, the change of dimension is quite small for the subsamples with , which concentrate most of the stellar mass. This quasi-uniformity of scaling properties, which reflects the phenomenon of concentration of mass in multifractal geometry, justify us to consider all the galaxies at once for some purposes and, in particular, to study the angular projection of the full apparent-magnitude sample, without concern about mixing galaxies in a broad range of luminosities.
5. Fractal Projections
The study of fractal projections has tradition in mathematics [34, 39]. In cosmology, the properties of the angular projection of a fractal set have been considered in regard to the possibility of a fractal universe with no transition to homogeneity [7, 9, 26]. Of course, the study of the angular projection of the distribution of galaxies, that is to say, of angular galaxy surveys, is old and predates the study of redshift surveys [4, 5]. In fact, Peebles [3] uses basic properties of the angular distribution of galaxies as an argument against a fractal universe with no transition to homogeneity. The argument is also applicable to the stellar mass distribution.
To study the angular distribution of stellar mass, we have to consider the generalization of the theory of projection of fractal sets to the theory of projection of fractal measures, which has been developed more recently [39]. This is a necessary step because the large-scale mass distribution has to be treated not as a set but as a measure, the measure being the mass. A useful concept is the support of a mass distribution, namely, the smallest closed set that contains all the mass. The large-scale mass distribution appears to be a multifractal mass distribution of nonlacunar type, that is to say, with support in the full space [15]. This property amounts to the absence of totally empty cosmic voids. These concepts are explained in detail by Gaite [23].
Regarding lacunarity, let us recall that Durrer et al. [26] argued that the angular projection of a fractal set can have a vanishing lacunarity, by noting that the projection of a fractal with dimension 2 or higher is nonfractal and appealing to galaxy properties (apparent sizes and opacities). If the three-dimensional stellar mass distribution is already nonlacunar, then the issue is no longer meaningful. Furthermore, the presence of a small lacunarity would also give rise to a nonlacunar projection, as a consequence of the almost certain fact that the fractal dimension of the support of the mass distribution in space is larger than two [15].
The study of multifractal projections is not reduced to the behaviour of lacunarity. The main question is how to characterize the behavior of dimensions under a projection. This question has been partially answered by Hunt and Kaloshin [40], in terms of the Rényi dimension spectrum . The answer is the natural generalization of the standard and intuitive result obtained for fractal sets [34, 39]: the dimension of the projection of a fractal set onto a plane (or a smooth surface), in particular, equals the dimension of the set if it is smaller than two and it is two otherwise (this statement has to be qualified for nonrandom fractals, which can have special projections along some directions). This result is still valid for the dimension spectrum of a mass distribution, with the restriction that [40]. The restriction is due to technical reasons and may not apply to every type of mass distribution. Anyway, we are mostly interested in that interval and, especially, in the correlation dimension . Since we expect that is in the range 1–1.5 for the stellar mass distribution (Sections 2 and 3.1), it has to be preserved under angular projections.
Here, we should notice that some authors find that (Table 1 of [25]). However, large values of often come from power-law fits of on scale ranges that extend too far and include a quasi-homogeneous range, thus being biased towards , as explained in Section 2. Moreover, all those values of do not refer to the stellar mass distribution but to the galaxy number distribution.
To apply coarse multifractal analysis to an angular projection, we partition the projected region in cells of solid angle and replace in (9) the length by . Now, the formula that is analogous to equation (13) and gives the moment iswhere is assumed that . The expected scaling exponent of is, in analogy with equation (14):
In particular, is expected to scale with exponent . Naturally, the exponent is only valid provided that , because 2 is the maximal fractal dimension of a two-dimensional projection.
Let us now recall and generalize some standard results about the angular two-point correlation function, in the light of the theory of fractal projections.
5.1. Angular Correlation of Galaxies
The reduced two-point correlation function of the angular positions in a flux-limited sample is denoted by , where is the angular distance between two points. This function can be expressed as an integral of the two-point correlation function of ordinary positions over the radial coordinates and . The integral can be simplified if is a power law, equation (1). Then, in the small angle approximation, (in radians), the integral giveswhere is a nondimensional constant that depends on and the radial selection function and is the characteristic sample depth [3]. The integration variable is . The integral over is left unevaluated to show that the present approximation fails if the integral is divergent, namely, if . In this case, the projected correlation function is dominated by pairs of points at large relative radial distances.
Of course, all the above is as applicable to the stellar mass distribution as to the galaxy number distribution. To the scaling exponent corresponds the fractal dimension , so the cases or correspond, respectively, to or , the latter being the case of nonfractal projection. If , then the projected correlation function is dominated by pairs of points at small relative radial distances, and the three-dimensional fractal structure is preserved in the projection.
It is also relevant that and the integral in equation (22), for and not too close to one, are factors of the order of unity, so the magnitude of is ruled by the quotient [3]. If the characteristic sample depth is much larger than , as is normal in deep surveys, the angular correlations are small, except at very small angles. This means that the projected fractal structure is only observable at very small , where and the mass fluctuations are large.
Notice that the fractal structure appears blurred and the cosmic web features are hardly perceptible in angular images of the galaxy distribution; for example, in the image of the Lick survey ([3]; p. 41). Images like this one show that , which is a proof of large-scale homogeneity, namely, of a small ratio . These images do not reveal the web features that are so conspicuous in slices of three-dimensional redshift surveys. In the angular projection of an ideal mathematical fractal, the relevant features must always appear on very small angles. But the galaxies have discrete nature. Although is integrated down to zero in equation (22), correlation functions are actually calculated as sums over pairs of points in a finite set. In particular, there is some pair at a minimal distance, while the density of projected far points on the angular location of that pair keeps growing as grows. It is evident that uncorrelated far points must blur the small scale features and obliterate them at some stage. In fact, the projection of uncorrelated pairs of points adds a sort of Poisson distributed component, thus making the average angular density grow without altering its fluctuations. Therefore, the variance gets uniformly depressed, with no change of the scaling behavior (if it exists).
Peebles ([3]; p. 220) shows log-log plots of for the Zwicky, Lick, and Jagellonian catalogs. In fact, only the Zwicky catalog, with limiting apparent magnitude , has a range of small with . However, the absolute value of the slope grows in that range and tends to two, the value that corresponds to , that is to say, to a distribution of isolated points. Something similar happens in the plot of for successive slices of the APM catalog with ([3]; p. 221). Analyses of modern and hence deeper catalogs seldom show any . For example, Wang et al. [31] only show a very small range of with .
In conclusion, any strong angular correlations are likely to be buried in a homogeneous background, as indeed happens in our case (Section 6).
6. Angular Analysis of the Stellar Mass Distribution
We employ the same set of galaxies as in Section 3, in the angular rectangle defined therein. The lattice of cells and coarse-graining system are the same as in Section 3 (from [15]), without the radial coordinate, that is to say, a basic lattice and a sequence of binary subdivisions. To cover larger solid angles, we add a coarser lattice, namely, a lattice. All the cells have an aspect ratio close to one, which is convenient.
According to equation (22), the average of over a small cell of solid angle is proportional to . Therefore, the expected scaling of the cell mass variance iswhere is the solid angle at which the projected distribution approaches homogeneity and can be expressed in terms of , and . If , then equation (23) is a particular case of the fractal scaling of -moments (20), with exponent (21), in this case, .
In Figure 3 is the log-log plot of ( is normalized to the total solid angle from now onwards, unless the unit “steradian” appears explicitly). The exponent of the power-law fit of in turns out to be somewhat large in absolute value, namely, , giving . In the fitted range, , while grows for small and tends to three, due to the effect of discreteness.
Therefore, in the range of where the projected stellar mass distribution can be considered as a continuous mass distribution, it is quite uniform, with small fluctuations that are power-law correlated. The mass fluctuations are reduced by the projection, as explained in Section 5.1. Presumably, we can have a better representation of the correlation function for small by suppressing the shot noise.
6.1. Shot Noise Suppression
The Poissonian fluctuations of an uncorrelated distribution of points in some volume are characterized by a number density variance equal to , where is the mean number of points. This term is also present when the points are correlated. In the case of a distribution of particles with different masses, such as galaxies, and assuming that masses are statistically independent of positions, the shot noise term becomes ([3]; 509–510):
In the coarse multifractal analysis of a sample of particles, when the cell size is so small that no cell contains more than one particle and therefore the sums over cells and over particles coincide, the calculation of the moment by equation (9) obtains the shot noise term:where the index of the sums runs over the set of particles. If we start with a cell size such that no cell contains more than one particle and let the size grow, then, at some point, some cells contain more than one particle and grows, because the restricted sum over these cells is larger than the sum over the particles in the cells [for a cell with two particles, ]. Therefore, the value given by (25) is the minimum value of , and while stays constant, , according to equations (10) and (11). At some larger scale, is definitely growing, and there is a crossover to a scaling with .
Peebles [3] says about the shot noise term that “in most applications to be discussed here, this shot noise term is subdominant and will be dropped” and advises that, where the shot noise term is appreciable, it is to be subtracted. We subtract the shot noise term (25) from and see the effect in Figure 4, to be compared with Figure 3. It is now possible to extend the power-law range to , with exponent (a long scaling range and a small error). This exponent gives , which agrees with the values found in Section 3. It also agrees with the results of Li and White [16] and of Wang et al. [31] (the latter for galaxy positions). Their procedure for calculating the correlation function automatically suppresses the shot noise.
6.2. Finding the Scale
The scale can be found from the projected distribution by first finding the angular scale of homogeneity in equation (23). According to equation (23), is the solution of the equation . The scale of homogeneity of the three-dimensional distribution is found by expressing the mass variance as the average of , which is in turn given by equation (22). The equation that relates , , , and follows from equations (22) and (23) and is written aswhere is the integral in equation (22), and the integral over and extends over a cell.
Given , the evaluation of and the integral over and are straightforward numerical computations (for this integral, we take into account the small size of the relevant cells to ignore the spherical geometry and carry out the computation in the plane). To calculate the constant , which is an integral that contains the luminosity function, we can use the Schechter model [3]. This model has two parameters, namely, the characteristic galaxy luminosity and the slope , but only depends on the latter. We take , according to Zandivarez and Martínez [41] and references therein.
To solve for we have to fix and . Taking the value of found in Section 6.1, we obtain that steradians. With , we compute that and that the integral over and is 2.33, while (the uncertainties in and the integral are negligible, in comparison). Employing equation (26) (and neglecting the uncertainty in ), we obtainwhere the uncertainty is small enough to assume that it is normally distributed.
To estimate the characteristic sample depth , which is proportional to the square root of , we take , the absolute magnitude corresponding to , to be
This value is based on the analysis of Zandivarez and Martínez [41] and is in the interval of magnitudes of the most representative galaxies by mass density (Table 1). We obtainwhere the uncertainty is again small enough to assume that it is normally distributed. Hence,
This is the clustering length of equation (1), in contrast with the (luminosity-dependent) homogeneity scale of Section 3, which is tailored to the sample geometry and coarse-graining method. The value of agrees with the value of Li and White [16].
7. Summary and Discussion
The standard scaling law of galaxy clustering is the power-law correlation function of galaxy positions, with canonical exponent and clustering length h−1 Mpc, which dates back over 50 years, since the early analyses of the galaxy-galaxy angular correlation function . We have shown that there are two important ingredients that are missing in most analyses of the galaxy distribution, namely, galaxy masses and relevant statistics other than the variance. Indeed, more important than the distribution of galaxy positions is the distribution of stellar mass, which is a proxy for the distribution of baryonic matter. Besides, statistics appropriate to describe strong galaxy clustering are the fractional moments, seldom employed. We have shown that a complete study of scaling laws of the distribution of stellar mass reveals novelties that are worth considering.
Our study is based on the theory of multifractal geometry. This theory assumes scaling laws for the growth of mass fluctuations on small scales, which define dimensions. The correlation dimension , which is only one of many dimensions , derives from the scaling of the second-order statistical moment of the mass probability function. This moment also yields the scale of homogeneity. It is usually defined as the clustering length , but a coarse-graining analysis obtains a volume (which depends somewhat on the method). In Section 2, we have seen that but that the mass fluctuations at the scale are still considerable. A state of quasi-homogeneity is reached only for length scales that are almost an order of magnitude larger.
Conversely, the strong clustering regime, where fractal geometry applies, takes place for length scales that are almost an order of magnitude smaller, close to 1 Mpc. Unfortunately, there is not even one galaxy, on average, in a volume of such diameter; especially, if the volume is extracted from a volume-limited sample of galaxies. This leads to strong discreteness effects, which appear as errors in the calculation of and . We have indeed shown that a strict definition of fractal scaling leads to considerable errors, namely, relative errors of 20% in and of 30% in .
Nevertheless, a definition of scaling that encompasses the quasi-homogeneous regime leads to more precise results. The cell mass variance decreases in this regime yet the scaling continues, allowing us to calculate a precise scaling exponent . Although is somewhat dependent on the volume-limited sample that we consider, its growth with luminosity is expected in a multifractal stellar mass distribution. Indeed, galaxy luminosity and stellar mass are strongly correlated and the stellar mass of a galaxy is correlated with the local dimension of the stellar mass distribution at its position. However, most stellar mass is concentrated in a restricted range of dimensions, namely, –1.3 (), which is representative of the full stellar mass distribution. The homogeneity scale also grows with luminosity, being about h−1 Mpc for the set of galaxies that concentrates most of the mass.
It is to be remarked that the concentration of stellar mass in a range of galaxy luminosities that contains a small fraction of the galaxy number density makes the analysis of the stellar mass distribution a subject in its own right, quite different from the usual analysis of the galaxy number distribution. Nevertheless, the values of and from both analyses are compatible (we recall our calculation of below). In this regard, we agree with the results of Li and White [16].
The mass variance is not a sufficient statistic in the fractal regime and must be complemented with fractional moments , for , especially, . The study of fractional moments confirms the multifractal nature of the stellar mass distribution but shows no way to redefine -moments (with ) that extends the multifractal scaling to the quasi-homogeneous regime. In this regime, fractional moments simply adopt the values that correspond to a Gaussian distribution, as shown using the lognormal model, which interpolates between strong clustering and homogeneity. The long scaling range of the mass variance, which goes across the transition to homogeneity into the quasi-homogeneous regime, is surely a consequence of its origin in the linear evolution of Gaussian initial conditions.
Even the comparatively long scaling range of the mass variance is limited to less than two orders of magnitude in length (), as obtained from our volume-limited samples. This limitation is due to discreteness errors. In an attempt to reduce the discreteness errors, we study the angular distribution of the full flux-limited sample, which avails the galaxies in our angular rectangle. However, in the calculation of , many pairs of galaxies in a given angular region correspond to galaxies at very large distances, which are uncorrelated. So to speak, we enhance both the signal and the noise. After suppressing the shot noise, the range of scaling of the angular stellar mass distribution grows to sr, with exponent . This range of , almost five orders of magnitude, is equivalent to a linear range of half of that, namely, two and half orders of magnitude.
From the angular scale of transition to homogeneity , it is possible to calculate by expressing in terms of and several calculable magnitudes. In particular, this calculation involves the characteristic parameters of the galaxy luminosity function, which are somewhat uncertain. We obtain a reasonable value of , in the range h−1 Mpc.
To recapitulate and conclude, we have carried out a multifractal analysis of the stellar mass distribution, combining a set of volume-limited samples with the angular distribution. The methods developed here constitute an appealing alternative to other methods, such as methods that only use the galaxy positions or the two-point correlation function. Given that the important cosmological parameter is theoretically defined in terms of the fluctuations of the full mass distribution, our methods can be relevant in the calculation of precision values of .
Data Availability
The galaxy survey data used to support the findings of this study have been downloaded from the New York University Value-Added Galaxy Catalog (http://sdss.physics.nyu.edu/vagc) and are available from the corresponding author upon request.
Disclosure
This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.
Conflicts of Interest
The author declares that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
The author thanks C.A. Chacón-Cardona for the SDSS-DR7 file.