Relativistic cosmology number densities and the luminosity function
^{1} Valongo Observatory, Federal University of Rio de Janeiro, Ladeira Pedro Antônio 43, Rio de Janeiro CEP 20080090, Brazil
email: iribarrem@astro.ufrj.br; amanda05@astro.ufrj.br
^{2} Physics Institute, Federal University of Rio de Janeiro, CxP 68532, Rio de Janeiro CEP 21941972, Brazil
email: mbr@if.ufrj.br
^{3} Vatican Observatory Research Group, Steward Observatory, University of Arizona, Tucson, AZ 85721, USA
email: wstoeger@email.arizona.edu
Received: 21 June 2011
Accepted: 15 December 2011
Aims. This paper studies the connection between the relativistic number density of galaxies down the past light cone in a FriedmannLemaîtreRobertsonWalker spacetime with nonvanishing cosmological constant and the galaxy luminosity function (LF) data. It extends the redshift range of previous results presented in Albani et al. (2007, ApJ, 657, 760), where the galaxy distribution was studied out to z = 1. Observational inhomogeneities were detected at this range. This research also searches for LF evolution in the context of the framework advanced by Ribeiro and Stoeger (2003, ApJ, 592, 1), further developing the theory linking relativistic cosmology theory and LF data.
Methods. Selection functions are obtained using the Schechter parameters and redshift parametrization of the galaxy LF obtained from an Iband selected dataset of the FORS deep field galaxy survey in the redshift range 0.5 ≤ z ≤ 5.0 for its blue bands and 0.75 ≤ z ≤ 3.0 for its red ones. Differential number counts, densities and other related observables are obtained, and then used with the calculated selection functions to study the empirical radial distribution of the galaxies in a fully relativistic framework.
Results. The redshift range of the dataset used in this work, which is up to five times larger than the one used in previous studies, shows an increased relevance of the relativistic effects of expansion when compared to the evolution of the LF at the higher redshifts. The results also agree with the preliminary ones presented in Albani et al., suggesting a powerlaw behavior of relativistic densities at high redshifts when they are defined in terms of the luminosity distance.
Key words: cosmology: observations / largescale structure of Universe / galaxies: luminosity function, mass function / cosmology: theory
© ESO, 2012
1. Introduction
The galaxy volume number density, that is, the number of galaxies enclosed in a given volume, is a very important quantity in cosmology. It gives information about the density of the massenergy in the universe and its evolution, allowing us to test various observational features of cosmological models like, for instance, the observational inhomogeneities of the number counts data down our past light cone (Albani et al. 2007). This important quantity can be observationally determined by a careful analysis of data from galaxy redshift surveys, as has been done in a systematic way by several authors who derived the galaxy luminosity function (LF) φ from such surveys. The LF of galaxies is a number density per unit of luminosity (Peacock 1999). The LF stemming from the galaxy distribution in a given dataset^{1} is commonly fitted using the profile due to Schechter (1976), (1)where l ≡ L/L_{∗}, L being the observed luminosity, L_{∗} the luminosity scale parameter, φ^{∗} a normalization parameter and α the faintend slope parameter. Various papers on the LF have used a variety of datasets in different wavelength regions: Lin et al. (1999), Fried et al. (2001), Blanton et al. (2003), Pozzetti et al. (2003), Bell et al. (2003), Norman et al. (2004), Wilmer et al. (2006), Ly et al. (2007), and Tzanavaris & Georgantopoulos (2008) obtained Schechtertype LF parameters for galaxies with redshift values out to z ~ 1.5. Poli et al. (2003), and Rudnick et al. (2003) did the same for galaxies with redshift values out to z ~ 3. More recently Bouwens et al. (2007), and Gabasch et al. (2008) obtained those parameters for galaxies with redshift values out to z ~ 6. It was found that the LF evolves with redshift in all galaxy redshift surveys.
Using the LF data to study observational features of cosmological models requires, however, some model linking relativistic cosmology number count theory with LF astronomical data and practice. One approach to this theoretical link was developed by Ribeiro and Stoeger (2003, hereafter RS03), who started from very general relativistic considerations and then specialized their theoretical results, applying them to the LF data provided by the CNOC2 redshift survey (Lin et al. 1999). This provides an example of how such a link is established, as well as facilitating various consistency tests between the data and the expected number count predictions from the assumed cosmological model. In a sequel Albani et al. (2007, hereafter A07) used the same CNOC2 dataset to actually obtain and analyze various types of observational relativistic densities, that is, the ones defined on the observer’s past light cone^{2}.
The aim of this work is to further study the observed number densities of galaxies down the observer’s past light cone. We want to focus on observational inhomogeneities of the number count data, analyzing these with differential density and integral density measures in order to separate out the effects of cosmic expansion from the effect of LF evolution. As has been discussed in detail in A07 and in Rangel Lemos & Ribeiro (2008), these measures are capable of revealing observational inhomogeneities even in spatially homogeneous cosmological models. This is because spatial homogeneity is defined on constant time hypersurfaces, and is therefore a builtin feature of Friedmann cosmologies, whereas observational homogeneity is defined by measures of massenergy densities remaining unchanged along the past light cone. We also seek to investigate if there is significant evidence of powerlaw behavior in the observed galaxy distribution at high values of redshift, since A07 detected such a behavior out to z ≈ 1. The importance of powerlaws comes from the fact that they are indicative of selforganized criticality in dynamical systems. In other words, scale invariant phenomena such as power laws are emergent features of a distribution, this being a mechanism by which complexity arises in nature through simple local interactions.
The study presented here expands the results presented in RS03 and A07, since the model specific equations obtained in RS03 were applicable only to LF data whose parameters assumed the Einsteinde Sitter (EdS) cosmology. And although the study presented in A07 was not limited to the EdS model, it was limited to the CNOC2 survey, which only probed out to z = 0.75 and extrapolated to z = 1. In addition, although both papers made some important theoretical connections between the relativistic cosmology equations and the LF data for a specific cosmology, they omitted important steps. These can be summarized as follows. The equations and methodology used in A07 for dealing with nonzero cosmological constant models are generalized, fully presented and discussed in this paper. An appendix containing a detailed algorithm on how to obtain the key results is also included. The role of the completeness function J connecting theory and observations, as advanced in RS03, is made explicit here, adding clarity to the results and presenting them in a more straightforward manner. Finally, a comparison of the geometrical versus evolutionary effects, the latter being empirically obtained from the LF, on the relativistic volume number densities is discussed, considering a broad spectral classification of late versus early type galaxies. So, this work further advances the treatment of both RS03 and A07 and presents a comprehensive and detailed discussion of how one can proceed in applying the very general, model independent, relativistic cosmology equations to specific cosmologies, thus linking any cosmological model to LF data. This is necessary if we want to use LF data to observationally test more general cosmological models (any more general models will be inhomogeneous and/or anisotropic), a task we intend to carry out in forthcoming papers. Table 1 summarizes the new contributions of this paper as compared with the results obtained in the previous ones.
Comparison with previous works.
In this paper we assume a FriedmannLemaîtreRobertsonWalker (FLRW) cosmology with nonzero cosmological constant Λ and derive model independent relativistic equations which enable us to solve numerically expressions which include LF data. The numerical scheme is presented in detail and applied to a galaxy redshift survey dataset that assumes the FLRW standard cosmological model in its volume and distance definitions: the Iband selected LF of Gabasch et al. (2004, 2006, hereafter G04, G06 respectively) for the FORS deep field (FDF). From the Schechter parameters and their redshift parametrization appearing in G04 and G06, we calculated the selection functions ψ in eight bandwidths and in equally spaced redshift intervals of 0.25 spanning the range 0.5 ≤ z ≤ 5.0, for the blue bands (G04), and 0.75 ≤ z ≤ 3.0 for the red bands (G06). Once the selection functions were calculated, we were able to obtain the observational number counts, differential number counts and differential and integral densities, as defined in RS03 and Ribeiro (2005) and further discussed in A07 and Rangel Lemos & Ribeiro (2008). These quantities were then used to study observational inhomogeneities in the relativistic radial distribution of galaxies belonging to that dataset. Therefore, the calculations of this paper go to much higher redshift than the ones presented in A07. The results reached here generally agree with theirs, also indicating that observational inhomogeneities in the relativistic galaxy distributions can arise due to geometrical, nullgeodesic effects, even in the context of spatially homogeneous universes (Ribeiro 1992, 1995, 2001, 2005; Rangel Lemos & Ribeiro 2008).
The plan of the paper is as follows. In Sect. 2 we briefly present the FLRW model, basically to set up notation and write essential results. Then we apply the general, model independent, equations of relativistic cosmology numbercount theory to this universe model in order to obtain the differential equations governing the evolution of the two basic quantities of our analysis, namely, the scale factor and the relativistic number counts. Section 3 describes in detail our numerical approach for solving these two differential equations, as well as our procedures for deriving all other quantities relevant in our analysis. In doing that we also explain how they connect to the the scale factor and the number counts. In Sect. 4 we further develop general equations connecting LF data with relativistic cosmology numbercount theory and obtain the selection functions for the FDF dataset. In Sect. 5 we study the behavior of radial statistics of the relativistic densities obtained in the context of the FLRW model with Λ ≠ 0 and out to z = 5.0. We present our conclusions in Sect. 6. The appendix presents an algorithm detailing how the methodology discussed in Sects. 2–4 can be applied to a given parametrized LF in order to obtain the observational relativistic densities discussed in Sect. 5.
2. Standard cosmology with nonzero cosmological constant
2.1. The scale factor
We begin by writing the FLRW line element as follows, (2)where the timedependent function S = S(t) is the cosmic scale factor, k is the curvature parameter (k = + 1,0, − 1) and c is the light speed. As is well known, the Einstein’s field equations with the metric corresponding to this line element yields the Friedmann equation, which, if the cosmological constant Λ is included, may be written (e.g. Roos 1994), (3)where ρ_{m} is the matter density and we have assumed the usual definition for the Hubble parameter, (4)Let us now define the vacuum energy density in terms of the cosmological constant, (5)Since the critical density at the present time is given by, (6)where H_{0} is the Hubble constant, the following relativetocritical density parameter relations hold, (7)We have used the zero index to indicate observable quantities at the present time. Notice that since Λ is a constant, then ρ_{Λ} = ρ_{Λ0}. Considering the definitions (7), we can rewrite the Friedmann Eq. (3) at the present time as follows, (8)In addition, from the law of conservation of energy applied to the zero pressure, matter dominated era, we know that, (9)which leads to, (10)since the matterdensity parameter can also be written in terms of the critical density as, (11)We can rewrite Eq. (3) as a first order ordinary differential equation in terms of the scale factor S(t), by using the results in Eqs. (4)–(11), yielding, (12)The problem we deal with in this paper requires the solution of the equation above as well as of another differential equation for the cumulative number count N (see Eq. (20) below). The simplest approach is to simultaneously solve these two differential equations by numerical means.
2.2. Relativistic number counts
Let us write the completely general, cosmological, modelindependent expression derived by Ellis (1971) for the number count of cosmological sources, which takes fully into account relativistic effects, as follows, (13)Here dN is the number of cosmological sources in a volume section at a point P down the null cone, n is the number density of radiating sources per unit of proper volume in a section of a bundle of light rays converging towards the observer and subtending a solid angle dΩ at the observer’s position, d_{A} is the area distance of this section from the observer’s viewpoint (also known as angular diameter distance, observer area distance and corrected luminosity distance), u^{a} is the observer’s 4velocity, k^{a} is the tangent vector along the light rays and y is the affine parameter distance down the light cone constituting the bundle (see RS03, Sect. 2.1, Fig. 1). The number density, that is, the number of cosmological sources per proper volume unit, can be related to the matter density ρ_{m} by means of, (14)where, ℳ_{g} is simply the average galaxy rest mass, dark matter included.
One should point out that the details of the galaxy mass function and how it evolves with the redshift are imprinted in the LF itself. As a consequence, it will be included in any selection functions stemming from the LF. For the empirical purposes of this paper this equation is correct to an order of magnitude and should be regarded only as an estimate. This means that Eq. (14) enables us to connect the theoretical relativistic quantities to the LF data and thus include the redshift evolution of the mass function empirically (more on this point at the end of Sect. 4.1 below). However, to actually extract from the LF its implicit galaxy mass function requires some assumption about the function ℳ_{g}(z). We shall deal with this problem in a forthcoming paper.
So, for the empirical approach of this paper it is enough to assume a constant average galaxy rest mass with the working value of ℳ_{g} ≈ 10^{11} M_{⊙} based on the estimate by Sparke & Gallagher (2000). Hence, if we use Eqs. (10) and (11) we can rewrite Eq. (14) as, (15)The past radial null geodesic in the geometry given by metric (2) may be written (16)Since both coordinates t and r are functions of the affine parameter y along the past null cone, we have that (17)Assuming now that both source and observer are comoving, then u^{a} = cδ_{0}^{a} and the following results hold, (18)In addition, remembering that the area distance d_{A} is defined by means of a relation between the intrinsically measured crosssectional area element dσ of the source and the observed solid angle dΩ_{0} (Ellis 1971, 2007; Plebański & Krasiński 2006), we have that (19)We note that for this metric the area distance, a strictly observational distance definition, equals the proper distance d_{Pr}, a relativistic one. Considering Eq. (8) and substituting Eqs. (15), (18) and (19) into Eq. (13), and remembering that here dΩ_{0} = 4π, we are then able to write the number of sources along the past light cone in terms of the radial coordinate r and in the FLRW model defined in Sect. 2.1. This quantity is given by the following expression, (20)
3. Numerical problem
The differential equations for the scale factor S and the cumulative number count N, Eqs. (12) and (20) respectively, provide the basic quantities necessary here, since the expressions which actually connect the relativistic theory to LF data are built upon them. However, before these two differential equations can be solved numerically, some further algebraic manipulations are necessary. Next we shall describe these steps and how other quantities relevant to our analysis are written in terms of these two functions.
The numerical problem can be summarized as follows. We take the radial coordinate r as the independent variable in order to numerically obtain the functions N(r) and S(r) along the past light cone. Then all other quantities of interest are written in terms of these two functions, meaning that numerical results for N(r) and S(r) allow us to straightforwardly obtain z(r) and dN/dz(r), as well as various cosmological distances, their derivatives, observational volumes and number densities. Thus, all quantities used in our analysis end up being written in terms of the radial coordinate r.
3.1. Scale factor
A simple and straightforward approach to solving the differential equations for the scale factor S and the cumulative number count N is to take the radial coordinate r as the independent variable and simultaneously integrate them by numerical means. To do so we shall proceed as follows.
Equation (12) has its time coordinate implicitly defined in the scale factor, therefore we need to rewrite that differential equation in such a way that the independent variable becomes explicit. Considering Eq. (8), the null geodesic (16) becomes, (21)Along the past light cone, we can write the scale factor in terms of the radial coordinate as (22)Thus, we are able to rewrite Eq. (12) in terms of the radial coordinate, (23)To find solutions for S(r) we must assume numerical values for Ω_{m0}, Ω_{Λ0} and H_{0}. In this paper we adopt the FLRW model with Ω_{m0} = 0.3, Ω_{Λ0} = 0.7 and H_{0} = 70 km s^{1} Mpc^{1}.
The two differential Eqs. (20) and (23) comprise our numerical problem. Solving them simultaneously enables us to generate tables for r, S and N. A computer code using the fourthorder RungeKutta method is good enough to successfully carry out the numerical tasks. The initial conditions r_{0} and N_{0} are set to zero, whereas S_{0} can be derived considering that as r → 0 the spacetime is approximately Euclidean, that is, k ≈ 0. This leads, from Eq. (16), to ct = −r as well as S_{0} = 1.
3.2. Differential number count
The redshift z can be written as (24)where it is clear that a numerical solution of the scale factor S(r) immediately gives us the numerical solution for z(r). We can derive the differential number counts dN/dz by means of the expression, (25)with the help of the useful relation (26)Numerically, since we build all our quantities using N(r) and S(r), any derivative with respect to the redshift we wish to evaluate will be similarly written in terms of the derivative of that quantity in terms of the radial coordinate. The derivatives in Eq. (26) can be taken from definition (24) and Eq. (23), enabling us to write, (27)This, together with Eq. (20), allows us to write Eq. (25) as, (28)
3.3. Proper and comoving volumes
As discussed in RS03, most cosmological densities obtained from astronomical observations assume comoving volumes, whereas densities derived from theory often assume the local, or proper, volumes. The LF, for instance, is nowadays always obtained from galaxy catalogues by assuming the comoving volume in its calculation. So, we are only able to compare those observationally derived parameters with theory if we carry out a conversion of volume units. From metric (2) it is obvious that, (29)This equation clearly defines the conversion factor between these two volume definitions.
3.4. Distance measures
So far we have used only the area distance d_{A} as distance definition. However, other cosmological distances can, and will, be used later. They can be easily obtained from the area distance by invoking Etherington’s reciprocity law (Etherington 1933; Ellis 1971, 2007), (30)Here d_{L} is the luminosity distance and d_{G} is the galaxyarea distance. The latter is also known as angulardiameter distance, transverse comoving distance or propermotion distance. A fourth distance will also be useful later, the redshift distance d_{z}, defined by the following equation, (31)Although the reciprocity law is independent of any cosmological model, the detailed calculations presented in the previous sections are for the FLRW cosmology only. This means that in deriving d_{A}, d_{L}, d_{G} and d_{z} we write the FLRW expression for d_{A} and derive the others using the general, model independent, expression (30), as well as Eq. (31). Thus, considering Eq. (19) for d_{A} and the reciprocity theorem (30), it is straightforward to write the other cosmological distances in terms of the scale factor, as follows, The derivatives of each distance with respect to the redshift still need to be determined since they will be necessary later. Starting with the area distance d_{A}, they can be easily obtained from Eqs. (19) and (24). Then it follows that, (35)Using Eq. (23), this expression may be rewritten as, (36)The other observational distances can be numerically calculated from Eq. (35) if we consider the reciprocity law (30). Therefore, we have, These two equations can also be rewritten in terms of the scale factor if we consider Eqs. (24) and (36), yielding, From now on we shall generically use d_{i} to indicate any observational distance, which can be any one of the four cosmological distances defined above (i = A, G, L, Z).
3.5. Differential and integral densities
The differential density γ_{i} gives the rate of growth in number counts, or more exactly in their density, as one moves along the observational distance d_{i}. It is defined by (Ribeiro 2005; Albani et al. 2007; Rangel Lemos & Ribeiro 2008), (41)whereas the integrated differential density, or simply integral density, gives the number of sources per unit of observational volume located inside the observer’s past light cone out to a distance d_{i}. Is is written as, (42)where V_{i} is the observational volume, (43)These quantities are useful in determining whether or not, and within what ranges, a spatially homogeneous cosmological model can or cannot be observationally homogeneous as well (Ribeiro 1995, 2005; Rangel Lemos & Ribeiro 2008). This is because these densities behave very differently depending on the distance measure used in their definitions – that is, they show a strong dependence on the cosmological distance adopted. Therefore, as discussed in Rangel Lemos & Ribeiro (2008), these measures are the ones we employ in this paper, because they are capable of probing the possible observational inhomogeneity of the number counts.
From a numerical viewpoint it is preferable to write the densities given by Eqs. (41) and (42) in terms of the redshift. Thus, the differential density (41) may be written as, (44)A final point still needs to be discussed. The integral density (42) is a result of integrating γ_{i} over an observational volume. The simplest way of numerically deriving it is shown in what follows. Let us differentiate in terms of the volume, so that, (45)Considering Eq. (44), this result leads to the following expression, (46)Similarly, it is simple to show that, (47)Finally, from these two equations above, as well as from the definitions of γ_{i} and , it is easy to conclude that the following expression holds, (48)So, the numerical solution of Eq. (20) together with the numerical determination of all distances, as given by Eqs. (19), (32), (33), (34), allow us to calculate the volume (43) and evaluate . These results fully determine the numerical problem for the cosmological model under study.
In conclusion, once N(r) and S(r) are calculated along the past light cone in terms of the radial coordinate r, all other quantities are straightforwardly obtained with the same functional dependence: z(r), dN/dz(r), d_{i} (r), d(d_{i})/dz(r), γ_{i}(r), V_{i}(r) and .
4. Observational quantities
We shall now specialize the equations above to obtain their observational counterparts based on the LF data from a specific galaxy catalog. As a direct consequence, we show how these results can be used to test the consistency of the number count theory detailed above.
4.1. General equations
Generally speaking, we shall assume that an observational quantity [T] _{obs} can be related to its theoretical counterpart T by means of a completeness function J, such that, (49)RS03 showed that such a completeness function can be obtained by relating the selection function ψ, which gives the number of galaxies with luminosity above a given threshold in a given comoving volume, to the radial number density n_{C} (z) in comoving volume units – that is, the number of sources in a given comoving volume. As already mentioned, it is common practice to adopt this volume definition to calculate the LF from galaxy datasets. The volume number density n obtained in Eq. (15) comes from the right hand side of Einstein’s field equations and, therefore, it is written in terms of the proper volume. Thus, Eq. (15) must have its volume units corrected by means of Eq. (29) in order for the former to be correctly related to a selection function ψ stemming from the LF.
The relationship allowing us to calculate the completeness function used in this work may be written, (50)Now, if we go back to Eq. (14) and the discussion in the paragraph below it, we can see that applying the completeness function by means of Eq. (49) actually means replacing the theoretical comoving number density n_{C}, based on a constant average galaxy rest mass, with its observed counterpart ψ obtained directly from the LF. In that sense, all the observational quantities obtained by applying J(z) to their theoretical counterparts will inherit that same empirical number count redshift evolution encoded in the parametrization of the LF. As the goal of this paper is to consider the effect of redshift evolution on the relativistic number densities, such an approach should suffice. Also, in Eq. (50) the completeness function is independent of volume units, therefore, if an observational quantity is obtained using J(z) by means of Eq. (49) its original volume units dependence is preserved.
The selection function may be written in terms of the LF, (51)where l(z) is the lower luminosity threshold below which the sources are not observed. The theoretical radial number density in comoving volume is given by, (52)Most observational quantities of interest studied in this paper require previous knowledge of the observed differential number counts [ dN] _{obs}. Therefore, linking this quantity to its theoretical counterpart is an essential step in order to apply Ellis’ Eq. (13). In view of this, to find [ dN] _{obs} we identify T as dN in Eq. (49), yielding, (53)In essence this expression describes the same number counts as in Eq. (50) and both are in agreement with Ellis’ (1971) key Eq. (13). Considering now both Eqs. (50) and (53), it easily follows that, (54)This is our key equation relating the relativistic theory to the observations. The number density of a proper volume is given by (55)whereas the number density of a comoving volume is given by Eq. (52). Thus, we may rewrite Eq. (54) as follows, (56)It is important to point out that the volume transformation in this equation merely reflects the fact that J(z) is independent of volume units, proper or comoving, and does not in itself change the volume units underlying both [dN/dz]_{obs} and dN/dz when they are derived by means of Ellis’ Eq. (13). It is the number density n in Eq. (13) that defines the volume units of both [dN/dz]_{obs} and dN/dz (see below).
A further step towards writing these results in terms of the underlying cosmology is taken if we remember that RS03 showed that Ellis’ (1971) differential number counts (13) can be rewritten (57)Thus, Eq. (56) becomes (58)Note that in this expression [dN/dz]_{obs} is given in terms of proper volume units because the number density n is written in terms of proper volume units. Also note that these expressions are general and independent of any cosmological model. The specific cosmology will appear once we specialize the terms inside the brackets on the right hand side.
Galaxy redshift surveys use bandwidth filters, and they sometimes include some sort of morphological classification. Therefore, a selection function ψ which considers a set of filters W, and morphological types v, may be written as follows (see A07, Eqs. (3) − (10)), (59)Here P_{v} is the fraction of galaxies in the dataset that were classified with the morphological type v (see RS03, Eq. (13)), ℳ_{v} is the typical local restmass value for a galaxy of the morphological type v, and a_{W} are constants introduced to avoid multiple counting of the same objects in the various filters, defined as follows (see RS03, Eq. (46)), (60)and (61)Equation (61) simply states that when there exists more than one observed waveband (W > 1), a_{W} gives the fraction of galaxies in waveband W > 1 that are not counted in wavebands 1,2,...,(W − 1). Considering these expressions, the completeness function defined in Eq. (50) can be rewritten (62)Substituting the expression (59) in Eq. (58) we obtain the following result, (63)This expression already appeared in A07. Here it has been reached through a more straightforward and simpler derivation and with the volume dependence explicitly shown. This means that volume conversions like the one discussed above in Sect. 3.3 may be needed if we are to calculate the observational differential number count [ dN/dz] _{obs} in a consistent way. The above equation is shown here only to make explicit how the methodology of previous work can be understood in relationship to that being developed in this paper. In addition, it is important to note that Eqs. (54), (56) and (63) are closely related. The first two are written in a more compact form, whereas the last one expands the morphological and bandwidth dependencies of the selection function and the relativistic features of the underlying fourdimensional spacetime, as indicated by the expression (57).
Equation (63) also allows us to see how a theoretical mass evolution does not affect the differential number counts constructed with LF parameters. This is based on the form of the summation over W, inasmuch as if we allow an evolution of the galactic mass by considering it dependent on the redshift, such a dependency will appear both in the numerator and the denominator of the summation and, therefore, cancels out, at least to first order. However, although mass evolution does not directly alter Eq. (63), it will indirectly change [ dN/dz] _{obs} as the LF parameters used in its construction will include some form of source evolution. See Sect. 2.2 of A07 for more details.
4.2. Selection functions
To calculate the completeness function and, as a consequence, other observational quantities, we need to compute the selection functions with respect to redshift in a given galaxy survey. G04 and G06 fitted the LF Schechter parameters over the redshifts of 5558 Iband selected galaxies in the FDF dataset, photometrically measured down to an apparent magnitude limit of I_{AB} = 26.8. G04 and G06 also showed that the selection in the Iband is expected to miss less than 10% of the objects detected in the Kband, given that the ABmagnitudes of the Iband are half a magnitude deeper than those of the Kband, out to redshift 6, beyond which no signal is detectable in the Iband due to the Lyman break. Heidt et al. (2001) also argue that the Iband selection minimizes biases like dust absorption. All galaxies in those studies were therefore selected in the Iband and then had their magnitudes for each of the five blue bands (1500 Å, 2800 Å, u′, g′ and B) and the three red ones (r′, i′ and z′) computed using the best fitting SED given by their authors’ photometric redshift code convolved with the associated filter function. Gabash et al. (2004, 2006) determined the photometric redshifts by fitting template spectra to the measured fluxes on the optical and near infrared images of the galaxies. Heidt et al. (2001) reported that approximately 80% of the objects in that sample were classified as Im/Irrtype, 15% as E/S0type or Sa/Sb/Sctype, and finally 5% as stars.
The redshift ranges of 0.75 ≤ z ≤ 3.0 for the red bands and 0.5 ≤ z ≤ 5.0 for the blue bands are large enough for checking possible observational inhomogeneities, since Rangel Lemos & Ribeiro (2008) showed that in the Einsteinde Sitter cosmology it is only theoretically possible, i.e., not observationally, to distinguish these two features when z ≳ 0.1, whereas for z > 0.5 such a distinction becomes significant. Therefore, using galaxies whose measured redshifts are mostly greater than z = 1 allows us to obtain data along the past light cone far enough from our present time hypersurface. One must remember that spatial homogeneity is defined on our constant time hypersurface whereas observational homogeneity occurs along our past light cone (Rangel Lemos & Ribeiro 2008; see also Ribeiro 1992, 1995, 2001, 2005).
Selection functions for the FDF blue UV bands.
To calculate the actual values of these functions, we begin by writing the result obtained in RS03 for the limited bandwidth version of the selection function of a given LF, fitted by a Schechter analytical profile, in terms of absolute magnitudes, (64)where, as discussed earlier, the index W indicates the bandwidth filter. The equations for the redshift evolution of the LF parameters adopted by G04 and G06 are, \arraycolsep1.75ptwith A^{W} and B^{W} being the evolution parameters fitted for the different W bands and , and α_{0} the local (z ≈ 0) values of the Schechter parameters as defined in G04 and G06. Since all galaxies were detected and selected in the Iband, we can write, (65)for a luminosity distance d_{L} given in Mpc. I_{lim} is the limiting apparent magnitude of the Iband of the FDF survey and equals to 26.8. Its reddening correction is A^{I} = 0.035 (in Heidt et al. 2001). We computed the selection functions for all eight bands of the dataset by means of simple numerical integrations at equally spaced values spanning the whole redshift interval. The errors were propagated quadratically, as mentioned before. The results are summarized in Tables 2 to 4. Regarding the blue band dataset of G04, two distinct patterns of the selection functions in the different bands are noticeable. The UV bands, 1500 Å, 2800 Å, and u′, evolve tightly with redshift, having values that are consistent with each other within the uncertainties, while at the same time assuming values outside the uncertainties of those in the blue optical bands g′ and B. Therefore we chose to use the selection functions of the combined UV bands and those of the combined blue optical bands separately. The selection functions in the redband dataset of G06, r′, i′ and z′ are also combined. Once the combined selection functions have been obtained, we can calculate [ dN/dz] _{obs} by means of Eqs. (54) and (62).
Selection functions for FDF blue optical bands.
Selection functions for the FDF red bands.
Differential number counts.
4.3. Observed number counts
Following Eq. (48), the best way to obtain an observed relativistic number density [γ_{i}(z)] _{obs} for a given distance definition d_{i} is by calculating the observed number counts [N(z)] _{obs}. These can be written as, (66)The uncertainty of the cumulative number counts [ N] _{obs} can be derived from the already determined uncertainty in the number counts [ dN/dz] _{obs}. In this regard, the discussion about error analysis in A07 (Appendix) showed that, (67)where we can write, (68)The derivative of the completeness function dJ(z)/dz can be obtained from Eq. (50) and is proportional to the derivative of the selection function ψ^{W} itself with respect to redshift. This can be directly obtained from the parametrization of the LF.
Next, we shall use the observed differential number counts [dN/dz] _{obs} to study the relativistic radial distribution of the galaxies in the G04 and G06 datasets in the same manner as in A07.
5. Relativistic radial statistics of the FDF survey
The observational values of the differential number counts [dN/dz] _{obs} given in Table 5 for all filters in the dataset permit us to evaluate the observational differential number densities [γ] _{obs} for the various cosmological distance definitions in the FLRW model we are assuming. What it takes to do so is replacing dN/dz with [dN/dz] _{obs} in Eq. (44), which is exactly the same as applying the completeness function J(z) to the theoretical values of the relativistic differential number densities γ_{i}. As discussed at the end of Sect. 4.1, such application of the completeness function naturally includes the observed evolution of the number counts in the selection functions of a given dataset.
A look at Eq. (44) allows us to understand the separate roles of number counts and geometry, encoded in the term inside the brackets, in the relativistic differential densities γ_{i}. It shows that the underlying cosmological spacetime affects such quantities and that their theoretical departure from homogeneity at high redshifts is ultimately a consequence of the Universe’s expansion, through the inclusion of the scale factor S in the distance definitions, Eqs. (19), (32), and (34). The absence of the scale factor in the galaxy area distance d_{G} (Eq. (33)) allows us to identify this distance with the radial coordinate r in this particular spacetime. As a consequence, the differential number densities γ_{G} obtained using this distance definition are not subject to the redshift evolution caused by expansion and will remain constant if a constant average galaxy rest mass is assumed.
The fact that such differential densities depend on the geometry of the Universe through the term inside the brackets of Eq. (44) causes them to vary with the redshift differently, and they do not necessarily remain constant. Such inhomogeneities, however, are to be understood simply as a consequence of using relativistic distances, defined on the observer’s past light cone. They are therefore easily reconciled with the Cosmological Principle, which requires largescale spatial homogeneity. As can be seen in Figs. 5 and 6 of A07, both differential and integral number densities can be inhomogeneous, even when a constant average galaxy mass is assumed in an otherwise spatially homogeneous FLRW metric. Moreover, both quantities calculated using d_{G} remain homogeneous, as expected, since this distance in this particular spacetime is not affected by expansion, as discussed before.
We can use this important result to further compare the relativistic effects of considering densities down the past light cone in an expanding FLRW spacetime and the observational evolution of the LF in that same hypersurface. The argument is simple: if we compute the ratio of the theoretical values of a given relativistic density at a given redshift – say the differential number density using the luminosity distance d_{L} – to the one using d_{G}, this indicates how much observational inhomogeneity is introduced through the use of d_{L} in the calculation of this relativistic density. Similarly, the ratio of the observational value of [ γ_{G}] _{obs} to its theoretical one, γ_{G}, indicates how much observational inhomogeneity the redshift evolution of the LF introduces. In fact, considering Eqs. (44) and (54) we can see that such a ratio is simply the completeness function J(z) in a particular dataset, that is, (69)Now, if there were no significant departure from homogeneity at a given redshift such ratios should be approximately unity. Therefore, if we subtract the values of those ratios at a given redshift from one, this number should express how much inhomogeneity that relativistic density shows at that redshift: the bigger the number, the bigger the inhomogeneity. Table 6 summarizes these results. At a fixed redshift entry one can compare the effects on the homogeneity of the relativistic differential densities caused by the LF redshift evolution of the different filters in the G04 and G06 datasets.
Spacetime expansion vs. galaxy evolution effects on relativistic homogeneity.
The results in Table 6 allow us to see that the evolution of the LF in the combined UV bands produces a bigger departure from homogeneity of the differential densities than the one in the combined optical bands, and this departure is even bigger than the one in the combined red ones. On the other hand, fixing a given ratio, say the observationaltotheoretical ratio for the UV bands, [γ_{G}] _{obs}/γ_{G}, one can observe that the departure from homogeneity of these relativistic differential densities is quite a common feature and tends to increase with redshift. Finally, by comparing the values of the purely theoretical γ_{L}/γ_{G} ratio to the various observationaltotheoretical ones, one can investigate which contribution dominates at each redshift. That is particularly interesting, since the redshift evolution of any of the other differential densities can be obtained by a simple combination of those two separate effects, that of the expansion of the geometry in the distance definition and that of the redshift evolution of the LF. For instance, let us consider the observational differential density in the combined UV filters calculated using the luminosity distance d_{L}, that is, [γ_{L}] _{UV}. Its departure from homogeneity can be investigated using the ratio [γ_{L}] _{UV}/γ_{G}, as discussed above. Remembering that, by construction [γ_{L}] _{UV} = J_{UV} γ_{L}, it follows, (70)which is simply the combination of the observationaltotheoretical ratio for that dataset, namely the LF of the combined UV bands with the purely theoretical ratio for that distance definition, in that case, d_{L}. The values in Table 6 indicate that the geometrical effect is comparable to the evolution of the LF in the whole redshift range of the G04 and G06 datasets.
At this point, one should note that determining the spacetime geometry of the Universe or its density parameters is not the goal of the present paper. We simply assume the ΛCDM metric, a necessary step in calculating the above mentioned geometrical terms, in order to see how the empirical evolution of the differential number counts [dN/dz] _{obs} affects the theoretical results discussed above.
Similarly, we can obtain the observational values of the integral densities by means of Eq. (46), or (48), because, by their very definitions, we have that [N] _{obs} = V_{C}ψ(z).
Figures 1 and 2 show graphs of the observational differential densities determined using the combined UV and optical bands in the G04 dataset plotted against the redshift, whereas Fig. 3 shows the results for those same densities versus the redshift in the combined red bands. The dependence of such relativistic densities on the distance definition, a known theoretical result (Ribeiro 2005; A07), can be easily observed in these three figures. In addition, the dependence of [γ_{G}] _{obs} with the redshift in these same figures is solely due to the evolution of the mass function, since its relativistic volume definition remains constant in the assumed FLRW spacetime. We note that values of [γ_{G}] _{obs} obtained from the LF in the combined red bands (Fig. 3) seem less dependent on the redshift than those in the UV and optical bands (Figs. 1 and 2).
Fig. 1 Observational relativistic differential densities versus redshift in the combined UV bands of the FDF dataset of G04. Symbols are as in the legend. 

