Stellar systems are broadly divided into collisional and non-collisional. The latter are large-N systems with long relaxation timescales and can be simulated disregarding two-body interactions, while either computationally expensive direct N-body simulations or approximate schemes are required to properly model the former. Large globular clusters and nuclear star clusters, with relaxation timescales of the order of a Hubble time, are small enough to display some collisional behaviour and big enough to be impossible to simulate with direct $N$-body codes and current hardware. We introduce a new method to simulate collisional stellar systems, and validate it by comparison with direct N-body codes on small-N simulations. The Multi-Particle collision for Dense stellar systems Code (MPCDSS) is a new code for evolving stellar systems with the Multi-Particle Collision method. Such method amounts to a stochastic collision rule that allows to conserve exactly the energy and momentum over a cluster of particles experiencing the collision. The code complexity scales with $N \log N$ in the number of particles. Unlike Monte-Carlo codes, MPCDSS can easily model asymmetric, non-homogeneous, unrelaxed and rotating systems, while allowing us to follow the orbits of individual stars. We evolve small ($N = 3.2 \times 10^4$) star clusters with MPCDSS and with the direct-summation code NBODY6, finding a similar evolution of key indicators. We then simulate different initial conditions in the $10^4 - 10^6$ star range. MPCDSS bridges the gap between small, collisional systems that can be simulated with direct $N$-body codes and large noncollisional systems. MPCDSS in principle allows us to simulate globular clusters such as Omega Cen and M54 and even the nuclear star cluster, beyond the limits of current direct N-body codes in terms of the number of particles.

We investigate the properties of the host galaxies of compact binary mergers across cosmic time. To this end, we combine population synthesis simulations together with galaxy catalogues from the hydrodynamical cosmological simulation EAGLE to derive the properties of the host galaxies of binary neutron star (BNS), black hole-neutron star (BHNS) and binary black hole (BBH) mergers. Within this framework, we derive the host galaxy probability, i.e., the probability that a galaxy hosts a compact binary coalescence as a function of its stellar mass, star formation rate, $K_s$ magnitude and $B$ magnitude. This quantity is particularly important for low-latency searches of gravitational wave (GW) sources as it provides a way to rank galaxies lying inside the credible region in the sky of a given GW detection, hence reducing the number of viable host candidates. Furthermore, even if no electromagnetic counterpart is detected, the proposed ranking criterion can still be used to classify the galaxies contained in the error box. Our results show that massive galaxies (or equivalently galaxies with a high luminosity in $K_s$ band) have a higher probability of hosting BNS, BHNS, and BBH mergers. We provide the probabilities in a suitable format to be implemented in future low-latency searches.

In globular clusters, blue straggler stars are heavier than the average star, so dynamical friction strongly affects them. The radial distribution of BSS, normalized to a reference population, appears bimodal in a fraction of Galactic GCs, with a density peak in the core, a prominent zone of avoidance at intermediate radii, and again higher density in the outskirts. The zone of avoidance appears to be located at larger radii the more relaxed the host cluster, acting as a sort of dynamical clock. We use a new method to compute the evolution of the BSS radial distribution under dynamical friction and diffusion. We evolve our BSS in the mean cluster potential under dynamical friction plus a random fluctuating force, solving the Langevin equation with the Mannella quasi symplectic scheme. This amounts to a new simulation method which is much faster and simpler than direct N-body codes but retains their main feature: diffusion powered by strong, if infrequent, kicks. We compute the radial distribution of initially unsegregated BSS normalized to a reference population as a function of time. We trace the evolution of its minimum, corresponding to the zone of avoidance. We compare the evolution under kicks extracted from a Gaussian distribution to that obtained using a Holtsmark distribution. The latter is a fat tailed distribution which correctly models the effects of close gravitational encounters. We find that the zone of avoidance moves outwards over time, as expected based on observations, only when using the Holtsmark distribution. Thus the correct representation of near encounters is crucial to reproduce the dynamics of the system. We confirm and extend earlier results that showed how the dynamical clock indicator depends both on dynamical friction and effective diffusion powered by dynamical encounters.

