The Kinematics of 30 Milky Way Globular Clusters and the Multiple Stellar Populations within

The spectroscopic and photometric classification of multiple stellar populations (MPs) in Galactic globular clusters (GCs) has enabled comparisons between contemporary observations and formation theories regarding the initial spatial configurations of the MPs. However, the kinematics of these MPs is an aspect that requires more attention. We investigated the 3D kinematics of 30 Galactic GCs, extending to 3-5 half-light radii, as well as their MPs, in order to uncover clues of the initial conditions of GCs and the MPs within. We have combined Hubble Space Telescope and Gaia DR3 proper motions together with a comprehensive set of line-of-sight velocities to determine the 3D rotation amplitudes, rotation axes, and anisotropy profiles of the clusters. We include radial velocities from new IFU observations of NGC 5024 and an analysis of archival MUSE data of NGC 6101. We compare our kinematic results with structural and orbital parameters of each cluster, reporting the most significant correlations and common features. We find significant rotation in 21 GCs, with no significant differences between the total rotational amplitudes of the MPs, except for NGC 104. We find no significant differences in the position angles or inclination angles. We find that the 3D rotational amplitude is strongly correlated with mass, relaxation time, enriched star fraction and concentration. We determine the anisotropy profiles of each cluster and the MPs where possible. We investigate correlations with the structural parameters, orbital parameters and accretion history of the clusters, finding that the dynamically young clusters with the highest central concentrations of primordial stars show radial anisotropy in their outer regions ($>2$ half-light radii). The dynamically young clusters with a central concentration of enriched stars show significant tangential anisotropy or isotropy in their outer regions.