## Using Gaia for studying Milky Way star clusters

## Eugene Vasiliev 1 2 ,

1 Institute of Astronomy, University of Cambridge, UK 2 Lebedev Physical Institute, Moscow, Russia email: eugvas@lpi.ru

Abstract. We review the implications of the Gaia Data Release 2 catalogue for studying the dynamics of Milky Way globular clusters, focusing on two separate topics.

The first one is the analysis of the full 6-dimensional phase-space distribution of the entire population of Milky Way globular clusters: their mean proper motions (PM) can be measured with an exquisite precision (down to 0.05 mas yr -1 , including systematic errors). Using these data, and a suitable ansatz for the steady-state distribution function (DF) of the cluster population, we then determine simultaneously the best-fit parameters of this DF and the total Milky Way potential. We also discuss possible correlated structures in the space of integrals of motion. The second topic addresses the internal dynamics of a few dozen of the closest and richest globular clusters, again using the Gaia PM to measure the velocity dispersion and internal rotation, with a proper treatment of spatially correlated systematic errors. Clear rotation signatures are detected in 10 clusters, and a few more show weaker signatures at a level glyph[greaterorsimilar] 0 05 . masyr -1 . PM dispersion profiles can be reliably measured down to 0.1 mas yr -1 , and agree well with the line-of-sight velocity dispersion profiles from the literature.

Keywords. proper motions - globular clusters: general - Galaxy: kinematics and dynamics

## 1. Introduction

The second data release (DR2) of the Gaia mission (Gaia Collaboration 2018a) provides proper motions (PM) and broad-band photometry for more than a billion stars down to magnitude G = 21, based on 22 months of observations. Despite some limitations (e.g., no special treatment of binary stars, and poor completeness in crowded fields), this dataset is a giant leap forward for studies of Milky Way and its constituents, and has been used in more than 1000 papers over the last year.

The mean PM of roughly half of the Milky Way globular clusters have been calculated in Gaia Collaboration (2018b). Subsequently, mean PM of almost the entire population of the ∼ 150 clusters (Harris 1996, 2010) have been independently measured by Vasiliev (2019a) and Baumgardt et al. (2019). The former study additionally considered the constraints on the Milky Way potential from the dynamics of globular clusters, while the latter also analyzed the internal PM for most of them. Vasiliev (2019b) further extended the analysis of internal kinematics by considering the impact of spatially correlated systematic errors in Gaia astrometry.

## 2. Membership determination

The first step in analysis is the selection of cluster members. We use only the positions and PM of stars in a circular area centered on each cluster, but not their magnitudes or colors, which may be affected by extinction and (in the case of colors) crowding. We

Figure 1. Left panel: correlation function of PM measurement errors as a function of angular separation between two sources.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00009v1-figure-4d7b9c70-d0bd-448b-b2d5-cca2f192ffe0.jpg)

Right panel: example of measurement of internal kinematics for the globular cluster NGC 7078. Shaded regions show the 68% and 95% confidence intervals for the rotation (red) and PM dispersion (blue) profiles, determined in Vasiliev (2019b); points with error bars show measurements from other studies.

select only stars with parallax values consistent with the inverse distance to a cluster at 3 σ level, satisfying certain quality criteria.

Instead of using hard cuts (e.g. removing 3 σ outliers in PM), we employ a probabilistic Gaussian mixture classification method, described in the appendix of Vasiliev (2019a); the code is provided at https://github.com/GalacticDynamics-Oxford/GaiaTools . First, we define a two-component model for the joint distribution of stars in the sky plane and PM space. The density of cluster stars is assumed to follow a particular functional form (Plummer profile) in the sky plane, and their PM have an isotropic Gaussian distribution with a spatially-varying width σ c ( R ), while the distribution of field stars is uniform on the sky plane and is described by a generic 2d Gaussian in the PM space. The fraction of cluster stars, their mean PM and the dispersion σ c (0), the scale radius of the Plummer profile, the mean and covariance of the Gaussian distribution of field stars, are all free parameters in the model. The likelihood of measuring a particular value of PM is given by the sum of the values of cluster and field distribution functions (DFs), convolved with measurement errors (different for each star). The parameters of the model are then optimized to maximize the log-likelihood of the entire dataset, which is a sum of log-likelihoods of individual stars. Finally, after the best-fit parameters have been determined, we evaluate the membership probability for each star.