Young star clusters are the most common birth-place of massive stars and are dynamically active environments. Here, we study the formation of black holes (BHs) and binary black holes (BBHs) in young star clusters, by means of 6000 N-body simulations coupled with binary population synthesis. We probe three different stellar metallicities (Z=0.02, 0.002 and 0.0002) and two initial density regimes (density at the half-mass radius $\rho_{\rm h}\ge{}3.4\times10^4$ and $\ge{1.5\times10^2}$ M$_\odot$ pc$^{-3}$ in dense and loose star clusters, respectively). Metal-poor clusters tend to form more massive BHs than metal-rich ones. We find $\sim{}6$, $\sim{}2$, and $<1$% of BHs with mass $m_{\rm BH}>60$ M$_\odot$ at Z=0.0002, 0.002 and 0.02, respectively. In metal-poor clusters, we form intermediate-mass BHs with mass up to $\sim{}320$ M$_\odot$. BBH mergers born via dynamical exchanges (exchanged BBHs) can be more massive than BBH mergers formed from binary evolution: the former (latter) reach total mass up to $\sim{}140$ M$_\odot$ ($\sim{}80$ M$_\odot$). The most massive BBH merger in our simulations has primary mass $\sim{}88$ M$_\odot$, inside the pair-instability mass gap, and a mass ratio of $\sim{}0.5$. Only BBHs born in young star clusters from metal-poor progenitors can match the masses of GW170729, the most massive event in O1 and O2, and those of GW190412, the first unequal-mass merger. We estimate a local BBH merger rate density $\sim{}110$ and $\sim{}55$ Gpc$^{-3}$ yr$^{-1}$, if we assume that all stars form in loose and dense star clusters, respectively.

We investigate the properties of the host galaxies of compact binary mergers across cosmic time, by means of population-synthesis simulations combined with galaxy catalogues from the EAGLE suite. We analyze the merger rate per galaxy of binary neutron stars (BNSs), black hole--neutron star binaries (BHNSs) and binary black holes (BBHs) from redshift zero up to redshift six. The binary merger rate per galaxy strongly correlates with the stellar mass of the host galaxy at any redshift considered here. This correlation is significantly steeper for BNSs than for both BHNSs and BBHs. Moreover, we find that the merger rate per galaxy depends also on host galaxy's star formation rate and metallicity. We derive a robust fitting formula that relates the merger rate per galaxy with galaxy's star formation rate, stellar mass and metallicity at different redshifts. The typical masses of the host galaxies increase significantly as redshift decreases, as a consequence of the interplay between delay time distribution of compact binaries and cosmic assembly of galaxies. Finally, we study the evolution of the merger rate density with redshift. At low redshift ($z\leq{}0.1$) early-type galaxies give a larger contribution to the merger rate density than late-type galaxies. This trend reverts at $z\ge{}1$.

Recent interpretation of the color$-$magnitude diagrams of the Milky Way (MW) bulge has suggested that the observed double red-clump feature can be a natural consequence of He-enhanced stellar populations in the MW bulge. This implies that globular clusters (GCs), where the He-enhanced second-generation (SG) stars can be efficiently created, are the most likely candidate contributors of He-rich stars to the MW bulge. We extend this idea to the Galactic inner halo and investigate the fraction of the SG stars as a function of the Galactocentric distance. We use bluer blue-horizontal branch (bBHB) stars, which are assumed to be originated from He-rich SG populations, as proxies of SG stars, and find that the fraction of bBHB stars increases with decreasing Galactocentric distance. Simulations of the GC evolution in the MW tidal field qualitatively support the observed trend of bBHB enhancement in the inner halo. In these simulations, the increasing tidal force with decreasing Galactocentric distance leads to stripping of stars not only from the outskirts but also from the central regions of GCs, where SG stars are more abundant. We discuss the implication and prospect of our findings concerning the formation history of the bulge and inner halo of the MW.

Sizeable number of stellar-mass black holes (BHs) in globular clusters (GCs) can strongly influence the dynamical evolution and observational properties of their host cluster. Using results from a large set of numerical simulations, we identify the key ingredients needed to sustain a sizeable population of BHs in GCs up to a Hubble time. We find that while BH natal kick prescriptions are essential in determining the initial retention fraction of BHs in GCs, the long-term survival of BHs is determined by the size, initial central density and half-mass relaxation time of the GC. Simulated GC models that contain many BHs are characterized by relatively low central surface brightness, large half-light and core radii values. We also discuss novel ways to compare simulated results with available observational data to identify GCs that are most likely to contain many BHs.

