Solvable model of a self-gravitating system

We introduce and discuss an effective model of a self-gravitating system whose equilibrium thermodynamics can be solved in both the microcanonical and the canonical ensemble, up to a maximization with respect to a single variable. Such a model can be derived from a model of self-gravitating particles confined on a ring, referred to as the self-gravitating ring (SGR) model, allowing a quantitative comparison between the thermodynamics of the two models. Despite the rather crude approximations involved in its derivation, the effective model compares quite well with the SGR model. Moreover, we discuss the relation between the effective model presented here and another model introduced by Thirring forty years ago. The two models are very similar and can be considered as examples of a class of minimal models of self-gravitating systems.Comment: 21 pages, 6 figures; submitted to JSTAT for the special issue on long-range interaction

Thermodynamics of the self-gravitating ring model

We present the phase diagram, in both the microcanonical and the canonical ensemble, of the selfgravitating- ring sSGRd model, which describes the motion of equal point masses constrained on a ring and subject to 3D gravitational attraction. If the interaction is regularized at short distances by the introduction of a softening parameter, a global entropy maximum always exists, and thermodynamics is well defined in the mean-field limit. However, ensembles are not equivalent and a phase of negative specific heat in the microcanonical ensemble appears in a wide intermediate energy region, if the softening parameter is small enough. The phase transition changes from second to first order at a tricritical point, whose location is not the same in the two ensembles. All these features make of the SGR model the best prototype of a self-gravitating system in one dimension. In order to obtain the stable stationary mass distribution, we apply an iterative method, inspired by a previous one used in 2D turbulence, which ensures entropy increase and, hence, convergence towards an equilibrium state.othe

Thermodynamics of rotating self-gravitating systems

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.Comment: 12 pages, 9 figure

Phase transitions in simplified models with long-range interactions

We study the origin of phase transitions in some simplified models with long range interactions. For the ring model, we show that a possible new phase transition predicted in a recent paper by Nardini and Casetti from an energy landscape analysis does not occur. Instead of such phase transitions we observe a sharp, although without any non-analiticity, change from a core-halo to an only core configuration in the spatial distribution functions for low energies. By introducing a new class of solvable simplified models without any critical points in the potential energy, we show that a similar behaviour to the ring model is obtained, with a first order phase transition from an almost homogeneous high energy phase to a clustered phase, and the same core-halo to core configuration transition at lower energies. We discuss the origin of these features of the simplified models, and show that the first order phase transition comes from the maximization of the entropy of the system as a function of energy an an order parameter, as previously discussed by Kastner, which seems to be the main mechanism causing phase transitions in long-range interacting systems

Solving the Vlasov equation for one-dimensional models with long range interactions on a GPU

We present a GPU parallel implementation of the numeric integration of the Vlasov equation in one spatial dimension based on a second order time-split algorithm with a local modified cubic-spline interpolation. We apply our approach to three different systems with long-range interactions: the Hamiltonian Mean Field, Ring and the self-gravitating sheet models. Speedups and accuracy for each model and different grid resolutions are presented