We present a study of the internal kinematics of two globular clusters, M10 (NGC 6254) and M71 (NGC 6838), using individual radial velocity (RV) measurements obtained from observations using the Hydra multi-object spectrograph on the WIYN 3.5 m telescope. We measured 120 RVs for stars in M10, of which 107 were determined to be cluster members. In M71, we measured 82 RVs and determined 78 of those measurements belonged to cluster members. Using the cluster members, we determine a mean RV of $75.9 \pm 4.0$ (s.d.) km s$^{-1}$ and $-22.9 \pm 2.2$ (s.d.) km s$^{-1}$ for M10 and M71, respectively. We combined the Hydra RV measurements with literature samples and performed a line-of-sight rotational analysis on both clusters. Our analysis has not revealed a statistically significant rotation in either of these clusters with the exception of the inner region (10\arcsec - 117\arcsec) of M10 for which we find hints of a marginally significant rotation with amplitude V$_{rot} = 1.14 \pm 0.18$ km s$^{-1}$. For M10, we calculate a central velocity dispersion of $\sigma_{0} = 5.44 \pm 0.61$ km s$^{-1}$, which gives a ratio of the amplitude of rotation to the central velocity dispersion $V_{rot}/\sigma_{0} = 0.21 \pm 0.04$. We also explored the rotation of the multiple stellar populations identified in M10 and M71 and found rotation (or lack thereof) in each population consistent with each other and the cluster global rotation signatures.

A number of scenarios for the formation of multiple populations in globular clusters (GCs) predict that second generation (2G) stars form in a compact and dense subsystem embedded in a more extended first-generation (1G) system. If these scenarios are accurate, a consequence of the denser 2G formation environment is that 2G binaries should be more significantly affected by stellar interactions and disrupted at a larger rate than 1G binaries. The fractions and properties of binary stars can thus provide a dynamical fingerprint of the formation epoch of multiple-population GCs and their subsequent dynamical evolution. We investigate the connection between binaries and multiple populations in five GCs, NGC 288, NGC 6121 (M 4), NGC 6352, NGC 6362, and NGC 6838 (M 71). To do this, we introduce a new method based on the comparison of Hubble Space Telescope observations of binaries in the F275W, F336W, F438W, F606W and F814W filters with a large number of simulated binaries. In the inner regions probed by our data we do not find large differences between the local 1G and the 2G binary incidences in four of the studied clusters, the only exception being M 4 where the 1G binary incidence is about three times larger than the 2G incidence. The results found are in general agreement with the results of simulations predicting significant differences in the global 1G and 2G incidences and in the local values in the clusters' outer regions but similar incidences in the inner regions. The significant difference found in M 4 is consistent with simulations with a larger fraction of wider binaries. Our analysis also provides the first evidence of mixed (1G-2G) binaries, a population predicted by numerical simulations to form in a cluster's inner regions as a result of stellar encounters during which one component of a binary is replaced by a star of a different population.

We have observed the CN features at $\sim$3800 ~\AA~ and 4120~ \AA~ as well as the CH band at $\sim$4300~\AA~ for 145 evolved stars in the Galactic globular cluster M71 using the multi-object spectrograph, Hydra, on the WIYN-3.5m telescope. We use these measurements to create two $\delta$CN indices finding that both distributions are best fit by two populations, a CN-enhanced and CN-normal. We find that 42 $\pm$ 4\% of the RGB stars in our sample are CN-enhanced. The percentage of CN-enhanced is 40 $\pm$ 13\% for the AGB and 33 $\pm$ 9\% for the HB stars, which suggests there are no missing second generation stars at these stages of stellar evolution. The two generations also separate in magnitude and color on the HB, which allows us to find the difference in He abundance between the two populations by fitting appropriate ZAHBs. The broad range of distances from the cluster's center covered by our sample allows us to study the dependence of the ratio of the number of first to second population stars on the distance from the cluster's center, and we find that this ratio does not vary radially and that the two populations are spatially mixed. Finally, we compare our identification of multiple populations with the classification based on the Na-O anti-correlation and the HST UV photometry, and we find good agreement with both methods.

In the present work we analyzed seven globular clusters selected from their location in the Galactic bulge and with metallicity values in the range $-1.30\lesssim\rm{[Fe/H]}\lesssim-0.50$. The aim of this work is first to derive cluster ages assuming single stellar populations, and secondly, to identify the stars from first (1G) and second generations (2G) from the main sequence, subgiant and red giant branches, and to derive their age differences. Based on a combination of UV and optical filters used in this project, we apply the Gaussian mixture models to distinguish the multiple stellar populations. Applying statistical isochrone fitting, we derive self-consistent ages, distances, metallicities, and reddening values for the sample clusters. An average of $12.3\pm0.4$ Gyr was obtained both using Dartmouth and BaSTI (accounting atomic diffusion effects) isochrones, without a clear distinction between the moderately metal-poor and the more metal-rich bulge clusters, except for NGC 6717 and the inner halo NGC 6362 with $\sim 13.5$ Gyr. We derived a weighted mean age difference between the multiple populations hosted by each globular cluster of $41\pm170$ Myr adopting canonical He abundances; whereas for higher He in 2G stars, this difference reduces to $17\pm170$ Myr, but with individual uncertainties of $500$ Myr.

We studied the spatial distributions of multiple stellar populations (MPs) in a sample of 20 globular clusters (GCs) spanning a broad range of dynamical ages. The differences between first-population (FP) and second-population (SP) stars were measured by means of the parameter $A^+$, defined as the area enclosed between their cumulative radial distributions. We provide the first purely observational evidence of the dynamical path followed by MPs from initial conditions toward a complete FP-SP spatial mixing. Less dynamically evolved clusters have SP stars more centrally concentrated than FPs, while in more dynamically evolved systems the spatial differences between FP and SP stars decrease and eventually disappear. By means of an appropriate comparison with a set of numerical simulations, we show that these observational results are consistent with the evolutionary sequence expected by the long-term dynamical evolution of clusters forming with an initially more centrally concentrated SP sub-system. This result is further supported by the evidence of a trend between $A^+$ and the stage of GC dynamical evolution inferred by the ratio between the present-day and the initial mass of the cluster.

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.

By means of grid-based, 3D hydrodynamical simulations we study the formation of second generation (SG) stars in a young globular cluster (GC) of mass 10^7 Msun, the possible progenitor of an old GC with a present mass ~(1-5) * 10^6 Msun. The cluster accretes external gas as its first generation (FG) asymptotic giant branch (AGB) stars release their ejecta and SG stars form. We consider two models characterised by different densities of the external gas. In both cases, we find that a very compact SG subsystem with central density >10^5M sun/pc^3 forms in the innermost regions of the cluster. The low-density model forms a population of extreme SG stars with high helium enhancement, followed by the formation of another SG group out of a mix of pristine gas and AGB ejecta and characterised by a modest helium enhancement. On the other hand, the high-density model forms in prevalence SG stars with modest helium enhancement. Our simulations illustrate the dynamical processes governing the formation of SG populations in GCs and shed light on the structural properties emerging at the end of this phase. The newly born SG groups have different concentrations, with more extreme SG stars more centrally concentrated than those with less extreme chemical abundances. The very high density of the SG subsystems implies that SG massive stars, if formed, might suffer frequent close encounters, collisions and gas stripping, thus possibly contributing further gas to the SG formation.

In this paper we report on the discovery of a double blue straggler star (BSS) sequence in the core of the core-collapsed cluster M15 (NGC 7078). We performed a detailed photometric analysis of the extremely dense core of the cluster using a set of images secured with the Advanced Camera for Survey in the High Resolution Channel mode on-board the Hubble Space Telescope. The proper combination of the large number of single frames in the near-UV (F220W), and blue (F435W) filters allowed us to perform a superb modeling of the Point Spread Function and an accurate deblending procedure. The Color-Magnitude diagram revealed the presence of two distinct parallel sequences of blue stragglers. In particular, the blue BSS sequence is characterized by the intriguing presence of two different branches. The first branch appears extremely narrow, it extends up to 2.5 magnitudes brighter than the cluster main-sequence turnoff (MS-TO) point, and it is nicely reproduced by a 2 Gyr-old collisional isochrone. The second branch extends up to 1.5 magnitudes from the MS-TO and it is reproduced by a 5.5 Gyr-old collisional isochrone. Our observations suggest that each of these branches is mainly constituted by a population of nearly coeval collisional BSS of different masses generated during two episodes of high collisional activity. We discuss the possibility that the oldest episode corresponds to the core-collapse event (occurred about 5.5 Gyr ago), while the most recent one (occurred about 2 Gyr ago) is associated with a core oscillation in the post-core collapse evolution. The discovery of these features provides further strong evidence in support of the connection between the BSS properties and GC dynamical evolution, and it opens new perspectives on the study of core-collapse and post core-collapse evolution.

The outer regions of globular clusters can enable us to answer many fundamental questions concerning issues ranging from the formation and evolution of clusters and their multiple stellar populations to the study of stars near and beyond the hydrogen-burning limit and to the dynamics of the Milky Way. The outskirts of globular clusters are still uncharted territories observationally. A very efficient way to explore them is through high-precision proper motions of low-mass stars over a large field of view. The Wide Field InfraRed Survey Telescope (WFIRST) combines all these characteristics in a single telescope, making it the best observational tool to uncover the wealth of information contained in the clusters' outermost regions.

We used high-resolution optical HST/WFC3 and multi-conjugate adaptive optics assisted GEMINI GeMS/GSAOI observations in the near-infrared to investigate the physical properties of the globular cluster NGC 6569 in the Galactic bulge. We have obtained the deepest purely NIR color-magnitude diagram published so far for this cluster using ground-based observations, reaching $K_{s}$ $\approx$ 21.0 mag (two magnitudes below the main-sequence turn-off point). By combining the two datasets secured at two different epochs, we determined relative proper motions for a large sample of individual stars in the center of NGC 6569, allowing a robust selection of cluster member stars. Our proper motion analysis solidly demonstrates that, despite its relatively high metal content, NGC 6569 hosts some blue horizontal branch stars. A differential reddening map has been derived in the direction of the system, revealing a maximum color excess variation of about $\delta E(B-V)$ $\sim$ 0.12 mag in the available field of view. The absolute age of NGC 6569 has been determined for the first time. In agreement with the other few bulge globular clusters with available age estimates, NGC 6569 turns out to be old, with an age of about 12.8 Gyr, and a typical uncertainty of 0.8-1.0 Gyr.

We present the analysis of the radial distributions and kinematic properties of the multiple stellar populations (mPOPs) hosted in the globular cluster (GC) NGC 6352 as part of the Hubble Space Telescope (HST) "UV Legacy Survey of Galactic Globular Clusters" program. NGC 6352 is one of the few GCs for which the mPOP tagging in appropriate color-magnitude diagrams is clear in all evolutionary sequences. We computed high-precision stellar proper motions for the stars from the cluster's core out to 75 arcsec ($\sim$1.5 core radii, or $\sim$0.6 half-light radii). We find that, in the region explored, first- and second-generation stars share the same radial distribution and kinematic properties. Velocity dispersions, anisotropy radial profiles, differential rotation, and level of energy equipartition, all suggest that NGC 6352 is probably in an advanced evolutionary stage, and any possible difference in the structural and kinematic properties of its mPOPs have been erased by dynamical processes in the core of the cluster. We also provide an estimate of the mass of blue stragglers and of main-sequence binaries through kinematics alone. In general, in order to build a complete dynamical picture of this and other GCs, it will be essential to extend the analyses presented in this paper to the GCs' outer regions where some memories of the initial differences in the mPOP properties, and those imprinted by dynamical processes, might still be present.

In the globular cluster NGC2808, a quasi-standard initial lithium abundance is derived for a red giant belonging to the `extreme' population, characterized by a large helium overabundance, and by abundances of proton capture elements typical of nuclear processing in gas at very high temperatures, where the initial lithium has been fully destroyed. The observations of lithium in such extreme cluster stars are important to test different models for the formation of multiple populations in old Globular Clusters. In the asymptotic giant branch (AGB) scenario, fresh lithium is synthetized during the initial phases of hot bottom burning which, afterwards, synthetize the other p-capture elements. We model the abundance of lithium in the ejecta of superAGB models, finding values consistent or larger than observed in the `extreme' giant; these same models describe correctly the magnesium depletion and silicon enrichment of the extreme population of NGC 2808, so the overall agreement provides further support to the AGB scenario. In the models involving massive or supermassive stars, the Lithium observed requires a mixture of the lithium-free ejecta of the polluting population with more than 40% of standard-lithium pristine gas. The extended chemical anomalies of NGC 2808 stars are then to be all explained within at most 60% of the possible dilution range, the initial helium mass fraction in the ejecta should be Y >= 0.5, to account for the Ye 0.38-0.40 of the extreme population, and further observations of p-process elements are needed to check the model.

