Showing 1–42 of 42 results for author: Varri, A L

Dynamics of star clusters with tangentially anisotropic velocity distribution

Authors: Václav Pavlík, Douglas C. Heggie, Anna Lisa Varri, Enrico Vesperini

Abstract: Recent high-precision observations with HST and Gaia enabled new investigations of the internal kinematics of star clusters (SCs) and the dependence of kinematic properties on the stellar mass. These studies raised new questions about the dynamical evolution of self-gravitating stellar systems. We aim to develop a more complete theoretical understanding of how various kinematical properties of sta… ▽ More Recent high-precision observations with HST and Gaia enabled new investigations of the internal kinematics of star clusters (SCs) and the dependence of kinematic properties on the stellar mass. These studies raised new questions about the dynamical evolution of self-gravitating stellar systems. We aim to develop a more complete theoretical understanding of how various kinematical properties of stars affect the global dynamical development of their host SCs. We perform N-body simulations of SCs with isotropic, radially anisotropic and tangentially anisotropic initial velocity distributions. We also study the effect of an external Galactic tidal field. We find three main results. First, compared to the conventional, isotropic case, the relaxation processes are accelerated in the tangentially anisotropic models and, in agreement with our previous investigations, slower in the radially anisotropic ones. This leads to, e.g., more rapid mass segregation in the central regions of the tangential models or their earlier core collapse. Second, although all SCs become isotropic in the inner regions after several relaxation times, we observe differences in the anisotropy profile evolution in the outer cluster regions - all tidally filling models gain tangential anisotropy there while the underfilling models become radially anisotropic. Third, we observe different rates of evolution towards energy equipartition (EEP). While all SCs evolve towards EEP in their inner regions (regardless of the filling factor), the outer regions of the tangentially anisotropic and isotropic models are evolving to an "inverted" EEP (i.e., the high-mass stars having higher velocity dispersion than the low-mass ones). The extent (both spatial and temporal) of this inversion can be attributed to the initial velocity anisotropy - it grows with increasing tangential anisotropy and decreases as the radial anisotropy rises. △ Less

Submitted 29 May, 2024; originally announced May 2024.

Comments: 13 pages, 13 figures, 1 table, 1 appendix. Accepted for publication in A&A

Ongoing hierarchical massive cluster assembly: the LISCA II structure in the Perseus complex

Authors: A. Della Croce, E. Dalessandro, A. Livernois, E. Vesperini, C. Fanelli, L. Origlia, M. Bellazzini, E. Oliva, N. Sanna, A. L. Varri

Abstract: We report on the identification of a massive ($\sim10^5$ M$_\odot$) sub-structured stellar system in the Galactic Perseus complex likely undergoing hierarchical cluster assembly. Such a system comprises nine star clusters (including the well-known clusters NGC 654 and NGC 663) and an extended and low-density stellar halo. Gaia-DR3 and available spectroscopic data show that all its components are p… ▽ More We report on the identification of a massive ($\sim10^5$ M$_\odot$) sub-structured stellar system in the Galactic Perseus complex likely undergoing hierarchical cluster assembly. Such a system comprises nine star clusters (including the well-known clusters NGC 654 and NGC 663) and an extended and low-density stellar halo. Gaia-DR3 and available spectroscopic data show that all its components are physically consistent in the 6D phase-space (position, parallax, and 3D motion), homogeneous in age (14 $-$ 44 Myr), and chemical content (half-solar metallicity). In addition, the system's global stellar density distribution is that of typical star clusters and shows clear evidence of mass segregation. We find that the hierarchical structure is mostly contracting towards the center with a speed of up to $\simeq4-5$ km s$^{-1}$, while the innermost regions expand at a lower rate (about $\simeq1$ km s$^{-1}$) and are dominated by random motions. Interestingly, this pattern is dominated by the kinematics of massive stars, while low-mass stars ($M<2$ M$_\odot$) are characterized by contraction across the entire cluster. Finally, the nine star clusters in the system are all characterized by a relatively flat velocity dispersion profile possibly resulting from ongoing interactions and tidal heating. We show that the observational results are generally consistent with those found in $N$-body simulations following the cluster violent relaxation phase strongly suggesting that the system is a massive cluster in the early assembly stages. This is the second structure with these properties identified in our Galaxy and, following the nomenclature of our previous work, we named it LISCA II. △ Less

Submitted 27 March, 2023; originally announced March 2023.

Comments: 21 pages, 24 figures, 1 table; accepted for publication in A&A

Journal ref: A&A 674, A93 (2023)

Linking the Internal Properties of Infant Globular Clusters to their Formation Environments

Authors: Frederika Phipps, Sadegh Khochfar, Anna Lisa Varri, Claudio Dalla Vecchia

Abstract: We investigate the formation of infant globular cluster (GC) candidates in high-resolution cosmological simulations from the First Billion Years (FiBY) project. By analysing the evolution of the systems in the energy and angular momentum plane, we identify the redshift at which the infant GCs first became gravitationally bound, and we find evidence of radial infall of their gaseous and stellar com… ▽ More We investigate the formation of infant globular cluster (GC) candidates in high-resolution cosmological simulations from the First Billion Years (FiBY) project. By analysing the evolution of the systems in the energy and angular momentum plane, we identify the redshift at which the infant GCs first became gravitationally bound, and we find evidence of radial infall of their gaseous and stellar components. The collapse appears to be driven by internal self-gravity, however, the initial trigger is sourced from the external environment. The phase space behaviour of the infant GCs also allows us to identify some characteristic groupings of objects. Such a classification based on internal properties appears to be reflected in the formation environment: GC candidates that belong to the same class are found in host galaxies of similar morphology, with the majority of the infant GCs located in clumpy, irregular proto-galaxies. Finally, through the inspection of two GC candidates that contain only stars by z = 6, we find that supernova feedback is the main physical mechanism behind their dearth of gas and that the systems subsequently respond with an approximately adiabatic expansion. Such infant GC candidates already resemble the GCs we currently observe in the local Universe. △ Less

Submitted 17 November, 2022; originally announced November 2022.

Comments: 17 pages, 12 Figures, Accepted by MNRAS

Long-term evolution of multimass rotating star clusters

Authors: Alexander R. Livernois, Enrico Vesperini, Anna Lisa Varri, Jongsuk Hong, Maria Tiongco

Abstract: We investigate the long-term dynamical evolution of the internal kinematics of multimass rotating star clusters. We have performed a set of N-body simulations to follow the internal evolution of clusters with different degrees of initial rotation and have explored the evolution of the rotational velocity, the degree of energy equipartition, and anisotropy in the velocity distribution. Our simulati… ▽ More We investigate the long-term dynamical evolution of the internal kinematics of multimass rotating star clusters. We have performed a set of N-body simulations to follow the internal evolution of clusters with different degrees of initial rotation and have explored the evolution of the rotational velocity, the degree of energy equipartition, and anisotropy in the velocity distribution. Our simulations show that: 1) as the cluster evolves, the rotational velocity develops a dependence on the stellar mass with more massive stars characterised by a more rapid rotation and a peak in the rotation curve closer to the cluster centre than low-mass stars; 2) the degree of energy equipartition in the cluster's intermediate and outer regions depends on the component of the velocity dispersion measured; for more rapidly rotating clusters, the evolution towards energy equipartition is more rapid in the direction of the rotational velocity; 3) the anisotropy in the velocity distribution is stronger for massive stars; 4) both the degree of mass segregation and energy equipartition are characterised by spatial anisotropy; they have a dependence on both $R$ and $z$, correlated with the flattening in the spatial variation of the cluster's density and velocity dispersion, as shown by 2D maps of the mass segregation and energy equipartition on the ($R$-$z$) meridional plane. △ Less

Submitted 13 April, 2022; originally announced April 2022.

