In this paper, we analyze the clustering of galaxies using a modified theory of gravity, in which the field content of general relativity has been be increased. This increasing in the field content of general relativity changes the large distance behavior of the theory, and in weak field approximation, it will also modify the large distance behavior of Newtonian potential. So, we will analyzing the clustering of multi-component system of galaxies interacting through this modified Newtonian potential. We will obtain the partition function for this multi-component system, and study the thermodynamics of this system. So, we will analyze the effects of the large distance modification to the Newtonian potential on Helmholtz free energy, internal energy, entropy, pressure and chemical potential of this system. We obtain also the modified distribution function and the modified clustering parameter for this system, and hence observe the effect of large distance modification of Newtonian potential on clustering of galaxies.

This paper is an attempt to study the thermodynamics of the structure formation in the large scale universe in the non local gravity using Boltzmann statistics and the Tsallis statistics. The partition function is obtained in both the approaches and the corresponding thermodynamics properties are evaluated. The important thing about the paper is that we surprisingly get the divergence free integrals and thus stress upon the fact that the nonlocal gravity is the singularity free model of gravity.

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.