Blue straggler stars are more massive than the average star in globular clusters, as they originate from the merger of two stars. Consequently, they experience dynamical friction, progressively sinking to the cluster center. Recently, several indicators of the degree of dynamical relaxation of a globular cluster have been proposed, based on the observed radial distribution of blue straggler stars. The most successful is the Alessandrini indicator, or A+ for short, which is the integral of the cumulative distribution of the blue straggler stars minus that of a lighter reference population. A+ correlates with the dynamical age of a cluster both in realistic simulations and in observations. Here I calculate the temporal dependence of the A+ indicator analytically in a simplified model of the evolution of the blue straggler star distribution under dynamical friction only.

The Alessandrini A$+$ indicator is a quantitative measure of star cluster dynamical evolution based on the mass-segregation of blue straggler stars. A$+$ is defined as the integral of the cumulative distribution of blue straggler stars with the radius measured in logarithmic scale, minus a term related to the reference population used. In a companion paper I introduced a simplified model of dynamical friction and calculated the A$+$ indicator analytically. Here I show further properties of the time evolution of A$+$, focusing in particular on the physical interpretation of its time derivative $d{A^+}/dt$. I find that $d{A^+}/dt$ is the mean of the reciprocal dyamical friction timescale, weighted by the density of blue straggler stars. I show that it is non-negative (as expected based on monotonicity) due to the density of blue-straggler stars being non-negative and that, for a dynamical friction timescale that is non-decreasing with radius, $d{A^+}/dt$ is also non-decreasing with time, making the A$+$ indicator a convex function.

In the massive globular cluster NGC 2808, RGB stars form at least five distinct groups in the so-called chromosome map photometric plane, arguably corresponding to different stellar populations. While a human expert can separate the groups by eye relatively easily, algorithmic approaches are desirable for reproducibility and for handling a larger sample of globular clusters. Unfortunately, cluster analysis algorithms often produced unsatisfactory results. Here we apply a range of non-parametric clustering algorithms to the NGC 2808 RGB dataset: partitioning (k-means, Partitioning Around Medoids - PAM), hierarchical (AGglomerative NESting - AGNES, DIvisive ANAlysis - DIANA), and density based (Density-Based Spatial Clustering of Applications with Noise - DBSCAN, Ordering Points To Identify the Clustering Struture - OPTICS). For each algorithm we discuss different choices of the relevant hyperparameters and their impact on the resulting clustering. We find that AGNES produces results that are most similar to the expectations of a human expert, depending on the prescription used for joining adjacent groups - linkage. Among the linkage prescriptions we tested, Ward's method performs best, and average linkage obtains comparable results only if outliers are removed beforehand. We recommend using AGNES with Ward's method or similar linkages in future studies to automatically identify stellar populations in the chromosome map plane.

Context: NGC 2682 is a nearby open cluster, approximately 3.5 Gyr old. Dynamically, most open clusters should dissolve on shorter timescales, of ~ 1 Gyr. Having survived until now, NGC 2682 was likely much more massive in the past, and is bound to have an interesting dynamical history. Aims: We investigate the spatial distribution of NGC 2682 stars to constrain its dynamical evolution, especially focusing on the marginally bound stars in the cluster outskirts. Methods: We use Gaia DR2 data to identify NGC 2682 members up to a distance of ~150 pc (10 degrees). Two methods (Clusterix and UPMASK) are applied to this end. We estimate distances to obtain three-dimensional stellar positions using a Bayesian approach to parallax inversion, with an appropriate prior for star clusters. We calculate the orbit of NGC 2682 using the GRAVPOT16 software. Results: The cluster extends up to 200 arcmin (50 pc) which implies that its size is at least twice as previously believed. This exceeds the cluster Hill sphere based on the Galactic potential at the distance of NGC 2682. Conclusions: The extra-tidal stars in NGC 2682 may originate from external perturbations such as disk shocking or dynamical evaporation from two-body relaxation. The former origin is plausible given the orbit of NGC 2682, which crossed the Galactic disk ~40 Myr ago.