We present the results of a study aimed at investigating the effects of dynamical evolution on the spatial distribution and mixing of primordial binary stars in multiple-population globular clusters. Multiple stellar population formation models predict that second-generation (SG) stars form segregated in the inner regions of a more extended first-generation (FG) cluster. Our study, based on the results of a survey of N-body simulations, shows that the spatial mixing process for binary stars is more complex than that of single stars since additional processes such as binary ionization, recoil and ejection following binary-single and binary-binary interactions play a key role in determining the spatial distribution of the population of surviving binaries. The efficiency and relative importance of these additional effects depends on the binary binding energy and determines the timescale of the spatial mixing of FG and SG binaries. Our simulations illustrate the role of ionization, recoil and ejection combined with the effects of mass segregation driven by two-body relaxation and show that the complex interplay of all these processes results in a significant extension of the time needed for the complete spatial mixing of FG and SG binaries compared to that of single stars. Clusters in which FG and SG single stars have already reached complete spatial mixing might be characterized by a significant radial gradient in the ratio of the FG-to-SG binary fraction. The implications of the delayed mixing of FG and SG binaries for the differences between the kinematics of the two populations are discussed.

Recent work, based on data from the Hubble Space Telescope (HST) UV Legacy Survey of Galactic Globular Clusters (GCs), has revealed that all the analyzed clusters host two groups of first- (1G) and second-generation (2G) stars. In most GCs, both 1G and 2G stars host sub-stellar populations with different chemical composition. We compare multi-wavelength HST photometry with synthetic spectra to determine for the first time the average helium difference between the 2G and 1G stars in a large sample of 57 GCs and the maximum helium variation within each of them. We find that in all clusters 2G stars are consistent with being enhanced in helium with respect to 1G. The maximum helium variation ranges from less than 0.01 to more than 0.10 in helium mass fraction and correlates with both the cluster mass and the color extension of the horizontal branch (HB). These findings demonstrate that the internal helium variation is one of the main (second) parameters governing the HB morphology.

The dense environment of globular clusters (GCs) can facilitate the formation of binary black holes (BBHs), some of which can merge with gravitational waves (GW) within the age of the Universe. We have performed a survey of Monte-Carlo simulations following the dynamical evolution of GCs with different masses, sizes and binary fractions and explored the impact of the host GC properties on the formation of BBH mergers. We find that the number of BBH mergers from GCs is determined by the GC's initial mass, size and primordial binary fraction. We identify two groups of BBH mergers: a primordial group whose formation does not depend on cluster's dynamics and a dynamical group whose formation is driven by the cluster's dynamical evolution. We show how the BBH origin affects the BBH mergers' main properties such as the chirp mass and merging time distributions. We provide analytic expressions for the dependence of the number of BBH mergers from individual GCs on the main cluster's structural properties and the time evolution of the merger rates of these BBHs. These expressions provide an essential ingredient for a general framework allowing to estimate the merger rate density. Using the relations found in our study, we find a local merger rate density of 0.18$-$1.8 ${\rm Gpc}^{-3}{\rm yr}^{-1}$ for primordial BBH mergers and 0.6$-$18 ${\rm Gpc}^{-3}{\rm yr}^{-1}$ for dynamical BBH mergers, depending on the GC mass and size distributions, initial binary fraction and the number density of GCs in the Universe.

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

We present a detailed analysis of the kinematic properties of the multiple populations (MPs) in the low-mass Galactic globular cluster NGC6362 based on a sample of about 500 member stars for which radial velocities (RVs), Fe and Na abundances have been homogeneously derived. At distances from the cluster center larger than about 0.5 $r_h$, we find that first (FG - Na-poor) and second generation (SG - Na-rich) stars show hints of different line-of-sight velocity dispersion profiles, with FG stars being dynamically hotter. This is the first time that differences in the velocity dispersion of MPs are detected by using only RVs. While kinematic differences between MPs in globular clusters are usually described in terms of anisotropy differences driven by the different radial distributions, this explanation seems hardly viable for NGC6362, where SG and FG stars are spatially mixed. We demonstrate that the observed difference in the velocity dispersion profiles can be accounted for by the effect of binary stars. In fact, thanks to our multi-epoch RV measurements we find that the binary fraction is significantly larger in the FG sample (f~14%) than in SG population (f<1%), and we show that such a difference can inflate the velocity dispersion of FG with respect to SG by the observed amount in the relevant radial range. Our results nicely match the predictions of state-of-the art $N$-body simulations of the co-evolution of MPs in globular clusters that include the effect of binaries.

Through a suite of direct N-body simulations, we explore how the structural and kinematic evolution of a star cluster located at the center of a dwarf galaxy is affected by the shape of its host's dark matter density profile. The stronger central tidal fields of cuspier halos minimize the cluster's ability to expand in response to mass loss due to stellar evolution during its early evolutionary stages and during its subsequent long-term evolution driven by two-body relaxation. Hence clusters evolving in cuspier dark matter halos are characterized by more compact sizes, higher velocity dispersions and remain approximately isotropic at all clustercentric distances. Conversely, clusters in cored halos can expand more and develop a velocity distribution profile that becomes increasingly radially anisotropic at larger clustercentric distances. Finally, the larger velocity dispersion of clusters evolving in cuspier dark matter profiles results in them having longer relaxation times. Hence clusters in cuspy galaxies relax at a slower rate and, consequently, they are both less mass segregated and farther from complete energy equipartition than cluster's in cored galaxies. Application of this work to observations allows for star clusters to be used as tools to measure the distribution of dark matter in dwarf galaxies and to distinguish isolated star clusters from ultra-faint dwarf galaxies.

We present an improved data-reduction technique to obtain high-precision proper motions (PMs) of globular clusters using Hubble Space Telescope data. The new reduction is superior to the one presented in the first paper of this series for the faintest sources in very crowded fields. We choose the globular cluster NGC 362 as a benchmark to test our new procedures. We measure PMs of 117 450 sources in the field, showing that we are able to obtain a PM precision better than 10 $\mu$as yr$^{-1}$ for bright stars. We make use of this new PM catalog of NGC 362 to study the cluster's internal kinematics. We investigate the velocity-dispersion profiles of the multiple stellar populations hosted by NGC 362 and find new pieces of information on the kinematics of first- and second-generation stars. We analyze the level of energy equipartition of the cluster and find direct evidence for its post-core-collapsed state from kinematic arguments alone. We refine the dynamical mass of the blue stragglers and study possible kinematic differences between blue stragglers formed by collisions and mass transfer. We also measure no significant cluster rotation in the plane of the sky. Finally, we measure the absolute PM of NGC 362 and of the background stars belonging to the Small Magellanic Cloud, finding a good agreement with previous estimates in the literature. We make the PM catalog publicly available.

We present CN and CH band measurements for 137 RGB and AGB stars in the Galactic globular cluster M10. Our measurements come from low resolution spectroscopy taken with the Hydra spectrograph on the WIYN-3.5m telescope. We use these measurements to identify two populations of stars within the cluster, a CN-normal and CN-enhanced, and find that in our sample 60% of stars are CN-enhanced. Our large sample allows us to conduct a detailed analysis on the carbon and nitrogen abundances and the radial distribution of each population separately. Our analysis of the radial dependence shows that each population has the same radial distribution in the cluster, which is likely due to the cluster being dynamically evolved. We also compare our results to other methods of classifying multiple populations in globular clusters such as the Na-O anti-correlation and the HST pseudo-color magnitude diagrams. We find that these three methods of identifying multiple populations are in good agreement with each other for M10 and all lead to an estimate of the fraction of second generation stars approximately equal to 60%. Among AGB stars, when classified by the CN band, there appears to be a lack of second generation stars when compared to the RGB stars. However, when classified by [N/Fe], we find a similar 60% of AGB stars in the second generation. Finally, we use the measured carbon and nitrogen abundances in RGB stars to study the change of each element with magnitude as stars evolve up the RGB, comparing the results to globular clusters of similar metallicity, M3 and M13.

The parameter A+, defined as the area enclosed between the cumulative radial distribution of blue straggler stars (BSSs) and that of a reference population, is a powerful indicator of the level of BSS central segregation. As part of the Hubble Space Telescope UV Legacy Survey of Galactic globular clusters (GCs), here we present the BSS population and the determination of A+ in 27 GCs observed out to about one half-mass radius. In combination with 21 additional clusters discussed in a previous paper this provides us with a global sample of 48 systems (corresponding to \sim 32\% of the Milky Way GC population), for which we find a strong correlation between A+ and the ratio of cluster age to the current central relaxation time. Tight relations have been found also with the core radius and the central luminosity density, which are expected to change with the long-term cluster dynamical evolution. An interesting relation is emerging between A+ and the ratio of the BSS velocity dispersion relative to that of main sequence turn-off stars, which measures the degree of energy equipartition experienced by BSSs in the cluster. These results provide further confirmation that BSSs are invaluable probes of GC internal dynamics and A+ is a powerful dynamical clock.

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

We present the first results of the Multi-Instrument Kinematic Survey of Galactic Globular Clusters, a project aimed at exploring the internal kinematics of a representative sample of Galactic globular clusters from the radial velocity of individual stars, covering the entire radial extension of each system. This is achieved by exploiting the formidable combination of multi-object and integral field unit spectroscopic facilities of the ESO Very Large Telescope. As a first step, here we discuss the results obtained for 11 clusters from high and medium resolution spectra acquired through a combination of FLAMES and KMOS observations. We provide the first kinematical characterization of NGC 1261 and NGC 6496. In all the surveyed systems, the velocity dispersion profile declines at increasing radii, in agreement with the expectation from the King model that best fits the density/luminosity profile. In the majority of the surveyed systems we find evidence of rotation within a few half-mass radii from the center. These results are in general overall agreement with the predictions of recent theoretical studies, suggesting that the detected signals could be the relic of significant internal rotation set at the epoch of the cluster's formation.

We present a detailed analysis of the radial distribution of light-element multiple populations (LE-MPs) in the massive and dense globular cluster M80 based on the combination of UV and optical Hubble Space Telescope data. Surprisingly, we find that first generation stars (FG) are significantly more centrally concentrated than extreme second generation ones (SG) out to $\sim 2.5 r_h$ from the cluster center. To understand the origin of such a peculiar behavior, we used a set of $N$-body simulations following the long-term dynamical evolution of LE-MPs. We find that, given the advanced dynamical state of the cluster, the observed difference does not depend on the primordial relative distributions of FG and SG stars. On the contrary, a difference of $\sim 0.05-0.10 M_{\odot}$ between the average masses of the two sub-populations is needed to account for the observed radial distributions. We argue that such a mass difference might be the result of the higher He abundance of SG stars (of the order of $\Delta Y\sim 0.05-0.06$) with respect to FG. Interestingly, we find that a similar He variation is necessary to reproduce the horizontal branch morphology of M80. These results demonstrate that differences in mass among LE-MPs, due to different He content, should be properly taken into account for a correct interpretation of their radial distribution, at least in dynamically evolved systems.

We present the results of a survey of N-body simulations aimed at studying the effects of the long-term dynamical evolution on the stellar mass function (MF) of multiple stellar populations in globular clusters. Our simulations show that if first-(1G) and second-generation (2G) stars have the same initial MF (IMF), the global MFs of the two populations are affected similarly by dynamical evolution and no significant differences between the 1G and the 2G MFs arise during the cluster's evolution. If the two populations have different IMFs, dynamical effects do not completely erase memory of the initial differences. Should observations find differences between the global 1G and 2G MF, these would reveal the fingerprints of differences in their IMFs. Irrespective of whether the 1G and 2G populations have the same global IMF or not, dynamical effects can produce differences between the local (measured at various distances from the cluster centre) 1G and 2G MFs; these differences are a manifestation of the process of mass segregation in populations with different initial structural properties. In dynamically old and spatially mixed clusters, however, differences between the local 1G and 2G MFs can reveal differences between the 1G and 2G global MFs. In general, for clusters with any dynamical age, large differences between the local 1G and 2G MFs are more likely to be associated with differences in the global MF. Our study also reveals a dependence of the spatial mixing rate on the stellar mass, another dynamical consequence of the multiscale nature of multiple-population clusters.