Open with DEXTER 
Fig. 2 Observational relativistic differential densities versus redshift in the combined optical bands of the FDF dataset of G04. Symbols are as in the legend. 

Open with DEXTER 
Fig. 3 Observational relativistic differential densities versus redshift in the combined red bands of the FDF dataset of G06. Symbols are as in the legend. 

Open with DEXTER 
Fig. 4 Observational relativistic integral densities versus redshift in the combined UV bands of the FDF dataset of G04. Symbols are as in the legend. 

Open with DEXTER 
Fig. 5 Observational relativistic integral densities versus redshift in the combined optical bands of the FDF dataset of G04. Symbols are as in the legend. 

Open with DEXTER 
Fig. 6 Observational relativistic integral densities versus redshift in the combined red bands of the FDF dataset of G06. Symbols are as in the legend. 

Open with DEXTER 
Fig. 7 Observational relativistic differential and integral densities versus luminosity distance in the combined UV bands of the FDF dataset of G04. The solid line was drawn just as reference, to show the apparent powerlaw behavior of both relativistic densities in that distance definition at high redshifts. Symbols are as in the legend. 

Open with DEXTER 
Fig. 8 Observational relativistic differential and integral densities versus luminosity distance in the combined optical bands of the FDF dataset of G04. The solid line is plotted just as a reference to the power law behavior stated in the legend. Symbols are as in the legend. 

