Based on thermodynamics, we discuss the galactic clustering of expanding Universe by assuming the gravitational interaction through the modified Newton's potential given by $f(R)$ gravity. We compute the corrected $N$-particle partition function analytically. The corrected partition function leads to more exact equations of states of the system. By assuming that system follows quasi-equilibrium, we derive the exact distribution function which exhibits the $f(R)$ correction. Moreover, we evaluate the critical temperature and discuss the stability of the system. We observe the effects of correction of $f(R)$ gravity on the power law behavior of particle-particle correlation function also. In order to check feasibility of an $f(R)$ gravity approach to the clustering of galaxies, we compare our results with an observational galaxy cluster catalog.

In this paper, we study thermodynamics of the cluster of galaxies under the effect of dynamical dark energy. We evaluate the configurational integral for interacting system of galaxies in an expanding universe by including the effects produced by the varying $\Lambda$. The gravitational partition function is obtained using this configuration integral. We obtain thermodynamics quantities in canonical ensemble which depend on time and investigate the second law of thermodynamics. We also calculate the distribution function in grand canonical ensemble. The time evolution of the clustering parameter of galaxies is investigated for the time dependent (dynamical) dark energy. We conclude that the second law of thermodynamics is valid for the total system of cluster of galaxies and dynamical dark energy. We calculate correlation function and show that our model is very close to Peebles's power law, in agreement with the N-body simulation. It is observed that thermodynamics quantities depend on the modified clustering parameter for this system of galaxies.

In this paper, we analyse the effect of the expansion of the universe on the clustering of galaxies. We evaluate the configurational integral for interacting system of galaxies in an expanding universe by including effects produced by the cosmological constant. The gravitational partition function is obtained using this configuration integral. Thermodynamic quantities, specifically, Helmholtz free energy, entropy, internal energy, pressure and chemical potential are also derived for this system. It is observed that they depend on the modified clustering parameter for this system of galaxies. It is also demonstrated that these thermodynamical quantities get corrected because of the cosmological constant.