## Imprints of Evolution on the Internal Kinematics of Globular Clusters

1 2 ,

Laura L. Watkins with Roeland P. van der Marel , Andrea Bellini , 3 3 Mattia Libralato 3 and Jay Anderson 3

1 Department of Astrophysics, University of Vienna T¨ urkenschanzstraße 17, 1180 Vienna, Austria email: lwatkins@eso.org

2

European Southern Observatory Karl-Schwarzschild-Str. 2, 85748 Garching bei M¨nchen, Germany u 3 Space Telescope Science Institute 3700 San Martin Drive, Baltimore MD 21218, USA

Abstract. Globular clusters are collisional systems, meaning that the stars inside them interact on timescales much shorter than the age of the Universe. These frequent interactions transfer energy between stars and set up observable trends that tell the story of a cluster's evolution. This contribution focuses on what we can learn by studying velocity anisotropy and energy equipartition in globular clusters with Hubble Space Telescope proper motions.

Keywords. astrometry, globular clusters: general, stars: kinematics, stellar dynamics

## 1. Introduction

Globular clusters are old, dense systems. The high density of stars means that the timescale for interactions between stars is short enough to be meaningful. The old age of the clusters means that the stars inside have had a long time in which to interact. Taken together, the consequence is that stars inside clusters have experienced many interactions. While any changes that result from a single interaction are small, after many interactions the cumulative effects are significant and measurable.

The history of interactions between stars in a globular cluster is imprinted on its kinematics and manifests in a number of different ways. In this contribution, I will focus only on two of these: 1) velocity anisotropy; and 2) kinetic energy equipartition.

Studying velocity anisotropy requires at least two components of motion, which necessitates proper motion data. Studying energy equipartition typically requires deep astrometry, which necessitates data from the Hubble Space Telescope ( HST ). As such, the HST proper motion catalogues for 22 Milky Way globular clusters presented in Bellini et al. (2014) are ideally suited for this work.

## 2. Anisotropy

Velocity anisotropy is defined as the ratio between the velocity dispersions in two orthogonal directions. These can be any two orthogonal directions, and the choice of directions can depend on the structural properties of the system and the effects we wish to measure. In globular clusters, it is common to consider the projected radial and projected tangential velocity dispersions σ R and σ T .

In Watkins et al. (2015), we measured the projected tangential-radial velocity anisotropy

Figure 1. Left panel : Anisotropy profile for NGC 2808. The black points show the binned profile measured with HST proper motions. The blue lines are fits to the data; in this case we fit straight lines to the data in linear-radius, which appear curved in log-radius as shown. The red, green and orange lines show the half-light radius R half , the core radius R core , and 1 10 R half , which indicate the size of the cluster and the extent of the data. R half and R core were taken from Harris (1996, 2010 edition). Right panel : Anisotropy as a function of relaxation time for our sample of 22 clusters. Red points show the core values and the blue points show the half-light values. The grey lines show a fit that is flat for relaxation times shorter than some characteristic time and linearly decreasing at longer times. The dashed line marks the characteristic time.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00789v1-figure-7dd20ac8-62e0-4d20-9150-c70d27c1920e.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00789v1-figure-6799d6b1-0546-44df-ad4f-413ec8e6ca4a.jpg)

σ T /σ R in all 22 GCs in the Bellini et al. (2014) sample. The anisotropy profile for a single cluster, NGC 2808, is shown in the left panel of Figure 1. We see that the centre is isotropic ( σ T /σ R ∼ 1), and that the anisotropy becomes mildly radial at the outer limits of the area covered by the data ( σ T /σ R &lt; 1). In fact, we see this trend across all 22 clusters: all clusters are isotropic at their centres; some stay isotropic out to the limits of the data, and some become mildly radial in the outer regions . Figures 5-10 of Watkins † et al. (2015) show the corresponding profiles for the other 21 clusters in the sample.

We used the fits to estimate the anisotropy at the core and half-light radii, and compared these values to the relaxation times estimated at the core and half-light radii, from Harris (1996, 2010 edition). These estimates are shown in the right panel of Figure 1. We see indeed that the clusters are almost all isotropic out to their core radii; some are isotropic even out to their half-light radii, while some show mild radial anisotropy at the half-light radii. The degree of anisotropy increases with relaxation time.

To quantify this, to these points we fit a simple model that assumes a constant anisotropy for relaxation times shorter than some characteristic value and a linearly decreasing anisotropy for time larger than that value. The constant anisotropy, characteristic time and slope of the decrease were all left free. We see that the best fit shows a constant anisotropy for short relaxation times very close to 1 (isotropic) and an increase in the anisotropy in the radial direction for longer relaxation times.

What can we learn about the evolution of clusters from this? We suppose that GCs in the past had some radial anisotropy then as the stars inside interact, the cluster relaxes and signals of anisotropy are washed out. The densest regions (in the cores) have the shortest relaxation times and so reach isotropy first; the intermediate regions (near the half-light radii) have longer relaxation times so have not yet reached isotropy in some clusters and we still see mild anisotropy present. The cause of the radial anisotropy, whether it was primordial or manifested later in the lifetime of the cluster, is beyond the scope of this work to say.

† All of our fields are central fields, so the limits of the data are not the limits of the cluster.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00789v1-figure-fa2d12e7-2ead-4aa8-b870-6a941fd13cc3.jpg)

/circledot

Figure 2. Left panel : Velocity dispersion as a function of stellar mass in NGC 104 (47 Tuc). The black points were measured using HST proper motions. The yellow line shows a fit to the data with the resulting η estimate given in the top-right. Right panel : Distribution of η estimates for all 11 clusters for which we have performed this analysis.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00789v1-figure-8eeda0bf-f7ff-4e68-a08a-c1b0c566ba4e.jpg)