Open with DEXTER 
Similar graphs for the observational integral densities in the UV, optical and red combined bands versus the redshift are respectively shown in Figs. 4–6, where their dependency with the distance definition can also be clearly seen. Although less pronounced in these smoothed out relativistic densities, the effect of the evolution of the mass function can still be observed in these three graphs, as is also unaffected by the geometrical effects of considering the densities on the past null light cone of an FLRW spacetime. In other words, the LF redshift evolution seems to further enhance the inhomogeneity of these densities, making observationally inhomogeneous even the theoretically homogeneous γ_{G} and . We observe that in Fig. 6, is also less dependent on the redshift than in the UV and optical bands shown in Figs. 4 and 5.
Figures 7–9 plot the observational differential and integral densities built with the luminosity distance in the combined UV, optical and red bands, respectively, versus the luminosity distance. A powerlaw pattern for both densities in the three combined bands can be observed when one compares the data points with the solid line drawn just as reference in all three plots. The slope of the reference lines is the same in all bands, which indicates that there is little change in the slope of the data points. This possibly indicates that such an effect is predominantly relativistic.
Fig. 9 Observational relativistic densities versus luminosity distance in the combined red bands of the FDF dataset of G06. The solid line is shown as a reference to the power law behavior stated in the legend. Symbols are as in the legend. 

