Abstract
We study the statistical mechanics of binary systems under the gravitational interaction of the Modified Newtonian Dynamics (MOND) in three-dimensional space. Considering the binary systems in the microcanonical and canonical ensembles, we show that in the microcanonical systems, unlike the Newtonian gravity, there is a sharp phase transition, with a high-temperature homogeneous phase and a low-temperature clumped binary one. Defining an order parameter in the canonical systems, we find a smoother phase transition and identify the corresponding critical temperature in terms of the physical parameters of the binary system.