Machine learning is a powerful technique, becoming increasingly popular in astrophysics. In this paper, we apply machine learning to more than a thousand globular cluster (GC) models simulated as part of the 'MOCCA-Survey Database I' project in order to correlate present-day observable properties with the presence of a subsystem of stellar mass black holes (BHs). The machine learning model is then applied to available observed parameters for Galactic GCs to identify which of them that are most likely to be hosting a sizeable number of BHs and reveal insights into what properties lead to the formation of BH subsystems. With our machine learning model, we were able to shortlist 21 Galactic GCs that are most likely to contain a BH subsystem. We show that the clusters shortlisted by the machine learning classifier include those in which BH candidates have been observed (M22, M10 and NGC 3201) and that our results line up well with independent simulations and previous studies that manually compared simulated GC models with observed properties of Galactic GCs. These results can be useful for observers searching for elusive stellar mass BH candidates in GCs and further our understanding of the role BHs play in GC evolution. In addition, we have released an online tool that allows one to get predictions from our model after they input observable properties.

Globular Clusters (GCs) have historically been subdivided in either two (disk/halo) or three (disk/inner-halo/outer-halo) groups based on their orbital, chemical and internal physical properties. The qualitative nature of this subdivision makes it impossible to determine whether the natural number of groups is actually two, three, or more. In this paper we use cluster analysis on the $(\log M, \log \sigma_0, \log R_e, [Fe/H], \log | Z |)$ space to show that the intrinsic number of GC groups is actually either $k=2$ or $k=3$, with the latter being favored albeit non-significantly. In the $k=2$ case, the Partitioning Around Medoids (PAM) clustering algorithm recovers a metal-poor halo GC group and a metal-rich disk GC group. With $k=3$ the three groups can be interpreted as disk/inner-halo/outer-halo families. For each group we obtain a medoid, i.e. a representative element (NGC $6352$, NGC $5986$, and NGC $5466$ for the disk, inner halo, and outer halo respectively), and a measure of how strongly each GC is associated to its group, the so-called silhouette width. Using the latter, we find a correlation with age for both disk and outer halo GCs where the stronger the association of a GC with the disk (outer halo) group, the younger (older) it is.

Searching for distinctive signatures, which characterize different formation channels of binary black holes (BBHs), is a crucial step towards the interpretation of current and future gravitational wave detections. Here, we investigate the demography of merging BBHs in young star clusters (SCs), which are the nursery of massive stars. We performed $4\times{} 10^3$ N-body simulations of SCs with metallicity $Z=0.002$, initial binary fraction $0.4$ and fractal initial conditions, to mimic the clumpiness of star forming regions. Our simulations include a novel population-synthesis approach based on the code MOBSE. We find that SC dynamics does not affect the merger rate significantly, but leaves a strong fingerprint on the properties of merging BBHs. More than 50 % of merging BBHs in young SCs form by dynamical exchanges in the first few Myr. Dynamically formed merging BBHs are significantly heavier than merging BBHs in isolated binaries: merging BBHs with total mass up to $\sim{}120$ M$_\odot$ form in young SCs, while the maximum total mass of merging BBHs in isolated binaries with the same metallicity is only $\sim{}70$ M$_\odot$. Merging BBHs born via dynamical exchanges tend to have smaller mass ratios than BBHs in isolated binaries. Furthermore, SC dynamics speeds up the merger: the delay time between star formation and coalescence is significantly shorter in young SCs. In our simulations, massive systems such as GW170729 form only via dynamical exchanges. Finally $\sim{}2$ % of merging BBHs in young SCs have mass in the pair-instability mass gap ($\sim{}60-120$ M$_\odot$). This represents a unique fingerprint of merging BBHs in SCs.

The shape of the radial distribution of Blue Straggler Stars (BSS), when normalized to a reference population of Horizontal Branch (HB) stars, has been found to be a powerful indicator of the dynamical evolution reached by a Globular Cluster (GC). In particular, observations suggest that the BSS distribution bimodality is modulated by the dynamical age of the host GC, with dynamically unrelaxed GCs showing a flat BSS distribution, and more relaxed GCs showing a minimum at a radius that increases for increasing dynamical age, resulting in a natural dynamical clock. While direct N-body simulations are able to reproduce the general trend, thus supporting its dynamical origin, the migration of the minimum of the distribution appears to be noisy and not well defined. Here we show that a simple unidimensional model based on dynamical friction (drift) and Brownian motion (diffusion) correctly reproduces the qualitative motion of the minimum, without adjustable parameters except for the BSS to HB stars mass-ratio. Differential dynamical friction effects combine with diffusion in creating a bimodality in the BSS distribution and determining its evolution, driving the migration of the minimum to larger radii over time. The diffusion coefficient is strongly constrained by the need to reproduce the migratory behaviour of the minimum, and the radial dependence of diffusion set by fundamental physical arguments automatically satisfies this constraint. Therefore, our model appears to capture the fluctuation-dissipation dynamics that underpins the dynamical clock.