Open with DEXTER 
Fig. 10 Observational relativistic differential and integral densities versus redshift distance in the combined UV bands of the FDF dataset of G04. The solid line is shown just as a reference to the power law behavior stated in the legend. Symbols are as in the legend. 

Open with DEXTER 
A similar pattern can be found in Figs. 10–12, which plot the observational differential and integral densities using the redshift distance against that same distance in the combined UV, optical and red bands respectively. We found a similar indication of a powerlaw pattern in both densities at high redshifts, but the actual power exponent of the distributions is different from the one using the luminosity distance. This comes from the fact that these two relativistic distances are affected differently by the expansion of the FLRW spacetime. Comparing the results shown in Fig. 11 with the ones in the UV band in Fig. 10 we can see that the slopes of the data points should agree with one another. Similarly as in Figs. 10 and 11, the power law exponent shown in Fig. 12 seems to be the same for both [ γ_{Z}] _{obs} and in the given distance definition.
Fig. 11 Observational relativistic differential and integral densities versus redshift distance in the combined optical bands of the FDF dataset of G04. The solid line is shown as a reference to the power law behavior stated in the legend. Symbols are as in the legend. 

Open with DEXTER 
Fig. 12 Observational relativistic differential and integral densities versus redshift distance in the combined red bands of the FDF dataset of G06. The solid line is plotted as a reference to the power law behavior stated in the legend. Symbols are as in the legend. 