Comments: 11 pages, 10 figures, Published in MNRAS

Journal ref: MNRAS Volume 512, Issue 2, May 2022, Pages 2584-2593

The ESO-VLT MIKiS survey reloaded: velocity dispersion profile and rotation curve of NGC 1904

Authors: S. Leanza, C. Pallanca, F. R. Ferraro, B. Lanzoni, E. Dalessandro, L. Origlia, A. Mucciarelli, E. Valenti, M. Tiongco, A. L. Varri, E. Vesperini

Abstract: We present an investigation of the internal kinematic properties of M79 (NGC 1904). Our study is based on radial velocity measurements obtained from the ESO-VLT Multi-Instrument Kinematic Survey (MIKiS) of Galactic globular clusters for more than 1700 individual stars distributed between $\sim 0.3^{\prime\prime}$ and $770^{\prime\prime}$ ($\sim14$ three-dimensional half-mass radii), from the cente… ▽ More We present an investigation of the internal kinematic properties of M79 (NGC 1904). Our study is based on radial velocity measurements obtained from the ESO-VLT Multi-Instrument Kinematic Survey (MIKiS) of Galactic globular clusters for more than 1700 individual stars distributed between $\sim 0.3^{\prime\prime}$ and $770^{\prime\prime}$ ($\sim14$ three-dimensional half-mass radii), from the center. Our analysis reveals the presence of ordered line-of-sight rotation with a rotation axis almost aligned along the East-West direction and a velocity peak of $1.5$ km s$^{-1}$ at $\sim 70^{\prime\prime}$ from the rotation axis. The velocity dispersion profile is well described by the same King model that best fits the projected density distribution, with a constant central plateau at $σ_0\sim 6$ km s$^{-1}$. To investigate the cluster rotation in the plane of the sky, we have analyzed the proper motions provided by the Gaia EDR3, finding a signature of rotation with a maximum amplitude of $\sim 2.0$ km s$^{-1}$ at $\sim 80^{\prime\prime}$ from the cluster center. Analyzing the three-dimensional velocity distribution, for a sub-sample of 130 stars, we confirm the presence of systemic rotation and find a rotation axis inclination angle of $37$° with respect to the line-of-sight. As a final result, the comparison of the observed rotation curves with the results of a representative N-body simulation of a rotating star cluster shows that the present-day kinematic properties of NGC 1904 are consistent with those of a dynamically old system that has lost a significant fraction of its initial angular momentum. △ Less

Submitted 14 March, 2022; originally announced March 2022.

Comments: 30 pages, 20 figures, 8 tables, accepted for publication in ApJ

Early dynamical evolution of rotating star clusters in a tidal field

Authors: Maria Tiongco, Enrico Vesperini, Anna Lisa Varri

Abstract: In order to explore how the early internal rotational properties of star clusters are affected by the external potential of their host galaxies, we have run a suite of $N$-body simulations following the early dynamical evolution and violent relaxation of rotating star clusters embedded in a tidal field. Our study focuses on models for which the cluster's rotation axis has a generic orientation rel… ▽ More In order to explore how the early internal rotational properties of star clusters are affected by the external potential of their host galaxies, we have run a suite of $N$-body simulations following the early dynamical evolution and violent relaxation of rotating star clusters embedded in a tidal field. Our study focuses on models for which the cluster's rotation axis has a generic orientation relative to the torque of the tidal field. The interaction between the violent relaxation process, angular momentum of the cluster, and the external torque creates a complex kinematic structure within the cluster, most prominently a radial variation in the position of the rotation axis along both the polar and azimuthal directions. We also examine the cluster's velocity dispersion anisotropy and show that the projected anisotropy may be affected by the variation of the rotation axis directions within the cluster; the combination of projection effects and the complex kinematical features may result in the measurement of tangential anisotropy in the cluster's inner regions. We also characterize the structural properties of our clusters as a function of their initial rotation and virial ratio and find that clusters may develop a triaxial morphology and a radial variation of the minor axis not necessarily aligned with the rotation axis. Finally, we examine the long-term evolution of these complex kinematic features. △ Less

Submitted 8 March, 2022; originally announced March 2022.

Comments: 15 pages, 12 figures, accepted by MNRAS

Forward and Back: Kinematics of the Palomar 5 Tidal Tails

Authors: Pete B. Kuzma, Annette M. N. Ferguson, Anna Lisa Varri, Michael J. Irwin, Edouard J. Bernard, Eline Tolstoy, Jorge Peñarrubia, Daniel B. Zucker

Abstract: The tidal tails of Palomar 5 (Pal 5) have been the focus of many spectroscopic studies in an attempt to identify individual stars lying along the stream and characterise their kinematics. The well-studied trailing tail has been explored out to a distance of 15^\text{o} from the cluster centre, while less than four degrees have been examined along the leading tail. In this paper, we present results… ▽ More The tidal tails of Palomar 5 (Pal 5) have been the focus of many spectroscopic studies in an attempt to identify individual stars lying along the stream and characterise their kinematics. The well-studied trailing tail has been explored out to a distance of 15^\text{o} from the cluster centre, while less than four degrees have been examined along the leading tail. In this paper, we present results of a spectroscopic study of two fields along the leading tail that we have observed with the AAOmega spectrograph on the Anglo-Australian telescope. One of these fields lies roughly 7^\text{o} along the leading tail, beyond what has been previously been explored spectroscopically. Combining our measurements of kinematics and line strengths with Pan-STARRS1 photometric data and Gaia EDR3 astrometry, we adopt a probabilistic approach to identify 16 stars with high probability of belonging to the Pal 5 stream. Eight of these stars lie in the outermost field and their sky positions confirm the presence of ``fanning'' in the leading arm. We also revisit previously-published radial velocity studies and incorporate Gaia EDR3 astrometry to remove interloping field stars. With a final sample of 109 {\it bona fide} Pal 5 cluster and tidal stream stars, we characterise the 3D kinematics along the the full extent of the system. We provide this catalogue for future modeling work. △ Less

Submitted 9 February, 2022; originally announced February 2022.

Comments: 13 pages, 8 figures, accepted for publication in MNRAS. The extended version of Table 4 is available as an ancillery file, and will be supplementary material in MNRAS

Early dynamics and violent relaxation of multi-mass rotating star clusters

Authors: A. R. Livernois, E. Vesperini, M. Tiongco, A. L. Varri, E. Dalessandro

Abstract: We present the results of a study aimed at exploring, by means of N-body simulations, the evolution of rotating multi-mass star clusters during the violent relaxation phase, in the presence of a weak external tidal field. We study the implications of the initial rotation and the presence of a mass spectrum for the violent relaxation dynamics and the final properties of the equilibria emerging at t… ▽ More We present the results of a study aimed at exploring, by means of N-body simulations, the evolution of rotating multi-mass star clusters during the violent relaxation phase, in the presence of a weak external tidal field. We study the implications of the initial rotation and the presence of a mass spectrum for the violent relaxation dynamics and the final properties of the equilibria emerging at the end of this stage. Our simulations show a clear manifestation of the evolution towards spatial mass segregation and evolution towards energy equipartition during and at the end of the violent relaxation phase. We study the final rotational kinematics and show that massive stars tend to rotate more rapidly than low-mass stars around the axis of cluster rotation. Our analysis also reveals that during the violent relaxation phase, massive stars tend to preferentially segregate into orbits with angular momentum aligned with the cluster's angular momentum, an effect previously found in the context of the long-term evolution of star clusters driven by two-body relaxation. △ Less

Submitted 23 July, 2021; originally announced July 2021.

Comments: 18 pages, 20 figures, 1 appendix. Accepted for publication in MNRAS

Central Dynamics of Multi-mass Rotating Star Clusters

Authors: Maria Tiongco, Angela Collier, Anna Lisa Varri

