2024A&A...691A..94D

We present the first 3D kinematic analysis of multiple stellar populations (MPs) in a representative sample of 16 Galactic globular clusters (GCs). For each GC in the sample, we studied the MP line-of-sight, plane-of-the-sky and 3D rotation, and velocity distribution anisotropy. The differences between first-population (FP) and second-population (SP) kinematic patterns were constrained by means of parameters specifically defined to provide a global measure of the relevant physical quantities and to enable a meaningful comparison among different clusters. Our analysis provides the first observational description of the MP kinematic properties and of the path they follow during their long-term dynamical evolution. In particular, we find evidence of differences between the rotation of MPs along all velocity components with the SP preferentially rotating faster than the FP. The difference between the rotation strength of MPs is anticorrelated with the cluster dynamical age. We also observe that FPs are characterized by isotropic velocity distributions at any dynamical age probed by our sample. On the contrary, the velocity distribution of SP stars is found to be radially anisotropic in dynamically young clusters and isotropic at later evolutionary stages. The comparison with a set of numerical simulations shows that these observational results are consistent with the long-term evolution of clusters forming with an initially more centrally concentrated and more rapidly rotating SP subsystem. We discuss the possible implications these findings have on our understanding of MP formation and early evolution.