Abstract
We numerically investigate the motion of stars on the meridional plane of an axially symmetric disk galaxy model, containing a central supermassive black hole, represented by the Paczyski-Wiita potential. By using this pseudo-Newtonian potential we can replicate important relativistic properties such as the existence of the Schwarzschild radius. After classifying extensive samples of initial conditions of trajectories, we managed to distinguish between collisional, ordered, and chaotic motion. Besides all starting conditions of regular orbits were further classified into families of regular orbits. Our results are presented via color-coded basin diagrams on several types of two-dimensional planes. Our analysis reveals that both the mass of the black hole (in direct relation with the Schwarzschild radius) as well as angular momentum play an important role in the character of the orbits of stars. More specifically, the trajectories of low angular momentum stars are highly affected by the mass of the black hole, while high angular momentum stars seem to be unaffected by the central black hole. A comparison with previous related outcomes, using Newtonian potentials for the central region of the galaxy, is also made.