Milky Way Subsystems from Globular Clusters Kinematics Using Gaia DR2 and HST Data A. T. Bajkova1, G. Carraro2, V. I. Korchagin3, N. O. Budanova3, V. V. Bobylev1 1Pulkovo Astronomical Observatory, St.-Petersburg, Russia 2Dipartimento di Fisica e Astronomia, Universit´adi Padova, Vicolo Osservatorio 3 I-35122, Padova, Italy 3Southern Federal University, Rostov-on-Don, Russia Abstract We employ Gaia DR2 proper motions for 151 Milky Way globular clusters from Vasiliev (2019) in tandem with distances and line-of-sight velocities to derive their kinematical properties. To assign clusters to the Milky Way thick disk, bulge, and halo we follow the approach of Posti et al. (2018) who distinguished among different Galactic stellar components using stars’s orbits. In particular, we use the ratio Lz/e, the Z projection of the angular momentum to the eccentricity, as population tracer, which we complement with chemical abundances extracted from the literature and Monte-Carlo simulations. We find that 20 globular clusters belong to the bar/bulge of the Milky Way, 35 exhibit disk properties, and 96 are members of the halo. Moreover, we find that halo globular clusters have close to zero rotational velocity with average value < Θ > =1± 4 km s−1. On the other hand, the sample of clusters that belong to the thick disk possesses a significant rotation with average rotational velocity 179 ± 6 kms−1. The twenty globular clusters orbiting within the bar/bulge region of the Milky Way galaxy have average rotational velocity of 49 ± 11 km s−1. Key words: Milky Way, globular clusters (GCs),thin disk, thick disk, stellar halo 1 INTRODUCTION The study of kinematical properties of the Milky Way globular clusters (GCs) has a long tradition. One of the most important step in understanding the nature of the globular cluster system was done by Zinn (1993) who suggested a cluster classification based on the horizontal branch morphology and chemical abundances. Zinn (1993) found that the metal rich globular clusters ([F e/H] > −0.8) are confined to the bar/bulge and to the disk of the Galaxy, while the metal-poor globular clusters ([F e/H] ≤ −0.8) are generally located in the Galactic halo. Later on, Mackey & Gilmore (2004) arXiv:2004.13597v1 [astro-ph.GA] 28 Apr 2020 used Zinn’s (1993) classification and compared the cluster properties with those in the Large and the Small Magellanic Clouds concluding that all young globular clusters and 15 -17 percent of the old ones are of external origin. More recent progresses in understanding the structural properties of the Milky Way galaxy come from the Gaia space mission. The recently released Gaia DR2 catalog Gaia Collaboration (Helmi et al., 2018) allowed to make significant progress in understanding the nature and the origin of the Milky Way globular cluster system. Gaia collaboration (Helmi et al., 2018) determined proper motions of seventy five globular clusters. In the subsequent paper Posti & Helmi (2019) used the Gaia Data Release 2 proper motions catalog of 75 globular clusters and measurements of the proper motions of another 20 distant clusters obtained with the Hubble Space Telescope, and determined the kinematical properties of the globular clusters in the Milky Way galaxy. Posti & Helmi (2019) followed up an earlier work by Binney & Wong (2017), and described the population of the Milky Way globular clusters by the two-component distribution function that describes a flat, rotating disk and a spherical extended halo. Compared to the paper of Binney & Wong (2017), 1 Posti & Helmi (2019) used larger data set, and varied the mass of the dark matter halo and its shape. Vasiliev (2019) published the proper motion catalog of 151 Milky Way globular clusters with estimated uncertainties smaller than 0.1 mas yr−1 for most objects. Catalog of Vasiliev (2019) complements the measurements provided by the Gaia collaboration (Helmi et al.,2018) and HST-based determinations Sohn et al. (2018), and is in excellent agreement with both. Vasiliev (2019) compared his results with the ground based measurements, and found that they are not only much less precise than the space-based ones, but often deviate from the space- based determinations by a far larger amount than implied by the quoted uncertainties. Combining proper motion measurements with the distance and and line-of-sight velocity measurements from the literature, Vasiliev (2019) found that within inner 10 kpc the population of the globular clusters has a mean rotational velocity of 50 –80 km s−1 with nearly isotropic velocity dispersion about 100 – 120 km s−1. Baumgardt et al. (2019) re-derived the mean proper motions and space velocities for 154 Milky Way globular clusters using a combination of the Gaia DR2 proper motions and ground-based line- of-sight velocities of individual member stars. Mean proper motions of Baumgardt et al. (2019) are in a good agreement with measurements of Gaia Collaboration (Helmi et al.,2018) and with Vasiliev’s (2019) measurements, and also are in a good agreement with the proper motions derived by Sohn et al. (2018) from HST data if the systematic error of 0.10 mas yr−1 in the Gaia proper motions is taken into account. Good quality kinematical data of nearly all known globular clusters of the Milky Way galaxy allow one to address the problem of assigning them to different Milky Way components. A number of recent papers have been devoted to the problem of classification of GSc and the understanding the their origin (Myeong et al. 2019; Massari et al., 2019; Piatti, 2019; Horta et al., 2020; Forbes, 2020; Perez-Villegas et al., 2020). The investigation of Massari et al. (2019) is particularly relevant here. This work essentially provides a solution to the problems of the separation of GCs into the subsystems of the Galaxy, and, most importantly, answers the question about the origin of the halo GCs. To separate the GCs between internal (in-situ) and external (ex-situ), the authors use various dynamical properties, namely, orbital parameters like the apocenter (apo), maximum height from the disc (Zmax), and circ the orbital circularity parameter Circ. The latter parameter is defined as Circ = Lz/Lz where Lz circ is the Z-component of the angular momentum of the cluster, and Lz is the angular momentum of the cluster placed on the circular orbit with the same total energy. Combination of these three parameters with the age-metallicity relation allowed to solve the problem of the separation. In this paper we readdress the problem of defining globular clusters membership to Milky Way different components, – the bar/bulge, the disk, and the halo only on the basis of the dynamical properties of the clusters. Basing on the clusters orbit parameters such as apo and Zmax, we first select those clusters that belong to the bar/bulge of the Milky Way and to its halo similar to how Massari et al. (2019) did. To separate the globular clusters within the thick disk region of the Milky Way galaxy into the disk and halo fractions, we use instead of the circular parameter Circ the distribution of the clusters on the ratio of the Z projection of the angular momentum of the clusters to their eccentricity e, Lz/e, and demonstrate that this quantity allows us to distinguish clusters that belong to the halo and to the thick disk effectively. In detail, the goal of this paper is to study the properties of the globular cluster subsystems based on the dynamical properties of the their orbital parameters and to test the proposed criterion using accepted models of the Galactic potential. The paper is organised as follows. In Section 2 we describe available kinematical data, and derive spatial velocity components of the globular clusters in the Galacto-centric coordinate system. In this Section we also describe the Galactic potential models (both axisymmetric and barred ones), and perform orbit integration of the clusters. Section 3 is devoted to a new method for separat- ing globular clusters into Galactic subsystems, an analysis of the kinematic characteristics of the 2 resulting subsystems, and a comparison of our results with the results of Massari et al. (2019). Section 4 is devoted to testing the separation algorithm against different potential models. Section 5 discusses results of the separation in a barred potential with different bar rotational velocities. In Section 6 we add metallicity properties of the globular clusters to support our findings. Finally, section 7 summarises our main results. The Appendix contains several tables with clusters’ orbital parameters. 2 Galactic potential model and Globular Clusters orbits Globular clusters are distributed over a wide range of distances. To integrate their orbits one needs therefore a model of the Milky Way potential reliable in a wide range of distances from the Galactic center out to the very distant halo. Bajkova & Bobylev (2016, 2017) specified the parameters of a few models of the Milky Way galaxy. They assumed that the potentials of the disk and of the bulge are described by the Miyamoto & Nagai (1975) and the Plummer models, respectively, while for the halo potential they adopted the models of Allen & Santillan (1991), Wilkinson & Evans (1999), Navarro, Frenk & White (1997), the spherical logarithmic potential of Binney (Binney & Tremaine 1987), the Plummer sphere, and the Hernquist (1990) potential. To fix the parameters of the models, Bajkova and Bobylev (2016, 2017) used observational data, covering galactocentric distances up to 200 kpc. Within 20 kpc, Bajkova and Bobylev additionally constrained the models by the line-of-sight velocities of the molecular clouds in their tangential points and by the kinematical data for 130 masers with known trigonometric parallaxes. At larger distances, Bajkova & Bobylev (2016, 2017) took into account the rotational velocities of the objects taken from the catalogue of Bhattacharjee et al.