Abstract: We investigate the evolutionary nexus between the morphology and internal kinematics of the central regions of collisional, rotating, multi-mass stellar systems, with special attention to the spatial characterisation of the process of mass segregation. We report results from idealized, purely $N$-body simulations that show multi-mass, rotating, and spherical systems rapidly form an oblate, spheroi… ▽ More We investigate the evolutionary nexus between the morphology and internal kinematics of the central regions of collisional, rotating, multi-mass stellar systems, with special attention to the spatial characterisation of the process of mass segregation. We report results from idealized, purely $N$-body simulations that show multi-mass, rotating, and spherical systems rapidly form an oblate, spheroidal massive core, unlike single-mass rotating or multi-mass non-rotating configurations with otherwise identical initial properties, indicating that this evolution is a result of the interplay between the presence of a mass spectrum and angular momentum. This feature appears to be long-lasting, preserving itself for several relaxation times. The degree of flattening experienced by the systems is directly proportional to the initial degree of internal rotation. In addition, this morphological effect has a clear characterisation in terms of orbital architecture, as it lowers the inclination of the orbits of massive stars. We offer an idealised dynamical interpretation that could explain the mechanism underpinning this effect and we highlight possible useful implications, from kinematic hysteresis to spatial distribution of dark remnants in dense stellar systems. △ Less

Submitted 7 July, 2021; originally announced July 2021.

Comments: 13 pages, 12 figures, Accepted for publication in MNRAS

The Kinematic Richness of Star Clusters - II. Stability of Spherical Anisotropic Models with Rotation

Authors: Philip G. Breen, Simon Rozier, Douglas C. Heggie, Anna Lisa Varri

Abstract: We study the bar instability in collisionless, rotating, anisotropic, stellar systems, using N-body simulations and also the matrix technique for calculation of modes with the perturbed collisionless Boltzmann equation. These methods are applied to spherical systems with an initial Plummer density distribution, but modified kinematically in two ways: the velocity distribution is tangentially aniso… ▽ More We study the bar instability in collisionless, rotating, anisotropic, stellar systems, using N-body simulations and also the matrix technique for calculation of modes with the perturbed collisionless Boltzmann equation. These methods are applied to spherical systems with an initial Plummer density distribution, but modified kinematically in two ways: the velocity distribution is tangentially anisotropic, using results of Dejonghe, and the system is set in rotation by reversing the velocities of a fraction of stars in various regions of phase space, a la Lynden-Bell. The aim of the N-body simulations is first to survey the parameter space, and, using those results, to identify regions of phase space (by radius and orbital inclination) which have the most important influence on the bar instability. The matrix method is then used to identify the resonant interactions in the system which have the greatest effect on the growth rate of a bar. Complementary series of N-body simulations examine these processes in relation to the evolving frequency distribution and the pattern speed. Finally, the results are synthesised with an existing theoretical framework, and used to consider the old question of constructing a stability criterion. △ Less

Submitted 8 February, 2021; originally announced February 2021.

Comments: 17 pages, 18 figures, accepted for publication in MNRAS

First phase space portrait of a hierarchical stellar structure in the Milky Way

Authors: E. Dalessandro, A. L. Varri, M. Tiongco, E. Vesperini, C. Fanelli, A. Mucciarelli, L. Origlia, M. Bellazzini, S. Saracino, E. Oliva, N. Sanna, M. Fabrizio, A. Livernois

Abstract: We present the first detailed observational picture of a possible ongoing massive cluster hierarchical assembly in the Galactic disk as revealed by the analysis of the stellar full phase-space (3D positions and kinematics and spectro-photometric properties) of an extended area ($6^{\circ}$ diameter) surrounding the well-known $\it h$ and $χ$ Persei double stellar cluster in the Perseus Arm. Gaia-E… ▽ More We present the first detailed observational picture of a possible ongoing massive cluster hierarchical assembly in the Galactic disk as revealed by the analysis of the stellar full phase-space (3D positions and kinematics and spectro-photometric properties) of an extended area ($6^{\circ}$ diameter) surrounding the well-known $\it h$ and $χ$ Persei double stellar cluster in the Perseus Arm. Gaia-EDR3 shows that the area is populated by seven co-moving clusters, three of which were previously unknown, and by an extended and quite massive ($M\sim10^5 M_{\odot}$) halo. All stars and clusters define a complex structure with evidence of possible mutual interactions in the form of intra-cluster over-densities and/or bridges. They share the same chemical abundances (half-solar metallicity) and age ($t\sim20$ Myr) within a small confidence interval and the stellar density distribution of the surrounding diffuse stellar halo resembles that of a cluster-like stellar system. The combination of these evidences suggests that stars distributed within a few degrees from $\it h$ and $χ$ Persei are part of a common, sub-structured stellar complex that we named LISCA I. Comparison with results obtained through direct $N$-body simulations suggest that LISCA I may be at an intermediate stage of an ongoing cluster assembly that can eventually evolve in a relatively massive (a few $10^5 M_{\odot}$) stellar system. We argue that such cluster formation mechanism may be quite efficient in the Milky Way and disk-like galaxies and, as a consequence, it has a relevant impact on our understanding of cluster formation efficiency as a function of the environment and redshift. △ Less

Submitted 11 January, 2021; originally announced January 2021.

Comments: 19 pages, 8 figures, 1 table; accepted for publication in ApJ

$l=1$: Weinberg's weakly damped mode in an $N$-body model of a spherical stellar system

Authors: Douglas C. Heggie, Philip G. Breen, Anna Lisa Varri

Abstract: Spherical stellar systems such as King models, in which the distribution function is a decreasing function of energy and depends on no other invariant, are stable in the sense of collisionless dynamics. But Weinberg showed, by a clever application of the matrix method of linear stability, that they may be nearly unstable, in the sense of possessing {\sl weakly} damped modes of oscillation. He also… ▽ More Spherical stellar systems such as King models, in which the distribution function is a decreasing function of energy and depends on no other invariant, are stable in the sense of collisionless dynamics. But Weinberg showed, by a clever application of the matrix method of linear stability, that they may be nearly unstable, in the sense of possessing {\sl weakly} damped modes of oscillation. He also demonstrated the presence of such a mode in an $N$-body model by endowing it with initial conditions generated from his perturbative solution. In the present paper we provide evidence for the presence of this same mode in $N$-body simulations of the King $W_0 = 5$ model, in which the initial conditions are generated by the usual Monte Carlo sampling of the King distribution function. It is shown that the oscillation of the density centre correlates with variations in the structure of the system out to a radius of about 1 virial radius, but anticorrelates with variations beyond that radius. Though the oscillations appear to be continually reexcited (presumably by the motions of the particles) we show by calculation of power spectra that Weinberg's estimate of the period (strictly, $2π$ divided by the real part of the eigenfrequency) lies within the range where the power is largest. In addition, however, the power spectrum displays another very prominent feature at shorter periods, around 5 crossing times. △ Less

Submitted 2 March, 2021; v1 submitted 3 February, 2020; originally announced February 2020.

Comments: 9 pages, 6 figures, footnote added in Sec.3.2, no other change to the version published in MNRAS

The First Billion Years Project: Finding Infant Globular Clusters at z=6

Authors: Frederika Phipps, Sadegh Khochfar, Anna Lisa Varri, Claudio Dalla Vecchia

