Reply to the recent comment by I.Ispolatov and M.Karttunen, cond-mat/0303564

We study the classical Antonov problem (of retrieving the statistical equilibrium properties of a self-gravitating gas of classical particles obeying Boltzmann statistics in space and confined in a spherical box) for a rotating system. It is shown that a critical angular momentum $\lambda_c$ (or, in the canonical language, a critical angular velocity $\omega_c$) exists, such that for $\lambda<\lambda_c$ the system's behaviour is qualitatively similar to that of a non-rotating gas, with a high energy disordered phase and a low energy collapsed phase ending with Antonov's limit, below which there is no equilibrium state. For $\lambda>\lambda_c$, instead, the low-energy phase is characterized by the formation of two dense clusters (a ``binary star''). Remarkably, no Antonov limit is found for $\lambda>\lambda_c$. The thermodynamics of the system (phase diagram, caloric curves, local stability) is analyzed and compared with the recently-obtained picture emerging from a different type of statistics which forbids particle overlapping.

We complete the analysis of the equilibrium shapes of rotating self-gravitating gases initiated in Phys. Rev. Lett. 89 031101 (2002) (cond-mat/0202140) by studying the formation and thermodynamics of asymmetric binaries and by introducing a family of order parameters to classify the possible equilibrium states.

We investigate the statistical equilibrium properties of a system of classical particles interacting via Newtonian gravity, enclosed in a three-dimensional spherical volume. Within a mean-field approximation, we derive an equation for the density profiles maximizing the microcanonical entropy and solve it numerically. At low angular momenta, i.e. for a slowly rotating system, the well-known gravitational collapse ``transition'' is recovered. At higher angular momenta, instead, rotational symmetry can spontaneously break down giving rise to more complex equilibrium configurations, such as double-clusters (``double stars''). We analyze the thermodynamics of the system and the stability of the different equilibrium configurations against rotational symmetry breaking, and provide the global phase diagram.

We derive the global phase diagram of a self-gravitating $N$-body system enclosed in a finite three-dimensional spherical volume $V$ as a function of total energy and angular momentum, employing a microcanonical mean-field approach. At low angular momenta (i.e. for slowly rotating systems) the known collapse from a gas cloud to a single dense cluster is recovered. At high angular momenta, instead, rotational symmetry can be spontaneously broken and rotationally asymmetric structures (double clusters) appear.