The Sagittarius dwarf spheroidal galaxy (Sgr dSph) is in an advanced stage of disruption but still hosts its nuclear star cluster (NSC), M54, at its center. In this paper, we present a detailed kinematic characterization of the three stellar populations present in M54: young metal-rich (YMR); intermediate-age metal-rich (IMR); and old metal-poor (OMP), based on the spectra of $\sim6500$ individual M54 member stars extracted from a large MUSE/VLT dataset. We find that the OMP population is slightly flattened with a low amount of rotation ($\sim0.8$ km s$^{-1}$) and with a velocity dispersion that follows a Plummer profile. The YMR population displays a high amount of rotation ($\sim5$ km s$^{-1}$) and a high degree of flattening, with a lower and flat velocity dispersion profile. The IMR population shows a high but flat velocity dispersion profile, with some degree of rotation ($\sim2$ km s$^{-1}$). We complement our MUSE data with information from \textit{Gaia DR2} and confirm that the stars from the OMP and YMR populations are comoving in 3D space, suggesting that they are dynamically bound. While dynamical evolutionary effects (e.g. energy equipartition) are able to explain the differences in velocity dispersion between the stellar populations, the strong differences in rotation indicate different formation paths for the populations, as supported by an $N$-body simulation tailored to emulate the YMR-OMP system. This study provides additional evidence for the M54 formation scenario proposed in our previous work, where this NSC formed via GC accretion (OMP) and in situ formation from gas accretion in a rotationally supported disc (YMR).

The outskirts of globular clusters (GCs) simultaneously retain crucial information about their formation mechanism and the properties of their host galaxy. Thanks to the advent of precision astrometry both their morphological and kinematic properties are now accessible. Here we present the first dynamical study of the outskirts of the retrograde GC NGC 3201 until twice its Jacobi radius (< 100 pc), using specifically-selected high-quality astrometric data from Gaia DR2. We report the discovery of a stellar overdensity along the South-East/North-West direction that we identify as tidal tails. The GC is characterized globally by radial anisotropy and a hint of isotropy in the outer parts, with an excess of tangential orbits around the lobes corresponding to the tidal tails, in qualitative agreement with an N-body simulation. Moreover, we measure flat velocity dispersion profiles, reaching values of $3.5\pm0.9$ km/s until beyond the Jacobi radius. While tidal tails could contribute to such a flattening, this high velocity dispersion value is in disagreement with the expectation from the sole presence of potential escapers. To explain this puzzling observation, we discuss the possibility of an accreted origin of the GC, the presence of a dark matter halo --leftover of its formation at high redshift -- and the possible effects of non-Newtonian dynamics. Our study uncovers a new path for the study of GC formation and of the properties of the Milky Way potential in the era of precision astrometry.

In weakly collisional stellar systems such as some globular clusters, partial energy equipartition and mass segregation are expected to develop as a result of the cumulative effect of stellar encounters even in systems initially characterized by star-mass independent density and energy distributions. In parallel, numerical simulations have demonstrated that radially-biased pressure anisotropy slowly builds up in realistic models of globular clusters from initial isotropic conditions, leading to anisotropy profiles that, to some extent, mimic those resulting from incomplete violent relaxation known to be relevant to elliptical galaxies. In this paper we consider a set of realistic simulations realized by means of Monte Carlo methods and analyze them by means of self-consistent two-component models. For the purpose, we refer to an underlying distribution function, originally conceived to describe elliptical galaxies, that has been recently truncated and adapted to the context of globular clusters. The two components are meant to represent light stars (combining all main sequence stars) and heavy stars (giants, dark remnants, and binaries). We show that this conceptually simple family of two-component truncated models provides a reasonable description of simulated density, velocity dispersion, and anisotropy profiles, especially for the most relaxed systems, with the capability to express quantitatively the attained levels of energy equipartition and mass segregation. In contrast, two-component isotropic models based on the King distribution function do not offer a comparably satisfactory representation of the simulated globular clusters. With this work we provide a new reliable diagnostic tool applicable to nonrotating globular clusters that are characterized by significant gradients in the local value of the mass-to-light ratio, beyond the commonly used one-component dynamical models.