Abstract: We explored a suite of high-resolution cosmological simulations from the First Billion Years Project (FiBY) at $z \geq 6$. All substructures within the simulations have been identified with the SUBFIND algorithm. From our analysis, two distinct groups of objects emerge. We hypothesise that the substructures in the first group, which appear to have a high baryon fraction ($f_{\rm b} \geq 0.95$), ar… ▽ More We explored a suite of high-resolution cosmological simulations from the First Billion Years Project (FiBY) at $z \geq 6$. All substructures within the simulations have been identified with the SUBFIND algorithm. From our analysis, two distinct groups of objects emerge. We hypothesise that the substructures in the first group, which appear to have a high baryon fraction ($f_{\rm b} \geq 0.95$), are possible infant GC candidates. Objects belonging to the second group have a high stellar fraction ($f_{\rm star} \geq 0.95$) and show a potential resemblance to infant ultra-faint dwarf galaxies. The high baryon fraction objects identified in this study are characterised by a stellar content similar to the one observed in present-day GCs, but they still contain a high gas fraction ($f_{\rm gas} \sim 0.95$) and a relatively low amount of dark matter. They are compact, dense systems. Their sizes are consistent with recent estimates based on the first observations of possible proto-GCs at high redshifts. These types of infant GC candidates appear to be more massive and more abundant in massive host galaxies, indicating that the assembly of galaxies via mergers may play an important role in building several GC-host scaling relations. Specifically, we express the relation between the mass of the most massive infant GC and its host stellar mass as $\log(M_{\rm cl}) = (0.31\pm0.15)\log(M_{\rm *,gal} + (4.17\pm1.06)$. We also report a new relation between the most massive infant GC and the parent specific star formation rate of the form $\log(M_{\rm cl}) = (0.85\pm0.30)\log(sSFR) + α$ that describes the data at both low and high redshift. Finally, we assess the present-day GC mass (GC number) -- halo mass relation offers a satisfactory description of the behaviour of our infant GC candidates at high redshift, suggesting that such a relation may be set at formation. △ Less

Submitted 3 July, 2020; v1 submitted 22 October, 2019; originally announced October 2019.

Comments: 17 pages, 11 figures, accepted by A&A

Journal ref: A&A 641, A132 (2020)

Kinematical evolution of multiple stellar populations in star clusters

Authors: Maria Tiongco, Enrico Vesperini, Anna Lisa Varri

Abstract: 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 mult… ▽ More 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. △ Less

Submitted 12 July, 2019; originally announced July 2019.

Comments: 14 pages, 12 figures, Accepted by MNRAS

Spherical models of star clusters with potential escapers

Authors: Ian Claydon, Mark Gieles, Anna Lisa Varri, Douglas C. Heggie, Alice Zocchi

Abstract: An increasing number of observations of the outer regions of globular clusters (GCs) have shown a flattening of the velocity dispersion profile and an extended surface density profile. Formation scenarios of GCs can lead to different explanations of these peculiarities, therefore the dynamics of stars in the outskirts of GCs are an important tool in tracing back the evolutionary history and format… ▽ More An increasing number of observations of the outer regions of globular clusters (GCs) have shown a flattening of the velocity dispersion profile and an extended surface density profile. Formation scenarios of GCs can lead to different explanations of these peculiarities, therefore the dynamics of stars in the outskirts of GCs are an important tool in tracing back the evolutionary history and formation of star clusters. One possible explanation for these features is that GCs are embedded in dark matter halos. Alternatively, these features are the result of a population of energetically unbound stars that can be spatially trapped within the cluster, known as potential escapers (PEs). We present a prescription for the contribution of these energetically unbound members to a family of self-consistent, distribution function-based models, which, for brevity, we call the Spherical Potential Escapers Stitched (SPES) models. We show that, when fitting to mock data of bound and unbound stars from an N-body model of a tidally-limited star cluster, the SPES models correctly reproduce the density and velocity dispersion profiles up to the Jacobi radius, and they are able to recover the value of the Jacobi radius itself to within 20%. We also provide a comparison to the number density and velocity dispersion profiles of the Galactic cluster 47 Tucanae. Such a case offers a proof of concept that an appropriate modelling of PEs is essential to accurately interpret Gaia data in the outskirts of GCs, and, in turn, to formulate meaningful present-day constraints for GC formation scenarios in the early universe. △ Less

Submitted 15 March, 2019; v1 submitted 14 March, 2019; originally announced March 2019.

Comments: 16 pages, 11 figures, submitted to MNRAS

Mapping the stability of stellar rotating spheres via linear response theory

Authors: Simon Rozier, Jean-Baptiste Fouvry, Philip G. Breen, Anna Lisa Varri, Christophe Pichon, Douglas C. Heggie

Abstract: Rotation is ubiquitous in the Universe, and recent kinematic surveys have shown that early type galaxies and globular clusters are no exception. Yet the linear response of spheroidal rotating stellar systems has seldom been studied. This paper takes a step in this direction by considering the behaviour of spherically symmetric systems with differential rotation. Specifically, the stability of seve… ▽ More Rotation is ubiquitous in the Universe, and recent kinematic surveys have shown that early type galaxies and globular clusters are no exception. Yet the linear response of spheroidal rotating stellar systems has seldom been studied. This paper takes a step in this direction by considering the behaviour of spherically symmetric systems with differential rotation. Specifically, the stability of several sequences of Plummer spheres is investigated, in which the total angular momentum, as well as the degree and flavour of anisotropy in the velocity space are varied. To that end, the response matrix method is customised to spherical rotating equilibria. The shapes, pattern speeds and growth rates of the systems' unstable modes are computed. Detailed comparisons to appropriate N-body measurements are also presented. The marginal stability boundary is charted in the parameter space of velocity anisotropy and rotation rate. When rotation is introduced, two sequences of growing modes are identified corresponding to radially and tangentially-biased anisotropic spheres respectively. For radially anisotropic spheres, growing modes occur on two intersecting surfaces (in the parameter space of anisotropy and rotation), which correspond to fast and slow modes, depending on the net rotation rate. Generalised, approximate stability criteria are finally presented. △ Less

Submitted 25 February, 2019; originally announced February 2019.

Comments: 18 pages, 20 figures. Submitted to MNRAS

A MODEST review

Authors: Anna Lisa Varri, Maxwell Xu Cai, Francisca Concha-Ramírez, František Dinnbier, Nora Lützgendorf, Václav Pavlík, Sara Rastello, Antonio Sollima, Long Wang, Alice Zocchi

Abstract: We present an account of the state of the art in the fields explored by the research community invested in 'Modeling and Observing DEnse STellar systems'. For this purpose, we take as a basis the activities of the MODEST-17 conference, which was held at Charles University, Prague, in September 2017. Reviewed topics include recent advances in fundamental stellar dynamics, numerical methods for the… ▽ More We present an account of the state of the art in the fields explored by the research community invested in 'Modeling and Observing DEnse STellar systems'. For this purpose, we take as a basis the activities of the MODEST-17 conference, which was held at Charles University, Prague, in September 2017. Reviewed topics include recent advances in fundamental stellar dynamics, numerical methods for the solution of the gravitational N-body problem, formation and evolution of young and old star clusters and galactic nuclei, their elusive stellar populations, planetary systems, and exotic compact objects, with timely attention to black holes of different classes of mass and their role as sources of gravitational waves. Such a breadth of topics reflects the growing role played by collisional stellar dynamics in numerous areas of modern astrophysics. Indeed, in the next decade, many revolutionary instruments will enable the derivation of positions and velocities of individual stars in the Milky Way and its satellites and will detect signals from a range of astrophysical sources in different portions of the electromagnetic and gravitational spectrum, with an unprecedented sensitivity. On the one hand, this wealth of data will allow us to address a number of long-standing open questions in star cluster studies; on the other hand, many unexpected properties of these systems will come to light, stimulating further progress of our understanding of their formation and evolution. △ Less

Submitted 17 October, 2018; originally announced October 2018.

Comments: 42 pages; accepted for publication in 'Computational Astrophysics and Cosmology'. We are much grateful to the organisers of the MODEST-17 conference (Charles University, Prague, September 2017). We acknowledge the input provided by all MODEST-17 participants, and, more generally, by the members of the MODEST community

The complex kinematics of rotating star clusters in a tidal field

Authors: Maria Tiongco, Enrico Vesperini, Anna Lisa Varri

Abstract: 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… ▽ More 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. △ Less

Submitted 18 January, 2018; originally announced January 2018.

Comments: 5 pages, 4 figures, accepted for publication in MNRAS

Three candidate double clusters in the LMC: truth or dare?

Authors: E. Dalessandro, A. Zocchi, A. L. Varri, A. Mucciarelli, M. Bellazzini, F. R. Ferraro, B. Lanzoni, E. Lapenna, L. Origlia

Abstract: The Large Magellanic Cloud (LMC) hosts a large number of candidate stellar cluster pairs. Binary stellar clusters provide important clues about cluster formation processes and the evolutionary history of the host galaxy. However, to properly extract and interpret this information, it is crucial to fully constrain the fraction of real binary systems and their physical properties. Here we present a… ▽ More The Large Magellanic Cloud (LMC) hosts a large number of candidate stellar cluster pairs. Binary stellar clusters provide important clues about cluster formation processes and the evolutionary history of the host galaxy. However, to properly extract and interpret this information, it is crucial to fully constrain the fraction of real binary systems and their physical properties. Here we present a detailed photometric analysis based on ESO-FORS2 images of three candidate cluster multiplets in the LMC, namely SL349-SL353, SL385-SL387-NGC1922 and NGC1836-BRHT4b-NGC1839. For each cluster we derived ages, structural parameters and morphological properties. We have also estimated the degree of filling of their Roche lobe, as an approximate tool to measure the strength of the tidal perturbations induced by the LMC. We find that the members of the possible pairs SL349-SL353 and BRHT4b-NGC1839 have a similar age ($t = 1.00 \pm 0.12$ Gyr and $t = 140 \pm 15$ Myr, respectively), thus possibly hinting to a common origin of their member systems We also find that all candidate pairs in our sample show evidence of intra-cluster overdensities that can be a possible indication of real binarity. Particularly interesting is the case of SL349-SL353. In fact, SL353 is relatively close to the condition of critical filling, thus suggesting that these systems might actually constitute an energetically bound pair. It is therefore key to pursue a detailed kinematic screening of such clusters, without which, at present, we do not dare making a conclusive statement about the true nature of this putative pair. △ Less

Submitted 7 November, 2017; originally announced November 2017.

Comments: 14 pages, 12 figures, 3 tables; Accepted for publication in MNRAS

The Kinematic Richness of Star Clusters I. Isolated Spherical Models with Primordial Anisotropy

Authors: Philip G. Breen, Anna Lisa Varri, Douglas C. Heggie

Abstract: We investigate the dynamical evolution of isolated equal-mass star cluster models by means of direct N-body simulations, primarily focusing on the effects of the presence of primordial anisotropy in the velocity space. We found evidence of the existence of a monotonic relationship between the moment of core collapse and the amount and flavour of anisotropy in the stellar system. Specifically, equi… ▽ More We investigate the dynamical evolution of isolated equal-mass star cluster models by means of direct N-body simulations, primarily focusing on the effects of the presence of primordial anisotropy in the velocity space. We found evidence of the existence of a monotonic relationship between the moment of core collapse and the amount and flavour of anisotropy in the stellar system. Specifically, equilibria characterised by the same initial structural properties (Plummer density profile) and with different degrees of tangentially-biased (radially-biased) anisotropy, reach core collapse earlier (later) than isotropic models. We interpret this result in light of an accelerated (delayed) phase of the early evolution of collisional stellar systems "anisotropic-response"), which we have characterised both in terms of the evolution of the velocity moments and of a fluid model of two-body relaxation. For the case of the most tangentially anisotropic model the initial phase of evolution involves a catastrophic collapse of the inner part of the system which continues until an isotropic velocity distribution is reached. This study represents a first step towards a comprehensive investigation of the role played by kinematic richness in the long-term dynamical evolution of collisional systems. △ Less