The split main sequences (MSs) and extended MS turnoffs (eMSTOs) detected in a few young clusters have demonstrated that these stellar systems host multiple populations differing in a number of properties such as rotation and, possibly, age.We analyze Hubble Space Telescope photometry for thirteen clusters with ages between ~40 and ~1000 Myrs and of different masses. Our goal is to investigate for the first time the occurrence of multiple populations in a large sample of young clusters. We find that all the clusters exhibit the eMSTO phenomenon and that MS stars more massive than ~1.6 solar masses define a blue and red MS, with the latter hosting the majority of MS stars. The comparison between the observations and isochrones suggests that the blue MSs are made of slow-rotating stars, while the red MSs host stars with rotational velocities close to the breakup value. About half of the bright MS stars in the youngest clusters are H-alpha emitters. These Be stars populate the red MS and the reddest part of the eMSTO thus supporting the idea that the red MS is made of fast rotators. We conclude that the split MS and the eMSTO are a common feature of young clusters in both Magellanic Clouds. The phenomena of a split MS and an eMSTO occur for stars that are more massive than a specific threshold which is independent of the host-cluster mass. As a by-product, we report the serendipitous discovery of a young SMC cluster, GALFOR1.

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

As part of the Hubble Space Telescope UV Legacy Survey of Galactic Globular Clusters, 110 parallel fields were observed with the Wide Field Channel of the Advanced Camera for Surveys, in the outskirts of 48 globular clusters, plus the open cluster NGC 6791. Totalling about $0.3$ square degrees of observed sky, this is the largest homogeneous Hubble Space Telescope photometric survey of Galalctic globular clusters outskirts to date. In particular, two distinct pointings have been obtained for each target on average, all centred at about $6.5$ arcmin from the cluster centre, thus covering a mean area of about $23\,{\rm arcmin^{2}}$ for each globular cluster. For each field, at least one exposure in both F475W and F814W filters was collected. In this work, we publicly release the astrometric and photometric catalogues and the astrometrised atlases for each of these fields.

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.

The Hubble Space Telescope UV Legacy Survey of Galactic Globular Clusters is providing a major breakthrough in our knowledge of Globular Clusters (GCs) and their stellar populations. Among the main results, we discovered that all the studied GCs host two main discrete groups consisting of first generation (1G) and second generation (2G) stars. We exploit the multiwavelength photometry from this project to investigate, for the first time, the Red Giant Branch Bump (RGBB) of the two generations in a large sample of GCs. We identified, with high statistical significance, the RGBB of 1G and 2G stars in 26 GCs and found that their magnitude separation as a function of the filter wavelength follows comparable trends. The comparison of observations to synthetic spectra reveals that the RGBB luminosity depends on the stellar chemical composition and that the 2G RGBB is consistent with stars enhanced in He and N and depleted in C and O with respect to 1G stars. For metal-poor GCs the 1G and 2G RGBB relative luminosity in optical bands mostly depends on helium content, Y. We used the RGBB observations in F606W and F814W bands to infer the relative helium abundance of 1G and 2G stars in 18 GCs, finding an average helium enhancement $\Delta$Y=0.011$\pm$0.002 of 2G stars with respect to 1G stars. This is the first determination of the average difference in helium abundance of multiple populations in a large number of clusters and provides a lower limit to the maximum internal variation of helium in GCs.

The color-magnitude diagrams of many Magellanic Cloud clusters (with ages up to 2 billion years) display extended turnoff regions where the stars leave the main sequence, suggesting the presence of multiple stellar populations with ages which may differ even by hundreds million years (Mackey et al. 2008, Milone et al. 2009, Girardi et al. 2011). A strongly debated question is whether such an extended turnoff is instead due to populations with different stellar rotations (Girardi et al. 2011, Goudfrooij et al. 2011, Rubele et al. 2013, Li et al. 2014). The recent discovery of a `split' main sequence in some younger clusters (about 80--400Myr) added another piece to this puzzle. The blue (red) side of the main sequence is consistent with slowly (rapidly) rotating stellar models (D'Antona et al. 2015, Milone et al. 2016, Correnti et al. 2017, Milone et al 2016), but a complete theoretical characterization of the observed color-magnitude diagram appeared to require also an age spread (Correnti et al. 2017). We show here that, in three clusters so far analyzed, if the blue main sequence stars are interpreted with models that have been always slowly rotating, they must be about 30% younger than the rest of the cluster. If they are instead interpreted as stars initially rapidly rotating, but that have later slowed down, the age difference disappears, and "braking" also helps to explain the apparent age differences of the extended turnoff. The age spreads in Magellanic Cloud clusters are a manifestation of rotational stellar evolution. Observational tests are suggested.

We measure how the slope $\alpha$ of the stellar mass function (MF) changes as a function of clustercentric distance $r$ in five Galactic globular clusters and compare $\alpha(r)$ to predictions from direct $N$-body star cluster simulations. Theoretical studies predict that $\alpha(r)$ (which traces the degree of mass segregation in a cluster) should steepen with time as a cluster undergoes two-body relaxation and that the amount by which the global MF can evolve from its initial state due to stellar escape is directly linked to $\alpha(r)$. We find that the amount of mass segregation in M10, NGC 6218, and NGC 6981 is consistent with their dynamical ages, but only the global MF of M10 is consistent with its degree of mass segregation as well. NGC 5466 and NGC 6101 on the other hand appear to be less segregated than their dynamical ages would indicate. Furthermore, despite the fact that the escape rate of stars in non-segregated clusters is independent of stellar mass, both NGC 5466 and NGC 6101 have near-flat MFs. We discuss various mechanisms which could produce non-segregated clusters with near-flat MFs, including higher mass-loss rates and black hole retention, but argue that for some clusters (NGC 5466 and NGC 6101) explaining the present-day properties might require either a non-universal IMF or a much more complex dynamical history.

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.

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.

We used data from the HST UV Legacy Survey of Galactic Globular Clusters to select the Blue Straggler Star (BSS) population in four intermediate/high density systems (namely NGC 2808, NGC 6388, NGC 6541 and NGC 7078) through a "UV-guided search". This procedure consists in using the F275W images in each cluster to construct the master list of detected sources, and then force it to the images acquired in the other filters. Such an approach optimizes the detection of relatively hot stars and allows the detection of complete sample of BSSs even in the central region of high-density clusters, because the light from the bright cool giants, which dominates the optical emission in old stellar systems, is sensibly reduced at UV wavelengths. Our UV-guided selections of BSSs have been compared to the samples obtained in previous, optical-driven surveys, clearly demonstrating the efficiency of the UV approach. In each cluster we also measured the parameter A+, defined as the area enclosed between the cumulative radial distribution of BSSs and that of a reference population, which traces the level of BSS central segregation and the level of dynamical evolution suffered by the system. The values measured for the four clusters studied in this paper nicely fall along the dynamical sequence recently presented for a sample of 25 clusters.

We present N-body simulations of star clusters that initially evolve within a strong compressive tidal field and then transition into extensive tidal fields of varying strengths. While subject to compressive tides, clusters can undergo significant heating due to two-body interactions and mass loss due to stellar evolution. When the cluster transitions into an extensive tidal field it is super-virialized, which leads to a rapid expansion and significant mass loss before the cluster reaches virial equilibrium. After the transition, clusters are significantly less massive, more extended and therefore more tidally filling than clusters which have spent their entire lifetime in a similar extensive tidal field.

The formation and evolution of X-ray sources in globular clusters is likely to be affected by the cluster internal dynamics and the stellar interactions in the cluster dense environment.Several observational studies have revealed a correlation between the number of X-ray sources and the stellar encounter rate and provided evidence of the role of dynamics in the formation of X-ray binaries. We have performed a survey of Monte-Carlo simulations aimed at exploring the connection between the dynamics and formation of cataclysmic variables (CVs) and the origin of the observed correlation between the number of these objects, $N_{\rm cv}$, and the stellar encounter rate, $\Gamma$.The results of our simulations show a correlation between $N_{\rm cv}$ and $\Gamma$ as found in observational data, illustrate the essential role played by dynamics, and shed light on the dynamical history behind this correlation. CVs in our simulations are more centrally concentrated than single stars with masses close to those of turn-off stars, although this trend is stronger for CVs formed from primordial binaries undergoing exchange encounters, which include a population of more massive CVs absent in the group of CVs formed from binaries not suffering any component exchange.

The Hubble Space Telescope UV Legacy survey of Galactic Globular Clusters (GC) is characterising many different aspects of their multiple stellar populations. The "Grundahl-jump" (G-jump) is a discontinuity in ultraviolet brightness of blue horizontal branch (HB) stars, signalling the onset of radiative metal levitation. The HB Legacy data confirmed that the G-jump is located at the same T$_{eff}$ ($\simeq$11,500 K) in nearly all clusters. The only exceptions are the metal-rich clusters NGC 6388 and NGC 6441, where the G-jump occurs at T$_{eff}\simeq$13-14,000K. We compute synthetic HB models based on new evolutionary tracks including the effect of helium diffusion, and approximately accounting for the effect of metal levitation in a stable atmosphere. Our models show that the G-jump location depends on the interplay between the timescale of diffusion and the timescale of the evolution in the T$_{eff}$ range 11,500 K$\lessapprox$T$_{eff}\lessapprox$14,000 K. The G-jump becomes hotter than 11,500 K only for stars that have, in this T$_{eff}$ range, a helium mass fraction Y>0.35. Similarly high Y values are also consistent with the modelling of the HB in NGC 6388 and NGC 6441. In these clusters we predict that a significant fraction of HB stars show helium in their spectra above 11,500 K, and full helium settling should only be found beyond the hotter G-jump.

We present the results of a set of N-body simulations aimed at exploring how the process of mass segregation (as traced by the spatial distribution of blue straggler stars, BSSs) is affected by the presence of a population of heavy dark remnants (as neutron stars and black holes). To this end, clusters characterized by different initial concentrations and different fractions of dark remnants have been modeled. We find that an increasing fraction of stellar-mass black holes significantly delays the mass segregation of BSSs and the visible stellar component. In order to trace the evolution of BSS segregation, we introduce a new parameter ($A^+$) that can be easily measured when the cumulative radial distribution of these stars and a reference population are available. Our simulations show that $A^+$ might also be used as an approximate indicator of the time remaining to the core collapse of the visible component.

We use high-precision photometry of red-giant-branch (RGB) stars in 57 Galactic globular clusters (GCs), mostly from the `Hubble Space Telescope (HST) UV Legacy Survey of Galactic globular clusters', to identify and characterize their multiple stellar populations. For each cluster the pseudo two-color diagram (or `chromosome map') is presented, built with a suitable combination of stellar magnitudes in the F275W, F336W, F438W and F814W filters that maximizes the separation between multiple populations. In the chromosome map of most GCs (Type I clusters), stars separate in two distinct groups that we identify with the first (1G) and the second generation (2G). This identification is further supported by noticing that 1G stars have primordial (oxygen-rich, sodium-poor) chemical composition, whereas 2G stars are enhanced in sodium and depleted in oxygen. This 1G-2G separation is not possible for a few GCs where the two sequences have apparently merged into an extended, continuous sequence. In some GCs (Type II clusters) the 1G and/or the 2G sequences appear to be split, hence displaying more complex chromosome maps. These clusters exhibit multiple SGBs also in purely optical color-magnitude diagrams, with the fainter SGB joining into a red RGB which is populated by stars with enhanced heavy-element abundance. We measure the RGB width by using appropriate colors and pseudo-colors. When the metallicity dependence is removed, the RGB width correlates with the cluster mass. The fraction of 1G stars ranges from ~8% to ~67% and anticorrelates with the cluster mass, indicating that incidence and complexity of the multiple population phenomenon both increase with cluster mass.

We make use of $N$-body simulations to determine the relationship between two observable parameters that are used to quantify mass segregation and energy equipartition in star clusters. Mass segregation can be quantified by measuring how the slope of a cluster's stellar mass function $\alpha$ changes with clustercentric distance r, and then calculating $\delta_\alpha = \frac{d \alpha(r)}{d ln(r/r_m)}$ where $r_m$ is the cluster's half-mass radius. The degree of energy equipartition in a cluster is quantified by $\eta$, which is a measure of how stellar velocity dispersion $\sigma$ depends on stellar mass m via $\sigma(m) \propto m^{-\eta}$. Through a suite of $N$-body star cluster simulations with a range of initial sizes, binary fractions, orbits, black hole retention fractions, and initial mass functions, we present the co-evolution of $\delta_\alpha$ and $\eta$. We find that measurements of the global $\eta$ are strongly affected by the radial dependence of $\sigma$ and mean stellar mass and the relationship between $\eta$ and $\delta_\alpha$ depends mainly on the cluster's initial conditions and the tidal field. Within $r_m$, where these effects are minimized, we find that $\eta$ and $\delta_\alpha$ initially share a linear relationship. However, once the degree of mass segregation increases such that the radial dependence of $\sigma$ and mean stellar mass become a factor within $r_m$, or the cluster undergoes core collapse, the relationship breaks down. We propose a method for determining $\eta$ within $r_m$ from an observational measurement of $\delta_\alpha$. In cases where $\eta$ and $\delta_\alpha$ can be measured independently, this new method offers a way of measuring the cluster's dynamical state.

