We have carried out a set of Monte Carlo simulations to study a number of fundamental aspects of the dynamical evolution of multiple stellar populations in globular clusters with different initial masses, fractions of second generation (2G) stars, and structural properties. Our simulations explore and elucidate: 1) the role of early and long-term dynamical processes and stellar escape in the evolution of the fraction of 2G stars and the link between the evolution of the fraction of 2G stars and various dynamical parameters; 2) the link between the fraction of 2G stars inside the cluster and in the population of escaping stars during a clusters dynamical evolution; 3) the dynamics of the spatial mixing of the first-generation (1G) and 2G stars and the details of the structural properties of the two populations as they evolve toward mixing; 4) the implications of the initial differences between the spatial distribution of 1G and 2G stars for the evolution of the anisotropy in the velocity distribution and the expected radial profile of the 1G and 2G anisotropy for clusters at different stages of their dynamical history; 5) the variation of the degree of energy equipartition of the 1G and the 2G populations as a function of the distance from the clusters centre and the clusters evolutionary phase.
Download