Submitted 10 July, 2017; originally announced July 2017.

Comments: 13 pages, 12 figures, accepted for publication in MNRAS

Hubble Space Telescope Proper Motion (HSTPROMO) Catalogs of Galactic Globular Clusters. V. The rapid rotation of 47 Tuc traced and modeled in three dimensions

Authors: A. Bellini, P. Bianchini, A. L. Varri, J. Anderson, G. Piotto, R. P. van der Marel, E. Vesperini, L. L. Watkins

Abstract: 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… ▽ More 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. △ Less

Submitted 27 June, 2017; originally announced June 2017.

Comments: 15 pages, 14 figures, 8 tables, accepted for publication in ApJ

Stellar envelopes of globular clusters embedded in dark mini-haloes

Authors: Jorge Peñarrubia, Anna Lisa Varri, Philip G. Breen, Annette M. N. Ferguson, Rubén Sánchez-Janssen

Abstract: We show that hard encounters in the central regions of globular clusters embedded in dark matter (DM) haloes necessarily lead to the formation of gravitationally-bound stellar envelopes that extend far beyond the nominal tidal radius of the system. Using statistical arguments and numerical techniques we derive the equilibrium distribution function of stars ejected from the centre of a non-divergen… ▽ More We show that hard encounters in the central regions of globular clusters embedded in dark matter (DM) haloes necessarily lead to the formation of gravitationally-bound stellar envelopes that extend far beyond the nominal tidal radius of the system. Using statistical arguments and numerical techniques we derive the equilibrium distribution function of stars ejected from the centre of a non-divergent spherical potential. Independently of the velocity distribution with which stars are ejected, GC envelopes have density profiles that approach asymptotically $ρ\sim r^{-4}$ at large distances and become isothermal towards the centre. Adding a DM halo component leaves two clear-cut observational signatures: (i) a flattening, or slightly increase of the projected velocity dispersion profile at large distances, and (ii) an outer surface density profile that is systematically shallower than in models with no dark matter. △ Less

Submitted 8 June, 2017; originally announced June 2017.

Comments: Accepted to MNRAS letters

Kinematical evolution of tidally limited star clusters: rotational properties

Authors: Maria A. Tiongco, Enrico Vesperini, Anna Lisa Varri

Abstract: 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, ultimate… ▽ More 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. △ Less

Submitted 19 April, 2017; originally announced April 2017.

Comments: 10 pages, 14 figures, accepted for publication in MNRAS

An Approximate Analytic Model of a Star Cluster with Potential Escapers

Authors: Kathryne J. Daniel, Douglas C. Heggie, Anna Lisa Varri

Abstract: In the context of a star cluster moving on a circular galactic orbit, a "potential escaper" is a cluster star that has orbital energy greater than the escape energy, and yet is confined within the Jacobi radius of the stellar system. On the other hand analytic models of stellar clusters typically have a truncation energy equal to the cluster escape energy, and therefore explicitly exclude these en… ▽ More In the context of a star cluster moving on a circular galactic orbit, a "potential escaper" is a cluster star that has orbital energy greater than the escape energy, and yet is confined within the Jacobi radius of the stellar system. On the other hand analytic models of stellar clusters typically have a truncation energy equal to the cluster escape energy, and therefore explicitly exclude these energetically unbound stars. Starting from the landmark analysis performed by Henon of periodic orbits of the circular Hill equations, we present a numerical exploration of the population of "non-escapers", defined here as those stars which remain within two Jacobi radii for several galactic periods, with energy above the escape energy. We show that they can be characterised by the Jacobi integral and two further approximate integrals, which are based on perturbation theory and ideas drawn from Lidov-Kozai theory. Finally we use these results to construct an approximate analytic model that includes a phase space description of a population resembling that of potential escapers, in addition to the usual bound population. △ Less

Submitted 13 March, 2017; originally announced March 2017.

Comments: 23 pages, 19 figures, 1 table

