Milky Way globular clusters on cosmological timescales

III. Interaction rates

¹ Main Astronomical Observatory, National Academy of Sciences of Ukraine, 27 Akademika Zabolotnoho St, 03143 Kyiv, Ukraine
e-mail: marina@mao.kiev.ua
² Fesenkov Astrophysical Institute, Observatory 23, 050020 Almaty, Kazakhstan
³ Astronomisches Rechen-Institut, Zentrum für Astronomie, University of Heidelberg, Mönchhofstrasse 12-14, 69120 Heidelberg, Germany
⁴ Nicolaus Copernicus Astronomical Centre, Polish Academy of Sciences, ul. Bartycka 18, 00-716 Warsaw, Poland
⁵ Konkoly Observatory, Research Centre for Astronomy and Earth Sciences, Eötvös Loránd Research Network (ELKH), MTA Centre of Excellence, Konkoly Thege Miklós út 15-17, 1121 Budapest, Hungary
⁶ Al-Farabi Kazakh National University, 71 al-Farabi Ave., 050040 Almaty, Kazakhstan
⁷ Energetic Cosmos Laboratory, Nazarbayev University, 53 Kabanbay Batyr Ave., 010000 Astana, Kazakhstan

Received: 31 March 2023
Accepted: 7 August 2023

Abstract

Aims. We carry out the self-consistent dynamic evolution of the orbital structure of Milky Way globular clusters. This allows us to estimate possible and probable close passages and even collisions of the clusters with each other.

Methods. We reproduced the orbits of 147 globular clusters in 10 Gyr lookback time using our own high-order N-body parallel dynamic φ-GPU code. The initial conditions (three coordinates and three velocities for the present time) were derived from the Gaia DR3 catalogue. For each of the GCs, 1000 initial conditions were additionally generated, taking the Gaia measurement errors into account. The galaxy is represented by five external potentials from the IllustrisTNG-100, whose masses and sizes of the disk and halo components are similar to the physical values of the Milky Way at present.

Results. We present a statistical analysis of the cumulative close passages rate: About ten close passages with relative distances shorter than 50 pc for every billion years for each of the five external potentials. We present the 22 most reliable collision pairs with a good probability. As an example: Terzan 4 vs. Terzan 2 (49%), Terzan 4 vs. NGC 6624 (44%), Terzan 4 vs. Terzan 5 (40%), Terzan 4 vs. NGC 6440 (40%), and Terzan 4 vs. Liller 1 (42%). The most active globular cluster in the collision sense is Terzan 4, which has 5.65 collision events on average (averaged over all individual 1000 initial condition realisations). Most collisions are located inside the Galactic disk and form two ring-like structures. The first ring-like structure has the highest collision number density at ∼1 kpc, and the second sturcture has a maximum at ∼2 kpc.

Conclusions. Based on our numerical simulations, we can conclude that the few dozen Milky Way globular clusters probably undergo some close encounters and even possible collisions during their lifetimes, which can significantly affect their individual dynamical evolution and possibly even their stellar content.