## 3. Systematic errors in astrometry

Because of the way that Gaia surveys the sky, measurement errors appear to be spatially correlated on angular scales glyph[lessorsimilar] 0 5 . ◦ . In the first approximation, the two-point correlation function depends only on the angular distance between sources (Figure 15 in Lindegren et al. 2018, or left panel of Figure 1), and a suitable functional form is given by Equation 1 in Vasiliev (2019b). The latter paper introduces the mathematical formalism for taking into account these spatially correlated PM errors when performing model fits (in particular, measuring the mean PM of a cluster, or its internal kinematics - rotation and velocity dispersion). The uncertainties on the mean PM of almost all clusters are dominated by the systematic errors, and are at the level 0 05 . -0 1 mas yr . -1 .

## 4. Internal kinematics of globular clusters

We use this formalism to fit more sophisticated models, in which the PM of cluster stars has a radially-varying rotational component and a radially-varying [isotropic] dispersion. By calibrating the method on mock datasets with full account of spatially correlated systematic errors, we find that the rotation can be reliably measured when its peak amplitude exceeds ∼ 0 05 mas yr . -1 , and the PM dispersion when it exceeds ∼ 0 1 mas yr . -1 ; these conditions are fulfilled for ∼ 60 clusters mostly within 10 kpc and having at least 100 member stars passing all quality cuts.

We find high-significance rotation signatures in 8 clusters, and a weaker signatures in another 10 clusters; the results agree well with independent studies of Bianchini et al. (2018); Sollima et al. (2019); Jindal et al. (2019), even though none of these papers took into account systematic errors. The PM dispersion profiles for 60 clusters are also in a good agreement with independent PM measurements of Baumgardt et al. (2019) based on Gaia data, HST -based PM measurements (Watkins et al. 2015) in the central parts of 22 clusters, and various studies of line-of-sight velocity dispersion profiles, although for several most crowded clusters our PM dispersions are likely unreliable in the central parts. An example is shown in the right panel of Figure 1.

## 5. Dynamics of the entire Milky Way globular cluster system

With our measurements of mean PM of clusters, and the distances and line-of-sight velocities taken from the literature, we analyze the distribution of [almost] the entire population of Milky Way globular clusters in the 6d phase space. We confirm the previously known overall rotation of clusters (especially metal-rich ones) in the inner part ( r glyph[lessorsimilar] 10 kpc) of the Milky Way, with the mean v φ reaching ∼ 70 -80 kms -1 . The velocity dispersion is nearly isotropic ( ∼ 100 kms -1 ) in the inner part, and becomes more radially anisotropic further out. The distribution of clusters in the space of integrals of motion (actions) reveals several possible kinematic subgroups (Figure 2), which may be remnants of ancient accretion events.

We use clusters as dynamical probes of the Milky Way potential. Assuming a particular functional form for the DF of clusters in the action space (similar to Binney &amp; Wong 2017 and Posti &amp; Helmi 2019), we explore the range of parameters of the Milky Way potential and the DF consistent with the measured 6d phase-space coordinates of clusters (again convolved with observational uncertainties, which are mostly dominated by distance errors). We find an approximately flat rotation curve ( v circ glyph[similarequal] 200 -240 kms -1 between 10 and 100 kpc), consistent with the best-fit potential of McMillan (2017), but higher than the one from Bovy (2015). The enclosed mass within 50 kpc is inferred to be (0 45 . -0 65) . × 10 12 M glyph[circledot] , and within 100 kpc - (0 65 . -1 2) . × 10 12 M glyph[circledot] ; this is mostly consistent with other recent measurements (e.g., Watkins et al. 2019), but higher than the results of Eadie &amp; Juric (2019) based on the same kinematic sample.