A number of recent observational studies of Galactic globular clusters have measured the variation in the slope of a cluster's stellar mass function $\alpha$ with clustercentric distance $r$. In order to gather a deeper understanding of the information contained in such observations, we have explored the evolution of $\alpha(r)$ for star clusters with a variety of initial conditions using a large suite of $N$-body simulations. We have specifically studied how the time evolution of $\alpha(r)$ is affected by initial size, mass, binary fraction, primordial mass segregation, black hole retention, an external tidal field, and the initial mass function itself. Previous studies have shown that the evolution of $\alpha_G$ is closely related to the amount of mass loss suffered by a cluster. Hence for each simulation we have also followed the evolution of the slope of the cluster's global stellar mass function, $\alpha_G$, and have shown that clusters follow a well-defined track in the $\alpha_G$-$d\alpha(r)/d(ln(r/r_m))$ plane. The location of a cluster on the $\alpha_G-d\alpha(r)/d(ln(r/r_m))$ plane can therefore constrain its dynamical history and, in particular, constrain possible variations in the stellar initial mass function. The $\alpha_G$-$d\alpha(r)/d(ln(r/r_m))$ plane thus serves as a key tool for fully exploiting the information contained in wide field studies of cluster stellar mass functions.

Due to their extreme helium abundance, the multiple stellar populations of the globular cluster NGC 2808 have been widely investigated from a photometric, spectroscopic, and kinematic perspective. The most striking feature of the color-magnitude diagram of NGC 2808 is the triple main sequence (MS), with the red MS corresponding to a stellar population with primordial helium, and the middle and the blue MS being enhanced in helium up to Y$\sim$0.32 and $\sim$0.38, respectively. A recent study has revealed that this massive cluster hosts at least five distinct stellar populations (A, B, C, D, and E). Among them populations A, B, and C correspond to the red MS, while populations C and D are connected to the middle and the blue MS. In this paper we exploit Hubble-Space-Telescope photometry to investigate the radial distribution of the red, the middle and the blue MS from the cluster center out to about 8.5 arcmin. Our analysis shows that the radial distribution of each of the three MSs is different. In particular, as predicted from multiple-population formation models, both the blue MS and the middle MS appears to be more concentrated than the red MS with a significance level for this result wich is above 3{\sigma}.

We present an observational estimate of the fraction and distribution of dark mass in the innermost region of the two Galactic globular clusters NGC 6218 (M12) and NGC 288. Such an assessment has been made by comparing the dynamical and luminous mass profiles derived from an accurate analysis of the most extensive spectroscopic and photometric surveys performed on these stellar systems. We find that non-luminous matter constitutes more than 60% of the total mass in the region probed by our data (R<1.6 arcmin~r_h) in both clusters. We have carefully analyzed the effects of binaries and tidal heating on our estimate and ruled out the possibility that our result is a spurious consequence of these effects. The dark component appears to be more concentrated than the most massive stars suggesting that it is likely composed of dark remnants segregated in the cluster core.

We study the interaction of the early spherical GC wind powered by Type II supernovae (SNe II) with the surrounding ambient medium consisting of the gaseous disk of a star forming galaxy at redshift z ~> 2. The bubble formed by the wind eventually breaks out of the disk, and most of the wind moves directly out of the galaxy and is definitively lost. The fraction of the wind moving nearly parallel to the galactic plane carves a hole in the disk which will contract after the end of the SN activity. During the interval of time between the end of the SN explosions and the "closure" of the hole, very O-poor stars (the Extreme population) can form out of the super-AGB (asymptotic giant branch) ejecta collected in the GC center. Once the hole contracts, the AGB ejecta mix with the pristine gas, allowing the formation of stars with an oxygen abundance intermediate between that of the very O-poor stars and that of the pristine gas. We show that this mechanism may explain why Extreme populations are present only in massive clusters, and can also produce a correlation between the spread in helium and the cluster mass. Finally, we also explore the possibility that our proposed mechanism can be extended to the case of multiple populations showing bimodality in the iron content, with the presence of two populations characterized by a small difference in [Fe/H]. Such a result can be obtained taking into account the contribution of delayed SN II.

