Kinematical evolution of tidally limited star clusters: the role of retrograde stellar orbits
The presence of an external tidal field often induces significant dynamical evolutionary effects on the internal kinematics of star clusters. Previous studies investigating the restricted three-body problem with applications to star cluster dynamics have shown that unbound stars on retrograde orbits (with respect to the direction of the cluster’s orbit) are more stable against escape than prograde orbits, and predicted that a star cluster might acquire retrograde rotation through preferential escape of stars on prograde orbits. In this study we present evidence of this prediction, but we also illustrate that there are additional effects that cannot be accounted for by the preferential escape of prograde orbits alone. Specifically, in the early evolution, initially underfilling models increase their fraction of retrograde stars without losing significant mass, and acquire a retrograde angular velocity. We attribute this effect to the development of preferentially eccentric/radial orbits in the outer regions of star clusters as they are expanding into their tidal limitation.
We explore the implications of the evolution of the fraction of prograde and retrograde stars for the evolution of the cluster internal rotation, and its dependence on the initial structural properties. Although all the systems studied here evolve towards an approximately solid-body internal rotation with angular velocity equal to about half of the angular velocity of the cluster orbital motion around the host galaxy, the evolutionary history of the radial profile of the cluster internal angular velocity depends on the cluster initial structure.
keywords:galaxies: star clusters: general; Galaxy: globular clusters: general; methods: numerical;
The combined effects of internal two-body relaxation and the external potential of the host galaxy play a major role in determining the structural and kinematical evolution of globular clusters, the rate at which clusters lose mass and their lifetimes. Early works in this regard include, for example, Vesperini & Heggie (1997), Giersz & Heggie (1997), and Baumgardt & Makino (2003), and a number of more recent investigations including Madrid et al. (2012), Webb et al. (2014), Haghi et al. (2014), and Zonoozi et al. (2016) have further explored these issues. The complex interplay between internal and external evolutionary processes alters in a non-trivial way the initial cluster properties and complicates any effort aimed at understanding which of the current cluster properties were imprinted at the epoch of its formation and which are instead the result of its evolutionary history (see e.g. Zonoozi et al., 2011, 2014; Webb & Leigh, 2015, for some studies attempting to reconstruct the cluster initial properties).
The observational characterization of the internal kinematics of several Galactic globular clusters will soon reach an unprecedented level of richness, thanks the synergy between the forthcoming astrometric data provided by Gaia and HST (see e.g. Richer et al., 2013; Bellini et al., 2014, 2015; Watkins et al., 2015; Sollima et al., 2015), and a number of ESO/VLT spectroscopic programs (see e.g. Lanzoni et al., 2013; Lapenna et al., 2015; Lardo et al., 2015). The knowledge of cluster internal kinematical properties opens the way to major advances in our understanding of cluster formation and dynamical evolution and of the role played by different internal and external evolutionary processes in shaping the current cluster properties.
In Tiongco et al. (2016), we have presented the results of a survey of simulations aimed at exploring the link between anisotropy in the velocity distribution and cluster initial structural properties, dynamical phase and mass loss history.
In this paper, we will focus on the characterization of the stellar orbital properties of a star cluster evolving in the external potential of the host galaxy, and we will discuss their implications on the global kinematics of the cluster.
A number of numerical studies (Keenan & Innanen, 1975; Fukushige & Heggie, 2000; Ernst et al., 2007; Zotos, 2015) have shown that unbound stars still residing within the cluster are more stable against escape if they are on retrograde orbits (with respect to the direction of the cluster’s orbit), with non-trivial implications on the definition of boundary of a tidally perturbed stellar system in coordinate and velocity space (see Read et al., 2006; Gajda & Łokas, 2016). As pointed out by Keenan & Innanen (1975), the preferential loss of prograde stars and the ensuing excess of retrograde stars implies that the process of mass loss may be accompanied by the development of an internal retrograde rotation. A study of the evolution of the velocity moments of -body models performed by Baumgardt & Makino (2003) has provided some evidence of the presence of such a behaviour, especially in the outer parts of tidally perturbed star clusters.
At the same time, the interaction with the external tidal field of the host galaxy might in part offset this effect and drive the system towards an internal rotation synchronous with that associated with the cluster orbital motion. Such a condition of “synchronization” of the internal and orbital motion (“tidal locking”) is indeed a standard assumption underlying many investigations of the dynamics of star clusters, which often rely on the classical formulation of the Hill’s problem (see e.g. Chandrasekhar, 1942; Heggie & Hut, 2003).
In this paper we present a detailed analysis of the long-term evolution of the fractions of the prograde and retrograde orbiting populations in -body models of star clusters evolving in an external tidal field. We study the evolution of clusters models for a range of different initial structural properties, and we explore how these properties affect the evolution of the fraction of prograde and retrograde stars and the implications for the cluster internal kinematics, especially its angular momentum. We explore in detail the dynamical processes behind the evolution of the fraction of prograde and retrograde stars, both for bound stars and for the population of unbound stars still residing within the cluster (“potential escapers”).
The outline of the paper is the following. In Section 2, we describe our method and the initial conditions adopted. We present our results in Section 3, and in Section 4 we summarize our conclusions.
2 Method and Initial Conditions
We have carried out a survey of -body simulations using the NBODY6 code (Aarseth, 2003) accelerated by a GPU (Nitadori & Aarseth, 2012); the simulations were run on the BIG RED II cluster at Indiana University.
For all the simulations, we have assumed clusters on circular orbits in the host galaxy external field modelled as that of a point-mass. The coordinate system in which equations of motion are solved is corotating with the cluster around the host galaxy with angular velocity (see e.g. Heggie & Hut, 2003). For the primary set of simulations, we assume that the system has no internal rotation in this reference frame; in a non-rotating reference frame, this corresponds to a rotation synchronous with that of the cluster around the host galaxy. Hereafter we shall refer to these set of models as ‘locked’ models. To explore the potential dependence of the evolution of the kinematics on whether or not the system is initially locked, we have also considered some models that are initially non-rotating in a non-rotating frame reference. In the rotating frame of reference described above these models have a retrograde solid-body rotation with angular velocity equal to . We will refer to these models as ‘unlocked’ models, and identify them by adding U to the model ID (see Table 1).
We consider systems with 16 384 equal-mass particles, and particles moving beyond a distance from the cluster centre equal to two times the Jacobi radius, , are removed from the simulation. All the simulations are run until 75% of the initial cluster mass is lost.
The initial conditions of our systems are -body realizations of isotropic, King (1966) models, with dimensionless central potential , and different values of the initial filling factor, defined as the ratio of the three-dimensional half-mass radius to the Jacobi radius or equivalently, the ratio of the model truncation radius to the Jacobi radius (see Table 1). These models were originally presented in the context of a study of the long-term evolution of the velocity anisotropy of tidally limited star clusters and we refer the reader to Tiongco et al. (2016) for a complete description (notice that there was a typo in reporting the of the model KF0125 in Table 1 of Tiongco et al. 2016. The correct value is shown in this paper).
By setting the value of according to the desired filling factor, is determined as
As increases for more underfilling models, decreases and this implies that the initial conditions and evolution of the unlocked versions of our underfilling models do not differ significantly from those of the corresponding locked versions. In consideration of this, we focus our attention only on unlocked models with relatively large filling factors, KF1U, KF075U and KF05U.
We emphasize that in the rest of the paper, whenever we discuss cluster rotation or refer to the cluster internal angular velocity, , we will refer to the angular velocity as measured in the co-rotating reference frame described above. In this reference frame, prograde and retrograde orbits are, respectively, those with larger or smaller than zero and an orbit with zero angular velocity is that of a star in corotation with the cluster around the host galaxy.
3.1 Evolution of the angular velocity profile
We approach the characterization of the kinematical properties of the models by studying the evolution of the radial profile of the internal angular velocity, .
The angular velocity is measured in the direction parallel to the - axis, binning the particles in equally thick cylindrical shells with heights encompassing the entire cluster at different projected radii, . We have combined 5 snapshots around the desired times. In the construction of the radial profiles, all particles within the Jacobi radius have been used.
The radial profiles of , normalized to the angular velocity of the cluster’s orbit around the host galaxy, , are illustrated in Fig. 1. Each panel offers a comparison between the profiles of all models evaluated at four different stages of their evolution, identified by the amount of mass that is retained by the system (), where denotes the total mass of all particles enclosed in the Jacobi radius.
All the models that are initially locked start with , while the unlocked models have, in the rotating frame, a solid body retrograde rotation (). In all cases, as the locked systems start to lose mass they all acquire a retrograde rotation in their outer regions (beyond approximately 0.1 ; see the last panel of Fig. 1).
Besides providing a quantitative measurement of the effect of the preferential loss of stars in prograde orbits on the cluster kinematics, this figure reveals several interesting and non-trivial results.
First of all, it appears that the early evolution of the radial profile of depends on the cluster initial filling factor for the locked models. Strongly underfilling models develop a strong retrograde angular velocity already at the beginning of their evolution (see the behaviour of KF0125 and KF025 in the first panel), while more filling models (i.e., KF075, KF1) are initially characterized by a more moderate retrograde rotation. As the evolution proceeds and the systems lose more mass, the behaviour of the angular velocity of the filling models converge to an approximately retrograde solid-body rotation such that , starting with the outermost regions and later extending to the inner regions. The model KF05 is an intermediate case, having a stronger retrograde angular velocity than the more filling models at = 0.95, but follows a similar behaviour by the time it has lost 50% of its mass. The most underfilling models, KF025 and KF0125, and the unlocked models take longer to converge but eventually at = 0.25 all models share the same approximately solid-body rotation with -0.5 at .
Both initially locked and unlocked systems evolve towards a common radial profile of but the relative importance of different dynamical processes in driving evolution towards the final solid-body rotation profile depends on the cluster initial properties. The more filling and initially locked systems lose their initial synchronous rotation as a result of the preferential loss of prograde stars (in this respect, the role of the so-called “potential escapers” will be discussed in detail in Section 3.4). Initially underfilling models, on the other hand, experience an early evolutionary phase during which their outer regions develop a strong retrograde rotation which may be interpreted in light of an additional mechanism not relying on mass loss (we will further discuss in Sections 3.2 and 3.3 the dynamics underpinning this mechanism).
In the subsequent evolution, while the interaction with the tidal field drives the systems towards a synchronization of the internal angular velocity with that of the cluster orbital motion around the host galaxy ( in the co-rotating reference frame), the preferential loss of particles on prograde orbits counteracts in part this effect. The net result is the retrograde approximate solid-body rotation shown in Fig. 1.
We wish to emphasize that for the purpose of the current investigation, we will focus exclusively on understanding the physical mechanisms that are responsible for the presence of an approximate solid-body rotation in the intermediate to outer parts of the models under consideration. The role played by the presence of non-vanishing total angular momentum in the dynamical evolution of collisional stellar systems is indeed a much broader topic (see e.g. Goodman, 1983; Einsel & Spurzem, 1999), and in particular, the long-term kinematical evolution of models with differential rotation in a tidal field has an added layer of complexity that has been only partially explored (see e.g. Ernst et al., 2007; Hong et al., 2013). We plan to fully address the interplay between internal rotation and external tidal field in a future study (Tiongco et al., in prep.).
3.2 Early evolution of the global fraction of retrograde orbits
The possible development of retrograde rotation has been often associated exclusively with the preferential loss of prograde orbiting stars. In this section, we delve further into this issue and present the results of our analysis aimed at quantifying the evolution of the populations of prograde and retrograde orbiting stars. In Fig. 2, we show the time evolution of , the global fraction of the total number of stars within the Jacobi radius on retrograde orbits.
We focus our attention first on the initially locked models; the time evolution of the global fraction of stars on retrograde orbits is characterized by an initial progressive growth, at a rate approximately independent of the initial filling factor of the system followed by a phase in which the fraction of retrograde orbits reaches a plateau and remains approximately constant and equal to about 0.52-0.54 until the end of the simulations. The fraction of retrograde orbits at the plateau is slightly larger and reached later in more underfilling systems (see e.g. in Fig. 2 the slight differences in evolution of the fraction of retrograde orbits in the KF05, KF025, and KF0125 models).
On the other hand, the unlocked models begin their evolution with a larger fraction of retrograde orbiting particles, and the time evolution of the total fraction of retrograde orbits is less significant, especially for the more filling models (see model KF1U in Fig. 2).
When the evolution of the global fraction of retrograde orbits is represented as a function of the fraction of the initial mass remaining in the system (see Fig. 3), some interesting aspects of the dynamics behind the evolution of the fraction of prograde and retrograde stars emerge. Initially locked, tidally filling models require a significant mass loss in order to show an appreciable growth of the fraction of retrograde orbits (see models KF1 and KF075, denoted with orange and green lines in Fig. 3), which therefore may be interpreted as associated with the loss of stars on prograde orbits.
But we emphasize here that in initially locked underfilling models, the initial evolution of the fraction of prograde and retrograde orbits is not associated with mass loss. This is particularly evident for the KF0125 and the KF025 models (depicted with black and red lines in Fig. 3), and indicates that for these systems there must be an additional internal mechanism affecting the fraction of prograde and retrograde stars before the preferential escape of prograde stars becomes important. We will now provide evidence that such a mechanism is indeed tied to the development of a halo of stars on eccentric/radial orbits driven by the effects of two-body relaxation.
As discussed in Tiongco et al. (2016) and in a number of previous studies (see e.g. Giersz & Heggie, 1994, 1997, 2011) initially isotropic star clusters that are either isolated or underfilling, develop during their long-term evolution an outer halo characterized by a radially anisotropic velocity distribution.
As a particle moves outwards on a highly eccentric orbit, angular momentum conservation (in a non-rotating frame of reference) implies that its rotational velocity around the z-axis slows down, and therefore the particle’s rotational velocity tends to become retrograde with respect to the rotating frame.
By inspecting a zoomed-in version of the global fraction of retrograde particles in, for example, the most underfilling system KF0125 (see Fig. 4, top panel), the growth of the fraction of retrograde particles is particularly evident, in the absence of significant mass loss. When compared to the growth of the radial velocity anisotropy in the outer regions (> 3, where is the 2D projected half-mass radius; see Fig. 4, bottom panel), the timescale of the increase appears to be consistent. To measure the degree of anisotropy we use the ratio of the projected tangential to the radial velocity dispersion, , measured on the cluster orbital plane.
We note that a similar interpretation applies also to unlocked systems, in which the difference between the behaviour of filling and underfilling models is less prominent (especially because, as discussed above, their initial fraction of retrograde orbits is initially higher compared to the corresponding locked systems), but still appreciable.
To provide a more quantitative characterization of the dynamics of the stars populating the outer regions of the star clusters, we explore the evolution of the models in the plane , where is the component of the angular momentum calculated in the rotating frame of reference and is projected distance from the cluster centre measured on the plane (the analysis performed on model KF05 is illustrated in Fig. 5). As particles start to populate the halo, they do so mostly with retrograde velocities and the outer regions are dominated by retrograde particles. A few trajectories on this plane determined by assuming angular momentum conservation are also plotted to further illustrate this effect.
3.3 The radial variation of the fraction of prograde and retrograde orbits and the origin of the plateau
While in the previous section we have explored the evolution of the fraction of retrograde orbits mostly from a global perspective, we now wish to investigate how the orbital distribution evolves within different regions of a tidally perturbed star cluster.
We first consider the filling system KF1 (see Fig. 6), by exploring the evolution of the fraction of retrograde stars within several radial shells, again using cylindrical shells with heights encompassing the entire cluster (notice that the fraction of retrograde particles for each shell is calculated with respect to the number of stars within each shell and not the total number of stars within the system). Indeed, the outermost shells show the strongest signature of asymmetry between the prograde and retrograde populations while in the innermost shells (), the fraction of retrograde particles remains approximately constant and equal to its initial value. The growth of the fraction retrograde orbits in the intermediate to outer shells is characterized by an early growth followed by a late flattening similar to that already seen in Fig. 2 for the evolution of the global fraction of retrograde orbits. For this model, the evolution of the fraction of retrograde particles is therefore dominated by the preferential loss of prograde particles although, as we will further discuss in Section 3.4, the additional mechanism tied to the development of radial anisotropy is affecting also this model.
As for star clusters in the underfilling regime, we perform the same analysis on the representative model KF025 (see Fig. 7). In this case the outermost regions near the cluster Jacobi radius are initially empty, and, as illustrated in Fig. 3, the fraction of retrograde particles starts to increase well before any significant mass loss. As discussed in Section 3.2, we ascribe this initial increase in the fraction of retrograde particles to the development of a halo of particles on highly-eccentric/radial orbits. Fig. 7 further illustrates how retrograde particles largely dominate the cluster outermost regions. Afterwards, when the system starts to lose significant mass, the preferential loss of prograde stars becomes important.
This exploration of the radial dependence of the fraction of prograde and retrograde orbits allows us also to provide an interpretation of the origin of the plateau observed in the second phase of the evolution of the global fraction of retrograde stars (see Figs. 2 and 3). We suggest that the origin of such plateau lies in the balancing act between the mass loss due to the expansion of the outer layers across the tidal boundary and the continued “production” of retrograde orbits through the mechanism associated with the ejection of stars from the core to the halo on highly eccentric/radial orbits (a process further enhanced in the post-core-collapse).
This interpretation seems to be well supported by the behaviour of the intermediate to outer shells of our models, which show a moderate increase of the fraction of retrograde orbits, followed by a stabilization or even a mild decline, especially for the last shell of initially underfilled systems (see purple, blue, and orange lines in Figs. 6 and 7). Additional evidence of the role of the tidal stripping of a cluster outer layers in the evolution of cluster anisotropy and structural properties has been discussed by Giersz & Heggie (1997).
3.4 Fraction of prograde and retrograde orbits for bound and unbound particles
As studied in detail by Fukushige & Heggie (2000), an unbound particle still within a cluster Jacobi radius can take many dynamical times before actually escaping the cluster. This can lead to the formation of a population of potential escapers residing within the cluster and affecting the cluster structural and kinematical properties (Küpper et al., 2010). In this section we further analyse our results in the context of energetically bound or unbound particles. In the frame co-rotating with a star cluster on a circular orbit in a point-mass galaxy, the energy of a star is given by (see e.g. Fukushige & Heggie, 2000)
where is the speed of the star and is the potential of the cluster. Following Fukushige & Heggie (2000) we define the critical energy as the potential at the Lagrangian points
and identify bound or unbound stars as those having an energy, respectively, below or above, .
We find that potential escapers may dominate the outermost regions of clusters. For example, Fig. 8 shows a snapshot of the plane for model KF05, with the potential escapers highlighted in green. It is evident that almost all of the outer cluster stars beyond approximately 0.5 are potential escapers. This result is in general agreement with the findings of Küpper et al. (2010), who first showed that the behaviour of the surface brightness and velocity dispersion radial profiles are significantly affected by the presence of an outer population of potential escapers.
The fraction of unbound particles as a function of the fraction of the mass remaining in the cluster for all our locked models is illustrated in Fig. 9. As discussed in Fukushige & Heggie (2000) and Baumgardt (2001), initial conditions set by sampling, for example, a King model may include a fraction of stars which are unbound when an external tidal field is added. This is the case for our KF1 model. For this model we show the evolution of the fraction of potential escapers both with and without this population of primordial potential escapers. All the other models do not have any primordial potential escapers.
Fig. 9 also shows that the fraction of potential escapers increases significantly over time reaching a fraction of the remaining mass equal to about 15% by the end of our simulations. Such a behaviour is qualitatively consistent with the findings reported by Baumgardt (2001), in which tidally perturbed models (initialised as critically filling King model , i.e. characterized by a much higher fraction of primordial potential escapers, about 15 per cent) possess approximately 20 per cent of potential escapers during the final stages of evolution.
We have examined the evolution of the fraction of prograde and retrograde particles within the populations of bound and unbound particles in models KF1 and KF025, as representative of the two filling regimes (see, Figs. 10 and 11, respectively). The larger fraction of retrograde particles in the population of unbound particles is the result of the two effects discussed in Section 3.2: the preferential escape of prograde particles and the fact that, in the transition from bound to unbound, some particles move on more eccentric/radial orbits becoming retrograde as they move outwards.
It is interesting that also the bound population shows an asymmetry in the fractions of prograde and retrograde particles. The growth in the fraction of retrograde particles in this case must be due to particles that develop more eccentric orbits while remaining bound to the cluster. The development of an excess of radially-biased orbits (as measured by the presence of radial anisotropy, see bottom panel of Fig. 4) in the bound population is therefore the cause of the increase in the fraction of retrograde particles in this population. As discussed in Tiongco et al. (2016), even the most filling King model (KF1) was characterized after core collapse by the development of radial anisotropy (although very weak) in the inner regions near the half-mass radius (see Fig. 10 in Tiongco et al. 2016) and indeed, Fig. 10 shows that for the KF1 model the slight increase in the fraction of retrograde particles starts only after core collapse at .
We have studied the effects of the host galaxy tidal field and star escape on the long-term evolution of a star cluster’s orbital distribution and internal kinematics. In particular, we have focused our attention on the evolution in the number of stars on prograde and retrograde orbits with respect to the cluster’s orbit about the galaxy and the implications for the evolution of the cluster internal rotation. We have carried out a number of -body simulations using isotropic King (1966) models varying the initial filling factor, and experimenting with models without rotation in the standard frame of reference co-rotating with an angular velocity equal to that of the cluster around the host galaxy, i.e. “locked", as well as with “unlocked” systems.
Previous studies investigating the restricted three-body problem of a star orbiting in the potential of the cluster and that of the cluster host galaxy have shown that unbound stars on retrograde orbits are more stable against escape, and predicted that a cluster might acquire retrograde rotation through preferential escape of stars on prograde orbits (Keenan & Innanen, 1975). In this study, we present evidence of this prediction in action and quantify its effect: the rotation curve of all our models evolves towards a retrograde rotational velocity beginning with the outer regions and moving inward, and the cluster’s fraction of retrograde orbiting particles slightly increases over time.
However, we have shown there are additional effects that cannot be accounted for by preferential escape of prograde orbits alone. In the early evolution, initially underfilling models increase their fraction of retrograde stars without losing significant mass and acquire a retrograde angular velocity. We attribute this effect to the development of preferentially eccentric/radial orbits in the outer regions of star clusters as they are expanding into their tidal limitation. Such an excess of elongated orbits can be measured in terms of the presence of a radially-biased anisotropy in the velocity space (see e.g. Giersz & Heggie, 1994, 1997, 2011; Tiongco et al., 2016; Zocchi et al., 2016). As a particle moves outwards on a highly eccentric orbit, angular momentum conservation (in a non-rotating frame of reference) implies that its rotational velocity around the -axis slows down and therefore that the particle rotational velocity tends to become retrograde with respect to the rotating frame. We have provided evidence that particles on these eccentric orbits populate the cluster halo mainly on retrograde orbits. We have also explored the time evolution of the radial variation of the fraction of retrograde particles and shown that, as the cluster evolves, the fraction of retrograde orbits increases with the distance from the cluster centre and the outer regions are increasingly populated by stars on retrograde orbits.
It is important to emphasize that once an underfilling cluster has expanded into its tidal limitation and is able to lose particles more rapidly, the preferential escape of prograde orbits becomes important and these models behave similarly to the initially filling models.
For all initially locked models, we observe that the rate at which the fraction of retrograde stars increases slows down until it stops and the fraction remains almost constant and equal to about 0.52-0.54 for the rest of the simulation. We have interpreted such a plateau as the result of the balance between the mass loss associated with the tidal stripping of the cluster outer layers expanding across the tidal boundary and the continued generation of retrograde orbits through the mechanism associated with the ejection of stars from the core to halo on preferentially highly eccentric/radial orbits.
The behaviour of the unlocked models considered in this study may be interpreted in terms of the same mechanisms illustrated above, although their evolution start with slightly larger fractions of retrograde orbits, which eventually converge to values similar to those reached in the simulations of non-rotating systems.
The interplay between the evolution of the fraction of prograde and retrograde stars and the effects of the tidal field leads to a rotational velocity radial profile with an angular velocity approximately constant and equal to about half of the angular velocity of the cluster orbital motion around the centre of the host galaxy. Although all the systems studied eventually reach this profile after losing a large fraction of their initial mass, the evolutionary history of the angular velocity radial profile depends on the cluster initial structural and kinematical properties.
Finally, we have explored the orbital properties of energetically bound and unbound particles. The fraction of unbound particles within the cluster (potential escapers) increases with time and by the end of our simulations (when the systems have lost 75% of their initial mass) they comprise almost 15% of the remaining mass (see also Baumgardt, 2001; Küpper et al., 2010) and dominate the outer regions beyond approximately 0.5 . Unbound particles have a higher fraction of retrograde particles when compared to the fraction calculated for bound particles or the whole system. Interestingly, also within the population of bound particles the retrograde fraction slightly increases with time; we ascribe this evolution to the increase of elongated orbits (i.e., an increase of radial anisotropy) in the bound population.
In the future we will investigate the behaviour of models with higher degree of rotation, and thus beginning with a stronger difference in the fraction of prograde and retrograde orbits.
Our investigation therefore further illustrates the complex effects of internal and external dynamical processes on the cluster phase-space density distribution. As observations from ESO/VLT, HST and the Gaia mission shed light on the internal kinematical properties of star clusters, theoretical efforts aimed at building a framework for the interpretation of these results must include the development of new families of dynamical models including features produced by dynamical processes such as different fractions of prograde/retrograde orbits and the presence of a significant fraction of unbound potential escapers (see Daniel, Heggie & Varri in prep).
- Aarseth (2003) Aarseth S. J., 2003, Gravitational N-Body Simulations. Cambridge Univ. Press, Cambridge
- Baumgardt (2001) Baumgardt H., 2001, MNRAS, 325, 1323
- Baumgardt & Makino (2003) Baumgardt H., Makino J., 2003, MNRAS, 340, 227
- Bellini et al. (2014) Bellini A., et al., 2014, ApJ, 797, 115
- Bellini et al. (2015) Bellini A., et al., 2015, ApJ, 810, L13
- Chandrasekhar (1942) Chandrasekhar S., 1942, Principles of stellar dynamics. Univ. Chicago Press, Chicago, IL
- Einsel & Spurzem (1999) Einsel C., Spurzem R., 1999, MNRAS, 302, 81
- Ernst et al. (2007) Ernst A., Glaschke P., Fiestas J., Just A., Spurzem R., 2007, MNRAS, 377, 465
- Fukushige & Heggie (2000) Fukushige T., Heggie D. C., 2000, MNRAS, 318, 753
- Gajda & Łokas (2016) Gajda G., Łokas E. L., 2016, ApJ, 819, 20
- Giersz & Heggie (1994) Giersz M., Heggie D. C., 1994, MNRAS, 268, 257
- Giersz & Heggie (1997) Giersz M., Heggie D. C., 1997, MNRAS, 286, 709
- Giersz & Heggie (2011) Giersz M., Heggie D. C., 2011, MNRAS, 410, 2698
- Goodman (1983) Goodman J. J., 1983, PhD thesis, Princeton Univ., NJ.
- Haghi et al. (2014) Haghi H., Hoseini-Rad S. M., Zonoozi A. H., Küpper A. H. W., 2014, MNRAS, 444, 3699
- Heggie & Hut (2003) Heggie D., Hut P., 2003, The Gravitational Million-Body Problem: A Multidisciplinary Approach to Star Cluster Dynamics. Cambridge Univ. Press, Cambridge
- Hong et al. (2013) Hong J., Kim E., Lee H. M., Spurzem R., 2013, MNRAS, 430, 2960
- Keenan & Innanen (1975) Keenan D. W., Innanen K. A., 1975, AJ, 80, 290
- King (1966) King I. R., 1966, AJ, 71, 64
- Küpper et al. (2010) Küpper A. H. W., Kroupa P., Baumgardt H., Heggie D. C., 2010, MNRAS, 407, 2241
- Lanzoni et al. (2013) Lanzoni B., et al., 2013, ApJ, 769, 107
- Lapenna et al. (2015) Lapenna E., Origlia L., Mucciarelli A., Lanzoni B., Ferraro F. R., Dalessandro E., Valenti E., Cirasuolo M., 2015, ApJ, 798, 23
- Lardo et al. (2015) Lardo C., et al., 2015, A&A, 573, A115
- Madrid et al. (2012) Madrid J. P., Hurley J. R., Sippel A. C., 2012, ApJ, 756, 167
- Nitadori & Aarseth (2012) Nitadori K., Aarseth S. J., 2012, MNRAS, 424, 545
- Read et al. (2006) Read J. I., Wilkinson M. I., Evans N. W., Gilmore G., Kleyna J. T., 2006, MNRAS, 367, 387
- Richer et al. (2013) Richer H. B., Heyl J., Anderson J., Kalirai J. S., Shara M. M., Dotter A., Fahlman G. G., Rich R. M., 2013, ApJ, 771, L15
- Sollima et al. (2015) Sollima A., Baumgardt H., Zocchi A., Balbinot E., Gieles M., Hénault-Brunet V., Varri A. L., 2015, MNRAS, 451, 2185
- Tiongco et al. (2016) Tiongco M. A., Vesperini E., Varri A. L., 2016, MNRAS, 455, 3693
- Vesperini & Heggie (1997) Vesperini E., Heggie D. C., 1997, MNRAS, 289, 898
- Watkins et al. (2015) Watkins L. L., van der Marel R. P., Bellini A., Anderson J., 2015, ApJ, 803, 29
- Webb & Leigh (2015) Webb J. J., Leigh N. W. C., 2015, MNRAS, 453, 3278
- Webb et al. (2014) Webb J. J., Leigh N., Sills A., Harris W. E., Hurley J. R., 2014, MNRAS, 442, 1569
- Zocchi et al. (2016) Zocchi A., Gieles M., Hénault-Brunet V., Varri A. L., 2016, MNRAS, in press
- Zonoozi et al. (2011) Zonoozi A. H., Küpper A. H. W., Baumgardt H., Haghi H., Kroupa P., Hilker M., 2011, MNRAS, 411, 1989
- Zonoozi et al. (2014) Zonoozi A. H., Haghi H., Küpper A. H. W., Baumgardt H., Frank M. J., Kroupa P., 2014, MNRAS, 440, 3172
- Zonoozi et al. (2016) Zonoozi A. H., Rabiee M., Haghi H., Küpper A. H. W., 2016, ApJ, 818, 58
- Zotos (2015) Zotos E. E., 2015, MNRAS, 446, 770
This research was supported in part by Lilly Endowment, Inc., through its support for the Indiana University Pervasive Technology Institute, and in part by the Indiana METACyt Initiative. The Indiana METACyt Initiative at IU is also supported in part by Lilly Endowment, Inc. ALV is grateful to Douglas Heggie for stimulating conversations, and to the EU Horizon 2020 program for support in form of a Marie Sklodowska-Curie Fellowship (MSCA-IF-EF-RI 658088). We wish to thank the referee for a careful report.