Motion in a double nucleus galactic dynamical model

We construct a 2D-dynamical model in order to study the motion in a galaxy with a double nucleus. Numerical calculations show that the majority of high energy stars, near the double nuclear region, are on chaotic orbits. On the contrary, low energy stars are on chaotic orbits only near massive nuclei while, for less massive nuclei, the motion is regular. Using a semi-numerical approach we explain the chaotic scattering of orbits near each nucleus. The high values of the velocities near the dense nuclei are also explained. Computation of the LCEs indicates a very high degree of chaos. The results derived from our double nucleus dynamical model are compared with models with single ones. Our outcomes are in satisfactory agreement with observational data for the double nucleus system $\mathit{NGC}6240$.