Our Galaxy and the nearby Andromeda Galaxy (M31) form a bound system, even though the relative velocity vector of M31 is currently not well constrained. Their orbital motion is highly dependent on the initial conditions, but all the reliable scenarios imply a first close approach in the next 3-5 Gyrs. In our study, we simulate this interaction via direct N-body integration, using the HiGPUs code. Our aim is to investigate the dependence of the time of the merger on the physical and dynamical properties of the system. Finally, we study the dynamical evolution of the two Supermassive Black Holes placed in the two galactic centers, with the future aim to achieve a proper resolution to follow their motion until they form a tight binary system.

According to our current knowledge about physical and dynamical properties of the Milky Way-M31 system, it seems likely that these two galaxies will collide and eventually merge in a time very sensitive to initial conditions. Using the HiGPUs code, we have performed several numerical simulations to study the dynamics of the system, trying to define the role of indeterminacy in the present day observed relative velocities of the two galaxies and the time of the merger. At the same time, we have followed the dynamics of the two massive black holes sitting in the galactic centers, to check (within the space and time resolution limits of our simulation) their relative motion upon the completion of the galaxies merger process.