Differences in the rotational properties of multiple stellar populations in M 13: a faster rotation for the "extreme" chemical subpopulation

Authors: M. J. Cordero, V. Hénault-Brunet, C. A. Pilachowski, E. Balbinot, C. I. Johnson, A. L. Varri

Abstract: We use radial velocities from spectra of giants obtained with the WIYN telescope, coupled with existing chemical abundance measurements of Na and O for the same stars, to probe the presence of kinematic differences among the multiple populations of the globular cluster (GC) M13. To characterise the kinematics of various chemical subsamples, we introduce a method using Bayesian inference along with… ▽ More We use radial velocities from spectra of giants obtained with the WIYN telescope, coupled with existing chemical abundance measurements of Na and O for the same stars, to probe the presence of kinematic differences among the multiple populations of the globular cluster (GC) M13. To characterise the kinematics of various chemical subsamples, we introduce a method using Bayesian inference along with an MCMC algorithm to fit a six-parameter kinematic model (including rotation) to these subsamples. We find that the so-called "extreme" population (Na-enhanced and extremely O-depleted) exhibits faster rotation around the centre of the cluster than the other cluster stars, in particular when compared to the dominant "intermediate" population (moderately Na-enhanced and O-depleted). The most likely difference between the rotational amplitude of this extreme population and that of the intermediate population is found to be $\sim$4 km s$^{-1}$, with a 98.4% probability that the rotational amplitude of the extreme population is larger than that of the intermediate population. We argue that the observed difference in rotational amplitudes, obtained when splitting subsamples according to their chemistry, is not a product of the long-term dynamical evolution of the cluster, but more likely a surviving feature imprinted early in the formation history of this GC and its multiple populations. We also find an agreement (within uncertainties) in the inferred position angle of the rotation axis of the different subpopulations considered. We discuss the constraints that these results may place on various formation scenarios. △ Less

Submitted 28 October, 2016; originally announced October 2016.

Comments: Accepted for publication in MNRAS

Kinematical evolution of tidally limited star clusters: the role of retrograde stellar orbits

Authors: Maria Tiongco, Enrico Vesperini, Anna Lisa Varri

Abstract: 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… ▽ More 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. △ Less

Submitted 21 June, 2016; originally announced June 2016.

Comments: 10 pages, 11 figures, accepted for publication in MNRAS

Spatial and kinematic segregation in star cluster merger remnants

Authors: David R. Cole, Victor P. Debattista, Anna Lisa Varri, Markus Hartmann, Anil C. Seth

Abstract: Globular clusters which exhibit chemical and dynamical complexity have been suggested to be the stripped nuclei of dwarf galaxies (e.g., M54, $ω$ Cen). We use $N$-body simulations of nuclear star clusters forming via the mergers of star clusters to explore the persistence of substructure in the phase space. We find that the observed level of differentiation is difficult to reconcile with the obser… ▽ More Globular clusters which exhibit chemical and dynamical complexity have been suggested to be the stripped nuclei of dwarf galaxies (e.g., M54, $ω$ Cen). We use $N$-body simulations of nuclear star clusters forming via the mergers of star clusters to explore the persistence of substructure in the phase space. We find that the observed level of differentiation is difficult to reconcile with the observed if nuclear clusters form wholly out of the mergers of star clusters. Only the star clusters that merged most recently retain sufficiently distinct kinematics to be distinguishable from the rest of the nuclear cluster though the critical factor is the number of merger events not the elapsed time. In situ star formation must therefore be included to explain the observed properties of nuclear star clusters, in good agreement with previous results. △ Less

Submitted 10 May, 2016; originally announced May 2016.

Comments: 16 pages, 17 figures, 3 tables

Testing lowered isothermal models with direct N-body simulations of globular clusters

Authors: Alice Zocchi, Mark Gieles, Vincent Hénault-Brunet, Anna Lisa Varri

Abstract: Several self-consistent models have been proposed, aiming at describing the phase space distribution of stars in globular clusters. This study explores the ability of the recently proposed LIMEPY models (Gieles & Zocchi) to reproduce the dynamical properties of direct N-body models of a cluster in a tidal field, during its entire evolution. These dynamical models include prescriptions for the trun… ▽ More Several self-consistent models have been proposed, aiming at describing the phase space distribution of stars in globular clusters. This study explores the ability of the recently proposed LIMEPY models (Gieles & Zocchi) to reproduce the dynamical properties of direct N-body models of a cluster in a tidal field, during its entire evolution. These dynamical models include prescriptions for the truncation and the degree of radially-biased anisotropy contained in the system, allowing us to explore the interplay between the role of anisotropy and tides in various stages of the life of star clusters. We show that the amount of anisotropy in an initially tidally underfilling cluster increases in the pre-collapse phase, and then decreases with time, due to the effect of the external tidal field on its spatial truncation. This is reflected in the correspondent model parameters, and the best-fit models reproduce the main properties of the cluster at all stages of its evolution, except for the phases immediately preceding and following core collapse. We also notice that the best-fit LIMEPY models are significantly different from isotropic King models, especially in the first part of the evolution of the cluster. Our results put limits on the amount of radial anisotropy that can be expected for clusters evolving in a tidal field, which is important to understand other factors that could give rise to similar observational signatures, such as the presence of an intermediate-mass black hole. △ Less

Submitted 6 May, 2016; originally announced May 2016.

Comments: 20 pages, 11 figures. Accepted for publication in MNRAS

A novel look at energy equipartition in globular clusters

Authors: P. Bianchini, G. van de Ven, M. A. Norris, E. Schinnerer, A. L. Varri

Abstract: 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… ▽ More 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 $σ(m)$ can be described by an exponential function $σ^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. △ Less

Submitted 2 March, 2016; originally announced March 2016.

Comments: 11 pages, 6 figures, 3 tables. Accepted for publication in MNRAS

Velocity anisotropy in tidally limited star clusters

Authors: Maria Tiongco, Enrico Vesperini, Anna Lisa Varri

Abstract: 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.… ▽ More 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. △ Less

Submitted 6 November, 2015; originally announced November 2015.

Comments: 9 pages, 13 figures. Accepted for publication in MNRAS

Evolution of star clusters on eccentric orbits

Authors: Maxwell Xu Cai, Mark Gieles, Douglas C. Heggie, Anna Lisa Varri

Abstract: We study the evolution of star clusters on circular and eccentric orbits using direct $N$-body simulations. We model clusters with initially $N=8{\rm k}$ and $N=16{\rm k}$ single stars of the same mass, orbiting around a point-mass galaxy. For each orbital eccentricity that we consider, we find the apogalactic radius at which the cluster has the same lifetime as the cluster with the same $N$ on a… ▽ More We study the evolution of star clusters on circular and eccentric orbits using direct $N$-body simulations. We model clusters with initially $N=8{\rm k}$ and $N=16{\rm k}$ single stars of the same mass, orbiting around a point-mass galaxy. For each orbital eccentricity that we consider, we find the apogalactic radius at which the cluster has the same lifetime as the cluster with the same $N$ on a circular orbit. We show that then, the evolution of bound particle number and half-mass radius is approximately independent of eccentricity. Secondly, when we scale our results to orbits with the same semi-major axis, we find that the lifetimes are, to first order, independent of eccentricity. When the results of Baumgardt and Makino for a singular isothermal halo are scaled in the same way, the lifetime is again independent of eccentricity to first order, suggesting that this result is independent of the Galactic mass profile. From both sets of simulations we empirically derive the higher order dependence of the lifetime on eccentricity. Our results serve as benchmark for theoretical studies of the escape rate from clusters on eccentric orbits. Finally, our results can be useful for generative models for cold streams and cluster evolution models that are confined to spherical symmetry and/or time-independent tides, such as Fokker-Planck models, Monte Carlo models, and (fast) semi-analytic models. △ Less