## 3. Energy Equipartition

When two stars pass close to one another, energy is transferred from the star with higher energy to the star with lower energy. Given enough interactions, it is expected that a stellar system will reach kinetic energy equipartition whereby all stars have the same kinetic energy. In complete equipartition we expect the velocity dispersion σ of the stars to depend on stellar mass m like σ ∝ m -1 2 . In practice, most clusters will be in a state of partial equipartition with dispersions that vary with stellar mass like σ ∝ m -η where η quantifies the degree of equipartition (0 for none and 0.5 for full).

Naively, we might expect that given long enough clusters will reach full equipartition. However Trenti &amp; van der Marel (2013) used a suite of N-body simulations to predict that clusters never reach full equipartition. The simulations started with no primordial equipartition and were left to evolve. The inner regions of the clusters reached a maximum η ∼ 0 2 and the degree of equipartition actually began to decrease again to level off at . η ∼ 0 1. The outer regions of the clusters increased steadily to . η ∼ 0 1 and leveled off . without first overshooting. (See their Figure 3.) We would like to test this predict with real data.

Whether partial or full, equipartition means that low mass stars tend to be moving faster than high mass stars, so to measure equipartition we need to measure velocities for stars over a wide range of stellar mass. The stars for which we can most easily measure kinematics (particularly those stars that have been most frequently targeted with spectrographs to obtain line-of-sight velocities) are the bright stars above the mainsequence turn-off. Unfortunately, the phases of stellar evolution above the turn-off are fast and so the range of stellar masses covered by these stars is small and not ideal for measuring equipartition (with the exception of blue stragglers stars, as discussed in Baldwin et al. 2016).

To get a decent range of stellar masses, we really need to measure velocities for stars along the main sequence. This is beyond the reach of existing spectrographs, except perhaps MUSE for the very closest clusters, and also beyond the reach of Gaia with sufficient accuracy, again except perhaps for one or two of the closest clusters. But this is achievable with HST , although this work is still challenging and pushes both the observatory and the data analysis techniques to their limits.

In Watkins et al. (in preparation) we have measured the velocity dispersion as a func-

00

125

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00789v1-figure-d4e01fe0-19f7-4849-b37c-d8322e94a07b.jpg)

150

Figure 3. Left panel : η profile for NGC 104 (47 Tuc). The green lines show the median η estimate in each radial bin; the dark shaded regions show the region spanning the 15.9 to 84.1 percentiles, and the light shaded regions the 2.5 to 97.5 percentiles. The green points mark the centre of each bin. The grey line with shaded region shows the average η estimated in the previous analysis. Right panel : The same for NGC6397. The dotted line marks η = 0, a theoretical minimum.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00789v1-figure-6a939b8c-56b1-487c-bc0c-1371fc29b6f4.jpg)

tion of stellar mass and estimated the degree of equipartition † for 11 of the 22 clusters in Bellini et al. (2014). The left panel of Figure 2 shows the results for NGC 104 (47 Tuc). We estimate η ∼ 0 2 in this cluster, in line with the theoretical expectations. The right . panel of Figure 2 shows a histogram of the η values obtain for all 11 clusters. We see that no cluster reaches full equipartition ( η ∼ 0 5), . consistent with Trenti &amp; van der Marel (2013), but we do see some values above η ∼ 0 2, inconsistent with their simulations. .

So far we have considered average η values, however the relaxation times in the dense cluster centres are shorter than the relaxation times in the less dense regions further out. We saw the consequence of this for anisotropy in Figure 1. So the next step is to measure how η changes as a function of radius. Figure 3 shows the η profiles for NGC 104 (left) and NGC6397 (right). NGC104 appears to have a flat η profile over the range of our dataset, while NGC 6397 shows a gradient in η that is highest at the centre and decreases with radius, consistent with expectations.

Recently Libralato et al. (2018) studied equipartition in NGC 362 using an improved proper motion catalogue that better controls systematics. They also saw a similar gradient in η across the range of the data. Their smaller uncertainties highlight the efficacy of the improved data reduction and proper-motion measurement techniques.

Acknowledgements : Support for this work was provided by grants for HST programs AR12845 and AR-15055, provided by the Space Telescope Science Institute, which is operated by AURA, Inc., under NASA contract NAS 5-26555. LLW acknowledges support from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme under grant agreement No 724857 (Consolidator Grant ArcheoDyn).

## References

Baldwin, A. T., Watkins, L. L., van der Marel, R. P., et al. 2016, ApJ , 827, 12 Bellini, A., Anderson, J., van der Marel, R. P., et al. 2014, ApJ , 797, 115 Harris, W. E. 1996, AJ , 112, 1487 Libralato, M., Bellini, A., van der Marel, R. P., et al. 2018, ApJ , 861, 99 Trenti, M., &amp; van der Marel, R. 2013, MNRAS , 435, 3272 Watkins, L. L., van der Marel, R. P., Bellini, A., &amp; Anderson, J. 2015, ApJ , 803, 29

† This was done using only stars within the 25th and 75th percentiles in radius to mitigate the effects of incompleteness in the dense centres and edge effects in the outer regions, and the fact that velocity dispersion also changes as a function of radius.
Laura L. Watkins Roeland P. van der Marel Andrea Bellini Mattia Libralato Jay Anderson
Imprints of Evolution on the Internal Kinematics of Globular Clusters
Globular clusters are collisional systems, meaning that the stars inside them
interact on timescales much shorter than the age of the Universe. These
frequent interactions transfer energy between stars and set up observable
trends that tell the story of a cluster's evolution. This contribution focuses
on what we can learn by studying velocity anisotropy and energy equipartition
in globular clusters with Hubble Space Telescope proper motions.