We exploit a new theory of gravity proposed by Finzi, which gives stronger interaction at large
scales, to study the thermodynamic description of galaxy clusters. We employ a statistical model
to deduce various thermodynamics equations of state. In addition, we analyze the behavior
of clustering parameter in comparison to its standard (Newtonian) counterpart. The general
distribution function and its behavior with varying strength of clustering parameter are also
studied. The possibility of phase transition is investigated and it is observed that a phase
transition is possible though hierarchically. We also analyze the model by comparing the results
with data available through SDSS-III, and obtain the parameters involved.