Disks of low-mass bodies on high-eccentricity orbits in near-Keplerian potentials can be dynamically unstable to buckling out of the plane. In this letter, we present $N$-body simulations of the long-term behavior of such a system, finding apsidal clustering of the orbits in the disk plane. The timescale over which the clustering is maintained increases with number of particles, suggesting that lopsided configurations are stable at large $N$. This discovery may explain the observed apsidal ($\varpi$) clustering of extreme trans-Neptunian Objects in the outer solar system.

We present the results of a suite of \Nbody simulations aimed at understanding the fundamental aspects of the long-term evolution of the internal kinematics of multiple stellar populations in globular clusters. Our models enable us to study the cooperative effects of internal, relaxation-driven processes and external, tidally-induced perturbations on the structural and kinematic properties of multiple-population globular clusters. To analyse the dynamical behaviour of the multiple stellar populations in a variety of spin-orbit coupling conditions, we have considered three reference cases in which the tidally perturbed star cluster rotates along an axis oriented in different directions with respect to the orbital angular momentum vector. We focus specifically on the characterisation of the evolution of the degree of differential rotation and anisotropy in the velocity space, and we quantify the process of spatial and kinematic mixing of the two populations. In light of recent and forthcoming explorations of the internal kinematics of this class of stellar systems by means of line-of sight and astrometric measurements, we also investigate the implications of projection effects and spatial distribution of the stars adopted as tracers. The kinematic and structural richness emerging from our models further emphasises the need and the importance of observational studies aimed at building a complete kinematical picture of the multiple population phenomenon.

As part of the ESO-VLT Multi-Instrument Kinematic Survey (MIKiS) of Galactic globular clusters, we present a detailed investigation of the internal kinematics of NGC 5986. The analysis is based on about 300 individual radial velocities of stars located at various distances from the cluster center, up to 300 arcseconds (about 4 half-mass radii). Our analysis reveals the presence of a solid-body rotation extending from the cluster center to the outermost regions probed by the data, and a velocity dispersion profile initially declining with the distance from the cluster's center, but flattening and staying constant at ~5 km/s for distances larger than about one half-mass radius. This is the first globular cluster for which evidence of the joint presence of solid-body rotation and flattening in the outer velocity dispersion profile is found. The combination of these two kinematical features provides a unique opportunity to shed light on fundamental aspects of globular cluster dynamics and probe the extent to which internal relaxation, star escape, angular momentum transport and loss, and the interaction with the Galaxy tidal field can affect a cluster's dynamical evolution and determine its current kinematical properties. We present the results of a series of N-body simulations illustrating the possible dynamical paths leading to kinematic features like those observed in this cluster and the fundamental dynamical processes that underpin them.

In the context of the ESO-VLT Multi-Instrument Kinematic Survey (MIKiS) of Galactic globular clusters, we present the line-of-sight rotation curve and velocity dispersion profile of M5 (NGC 5904), as determined from the radial velocity of more than 800 individual stars observed out to 700" (~ 5 half-mass radii) from the center. We find one of the cleanest and most coherent rotation patterns ever observed for globular clusters, with a very stable rotation axis (having constant position angle of 145^o at all surveyed radii) and a well-defined rotation curve. The density distribution turns out to be flattened in the direction perpendicular to the rotation axis, with a maximum ellipticity of 0.15. The rotation velocity peak (~3 km/s in projection) is observed at ~0.6 half-mass radii, and its ratio with respect to the central velocity dispersion (~0.3-0.4 at 4 projected half-mass radii) indicates that ordered motions play a significant dynamical role. This result strengthens the growing empirical evidence of the kinematic complexity of Galactic globular clusters and motivates the need of fundamental investigations of the role of angular momentum in collisional stellar dynamics.

We broaden the investigation of the dynamical properties of tidally perturbed, rotating star clusters by relaxing the traditional assumptions of coplanarity, alignment, and synchronicity between the internal and orbital angular velocity vector of their initial conditions. We show that the interplay between the internal evolution of these systems and their interaction with the external tidal field naturally leads to the development of a number of evolutionary features in their three-dimensional velocity space, including a precession and nutation of the global rotation axis and a variation of its orientation with the distance from the cluster centre. In some cases, such a radial variation may manifest itself as a counter-rotation of the outermost regions relative to the inner ones. The projected morphology of these systems is characterized by a non-monotonic ellipticity profile and, depending on the initial inclination of the rotation axis, it may also show a twisting of the projected isodensity contours. These results provide guidance in the identification of non-trivial features which may emerge in upcoming investigations of star cluster kinematics and a dynamical framework to understand some of the complexities already hinted by recent observational studies.

We present the results of a set of N-body simulations following the long-term evolution of the rotational properties of star cluster models evolving in the external tidal field of their host galaxy, after an initial phase of violent relaxation. The effects of two-body relaxation and escape of stars lead to a redistribution of the ordered kinetic energy from the inner to the outer regions, ultimately determining a progressive general loss of angular momentum; these effects are reflected in the overall decline the rotation curve as the cluster evolves and loses stars. We show that all of our models share the same dependence of the remaining fraction of the initial rotation on the fraction of the initial mass lost. As the cluster evolves and loses part of its initial angular momentum, it becomes increasingly dominated by random motions, but even after several tens of relaxation times, and losing a significant fraction of its initial mass, a cluster can still be characterized by a non-negligible ratio of the rotational velocity to the velocity dispersion. This result is in qualitative agreement with the recently observed kinematical complexity which characterizes several Galactic globular clusters.

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.

We explore the long-term evolution of the anisotropy in the velocity space of star clusters starting with different structural and kinematical properties. We show that the evolution of the radial anisotropy strength and its radial variation within a cluster contain distinct imprints of the cluster initial structural properties, dynamical history, and of the external tidal field of its host galaxy. Initially isotropic and compact clusters with small initial values of the ratio of the half-mass to Jacobi radius, $r_h/r_J$, develop a strong radial anisotropy during their long-term dynamical evolution. Many clusters, if formed with small values of $r_h/r_J$, should now be characterized by a significant radial anisotropy increasing with the distance from the cluster centre, reaching its maximum at a distance between 0.2 $r_J$ and 0.4 $r_J$, and then becoming more isotropic or mildly tangentially anisotropic in the outermost regions. A similar radial variation of the anisotropy can also result from an early violent relaxation phase. In both cases, as a cluster continues its evolution and loses mass, the anisotropy eventually starts to decrease and the system evolves toward an isotropic velocity distribution. However, in order to completely erase the strong anisotropy developed by these compact systems during their evolution, they must be in the advanced stages of their evolution and lose a large fraction of their initial mass. Clusters that are initially isotropic and characterized by larger initial values of $r_h/r_J$, on the other hand, never develop a significant radial anisotropy.