The presence of an external tidal field often induces significant dynamical evolutionary effects on the internal kinematics of star clusters. Previous studies investigating the restricted three-body problem with applications to star cluster dynamics have shown that unbound stars on retrograde orbits (with respect to the direction of the cluster's orbit) are more stable against escape than prograde orbits, and predicted that a star cluster might acquire retrograde rotation through preferential escape of stars on prograde orbits. In this study we present evidence of this prediction, but we also illustrate that there are additional effects that cannot be accounted for by the preferential escape of prograde orbits alone. Specifically, in the early evolution, initially underfilling models increase their fraction of retrograde stars without losing significant mass, and acquire a retrograde angular velocity. We attribute this effect to the development of preferentially eccentric/radial orbits in the outer regions of star clusters as they are expanding into their tidal limitation. We explore the implications of the evolution of the fraction of prograde and retrograde stars for the evolution of the cluster internal rotation, and its dependence on the initial structural properties. Although all the systems studied here evolve towards an approximately solid-body internal rotation with angular velocity equal to about half of the angular velocity of the cluster orbital motion around the host galaxy, the evolutionary history of the radial profile of the cluster internal angular velocity depends on the cluster initial structure.

We present the results of a survey of N-body simulations aimed at exploring the evolution of compact binaries in multiple-population globular clusters.We show that as a consequence of the initial differences in the structural properties of the first-generation (FG) and the second-generation (SG) populations and the effects of dynamical processes on binary stars, the SG binary fraction decreases more rapidly than that of the FG population. The difference between the FG and SG binary fraction is qualitatively similar to but quantitatively smaller than that found for wider binaries in our previous investigations.The evolution of the radial variation of the binary fraction is driven by the interplay between binary segregation, ionization and ejection. Ionization and ejection counteract in part the effects of mass segregation but for compact binaries the effects of segregation dominate and the inner binary fraction increases during the cluster evolution. We explore the variation of the difference between the FG and the SG binary fraction with the distance from the cluster centre and its dependence on the binary binding energy and cluster structural parameters. The difference between the binary fraction in the FG and the SG populations found in our simulations is consistent with the results of observational studies finding a smaller binary fraction in the SG population.

We explain the multiple populations recently found in the 'prototype' Globular Cluster (GC) NGC 2808 in the framework of the asymptotic giant branch (AGB) scenario. The chemistry of the five -or more- populations is approximately consistent with a sequence of star formation events, starting after the supernovae type II epoch, lasting approximately until the time when the third dredge up affects the AGB evolution (age ~90-120Myr), and ending when the type Ia supernovae begin exploding in the cluster, eventually clearing it from the gas. The formation of the different populations requires episodes of star formation in AGB gas diluted with different amounts of pristine gas. In the nitrogen-rich, helium-normal population identified in NGC 2808 by the UV Legacy Survey of GCs, the nitrogen increase is due to the third dredge up in the smallest mass AGB ejecta involved in the star formation of this population. The possibly-iron-rich small population in NGC 2808 may be a result of contamination by a single type Ia supernova. The NGC 2808 case is used to build a general framework to understand the variety of 'second generation' stars observed in GCs. Cluster-to-cluster variations are ascribed to differences in the effects of the many processes and gas sources which may be involved in the formation of the second generation. We discuss an evolutionary scheme, based on pollution by delayed type II supernovae, which accounts for the properties of s-Fe-anomalous clusters.

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.

Binary fraction and orbital characteristics provide indications on the conditions of star formation, as they shed light on the environment they were born in. Multiple systems are more common in low density environments rather than in higher density ones. In the current debate about the formation of Globular Clusters and their multiple populations, studying the binary incidence in the populations they host offers a crucial piece of information on the environment of their birth and their subsequent dynamical evolution. Through a multi-year observational campaign using FLAMES at VLT, we monitored the radial velocity of 968 Red-Giant Branch stars located around the half-light radii in a sample of 10 Galactic Globular Clusters. We found a total of 21 radial velocity variables identified as {\it bona fide} binary stars, for a binary fraction of 2.2%$\pm$0.5%. When separating the sample into first generation and second generation stars, we find a binary fraction of 4.9%$\pm$1.3% and 1.2%$\pm$0.4% respectively. Through simulations that take into account possible sources of bias in detecting radial velocity variations in the two populations, we show that the difference is significant and only marginally affected by such effects. Such a different binary fraction strongly suggests different conditions in the environment of formation and evolution of first and second generations stars, with the latter being born in a much denser environment. Our result hence strongly supports the idea that the second generation forms in a dense subsystem at the center of the loosely distributed first generation, where (loose) binaries are efficiently destroyed.

Multiple or extended turnoffs in young clusters in the Magellanic Clouds have recently received large attention. A number of studies have shown that they may be interpreted as the result of a significant age spread (several 10^8yr in clusters aged 1--2 Gyr), while others attribute them to a spread in stellar rotation. We focus on the cluster NGC 1856, showing a splitting in the upper part of the main sequence, well visible in the color m_{F336W}-m_{F555W}$, and a very wide turnoff region. Using population synthesis available from the Geneva stellar models, we show that the cluster data can be interpreted as superposition of two main populations having the same age (~350Myr), composed for 2/3 of very rapidly rotating stars, defining the upper turnoff region and the redder main sequence, and for 1/3 of slowly/non-rotating stars. Since rapid rotation is a common property of the B-A type stars, the main question raised by this model concerns the origin of the slowly/non-rotating component. Binary synchronization is a possible process behind the slowly/non-rotating population; in this case, many slowly/non-rotating stars should still be part of binary systems with orbital periods in the range from 4 to 500 days. Such periods imply that Roche lobe overflow occurs, during the evolution of the primary off the main sequence, so most primaries may not be able to ignite core helium burning, consistently why the lack of a red clump progeny of the slowly rotating population.

Numerous observational studies have revealed the ubiquitous presence of multiple stellar populations in globular clusters and cast many hard challenges for the study of the formation and dynamical history of these stellar systems. In this Letter we present the results of a study of the kinematic properties of multiple populations in NGC 2808 based on high-precision Hubble Space Telescope proper-motion measurements. In a recent study, Milone et al. have identified five distinct populations (A, B, C, D, and E) in NGC 2808. Populations D and E coincide with the helium-enhanced populations in the middle and the blue main sequences (mMS and bMS) previously discovered by Piotto et al.; populations A, B, and C correspond to the redder main sequence (rMS) that in the Piotto et al. was associated with the primordial stellar population. Our analysis shows that, in the outermost regions probed (between about 1.5 and 2 times the cluster half-light radius), the velocity distribution of populations D and E is radially anisotropic (the deviation from an isotropic distribution is significant at the ~3.5-sigma level). Stars of populations D and E have a smaller tangential velocity dispersion than those of populations A, B, and C, while no significant differences are found in the radial-velocity dispersion. We present the results of a numerical simulation showing that the observed differences between the kinematics of these stellar populations are consistent with the expected kinematic fingerprint of the diffusion towards the cluster outer regions of stellar populations initially more centrally concentrated.

Horizontal Branch stars belong to an advanced stage in the evolution of the oldest stellar galactic population, occurring either as field halo stars or grouped in globular clusters. The discovery of multiple populations in these clusters, that were previously believed to have single populations gave rise to the currently accepted theory that the hottest horizontal branch members (the blue hook stars, which had late helium-core flash ignition, followed by deep mixing) are the progeny of a helium-rich "second generation" of stars. It is not known why such a supposedly rare event (a late flash followed by mixing) is so common that the blue hook of {\omega} Cen contains \sim 30% of horizontal branch stars 10 , or why the blue hook luminosity range in this massive cluster cannot be reproduced by models. Here we report that the presence of helium core masses up to \sim 0.04 solar masses larger than the core mass resulting from evolution is required to solve the luminosity range problem. We model this by taking into account the dispersion in rotation rates achieved by the progenitors, whose premain sequence accretion disc suffered an early disruption in the dense environment of the cluster's central regions where second-generation stars form. Rotation may also account for frequent late-flash-mixing events in massive globular clusters.

In this paper we combine WFC3/UVIS F275W, F336W, and F438W data from the "UV Legacy Survey of Galactic Globular Clusters: Shedding Light on Their Populations and Formation" (GO-13297) HST Treasury program with F606W, F625W, F658N, and F814W ACS archive data for a multi-wavelength study of the globular cluster NGC 6352. In the color-magnitude and two-color diagrams obtained with appropriate combination of the photometry in the different bands we separate two distinct stellar populations and trace them from the main sequence to the subgiant, red giant, horizontal and asymptotic giant branches. We infer that the two populations differ in He by Delta Y=0.029+/-0.006. With a new method, we also estimate the age difference between the two sequences. Assuming no difference in [Fe/H] and [alpha/Fe], and the uncertainties on Delta Y, we found a difference in age between the two populations of 10+/-120 Myr. If we assume [Fe/H] and [alpha/Fe] differences of 0.02 dex (well within the uncertainties of spectroscopic measurements), the total uncertainty in the relative age rises to ~300 Myr.

NGC 5053 provides a rich environment to test our understanding of the complex evolution of globular clusters (GCs). Recent studies have found that this cluster has interesting morphological features beyond the typical spherical distribution of GCs, suggesting that external tidal effects have played an important role in its evolution and current properties. Additionally, simulations have shown that NGC 5053 could be a likely candidate to belong to the Sagittarius dwarf galaxy (Sgr dSph) stream. Using the Wisconsin-Indiana-Yale-NOAO-Hydra multi-object spectrograph, we have collected high quality (signal-to-noise ratio $\sim$ 75-90), medium-resolution spectra for red giant branch stars in NGC 5053. Using these spectra we have measured the Fe, Ca, Ti, Ni, Ba, Na, and O abundances in the cluster. We measure an average cluster [Fe/H] abundance of -2.45 with a standard deviation of 0.04 dex, making NGC 5053 one of the most metal-poor GCs in the Milky Way (MW). The [Ca/Fe], [Ti/Fe], and [Ba/Fe] we measure are consistent with the abundances of MW halo stars at a similar metallicity, with alpha-enhanced ratios and slightly depleted [Ba/Fe]. The Na and O abundances show the Na-O anti-correlation found in most GCs. From our abundance analysis it appears that NGC 5053 is at least chemically similar to other GCs found in the MW. This does not, however, rule out NGC 5053 being associated with the Sgr dSph stream.

The discovery of multiple stellar populations in globular clusters has implications for all the aspects of the study of these stellar systems. In this paper, by means of N-body simulations, we study the evolution of binary stars in multiple-population clusters and explore the implications of the initial differences in the spatial distribution of different stellar populations for the evolution and survival of their binary stars. Our simulations show that initial differences between the spatial distribution of first-generation (FG) and second-generation (SG) stars can leave a fingerprint in the current properties of the binary population. SG binaries are disrupted more efficiently than those of the FG population resulting in a global SG binary fraction smaller than that of the FG. As for surviving binaries, dynamical evolution produces a difference between the SG and the FG binary binding energy distribution with the SG population characterized by a larger fraction of high binding energy (more bound) binaries. We have also studied the dependence of the binary properties on the distance from the cluster centre. Although the global binary fraction decreases more rapidly for the SG population, the local binary fraction measured in the cluster inner regions may still be dominated by SG binaries. The extent of the differences between the surviving FG and SG binary binding energy distribution also varies radially within the cluster and is larger in the cluster inner regions.

Aluminum is the heaviest light element displaying large star--to--star variations in Galactic globular clusters (GCs). This element may provide additional insight into the origin of the multiple populations now known to be a common place in GCs, and also the nature of the first-generation stars responsible for a cluster's chemical inhomogeneities. In a previous analysis we found that, unlike more metal-poor GCs, 47 Tuc did not exhibit a strong Na-Al correlation, which motivates a careful study of the similar metallicity but less massive GC M71. We present chemical abundances of O, Na, Al, and Fe for 33 giants in M71 using spectra obtained with the WIYN-Hydra spectrograph. Our spectroscopic analysis finds that, similar to 47 Tuc and in contrast with more metal-poor GCs, M71 stars do not exhibit a strong Na-Al correlation and span a relatively narrow range in [Al/Fe], which are characteristics that GC formation models must reproduce.

We present high-precision multi-band photometry for the globular cluster (GC) M2. We combine the analysis of the photometric data obtained from the Hubble Space Telescope UV Legacy Survey of Galactic GCs GO-13297, with chemical abundances by Yong et al.(2014), and compare the photometry with models in order to analyze the multiple stellar sequences we identified in the color-magnitude diagram (CMD). We find three main stellar components, composed of metal-poor, metal-intermediate, and metal-rich stars (hereafter referred to as population A, B, and C, respectively). The components A and B include stars with different $s$-process element abundances. They host six sub-populations with different light-element abundances, and exhibit an internal variation in helium up to Delta Y~0.07 dex. In contrast with M22, another cluster characterized by the presence of populations with different metallicities, M2 contains a third stellar component, C, which shows neither evidence for sub-populations nor an internal spread in light-elements. Population C does not exhibit the typical photometric signatures that are associated with abundance variations of light elements produced by hydrogen burning at hot temperatures. We compare M2 with other GCs with intrinsic heavy-element variations and conclude that M2 resembles M22, but it includes an additional stellar component that makes it more similar to the central region of the Sagittarius galaxy, which hosts a GC (M54) and the nucleus of the Sagittarius galaxy itself.

We present semi-analytical models and simplified $N$-body simulations with $10^4$ and $10^5$ particles aimed at probing the role of dynamical friction (DF) in determining the radial distribution of Blue Straggler Stars (BSSs) in globular clusters. The semi-analytical models show that DF (which is the only evolutionary mechanism at work) is responsible for the formation of a bimodal distribution with a dip progressively moving toward the external regions of the cluster. However, these models fail to reproduce the formation of the long-lived central peak observed in all dynamically evolved clusters. The results of $N$-body simulations confirm the formation of a sharp central peak, which remains as a stable feature over the time regardless of the initial concentration of the system. In spite of a noisy behavior, a bimodal distribution forms in many cases, with the size of the dip increasing as a function of time. In the most advanced stages the distribution becomes monotonic. These results are in agreement with the observations. Also the shape of the peak and the location of the minimum (which in most of the cases is within 10 core radii) turn out to be consistent with observational results. For a more detailed and close comparison with observations, including a proper calibration of the timescales of the dynamical processes driving the evolution of the BSS spatial distribution, more realistic simulations will be necessary.

In this paper we describe a new UV-initiative HST project (GO-13297) that will complement the existing F606W and F814W database of the ACS Globular Cluster (GC) Treasury by imaging most of its clusters through UV/blue WFC3/UVIS filters F275W, F336W and F438W. This "magic trio" of filters has shown an uncanny ability to disentangle and characterize multiple-population (MP) patterns in GCs in a way that is exquisitely sensitive to C, N, and O abundance variations. Combination of these passbands with those in the optical also gives the best leverage for measuring helium enrichment. The dozen clusters that had previously been observed in these bands exhibit a bewildering variety of MP patterns, and the new survey will map the full variance of the phenomenon. The ubiquity of multiple stellar generations in GCs has made the formation of these cornerstone objects more intriguing than ever; GC formation and the origin of their MPs have now become one and the same problem. In the present paper we will describe the data base and our data reduction strategy, as well as the uses we intend to make of the final photometry, astrometry, and proper motions. We will also present preliminary color-magnitude diagrams from the data so far collected. These diagrams also draw on data from GO-12605 and GO-12311, which served as a pilot project for the present GO-13297.

We investigate the viability of a model in which the chemical anomalies among Globular Cluster stars are due to accretion of gas onto the protostellar discs of low mass stars. This model has been suggested as a way to reduce the large initial cluster masses required by other models for the formation of multiple stellar generations. We numerically follow the evolution of the accreting stars, and we show that the structure of the seed star does not remain fully convective for the whole duration of the accretion phase. Stellar populations showing discrete abundances of helium in the core, that seem to be present in some clusters, might be formed with this mechanism only if accretion occurs before the core of the stars become radiative (within 2-3Myr) or if a thermohaline instability is triggered, to achieve full mixing after the accretion phase ends. We also show that the lithium abundances in accreted structures may vary by orders of magnitude in equal masses obtained by accreting different masses. In addition, the same thermohaline mixing which could provide a homogeneous helium distribution down to the stellar center, would destroy any lithium surviving in the envelope, so that both helium homogeneity and lithium survival require that the accretion phase be limited to the first couple of million years of the cluster evolution. Such a short accretion phase strongly reduces the amount of processed matter available, and reintroduces the requirement of an extremely large initial mass for the protocluster.

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.

We use the combination of photometric and spectroscopic data of 47 Tuc stars to reconstruct the possible formation of a second generation of stars in the central regions of the cluster, from matter ejected from massive Asymptotic Giant Branch stars, diluted with pristine gas. The yields from massive AGB stars with the appropriate metallicity (Z=0.004, i.e. [Fe/H]=-0.75) are compatible with the observations, in terms of extension and slope of the patterns observed, involving oxygen, nitrogen, sodium and aluminium. Based on the constraints on the maximum helium of 47 Tuc stars provided by photometric investigations, and on the helium content of the ejecta, we estimate that the gas out of which second generation stars formed was composed of about one-third of gas from intermediate mass stars, with M>= 5Mo and about two-thirds of pristine gas. We tentatively identify the few stars whose Na, Al and O abundances resemble the undiluted AGB yields with the small fraction of 47 Tuc stars populating the faint subgiant branch. From the relative fraction of first and second generation stars currently observed, we estimate that the initial FG population in 47 Tuc was about 7.5 times more massive than the cluster current total mass.

We perform a series of simulations of evolving star clusters using AMUSE (the Astrophysical Multipurpose Software Environment), a new community-based multi-physics simulation package, and compare our results to existing work. These simulations model a star cluster beginning with a King model distribution and a selection of power-law initial mass functions, and contain a tidal cut-off. They are evolved using collisional stellar dynamics and include mass loss due to stellar evolution. After determining that the differences between AMUSE results and prior publications are understood, we explored the variation in cluster lifetimes due to the random realization noise introduced by transforming a King model to specific initial conditions. This random realization noise can affect the lifetime of a simulated star cluster by up to 30%. Two modes of star cluster dissolution were identified: a mass evolution curve that contains a run-away cluster dissolution with a sudden loss of mass, and a dissolution mode that does not contain this feature. We refer to these dissolution modes as "dynamical" and "relaxation" dominated respectively. For Salpeter-like initial mass functions, we determined the boundary between these two modes in terms of the dynamical and relaxation time scales.

We used a proper combination of high-resolution HST observations and wide-field ground based data to derive the radial star density profile of 26 Galactic globular clusters from resolved star counts (which can be all freely downloaded on-line). With respect to surface brightness (SB) profiles (which can be biased by the presence of sparse, bright stars), star counts are considered to be the most robust and reliable tool to derive cluster structural parameters. For each system a detailed comparison with both King and Wilson models has been performed and the most relevant best-fit parameters have been obtained. This is the largest homogeneous catalog collected so far of star count profiles and structural parameters derived therefrom. The analysis of the data of our catalog has shown that: (1) the presence of the central cusps previously detected in the SB profiles of NGC 1851, M13 and M62 is not confirmed; (2) the majority of clusters in our sample are fitted equally well by the King and the Wilson models; (3) we confirm the known relationship between cluster size (as measured by the effective radius) and galactocentric distances; (4) the ratio between the core and the effective radii shows a bimodal distribution, with a peak at ~ 0.3 for about 80% of the clusters, and a secondary peak at ~ 0.6 for the remaining 20%. Interestingly, the main peak turns out to be in agreement with what expected from simulations of cluster dynamical evolution and the ratio between these two radii well correlates with an empirical dynamical age indicator recently defined from the observed shape of blue straggler star radial distribution, thus suggesting that no exotic mechanisms of energy generation are needed in the cores of the analyzed clusters.

We examine the photometric data for Fornax clusters, focussing our attention on their horizontal branch color distribution and, when available, on the RR Lyr variables fraction and period distribution. Based on our understanding of the HB morphology in terms of varying helium content in the context of multiple stellar generations, we show that clusters F2, F3 and F5 must contain substantial fractions of second generation stars (~54-65%). On the basis of a simple chemical evolution model we show that the helium distribution in these clusters can be reproduced by models with cluster initial masses ranging from values equal to ~4 to ~10 times larger than the current masses. Models with a very short second generation star formation episode can also reproduce the observed helium distribution but require larger initial masses up to about twenty times the current mass. While the lower limit of this range of possible initial GC masses is consistent with those suggested by the observations of the low metallicity field stars, we also discuss the possibility that the metallicity scale of field stars (based on CaII triplet spectroscopy) and the metallicities derived for the clusters in Fornax may not be consistent with each other. The reproduction of the HB morphology in F2,F3,F5 requires two interesting hypotheses: 1) the first generation HB stars lie all at "red" colours. According to this interpretation, the low metallicity stars in the field of Fornax, populating the HB at colours bluer than the blue side ((V-I)o<=0.3 or (B-V)o<=0.2) of the RR Lyrs, should be second generation stars born in the clusters;a preliminary analysis of available colour surveys of Fornax field provides a fraction ~20% of blue HB stars, in the low metallicity range; 2) the mass loss from individual second generation red giants is a few percent of a solar mass larger than the mass loss from first generation stars.