The center of the Sagittarius dwarf spheroidal galaxy (Sgr dSph) hosts a nuclear star cluster, M54, which is the only galaxy nucleus that can be resolved into individual stars at optical wavelengths. It is thus a key target for understanding the formation of nuclear star clusters and their relation to globular clusters. We present a large Multi-Unit Spectroscopic Explorer (MUSE) data set that covers M54 out to $\sim$2.5 half-light radius, from which we extracted the spectra of $\sim$6600 cluster member stars. We use these data in combination with HST photometry to derive age and metallicity for each star. The stellar populations show a well defined age-metallicity relation, implying an extended formation history for the central region of Sgr dSph. We classify these populations into three groups, all with the same systemic velocity: young metal-rich (YMR; 2.2\,Gyr, \mbox{[Fe/H]$=-0.04$}); intermediate-age metal-rich (IMR; 4.3\,Gyr, \mbox{[Fe/H]$=-0.29$}); and old metal-poor (OMP; 12.2\,Gyr, \mbox{[Fe/H]$=-1.41$}). The YMR and OMP populations are more centrally concentrated than the IMR population, which are likely stars of the Sgr dSph. We suggest the OMP population is the result of accretion and merging of two or more old and metal-poor globular clusters dragged to the center by dynamical friction. The YMR is consistent with being formed by in situ star formation in the nucleus. The ages of the YMR population suggest that it may have been triggered into forming when the Sgr dSph began losing its gas during the most recent interaction with the Milky Way, $\sim$3\,Gyr ago.

The study of the kinematics of globular clusters (GCs) offers the possibility of unveiling their long term evolution and uncovering their yet unknown formation mechanism. Gaia DR2 has strongly revitalized this field and enabled the exploration of the 6D phase-space properties of Milky Way GCs, thanks to precision astrometry. However, to fully leverage on the power of precision astrometry, a thorough investigations of the data is required. In this contribution, we show that the study of the mean radial proper motion profiles of GCs offers an ideal benchmark to assess the presence of systematics in crowded fields. Our work demonstrates that systematics in Gaia DR2 for the closest 14 GCs are below the random measurement errors, reaching a precision of ~0.015 mas/yr for mean proper motion measurements. Finally, through the analysis of the tangential component of proper motions, we report the detection of internal rotation in a sample of ~50 GCs, and outline the implications of the presence of angular momentum for the formation mechanism of proto-GC. This result gives the first taste of the unparalleled power of Gaia DR2 for GCs science, in preparation for the subsequent data releases.

As a result of the slow action of two-body encounters, globular clusters develop mass segregation and attain a condition of only partial energy equipartition even in their central, most relaxed regions. Realistic numerical simulations show that, during the process, a radially-biased anisotropy profile slowly builds up, mimicking that resulting from incomplete violent relaxation. Commonly used dynamical models, such as the one-component King models, cannot describe these properties. Here we show that simple two-component models based on a distribution function originally conceived to describe elliptical galaxies, recently truncated and adapted to the context of globular clusters, can describe in detail what is observed in complex and realistic numerical simulations.