Submitted 5 October, 2015; originally announced October 2015.

Comments: 7 pages, 8 figures, accepted for publication in the Monthly Notices of the Royal Astronomical Society

Biases in the determination of dynamical parameters of star clusters: today and in the Gaia era

Authors: A. Sollima, H. Baumgardt, A. Zocchi, E. Balbinot, M. Gieles, V. Henault-Brunet, A. L. Varri

Abstract: The structural and dynamical properties of star clusters are generally derived by means of the comparison between steady-state analytic models and the available observables. With the aim of studying the biases of this approach, we fitted different analytic models to simulated observations obtained from a suite of direct N-body simulations of star clusters in different stages of their evolution and… ▽ More The structural and dynamical properties of star clusters are generally derived by means of the comparison between steady-state analytic models and the available observables. With the aim of studying the biases of this approach, we fitted different analytic models to simulated observations obtained from a suite of direct N-body simulations of star clusters in different stages of their evolution and under different levels of tidal stress to derive mass, mass function and degree of anisotropy. We find that masses can be under/over-estimated up to 50% depending on the degree of relaxation reached by the cluster, the available range of observed masses and distances of radial velocity measures from the cluster center and the strength of the tidal field. The mass function slope appears to be better constrainable and less sensitive to model inadequacies unless strongly dynamically evolved clusters and a non-optimal location of the measured luminosity function are considered. The degree and the characteristics of the anisotropy developed in the N-body simulations are not adequately reproduced by popular analytic models and can be detected only if accurate proper motions are available. We show how to reduce the uncertainties in the mass, mass-function and anisotropy estimation and provide predictions for the improvements expected when Gaia proper motions will be available in the near future. △ Less

Submitted 11 May, 2015; originally announced May 2015.

Comments: 14 pages, 8 figures, accepted for publication by MNRAS

Neutron star natal kicks and the long-term survival of star clusters

Authors: Filippo Contenta, Anna Lisa Varri, Douglas C. Heggie

Abstract: We investigate the dynamical evolution of a star cluster in an external tidal field by using N-body simulations, with focus on the effects of the presence or absence of neutron star natal velocity kicks.We show that, even if neutron stars typically represent less than 2% of the total bound mass of a star cluster, their primordial kinematic properties may affect the lifetime of the system by up to… ▽ More We investigate the dynamical evolution of a star cluster in an external tidal field by using N-body simulations, with focus on the effects of the presence or absence of neutron star natal velocity kicks.We show that, even if neutron stars typically represent less than 2% of the total bound mass of a star cluster, their primordial kinematic properties may affect the lifetime of the system by up to almost a factor of four. We interpret this result in the light of two known modes of star cluster dissolution, dominated by either early stellar evolution mass loss or two-body relaxation. The competition between these effects shapes the mass loss profile of star clusters, which may either dissolve abruptly ("jumping"), in the pre-core-collapse phase, or gradually ("skiing"), after having reached core collapse. △ Less

Submitted 9 February, 2015; originally announced February 2015.

Comments: 5 pages, 5 figures, 1 table. Accepted for publication in MNRAS Letters

The inefficiency of satellite accretion in forming extended star clusters

Authors: Paolo Bianchini, Florent Renaud, Mark Gieles, Anna Lisa Varri

Abstract: 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… ▽ More 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. △ Less

Submitted 4 November, 2014; originally announced November 2014.

Comments: 5 pages, 4 figures, 1 table. Accepted for publication in MNRAS Letters

Kinematical fingerprints of star cluster early dynamical evolution

Authors: Enrico Vesperini, Anna Lisa Varri, Stephen L. W. McMillan, Stephen E. Zepf

Abstract: We study the effects of the external tidal field on the violent relaxation phase of star clusters dynamical evolution, with particular attention to the kinematical properties of the equilibrium configurations emerging at the end of this phase.We show that star clusters undergoing the process of violent relaxation in the tidal field of their host galaxy can acquire significant internal differential… ▽ More We study the effects of the external tidal field on the violent relaxation phase of star clusters dynamical evolution, with particular attention to the kinematical properties of the equilibrium configurations emerging at the end of this phase.We show that star clusters undergoing the process of violent relaxation in the tidal field of their host galaxy can acquire significant internal differential rotation and are characterized by a distinctive radial variation of the velocity anisotropy. These kinematical properties are the result of the symmetry breaking introduced by the external tidal field in the collapse phase and of the action of the Coriolis force on the orbit of the stars. The resulting equilibrium configurations are characterized by differential rotation, with a peak located between one and two half-mass radii. As for the anisotropy, similar to clusters evolving in isolation, the systems explored in this Letter are characterized by an inner isotropic core, followed by a region of increasing radial anisotropy. However for systems evolving in an external tidal field the degree of radial anisotropy reaches a maximum in the cluster intermediate regions and then progressively decreases, with the cluster outermost regions being characterized by isotropy or a mild tangential anisotropy. Young or old but less-relaxed dynamically young star clusters may keep memory of these kinematical fingerprints of their early dynamical evolution. △ Less

Submitted 13 June, 2014; originally announced June 2014.

Comments: 5 pages 4 figures. Accepted for publication in MNRAS Letters

Rotating Globular Clusters

Authors: P. Bianchini, A. L Varri, G. Bertin, A. Zocchi

Abstract: 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-consiste… ▽ More 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. △ Less

Submitted 26 May, 2013; originally announced May 2013.

Comments: 19 pages, 14 figures, accepted for publication in the Astrophysical Journal

On the density profile of the globular cluster M92

Authors: A. Di Cecco, A. Zocchi, A. L. Varri, M. Monelli, G. Bertin, G. Bono, P. B. Stetson, M. Nonino, R. Buonanno, I. Ferraro, G. Iannicola, A. Kunder, A. R. Walker

Abstract: We present new number density and surface brightness profiles for the globular cluster M92 (NGC 6341). These profiles are calculated from optical images collected with the CCD mosaic camera MegaCam at the Canada-France-Hawaii-Telescope and with the Advanced Camera for Surveys on the Hubble Space Telescope. The ground-based data were supplemented with the Sloan Digital Sky Survey photometric catalo… ▽ More We present new number density and surface brightness profiles for the globular cluster M92 (NGC 6341). These profiles are calculated from optical images collected with the CCD mosaic camera MegaCam at the Canada-France-Hawaii-Telescope and with the Advanced Camera for Surveys on the Hubble Space Telescope. The ground-based data were supplemented with the Sloan Digital Sky Survey photometric catalog. Special care was taken to discriminate candidate cluster stars from field stars and to subtract the background contamination from both profiles. By examining the contour levels of the number density, we found that the stellar distribution becomes clumpy at radial distances larger than about 13 arcminutes, and there is no preferred orientation of contours in space. We performed detailed fits of King and Wilson models to the observed profiles. The best-fit models underestimate the number density inside the core radius. Wilson models better represent the observations, in particular in the outermost cluster regions: the good global agreement of these models with the observations suggests that there is no need to introduce an extra-tidal halo to explain the radial distribution of stars at large radial distances. The best-fit models for the number density and the surface brightness profiles are different, even though they are based on the same observations. Additional tests support the evidence that this fact reflects the difference in the radial distribution of the stellar tracers that determine the observed profiles (main sequence stars for the number density, bright evolved stars for the surface brightness). △ Less

Submitted 7 February, 2013; originally announced February 2013.

Comments: 18 pages, 10 figures, Accepted by AJ

The flattening of globular clusters: internal rotation or velocity anisotropy?

Authors: Paolo Bianchini, Anna Lisa Varri, Giuseppe Bertin, Alice Zocchi

Abstract: 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 des… ▽ More 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. △ Less

Submitted 4 February, 2013; originally announced February 2013.