We describe the work of the Meta-Institute for Computational Astrophysics (MICA), the first professional scientific organization based in virtual worlds. MICA was an experiment in the use of this technology for science and scholarship, lasting from the early 2008 to June 2012, mainly using the Second Life and OpenSimulator as platforms. We describe its goals and activities, and our future plans. We conducted scientific collaboration meetings, professional seminars, a workshop, classroom instruction, public lectures, informal discussions and gatherings, and experiments in immersive, interactive visualization of high-dimensional scientific data. Perhaps the most successful of these was our program of popular science lectures, illustrating yet again the great potential of immersive VR as an educational and outreach platform. While the members of our research groups and some collaborators found the use of immersive VR as a professional telepresence tool to be very effective, we did not convince a broader astrophysics community to adopt it at this time, despite some efforts; we discuss some possible reasons for this non-uptake. On the whole, we conclude that immersive VR has a great potential as a scientific and educational platform, as the technology matures and becomes more broadly available and accepted.

Globular star clusters that formed at the same cosmic time may have evolved rather differently from a dynamical point of view (because that evolution depends on the internal environment) through a variety of processes that tend progressively to segregate stars more massive than the average towards the cluster centre. Therefore clusters with the same chronological age may have reached quite different stages of their dynamical history (that is, they may have different dynamical ages). Blue straggler stars have masses greater than those at the turn-off point on the main sequence and therefore must be the result of either a collision or a mass-transfer event. Because they are among the most massive and luminous objects in old clusters, they can be used as test particles with which to probe dynamical evolution. Here we report that globular clusters can be grouped into a few distinct families on the basis of the radial distribution of blue stragglers. This grouping corresponds well to an effective ranking of the dynamical stage reached by stellar systems, thereby permitting a direct measure of the cluster dynamical age purely from observed properties.

In this paper we study the long-term dynamical evolution of multiple-population clusters, focusing on the evolution of the spatial distributions of the first- (FG) and second-generation (SG) stars.In previous studies we have suggested that SG stars formed from the ejecta of FG AGB stars are expected initially to be concentrated in the cluster inner regions. Here, by means of N-body simulations, we explore the time scales and the dynamics of the spatial mixing of the FG and the SG populations and their dependence on the SG initial concentration.Our simulations show that, as the evolution proceeds, the radial profile of the SG/FG number ratio, NSG/NFG, is characterized by three regions: 1) a flat inner part; 2) a declining part in which FG stars are increasingly dominant; and 3) an outer region where the NSG/NFG profile flattens again (the NSG/NFG profile may rise slightly again in the outermost cluster regions). The radial variation of NSG/NFG implies that the fraction of SG stars determined by observations covering a limited range of radial distances is not, in general, equal to the SG global fraction, (NSG/NFG)glob. The distance at which NSG/NFG equals (NSG/NFG)glob is approximately between 1 and 2 cluster half-mass radii. The results of our simulations suggest that in many Galactic globular clusters the SG should still be more spatially concentrated than the FG.[abridged]

We follow the scenario of formation of second generation stars in globular clusters by matter processed by hot bottom burning (HBB) in massive asymptotic giant branch (AGB) stars and Super AGB stars (SAGB). In the cluster NGC 2419 we assume the presence of an extreme population directly formed from the AGB and SAGB ejecta, so we can directly compare the yields for a metallicity Z=0.0003 with the chemical inventory of the cluster NGC 2419. At such a low metallicity, the HBB temperatures (well above 108K) allow a very advanced nucleosynthesis. Masses of about 6Mo deplete Mg and synthesize Si, going beyond Al, so this latter element results only moderately enhanced; sodium can not be enhanced. The models are consistent with the observations, although the predicted Mg depletion is not as strong as in the observed stars. We predict that the oxygen abundance must be depleted by a huge factor in the Mg poor stars. The HBB temperatures are close to the region where other p capture reactions on heavier nuclei become possible. We show that high potassium abundance found in Mg poor stars can be achieved during HBB, by p captures on the argon nuclei, if the relevant cross section(s) are larger than listed in the literature or if the HBB temperature is higher. Finally, we speculate that some calcium production is occurring owing to proton capture on potassium. We emphasize the importance of a strong effort to measure a larger sample of abundances in this cluster.

Several dynamical scenarios have been proposed that can lead to prompt mass segregation on the crossing time scale of a young cluster. They generally rely on cool and/or clumpy initial conditions, and are most relevant to small systems. As a counterpoint, we present a novel dynamical mechanism that can operate in relatively large, homogeneous, cool or cold systems. This mechanism may be important in understanding the assembly of large mass-segregated clusters from smaller clumps.

Star formation plays a central role in the evolution of galaxies and of the Universe as a whole. Studies of star-forming regions in the local universe have shown that star formation typically occurs in a clustered fashion. Building a coherent picture of how star clusters form and evolve is therefore critical to our overall understanding of the star formation process. Most clusters disrupt after they form, thus contributing to the field star population. However, the most massive and dense clusters remain bound and survive for a Hubble time. These globular clusters provide unique observational probes of the formation history of their host galaxies. In particular, the age and metallicity can be determined for each globular cluster individually, allowing the distribution of ages and metallicities within host galaxies to be constrained. We show how space-based UV-to-near-IR imaging covering a wide field of view (>= 20 arcmin per axis) and deep UV/Optical multi-object spectroscopy of globular cluster systems in nearby galaxies would allow one to place important new constraints on the formation history of early-type galaxies and their structural subcomponents (e.g., bulge, halo).

Globular cluster stars show chemical abundance patterns typical of hot-CNO processing. Lithium is easily destroyed by proton capture in stellar environments, so its abundance may be crucial to discriminate among different models proposed to account for multiple populations. In order to reproduce the observed O-Na anticorrelation and other patterns typical of multiple populations, the formation of second generation stars must occur from the nuclearly processed stellar ejecta, responsible of the chemical anomalies, diluted with pristine gas having the composition of first generation stars. The lithium abundance in the unprocessed gas -which is very likely to be equal to the lithium abundance emerging from the Big Bang- affects the lithium chemical patterns among the cluster stars. This paper focuses on a scenario in which processed gas is provided by asymptotic giant branch (AGB) stars. We examine the predictions of this scenario for the lithium abundances of multiple populations. We study the role of the non-negligible lithium abundance in the ejecta of massive AGB (A(Li)~2), and, at the same time, we explore how our models can constrain the extremely large ---and very model dependent--- lithium yields predicted by recent super--AGB models. We show that the super--AGB yields may be tested by examining the lithium abundances in a large set of blue main sequence stars in wCen and/or NGC2808. In addition, we examine the different model results obtained by assuming for the pristine gas either the Big Bang abundance predicted by the standard models (A(Li)=2.6-2.7), or the abundance detected at the surface of population II stars (A(Li)=2.2-2.3). Once a chemical model is well constrained, the O--Li distribution could perhaps be used to shed light on the primordial lithium abundance.

In order to account for the chemical composition of a stellar second generation (SG), Globular Clusters (GCs) evolution models based on the asymptotic giant branch (AGB) scenario so far included only the yields available for the massive AGB stars, while the possible role of super-AGB ejecta was either extrapolated or not considered. In this work, we explore the role of super-AGB ejecta using yields recently calculated by Ventura and D'Antona. Models of clusters showing extended Na-O anticorrelations, like NGC 2808, indicate that a SG formation history similar to that outlined in our previous work is required: formation of an Extreme population with very large helium content from the pure ejecta of super-AGB stars, followed by formation of an Intermediate population by dilution of stellar ejecta with pristine gas. The very O-poor Na-rich Extreme stars can be accounted for once deep-mixing is assumed in SG giants forming in a gas with helium abundance Y> 0.34, which significantly reduces the atmospheric oxygen content, while preserving the sodium abundance. On the other hand, for clusters showing a mild O-Na anticorrelation, like M 4, the use of the new yields broadens the range of SG formation routes leading to abundance patterns consistent with observations. It is shown that models in which SG stars form only from super-AGB ejecta promptly diluted with pristine gas can reproduce the observations. We discuss the variety of small helium variations occurring in this model and its relevance for the horizontal branch morphology. In some of these models the duration of the SG formation episode can be as short as \sim10 Myr; the formation time of the SG is thus compatible with the survival of a cooling flow in the GC, previous to the explosion of the SG core collapse supernovae. We also explore models with formation of multiple populations in individual bursts, each lasting no longer than \sim10 Myr.

The discovery of multiple stellar populations in globular clusters raises fundamental questions concerning the formation and dynamical history of these systems. In a previous study aimed at exploring the formation of second-generation (SG) stars from the ejecta of first-generation (FG) AGB stars, and the subsequent dynamical evolution of the cluster, we showed that SG stars are expected to form in a dense subsystem concentrated in the inner regions of the FG cluster. In this paper we explore the implications of the structural properties of multiple-population clusters, and in particular the presence of the inner SG subsystem, for the disruption of binary stars. We quantify the enhancement of the binary disruption rate due to the presence of the central SG subsystem for a number of different initial conditions. Our calculations show that SG binaries, which are assumed to be more concentrated in the cluster inner regions, are disrupted at a substantially larger rate than FG binaries. Assuming a similar initial fraction of FG and SG binaries, our dynamical study indicates that the SG population is now expected to contain a significantly smaller binary fraction than the FG population.

Recent exam of large samples of omega Cen giants shows that it shares with mono-metallic globular clusters the presence of the sodium versus oxygen anticorrelation, within each subset of stars with iron content in the range -1.9<~[Fe/H]<~-1.3. These findings suggest that, while the second generation formation history in omega Cen is more complex than that of mono-metallic clusters, it shares some key steps with those simpler cluster. In addition, the giants in the range -1.3<[Fe/H]<~-0.7 show a direct O--Na correlation, at moderately low O, but Na up to 20 times solar. These peculiar Na abundances are not shared by stars in other environments often assumed to undergo a similar chemical evolution, such as in the field of the Sagittarius dwarf galaxy. These O and Na abundances match well the yields of the massive asymptotic giant branch stars (AGB) in the same range of metallicity, suggesting that the stars at [Fe/H]>-1.3 in omega Cen are likely to have formed directly from the pure ejecta of massive AGBs of the same metallicities. This is possible if the massive AGBs of [Fe/H]>-1.3 in the progenitor system evolve when all the pristine gas surrounding the cluster has been exhausted by the previous star formation events, or the proto--cluster interaction with the Galaxy caused the loss of a significant fraction of its mass, or of its dark matter halo, and the supernova ejecta have been able to clear the gas out of the system. The absence of dilution in the metal richer populations lends further support to a scenario of the formation of second generation stars in cooling flows from massive AGB progenitors. We suggest that the entire formation of omega Cen took place in a few 10^8yr, and discuss the problem of a prompt formation of s--process elements.

The star-to-star differences in the abundance of light elements observed in the globular clusters (GCs) can be explained assuming that a second generation (SG) of stars form in the gas ejected by the asymptotic giant branch (AGB) stars belonging to a first stellar generation. However, while Na and O appear to be anticorrelated in the cluster stars, from the stellar models they turn out to be correlated into the AGB ejecta. In order to reconcile the stellar theory with the observational findings, all the GC models invoke an early dilution of AGB ejecta with pristine gas occurring during the SG formation. Despite a vast consensus about the occurrence of such dilution, the physical process behind it is still unknown. In the present paper we set some general constraints on the pristine gas dynamics and on the possible amount of pristine gas involved in the SG formation, making use of a one zone chemical model. We find that such a dilution is a necessary ingredient in the SG star formation to explain the observed abundance patterns. We confirm the conclusion of our previous works showing that clusters must have been initially much more massive. We also show that models assuming that clusters had an initial mass similar to their current one, and adopting a large fraction of pristine gas to form SG stars, fail to reproduce the observed Na-O anticorrelation and are not viable. We finally show that the dilution event should be restricted in time, rather than extended for the all duration of the SG formation.

We analyse complementary HST and SUBARU data for the globular cluster NGC 2419. We make a detailed analysis of the horizontal branch (HB), that appears composed by two main groups of stars: the luminous blue HB stars ---that extend by evolution into the RR Lyrae and red HB region--- and a fainter, extremely blue population. We examine the possible models for this latter group and conclude that a plausible explanation is that they correspond to a significant (~30 %) extreme second generation with a strong helium enhancement (Y~0.4). We also show that the color dispersion of the red giant branch is consistent with this hypothesis, while the main sequence data are compatible with it, although the large observational error blurs the possible underlying splitting. While it is common to find an even larger (50 -- 80) percentage of second generation in a globular cluster, the presence of a substantial and extreme fraction of these stars in NGC 2419 might be surprising, as the cluster is at present well inside the radius beyond which the galactic tidal field would be dominant. If a similar situation had been present in the first stages of the cluster life, the cluster would have retained its initial mass, and the percentage of second generation stars should have been quite small (up to ~10 %). Such a large fraction of extreme second generation stars implies that the system must have been initially much more massive and in different dynamical conditions than today. We discuss this issue in the light of existing models of the formation of multiple populations in globular clusters.

