Exploring the dynamical evolution of binary stars in multiple-population globular clusters

The presence of multiple stellar populations in globular clusters leads to a complex dynamical environment that significantly influences the evolution of binary stars, which in turn impacts the evolution of the cluster itself. For this study, we used a series of Monte Carlo simulations run with the MOCCA code to investigate the long-term dynamical evolution of binary stars in globular clusters hosting two distinct stellar populations. We explored how global binary properties such as incidence, fraction, and spatial distribution evolve over time due to the unique dynamical environment associated with each population. Our results show how binaries in the more centrally concentrated second population (P2) experience increased rates of hardening and disruption relative to the first population (P1), leading to distinct radial profiles in binary incidence and fraction. We also demonstrate the difference in spatial mixing timescales for binaries compared to single stars, where binary stars in each population retain some memory of their initial configurations even after complete single star mixing. Additionally, we investigated the formation and evolution of mixed binaries (binaries composed of a P1 component and a P2 component), which form primarily within the core through dynamical interactions. Finally, we studied main sequence--white dwarf binaries and find that they represent a larger fraction of binaries in P1 compared to P2. The results of this paper highlight the interplay between cluster dynamics and the evolution of binary stars and how binaries can act as tracers of the cluster's initial conditions and dynamical evolution.