Comments: 4 pages, 3 figures. To appear in the conference proceedings "Reading the book of globular clusters with the lens of stellar evolution", Mem. S. A. It. Eds. P. Ventura, C. Charbonnel, M. Castellani and M. Di Criscienzo

Self-consistent models of quasi-relaxed rotating stellar systems

Authors: A. L. Varri, G. Bertin

Abstract: Two new families of self-consistent axisymmetric truncated equilibrium models for the description of quasi-relaxed rotating stellar systems are presented. The first extends the spherical King models to the case of solid-body rotation. The second is characterized by differential rotation, designed to be rigid in the central regions and to vanish in the outer parts, where the energy truncation becom… ▽ More Two new families of self-consistent axisymmetric truncated equilibrium models for the description of quasi-relaxed rotating stellar systems are presented. The first extends the spherical King models to the case of solid-body rotation. The second is characterized by differential rotation, designed to be rigid in the central regions and to vanish in the outer parts, where the energy truncation becomes effective. The models are constructed by solving the nonlinear Poisson equation for the self-consistent mean-field potential. For rigidly rotating configurations, the solutions are obtained by an asymptotic expansion on the rotation strength parameter. The differentially rotating models are constructed by means of an iterative approach based on a Legendre series expansion of the density and the potential. The two classes of models exhibit complementary properties. The rigidly rotating configurations are flattened toward the equatorial plane, with deviations from spherical symmetry that increase with the distance from the center. For models of the second family, the deviations from spherical symmetry are strongest in the central region, whereas the outer parts tend to be quasi-spherical. The relevant parameter spaces are explored and the intrinsic and projected structural properties are described. Special attention is given to the effect of different options for the truncation of the distribution function in phase space. Models in the moderate rotation regime are best suited to applications to globular clusters. For general interest in stellar dynamics, at high values of the rotation strength the differentially rotating models exhibit a toroidal core embedded in a quasi-spherical configuration. Physically simple analytical models of the kind presented here provide insights into dynamical mechanisms and may be a basis for more realistic investigations with the help of N-body simulations. △ Less

Submitted 9 January, 2012; originally announced January 2012.

Comments: 24 pages, 26 figures. To appear in Astronomy and Astrophysics

A dynamical study of Galactic globular clusters under different relaxation conditions

Authors: A. Zocchi, G. Bertin, A. L. Varri

Abstract: We perform a systematic combined photometric and kinematic analysis of a sample of globular clusters under different relaxation conditions, based on their core relaxation time (as listed in available catalogs), by means of two well-known families of spherical stellar dynamical models. Systems characterized by shorter relaxation time scales are expected to be better described by isotropic King mode… ▽ More We perform a systematic combined photometric and kinematic analysis of a sample of globular clusters under different relaxation conditions, based on their core relaxation time (as listed in available catalogs), by means of two well-known families of spherical stellar dynamical models. Systems characterized by shorter relaxation time scales are expected to be better described by isotropic King models, while less relaxed systems might be interpreted by means of non-truncated, radially-biased anisotropic f^(ν) models, originally designed to represent stellar systems produced by a violent relaxation formation process and applied here for the first time to the study of globular clusters. The comparison between dynamical models and observations is performed by fitting simultaneously surface brightness and velocity dispersion profiles. For each globular cluster, the best-fit model in each family is identified, along with a full error analysis on the relevant parameters. Detailed structural properties and mass-to-light ratios are also explicitly derived. We find that King models usually offer a good representation of the observed photometric profiles, but often lead to less satisfactory fits to the kinematic profiles, independently of the relaxation condition of the systems. For some less relaxed clusters, f^(ν) models provide a good description of both observed profiles. Some derived structural characteristics, such as the total mass or the half-mass radius, turn out to be significantly model-dependent. The analysis confirms that, to answer some important dynamical questions that bear on the formation and evolution of globular clusters, it would be highly desirable to acquire larger numbers of accurate kinematic data-points, well distributed over the cluster field. △ Less

Submitted 6 January, 2012; originally announced January 2012.

Comments: 18 pages, 7 figures. Accepted for publication in Astronomy & Astrophysics

Properties of quasi-relaxed stellar systems in an external tidal field

Authors: A. L. Varri, G. Bertin

Abstract: In a previous paper, we have constructed a family of self-consistent triaxial models of quasi-relaxed stellar systems, shaped by the tidal field of the hosting galaxy, as an extension of the well-known spherical King models. For a given tidal field, the models are characterized by two physical scales (such as total mass and central velocity dispersion) and two dimensionless parameters (the conce… ▽ More In a previous paper, we have constructed a family of self-consistent triaxial models of quasi-relaxed stellar systems, shaped by the tidal field of the hosting galaxy, as an extension of the well-known spherical King models. For a given tidal field, the models are characterized by two physical scales (such as total mass and central velocity dispersion) and two dimensionless parameters (the concentration parameter and the tidal strength). The most significant departure from spherical symmetry occurs when the truncation radius of the corresponding spherical King model is of the order of the tidal radius, which, for a given tidal strength, is set by the maximum concentration value admitted. For such maximally extended (or "critical") models the outer boundary has a generally triaxial shape, given by the zero-velocity surface of the relevant Jacobi integral, which is basically independent of the concentration parameter. In turn, the external tidal field can give rise to significant global departures from spherical symmetry (as measured, for example, by the quadrupole of the mass distribution of the stellar system) only for low-concentration models, for which the allowed maximal value of the tidal strength can be relatively high. In this paper we describe in detail the intrinsic and the projected structure and kinematics of the models, covering the entire parameter space, from the case of sub-critical (characterized by "underfilling" of the relevant Roche volume) to that of critical models. The intrinsic properties can be a useful starting point for numerical simulations and other investigations that require initialization of a stellar system in dynamical equilibrium. The projected properties are a key step in the direction of a comparison with observed globular clusters and other candidate stellar systems. △ Less

Submitted 17 August, 2009; originally announced August 2009.

Comments: 13 pages, 11 figures, uses emulateapj.cls with apjfonts.sty. Accepted for publication in The Astrophysical Journal

Journal ref: Astrophys. J. 703: 1911-1922, 2009

The construction of non-spherical models of quasi-relaxed stellar systems

Authors: G. Bertin, A. L. Varri

Abstract: Spherical models of collisionless but quasi-relaxed stellar systems have long been studied as a natural framework for the description of globular clusters. Here we consider the construction of self-consistent models under the same physical conditions, but including explicitly the ingredients that lead to departures from spherical symmetry. In particular, we focus on the effects of the tidal fiel… ▽ More Spherical models of collisionless but quasi-relaxed stellar systems have long been studied as a natural framework for the description of globular clusters. Here we consider the construction of self-consistent models under the same physical conditions, but including explicitly the ingredients that lead to departures from spherical symmetry. In particular, we focus on the effects of the tidal field associated with the hosting galaxy. We then take a stellar system on a circular orbit inside a galaxy represented as a "frozen" external field. The equilibrium distribution function is obtained from the one describing the spherical case by replacing the energy integral with the relevant Jacobi integral in the presence of the external tidal field. Then the construction of the model requires the investigation of a singular perturbation problem for an elliptic partial differential equation with a free boundary, for which we provide a method of solution to any desired order, with explicit solutions to two orders. We outline the relevant parameter space, thus opening the way to a systematic study of the properties of a two-parameter family of physically justified non-spherical models of quasi-relaxed stellar systems. The general method developed here can also be used to construct models for which the non-spherical shape is due to internal rotation. Eventually, the models will be a useful tool to investigate whether the shapes of globular clusters are primarily determined by internal rotation, by external tides, or by pressure anisotropy. △ Less

Submitted 18 August, 2008; originally announced August 2008.

Comments: AASTeX v5.2, 37 pages with 2 figures, accepted for publication in The Astrophysical Journal

Journal ref: Astrophys.J.689:1005-1019,2008