Logarithmic corrections to Newtonian gravity and Large Scale Structure

Effects from nonstandard corrections to Newtonian gravity, at large scale, can be investigated using the cosmological structure formation. In particular, it is possible to show if and how a logarithmic correction (as that induced from nonlocal gravity) modifies the clustering properties of galaxies and of clusters of galaxies. The thermodynamics of such systems can be used to obtain important information about the effects of such modification on clustering. We will compare its effects with observational data and it will be demonstrated that the observations seem to point to a characteristic scale where such a logarithmic correction might be in play at galactic scales. However, at larger scales such statistical inferences are much weaker, so that a fully reliable statistical evidence for this kind of corrections cannot be stated without further investigations and the use of more varied and precise cosmological and astrophysical probes.


Introduction
The clustering of galaxies is the main mechanism to address the large scale structure of the Universe. Such clustering mechanism can be studied confronting nua e-mail: capozziello@unina.it b e-mail: mirfaizalmir@googlemail.com c e-mail: hme123eda@gmail.com d e-mail: b.pourhassan@du.ac.ir e e-mail: vincenzo.salzano@usz.edu.pl merical simulations [1,2] with observations [3,4], using the local matter distributions of galaxies organized in groups, filaments and clusters. However, to analyze numerically the large scale structure formation, it is possible to approximate galaxies as point-like particles, and then study the clustering adopting the standard formalism of statistical mechanics. This approximation is valid because the distance between galaxies is much larger than their proper size. In fact, it has been demonstrated that the clustering mechanism can be studied considering a quasi-equilibrium description because the macroscopic quantities change very slowly as compared with the local relaxation time scales [5][6][7][8]. Pressure, average density and average temperature of clusters are the macroscopic quantities to be taken into account. Here the temperature is the effective one obtained by the kinetic theory of gases, for a gas of galaxies (each galaxy being approximated as a point particle). So, by this quasi-equilibrium description, the clustering of galaxies can be dealt as a thermodynamic system [9, 10]. Specifically, one has to adopt a gravitational partition function [11,10] where point-like galaxies gravitationally interact. Such a partition function can be also used to study phase transitions for systems of interacting galaxies [12][13][14].
It has to be pointed out that such a gravitational partition function diverges because the extended structure of galaxies has been neglected in this approximation. However, these divergences can be removed by adding a softening parameter, which accounts for the extended structure of galaxies as reported in [15]. In a recent study it has also been proved that the nonlocal gravity is divergence free [16]. Thermodynamics for systems of galaxies is obtained from the partition function (regularized with a softening parameter), and this can in turn be used to study the clustering properties [17]. It has been observed that any modification of gravitational partition function by a softening parameter, in turn, modifies the thermodynamic quantities of a given system. Vice versa, because these thermodynamic quantities are related to the clustering parameter, it also changes the clustering parameter itself. In fact, it can be demonstrated that the modification to the skewness and kurtosis of clustering systems occur due to such a softening parameter [18]. This formalism has also been used to demonstrated that galaxy clusters are surrounded by halos, and this has been done adopting a wide range of samples [19,20]. Furthermore the isothermal compressibility gives also information on the clustering of galaxies [21].
Specifically, the clustering of galaxies occurs due to the gravitational force between them. However, it is now well-known by a wide range of cosmological observations (e.g. Type Ia Supernovae, Baryon Acoustic Oscillations and Cosmic Microwave Background) that the Universe is undergoing an accelerating expansion [22][23][24]. This cosmic expansion is pulling the galaxies away from each other. So, the cosmic expansion seems to oppose the gravitational force, and thus it is important to include its effect in the gravitational partition function. It has been argued that the cosmological constant related to the accelerated expansion can produce important local effects on the scale of galaxy clusters [25][26][27][28][29]. These effects can be investigated using model of interacting dark matter and dark energy [25,26]. It is also possible to constrain the dark energy from the clustering of galaxies [30]. The cosmic expansion can also change the dynamics at the scale of galaxy clusters [31]. Such a modification of cluster dynamics has also been investigated using different groups of galaxies [32]. Finally, one can state that the cosmological constant has consequences for the formation of galaxy clusters [33]. In summary, it is important to study the effects of cosmic expansion at the scale of clusters. From a more technical point of view, it has to be noted that the Fisher matrix formalism has been used to analyze the effects of cosmic expansion on clusters [34].
Thus, it is possible to obtain information about the cosmological constant, and, in general, about dark energy dynamics, from the local clustering of galaxies. So, it is important to incorporate the effects of cosmic expansion in the statistical mechanical description of the clustering of galaxies. In fact, an accurate measurement of the statistics of galaxies can be used to constraint the value of the cosmological constant [35], and the distribution function of galaxies can be used to constraint the amount of the dark energy [36].
As such a statistical distribution function can be obtained from the gravitational partition function [9, 10], it is important to incorporate the effects of cosmic expansion in the partition function. So, the modification of the gravitational partition function by the cosmological constant term has also to be considered, and this modified gravitational partition function can be used to study the clustering process in an expanding Universe [37]. This can be technically done by first analyzing the Helmholtz free energy and entropy for a system of galaxies. Then the clustering parameter for this system can be obtained, and it has been observed that this clustering parameter depends on the value of the cosmological constant. The modification of the gravitational partition function with a time dependent cosmological constant can be also used to study the effects of dynamical dark energy on clustering of galaxies [38]. It has been observed that the correlation function for this system are consistent with observations. However, it is possible to obtain cosmic expansion without considering dark energy related to some new fundamental particle. In fact, large scale nonlocality related to some effective gravitational potential seems a realistic effect of gravitational field at infrared scales as recently reported in several papers, see for example [39][40][41][42][43]. It has been argued that such nonlocalities can occur due to some quantum gravitational effects [39,40,44]. They can be relevant also in gravitational waves modes and polarizations [45,46]. It has been demonstrated that string field theory produces nonlocal actions for the component fields, and these nonlocal actions can also be used to construct nonlocal cosmological models [47,48]. It is also possible to construct models of nonlocal cosmology from loop quantum gravity [49,50]. This can be done adopting a sort of condensation of states in loop quantum gravity. In summary, it is possible to obtain nonlocal cosmologies using quantum gravitational effects at ultraviolet scales that propagate up to infrared scales.
In other words, as nonlocal cosmology can explain the cosmic expansion without dark energy [39][40][41][42][43], such nonlocal corrections can also be used to incorporate the effects of cosmic expansion in the gravitational partition function. This one can then be used to analyze the effects of nonlocality on clustering of galaxies. It may be noted that even though general relativity is strongly constrained at the Solar System scales by several experiments and observations [51][52][53], nonlocal modification to general relativity can occur at large scales starting 3 from galactic one [54]. At the same time, nonlocal gravity is safe from possible issues related to oscillating features on small scales, falling in the regime where such oscillations are subdominant [55]. It has also been proposed that such nonlocality could be constrained from the clustering of galaxies [56]. Thus, it is important to analyze the effect of nonlocal modification of gravity on the gravitational partition function and point out how it affects the clustering process. Such a modification can be analyzed using nonlocal corrections to the gravitational potential [54,[57][58][59][60]. In general, the nonlocal gravitational partition function can be used to analyze the effects of nonlocality on clustering of galaxies. Although our main motivation for such a study is related to testing nonlocal gravity, we must note here that the logarithmic correction induced by such a theory is shared also by other alternative models of gravity, e.g. MOND [61,62], Randall-Sundrum brane scenario [63], or other string-theoretical motivations [64] which have been also tested with internal dynamics of spiral galaxies [65]. It may be noted that there is an intrinsic non-locality in string theory [66,67]. Furthermore, even MOND has been related to nonlocality [68][69][70][71][72]. This could be due to the fact that nonlocality [39,40], and several of these modifications, are motivated from quantum corrections to the original Einstein-Hilbert action. Even MOND can be obtained as a quantum correction in the Verlinde formalism [73]. Thus, we expect some kind of universal form of correction to occur in most of them. In this perspective, our results will be more general and not strictly related to one particular scenario.
In Sec. 2, we discuss the clustering mechanism. After defining the nonlocal gravitational partition function, we consider the related thermodynamics and spatial correlation function. Observational data analysis is described in detail in Sec. 3. Specifically, we take into account data from galaxies and clusters and study the effects of nonlocality on the clustering mechanism. Results are discussed in Sec. 4, while conclusions are drawn in Sec. 5.

Modeling Clusters of galaxies
In general, a nonlocal modification of the Hilbert-Einstein action can have the following form [39,40] where R is the Ricci scalar and f ( −1 R) is an arbitrary function, called distortion function, of the nonlocal term −1 R, which is explicitly given by the retarder Green's function (2) Setting f ( −1 R) = 0, the above action is equivalent to the Einstein-Hilbert one. The nonlocality is introduced by the inverse of the d'Alembert operator. In order to "localize" the action, an auxiliary scalar field can be introduced so that −1 R = φ and then formally φ = R. As a consequence, characteristic lengths related to the nonlocal terms naturally come out [54]. The weak-field limit of such a dynamics gives rise to corrections to the Newtonian potential which are interesting at large scales. These corrections can be polynomial or logarithmic and, in general, introduce characteristic scales. For a detailed discussion on this point see [54].
In the framework of this theory, let us now model the clustering of galaxies using the gravitational partition function which we are going to define. We will explicitly adopt the nonlocal gravity, as discussed in [58][59][60], to obtain the corrections to such a gravitational partition function. After analyzing the effects of nonlocality, the nonlocal gravitational partition function will be used to analyze the clustering of galaxies. This can be done by considering the thermodynamics of this system, and then relating it to the clustering parameter. Then we will calculate the spatial correlation function.

Nonlocal Gravitational Partition Function
Let us take into account a system with a large number of galaxies. This system can be analyzed using a quasiequilibrium description, where the change in macroscopic quantities is slower than the local relaxation time scales [5][6][7][8]. It is possible analyze such a system adopting an ensemble of cells, with same volume V , or radius R 1 and average densityρ. As the number of galaxies and their total energy can vary between these cells, the system can be analyzed in the framework of the grand canonical ensemble. Thus, it is possible to define a gravitational partition function. In this picture, a system of N galaxies of mass m interacting with a nonlocal gravitational potential can be written as [11,10] where p i are the momenta of different galaxies and Φ nl is the nonlocal gravitational potential energy. Here T is 4 the average temperature, which is obtained from the kinetic theory of gases where each particle is represented by a galaxy. Now integrating on the momentum space, we can write this gravitational partition function as where Q N (T, V ), the configuration integral, can be expressed as The nonlocal gravitational potential energy Φ nl (r 1 , . . . , r N ) is a function of the relative position vector r ij = |r i −r j |, as both the local and nonlocal terms represent central forces. The total potential energy can be obtained by summing up the potential energies. So, we can write the total nonlocal potential energy Φ nl (r 1 , r 2 , . . . , r N ) as We can write such a nonlocal gravitational potential energy for galaxies as [58][59][60] Since galaxies have an extended structure, we need to modify the gravitational potential energy by a softening parameter ǫ [15,21]. Thus, we can write Now we can use a nonlocal two-point interaction function for galaxies f ij = e −(Φij) nl /T − 1 to analyze this system. It vanishes in absence of interactions, and it is non-zero only for interacting galaxies. Thus, we can now express the nonlocal configuration integral as where we can express the nonlocal two-point interaction Mayer function as It is worth noticing that such an expression for clustering integral has been discussed for local gravitational interactions [18,21]. Here we have obtained the contributions to such Mayer function from nonlocal gravitational interactions [58][59][60]. We can expand this nonlocal two-point interaction Mayer function as a power series. So, we can write f ij as Now using l = Gm 2 /(T λ), we can write Q 2 (T, V ) as where α nl is given by It may be noted here that F (3/2, n + l/2; 5/2; −R 2 1 /ǫ 2 ) is the hypergeometric function. Now using R 1 ∼ ρ −1/3 , we observe that 3Gm 2 /2R 1 T = 3Gm 2 /2ρ −1/3 T = 3G m 2 ρ 1/3 T −1 /2. So, using the scale invariance, ρ → λ −3 ρ, T → λ −1 T and r → λr, we can define x as where β = 3(Gm 2 ) 3 /2. From this expression for x, we can obtain any general configuration integral. Using the so called dilute approximation, we can assume (Φ i,j ) nl /T very small and take only the first term of the exponential expansion, Thus, calculating the general configuration integral, we obtain Now the gravitational partition function for this dilute gravitating system can be written as This expression for the gravitational partition function can be used to analyze the effects of nonlocality on the clustering of galaxies.

Thermodynamics
It has to be noted that the gravitational partition function has already been used to study the thermodynamics of systems of galaxies, see [9] and [10]. We can thus adopt the nonlocal gravitational partition function to analyze the thermodynamics of this nonlocal system. The Helmholtz free energy In the plots of Fig. 1, we analyze the dependence of the Helmholtz free energy on various parameters. The dependence of Helmholtz free energy on x is represented in Fig. 1 (a). We can observe that the Helmholtz free energy is a decreasing function of x for large value of the nonlocality parameter λ and an increasing function of x for smaller value of the nonlocality parameter λ. So, the value of nonlocality parameter changes the dependence of the Helmholtz free energy on x. In Fig. 1 (b), the Helmholtz free energy is plotted in terms of λ. We can see how there is a minimum for the Helmholtz free energy, and this minimum seems to occur at a specific value of λ. This minimum does not seem to change with the temperatures of the system. This observation can explain the change in the behavior of the Helmholtz free energy with the change in the value of λ, as was observed in Fig. 1 (a). In Fig. 1 (c), we plot the dependence of the Helmholtz free energy on the temperature of the system. Here the temperature corresponds to the temperature of a gas of galaxies, with each galaxy acting as a particle analog. We observe that generally the Helmholtz free energy increases with the temperature, but the rate of this increase depends on the number of galaxies. As we have plotted the dependence on temperature by considering the system as a gas of galaxies, we also analyze the dependence of the Helmholtz free energy on the number of galaxies in Fig. 1 (d). We observe that the Helmholtz free energy first decreases for a relatively small number of galaxies, becoming negative, and then increases as the number of galaxies becomes more than a certain value.
The entropy S can now be calculated from this Helmholtz free energy, as Now, for large N , using N − 1 ≈ N , we obtain This is the general clustering parameter for a system of galaxies interacting in presence of nonlocal gravity. The internal energy U = F + T S of a system of galaxies can now be written as We have plotted the internal energy of a system of galaxies in Fig. 2. Its behavior is analyzed by investigating its dependence on various parameters. In Fig. 2 (a), we investigate how the internal energy of the system varies with x. We observe that it is an increasing function of x for smaller values of λ, and a decreasing function of x for larger values of λ. In Fig. 2 (b), we analyze the dependence of the internal energy on the nonlocality parameter λ. We observe that there is a minimum for internal energy, and again this minimum does not seem to depend on the temperature of the system. We also investigate the dependence of the internal energy on the temperate of the system in Fig.  2 (c). It is observed that there is a discontinuity in the dependence of internal energy on temperature, which occurs independently of the number of galaxies. In Fig.  2 (d), the dependence of internal energy on the number of galaxies is plotted. It is observed that internal energy increases with the increase in the number of galaxies. However, the rate of increase is larger for larger values of x.
Similarly, we can write the pressure P and chemical potential µ for galaxies interacting through a nonlocal gravitational potential as The probability of finding N galaxies can be written as where Z G is the grand-partition function defined by and z is the activity. Thus, for a system of gravitationally interacting galaxies in nonlocal gravity , we can write Now the grand-partition function can be written as and the distribution function can be expressed as In Fig. 3, we analyze the behavior of the distribution function for galaxies, and its dependence on different parameters used here. In general, the distribution function has the universal feature of having a maximum, while its functional shape changes, and depends on the considered parameters.

Spatial correlation function
It is well know that the interaction between different galaxies will cause correlation in their positions. The integral of the correlation function over a certain volume can be expressed in terms of the mean square fluctuation of the total number of galaxies in a given volume. Thus we can write This can be further represented in terms of thermodynamic quantities as Taking the volume derivative and using the equation of state for pressure P , we get The pressure can be written in terms of the clustering parameter as We can write the change of volume with pressure at constant temperature The expression for ξ is where The volume integral of correlation function is the quantity which can be compared with observational data in order to detect possible effects of nonlocal gravity.

Observational Data Analysis
We will test the nonlocal gravitational clustering with the cluster catalog in [74], containing 132, 684 clusters of galaxies in the redshift range 0.05 ≤ z < 0.8 from the Sloan Digital Sky Survey III (SDSS-III). We will take advantage of this catalog because it has all the needed ingredients to analyze clustering properties of both galaxies and clusters (although an alternative more recent version is in [75]), using two different approaches.

Galaxies
The analysis of the clustering of galaxies will be conducted with a sort of "smart" version of the standard counts-in-cell procedure. In general, one would first need to define cells/volumes with a given size, and then to count the gravitational structures of interest within them. Here, our "smart" volumes will be the clusters identified in the SDSS-III survey, and we will count the galaxies which are within each of them. From this procedure, we will be able to calculate the observed distribution function F (N ) and to compare it with the theoretical expectation in nonlocal gravity as given by Eq. (29). The "smartness" in this procedure is that the cells/volumes we are going to analyze are physically realized systems, structures which exist and have been identified, contrarily to the standard way of proceeding, where the existence of clusters and interacting systems of galaxies is not assured neither verified. On the other hand, performing the count in this way, we are going to inevitably miss voids, which instead indirectly retain some information about the clustering properties. All the needed data are provided by the catalogue in [74]: the volume of the cells or, equivalently, the radius of the clusters, r 200 in Mpc, is defined as the radius within which the mean density of a cluster is 200 times the critical density of the Universe at the same redshift and N 200 is the number of member galaxy candidates within r 200 of each cluster. The only caveat to take in mind is that the radii of the cells we are considering, i.e. the r 200 , are much smaller than the scales where the quasi-equilibrium clustering should be more effective [76, 77]; we will discuss this point when presenting our results.
We have divided the full catalogue in groups by redshift, with bin widths of ∆z = 0.05, and by radius, with bin widths of ∆r 200 = 0.1 Mpc, and we have eventually selected only the groups with a sufficiently large number of clusters to enable a strong and reliable statisti-cal analysis. The final groups with which we will work will have radii in the three ranges, 0.8 < r 200 < 0.9, 0.9 < r 200 < 1.0 and 1.0 < r 200 < 1.1 Mpc, and redshifts in the interval 0.05 ≤ z ≤ 0.65. We need to remind here that the catalogue is complete only up to a redshift z ∼ 0.42, for clusters with an estimated mass M 200 > 1 · 10 14 M ⊙ ; at higher redshift, a bias toward larger clusters (larger masses) and possibly to higher counts-in-cell is possible [74].
In Table 1, we show the results derived from using Eq. (29) with the chosen data. The parameters involved in the analysis of Eq. (29) have been rewritten as: -R = r 200 , the radius of the cell (cluster), which enters Eq. (29) through Eq. (17). We have assumed it to be constant within the three bins we have identified; in the following, it will act as a scaling factor for some parameters; -N = N 200 , the number of galaxies within each cell (cluster); -the clustering parameter b, defined as and which enters Eq. (29) through: with α given by Eq. (17); -the dimensionless softening parameter,ǫ = ǫ/R; -the dimensionless nonlocal characteristic length,λ = λ/R.
The only parameter which fully characterizes the nonlocal gravity model we are focusing on is thusλ. All other parameters are standard in the sense that they would appear also whether standard GR would be considered. In Table 1 we also report the number of clusters in each redshift bin, N cl , although this is not a fitting parameter.
When it comes to fit the data, we have defined and compared results from two different statistical tools. We have performed a least-square minimization using the χ 2 defined as where F obs (N i ) is the real counts-in-cell distribution extracted from the data, and F theo (N i ) is the theoretical counts-in-cell calculated from the nonlocal gravity model by using Eq. (29). The index i derives from the fact that in each radius bin, the clusters have a variable amount of galaxies within them, ranging from some minimal number N min to a maximal one N max . Thus, the index i selects the finite natural values of this range, where we evaluate the distribution function. We have also defined a second χ 2 as where the σ i are the errors on F obs (N i ) which we have derived from a jackknife-like procedure as exposed below.
For each redshift and radius bin: 1. we cut a variable fraction F (randomly selected in the range [10, 90]%) from the total bin population; 2. for each cut sample we derive the counts-in-cell distribution function ∼ 50 times; 3. we thus obtain a distribution of F obs (N i ) from the ∼ 50 sets of N i ; 4. from each distribution, which looks very close to a standard Gaussian, we derive the standard deviation and assume such value as the error σ i .
The data points and the errors on the distribution function F obs (N i ) which we finally obtain are shown respectively as black dots and bars in Figs. (5) -(6) -(7). The best fitting F (N ) distributions obtained from the minimization of χ 2 ls are shown in red; those ones derived from χ 2 jk are in green.

Clusters
Concerning the count-in-cells where the counted objects within volumes are the clusters of galaxies, we follow a procedure which is very similar to what is described in [76]. In order to be as much as clear as possible, we are going to enumerate all the steps in the following list.
The initial set up consists in the coming steps: and Ω m = 0.3121, the volume is ≈ 6.5 · 10 9 Mpc 3 . The number of clusters from the SDSS catalog falling into each group are, respectively, 37686, 24208 and 20634 (the total sums up to 82528, because we are taking only clusters with z ≤ 0.42, for which the catalogue is complete); (3a) all clusters in the catalog are provided with redshift z and celestial coordinates, the J2000.0 right ascension β and declination δ. These quantities are converted into physical Cartesian coordinates using the standard formulae, x = D C (z) cos δ cos β , where D C is the comoving distance calculated assuming the same fiducial cosmology as in the previous step. In Fig. 4 we show both the angular projection of the SDSS survey area (top panel), and the "physical length" distribution of the catalogue, with the three bins shown in different colors (bottom panel); (4a) we define a Cartesian grid covering the SDSS survey area. Each point in this grid will be the center of a cell which will be used for the count-in-cell procedure. In order to guarantee a full coverage of the area, considering that the volumes are spheres, the distance between each point (i.e. each center of each cell) is √ 2R; (4a) the SDSS survey does not cover the sky in a uniform way; thus, we are not going to include in our analysis all the points from the previously defined Cartesian grid. Actually, converting the Cartesian coordinates of each point back to celestial coordinates, we only select those ones which fall in the range defined by the initial cluster catalog. This step will acquire more significance as further steps in our procedure will be highlighted in the next points.
Having defined this initial set up, we have performed the following operations on/with it: (1b) we randomize the grid defined in the previous point (3a), by shifting its origin by a random number in the interval [0, √ 2R]. Actually, we apply this shift independently to each Cartesian axis of the grid; (2b) after each shift, we only retain cells which satisfy the condition described in the previous point (5a). This assures us that any cell which should fall out of the SDSS survey will not be considered in the following analysis; (3b) additionally, we apply a jackknife cut of 10% of cells selected randomly; (4b) the previous steps, After all these steps are performed, we eventually have at our disposal the grid of centers of cells/volumes, for each cell size R provided at point (1a), and divided in the three subgroups following the redshift criterion described at point (2a), which will serve to perform the count-in-cell with the cluster catalog data. Basically, for each point in the grid, we will retain only the clusters whose distance from a given grid point (i.e. cell center) will be lower than the size R taken into account. Because of point (4b) above, we basically have 100 different grid sets to compare with the original catalogue, namely, we will have 100 F (N ) spanning the ranges [N min , N max ] where, as before, N min and N max are the minimum and maximum numbers of objects (clusters) within each cell/volume. The median of such ensemble of 100 values for F (N ) will serve as observational estimation for this quantity in the χ ls analysis, as defined in Eq. (39). Instead, the median and the standard deviation, assumed as error σ i , will be used when performing the χ jk analysis, as defined in Eq. (40).

Results
The complete outcomes of our statistical analysis are presented in Table 1 for galaxies and Table 2 for clusters of galaxies. As a first comment, valid for both cases, let us notice that the softening parameterǫ is not present in the tables because it is basically unconstrained, exhibiting a uniform distribution all over the given prior range [0, 1]. Thus, we can infer that its role in the fit is marginal.

Galaxies
For what concerns galaxies, from Table 1 we can see how the parameterN is very well constrained, and there is no difference between using χ 2 ls or χ 2 jk . This statement is valid whatever radius R and redshift bin is considered. Actually, we must remember that the catalog is complete only up to z = 0.42, after which we do have a bias toward smaller clusters, namely, clusters with a lowerN . And this is effectively seen in the tables, wherē N rises slightly with redshift, becomes practically constant, and then starts to decrease for z > 0. 45.
The same consideration is basically true also for the clustering parameter b, even if in this case we can see how there is a trend with varying R, while no peculiar change with redshift z can be detected. In general, this parameter lies in the range [0.3, 0.4] for all cases; however, while for the smaller R we can see how it is very well constrained, for growing R, we are able to put only upper limits on it, which means that also small values of b are compatible with data. As described in [77], we know that smaller values of b correspond to clusters which have not yet fully virialized. It is worth noticing that such trend is on the opposite side with respect to [19,20,77] where smaller cells volumes correspond to smaller values of b; although we must also stress that in [19,20,77], the considered physical volumes are much larger than those we have defined here as well as the redshift bins, and the number of objects in each bin is much lower. In our case the bins are much finer, both in radius and in redshift.
For what concerns the only fully characterizing parameter of the nonlocal gravitational potential, i.e.λ, we can see how it is remarkably very well constrained in a specific and well defined range. Statistically speaking, we can see how it approximately lies in the range [0.4, 0.5], which, in physical units, corresponds roughly to [0.32, 0.55] Mpc. We can detect a slight decrease for varying redshift, with smaller values at larger z; while it seems to be more pronounced a trend with varying R, with largerλ for larger clusters. It is quite striking to have such values for this parameter, which are actually totally consistent with the size of gravitational structures/volumes under scrutiny. At the same time, they are also very different from a stable GR limit, which would be achieved when λ → ∞. Actually, from Eq. 7, one can see how it would be possible to recover an unstable GR-like limit even for λ → r ij .
Finally, from Figs. 5 -6 -7, we can also visually check the quality of the fits: the larger is the cell volume, the better is the fit.
We must stress here an important point: these results just show that a logarithmic correction to the Newtonian gravity might be consistent with the clustering of galaxies. But we cannot state at all that this correction is consistent with the full clustering process at these same scales, because, from the catalogue we have used, we have no information on voids, which should be taken into account. This is a further check which should be performed in the future.

Clusters
What we can check is how the full clustering process on larger scales behaves. As we have stated in the previous sections, if we move to consider the clustering of clusters of galaxies, the catalogue in hand does allow us to properly take into account voids. Although, in this case, results are a bit more fuzzy.
First of all, we focus onN . We can see here a trend both in redshift and in volume size R:N decreases with redshift and increases with R. The latter correspondence is of course expected, as large cells can contain more clusters; the former one is also somehow expected, because we might predict an intrinsic lack of clusters at higher redshifts, as they are caught in an epoch of formation, while a larger number of them is observable at smaller redshifts. Notice that our bins all have the same comoving volumes, so that any geometry-related issue should not be effective in this case.
The behavior of the clustering parameter b is now fully consistent with literature [19,20,77], with smaller values for smaller volumes, as well as for higher redshifts.
The nonlocal characteristic lengthλ, instead, in this case is quite unconstrained. We are unable to put strong boundaries on it and in Table 2 we can only assess upper limits on it. Although we must warn that such results are statistically weak, with the posteriors being quite irregular in most of the cases. This result can only be interpreted as a substantial negligible role of the parameter λ in the fitting procedure. Thus, we should conclude that on such large scales there is no evidence for a detectable logarithmic correction to the standard Newtonian potential.

Discussion and Conclusions
Nonlocality effects can occur at large distances as the consequence of modified gravitational potential. Specifically, nonlocal gravity, besides the effects at cosmological scale as a possible engine for accelerated expansion [39,42], is relevant also for large scale formation. In particular, it is capable of triggering the clustering of galaxies and clusters of galaxies.
In view of this statement, we have analyzed the effects of nonlocal gravitational potential on the structure formation comparing the model with observational data. This has been done by treating both galaxies and clusters of galaxies (in two distinct analysis) as points in a statistical mechanical system, and then analyzing their clustering using a gravitational partition function. This function has been modified according to the nonlocal gravitational potential and its influence on the properties of the gravitational clustering has been scrutinized. Finally, the model has also been compared with observational data, and we have demonstrated that a logarithmic correction to standard gravity seems to be consistent with the observations, at least with those concerning the galactic scale, i.e. within volumes corresponding to clusters of galaxies. Thus, it is possible to conclude that the nonlocal gravitational effects might be revealed at such scales. On the other hand, we must stress that the clustering of larger volumes of space, in which the points are the clusters of galaxies, does not show any evidence in favor of a such modification to gravity so that any inference should be suspended until further (maybe more precise) probes are checked.
For a more general discussion, it has be noted that other modifications of gravitational potential, explaining dark matter or dark energy effects, have been studied considering modifications of the gravitational partition function. For example, it is well-known that f (R) gravity modifies the large distance behavior of gravitational potential [78][79][80][81]. Thus, the clustering of galaxies, interacting through effective potentials derived from f (R) gravity, has also been studied using a modified gravitational partition function, and it has been observed that also these models are consistent with ob-servations for large samples of galaxies and clusters [82,83], although the corresponding statistical inferenec is a bit weaker that the nonlocal scenario we have considered here. It is worth noticing that the clustering of galaxies in f (R) gravity is also consistent with the constraints coming from the Planck data [84].
The MOND gravitational partition function has also been used to study the thermodynamics and the clustering properties of systems of galaxies [13,85].
Furthermore, the gravitational phase transition can also be analyzed using the gravitational partition function for systems of galaxies [14]. In fact, a first order phase transition occurs due to the clustering of galaxies from a homogeneous phase. In this framework, it is possible to take into account the Yang-Lee theory for systems of galaxies, and use the complex fugacity to analyze the phase transition [86]. A forthcoming step will be to extend the Yang-Lee theory for nonlocal gravitational potential, and analyze the phase transition for the related gravitational partition function.
In particular, it has to be noted that the clustering of galaxies has also been investigated using the cosmic energy equation [87]. As the cosmic energy equation is derived by approximating galaxies as points in a statistical mechanical system, a softening parameter can modify the cosmic energy equation [88]. It has been demonstrated that a large scale modification of the gravitational potential also modifies the cosmic energy equation [89]. Table 1 Galaxies: results from fitting Eq. (29) with data from the cluster catalog from [74]. Data and fits are divided in redshift bins (first column) and cell radius bin (as indicated by R intervals). N cl is not a fitting parameter, but the number of clusters in each redshift and radius bin. All the fitting parameters are described in the text. 0.8 < R < 0.9 Mpc 0.9 < R < 1.0 Mpc 1.0 < R < 1.   (29) with data from the cluster catalog from [74]. Data and fits are divided in redshift bins (first column) and cell radius bin (as indicated by R intervals). N cl is not a fitting parameter, but the number of clusters in each redshift and radius bin. All the fitting parameters are described in the text.   Fig. 8 Comparison between data and theoretical expectation Eq. (29) for cell size bin 1.0 < R 1 < 1.1 M pc. Black points: data; black bars: jackknife-like observational errors. Solid red line: best fit from minimization of χ 2 ls ; solid green line: best fit from minimization of χ 2 jk .