Surface brightness profiles of globular clusters with shallow central cusps (Sigma ~ R^v with -0.3<~ v <~ -0.05) have been associated by several recent studies with the presence of a central intermediate mass black hole (IMBH). Such shallow slopes are observed in several globular clusters thanks to the high angular resolution of Hubble Space Telescope imaging. In this Letter we evaluate whether shallow cusps are a unique signature of a central IMBH by analyzing a sample of direct N-body simulations of star clusters with and without a central IMBH. We ``observe'' the simulations as if they were HST images. Shallow cusps are common in our simulation sample: star clusters without an IMBH have v >~ -0.3 in the pre-core-collapse and core-collapse phases. Post-core-collapse clusters without an IMBH transition to steeper cusps, -0.7<~ v <~ -0.4, only if the primordial binary fraction is very small, f_{bin}< 3 per cent, and if there are few stellar-mass black holes remaining. Otherwise v values overlap the range usually ascribed to the presence of an IMBH throughout the entire duration of the simulations. In addition, measuring v is intrinsically prone to significant uncertainty, therefore typical measurement errors may lead to v > -0.3 even when <v> <~ -0.4. Overall our analysis shows that a shallow cusp is not an unequivocal signature of a central IMBH and casts serious doubts on the usefulness of measuring v in the context of the hunt for IMBHs in globular clusters.

Many observational studies have revealed the presence of multiple stellar generations in Galactic globular clusters. These studies suggest that second-generation stars make up a significant fraction of the current mass of globular clusters, with the second-generation mass fraction ranging from ~50 to 80 per cent in individual clusters. In this Letter we carry out hydrodynamical simulations to explore the dependence of the mass of second-generation stars on the initial mass and structural parameters and stellar initial mass function (IMF) of the parent cluster. We then use the results of these simulations to estimate the fraction, f_{SG,H}, of the mass of the Galactic stellar halo composed of second-generation stars that originated in globular clusters. We study the dependence of f_{SG,H} on the parameters of the initial mass function of the Galactic globular cluster system. For a broad range of initial conditions, we find that the fraction of mass of the Galactic stellar halo in second-generation stars is always small, f_{SG,H}<4-6 per cent for a Kroupa-1993 IMF and f_{SG,H}<7-9 per cent for a Kroupa-2001 IMF.

A large number of spectroscopic studies have provided evidence of the presence of multiple populations in globular clusters by revealing patterns in the stellar chemical abundances. This paper is aimed at studying the origin of these abundance patterns. We explore a model in which second generation (SG) stars form out of a mix of pristine gas and ejecta of the first generation of asymptotic giant branch stars. We first study the constraints imposed by the spectroscopic data of SG stars in globular clusters on the chemical properties of the asymptotic and super asymptotic giant branch ejecta. With a simple one-zone chemical model, we then explore the formation of the SG population abundance patterns focussing our attention on the Na-O, Al-Mg anticorrelations and on the helium distribution function. We carry out a survey of models and explore the dependence of the final SG chemical properties on the key parameters affecting the gas dynamics and the SG formation process. Finally, we use our chemical evolution framework to build specific models for NGC 2808 and M4, two Galactic globular clusters which show different patterns in the Na-O and Mg-Al anticorrelation and have different helium distributions. We find that the amount of pristine gas involved in the formation of SG stars is a key parameter to fit the observed O-Na and Mg-Al patterns. The helium distribution function for these models is in general good agreement with the observed one. Our models, by shedding light on the role of different parameters and their interplay in determining the final SG chemical properties, illustrate the basic ingredients, constraints and problems encountered in this self-enrichment scenario which must be addressed by more sophisticated chemical and hydrodynamic simulations.

We consider the horizontal branch (HB) of the Globular Cluster Terzan 5, recently shown to be split into two parts, the fainter one (delta M_K ~ 0.3mag) having a lower metallicity than the more luminous. Both features show that it contains at least two stellar populations. The separation in magnitude has been ascribed to an age difference of ~6 Gyr and interpreted as the result of an atypical evolutionary history for this cluster. We show that the observed HB morphology is also consistent with a model in which the bright HB is composed of second generation stars that are metal enriched and with a helium mass fraction larger (by delta Y ~ 0.07) than that of first generation stars populating the fainter part of the HB. Terzan 5 would therefore be anomalous, compared to most "normal" clusters hosting multiple populations, only because its second generation is strongly contaminated by supernova ejecta; the previously proposed prolonged period of star formation, however, is not required. The iron enrichment of the bright HB can be ascribed either to contamination from Type Ia supernova ejecta of the low-iron, helium rich, ejecta of the massive asympotic giant branch stars of the cluster, or to its mixing with gas, accreting on the cluster from the environment, that has been subject to fast metal enrichment due to its proximity with the galactic bulge. The model here proposed requires only a small age difference, of ~100Myr.

Spectroscopic and photometric observations show that many globular clusters host multiple stellar populations, challenging the common paradigm that globular clusters are "simple stellar populations" composed of stars of uniform age and chemical composition. The chemical abundances of second-generation (SG) stars constrain the sources of gas out of which these stars must have formed, indicating that the gas must contain matter processed through the high-temperature CNO cycle. First-generation massive Asymptotic Giant Branch (AGB) stars have been proposed as the source of this gas. In a previous study, by means of hydrodynamical and N-body simulations, we have shown that the AGB ejecta collect in a cooling flow in the cluster core, where the gas reaches high densities, ultimately forming a centrally concentrated subsystem of SG stars. In this Letter we show that the high gas density can also lead to significant accretion onto a pre-existing seed black hole. We show that gas accretion can increase the black hole mass by up to a factor of 100. The details of the gas dynamics are important in determining the actual black hole growth. Assuming a near-universal seed black hole mass and small cluster-to-cluster variations in the duration of the SG formation phase, the outcome of our scenario is one in which the present intermediate-mass black hole (IMBH) mass may have only a weak dependence on the current cluster properties. The scenario presented provides a natural mechanism for the formation of an IMBH at the cluster center during the SG star-formation phase.

[abridged] We present a unified picture for the evolution of star clusters on the two-body relaxation timescale. We use direct N-body simulations of star clusters in a galactic tidal field starting from different multi-mass King models, up to 10% of primordial binaries and up to Ntot=65536 particles. An additional run also includes a central Intermediate Mass Black Hole. We find that for the broad range of initial conditions we have studied the stellar mass function of these systems presents a universal evolution which depends only on the fractional mass loss. The structure of the system, as measured by the core to half mass radius ratio, also evolves toward a universal state, which is set by the efficiency of heating on the visible population of stars induced by dynamical interactions in the core of the system. Interactions with dark remnants are dominant over the heating induced by a moderate population of primordial binaries (3-5%), especially under the assumption that most of the neutron stars and black holes are retained in the system. All our models without primordial binaries undergo a deep gravothermal collapse in the radial mass profile. However their projected light distribution can be well fitted by medium concentration King models (with parameter W0 ~ 8), even though there tends to be an excess over the best fit for the innermost points of the surface brightness. This excess is consistent with a shallow cusp in the surface brightness (mu(R) ~ R^{-v} with v ~ 0.4-0.7), like it has been observed for many globular clusters from high-resolution HST imaging. Classification of core-collapsed globular clusters based on their surface brightness profile is likely to fail in systems that have already bounced back to lower concentrations.

Dynamical evolution plays a key role in shaping the current properties of star clusters and star cluster systems. A detailed understanding of the effects of evolutionary processes is essential to be able to disentangle the properties which result from dynamical evolution from those imprinted at the time of cluster formation. In this review, we focus our attention on globular clusters and review the main physical ingredients driving their early and long-term evolution, describe the possible evolutionary routes and show how cluster structure and stellar content are affected by dynamical evolution.

We investigate the ratio between the half-mass radii r_h of Galactic globular clusters and their Jacobi radii r_J given by the potential of the Milky Way and show that clusters with galactocentric distances R_{GC}>8 kpc fall into two distinct groups: one group of compact, tidally-underfilling clusters with r_h/r_J<0.05 and another group of tidally filling clusters which have 0.1 < r_h/r_J<0.3. We find no correlation between the membership of a particular cluster to one of these groups and its membership in the old or younger halo population. Based on the relaxation times and orbits of the clusters, we argue that compact clusters and most clusters in the inner Milky Way were born compact with half-mass radii r_h < 1 pc. Some of the tidally-filling clusters might have formed compact as well, but the majority likely formed with large half-mass radii. Galactic globular clusters therefore show a similar dichotomy as was recently found for globular clusters in dwarf galaxies and for young star clusters in the Milky Way. It seems likely that some of the tidally-filling clusters are evolving along the main sequence line of clusters recently discovered by Kuepper et al. (2008) and are in the process of dissolution.

We describe the Meta-Institute for Computational Astrophysics (MICA), the first professional scientific organization based exclusively in virtual worlds (VWs). The goals of MICA are to explore the utility of the emerging VR and VWs technologies for scientific and scholarly work in general, and to facilitate and accelerate their adoption by the scientific research community. MICA itself is an experiment in academic and scientific practices enabled by the immersive VR technologies. We describe the current and planned activities and research directions of MICA, and offer some thoughts as to what the future developments in this arena may be.

In this paper we use N-body simulations to study the effects of primordial mass segregation on the early and long-term evolution of star clusters. Our simulations show that in segregated clusters early mass loss due to stellar evolution triggers a stronger expansion than for unsegregated clusters. Tidally limited, strongly segregated clusters may dissolve rapidly as a consequence of this early expansion, while segregated clusters initially underfilling their Roche lobe can survive the early expansion and have a lifetime similar to that of unsegregated clusters. Long-lived initially segregated clusters tend to have looser structure and reach core collapse later in their evolution than initially unsegregated clusters. We have also compared the effects of dynamical evolution on the global stellar mass function (MF) of low-mass main sequence stars. In all cases the MF flattens as the cluster loses stars. The amount of MF flattening induced by a given amount of mass loss in a rapidly dissolving initially segregated cluster is less than for an unsegregated cluster. The evolution of the MF of a long-lived segregated cluster, on the other hand, is very similar to that of an initially unsegregated cluster.

Some globular clusters are observed to host a population of second generation (SG) stars which show chemical anomalies and must have formed from gas containing matter processed in the envelopes of first generation (FG) cluster stars. We study the SG formation process by means of 1D hydrodynamical simulations, assuming that the SG is formed by the gas ejected by AGB stars. This gas collects in a cooling flow into the cluster core, where it forms a SG star subsystem strongly concentrated in the cluster innermost regions with structural properties largely independent of the FG initial properties. We also present the results of a model in which pristine gas contributes to the SG formation. In this model a very helium-rich SG population and one with a moderate helium enrichment form; the resulting SG bimodal helium distribution resembles that observed for SG stars in NGC 2808. By means of N-body simulations, we study the two-population cluster dynamical evolution. In our simulations, a large fraction of FG stars are lost early in the cluster evolution due to the expansion and stripping of the cluster outer layers resulting from early mass loss associated with FG SN ejecta. The SG population is largely unscathed by this early mass loss, and this early evolution leads to values of the SG to FG number ratio consistent with observations. We also demonstrate possible evolutionary routes leading to the loss of most of the FG population, leaving an SG-dominated cluster. Until mixing of the two populations is complete, the radial profile of the SG to FG number ratio is characterized by a flat inner part and a declining portion in the outer cluster regions. (abridged)

The advent of affordable parallel computers such as Beowulf PC clusters and, more recently, of multi-core PCs has been highly beneficial for a large number of scientists and smaller institutions that might not otherwise have access to substantial computing facilities. However, there has not been an analogous progress in the development and dissemination of parallel software: scientists need the expertise to develop parallel codes and have to invest a significant amount of time in the development of tools even for the most common data analysis tasks. We describe the Beowulf Analysis Symbolic INterface (BASIN) a multi-user parallel data analysis and visualization framework. BASIN is aimed at providing scientists with a suite of parallel libraries for astrophysical data analysis along with general tools for data distribution and parallel operations on distributed data to allow them to easily develop new parallel libraries for their specific tasks.

Some young star clusters show a degree of mass segregation that is inconsistent with the effects of standard two-body relaxation from an initially unsegregated system without substructure, in virial equilibrium, and it is unclear whether current cluster formation models can account for this degree of initial segregation in clusters of significant mass. In this Letter we demonstrate that mergers of small clumps that are either initially mass segregated, or in which mass segregation can be produced by two-body relaxation before they merge, generically lead to larger systems which inherit the progenitor clumps' segregation. We conclude that clusters formed in this way are naturally mass segregated, accounting for the anomalous observations and suggesting that this process of prompt mass segregation due to initial clumping should be taken fully into account in constructing cluster dynamical models.