This work uses the data from the European Space Agency mission Gaia ( https:// www.cosmos.esa.int/gaia ). This work was supported by the European Research council under the 7th Framework programme (grant No. 308024).

## References

Baumgardt H., Hilker M., Sollima A., Bellini A., 2019, MNRAS , 482, 5138 Bianchini P., van der Marel R., del Pino A., et al., 2018, MNRAS , 481, 2125 Binney J., Wong L.K., 2017, MNRAS , 467, 2446

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00009v1-figure-829eafff-ecb4-4c46-8156-3443b2e26bdc.jpg)

radial (

J φ

=0

)

Figure 2. Action-space map of the outer Milky Way globular clusters (a modified version of Figure 5 of Vasiliev 2019a). The horizontal and vertical coordinates show the actions J φ and J z -J r , normalized by the angular momentum L circ ( E ) of a circular orbit with the given energy; they roughly correspond to eccentricity and inclination. The color shows the energy, expressed in terms of the radius of a circular orbit. Individual clusters are shown by clouds of Monte Carlo samples representing the measurement uncertainties. Encircled are possible distinct kinematic subgroups identified by Law &amp; Majewski (2010); Myeong et al. (2018, 2019); Helmi et al. (2018).

Bovy J., 2015, ApJS , 216, 29

Eadie G., Juri´ M., 2019, c ApJ , 875, 159 Harris W., 1996, AJ , 112, 1487; arXiv:1012.3224 Gaia Collaboration (Brown et al.), 2018a, A&amp;A , 616, 1

Gaia Collaboration (Helmi et al.), 2018b, A&amp;A , 616, 12

Helmi A., Babusiaux C., Koppelman H., et al., 2018, Nature , 563, 85

Jindal A., Webb J., Bovy J., 2019, MNRAS , 487, 3693

Law D., Majewski S., 2010, ApJ , 718, 1128

Lindegren L., Hernandez J., Bombrun A., et al. ( Gaia Collaboration), 2018, A&amp;A , 616, 2

McMillan P., 2017, MNRAS , 465, 76

Myeong G.C., Evans N.W., Belokurov V., Sanders J., Koposov S., 2018, ApJL , 863, L28

Myeong G.C., Vasiliev E., Iorio G., Evans N.W., Belokurov V., 2019, MNRAS , 488, 1235

Posti L., Helmi A., 2019, A&amp;A , 621, 56

Sollima A., Baumgardt H., Hilker M., 2019, MNRAS , 485, 1460

Vasiliev E., 2019a, MNRAS , 484, 2832

Vasiliev E., 2019b, MNRAS , in press; arXiv:1811.05345

Watkins L., van der Marel R., Bellini A., Anderson J., 2015, ApJ , 803, 29

Watkins L., van der Marel R., Sohn S.T., Evans N.W., 2019, ApJ , 873, 118
Eugene Vasiliev
Using Gaia for studying Milky Way star clusters
We review the implications of the Gaia Data Release 2 catalogue for studying
the dynamics of Milky Way globular clusters, focusing on two separate topics.
The first one is the analysis of the full 6-dimensional phase-space
distribution of the entire population of Milky Way globular clusters: their
mean proper motions (PM) can be measured with an exquisite precision (down to
0.05 mas/yr, including systematic errors). Using these data, and a suitable
ansatz for the steady-state distribution function (DF) of the cluster
population, we then determine simultaneously the best-fit parameters of this DF
and the total Milky Way potential. We also discuss possible correlated
structures in the space of integrals of motion. The second topic addresses the
internal dynamics of a few dozen of the closest and richest globular clusters,
again using the Gaia PM to measure the velocity dispersion and internal
rotation, with a proper treatment of spatially correlated systematic errors.
Clear rotation signatures are detected in 10 clusters, and a few more show
weaker signatures at a level >=0.05$ mas/yr. PM dispersion profiles can be
reliably measured down to 0.1 mas/yr, and agree well with the line-of-sight
velocity dispersion profiles from the literature.