An improved method for constructing models of self-gravitating tori   around black holes
  General-relativistic models of self-gravitating tori around black holes are constructed with a self-consistent-field method in compactified coordinates. The numerical code is highly accurate and robust, allowing for the construction of models that exactly fill their Roche lobe, when a cusp exists. As a first application, we focus on self-consistent models with cusp, having different values of constant specific angular momentum. Scaling all results with the mass of the black hole, we find evidence that models with constant specific angular momentum that can fill their Roche lobe are still limited by $l<4M_{\rm BH}$ (as is the case for models constructed in a fixed background metric) even for heavy tori. 
Nikolaos Stergioulas
gr-qc