Omega Centauri is the most massive globular cluster of the Milky Way, and possesses many peculiar properties. In particular, the cluster contains distinct multiple stellar populations, with a large spread in metallicity and different kinematics as a function of light elements abundance, implying that it formed over an extended period of time. This has lead to the suggestion that $\omega$ Cen is the remnant core of an accreted dwarf galaxy. If this scenario is correct, $\omega$ Cen should be tidally limited, and one should expect to find tidal debris spread along its orbit. Here we use N-body simulations to show that the recently-discovered `Fimbulthul' structure, identified in the second data release (DR2) of the Gaia mission, is the long sought-for tidal stream of $\omega$ Cen, extending up to $28\deg$ from the cluster. Follow-up high-resolution spectroscopy of 5 stars in the stream show that they are closely-grouped in velocity, and have metallicities consistent with having originated in that cluster. Guided by our N-body simulations, we devise a selection filter that we apply to Gaia data to also uncover the portion of the stream in the highly-contaminated and crowded field within $10\deg$ of $\omega$ Cen. Further modelling of the stream may help to constrain the dynamical history of the dwarf galaxy progenitor of this disrupting system and guide future searches for its remnant stars in the Milky Way.

Line-of-sight kinematic studies indicate that many Galactic globular clusters have a significant degree of internal rotation. However, three-dimensional kinematics from a combination of proper motions and line-of-sight velocities are needed to unveil the role of angular momentum in the formation and evolution of these old stellar systems. Here we present the first quantitative study of internal rotation on the plane-of-the-sky for a large sample of globular clusters using proper motions from Gaia DR2. We detect signatures of rotation in the tangential component of proper motions for 11 out of 51 clusters at a $>$3-sigma confidence level, confirming the detection reported in Gaia collaboration et al. (2018) for 8 clusters, and additionally identify 11 GCs with a 2-sigma rotation detection. For the clusters with a detected global rotation, we construct the two-dimensional rotation maps and proper motion rotation curves, and we assess the relevance of rotation with respect to random motions ($V/\sigma\sim0.08-0.51$). We find evidence of a correlation between the degree of internal rotation and relaxation time, highlighting the importance of long-term dynamical evolution in shaping the clusters current properties. This is a strong indication that angular momentum must have played a fundamental role in the earliest phases of cluster formation. Finally, exploiting the spatial information of the rotation maps and a comparison with line-of-sight data, we provide an estimate of the inclination of the rotation axis for a subset of 8 clusters. Our work demonstrates the potential of Gaia data for internal kinematic studies of globular clusters and provides the first step to reconstruct their intrinsic three-dimensional structure.

Dynamical evolution drives globular clusters toward core collapse, which strongly shapes their internal properties. Diagnostics of core collapse have so far been based on photometry only, namely on the study of the concentration of the density profiles. Here we present a new method to robustly identify core-collapsed clusters based on the study of their stellar kinematics. We introduce the \textit{kinematic concentration} parameter, $c_k$, the ratio between the global and local degree of energy equipartition reached by a cluster, and show through extensive direct $N$-body simulations that clusters approaching core collapse and in the post-core collapse phase are strictly characterized by $c_k>1$. The kinematic concentration provides a suitable diagnostic to identify core-collapsed clusters, independent from any other previous methods based on photometry. We also explore the effects of incomplete radial and stellar mass coverage on the calculation of $c_k$ and find that our method can be applied to state-of-art kinematic datasets.

Galactic globular clusters (GCs) are believed to have formed in-situ in the Galaxy as well as in dwarf galaxies later accreted onto the Milky Way. However, to date, there is no unambiguous signature to distinguish accreted GCs. Using specifically designed $N$-body simulations of GCs evolving in a variety of time-dependent tidal fields (describing the potential of a dwarf galaxy-Milky Way merger), we analyze the effects imprinted to the internal kinematics of an accreted GC. In particular, we look at the evolution of the velocity anisotropy. Our simulations show that at early phases, the velocity anisotropy is determined by the tidal field of the dwarf galaxy and subsequently the clusters will adapt to the new tidal environment, losing any signature of their original environment in a few relaxation times. At 10 Gyr, GCs exhibit a variety of velocity anisotropy profiles, namely, isotropic velocity distribution in the inner regions and either isotropy or radial/tangential anisotropy in the intermediate and outer regions. Independently of an accreted origin, the velocity anisotropy primarily depends on the strength of the tidal field cumulatively experienced by a cluster. Tangentially anisotropic clusters correspond to systems that have experienced stronger tidal fields and are characterized by higher tidal filling factor, $r_{50}/r_j\gtrsim0.17$, higher mass loss $\gtrsim60\%$ and relaxation times $t_{rel}\lesssim10^9$ Gyr. Interestingly, we demonstrate that the presence of tidal tails can significantly contaminate the measurements of velocity anisotropy when a cluster is observed in projection. Our characterization of the velocity anisotropy profiles in different tidal environments provides a theoretical benchmark for the interpretation of the unprecedented amount of three-dimensional kinematic data progressively available for Galactic GCs.

High-precision proper motions of the globular cluster 47 Tuc have allowed us to measure for the first time the cluster rotation in the plane of the sky and the velocity anisotropy profile from the cluster core out to about 13'. These profiles are coupled with prior measurements along the line of sight and the surface-brightness profile, and fit all together with self-consistent models specifically constructed to describe quasi-relaxed stellar systems with realistic differential rotation, axisymmetry and pressure anisotropy. The best-fit model provides an inclination angle i between the rotation axis and the line-of-sight direction of 30 deg, and is able to simultaneously reproduce the full three-dimensional kinematics and structure of the cluster, while preserving a good agreement with the projected morphology. Literature models based solely on line-of-sight measurements imply a significantly different inclination angle (i=45 deg), demonstrating that proper motions play a key role in constraining the intrinsic structure of 47 Tuc. Our best-fit global dynamical model implies an internal rotation higher than previous studies have shown, and suggests a peak of the intrinsic V/sigma ratio of ~0.9 at around two half-light radii, with a non-monotonic intrinsic ellipticity profile reaching values up to 0.45. Our study unveils a new degree of dynamical complexity in 47 Tuc, which may be leveraged to provide new insights into the formation and evolution of globular clusters.

The internal dynamics of globular clusters (GCs) is strongly affected by two-body interactions that bring the systems to a state of partial energy equipartition. Using a set of Monte Carlo clusters simulations, we investigate the role of the onset of energy equipartition in shaping the mass-to-light ratio (M/L) in GCs. Our simulations show that the M/L profiles cannot be considered constant and their specific shape strongly depends on the dynamical age of the clusters. Dynamically younger clusters display a central peak up to M/L $\simeq25$ $M_\odot/L_\odot$ caused by the retention of dark remnants; this peak flattens out for dynamically older clusters. Moreover, we find that also the global values of M/L correlate with the dynamical state of a cluster quantified as either the number of relaxation times a system has experienced $n_{rel}$ or the equipartition parameter $m_{eq}$: clusters closer to full equipartition (higher $n_{rel}$ or lower $m_{eq}$) display a lower M/L. We show that the decrease of M/L is primarily driven by the dynamical ejection of dark remnants, rather than by the escape of low-mass stars. The predictions of our models are in good agreement with observations of GCs in the Milky Way and M31, indicating that differences in relaxation state alone can explain variations of M/L up to a factor of $\simeq3$. Our characterization of the M/L as a function of relaxation state is of primary relevance for the application and interpretation of dynamical models.

The detection of intermediate mass black holes (IMBHs) in Galactic globular clusters (GCs) has so far been controversial. In order to characterize the effectiveness of integrated-light spectroscopy through integral field units, we analyze realistic mock data generated from state-of-the-art Monte Carlo simulations of GCs with a central IMBH, considering different setups and conditions varying IMBH mass, cluster distance, and accuracy in determination of the center. The mock observations are modeled with isotropic Jeans models to assess the success rate in identifying the IMBH presence, which we find to be primarily dependent on IMBH mass. However, even for a IMBH of considerable mass (3% of the total GC mass), the analysis does not yield conclusive results in 1 out of 5 cases, because of shot noise due to bright stars close to the IMBH line-of-sight. This stochastic variability in the modeling outcome grows with decreasing BH mass, with approximately 3 failures out of 4 for IMBHs with 0.1% of total GC mass. Finally, we find that our analysis is generally unable to exclude at 68% confidence an IMBH with mass of $10^3~M_\odot$ in snapshots without a central BH. Interestingly, our results are not sensitive to GC distance within 5-20 kpc, nor to mis-identification of the GC center by less than 2'' (<20% of the core radius). These findings highlight the value of ground-based integral field spectroscopy for large GC surveys, where systematic failures can be accounted for, but stress the importance of discrete kinematic measurements that are less affected by stochasticity induced by bright stars.

Proper motions (PMs) are crucial to fully understand the internal dynamics of globular clusters (GCs). To that end, the Hubble Space Telescope (HST) Proper Motion (HSTPROMO) collaboration has constructed large, high-quality PM catalogues for 22 Galactic GCs. We highlight some of our exciting recent results: the first directly-measured radial anisotropy profiles for a large sample of GCs; the first dynamical distance and mass-to-light (M/L) ratio estimates for a large sample of GCs; and the first dynamically-determined masses for hundreds of blue-straggler stars (BSSs) across a large GC sample.

We describe the dynamical evolution of a unique type of dark star cluster model in which the majority of the cluster mass at Hubble time is dominated by an intermediate-mass black hole (IMBH). We analyzed results from about 2000 star cluster models (Survey Database I) simulated using the Monte-Carlo code MOCCA and identified these dark star cluster models. Taking one of these models, we apply the method of simulating realistic "mock observations" by utilizing the COCOA and SISCO codes to obtain the photometric and kinematic observational properties of the dark star cluster model at 12 Gyr. We find that the perplexing Galactic globular cluster NGC 6535 closely matches the observational photometric and kinematic properties of the dark star cluster model presented in this paper. Based on our analysis and currently observed properties of NGC 6535, we suggest that this globular cluster could potentially harbour an IMBH. If it exists, the presence of this IMBH can be detected robustly with proposed kinematic observations of NGC 6535.

Nuclear Star Clusters (NSCs) are commonly observed in the centres of most galactic nuclei, including our own Milky Way. While their study can reveal important information about the build-up of the innermost regions of galaxies, the physical processes that regulate their formation are still poorly understood. NSCs might have been formed through gas infall and subsequent in situ star formation, and/or through the infall and merging of multiple star clusters into the centre of the galaxy. Here, we investigate the viability of the latter, by studying direct N-body simulations of inspiralling clusters to the centre of a Milky-Way-like nuclear bulge that hosts a massive black hole. We find that the NSC that forms through this process can show both morphological and kinematical properties that make it comparable with observations of the Milky Way NSC, including significant rotation- a fact that has been so far attributed mainly to gas infall. We explore its kinematic evolution, to see if and how the merger history can imprint fossil records on its dynamical structure. Moreover, we study the effect of stellar foreground contamination in the line-of-sight kinematics of the NSC. Our study shows that no fine tuning of the orientation of the infalling globular clusters is necessary to result in a rotating NSC. We suggest that cluster-inspiral is a viable mechanism for the formation of rotating NSCs.

We make use of the Hubble Space Telescope proper-motion catalogs derived by Bellini et al. (2014) to produce the first radial velocity-dispersion profiles sigma(R) for blue straggler stars (BSSs) in Galactic globular clusters (GCs), as well as the first dynamical estimates for the average mass of the entire BSS population. We show that BSSs typically have lower velocity dispersions than stars with mass equal to the main-sequence turnoff mass, as one would expect for a more massive population of stars. Since GCs are expected to experience some degree of energy equipartition, we use the relation sigma~M^-eta, where eta is related to the degree of energy equipartition, along with our velocity-dispersion profiles to estimate BSS masses. We estimate eta as a function of cluster relaxation from recent Monte Carlo cluster simulations by Bianchini et al. (2016b) and then derive an average mass ratio M_BSS /M_MSTO=1.50+/-0.14 and an average mass M_BSS=1.22+/-0.12 M_Sun from 598 BSSs across 19 GCs. The final error bars include any systematic errors that are random between different clusters, but not any potential biases inherent to our methodology. Our results are in good agreement with the average mass of M_BSS=1.22+/-0.06 M_Sun for the 35 BSSs in Galactic GCs in the literature with properties that have allowed individual mass determination.

High-precision kinematic studies of globular clusters require an accurate knowledge of all possible sources of contamination. Amongst other sources, binary stars can introduce systematic biases in the kinematics. Using a set of Monte Carlo cluster simulations with different concentrations and binary fractions, we investigate the effect of unresolved binaries on proper-motion dispersion profiles, treating the simulations like HST proper-motion samples. Since globular clusters evolve towards a state of partial energy equipartition, more massive stars lose energy and decrease their velocity dispersion. As a consequence, on average, binaries have a lower velocity dispersion, since they are more massive kinematic tracers. We show that, in the case of clusters with high binary fraction (initial binary fraction of 50%) and high concentration (i.e., closer to energy equipartition), unresolved binaries introduce a color-dependent bias in the velocity dispersion of main-sequence stars of the order of 0.1-0.3 km s$^{-1}$ (corresponding to $1-6$% of the velocity dispersion), with the reddest stars having a lower velocity dispersion, due to the higher fraction of contaminating binaries. This bias depends on the ability to distinguish binaries from single stars, on the details of the color-magnitude diagram and the photometric errors. We apply our analysis to the HSTPROMO data set of NGC 7078 (M15) and show that no effect ascribable to binaries is observed, consistent with the low binary fraction of the cluster. Our work indicates that binaries do not significantly bias proper-motion velocity-dispersion profiles, but should be taken into account in the error budget of kinematic analyses.

Two-body interactions play a major role in shaping the structural and dynamical properties of globular clusters (GCs) over their long-term evolution. In particular, GCs evolve toward a state of partial energy equipartition that induces a mass-dependence in their kinematics. By using a set of Monte Carlo cluster simulations evolved in quasi-isolation, we show that the stellar mass dependence of the velocity dispersion $\sigma(m)$ can be described by an exponential function $\sigma^2\propto \exp(-m/m_\mathrm{eq})$, with the parameter $m_\mathrm{eq}$ quantifying the degree of partial energy equipartition of the systems. This simple parametrization successfully captures the behaviour of the velocity dispersion at lower as well as higher stellar masses, that is, the regime where the system is expected to approach full equipartition. We find a tight correlation between the degree of equipartition reached by a GC and its dynamical state, indicating that clusters that are more than about 20 core relaxation times old, have reached a maximum degree of equipartition. This equipartition$-$dynamical state relation can be used as a tool to characterize the relaxation condition of a cluster with a kinematic measure of the $m_\mathrm{eq}$ parameter. Vice versa, the mass-dependence of the kinematics can be predicted knowing the relaxation time solely on the basis of photometric measurements. Moreover, any deviations from this tight relation could be used as a probe of a peculiar dynamical history of a cluster. Finally, our novel approach is important for the interpretation of state-of-the-art Hubble Space Telescope proper motion data, for which the mass dependence of kinematics can now be measured, and for the application of modeling techniques which take into consideration multi-mass components and mass segregation.

The detection of intermediate mass black holes in the centres of globular clusters is highly controversial, as complementary observational methods often deliver significantly different results. In order to understand these discrepancies, we develop a procedure to simulate integral field unit (IFU) observations of globular clusters: Simulating IFU Star Cluster Observations (SISCO). The input of our software are realistic dynamical models of globular clusters that are then converted in a spectral data cube. We apply SISCO to Monte Carlo cluster simulations from Downing et al. (2010), with a realistic number of stars and concentrations. Using independent realisations of a given simulation we are able to quantify the stochasticity intrinsic to the problem of observing a partially resolved stellar population with integrated-light spectroscopy. We show that the luminosity-weighted IFU observations can be strongly biased by the presence of a few bright stars that introduce a scatter in the velocity dispersion measurements up to $\simeq$40% around the expected value, preventing any sound assessment of the central kinematic and a sensible interpretation of the presence/absence of an intermediate mass black hole. Moreover, we illustrate that, in our mock IFU observations, the average kinematic tracer has a mass of $\simeq$0.75 solar masses, only slightly lower than the mass of the typical stars examined in studies of resolved line-of-sight velocities of giant stars. Finally, in order to recover unbiased kinematic measurements we test different masking techniques that allow us to remove the spaxels dominated by bright stars, bringing the scatter down to a level of only a few percent. The application of SISCO will allow to investigate state-of-the-art simulations as realistic observations.

The detection of intermediate mass black holes (IMBHs) in globular clusters has been hotly debated, with different observational methods delivering different outcomes for the same object. In order to understand these discrepancies, we construct detailed mock integral field spectroscopy (IFU) observations of globular clusters, starting from realistic Monte Carlo cluster simulations. The output is a data cube of spectra in a given field-of-view that can be analyzed in the same manner as real observations and compared to other (resolved) kinematic measurement methods. We show that the main discrepancies arise because the luminosity-weighted IFU observations can be strongly biased by the presence of a few bright stars that introduce a scatter in velocity dispersion measurements of several km/s. We show that this intrinsic scatter can prevent a sound assessment of the central kinematics, and therefore should be fully taken into account to correctly interpret the signature of an IMBH.

The distinction between globular clusters and dwarf galaxies has been progressively blurred by the recent discoveries of several extended star clusters, with size (20-30 pc) and luminosity (-6 < Mv < -2) comparable to the one of faint dwarf spheroidals. In order to explain their sparse structure, it has been suggested that they formed as star clusters in dwarf galaxy satellites that later accreted onto the Milky Way. If these clusters form in the centre of dwarf galaxies, they evolve in a tidally-compressive environment where the contribution of the tides to the virial balance can become significant, and lead to a super-virial state and subsequent expansion of the cluster, once removed. Using N-body simulations, we show that a cluster formed in such an extreme environment undergoes a sizable expansion, during the drastic variation of the external tidal field due to the accretion process. However, we show that the expansion due to the removal of the compressive tides is not enough to explain the observed extended structure, since the stellar systems resulting from this process are always more compact than the corresponding clusters that expand in isolation due to two-body relaxation. We conclude that an accreted origin of extended globular clusters is unlikely to explain their large spatial extent, and rather favor the hypothesis that such clusters are already extended at the stage of their formation.

We present the first study of high-precision internal proper motions (PMs) in a large sample of globular clusters, based on Hubble Space Telescope (HST) data obtained over the past decade with the ACS/WFC, ACS/HRC, and WFC3/UVIS instruments. We determine PMs for over 1.3 million stars in the central regions of 22 clusters, with a median number of ~60,000 stars per cluster. These PMs have the potential to significantly advance our understanding of the internal kinematics of globular clusters by extending past line-of-sight (LOS) velocity measurements to two- or three-dimensional velocities, lower stellar masses, and larger sample sizes. We describe the reduction pipeline that we developed to derive homogeneous PMs from the very heterogeneous archival data. We demonstrate the quality of the measurements through extensive Monte-Carlo simulations. We also discuss the PM errors introduced by various systematic effects, and the techniques that we have developed to correct or remove them to the extent possible. We provide in electronic form the catalog for NGC 7078 (M 15), which consists of 77,837 stars in the central 2.4 arcmin. We validate the catalog by comparison with existing PM measurements and LOS velocities, and use it to study the dependence of the velocity dispersion on radius, stellar magnitude (or mass) along the main sequence, and direction in the plane of the sky (radial/tangential). Subsequent papers in this series will explore a range of applications in globular-cluster science, and will also present the PM catalogs for the other sample clusters.

Aims: We present the first in-depth study of the kinematic properties and derive the structural parameters of NGC 4372 based on the fit of a Plummer profile and a rotating, physical model. We explore the link between internal rotation to different cluster properties and together with similar studies of more GCs, we put these in the context of globular cluster formation and evolution. Methods: We present radial velocities for 131 cluster member stars measured from high-resolution FLAMES/GIRAFFE observations. Their membership to the GC is additionally confirmed from precise metallicity estimates. Using this kinematic data set we build a velocity dispersion profile and a systemic rotation curve. Additionally, we obtain an elliptical number density profile of NGC 4372 based on optical images using a MCMC fitting algorithm. From this we derive the cluster's half-light radius and ellipticity as r_h=3.4'+/-0.04' and e=0.08+/-0.01. Finally, we give a physical interpretation of the observed morphological and kinematic properties of this GC by fitting an axisymmetric, differentially rotating, dynamical model. Results: Our results show that NGC 4372 has an unusually high ratio of rotation amplitude to velocity dispersion (1.2 vs. 4.5 km/s) for its metallicity. This, however, puts it in line with two other exceptional, very metal-poor GCs - M 15 and NGC 4590. We also find a mild flattening of NGC 4372 in the direction of its rotation. Given its old age, this suggests that the flattening is indeed caused by the systemic rotation rather than tidal interactions with the Galaxy. Additionally, we estimate the dynamical mass of the GC M_dyn=2.0+/-0.5 x 10^5 M_Sun based on the dynamical model, which constrains the mass-to-light ratio of NGC 4372 between 1.4 and 2.3 M_Sun/L_Sun, representative of an old, purely stellar population.

The Hubble Space Telescope (HST) has proven to be uniquely suited for the measurement of proper motions (PMs) of stars and galaxies in the nearby Universe. Here we summarize the main results and ongoing studies of the HSTPROMO collaboration, which over the past decade has executed some two dozen observational and theoretical HST projects on this topic. This is continuing to revolutionize our dynamical understanding of many objects, including: globular clusters; young star clusters; stars and stellar streams in the Milky Way halo; Local Group galaxies, including dwarf satellite galaxies, the Magellanic Clouds, and the Andromeda galaxy; and AGN Black Hole Jets.

Internal rotation is considered to play a major role in the dynamics of some globular clusters. However, in only few cases it has been studied by quantitative application of realistic and physically justified global models. Here we present a dynamical analysis of the photometry and three-dimensional kinematics of omega Cen, 47 Tuc, and M15, by means of a recently introduced family of self-consistent axisymmetric rotating models. The three clusters, characterized by different relaxation conditions, show evidence of differential rotation and deviations from sphericity. The combination of line-of-sight velocities and proper motions allows us to determine their internal dynamics, predict their morphology, and estimate their dynamical distance. The well-relaxed cluster 47 Tuc is very well interpreted by our model; internal rotation is found to explain the observed morphology. For M15, we provide a global model in good agreement with the data, including the central behavior of the rotation profile and the shape of the ellipticity profile. For the partially relaxed cluster omega Cen, the selected model reproduces the complex three-dimensional kinematics; in particular the observed anisotropy profile, characterized by a transition from isotropy, to weakly-radial anisotropy, and then to tangential anisotropy in the outer parts. The discrepancy found for the steep central gradient in the observed line-of-sight velocity dispersion profile and for the ellipticity profile is ascribed to the condition of only partial relaxation of this cluster and the interplay between rotation and radial anisotropy.

Internal rotation is considered to play a major role in determining the structure and dynamics of some globular clusters. We present a dynamical analysis of the photometry and three-dimensional kinematics of 47 Tuc and omega Cen, by means of a new family of self-consistent axisymmetric rotating models. The combined use of line-of-sight velocities and proper motions allows us to obtain a global description of the internal dynamical structure of the objects together with an estimate of their dynamical distances. The well-relaxed cluster 47 Tuc is very well interpreted by our dynamical models; in particular, internal rotation is found to explain the observed morphology. For the partially relaxed cluster omega Cen, the selected model provides a good representation of its complex three-dimensional kinematics, in general qualitative agreement with the observed anisotropy profile, which is characterized by tangential anisotropy in the outer parts; discrepancies are found between the observed and the expected ellipticity profile and are ascribed to the presence of a high degree of radial anisotropy in the intermediate region and to its interplay with rotation.