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

We use the Chandrasekhar formalism and direct N-body simulations to study the effect of dynamical friction on a test object only slightly more massive than the field stars, orbiting a spherically symmetric background of particles with a mass spectrum. The main goal is to verify whether the dynamical friction time (t_DF) develops a non-monotonic radial-dependence that could explain the bimodality of the Blue Straggler radial distributions observed in globular clusters. In these systems, in fact, relaxation effects lead to a mass and velocity radial segregation of the different mass components, so that mass-spectrum effects on t_DF are expected to be dependent on radius. We find that, in spite of the presence of different masses, t_DF is always a monotonic function of radius, at all evolutionary times and independently of the initial concentration of the simulated cluster. This because the radial dependence of t_DF is largely dominated by the total mass density profile of the background stars (which is monotonically decreasing with radius). Hence, a progressive temporal erosion of the BSS population at larger and larger distances from the cluster center remains the simplest and the most likely explanation of the shape of the observed BSS radial distributions, as suggested in previous works. We also confirm the theoretical expectation that approximating a multi-mass globular cluster as made of (averaged) equal-mass stars can lead to significant overestimates of t_DF within the half-mass radius.