We present new Hubble Space Telescope (HST)-NIC3, near-infrared H-band photometry of globular clusters (GC) around NGC 4365 and NGC 1399 in combination with archival HST-WCPC2 and ACS optical data. We find that NGC 4365 has a number of globular clusters with bluer optical colors than expected for their red optical to near-infrared colors and an old age. The only known way to explain these colors is with a significant population of intermediate-age (2-8 Gyr) clusters in this elliptical galaxy. In contrast, NGC 1399 reveals no such population. Our result for NGC 1399 is in agreement with previous spectroscopic work that suggests that its clusters have a large metallicity spread and are nearly all old. In the literature, there are various results from spectroscopic studies of modest samples of NGC 4365 globular clusters. The spectroscopic data allow for either the presence or absence of a significant population of intermediate-age clusters, given the index uncertainties indicated by comparing objects in common between these studies and the few spectroscopic candidates with optical to near-IR colors indicative of intermediate ages. Our new near-IR data of the NGC 4365 GC system with much higher signal-to-noise agrees well with earlier published photometry and both give strong evidence of a significant intermediate-age component. The agreement between the photometric and spectroscopic results for NGC 1399 and other systems lends further confidence to this conclusion, and to the effectiveness of the near-IR technique.

We investigate the role of dissolution of low-concentration clusters due to mass loss through stellar evolution on the evolution of the properties of globular cluster systems (GCSs) in elliptical galaxies. Our simulations show that, for an initial mass-concentration relationship based on that inferred from Galactic globular clusters, dissolution of low-concentration clusters leads to the disruption of a large number of clusters. A power-law initial globular cluster system mass function (GCMF) similar to that observed in young cluster systems in merging galaxies is transformed by this dissolution into a bell-shaped GCMF with a mean mass similar to that of old GCSs for all the galaxies investigated. Two-body relaxation and dynamical friction, which are also included in our simulations, subsequently lead to an additional significant evolution and disruption of the population of clusters. As shown previously, when these processes act on a bell-shaped GCMF with a mean mass similar to that of old GCS, they do not significantly alter the value of the mean mass. The final GCMFs are bell-shaped with similar peaks at different radii within galaxies and between different galaxies, in agreement with current observations.

We study the dynamical evolution of the M87 globular cluster system (GCS) with a number of numerical simulations. We explore a range of different initial conditions for the GCS mass function (GCMF), for the GCS spatial distribution and for the GCS velocity distribution. We confirm that an initial power-law GCMF like that observed in young cluster systems can be readily transformed through dynamical processes into a bell-shaped GCMF. However,only models with initial velocity distributions characterized by a strong radial anisotropy increasing with the galactocentric distance are able to reproduce the observed constancy of the GCMF at all radii.We show that such strongly radial orbital distributions are inconsistent with the observed kinematics of the M87 GCS. The evolution of models with a bell-shaped GCMF with a turnover similar to that currently observed in old GCS is also investigated. We show that models with this initial GCMF can satisfy all the observational constraints currently available on the GCS spatial distribution,the GCS velocity distribution and on the GCMF properties.In particular these models successfully reproduce both the lack of a radial gradient of the GCS mean mass recently found in an analysis of HST images of M87 at multiple locations, and the observed kinematics of the M87 GCS.Our simulations also show that evolutionary processes significantly affect the initial GCS properties by leading to the disruption of many clusters and changing the masses of those which survive.The preferential disruption of inner clusters flattens the initial GCS number density profile and it can explain the rising specific frequency with radius; we show that the inner flattening observed in the M87 GCS spatial distribution can be the result of the effects of dynamical evolution on an initially steep density profile. (abridged)

We have studied the evolution of globular cluster systems (GCS) in elliptical galaxies with a power-law initial GCS mass function (GCMF) ($f(M) \propto M^{-\alpha}$) similar to that predicted by some theoretical studies of globular clusters formation and to that of young cluster systems observed in merging galaxies. We have carried out a survey over a large number of different host galaxies and we have considered different values for the index, $\alpha$, of the initial power-law GCMF ($\alpha=1.5,1.8,2.0$); we show the dependence of the main GCS final properties (mean mass and dispersion of the final GCMF, fraction of surviving clusters, radial gradient of the GCMF parameters) on the structure of the host galaxy and on the slope of the initial GCMF. For a subsample of host galaxies with values of effective masses and radii equal to those determined using observational data for a number of giant, normal and dwarf galaxies our results show that the relation between the final GCMF properties and those of the host galaxies as well as the dependence of the final GCMF parameters on the galactocentric distance within individual galaxies differ from those observed in old GCS: the values of the final GCS mean mass are in general smaller (4.2 < log M_final < 5.0) than those observed and the galaxy-to-galaxy dispersion of logM_final is larger than that reported by observational analyses. The results are compared with those of a companion paper in which we investigated the evolution of GCS with a log-normal initial GCMF and in which the final GCS properties were perfectly consistent with observations.

We study the evolution of globular cluster systems (GCS) in elliptical galaxies and explore the dependence of their main properties on the mass and the size of the host galaxy.The dependence of the evolution of the GCS mass function (GCMF), of the fraction of surviving clusters and of the ratio of the final to initial mass in clusters on the structure of the host galaxy as well as their variation with the galactocentric distance inside individual host galaxies has been thoroughly investigated.After a survey over a large number of different host galaxies we have restricted our attention to a sample of galaxies with effective masses and radii equal to those observed for dwarf,normal and giant ellipticals. We show that, in spite of large differences in the fraction of surviving clusters, the final mean masses of the GCMF in massive galaxies are very similar to each other with a small galaxy-to-galaxy dispersion;low-mass compact galaxies tend to have smaller values of the final mean mass and a larger galaxy-to-galaxy dispersion. These findings are in agreement with those of recent observational analyses. The fraction of surviving clusters increases with the mass of the host galaxy. We show that a small difference between the initial and the final mean mass and dispersion of the GCMF and the lack of a significant radial dependence of the mean mass inside individual galaxies do not necessarily imply that evolutionary processes have been unimportant in the evolution of the initial population of clusters. For giant galaxies most disruption occurs within the effective radius while for low-mass galaxies a significant disruption of clusters takes place also at larger galactocentric distances. The dependence of the results obtained on the initial mean mass of the GCMF is investigated. (abridged)

We study the internal response of a galaxy to an unbound encounter and present a survey of orbital parameters covering typical encounters in different galactic environments. Overall, we conclude that relatively weak encounters by low-mass interloping galaxies can cause observable distortions in the primaries. The resulting asymmetries may persist long after the interloper is evident. We focus our attention on the production of structure in dark halos and in cluster ellipticals. Any distortion produced in a dark halo can distort the embedded stellar disk, possibly leading to the formation of lopsided and warped disks. Features excited in an elliptical are directly observable and we predict that asymmetries in the morphologies of these systems can be produced by relatively small perturbers. We use these detailed results to predict the distribution of the $A$ parameter defined by Abraham et al. and the shift between the center of mass of the primary system and the position of the peak of density for a range of environments. We find that high-density, low-velocity dispersion environments are more likely to host galaxies with significant asymmetries. Our distribution for the $A$ parameter is in good agreement with the range spanned by the observed values for local galaxy clusters and for distant galaxies in the MDS and in the HDF.Finally, we propose a generalized asymmetry parameter $A(r)$ which provides detailed information on the radial structure of the asymmetry produced by the mechanism explored in our work. (abridged)

In this paper we investigate the gravothermal instability of spherical stellar systems endowed with a radially anisotropic velocity distribution. We focus our attention on the effects of anisotropy on the conditions for the onset of the instability and in particular we study the dependence of the spatial structure of critical models on the amount of anisotropy present in a system. The investigation has been carried out by the method of linear series which has already been used in the past to study the gravothermal instability of isotropic systems. We consider models described by King, Wilson and Woolley-Dickens distribution functions. In the case of King and Woolley-Dickens models, our results show that, for quite a wide range of amount of anisotropy in the system, the critical value of the concentration of the system (defined as the ratio of the tidal to the King core radius of the system) is approximately constant and equal to the corresponding value for isotropic systems. Only for very anisotropic systems the critical value of the concentration starts to change and it decreases significantly as the anisotropy increases and penetrates the inner parts of the system. For Wilson models the decrease of the concentration of critical models is preceded by an intermediate regime in which critical concentration increases, it reaches a maximum and then it starts to decrease. The critical value of the central potential always decreases as the anisotropy increases.

We investigate the evolution of the mass function of the Galactic globular cluster system (GCMF) taking into account the effects of stellar evolution, two-body relaxation, disk shocking and dynamical friction on the evolution of individual globular clusters.We have adopted a log-normal initial GCMF and we have studied in detail the dependence on the initial dispersion and mean value of the GCMF of the fraction of the initial number of clusters surviving after one Hubble time, of the final GCMF, and of the difference between the properties of the GCMF of clusters closer to the Galactic center and of those in the outer regions of the Galaxy.The initial log-normal shape of the GCMF is preserved in most cases and if a power-law in M is adopted for the initial GCMF, evolutionary processes tend to modify it into a log-normal GCMF. The difference between initial and final values of dispersion and mean value of the GCMF as well as the difference between the final values of these parameters for inner and outer clusters can be positive or negative depending on initial conditions.A significant effect of evolutionary processes does not necessarily give rise to a strong trend of the mean value of the GCMF with the galactocentric distance.The existence of a particular initial GCMF able to keep its initial shape and parameters unaltered during the entire evolution through a subtle balance between disruption of clusters and evolution of the masses of those which survive,suggested in Vesperini(1997),is confirmed. (abridged)

In this paper we show the results of a large set of N-body simulations modelling the evolution of globular clusters driven by relaxation,stellar evolution,disk shocking and including the effects of the tidal field of the Galaxy. We investigate the evolution of multi-mass models with a power-law initial mass function (IMF) starting with different initial masses, concentrations, slopes of the IMF and located at different galactocentric distances. We show to what extent the effects of the various evolutionary processes alter the shape of the IMF and to what extent these changes depend on the position of the cluster in the Galaxy. Both the changes in the global mass function and in the local one (measured at different distances from the cluster center) are investigated showing whether and where the local mass function keeps memory of the IMF and where it provides a good indication of the current global mass function. The evolution of the population of white dwarfs is also followed in detail and we supply an estimate of the fraction of the current value of the total mass expected to be in white dwarfs depending on the main initial conditions for the cluster (mass and position in the Galaxy).Simple analytical expression by which it is possible to calculate the main quantities of interest (total mass, fraction of white dwarfs, slope of the mass function) at any time t for a larger number of different initial conditions than those investigated numerically have been derived.

In this paper we address the issue of the origin of some observational properties of the galactic globular cluster system. After a preliminary study of some general properties of the main evolutionary processes, we investigate the evolution of systems of globular clusters located in a model of the Milky Way starting from different initial conditions. We study the role of the evolutionary processes in changing the spatial distribution and mass function of the cluster system, in establishing and/or preserving some of the observed correlations and trends between internal properties of globular clusters and between internal properties and location inside the host galaxy and we provide an estimate for the rates of core collapse and disruption of globular clusters. The initial mass function and spatial distribution of the cluster system evolve quite significantly in one Hubble time and the evolution is toward a final state similar to the observed one. If the mass function is initially taken to be a log-normal distribution similar to the one currently observed in our galaxy, its shape is not significantly altered during the entire simulation even though a significant number of clusters are disrupted before one Hubble time, which suggests that the present mass function might represent a sort of 'quasi-equilibrium' distribution.

We investigate a population of primordial binaries during the initial stage of evolution of a star cluster. For our calculations we assume that equal mass stars form rapidly in a tidally truncated gas cloud, that $\sim10$\% of the stars are in binaries and that the resulting star cluster undergoes an epoch of violent relaxation. We study the collisional interaction of the binaries and single stars, in particular, the ionization of the binaries and the energy exchange between binaries and single stars. We find that for large $N$ systems ($N > 10^3$) even the most violent beginning leaves the binary distribution function largely intact. Hence, the binding energy originally tied up in the cloud's protostellar pairs is preserved during the relaxation process and the binaries are available to interact at later times within the virialized cluster.

The Halo Galactic Globular Clusters (GCs) lying inside the solar circle (\inn, \r $< 8 Kpc$, \feh$<-0.8$) are shown to keep no record of the strong correlation between GC absolute integrated magnitudes and core parameters (concentration and central density) clearly detected for the Halo GCs lying outside this circle (\out) as well as (at lower statistical significance) for those belonging to the \dis population ($\feh >-0.8$). We present here both observational evidences and theoretical simulations supporting the hypothesis of a {\it primordial} origin for such correlations and of an {\it environment-induced} destruction of them if absent. In particular, this hypothesis lead us to argue that \inn Halo GCs have undergone a chaotic, environment-driven, dynamical evolution not shared by the \dis GCs though the radial zone of the Galaxy occupied by the two groups is nearly the same.