The construction of non-spherical models of quasi-relaxed stellar systems

Spherical models of collisionless but quasi-relaxed stellar systems have long been studied as a natural framework for the description of globular clusters. Here we consider the construction of self-consistent models under the same physical conditions, but including explicitly the ingredients that lead to departures from spherical symmetry. In particular, we focus on the effects of the tidal field associated with the hosting galaxy. We then take a stellar system on a circular orbit inside a galaxy represented as a "frozen" external field. The equilibrium distribution function is obtained from the one describing the spherical case by replacing the energy integral with the relevant Jacobi integral in the presence of the external tidal field. Then the construction of the model requires the investigation of a singular perturbation problem for an elliptic partial differential equation with a free boundary, for which we provide a method of solution to any desired order, with explicit solutions to two orders. We outline the relevant parameter space, thus opening the way to a systematic study of the properties of a two-parameter family of physically justified non-spherical models of quasi-relaxed stellar systems. The general method developed here can also be used to construct models for which the non-spherical shape is due to internal rotation. Eventually, the models will be a useful tool to investigate whether the shapes of globular clusters are primarily determined by internal rotation, by external tides, or by pressure anisotropy.