Before to start we briefly review the Hubble law, after  , which is a linear relationship between expansion velocity and distance in
where is the Hubble constant , with when h is not specified, D is the distance in Mpc, c is the light velocity and z is the redshift. As an example a recent evaluation, see  , quotes
The original Hubble’s law was based on ≈20 galaxies of which was known the distance and the velocity; conversely in the last years the number of catalogs for galaxies with redshift available for public downloading has progressively grown. This huge amount of data allows setting up critical tests between different theoretical models over the various aspects of the Large Scale Structures. The layout of the paper is as follows. In Section 2, we describe some representative catalogs of galaxies. In Section 3, we introduce the adopted luminosity function for galaxies. In Section 4, we model the maximum in the number of galaxies as function of the distance. In Section 5, we evaluate the gravitational field as produced by the visible galaxies of Cosmicflows-2 catalog. We conclude in Section 6.
2. The Catalogs of Galaxies
A first kind of catalog for galaxies is represented by those focused as slice as the two-degree Field Galaxy Redshift Survey, in the following 2dFGRS, see  , or the Sloan Digital Sky Survey (SDSS), see  . A second classification is about the all-sky catalogs such as the 2MASS Redshift Survey (2MRS), see  , or the Cosmicflows-2, see  . As a first example a strip of the 2dFGRS is shown in Figure 1. The below Figure, which is a 3˚ slice, shows that the galaxies resides on filaments rather to be distributed in a uniform 2D way. A first extrapolation allows to state that the galaxies are disposed on the surface of bubbles rather than to fill uniformly the 3D space. The statistics of the voids is a topic of research and an important parameter is the average radius of the voids, :  quotes and Table 1 in  quotes a variable average radius ranging from to according to the selected sample. As second example is represented by the overall Cosmicflows-2 catalog, data available at http://vizier.u-strasbg.fr/viz-bin/VizieR, and Figure 2 reports the sky distribution of such a catalog in Galactic coordinates.
Also in this case the voids are can be visualized selecting a thin (20 Mpc) spherical shell, see Figure 3.
3. Luminosity Function for Galaxies
The distance modulus is
Figure 1. Cone-diagram of the galaxies in the 2dFGRS with distance Î [4 Mpc, 470 Mpc]. This plot contains 117,293 galaxies.
Figure 2. Sky distribution of 4970 Cosmicflows-2’s galaxies in galactic coordinates projected using the Mollweide projection.
Figure 3. Sky distribution of 713 Cosmicflows-2’s galaxies with distance Î [63.8 Mpc, 83.8 Mpc].
where m is the apparent magnitude, M is the absolute magnitude and d is the distance in pc.
Let L, the luminosity of a galaxy, be defined in . The Schechter LF of galaxies, , see  , is
where sets the slope for low values of , is the characteristic luminosity, and represents the number of galaxies per Mpc3. The normalization is
is the Gamma function. The average luminosity, , is
An equivalent form in absolute magnitude of the Schechter LF is
where is the characteristic magnitude.
A typical result of the Schechter LF in the case of Cosmicflows-2 is reported in Figure 4.
4. The Photometric Maximum
The flux, f, is
where r is the distance and L the luminosity of the galaxy. The joint distribution in distance, r, and flux, f, for the number of galaxies is
were the factor ( ) converts the number density into density for solid angle
and the Dirac delta function selects the required flux. In the case of Schechter LF of galaxies the number of galaxies as function of the distance is
We now introduce the critical radius
Figure 4. The observed Schechter LF for galaxies of Cosmicflows-2, empty stars with error bar, and the fit by the Schechter LF when , and .
Therefore the joint distribution in distance and flux becomes
The above number of galaxies has a maximum at :
and the average distance of the galaxies, , is
Figure 5 presents the number of galaxies observed in Cosmicflows-2 as a function of the distance for a given window in flux, as well as the theoretical curve. The luminosity of a galaxy is produced and attenuated according to the electromagnetism. We now deal with the mass of a galaxy. A transformation of the luminosity of a galaxy, L, by the mass, , is given by the nonlinear Formula (9) in 
Figure 5. The galaxies of Cosmicflows-2 with are organized by frequency versus distance, (empty circles); the error bar is given by the square root of the frequency. The maximum frequency of the observed galaxies is at Mpc. The full line is the theoretical curve generated by as given by the application of the Schechter LF which is Equation (11) and the theoretical maximum is at Mpc. The parameters are and .
We pay careful attention to the interval of existence in absolute magnitude versus distance for Cosmicflows-2 catalog, see Figure 6.
The progressive decrease of the interval of existence for the absolute magnitude is known as Malmquist bias  . Conversely the mass of a galaxy does not disappear with distance and therefore produce the gravitational field even if the galaxy is not visible due to the instrumental limitations.
5. Gravitational Forces
The shell theorem, according to proposition LXX, theorem XXX, in Principia  states that “If to every point of a spherical surface there tend equal centripetal forces decreasing in the square of the distances from those points, I say, that a corpuscle placed within that superficies will not be attracted by those forces any way” or more simply “The gravitational field inside a uniform spherical shell is zero”    . The galaxies are thought to reside on irregular shells which are approximated by a given radius. Due to the fact that the shapes of the above shells are neither exactly spherical nor uniformly populated by galaxies we search for a low value of the gravitational field rather than the theoretical zero. We therefore analyze a 2D box with size of 30 Mpc and masses as given by formula 16. The forces in every point of the considered box are evaluated according to the Newtonian force where G, the Newtonian constant of gravitation, is
Figure 6. The absolute magnitude M of 4970 galaxies belonging to Cosmicflows-2 when the absolute bolometric magnitude is 3.39 and (green points). The upper theoretical curve, Equation (3), is reported as a red thick line when .
where the length is in Mpc, mass is in which are and yr 8 are 108 yr, see more details in  . The masses do not disappear with distance and the above value of the box allows processing a complete sample of galaxies. Figure 7 reports the values of the gravitational field as a cut in the middle of the 2D box. Figure 8 reports the values of the gravitational field organized as a two color map.
Figure 7. Cut-line of the 2D gravitational forces expressed in (decimal logarithm).
Figure 8. Color scheme for the values of the gravitational field evaluated in 500 × 500 points. The red zone has values of gravitational field in the interval (the zones near the galaxies) and the green zone has values of gravitational field in the interval (the zones of the voids).
The observed “great wall” is theoretically explained by a maximum in the number of galaxies as function of the distance in Mpc, see Figure 5. The two concepts of repeller and attractor, see  , are not used in our analysis. The presence of voids approximated by spheres in the spatial distribution of galaxies allows analyzing the gravitational field at the light of the shell theorem. Due to the facts that the number of galaxies is finite and the masses of galaxies are distributed according to a gamma probability density function the gravitational field at the center of the voids turns out to be rather than zero.
This research has made use of the VizieR catalogue access tool, CDS, Strasbourg, France.