The high velocity gradient observed in the compact cloud CO-0.40-0.22, at a projected distance of 60 pc from the centre of the Milky Way, has led its discoverers to identify the closeby mm continuum emitter, CO-0.40-0.22*, with an intermediate mass black hole (IMBH) candidate. We describe the interaction between CO-0.40-0.22 and the IMBH, by means of a simple analytical model and of hydrodynamical simulations. Through such calculation, we obtain a lower limit to the mass of CO-0.40-0.22* of few $10^4 \times \; M_{\odot}$. This result tends to exclude the formation of such massive black hole in the proximity of the Galactic Centre. On the other hand, CO-0.40-0.22* might have been brought to such distances in cosmological timescales, if it was born in a dark matter halo or globular cluster around the Milky Way.

e-ASTROGAM (enhanced ASTROGAM) is a breakthrough Observatory space mission, with a detector composed by a Silicon tracker, a calorimeter, and an anticoincidence system, dedicated to the study of the non-thermal Universe in the photon energy range from 0.3 MeV to 3 GeV - the lower energy limit can be pushed to energies as low as 150 keV for the tracker, and to 30 keV for calorimetric detection. The mission is based on an advanced space-proven detector technology, with unprecedented sensitivity, angular and energy resolution, combined with polarimetric capability. Thanks to its performance in the MeV-GeV domain, substantially improving its predecessors, e-ASTROGAM will open a new window on the non-thermal Universe, making pioneering observations of the most powerful Galactic and extragalactic sources, elucidating the nature of their relativistic outflows and their effects on the surroundings. With a line sensitivity in the MeV energy range one to two orders of magnitude better than previous generation instruments, e-ASTROGAM will determine the origin of key isotopes fundamental for the understanding of supernova explosion and the chemical evolution of our Galaxy. The mission will provide unique data of significant interest to a broad astronomical community, complementary to powerful observatories such as LIGO-Virgo-GEO600-KAGRA, SKA, ALMA, E-ELT, TMT, LSST, JWST, Athena, CTA, IceCube, KM3NeT, and LISA.

Most globular clusters (GCs) are now known to host multiple stellar populations with different light element abundances. Here we use narrow-band photometry and low-resolution spectroscopy for NGC 362 and NGC 6723 to investigate their chemical properties and radial distributions of subpopulations. We confirm that NGC 362 and NGC 6723 are among the GCs with multiple populations showing bimodal CN distribution and CN-CH anti-correlation without a significant spread in calcium abundance. These two GCs show more centrally concentrated CN-weak earlier generation stars compared to the later generation CN-strong stars. These trends are reversed with respect to those found in previous studies for many other GCs. Our findings, therefore, seem contradictory to the current scenario for the formation of multiple stellar populations, but mass segregation acting on the two subpopulations might be a possible solution to explain this reversed radial trend.

Mergers of stellar-mass black holes were recently observed in the gravitational wave window opened by LIGO. This puts the spotlight on dense stellar systems and their ability to create intermediate-mass black holes (IMBHs) through repeated merging. Unfortunately, attempts at direct and indirect IMBH detection in star clusters in the nearby universe have proven inconclusive as of now. Indirect detection methods attempt to constrain IMBHs through their effect on star cluster photometric and kinematic observables. They are usually based on looking for a specific, physically motivated signature. While this approach is justified, it may be suboptimal in its usage of the available data. Here I present a new indirect detection method, based on machine learning, that is unaffected by these restrictions. I reduce the scientific question whether a star cluster hosts an IMBH to a classification problem in the machine learning framework. I present preliminary results to illustrate how machine learning models are trained on simulated datasets and measure their performance on previously unseen, simulated data.

