Data for the paper: Multiple Stellar Populations in MOCCA Models – Transient Spatial Over-Concentration of Pristine RGB Stars

This dataset contains MOCCA simulation data for a globular cluster model used in the publication:
"Multiple stellar populations in MOCCA globular cluster models: Transient spatial over-concentration of pristine red giant stars driven by strong dynamical encounters." Preprint for the paper is available at: https://arxiv.org/abs/2502.17162

The simulation was performed using the MOCCA code version described in Giersz et al. (2024):
https://arxiv.org/abs/2411.06421

Included Files

This dataset contains three files from the MOCCA simulation:

1. mocca.ini – The input configuration file used for the MOCCA simulation.
2. system.dat – A file containing detailed record of the cluster evolution, storing over 500 parameters related to:
   • Global cluster properties
   • Stellar and binary system properties
   • Technical parameters
   • Dynamical interactions
     The first 2450 lines in this file contain header descriptions for all columns.
     Please Note: Changes in the etot column during the first 100 Myr are due to the gradual accretion of mass onto gas particles in the initial model. This process leads to a time-dependent variation in total energy, which should be considered when analyzing energy evolution.
3. snap13.dat – The non-projected 13 Gyr snapshot of the simulation, containing 35 columns (described below).

Column Descriptions for snap13.dat (13 Gyr Snapshot)

1. im – Star/binary name (inames/inameb)
2. r – Position in the cluster (pc)
3. vr – Radial velocity (km/s)
4. vt – Tangential velocity (km/s)
5. u – Potential energy [MC units]
6. ikb – Binary type (0 = single, 1 = binary, 2 = semi-detached binary)
7. a – Semi-major axis of a binary (solar radii, R☉), 0.0 for single stars
8. e – Eccentricity of the binary, 0.0 for single stars
9. ik1 – SSE/BSE stellar type of the first component (see table below)
10. ik2 – SSE/BSE stellar type of the second component
11. sm1 – Mass of the first component (solar mass, M☉)
12. sm2 – Mass of the second component (solar mass, M☉), 0.0 for single stars
13. popId1 – Population ID of the first component (1 = pristine, 2 = enriched, >=100 for mixed pop)
14. popId2 – Population ID of the second component
15. timenr – Timestep (integer)
16. idd1 – Star/binary ID for the first component
17. idd2 – Star/binary ID for the second component
18. lum1 – Luminosity of the first component (bolometric solar luminosity, L☉)
19. lum2 – Luminosity of the second component (bolometric solar luminosity, L☉)
20. rad1 – Radius of the first component (solar radii, R☉)
21. rad2 – Radius of the second component (solar radii, R☉)
22. hist1 – Compact history of the first star (dynamical events, mergers, etc.)
23. hist2 – Compact history of the second star
24. spin1 – Spin of the first component (SSE/BSE)
25. spin2 – Spin of the second component (SSE/BSE)
26. mu1 – U-band magnitude for the first component
27. mu2 – U-band magnitude for the second component
28. mv1 – V-band magnitude for the first component
29. mv2 – V-band magnitude for the second component
30. mb1 – B-band magnitude for the first component
31. mb2 – B-band magnitude for the second component
32. mi1 – I-band magnitude for the first component
33. mi2 – I-band magnitude for the second component
34. mtr1 – Mass transfer rate of the first component (solar mass per year, M☉/yr)
35. mtr2 – Mass transfer rate of the second component (solar mass per year, M☉/yr)

The compact history columns (hist1, hist2) above encode the short history of the events with the star

• 0000000000000000001 - had dynamical collision (inter-coll)
• 0000000000000000010 - had binary merger (stellar evolution)
• 0000000000000000100 - had mass transfer event in binary evolution (increased the mass)
• 0000000000000001000 - had interaction (bin-bin, or bin-sin)
• 0000000000000010000 - had exchange in interaction
• 0000000000000100000 - went to a binary (inter-binform)
• 0000000000001000000 - went to single (dissolution in stellar evolution)
• 0000000000010000000 - went to single (dissolution in interaction)
• 0000000000100000000 - separation changed by 1% (in interaction)
• 0000000001000000000 - separation changed by 10% (in interaction)
• 0000000010000000000 - eccentricity changed by 1% (in interaction)
• 0000000100000000000 - eccentricity changed by 10% (in interaction)
• 0000001000000000000 - dissolution in bin-sin interaction
• 0000010000000000000 - dissolution in bin-bin interaction
• 0000100000000000000 - dissolution in binary evolution
• 0001000000000000000 - exchange in bin-sin interaction
• 0010000000000000000 - exchange in bin-bin interaction
• 0100000000000000000 - dynamical collision in bin-sin
• 1000000000000000000 - dynamical collision in bin-bin (or hierarchical)
  E.g. if the star had two dynamical collisions and 5 mass transfer events in binary evolution then it has hist column value of: 0000000052

The values in ik1 and ik2 correspond to SSE/BSE stellar types from Hurley et al. (2000; 2002):

0 - Fully convective low-mass MS star (M < 0.7 M☉)
1 - Main Sequence star (M > 0.7 M☉)
2 - Hertzsprung Gap (subgiant star)
3 - First Giant Branch
4 - Core Helium Burning
5 - First Asymptotic Giant Branch
6 - Second Asymptotic Giant Branch
7 - Main Sequence Naked Helium star
8 - Hertzsprung Gap Naked Helium star
9 - Giant Branch Naked Helium star
10 - Helium White Dwarf
11 - Carbon/Oxygen White Dwarf
12 - Oxygen/Neon White Dwarf
13 - Neutron Star
14 - Black Hole

If you have any questions about the data or require more extensive data then please feel free to contact us.