2022MNRAS.517.4768H

We present an upgraded version of the MOCCA code for the study of dynamical evolution of globular clusters (GCs) and its first application to the study of evolution of multiple stellar populations. We explore initial conditions spanning different structural parameters for the first (FG) and second generation of stars (SG) and we analyse their effect on the binary dynamics and survival. Here, we focus on the number ratio of FG and SG binaries, their spatial variation, and the way their abundances are affected by various cluster initial properties. We find that present-day SG stars are more abundant in clusters that were initially tidally filling. Conversely, FG stars stay more abundant in clusters that were initially tidally underfilling. We find that the ratio between binary fractions is not affected by the way we calculate these fractions [e.g. only main-sequence binaries (MS) or observational binaries, i.e. MS stars >0.4 M_⊙ mass ratios >0.5]. This implies that the MS stars themselves are a very good proxy for probing entire populations of FG and SG. We also discuss how it relates to the observations of Milky Way GCs. We show that MOCCA models are able to reproduce the observed range of SG fractions for Milky Way GCs for which we know these fractions. We show how the SG fractions depend on the initial conditions and provide some constraints for the initial conditions to have more numerous FG or SG stars at the Hubble time.