Recently, both stellar mass-segregation and binary-fractions were uniformly measured on relatively large samples of Galactic Globular Clusters (GCs). Simulations show that both sizeable binary-star populations and Intermediate-Mass Black Holes (IMBHs) quench mass-segregation in relaxed GCs. Thus mass-segregation in GCs with a reliable binary-fraction measurement is a valuable probe to constrain IMBHs. In this paper we combine mass-segregation and binary-fraction measurements from the literature to build a sample of 33 GCs (with measured core-binary fractions), and a sample of 43 GCs (with a binary fraction measurement in the area between the core radius and the half-mass radius). Within both samples we try to identify IMBH-host candidates. These should have relatively low mass-segregation, a low binary fraction (< 5%), and short (< 1 Gyr) relaxation time. Considering the core binary fraction sample, no suitable candidates emerge. If the binary fraction between the core and the half-mass radius is considered, two candidates are found, but this is likely due to statistical fluctuations. We also consider a larger sample of 54 GCs where we obtained an estimate of the core binary fraction using a predictive relation based on metallicity and integrated absolute magnitude. Also in this case no suitable candidates are found. Finally, we consider the GC core- to half-mass radius ratio, that is expected to be larger for GCs containing either an IMBH or binaries. We find that GCs with large core- to half-mass radius ratios are less mass-segregated (and show a larger binary fraction), confirming the theoretical expectation that the energy sources responsible for the large core are also quenching mass-segregation

Context. Machine-Learning (ML) solves problems by learning patterns from data, with limited or no human guidance. In Astronomy, it is mainly applied to large observational datasets, e.g. for morphological galaxy classification. Aims. We apply ML to gravitational N-body simulations of star clusters that are either formed by merging two progenitors or evolved in isolation, planning to later identify Globular Clusters (GCs) that may have a history of merging from observational data. Methods. We create mock-observations from simulated GCs, from which we measure a set of parameters (also called features in the machine-learning field). After dimensionality reduction on the feature space, the resulting datapoints are fed to various classification algorithms. Using repeated random subsampling validation we check whether the groups identified by the algorithms correspond to the underlying physical distinction between mergers and monolithically evolved simulations. Results. The three algorithms we considered (C5.0 trees, k-nearest neighbour, and support-vector machines) all achieve a test misclassification rate of about 10% without parameter tuning, with support-vector machines slightly outperforming the others. The first principal component of feature space correlates with cluster concentration. If we exclude it from the regression, the performance of the algorithms is only slightly reduced.

The formation of the Milky Way stellar halo is thought to be the result of merging and accretion of building blocks such as dwarf galaxies and massive globular clusters. Recently, Deason et al. (2015) suggested that the Milky Way outer halo formed mostly from big building blocks, such as dwarf spheroidal galaxies, based on the similar number ratio of blue straggler (BS) stars to blue horizontal-branch (BHB) stars. Here we demonstrate, however, that this result is seriously biased by not taking into detailed consideration on the formation mechanism of BHB stars from helium enhanced second-generation population. In particular, the high BS-to-BHB ratio observed in the outer halo fields is most likely due to a small number of BHB stars provided by GCs rather than to a large number of BS stars. This is supported by our dynamical evolution model of GCs which shows preferential removal of first generation stars in GCs. Moreover, there are a sufficient number of outer halo GCs which show very high BS-to-BHB ratio. Therefore, the BS-to-BHB number ratio is not a good indicator to use in arguing that more massive dwarf galaxies are the main building blocks of the Milky Way outer halo. Several lines of evidence still suggest that GCs can contribute a significant fraction of the outer halo stars.

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.

Globular Cluster (GC) Color-Magnitude Diagrams (CMDs) are reasonably well understood in terms of standard stellar-evolution. However, there are still some open issues, such as fully accounting for the Horizontal Branch (HB) morphology in terms of chemical and dynamical parameters. Mass-loss on the Red Giant Branch (RGB) shapes the mass-distribution of the HB stars, and the color distribution in turn. The physical mechanisms driving mass-loss are still unclear, as direct observations fail to reveal a clear correlation between mass-loss rate and stellar properties. The horizontal-branch mass-distribution is further complicated by Helium-enhanced multiple stellar populations, because of differences in the evolving mass along the HB. We present a simple analytical mass-loss model, based on tidal stripping through Roche-Lobe OverFlow (RLOF) during stellar encounters. Our model naturally results in a non-gaussian mass-loss distribution, with high skewness, and contains only two free parameters. We fit it to the HB mass distribution of four Galactic GCs, as obtained from fitting the CMD with Zero Age HB (ZAHB) models. The best-fit model accurately reproduces the observed mass-distribution. If confirmed on a wider sample of GCs, our results would account for the effects of dynamics in RGB mass-loss processes and provide a physically motivated procedure for synthetic CMDs of GCs. Our physical modeling of mass-loss may result in the ability to disentangle the effects of dynamics and helium-enhanced multiple-populations on the HB morphology and is instrumental in making HB morphology a probe of the dynamical state of GCs, leading to an improved understanding of their evolution.

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.

