Star accretion onto supermassive black holes in axisymmetric galactic nuclei

Tidal Disruption (TD) of stars by supermassive central black holes from dense rotating star clusters is modeled by high-accuracy direct N-body simulation. We study the time evolution of the stellar tidal disruption rate and the origin of tidally disrupted stars. Compared with that in spherical systems, we found a higher TD rate in axisymmetric systems. The enhancement can be explained by an enlarged loss-cone in phase space which is raised from the fact that total angular momentum J is not conserved. As in the case of spherical systems, the distribution of the last apocenter distance of tidally accreted stars peaks at the classical critical radius. However, the angular distribution of the origin of the accreted stars reveals bimodal features. We show that the bimodal structure can be explained by the presence of two families of regular orbits, namely short axis tube and saucer orbits.