Open with DEXTER 
One should notice that the slope of the straight lines used as references for the data points in Figs. 7 to 12 do not vary drastically over all bands of the FDF datasets, both for luminosity and redshift distance definitions. This possibly indicates that such an effect is predominantly relativistic.
The differential density constructed with the area distance d_{A} becomes discontinuous at z ≈ 1.5 due to the fact that, by definition, d_{A} has a maximum at that redshift and, therefore, its derivative with respect to z becomes zero, yielding an undefined [γ_{A}] _{obs} at high redshifts. In addition, differently from the decrease of and at higher redshifts, the integral density increases with z. Such pathologies have already been previously detected (see Figs. 5 and 6 of A07). They seem to render the number densities defined by d_{A} as unphysical.
The differential and integral densities defined by d_{G} lack the geometrical effect of expansion in the particular spacetime under consideration in this paper and, therefore, are understood to be unsuitable if one wants to fully characterize the relativistic powerlaw patterns arising from the combination of both spacetime expansion and LF evolution.
Finally, it is apparent that the differential densities [γ_{L}] _{obs} and [γ_{z}] _{obs} do not appear to exhibit a powerlaw decay as the integral ones do, especially at redshift values nearing the sample’s depth. This is very similar to what found in A07. This behavior is possibly a result of noisier data at the redshift limit of the samples. By definition, the differential densities measure the rate of growth in number counts, rendering them more sensitive to local fluctuations, whereas the integral densities indicate the change in number counts for entire observational volumes, rendering them less sensitive to the same fluctuations.
The indication of powerlaw behavior may be the result of observational, not necessarily spatial, inhomogeneity as we look down our past light cone, combined with the incompleteness of galaxy counts at higher redshifts. It may also be a result of other causes which we are investigating.
6. Conclusion
In this paper we have detailed a generalization of the framework connecting the relativistic cosmology number count theory with the astronomical data extracted from the galaxy luminosity function (LF), as proposed by Ribeiro & Stoeger (2003). This framework was used by Albani et al. (2007) to extract the observed differential number counts from the LF data and to study the observational inhomogeneities in the relativistic radial densities using different distance definitions. Here we have carried out an analysis of number counts using the LF dataset of Gabasch et al. (2004) in the redshift range 0.5 ≤ z ≤ 5.0, and of Gabasch et al. (2006) in the redshift range 0.75 ≤ z ≤ 3.0. In doing so we focused on the observational inhomogeneities in the relativistic radial statistics of the distribution of the galaxies by means of the differential densities γ_{i} (i = A, G, L, Z) and integral densities , both using the various cosmological distance definitions d_{i}, as well as the selection functions and the differential number counts [ dN/dz] _{obs} obtained from the LF of G04 and G06.
We confirmed the dependence of such empirically based relativistic densities on the distance definition, a geometrical consequence of the expansion of the Universe, which also is the cause of the observational inhomogeneities present even in theoretical analysis of the spatially homogeneous FLRW spacetime. Furthermore, as expected, the effect of the redshift evolution of the LF seems to increase the observational inhomogeneity of these densities.
We found evidence of a powerlaw pattern in the behavior of [γ_{L}] _{obs} and [γ_{z}] _{obs}, relative to their respective cosmological distance measures d_{L} and d_{G}, similar to that previously detected in A07. However, this pattern is somewhat less pronounced than the results of A07. This difference could be due to the use of a different, and somewhat more complete, redshift samples used in this paper, especially the red galaxies of G06.
Finally, we should note the dependence of these results and their interpretation on the assumed FLRW cosmology. Therefore, we look forward to finding out what an analysis of this data with reference to an inhomogeneous cosmological model (e.g., a LemaîtreTolmanBondi model) would reveal. If the universe is not almostFLRW or if the data are not probing almostFLRW scales, then hidden in the results are the effects of having chosen the wrong cosmological model. In addition, what we interpret as LF evolution, which is the same as number evolution, but with some indirect luminosity evolution effects, undoubtedly also reflects a lot of unobserved, and therefore uncounted, galaxies. Strictly speaking, if the differential and integral densities are comoving densities, then γ_{G} and should increase with redshift, as there should be more galaxies in a comoving volume in the past that there are now due to mergers (more, but smaller, galaxies in the past). So we are obviously not observing as many galaxies as we go farther out. Of course, that would mean that the average mass per galaxy decreases with redshift, an effect that should be implicit in the empirical LF, but which requires further theoretical considerations to be made explicit in our model. Including these points in our general approach to the problem of connecting relativistic cosmology theory and the empirically determined LF is indeed necessary in order to make the present study even more realistic. We intend to deal with these issues, as well as others like the possible meaning of power law patterns in the differential and integral densities, in forthcoming papers.
In this paper the term LF refers to the observationally determined quantity. It must be noted, however, that the LF does not depend only on observations, but also on an assumed cosmology. Thus, although it is a way of presenting observations, it also contains theoretical information. It could be considered to be fully observational if the assumed cosmology is observationally wellsubstantiated.
Usually a density is just a density, and the adjective “relativistic” does not apply. However, in this paper this expression has to do with the different ways the density can be defined, depending on the different measures of distance we use, since in relativistic cosmology a distance can be defined by different measures. In this context the adjective “cosmological” could also be used. See Sects. 3.4 and 3.5 below.
Acknowledgments
Thanks go to A. Gabasch for kindly providing the absolute magnitudes of their Iband selected dataset and thoroughly helping us with its observational details. We are also grateful to the referee for providing very helpful criticisms, recommendations, and remarks, which substantially improved the paper. Finally, we thank M. Giavalisco for some clarifying discussions.
References
 Albani, V. V. L., Iribarrem, A. S., Ribeiro, M. B., & Stoeger, W. R. 2007, ApJ, 657, 760 (A07) [NASA ADS] [CrossRef] (In the text)
 Bell, E. F., Mcintosh, D. H., Katz, N., & Weinberg, M. D. 2003, ApJS, 149, 289 [NASA ADS] [CrossRef] (In the text)
 Blanton, M. R., Hogg, D. W., Bahcall, N. A., et al. 2003, ApJ, 592, 819 [NASA ADS] [CrossRef] (In the text)
 Bouwens, R. J., Illingworth, G. D., Franx, M., & Ford, H. 2007, ApJ, 670, 928 [NASA ADS] [CrossRef] (In the text)
 Ellis, G. F. R. 1971, in General Relativity and Cosmology, ed. R. K. Sachs, Proc. Int. School Phys. Enrico Fermi (New York: Academic Press), reprinted in Gen. Rel. Grav., 41, 581, 2009 (In the text)
 Ellis, G. F. R. 2007, Gen. Rel. Grav., 39, 1047 [NASA ADS] [CrossRef] (In the text)
 Etherington, I. M. H. 1933, Phil. Mag. ser. 7, 15, 761; reprinted in Gen. Rel. Grav. 39, 1055, 2007 [NASA ADS] [CrossRef] (In the text)
 Fried, J. W., von Kuhlmann, B., Meisenheimer, K., et al. 2001, A&A, 367, 788 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Gabasch, A., Bender, R., Seitz, S., et al. 2004, A&A, 421, 41 (G04) [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Gabasch, A., Hopp, U., Feulner, G., et al. 2006, A&A, 448, 101 (G06) [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Gabasch, A., Goranova, Y., Hopp, U., et al. 2008, MNRAS, 383, 1319 [NASA ADS] [CrossRef] (In the text)
 Heidt, J., Appenzeller, I., Bender, R., et al. 2001, The FORS Deep Field, Reviews in Modern Astronomy, ed. R. E. Schielicke (Astronomische Gesellschaft), 14 (In the text)
 Lin, H., Yee, H. K. C., Carlberg, R. G., et al. 1999, ApJ, 518, 533 [NASA ADS] [CrossRef] (In the text)
 Ly, C., Malkan, M. A., Kashikawa, N., et al. 2007, ApJ, 657, 738 [NASA ADS] [CrossRef] (In the text)
 Norman, C., Ptak, A., Hornschemeier, A., et al. 2004, ApJ, 607, 721 [NASA ADS] [CrossRef] (In the text)
 Peacock, J. A. 1999, Cosmological Physics (Cambridge: Cambridge Univ. Press) (In the text)
 Plebański, J., & Krasiński, A. 2006, An Introduction to General Relativity and Cosmology (Cambridge University Press) (In the text)
 Poli, F., Giallongo, E., Fontana, A., et al. 2003, ApJ, 593, L1 [NASA ADS] [CrossRef] (In the text)
 Pozzetti, L., Cimatti, A., Zamorani, G., et al. 2003, A&A, 402, 837 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Rangel Lemos, L. J., & Ribeiro, M. B. 2008, A&A, 488, 55 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Ribeiro, M. B. 1992, ApJ, 395, 29 [NASA ADS] [CrossRef] (In the text)
 Ribeiro, M. B. 1995, ApJ, 441, 477 [NASA ADS] [CrossRef] (In the text)
 Ribeiro, M. B. 2001, Gen. Rel. Grav., 33, 1699 [NASA ADS] [CrossRef] (In the text)
 Ribeiro, M. B. 2005, A&A, 429, 65 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Ribeiro, M. B., & Stoeger, W. R. 2003, ApJ, 592, 1 (RS03) [NASA ADS] [CrossRef] (In the text)
 Roos, M. 1994, Introduction to Cosmology (Chichester: Wiley) (In the text)
 Rudnick, G., Rix, H.W., Franx, M., et al. 2003, ApJ, 599, 847 [NASA ADS] [CrossRef] (In the text)
 Schechter, P. 1976, ApJ, 203, 297 [NASA ADS] [CrossRef] (In the text)
 Sparke, L. S., & Gallagher, J. S. 2000, Galaxies in the Universe (Cambridge University Press) (In the text)
 Stoeger, W. R., Ellis, G. F. R., & Nel, S. D. 1992, Class. Q. Grav., 9, 509 [NASA ADS] [CrossRef] (In the text)
 Tzanavaris, P., & Georgantopoulos, I. 2008, A&A, 480, 663 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Wilmer, C. N. A., Faber, S. M., Koo, D. C., et al. 2006, ApJ, 647, 853 [NASA ADS] [CrossRef] (In the text)
Appendix A: Numerical algorithm
This appendix describes step by step how to apply the theoretical and numerical methodologies described in sections Sects. 2–4 in order to obtain the results used in this paper. Steps 1 to 6 involve the calculation of the theoretical quantities, as discussed in Sects. 2 and 3. Steps 7 to 10 show how to obtain the corresponding observational quantities for a given LF, as discussed in Sect. 4.
Step 1 – Coupled differential equations. Obtain the scale factor S(r) and the cumulative number count N(r) as a function of the radial coordinate, or comoving distance r, by solving the coupled Eqs. (20) and (23). We used the fourth order RungeKutta method to obtain values of S(r) and N(r) at fixed r increments, assuming as initial conditions N(0) = 0 and S(0) = 1.
Step 2 – Auxiliary quantities. Obtain the redshift z(r) for each increment of r using Eq. (24) and the corresponding value of the scale factor S(r). Do the same for the area distance d_{A}(r) using Eq. (19), the differential number count dN/dz(r) using Eq. (28), and the comoving volume number density n_{C}(r) using Eq. (52).
Step 3 – Cosmological distances. Having the values of z(r) and d_{A}(r) calculated for each r increment, obtain the values of the galaxy area distance d_{G}(r), and the luminosity distance d_{L}(r) by direct application of Eq. (30), Etherington’s reciprocity law itself. The values of the redshift distance d_{z}(r) can be easily obtained for each r increment using Eq. (31).
Step 4 – Derivatives of the distances with respect to redshift. For each r increment, calculate the values of the derivatives with respect to redshift of the distances obtained in the previous step by means of Eqs. (36), (39) and (40).
Step 5 – Theoretical differential densities. Obtain the values of the theoretical differential densities γ_{i}(r) for each of the distance definitions d_{i} (i = A,G,L,Z) calculated previously for each r increment, by means of Eq. (44).
Step 6 – Theoretical integral densities. Calculate for each distance d_{i}(r) entry its corresponding spherical volume V_{i}(r), by means of Eqs. (43), then use those values, together with the ones for the cumulative number count N(r), to obtain the values of the integrated differential densities , by means of Eq. (48).
Step 7 – Selection functions. Select a bin containing a group of redshift values in a certain interval [z_{a → b}] , from a point z_{a} to a point z_{b}, taken from the ones obtained in step 2, and obtain the corresponding values of the selection functions ψ^{W}(z_{a → b}) for each filter W using Eq. (64), with the appropriate values for the LF’s Schechter parameters φ^{∗}(z_{a → b}), M^{∗}(z_{a → b}), and α(z_{a → b}), and different absolute magnitude limits , given by Eq. (65).
Step 8 – Completeness functions. Calculate the values of the completeness function J^{W}(z_{a → b}) by means of Eq. (50) using the selection functions ψ^{W}(z_{a → b}) calculated in step 7 and the corresponding values of the comoving volume number density n_{C} [ z_{a → b}(r)] obtained in step 2.
Step 9 – Differential and cumulative number counts. Obtain the observational differential number count in each filter W by means of Eq. (54) and the observational cumulative number count , by means of Eq. (66).
Step 10 – Observational differential and integral densities. Obtain the observational differential densities , and the integrated differential densities , by using the values of and in place of their theoretical counterparts dN/dz and N in Eqs. (44) and (48), respectively.
All Tables
All Figures
Fig. 1 Observational relativistic differential densities versus redshift in the combined UV bands of the FDF dataset of G04. Symbols are as in the legend. 

Open with DEXTER  
In the text 
Fig. 2 Observational relativistic differential densities versus redshift in the combined optical bands of the FDF dataset of G04. Symbols are as in the legend. 

Open with DEXTER  
In the text 
Fig. 3 Observational relativistic differential densities versus redshift in the combined red bands of the FDF dataset of G06. Symbols are as in the legend. 

Open with DEXTER  
In the text 
Fig. 4 Observational relativistic integral densities versus redshift in the combined UV bands of the FDF dataset of G04. Symbols are as in the legend. 

Open with DEXTER  
In the text 
Fig. 5 Observational relativistic integral densities versus redshift in the combined optical bands of the FDF dataset of G04. Symbols are as in the legend. 

Open with DEXTER  
In the text 
Fig. 6 Observational relativistic integral densities versus redshift in the combined red bands of the FDF dataset of G06. Symbols are as in the legend. 

Open with DEXTER  
In the text 
Fig. 7 Observational relativistic differential and integral densities versus luminosity distance in the combined UV bands of the FDF dataset of G04. The solid line was drawn just as reference, to show the apparent powerlaw behavior of both relativistic densities in that distance definition at high redshifts. Symbols are as in the legend. 

Open with DEXTER  
In the text 
Fig. 8 Observational relativistic differential and integral densities versus luminosity distance in the combined optical bands of the FDF dataset of G04. The solid line is plotted just as a reference to the power law behavior stated in the legend. Symbols are as in the legend. 

Open with DEXTER  
In the text 
Fig. 9 Observational relativistic densities versus luminosity distance in the combined red bands of the FDF dataset of G06. The solid line is shown as a reference to the power law behavior stated in the legend. Symbols are as in the legend. 

Open with DEXTER  
In the text 
Fig. 10 Observational relativistic differential and integral densities versus redshift distance in the combined UV bands of the FDF dataset of G04. The solid line is shown just as a reference to the power law behavior stated in the legend. Symbols are as in the legend. 

Open with DEXTER  
In the text 
Fig. 11 Observational relativistic differential and integral densities versus redshift distance in the combined optical bands of the FDF dataset of G04. The solid line is shown as a reference to the power law behavior stated in the legend. Symbols are as in the legend. 

Open with DEXTER  
In the text 
Fig. 12 Observational relativistic differential and integral densities versus redshift distance in the combined red bands of the FDF dataset of G06. The solid line is plotted as a reference to the power law behavior stated in the legend. Symbols are as in the legend. 

Open with DEXTER  
In the text 