Context. Stellar collision rates in globular clusters (GCs) do not appear to correlate with horizontal branch (HB) morphology, sug- gesting that dynamics does not play a role in the second-parameter problem. However, core densities and collision rates derived from surface-brightness may be significantly underestimated as the surface-brightness profile of GCs is not necessarily a good indicator of the dynamical state of GC cores. Core-collapse may go unnoticed if high central densities of dark remnants are present. Aims. We test whether GC HB morphology data supports a dynamical contribution to the so-called second-parameter effect. Methods. To remove first-parameter dependence we fitted the maximum effective temperature along the HB as a function of metal- licity in a sample of 54 Milky Way GCs. We plotted the residuals to the fit as a function of second-parameter candidates, namely dynamical age and total luminosity. We considered dynamical age (i.e. the ratio between age and half-light relaxation time) among possible second-parameters. We used a set of direct N-body simulations, including ones with dark remnants to illustrate how core density peaks, due to core collapse, in a dynamical-age range similar to that in which blue HBs are overabundant with respect to the metallicity expectation, especially for low-concentration initial conditions. Results. GC total luminosity shows nonlinear behavior compatible with the self-enrichment picture. However, the data are amenable to a different interpretation based on a dynamical origin of the second-parameter effect. Enhanced mass-stripping in the late red-giant- branch phase due to stellar interactions in collapsing cores is a viable candidate mechanism. In this picture, GCs with HBs bluer than expected based on metallicity are those undergoing core-collapse.

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.

We study the binary fraction of the globular cluster M10 (NGC 6254) as a function of radius from the cluster core to the outskirts, by means of a quan- titative analysis of the color distribution of stars relative to the fiducial main sequence. By taking advantage of two data-sets, acquired with the Advanced Camera for Survey and the Wide Field Planetary Camera 2 on board the Hubble Space Telescope, we have studied both the core and the external regions of the cluster. The binary fraction is found to decrease from 14% within the core, to 1.5% in a region between 1 and 2 half-mass radii from the cluster centre. Such a trend and the derived values are in agreement with previous results ob- tained in clusters of comparable total magnitude. The estimated binary fraction is sufficient to account for the suppression of mass segregation observed in M10, without any need to invoke the presence of an intermediate-mass black hole in its centre.

Black holes are fascinating objects. As a class of solutions to the Einstein equations they have been studied a great deal, yielding a wealth of theoretical results. But do they really exist? What do astronomers really mean when they claim to have observational evidence of their existence? To answer these questions, I will focus on a particular range of black-hole masses, approximately from 100 to 10000 solar masses. Black holes of this size are named Intermediate Mass Black Holes (IMBHs) and their existence is still heavily disputed, so they will be perfect for illustrating the observational challenges faced by a black hole hunter

Studying the radial variation of the stellar mass function in globular clusters (GCs) has proved a valuable tool to explore the collisional dynamics leading to mass segregation and core collapse. In order to study the radial dependence of the luminosity and mass function of M 10, we used ACS/HST deep high resolution archival images, reaching out to approximately the cluster's half-mass radius (rhm), combined with deep WFPC2 images that extend our radial coverage to more than 2 rhm. From our photometry, we derived a radial mass segregation profile and a global mass function that we compared with those of simulated clusters containing different energy sources (namely hard binaries and/or an IMBH) able to halt core collapse and to quench mass segregation. A set of direct N-body simulations of GCs, with and without an IMBH of mass 1% of the total cluster mass, comprising different initial mass functions (IMFs) and primordial binary fractions, was used to predict the observed mass segregation profile and mass function. The mass segregation profile of M 10 is not compatible with cluster models without either an IMBH or primordial binaries, as a source of energy appears to be moderately quenching mass segregation in the cluster. Unfortunately, the present observational uncertainty on the binary fraction in M10 does not allow us to confirm the presence of an IMBH in the cluster, since an IMBH, a dynamically non-negligible binary fraction (~ 5%), or both can equally well explain the radial dependence of the cluster mass function.

