# Scaling Laws in the Stellar Mass Distribution and the Transition to Homogeneity

Authors: José Gaite
Journal: Advances in Astronomy
 (2021)

---

## 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γ=1.71−1.82, 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 γ=1.83 and a clustering length r0=5.8−7.0h−1Mpc. 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 10h−1Mpc.


---

## Body

## 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 8 h−1 Mpc, named σ8. 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 R, σR, defines a scale such that σR=1. 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 distancer0 such that ξr0=1, where ξr is the reduced galaxy-galaxy correlation function. Actually, this function is well approximated by the power-law:(1)ξr=r0rγ,for distances not much larger than r0 [3]. The scale such that the rms dispersion of galaxy counts is equal to one is related to r0 by a γ-dependent factor [3]. Historically, the power-law form of ξr was deduced from the angular positions of galaxies, yielding γ=1.8 and r0=4.7 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 wθ is used, in combination with other data, to determine precision values of σ8 [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 r≪r0, for any magnitude of r0.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 thanr0 [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 σ8 are assigned a relative error of about one percent [2]. This precision is surprisingly good, in comparison, given that σ8 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 galaxypositions. 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 wrp, on scales rp<30 h−1Mpc (rp is the separation perpendicular to the line of sight [17]). Li and White find that wrp is very well represented, over the range 10 h−1 kpc <rp<10 h−1 Mpc, by a power law that corresponds to equation (1) with γ=1.84 and r0=6.1 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 produceswrp is performed on the correlation function ξr 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 Section2 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 r0 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 ξr is not small, that is to say, for r not much larger than r0 (but it could also be valid for r≫r0). The coarse-grained mass fluctuation δMR/MR, namely, the scale-dependent rms fluctuation of the stellar mass in a cell of linear dimension R, 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,(2)δMRMR2=Cγr0Rγ,where the form of Cγ depends on the shape of the cell [3].Regardless of any scaling property,δMR/MR must decrease with R and the scale of homogeneity can be defined just by the condition δMR/MR=1. 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 δMR/MR2=0.1, 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 δMR/MR2=0.1 is close to the one defined by δMR/MR=1 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 radiusR, Cγ in equation (2) is given by Peebles [3]. The scale R such that δMR/MR equals a given number is a definite function of γ. For δMR/MR=1, R is only a little larger than r0, whereas, for δMR/MR2=0.1, it is between 4.7r0 and 12r0 (the larger, the smaller γ is). In conclusion, the scale such that δMR/MR=1 is roughly equivalent to r0, but the scale where the probability distribution is approximately Gaussian can be several times larger. This observation reveals that the values r0=5.4 h−1 Mpc and γ=1.8 are compatible with the presence of relatively large structures, for example, cosmic voids of size ≃30 h−1 Mpc [3]. Actually, even larger structures can be observed, provided that ξr does not fall too rapidly for r≫r0 [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 spectrumfα or, alternatively, of the Rényi dimension spectrum Dq [33, 34]. The multifractal spectrum fα 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 fα and Dq are related by a Legendre transform in terms of the variables α and q. The Rényi dimension Dq expresses the scaling of the statistical q-moment of the coarse-grained mass distribution and is, in general, a nonincreasing function of q. The second-order moment has an intuitive meaning, and the corresponding Rényi dimension is D2=3−γ. General q-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 toasymptotic 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 q-moments Mq 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 r from an occupied point [7, 9, 20]. It can be expressed as(3)Γr=ϱrϱ0ϱ=ϱξr+1.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Γr of the galaxy distribution (or of its integral in the ball of radius r). 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 ξr. If ϱ≠0, then the conditional density Γr is conceptually useful but can be expressed in terms of ξr, by equation (3), and it scales as ξr+1. Naturally, the condition for fractal scaling is ξ≫1, which makes ξ and ξ+1 equivalent. However, the condition is never fulfilled in a sufficiently long range of r. Therefore, the values of r0 and γ obtained from power-law fits of either function differ, namely, r0 is larger when calculated from ξ+1. 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 v. The second moment of the coarse-grained density is(4)μ2v=ϱv2ϱv2,and the density variance (mass variance) is(5)μ2v−1=δϱv2ϱv2=δMv2Mv2.The scaling law (1) is equivalent to the scaling law for the mass variance (2). Likewise, a scaling law for ξr+1 is equivalent to a scaling law for μ2v. An example of both the scaling laws for μ2v and μ2v−1 is displayed in Figure 1, which we now explain.Figure 1
Fits to scaling for sample VLS1: momentsμ2v and μ2v−1 (v in Mpc3/h3), with fits to the scaling v/v0−γ/3, represented as dashed lines. The crossing of each dashed line with the horizontal gray line gives v0.Figure1 refers to a volume-limited sample of 1765 galaxies in the redshift range 0.003,0.013, 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 D2. The solid blue lines in Figure 1 are the graphs of μ2 and μ2−1. The fit to the power-law(6)μ2=vv0−γ/3,in the interval v∈3,190Mpc3/h3 yields v0=1960±180Mpc3/h3 and γ=1.57±0.01, that is to say, D2=3−γ=1.43±0.01 (the power-law fit is a least-squares linear fit of logμ2 versus logv). An analogous fit to the power-law,(7)μ2−1=vv0−γ/3,yields v0=850±300Mpc3/h3 and D2=1.20±0.06 (γ=1.80±0.06). Both fits are represented by dashed red lines in Figure 1. The latter fit, with γ=1.80 and v01/3=9.5 Mpc/h, basically agrees with the canonical values of γ and r0 (from the reduced two-point correlation function of galaxy positions). However, both fits are questionable for calculating D2, because the fractal regime demands μ2−1≈μ2≫1, but the fitted range extends beyond the small scales where this condition holds. Therefore, some intermediate values of v0 and γ should be more appropriate. These values are uncertain because we have too small a range of scales with both μ2≫1 and negligible discreteness effects. These effects make D2 shrink towards D2=0, as the steeper left ends of the graphs reflect.The uncertainties of 20% inD2 and 30% in v01/3 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 Dq with q<0 are affected. Therefore, the evaluation of D2 should not be affected. To confirm it, we calculate now the uncertainty of μ2v 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 μ2 for each value of v, and we compute for each v the standard deviation of those ten values. The result is that the relative error of μ2 is always smaller than 4.2% (the corresponding error bars in Figure 1 would hardly be visible).Figure 2
Varianceμ2v−1 of the stellar mass distribution (v in Mpc3/h3), for three volume-limited samples: Mr∈−19,−20 (a), Mr∈−20,−21 (b), and Mr∈−21,−22 (c). The linear fits yield the values of γ and v01/3 in Table 2.
(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-z samples are also less useful to directly calculate D2 as above. The higher z 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 μ2−1≈μ2≫1.However, we may consider the scaling of justμ2v−1, in the longest possible interval, without demanding that μ2v−1>1. For example, one can dismiss the upper graph in Figure 1, corresponding to μ2, and think of extending the scaling of μ2−1 below the line μ2−1=1. Likewise, equation (1) usually extends to ξ<1 (weak correlations). Besides, we remark that the frequently employed projected correlation function wrp (e.g., the one in [16] or [18]) is a dimensional function that renders obscure the strength of correlations. The range with ξ<1 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 ξ≫1 to the large scales where ξ≪1, 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 Mr, not to be confused with our notation of M 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 12.5<mr<17.77. The angular ranges are defined in terms of equal-area coordinates sl=sinλ and f=η+0.567 radians, where the angles λ and η constitute the SDSS coordinate system. The available ranges of our coordinates are −0.760<sl<0.793 and −0.0227<f<1.200 (radians), providing a total solid angle Ω=1.554⋅1.222=1.899 steradians. We further restrict the sample to z<0.34, obtaining 529655 galaxies.The main characteristics of our volume-limited samples are reported in Table1. 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 Mr∈−21,−22, with N=101350. 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 Mr∈−21,−22 and number density 0.00159 has a volume per galaxy of 629 Mpc3/h3, that is to say, the equivalent length is 8.6 Mpc/h, as large as the expected value of r0.Table 1
Characteristics of the volume-limited samples: absolute magnitude, redshift, number of galaxies, number density (h3Mpc−3), fraction of the total mass density, and galaxy mass (M⊙).
MrzNnDens.Mass−14,−150.0011−0.00735540.085500.0092.7×104,1.4×108−15,−160.0017−0.01158850.034720.0104.1×105,4.0×108−16,−170.0026−0.018117080.017170.0129.4×106,9.2×108−17,−180.0042−0.028466250.017280.0391.7×106,2.8×109−18,−190.0066−0.0444174160.012020.0891.2×107,6.1×109−19,−200.0105−0.0689462040.008700.2032.1×107,2.3×1010−20,−210.0165−0.1058978470.005220.3331.7×106,6.6×1010−21,−220.0260−0.16131013500.001590.2471.2×107,1.7×1011−22,−230.0406−0.2444306800.000150.0567.5×108,4.4×1011−23,−240.0632−0.359714602.4×10−60.0027.7×109,7.9×1011As regards number density, the best sample is the first one, withMr∈−14,−15 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 Mr∈−14,−17, in spite of their abundance.The last variable in Table1 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, α=−logm/lnl, ([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 Dq 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 volumev, 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 Mr∈−19,−22, 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μ2v−1 appear in Table 2 and Figure 2. To wit, in Figure 2 are the log-log plots of μ2v−1 for the three samples with Mr∈−19,−22, 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 v, which crosses over, for small v, to the trivial scaling corresponding to isolated points (D2=0). 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.Table 2
Results for scaling in the five volume-limited samples with higher mass density (v01/3 in Mpc/h).
Mrγv01/3−18,−191.52±0.0713.5±1.6−19,−201.71±0.0412.4±0.8−20,−211.616±0.01314.2±0.3−21,−221.818±0.00716.9±0.2−22,−231.93±0.0229.4±0.6Regarding the values of the scaling exponentγ in Table 2, we can observe a definite growing trend, except in the sample with Mr∈−20,−21 (an anomaly in this range is also noted by [18]). The length scale v01/3, 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, withD2=3−γ∈1,1.5, 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 Mr∈−19,−22, which concentrate most of the mass, correspond to a narrow interval of the scaling exponent, namely, γ∈1.75,1.83, and hence to a narrow interval of D2. 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.


## 3.1. Scaling of Cell Mass Variances

Here, we examine if cell mass variances are power-law functions of the cell volumev, 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 Mr∈−19,−22, 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μ2v−1 appear in Table 2 and Figure 2. To wit, in Figure 2 are the log-log plots of μ2v−1 for the three samples with Mr∈−19,−22, 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 v, which crosses over, for small v, to the trivial scaling corresponding to isolated points (D2=0). 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.Table 2
Results for scaling in the five volume-limited samples with higher mass density (v01/3 in Mpc/h).
Mrγv01/3−18,−191.52±0.0713.5±1.6−19,−201.71±0.0412.4±0.8−20,−211.616±0.01314.2±0.3−21,−221.818±0.00716.9±0.2−22,−231.93±0.0229.4±0.6Regarding the values of the scaling exponentγ in Table 2, we can observe a definite growing trend, except in the sample with Mr∈−20,−21 (an anomaly in this range is also noted by [18]). The length scale v01/3, 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, withD2=3−γ∈1,1.5, 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 Mr∈−19,−22, which concentrate most of the mass, correspond to a narrow interval of the scaling exponent, namely, γ∈1.75,1.83, and hence to a narrow interval of D2. 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μ2v and, specifically, on the variance μ2v−1. 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μk≥μ2k−1,k∈ℕ ([36]; p. 55), we deduce a lower bound to the skewness:(8)δMv3Mv3≥μ2−1μ2−2.Hence,μ2≫1 implies that the probability distribution function is very skewed. The lognormal distribution is an example that goes from nearly normal for μ2≳1 to very skewed for μ2≫1. 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 fractionalq-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 l. Using this partition, fractional statistical moments are defined as(9)ℳql=∑imiMq,q∈ℝ,where the index i runs over the set of nonempty cells, mi is the mass in the cell i, and M=∑imi is the total mass. Then, multifractal behavior is given, for l⟶0, by(10)ℳql∼lτq,and the Rényi dimension spectrum is given by(11)Dq=τqq−1,(in particular, D2=τ2).The fractional moments of the mass distribution, namely,(12)μqv=ϱvqϱvq=MvqMvq,are related to ℳq by(13)μqv=ℳqvvq−1,where v is the cell volume (e.g., [23]). We restrict ourselves to moments with q>0, which are more interesting and easier to calculate from data. From equations (10)–(13), the scaling exponent of μqv is(14)τq3−q−1=q−1Dq−33.Let us remark that this exponent is negative forq>1, so that μqv grows as v shrinks, and the probability function becomes more skewed.Conversely,μqv⟶1 as v⟶∞. 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 μq with some moment that vanishes at homogeneity, we may consider the central absolute moments, defined as(15)δMvqMvq,q∈ℝ,or consider just μqv−1, because both the quantities generalize the variance to q≠2 (but are not equal).Actually, there is more information about the multifractal properties of a mass distribution in itsq-moments with q noninteger and close to one than in its q-moments with integer q>1. This is because the mass is concentrated on the singularities that fulfill D1=α1=fα1, where α1 is the local dimension corresponding to q=1 [33, 34]. However, μ1=1 (at any scale) and μq is close to one for q close to one. Therefore, we have to take q not too close to one to have a measurable change with scale. Let us see what happens in one example.Takingq=1.5 and employing the sample with absolute magnitudes Mr∈−21,−22, for example, we find that both μ1.5v−1 and the corresponding central absolute moment do have scaling ranges that are almost as large as possible, with exponents −0.54 and −0.41, respectively. However, on closer inspection, the scaling with exponent −0.54 splits into one part in the strong clustering regime, with exponent −0.51, and another in the quasi-homogeneous regime, with exponent −0.57. Both these exponents can be explained. The first value, using equation (14) for the exponent of μ1.5v≈μ1.5v−1, gives D1.5=−0.06, that is to say, it corresponds to the vanishing dimension of a set of isolated points. The exponent −0.57 should give, according to equation (14), D1.5=−0.42 (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(16)μq=eqq−1σ2/2.For smallσ2, μq≳1, and we obtain the ratio(17)μq−1μ2−1≈qq−12.Indeed, the exponent−0.57 of μ1.5v−1 is close to −0.61, which is the exponent of μ2v−1 (equal to −γ/3, with γ from Table 2). Moreover, if we compute the ratio μ1.5v−1/μ2v−1 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 of(18)δMv1.5Mv1.5,with exponent −0.41, we have no interpretation, but it is evident that it cannot correspond to a sensible value of D1.5: to see this, let us use expression (14), recalling that Dq is a nonincreasing function of q and D2>1.Sinceμ1.5v−1 or the central absolute moment have no useful scaling in the quasi-homogeneous regime, we may try to find the scaling of just μ1.5v in the strong clustering regime. Here and henceforth, we employ the sample with absolute magnitudes Mr∈−21,−22. A fit in the moderately strong clustering regime yields the exponent 0.19±0.04 and hence D1.5=1.8±0.2. This value is not very precise but is in accord with the also quite imprecise value deduced from the multifractal spectrum fα found by Gaite [15].The dimension of the mass concentrate,D1, cannot be calculated as the Dq with q≠1, because ℳ1=μ1=1. However, the limit q⟶1 is easily taken in equation (11), using equations (9) and (10), and it yields the so-called entropy dimension [33](19)D1=liml⟶0∑imi/Mlogmi/Mlogl.With a linear fit of the numerator of this formula versus logv, again in the moderately strong clustering regime, we obtain D1=1.9±0.2. It is somewhat smaller than the value found by Gaite [15], but it is compatible with it.Statistical moments with integerq>2 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, μ3v−1 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 withMr∈−21,−22, we have conducted an exploration of general scaling properties, for other samples and values of q, finding that the slight rise of scaling exponent with luminosity already seen for q=2 in Section 3.1 holds in general. This general rule confirms the expected multifractal behavior, namely, the smaller values of Dq for more luminous galaxy subsamples (Section 3). At any rate, the change of dimension is quite small for the subsamples with Mr∈−19,−22, 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 fractalmeasures, 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 Dq. 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 Dq of a mass distribution, with the restriction that 1<q≤2 [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 D2. Since we expect that D2 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 thatD2≥2 (Table 1 of [25]). However, large values of D2 often come from power-law fits of ξr+1 on scale ranges that extend too far and include a quasi-homogeneous range, thus being biased towards D2=3, as explained in Section 2. Moreover, all those values of D2 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 l by Ω1/2. Now, the formula that is analogous to equation (13) and gives the moment μqΩ is(20)μqΩ=ϱΩqϱΩq=ℳqΩΩq−1,where is assumed that q>0. The expected scaling exponent of μqΩ is, in analogy with equation (14):(21)τq2−q−1=q−1Dq2−1.In particular,μ2Ω is expected to scale with exponent D2/2−1. Naturally, the exponent q−1Dq/2−1 is only valid provided that Dq<2, 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 bywθ12, where θ12 is the angular distance between two points. This function can be expressed as an integral of the two-point correlation function of ordinary positions ξr12 over the radial coordinates r1 and r2. The integral can be simplified if ξr is a power law, equation (1). Then, in the small angle approximation, θ≪1 (in radians), the integral gives(22)wθ=Kθ1−γr0d∗γ∫−∞∞dx1+x2γ/2,where K is a nondimensional constant that depends on γ and the radial selection function and d∗ is the characteristic sample depth [3]. The integration variable is x=r2−r1/r1θ. The integral over x is left unevaluated to show that the present approximation fails if the integral is divergent, namely, if γ≤1. 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 D2=3−γ, so the cases γ>1 or γ≤1 correspond, respectively, to D2<2 or D2≥2, the latter being the case of nonfractal projection. If γ>1, 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 thatK and the integral in equation (22), for γ>1 and not too close to one, are factors of the order of unity, so the magnitude of wθ is ruled by the quotient r0/d∗ [3]. If the characteristic sample depth is much larger than r0, 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 wθ≫1 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 w≪1, which is a proof of large-scale homogeneity, namely, of a small ratio r0/d∗. 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 x 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 d∗ 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 δϱΩ/ϱΩ2 gets uniformly depressed, with no change of the scaling behavior (if it exists).Peebles ([3]; p. 220) shows log-log plots of wθ for the Zwicky, Lick, and Jagellonian catalogs. In fact, only the Zwicky catalog, with limiting apparent magnitude m=15, has a range of small θ with wθ>1. However, the absolute value of the slope grows in that range and tends to two, the value that corresponds to D2=0, that is to say, to a distribution of isolated points. Something similar happens in the plot of wθ for successive slices of the APM catalog with Δm=0.5 ([3]; p. 221). Analyses of modern and hence deeper catalogs seldom show any w>1. For example, Wang et al. [31] only show a very small range of θ with wθ>1.In conclusion, any strong angular correlations are likely to be buried in a homogeneous background, as indeed happens in our case (Section 6).


## 5.1. Angular Correlation of Galaxies

The reduced two-point correlation function of the angular positions in a flux-limited sample is denoted bywθ12, where θ12 is the angular distance between two points. This function can be expressed as an integral of the two-point correlation function of ordinary positions ξr12 over the radial coordinates r1 and r2. The integral can be simplified if ξr is a power law, equation (1). Then, in the small angle approximation, θ≪1 (in radians), the integral gives(22)wθ=Kθ1−γr0d∗γ∫−∞∞dx1+x2γ/2,where K is a nondimensional constant that depends on γ and the radial selection function and d∗ is the characteristic sample depth [3]. The integration variable is x=r2−r1/r1θ. The integral over x is left unevaluated to show that the present approximation fails if the integral is divergent, namely, if γ≤1. 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 D2=3−γ, so the cases γ>1 or γ≤1 correspond, respectively, to D2<2 or D2≥2, the latter being the case of nonfractal projection. If γ>1, 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 thatK and the integral in equation (22), for γ>1 and not too close to one, are factors of the order of unity, so the magnitude of wθ is ruled by the quotient r0/d∗ [3]. If the characteristic sample depth is much larger than r0, 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 wθ≫1 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 w≪1, which is a proof of large-scale homogeneity, namely, of a small ratio r0/d∗. 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 x 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 d∗ 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 δϱΩ/ϱΩ2 gets uniformly depressed, with no change of the scaling behavior (if it exists).Peebles ([3]; p. 220) shows log-log plots of wθ for the Zwicky, Lick, and Jagellonian catalogs. In fact, only the Zwicky catalog, with limiting apparent magnitude m=15, has a range of small θ with wθ>1. However, the absolute value of the slope grows in that range and tends to two, the value that corresponds to D2=0, that is to say, to a distribution of isolated points. Something similar happens in the plot of wθ for successive slices of the APM catalog with Δm=0.5 ([3]; p. 221). Analyses of modern and hence deeper catalogs seldom show any w>1. For example, Wang et al. [31] only show a very small range of θ with wθ>1.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 of529655 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 4×3 basic lattice and a sequence of binary subdivisions. To cover larger solid angles, we add a coarser lattice, namely, a 3×2 lattice. All the cells have an aspect ratio close to one, which is convenient.According to equation (22), the average of wθ over a small cell of solid angle Ω is proportional to Ω1−γ/2. Therefore, the expected scaling of the cell mass variance is(23)μ2Ω−1=δMΩ2MΩ2=ΩΩ0−γ−1/2,where Ω0 is the solid angle at which the projected distribution approaches homogeneity and can be expressed in terms of K,r0/d∗, and γ. If μ2≫1, then equation (23) is a particular case of the fractal scaling of q-moments (20), with exponent (21), in this case, D2/2−1=3−γ/2−1.In Figure3 is the log-log plot of μ2Ω−1 (Ω is normalized to the total solid angle from now onwards, unless the unit “steradian” appears explicitly). The exponent of the power-law fit of μ2Ω−1 in Ω∈8⋅10−5,8⋅10−2 turns out to be somewhat large in absolute value, namely, γ−1/2=0.486±0.008, giving γ=1.97±0.02. In the fitted range, μ2−1≪1, while μ2 grows for small Ω and γ tends to three, due to the effect of discreteness.Figure 3
Scale-dependent variance of the projected stellar mass, with a power-law fit givingγ=1.97.Figure 4
Values ofμ2Ω−1 after shot noise suppression, for the distributions of stellar mass, with linear fit (dashed line). It gives γ=1.83.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 to1/N, where N 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):(24)m2m2N.In the coarse multifractal analysis of a sample ofN 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 ℳ2 by equation (9) obtains the shot noise term:(25)∑mi2∑mi2=m2m2N,where the index i 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 ℳ2 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, m1+m22>m12+m22]. Therefore, the value given by (25) is the minimum value of ℳ2, and while ℳ2 stays constant, D2=0, according to equations (10) and (11). At some larger scale, ℳ2 is definitely growing, and there is a crossover to a scaling with D2>0.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 ℳ2Ω and see the effect in Figure 4, to be compared with Figure 3. It is now possible to extend the power-law range to Ω∈1.3⋅10−6,8⋅10−2, with exponent γ−1/2=0.415±0.007 (a long scaling range and a small error). This exponent gives γ=1.83, 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 Scaler0

The scaler0 can be found from the projected distribution by first finding the angular scale of homogeneity Ω0 in equation (23). According to equation (23), Ω0 is the solution of the equation μ2Ω0−1=1. The scale of homogeneity of the three-dimensional distribution is found by expressing the mass variance μ2Ω−1 as the average of wθ, which is in turn given by equation (22). The equation that relates Ω0, r0/d∗, K, and γ follows from equations (22) and (23) and is written as(26)Ω0γ−1/2=Ωγ−1/2∫dΩ1dΩ2Ω2wθ12=Kr0d∗γkγ∫dΩ1dΩ2Ω2θ12Ω1/21−γ,where kγ is the integral in equation (22), and the integral over Ω1 and Ω2 extends over a cell.Givenγ, the evaluation of kγ and the integral over Ω1 and Ω2 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 K, 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 L∗ and the slope α, but K only depends on the latter. We take α∈−1.2,−0.8, according to Zandivarez and Martínez [41] and references therein.To solve forr0/d∗ we have to fix γ and Ω0. Taking the value of γ found in Section 6.1, we obtain that Ω0=9.53⋅10−7 steradians. With γ=1.83, we compute that kγ=3.58 and that the integral over and Ω2 is 2.33, while K∈0.89,1.17 (the uncertainties in kγ and the integral are negligible, in comparison). Employing equation (26) (and neglecting the uncertainty in Ω0), we obtain(27)r0d∗=0.0134±0.0010,where the uncertainty is small enough to assume that it is normally distributed.To estimate the characteristic sample depthd∗, which is proportional to the square root of L∗, we take M∗, the absolute magnitude corresponding to L∗, to be(28)M∗−5log10h=−20.6±0.2.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 obtain(29)d∗=10pc×100.217.77+20.6±0.2−5log10h=470±40Mpc/h,where the uncertainty is again small enough to assume that it is normally distributed. Hence,(30)r0=6.3±0.7Mpc/h.This is the clustering length of equation (1), in contrast with the (luminosity-dependent) homogeneity scale v01/3 of Section 3, which is tailored to the sample geometry and coarse-graining method. The value of r0 agrees with the value of Li and White [16].


## 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 to1/N, where N 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):(24)m2m2N.In the coarse multifractal analysis of a sample ofN 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 ℳ2 by equation (9) obtains the shot noise term:(25)∑mi2∑mi2=m2m2N,where the index i 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 ℳ2 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, m1+m22>m12+m22]. Therefore, the value given by (25) is the minimum value of ℳ2, and while ℳ2 stays constant, D2=0, according to equations (10) and (11). At some larger scale, ℳ2 is definitely growing, and there is a crossover to a scaling with D2>0.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 ℳ2Ω and see the effect in Figure 4, to be compared with Figure 3. It is now possible to extend the power-law range to Ω∈1.3⋅10−6,8⋅10−2, with exponent γ−1/2=0.415±0.007 (a long scaling range and a small error). This exponent gives γ=1.83, 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 Scaler0

The scaler0 can be found from the projected distribution by first finding the angular scale of homogeneity Ω0 in equation (23). According to equation (23), Ω0 is the solution of the equation μ2Ω0−1=1. The scale of homogeneity of the three-dimensional distribution is found by expressing the mass variance μ2Ω−1 as the average of wθ, which is in turn given by equation (22). The equation that relates Ω0, r0/d∗, K, and γ follows from equations (22) and (23) and is written as(26)Ω0γ−1/2=Ωγ−1/2∫dΩ1dΩ2Ω2wθ12=Kr0d∗γkγ∫dΩ1dΩ2Ω2θ12Ω1/21−γ,where kγ is the integral in equation (22), and the integral over Ω1 and Ω2 extends over a cell.Givenγ, the evaluation of kγ and the integral over Ω1 and Ω2 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 K, 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 L∗ and the slope α, but K only depends on the latter. We take α∈−1.2,−0.8, according to Zandivarez and Martínez [41] and references therein.To solve forr0/d∗ we have to fix γ and Ω0. Taking the value of γ found in Section 6.1, we obtain that Ω0=9.53⋅10−7 steradians. With γ=1.83, we compute that kγ=3.58 and that the integral over and Ω2 is 2.33, while K∈0.89,1.17 (the uncertainties in kγ and the integral are negligible, in comparison). Employing equation (26) (and neglecting the uncertainty in Ω0), we obtain(27)r0d∗=0.0134±0.0010,where the uncertainty is small enough to assume that it is normally distributed.To estimate the characteristic sample depthd∗, which is proportional to the square root of L∗, we take M∗, the absolute magnitude corresponding to L∗, to be(28)M∗−5log10h=−20.6±0.2.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 obtain(29)d∗=10pc×100.217.77+20.6±0.2−5log10h=470±40Mpc/h,where the uncertainty is again small enough to assume that it is normally distributed. Hence,(30)r0=6.3±0.7Mpc/h.This is the clustering length of equation (1), in contrast with the (luminosity-dependent) homogeneity scale v01/3 of Section 3, which is tailored to the sample geometry and coarse-graining method. The value of r0 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γ=1.8 and clustering length r0≃5 h−1 Mpc, which dates back over 50 years, since the early analyses of the galaxy-galaxy angular correlation function wθ. 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 dimensionD2, which is only one of many dimensions Dq, 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 r0, but a coarse-graining analysis obtains a volume v0 (which depends somewhat on the method). In Section 2, we have seen that v01/3≃r0 but that the mass fluctuations at the scale v0 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 ofD2 and v0. We have indeed shown that a strict definition of fractal scaling leads to considerable errors, namely, relative errors of 20% in D2 and of 30% in v01/3.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, D2=1.2–1.3 (γ=1.7−1.8), which is representative of the full stellar mass distribution. The homogeneity scale v01/3 also grows with luminosity, being about 15 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 r0 from both analyses are compatible (we recall our calculation of r0 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μq, for q∈ℝ, especially, q≳1. The study of fractional moments confirms the multifractal nature of the stellar mass distribution but shows no way to redefine q-moments (with q≠2) 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 (v1/3), 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 529967 galaxies in our angular rectangle. However, in the calculation of μ2Ω−1, 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 Ω∈2.5⋅10−6,1.6⋅10−1 sr, with exponent γ=1.83. 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Ω0, it is possible to calculate r0 by expressing Ω0 in terms of r0 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 r0, in the range 5.8−7.0 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σ8 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 σ8.



---

*Source: 6680938-2021-02-08.xml*