Chaotic orbits in a 3D galactic dynamical model with a double nucleus

A 3D-dynamical model is constructed for the study of motion in the central regions of a disk galaxy with a double nucleus. Using the results of the 2D model, we find the regions of initial conditions in the (x,px,z,py)=EJ,(y=pz=0) phase space producing regular or chaotic orbits. The majority of stars are on chaotic orbits. All chaotic orbits come arbitrarily close to one or to both nuclei. Regular orbits are those starting near the stable periodic orbits of the 2D-system with small values of z. All regular orbits circulate around only one of the two nuclei.