Instrumentation and Methods for Astrophysics 5 JAN, 2015

MOCCA code for star cluster simulations - I. Blue Stragglers, first results

MOCCA code for star cluster simulations - I. Blue Stragglers, first results

We introduce an improved code for simulations of star clusters, called MOCCA. It combines the Monte Carlo method for star cluster evolution and the Fewbody code to perform scattering experiments. The Fewbody was added in order to track more precisely dynamical interactions between objects which can lead to creations of various exotic objects observed in the star clusters, like Blue Stragglers Stars (BSS). The MOCCA code is currently one of the most advanced codes for simulating real size star clusters. It follows the star cluster evolution closely to N-body codes but is much faster. We show that the MOCCA code is able to follow the evolution of BSS with details. It is a suitable tool to perform full scale evolution of real star clusters and detail comparison with observations of exotic star cluster objects like BSS. This paper is the first one of the series of papers about properties of BSS in star clusters. This type of stars is particularly interesting today, because by studying them one can get important constrains on a link between the stellar and dynamical evolution of star clusters. We discuss here first results concerning BSS for an arbitrary chosen test model. We investigate properties of BSS which characterize different channels of formation like masses, semi-major axes, eccentricities, and orbital periods. We show how BSS from different channels change their types, and discuss initial and final positions of BSS, their bimodal distribution in the star cluster, lifetimes and more.

💡 Research Summary

The paper presents MOCCA, a state‑of‑the‑art simulation framework that couples a Monte‑Carlo star‑cluster evolution engine with the Fewbody direct‑integration code for small‑N dynamical encounters. By embedding Fewbody, MOCCA can resolve binary–binary, binary–single, and three‑body interactions on the fly, thereby achieving N‑body‑level fidelity while retaining the speed advantage of Monte‑Carlo methods. The authors first validate the code against direct N‑body simulations, showing excellent agreement in global cluster properties such as mass loss, core contraction, half‑mass radius evolution, and escape rates.

The main scientific focus is the formation and evolution of Blue Straggler Stars (BSS). Two principal formation channels are examined: (1) mass‑transfer (MT) in primordial or dynamically formed binaries, and (2) direct collisions (COL) arising from close encounters or resonant three‑body interactions. MOCCA tracks each BSS from its birth time, recording the parent objects, orbital parameters (semi‑major axis, eccentricity, period), and spatial location within the cluster.

Key findings include:

  • Channel‑specific orbital characteristics – MT‑formed BSS typically originate in binaries with low eccentricities (e < 0.2) and wide separations (several thousand AU), often located at intermediate radii (0.5–1.0 r_h). COL‑formed BSS arise deep in the core (r < 0.1 r_h) from highly eccentric (e > 0.5) encounters, and they inherit compact orbits immediately after formation.

  • Channel transition – A non‑negligible fraction of BSS experience a “channel switch.” For example, a star that first gains mass via MT can later be involved in a resonant three‑body encounter that leads to a collision, effectively converting an MT‑type BSS into a COL‑type BSS. This demonstrates the fluid nature of BSS evolution in dense stellar environments.

  • Spatial distribution and bimodality – The simulated cluster exhibits a clear bimodal radial distribution of BSS: a central peak driven by collision‑induced BSS and an outer peak populated by MT‑induced BSS. This mirrors observational studies of Galactic globular clusters, supporting the idea that the two peaks trace distinct dynamical regimes.

  • Lifetimes – Collision‑formed BSS have relatively short lifetimes (~0.5 Gyr) due to rapid subsequent interactions and higher mass‑loss rates, whereas MT‑formed BSS survive longer (~2 Gyr), reflecting the more stable nature of mass‑transfer binaries.

  • Migration – After formation, most BSS migrate toward either the core or the outer peak within 1–2 Gyr, driven by mass segregation and dynamical friction. The authors quantify the typical radial drift and show that the final positions correlate strongly with the formation channel.

Overall, the study demonstrates that MOCCA can reproduce detailed BSS demographics—including masses, orbital elements, radial distributions, and evolutionary pathways—with a fidelity comparable to direct N‑body models but at a fraction of the computational cost. The authors argue that MOCCA is therefore an ideal platform for systematic investigations of exotic stellar populations (e.g., BSS, X‑ray binaries, stellar‑mass black holes) across a wide range of initial cluster conditions.

Future papers in the series will explore how varying initial binary fractions, metallicities, and cluster concentration parameters affect BSS statistics, aiming to place tighter constraints on the interplay between stellar evolution and cluster dynamics.