One-dimensional toy model of globular clusters

We introduce a one-dimensional toy model of globular clusters. The model is a version of the well-known gravitational sheets system, where we take additionally into account mass and energy loss by evaporation of stars at the boundaries. Numerical integration by the "exact" event-driven dynamics is performed, for initial uniform density and Gaussian random velocities. Two distinct quasi-stationary asymptotic regimes are attained, depending on the initial energy of the system. We guess the forms of the density and velocity profiles which fit numerical data extremely well and allow to perform an independent calculation of the self-consistent gravitational potential. Some power-laws for the asymptotic number of stars and for the collision times are suggested.