This dataset also contains MOCCA simulation data for the star cluster model simulated in https://arxiv.org/abs/2602.15953. 

This dataset contains 6 files from the MOCCA simulation:

1) 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 2497 lines in this file contain header descriptions for all columns.

and 5 snapshot files:

- snap1.18.dat (at 1180 Myr)
- snap1.19.dat (at 1190 Myr)
- snap1.20.dat (at 1200 Myr)
- snap1.21.dat (at 1210 Myr)
- snap1.22.dat (at 1220 Myr)

Column Descriptions the snapshot files are given below:

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.