In a previous paper we found that the Globular Clusters of our Galaxy lie around a line in the log(Re), SBe, log(sigma) parameter space, with a moderate degree of scatter and remarkable axi-symmetry. This implies the existence of a purely photometric scaling law obtained by projecting such a line onto the log(Re), SBe plane. Such photometric quantities are readily available for large samples of clusters, as opposed to stellar velocity dispersion data. We study a sample of 129 Galactic and extragalactic clusters on such photometric plane in the V-band. We look for a linear relation between SBe and log(Re) and study how the scatter around it is influenced by age and dynamical environment. We interpret our results as a test on the evolutionary versus primordial origin of the Fundamental Line. We perform a detailed analysis of surface brightness profiles, which allows us to present a catalogue of structural properties, without relying on a given dynamical model. We find a linear relation between SBe and log(Re), in the form SBe = (5.25 +- 0.44) log(Re) + (15.58 +- 0.28), where SBe is measured in mag/arcsec^2 and Re in parsec. Both young and old clusters lie on the scaling law, with a scatter of approximately 1 mag in SBe. However, young clusters display more scatter and a clear trend of such scatter with age, which old clusters do not. Such trend becomes tighter if cluster age is measured in units of the cluster half-light relaxation time. Two-body relaxation therefore plays a major role, together with passive stellar population evolution, in shaping the relation between SBe, log(Re), and cluster age. We argue that the log(Re)-SBe relation and hence the Fundamental Line scaling law is not primordially set at cluster formation, but rather is the result of combined stellar evolution and collisional dynamical evolution.

[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.

[abridged] Theoretical investigations have suggested the presence of Intermediate Mass Black Holes (IMBHs, with masses in the 100-10000 Msun range) in the cores of some Globular Clusters (GCs). In this paper we present the first application of a new technique to determine the presence or absence of a central IMBH in globular clusters that have reached energy equipartition via two-body relaxation. The method is based on the measurement of the radial profile for the average mass of stars in the system, using the fact that a quenching of mass segregation is expected when an IMBH is present. Here we measure the radial profile of mass segregation using main-sequence stars for the globular cluster NGC 2298 from resolved source photometry based on HST-ACS data. The observations are compared to expectations from direct N-body simulations of the dynamics of star clusters with and without an IMBH. The mass segregation profile for NGC 2298 is quantitatively matched to that inferred from simulations without a central massive object over all the radial range probed by the observations, that is from the center to about two half-mass radii. Profiles from simulations containing an IMBH more massive than ~ 300-500 Msun (depending on the assumed total mass of NGC 2298) are instead inconsistent with the data at about 3 sigma confidence, irrespective of the IMF and binary fraction chosen for these runs. While providing a null result in the quest of detecting a central black hole in globular clusters, the data-model comparison carried out here demonstrates the feasibility of the method which can also be applied to other globular clusters with resolved photometry in their cores.

The globular clusters (GC) of our Galaxy have been found to lie close to a plane in the log(R_e), log(sigma), SB_e space, on the continuation of the Fundamental Plane (FP) known to characterize the properties of early-type galaxies. We reexamine the issue on a sample of 48 GCs selected in terms of homogeneity criteria for the photometric data available from the literature and perform a model-independent analysis of surface brightness profiles (SBP) and distance moduli, estimating error bars and studying selection effects with non-parametric statistical tests. We determine the coefficients of the FP and their error bars. The scatter from the FP relation is likely to be intrinsic, i.e. not due to measurement errors only. We find that in the standard FP coordinates our sample occupies a slim, axisymmetric region of parameter space, suggesting that the scaling relation might be around a Fundamental Line, rather than a plane, confirming a result noted earlier. This is likely to be the origin of the difficulties in the fit by a plane mentioned in previous investigations. Such FL relation would imply a pure photometric scaling law, which might be tested on wider samples and on extra-galactic GC systems. We find a correlation of the residuals from the FP relation with the central slope of the SBP. No other correlations are found. Finally, we reconstruct the distribution of the values of the quantity log(K_V/(M/L)) (virial coefficient divided by the mass-to-light ratio) through kernel density estimation and find evidence for bimodality, which suggests that the galactic GC system may be composed of at least two dynamically different populations. Yet, these populations do not reflect the standard dichotomy between disk and halo clusters. (abridged).