Issue 
A&A
Volume 638, June 2020



Article Number  A155  
Number of page(s)  11  
Section  Galactic structure, stellar clusters and populations  
DOI  https://doi.org/10.1051/00046361/202037490  
Published online  29 June 2020 
Primordial mass segregation of star clusters with primordial binaries
^{1}
Astronomical Institute of the Czech Academy of Sciences, Boční II 1401, 141 31 Prague 4, Czech Republic
email: pavlik@asu.cas.cz
^{2}
Astronomical Institute of Charles University, V Holešovičkách 2, 180 00 Prague 8, Czech Republic
^{3}
Observatory and Planetarium of Prague, Královská obora 233, 170 21 Prague 7, Czech Republic
Received:
13
January
2020
Accepted:
29
April
2020
Context. Observations of young starforming regions suggest that star clusters are born completely mass segregated. These initial conditions are, however, gradually lost as the star cluster evolves dynamically. For star clusters with single stars only and a canonical initial mass function, it has been suggested that traces of these initial conditions vanish at a time τ_{v} between 3 and 3.5 t_{rh} (initial halfmass relaxation times).
Aims. Since a significant fraction of stars are observed in binary systems and it is widely accepted that most stars are born in binary systems, we aim to investigate what role a primordial binary population (even up to 100% binaries) plays in the loss of primordial mass segregation of young star clusters.
Methods. We used numerical Nbody models similar in size to the Orion Nebula Cluster (ONC) – a representative of young open clusters – integrated over several relaxation times to draw conclusions on the evolution of its mass segregation. We also compared our models to the observed ONC.
Results. We found that τ_{v} depends on the binary star fraction and the distribution of initial binary parameters that include a semimajor axis, eccentricity, and mass ratio. For instance, in the models with 50% binaries, we find τ_{v} = (2.7 ± 0.8) t_{rh}, while for 100% binary fraction, we find a lower value τ_{v} = (2.1 ± 0.6) t_{rh}. We also conclude that the initially completely mass segregated clusters, even with binaries, are more compatible with the presentday ONC than the nonsegregated ones.
Key words: methods: numerical / methods: data analysis / open clusters and associations: individual: ONC / stars: formation / binaries: general
© ESO 2020
1. Introduction
As a star cluster dynamically evolves, mass segregation establishes, and we observe a general tendency of clusters to evolve towards an even higher degree of mass segregation (e.g. Chandrasekhar 1943; Chandrasekhar & von Neumann 1942, 1943). Recent observations of the Serpens South starforming region performed with ALMA by Plunkett et al. (2018), however, suggest that young clusters are already born completely mass segregated. In Pavlík et al. (2019a, hereafter Paper I), we investigated this primordial mass segregation and were the first to point out that it is unexpectedly gradually lost because of twobody encounters that lead to energy equipartition. We then concluded that the degree of mass segregation of a nonsegregated and a completely segregated system should gradually settle at a similar level. Hence, both initial conditions become observationally indistinguishable (i.e. the primordial difference vanished) after some time designated as τ_{v}. In Nbody models with single stars only, we estimated this time to be 3 < τ_{v}/t_{rh} < 3.5 (where t_{rh} stands for the halfmass twobody relaxation time; cf. Spitzer & Hart 1971).
It is argued, for example, by Kroupa (1995a) and Goodwin & Kroupa (2005) that most (perhaps all) stars are preferentially born in pairs. Observations report 57% of Gdwarfs (Duquennoy & Mayor 1991; Raghavan et al. 2010), 45% of Kdwarfs (Mayor et al. 1992), and 42% of field, that is old, Mdwarfs (Fischer & Marcy 1992) are present in binaries. And the binary fraction increases with the stellar mass. This dependency on primary mass results if all stars form as binaries, which are dynamically processed in embedded clusters of various masses (Marks & Kroupa 2011; Thies et al. 2015, and their Fig. 6). Hence, we extended the work presented in Paper II by studying the evolution of mass segregation in star clusters that include primordial binaries in Pavlík (2020, hereafter Paper II), using several realisations of a model with 50% binaries. In this work, we go even further with higher number of realisations for better statistical description, different binary star distributions, and also binary fractions up to 100%.
2. Models
We prepared an ensemble of isolated Nbody models with 2.4k stars in the mass range from 0.08 to 50 M_{⊙} (comparable to the sources in the Orion Nebula Cluster, ONC, Pavlík et al. 2019b) with the optimally sampled canonical IMF (Kroupa 2001; Kroupa et al. 2013), where the highest stellar mass is determined from the maximum stellarmasstoclustermass relation (Weidner & Kroupa 2006; PflammAltenburg & Kroupa 2007)^{1}. All clusters started as Plummer models (Plummer 1911; André et al. 2014), where the halfmass radius was determined from the birth radius to embedded cluster mass relation (Marks & Kroupa 2012). We evolved them for several relaxation times using the stateoftheart numerical integrator nbody6 (Aarseth 2003)^{2}.
Similarly to Paper I, we set up two extreme degrees of primordial mass segregation according to a method of Baumgardt et al. (2008) – none or complete, for instance, the mass segregation parameter is S = 0 or S = 1, respectively. The primordially masssegregated models and their parameters are labelled seg, the nonsegregated ones are labelled non. The models’ physical properties are summarised in Table 1.
Initial parameters of the star cluster models.
In this work, two initial binary populations were analysed: (i) a conservative 50% (i.e. 601 binary stars in total) and (ii) 100% stars in binaries. In the models with 50% binaries, two pairing methods were employed (cf. Küpper et al. 2011)^{3}: (a) random pairing across the whole mass range (labelled P:rnd in the following plots), and (b) pairing based on a uniform mass ratio distribution (0.1 < q < 1.0) for all stars above 5 M_{⊙} and random pairing of the remaining stars up to the desired percentage (labelled P:uni, as well as in Paper II). In the case of the models with 100% binaries, only the uniform pairing method (described above) was used. The semimajor axes were distributed according to Sana et al. (2012) and Oh et al. (2015) distributions for stars with m > 5 M_{⊙} (with the periods P in days such that log_{10}P ≥ 0.15) and according to Kroupa (1995a) for lower mass stars (with log_{10}P ≥ 1). Eccentricity distribution of highmass systems was taken from Sana et al. (2011) and was thermal for lowmass stars (cf. Heggie 1975; Duquennoy & Mayor 1991; Kroupa et al. 2008). The lowmass shortperiod binaries were initialised in two different ways: (i) in the models labelled eig, we adopted Kroupa (1995b) premain sequence eigenevolution, and (ii) in the models without such a label, we did not. Results of both of these initial conditions are also compared with respect to the cluster’s mass segregation. For reference, the initial binary population parameters (semimajor axis, eccentricity, and mass ratio) are plotted in Figs. 1 and 2 for all models.
Fig. 1. Initial semimajor axis and eccentricity of the binary populations used in this work, distinguished by the primary component mass – either m_{1} < 5 M_{⊙} (grey) or m_{1} ≥ 5 M_{⊙} (black). Twenty realisations from models with and without eigenevolution are used, while segregated and nonsegregated models are combined in these plots. 
Fig. 2. Initial distributions of the semimajor axis, eccentricity and mass ratio of the binary populations used in this work. Two mass groups are distinguished according to the primary component mass: m_{1} < 5 M_{⊙} (grey) or m_{1} ≥ 5 M_{⊙} (black). Models with (solid line) and without (dashed line) eigenevolution are compared. Models with and without primordial mass segregation (ten realisations each) are combined in each curve. 
Each model in Table 1 is simulated many times (see the bottom two rows) for each combination of the initial parameters using different random seeds to acquire good statistics. The models presented are still idealised in the sense that the stellar evolution is suppressed, the interstellar gas (and therefore gas expulsion) is not considered, and the star clusters are isolated. We discuss the effects of these features in Appendix A.
The initial comparison of different pairing methods with the same fraction of binaries (i.e. models P:rnd and P:uni) in the highmass range (primary component has m_{1} > 5 M_{⊙}; i.e. black lines in Fig. 2) reveals that the P:rnd model has no equalmass binary stars and instead produces ≈80% of systems (≈70% after eigenevolution) where the primary component is at least 20 times more massive than the secondary (q < 0.05). Generally, when primordial binaries are included in the models, they are able to speed up mass segregation in the nonsegregated clusters because they are viewed by the rest of the stars as more massive bodies. In this particular case, where the random pairing of highmass stars with very lowmass stars predominantly occurs, the dynamical impact of binaries with highmass primary components is not going to be much higher than if these components were just single stars.
Concerning the models with and without eigenevolution, there is almost no difference in the initial distribution of the highmass binaries. On the other hand, the lowermass binaries (where the primary component has m_{1} < 5 M_{⊙}) do change their initial properties after the premain sequence eigenevolution. (i) There is a higher number of close binary systems (compare the plots in the top and bottom rows of Fig. 1, and the plotted dashed and solid grey lines in the left column of Fig. 2). (ii) The eccentricities become more uniformly distributed between 0 and 1 (see the middle column of plots in Fig. 2). (iii) The mass ratios shift towards unity, thus, more equalmass binaries are created (see the spike at q = 1 in the right column of plots in Fig. 2). These binaries could potentially speed up the mass segregation in the nonsegregated models.
3. Results
3.1. Mass segregation
Our models with binaries (both 50% and 100%) evolve in a similar way to the single star models presented in Paper I. The difference in the seg and non clusters is very high initially. The primordially fully masssegregated clusters, however, start to lose their initial ordering gradually due to random twobody encounters between massive stars, especially in the core region, which can efficiently eject stars out of the cluster (cf. Oh et al. 2015; Oh & Kroupa 2016; Kroupa et al. 2018; Wang et al. 2019). At the same time, the clusters without initial mass segregation do quickly establish it dynamically. After a period of time, both models settle almost at the same level of mass segregation. As in Paper I, we investigated the evolution of mass segregation using the spatial distribution of mass. We read the mean mass comprised in concentric spheres around the centre of the cluster (the density centre according to a method of Casertano & Hut 1985, implemented in nbody6, is taken), as plotted in Fig. 3 for each model, and calculated an integral parameter
with the mean mass in the kth bin at radius r_{k}
where n_{i} and m_{i} are the number of stars, and their total mass in an ith bin, respectively. Δr_{k} is the width of the kth bin, thus the thickness of the spherical shell that is being added, and f_{bin} is the total number of bins. The bins here are logarithmically equidistant, and, in particular, we chose r_{1} = 0.1 pc, r_{fbin} = 10 pc and f_{bin} = 50; the lower boundary is approximately the core radius, while the upper boundary initially contains almost the whole cluster (see Fig. A.2) and would be equivalent to the tidal radius (see Appendix A and Fig. A.3).
Fig. 3. Evolution of radial distribution of the mean stellar mass in the cluster models with eigenevolution (red lines are for nonsegregated models and blue lines for segregated models). Each row of plots corresponds to one model (top and middle have 50% initial binaries but different pairing methods; the bottom row is with 100% initial binaries). Time increases from left to right in the units of the initial halfmass relaxation time. 
The time when the difference of the initial conditions stabilises is estimated in each model from the ratio A_{seg}/A_{non}. We fitted this ratio by a linear function over a moving window of various lengths from 0.5 to 1.5 t_{rh} and looked for the moment when the slope first became horizontal (due to the noisiness of the data, we allowed for the slope to incline up to ±10^{−4} from zero). The median value of these times from all the used windows is then used as τ_{v}. Without any constraints on the vertical position of the fit, the plateau around τ_{v} appears always near A_{seg}/A_{non} ≈ 1, which signifies that the primordial mass segregation truly vanished around that time, and these clusters would be observationally indistinguishable from each other. Since in all models (with or without eigenevolution) the slopes after ⟨τ_{v}⟩ are near flat or have only a shallow slope (for example, in comparison to the first relaxation time) the clusters seem to evolve synchronously. This may also be verified with Fig. 3, where the evolution of each cluster radial profile is plotted (only models with eigenevolution are shown, because the models without it are very similar, which can also be verified from the evolution plotted in Fig. 4).
Fig. 4. Evolution of the ratio given by Eq. (1) in time of the models used. The plots on the righthand side are with shortperiod binary eigenevolution. The black vertical dashed line and the value ⟨τ_{v}⟩ represent the mean time when the slope of the data points became flat, the thinner vertical lines are the calculated 1σ uncertainties. The horizontal fit in the corresponding time window is indicated by a grey dashed line. 
Systems with 50% binaries and random pairing (P:rnd) yield ⟨τ_{v}⟩≈2.7 t_{rh} with or without the premain sequence eigenevolution of the shortperiod binaries (see the top panels of Fig. 4). The core of the primordially masssegregated cluster is populated by binaries composed of stars with more than 5 M_{⊙} paired with 10 to 20 times less massive companions (for the mass ratios, see the top right plot in Fig. 2; for the initial core population, compare the blue curves at t = 0 t_{rh} in Fig. 3). These binaries, which are effectively “single" stars, are not as efficient in ejecting themselves out of the core as the equalmass highmass binaries would be, so they reside there longer. From the point of view of the primordially nonsegregated model, these “single” stars are scattered all over the cluster, but also mainly around the core region. As they are still among the most massive bodies, they tend to segregate first, and because of their inefficiency to eject themselves out of the cluster afterwards, they partially halt further mass segregation of lower mass binaries. For the exact time needed to eject the most massive stars and how it compares to the other models, see the curves in Fig. A.1; and also compare the curves in Fig. 3, for example, at t = 5 t_{rh} (it stays higher in the P:rnd model than in the P:uni, meaning that these massive stars are still present in the core). Therefore, even though the premain sequence eigenevolution produces an excess of equalmass lowmass binaries that are able to speed up the mass segregation of the nonsegregated clusters (as they do in the following P:uni model), their effect is suppressed in the P:rnd model. Consequently, we cannot see any difference in the evolution of mass segregation between the P:rnd models with and without eigenevolution.
The systems with 50% binaries and P:uni pairing seem to evolve from their initial mass segregation faster when the eigenevolution is accounted for. This could be due to the early ejection of highmass binaries (see, e.g. Fig. A.1) and a higher number of lowmass equalmass binaries that are present in the eigenevolved models. Such systems act in the same way as twice as massive single stars and, hence, they help the nonsegregated clusters segregate faster. Yet, within their 1σ uncertainties, the resulting times τ_{v} still seem to be comparable (see the middle panels of Fig. 4). We note that the value of τ_{v} obtained in Paper II for P:uni model with only a handful of realisations holds even here, where more realisations are analysed.
When comparing both P:rnd and P:uni models with or without eigenevolution, Fig. 4 shows a much milder and longer decrease of the ratio A_{seg}/A_{non} in the P:rnd model, although its fitted value of τ_{v} happens to be smaller (at least in the model without eigenevolution). The steeper decrease towards A_{seg}/A_{non} ≈ 1 in the P:uni model may be caused again by a faster dynamical evolution of the nonsegregated clusters due to an early dynamical influence of the equalmass lowmass binaries.
In the case of the clusters with 100% binaries, the time of vanishing of the initial mass segregation is consistent within the uncertainties with the models containing just 50% binaries. In particular ⟨τ_{v}⟩≈2.1 t_{rh} (with eigenevolution) and ⟨τ_{v}⟩≈2.7 t_{rh} (without eigenevolution), see the bottom panels of Fig. 4 and Table 2.
Resulting time τ_{v}/t_{rh}.
The values of τ_{v}, which are summarised in Tab. 2, are generally smaller in comparison to the singlestar models in Paper I (i.e. 3 < τ_{v}/t_{rh} < 3.5), but not significantly so, with the exception of the eig model containing 100% binaries. Two factors play a role here. While in the clusters without primordial binary population binaries start to form as late as during the core collapse (cf. Fujii & Portegies Zwart 2014; Pavlík & Šubr 2018), here we already injected them in our models at the beginning of the integration. Binary stars then lead to a faster dynamical evolution, as they act as more massive bodies. The opposing factor is that most of the binaries that are present are initially soft (especially in the models without eigenevolution) so they very quickly disrupt, leaving their components as single stars behind. Therefore, the time τ_{v} is expected to be smaller but not significantly so, for example, approximately by the time the single star models from Paper I needed to reach the core collapse (i.e. 0.5 t_{rh}), which is what our results indicate. We also note that the core collapse of the models presented here happens only a fraction of the relaxation time sooner (i.e. t_{cc} ≈ 0.25 t_{rh}, estimated from the minimum of the inner Lagrangian radius, see Fig. A.2).
To summarise, if we compare different pairing schemes for the same percentage of binaries (i.e. P:rnd and P:uni models), the deduced values of τ_{v} are almost the same. The biggest difference is visible in the models with 50% P:uni or 100% binaries where the value of τ_{v} is on average about 20% smaller if the premain sequence eigenevolution is accounted for. Higher numbers of close binary systems and more equalmass binaries are generated initially in the eigenevolved models (compare the top and bottom panels in Fig. 1, and the left and right column of plots in Fig. 2). Such binaries are harder to disrupt and could potentially speed up the mass segregation of the initially nonsegregated clusters, thus shortening the time τ_{v}. However, based on the uncertainties of the fitted value τ_{v}, this difference may not have any statistical importance. Studying models with an even higher fraction of close binaries with the mass ratio biased towards unity could resolve this.
3.2. The Orion Nebula Cluster
As in Paper I and Paper II, we also tested which initial conditions are more compatible with the observed ONC. The present day ONC possesses about 8.5% binaries (Reipurth et al. 2007), but according to the simulations of King et al. (2012), it could even have had 73% binaries initially.
From now on, we only use the models with premain sequence eigenevolution, as they give a more realistic population of binary stars (e.g. in the models without eigenevolution, the lowmass binaries with small semimajor axes are missing, and an artificial truncation is present). The radial distribution of stars in two mass bins (divided by 5 M_{⊙}) was calculated at different times and then compared to the same mass bins derived from the observational data in Pavlík et al. (2019b). At each time, we provided the pvalue of the Kolmogorov–Smirnov (KS) test. The distribution of stars of mass m > 5 M_{⊙} is the same in the models as in the observed ONC. Lowmass stars, however, have significant differences – taking the model’s raw data showed no similarity with the data. This is mainly due to interstellar extinction (Scandariato et al. 2011) and geometry of the cluster (elongation beyond 0.5 pc from the core along the gaseous filament in which the ONC resides; cf. Hillenbrand & Hartmann 1998). By modifying the data to account for both of these features (as in Paper I, Sect. 4), the primordially masssegregated models reached the KS test p > 0.05 between 1 and 3 Myr, which is equivalent to the current age of the ONC (2.5 Myr; Hillenbrand 1997; Palla & Stahler 1999). In the case of the initially nonsegregated models, none are compatible with the presentday ONC, not even with scaling and extinction included, although the model with 100% binaries was closely approaching the desired p value at 3 Myr. At three time frames, we plot the average radial distributions of stars in Fig. 5 for visual comparison. The time evolution of the pvalue of these modified models is plotted in Fig. 6, the unmodified models (seg and non), which yield zero or almost zero everywhere in the lowmass range, are not included. These results suggest that the ONC was likely completely masssegregated at birth.
Fig. 5. Comparison of average radial profiles of the projected models of segregated (blue) and nonsegregated clusters (red) with the real ONC (black). Only the models with premain sequence eigenevolution are plotted. Three time frames are shown, and two mass groups are compared in each model (m ≥ 5 M_{⊙} with a solid line and m < 5 M_{⊙} with a dashed line). The pvalues for these groups are plotted in Fig. 6. All three models are modified by artificial extinction and scaling of the lowmass stars along the yaxis, as discussed in Paper I. 
Fig. 6. Results of KS test between the ONC data and models with 50% binaries (top and middle panels) and with 100% binaries (bottom panel). Only models with scaling and extinction are shown, thus, the time is also in scaled units. 
4. Conclusions
In this paper, we extend the works of Pavlík et al. (2019a, Paper I) and Pavlík (2020, Paper II). We investigated the role of the primordial binary star population on mass segregation in star clusters of the size of the ONC using models with 50% and 100% initial binary stars.
In the models with 50% initial binaries, we see no difference in the initial pairing when studying the evolution of mass segregation. We may also conclude that even with binaries, the models evolve in a similar fashion to the models with only single stars, yet they do so somewhat faster, which could be due to a faster dynamical evolution caused by the presence of hard binaries straight from the beginning. The loss of primordial mass segregation was milder in the models where pairing was random for all stars (P:rnd) in comparison to the models with a uniform distribution of mass ratio for binaries with stars above 5 M_{⊙}, and random only in the lowmass range (P:uni). Yet when the models were initialised with a premain sequence eigenevolution of shortperiod binaries, this apparent difference in the initial evolution partially faded. The time when the primordially masssegregated and nonsegregated models became observationally indistinguishable (i.e. the difference between different initial conditions vanishes) is ⟨τ_{v}⟩ = 2.7 t_{rh} (P:rnd,eig) or ⟨τ_{v}⟩ = 2.6 t_{rh} (P:uni,eig) for the models with premain sequence binary eigenevolution, and ⟨τ_{v}⟩ = 2.7 t_{rh} (P:rnd) or ⟨τ_{v}⟩ = 3.2 t_{rh} (P:uni) for those without it. However, we are still within the 1σ uncertainties of the singlestar model from Paper I where the value was ⟨τ_{v}⟩ = 3.3 t_{rh}.
In the models with 100% binaries, we can also see that the initial difference of mass segregation disappeared. Again, the time needed was shorter than in the singlestar clusters (⟨τ_{v}⟩ = 2.1 t_{rh} with eigenevolution and ⟨τ_{v}⟩ = 2.7 t_{rh} without it).
All models with eigenevolution, which undoubtedly better represent the cluster binary population, also ended up better erasing the differences of primordial mass segregation than models without it. The fitted values of τ_{v} from all models are mutually within their 1σ uncertainties, thus, there may be no statistical significance between the evolution of mass segregation in the presented models. However, on average, we can see that the models with premain sequence eigenevolution tend to erase the initial differences in mass segregation faster, for instance, the time τ_{v} is approximately 20% smaller than in the models without eigenevolution. Such a behaviour would be expected, because the premain sequence eigenevolution produces more equalmass lowmass binaries that can speed up the mass segregation of the initially nonsegregated clusters once the heavier binaries eject themselves from the cluster core. We note that models containing larger numbers of stars, and higher abundances of hard and equalmass binaries (both highmass and lowmass) are needed to fully determine this. Since the distribution of the initial binary parameters could have a stronger influence on the final results, it would also be important to follow the constrains from observations of semimajor axes, mass ratios, or metallicity (see e.g. Moe & Di Stefano 2017; Belloni et al. 2017; Di Carlo et al. 2020).
We also compared our models to the presentday Orion Nebula Cluster. When we also account for the interstellar extinction and physical elongation of the cluster, our results suggest that the ONC was probably primordially mass segregated. Models without primordial mass segregation are only compatible with the ONC for the highmass stars, but never in the lowmass range. Further modelling including stellar evolution, elongation along the molecular cloud, and gas expulsion are, however, still needed to support this.
Public version from Jul 6, 2018. https://github.com/ahwkuepper/mcluster
Acknowledgments
This study was supported by Charles University (Grant SVV260441), by the Czech Academy of Sciences (project RVO:67985815) and by the Czech Science Foundation (project of Excellence No. 1820083S). Computational resources were provided by the CESNET LM2015042 and the CERIT Scientific Cloud LM2015085, under the programme “Projects of Large Research, Development, and Innovations Infrastructures”. The author also appreciates discussion with the participants of the “(Ms4)^{th} Aarseth Nbody Meeting” and especially the comments from Pavel Kroupa. The author feels indebted to an anonymous referee for their comments and suggestions.
References
 Aarseth, S. J. 2003, Gravitational NBody Simulations (Cambridge, UK: Cambridge University Press) [Google Scholar]
 André, P., Di Francesco, J., WardThompson, D., et al. 2014, Protostars Planets VI, 27, [Google Scholar]
 Baumgardt, H., & Sollima, S. 2017, MNRAS, 472, 744 [NASA ADS] [CrossRef] [Google Scholar]
 Baumgardt, H., De Marchi, G., & Kroupa, P. 2008, ApJ, 685, 247 [NASA ADS] [CrossRef] [Google Scholar]
 Belloni, D., Askar, A., Giersz, M., Kroupa, P., & RochaPinto, H. J. 2017, MNRAS, 471, 2812 [NASA ADS] [CrossRef] [Google Scholar]
 Binney, J., & Tremaine, S. 1994, Galactic Dynamics (Princeton, New Jersey: Princeton University Press) [Google Scholar]
 BlandHawthorn, J., & Gerhard, O. 2016, ARA&A, 54, 529 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Casertano, S., & Hut, P. 1985, ApJ, 298, 80 [Google Scholar]
 Chandrasekhar, S. 1943, ApJ, 97, 255 [NASA ADS] [CrossRef] [Google Scholar]
 Chandrasekhar, S., & von Neumann, J. 1942, ApJ, 95, 489 [Google Scholar]
 Chandrasekhar, S., & von Neumann, J. 1943, ApJ, 97, 1 [NASA ADS] [CrossRef] [Google Scholar]
 De Donder, E., & Vanbeveren, D. 2003, New A, 8, 817 [NASA ADS] [CrossRef] [Google Scholar]
 Di Carlo, U. N., Giacobbo, N., Mapelli, M., et al. 2019, MNRAS, 487, 2947 [NASA ADS] [CrossRef] [Google Scholar]
 Di Carlo, U. N., Mapelli, M., Giacobbo, N., et al. 2020, MNRAS, submitted [arXiv: 2004.09525] [Google Scholar]
 Duquennoy, A., & Mayor, M. 1991, A&A, 500, 337 [NASA ADS] [Google Scholar]
 Faber, S. M., & Gallagher, J. S. 1979, ARA&A, 17, 135 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Fischer, D. A., & Marcy, G. W. 1992, ApJ, 396, 178 [NASA ADS] [CrossRef] [Google Scholar]
 Fujii, M. S., & Portegies Zwart, S. 2014, MNRAS, 439, 1003 [NASA ADS] [CrossRef] [Google Scholar]
 Goodwin, S. P., & Kroupa, P. 2005, A&A, 439, 565 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Hansen, B. M. S., & Phinney, E. S. 1997, MNRAS, 291, 569 [NASA ADS] [Google Scholar]
 Heggie, D. C. 1975, MNRAS, 173, 729 [NASA ADS] [CrossRef] [Google Scholar]
 Hillenbrand, L. A. 1997, AJ, 113, 1733 [NASA ADS] [CrossRef] [Google Scholar]
 Hillenbrand, L. A., & Hartmann, L. W. 1998, ApJ, 492, 540 [NASA ADS] [CrossRef] [Google Scholar]
 Hurley, J. R., Pols, O. R., & Tout, C. A. 2000, MNRAS, 315, 543 [NASA ADS] [CrossRef] [Google Scholar]
 Hurley, J. R., Tout, C. A., & Pols, O. R. 2002, MNRAS, 329, 897 [NASA ADS] [CrossRef] [Google Scholar]
 Jonker, P. G., & Nelemans, G. 2004, MNRAS, 354, 355 [NASA ADS] [CrossRef] [Google Scholar]
 King, R. R., Parker, R. J., Patience, J., & Goodwin, S. P. 2012, MNRAS, 421, 2025 [NASA ADS] [CrossRef] [Google Scholar]
 Kroupa, P. 1995a, MNRAS, 277, 1491 [NASA ADS] [CrossRef] [Google Scholar]
 Kroupa, P. 1995b, MNRAS, 277, 1507 [NASA ADS] [Google Scholar]
 Kroupa, P. 2001, MNRAS, 322, 231 [NASA ADS] [CrossRef] [Google Scholar]
 Kroupa, P. 2008, in The Cambridge NBody Lectures, eds. S. J. Aarseth, C. A. Tout, & R. A. Mardling (Berlin Springer Verlag), Lect. Notes Phys., 760, 181 [NASA ADS] [CrossRef] [Google Scholar]
 Kroupa, P., Weidner, C., PflammAltenburg, J., et al. 2013, in The Stellar and SubStellar Initial Mass Function of Simple and Composite Populations, eds. T. D. Oswalt, G. Gilmore, et al., 115 [Google Scholar]
 Kroupa, P., Jeřábková, T., Dinnbier, F., Beccari, G., & Yan, Z. 2018, A&A, 612, A74 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Küpper, A. H. W., Maschberger, T., Kroupa, P., & Baumgardt, H. 2011, MNRAS, 417, 2300 [NASA ADS] [CrossRef] [Google Scholar]
 Marks, M., & Kroupa, P. 2011, MNRAS, 417, 1702 [NASA ADS] [CrossRef] [Google Scholar]
 Marks, M., & Kroupa, P. 2012, A&A, 543, A8 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mayor, M., Duquennoy, A., Halbwachs, J. L., & Mermilliod, J. C. 1992, in IAU Colloq. 135: Complementary Approaches to Double and Multiple Star Research, eds. H. A. McAlister, & W. I. Hartkopf, ASP Conf. Ser., 32, 73 [Google Scholar]
 Moe, M., & Di Stefano, R. 2017, ApJS, 230, 15 [Google Scholar]
 Oh, S., & Kroupa, P. 2016, A&A, 590, A107 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Oh, S., Kroupa, P., & PflammAltenburg, J. 2015, ApJ, 805, 92 [NASA ADS] [CrossRef] [Google Scholar]
 Palla, F., & Stahler, S. W. 1999, ApJ, 525, 772 [NASA ADS] [CrossRef] [Google Scholar]
 Pavlík, V. 2020, Contributions of the Astronomical Observatory Skalnate Pleso, 50, 456 [Google Scholar]
 Pavlík, V., & Šubr, L. 2018, A&A, 620, A70 [CrossRef] [EDP Sciences] [Google Scholar]
 Pavlík, V., Jeřábková, T., Kroupa, P., & Baumgardt, H. 2018, A&A, 617, A69 [CrossRef] [EDP Sciences] [Google Scholar]
 Pavlík, V., Kroupa, P., & Šubr, L. 2019a, A&A, 626, A79 [CrossRef] [EDP Sciences] [Google Scholar]
 Pavlík, V., Kroupa, P., & Šubr, L. 2019b, VizieR Online Data Catalog: J/A+A/626/A79 [Google Scholar]
 PflammAltenburg, J., & Kroupa, P. 2007, MNRAS, 375, 855 [NASA ADS] [CrossRef] [Google Scholar]
 Plummer, H. C. 1911, MNRAS, 71, 460 [CrossRef] [Google Scholar]
 Plunkett, A. L., FernándezLópez, M., Arce, H. G., et al. 2018, A&A, 615, A9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Podsiadlowski, P., Joss, P. C., & Hsu, J. J. L. 1992, ApJ, 391, 246 [NASA ADS] [CrossRef] [Google Scholar]
 Raghavan, D., McAlister, H. A., Henry, T. J., et al. 2010, ApJS, 190, 1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Reipurth, B., Guimarães, M. M., Connelley, M. S., & Bally, J. 2007, AJ, 134, 2272 [NASA ADS] [CrossRef] [Google Scholar]
 Sana, H., & Evans, C. J. 2011, in Active OB Stars: Structure, Evolution, Mass Loss, and Critical Limits, eds. C. Neiner, G. Wade, G. Meynet, & G. Peters, IAU Symp., 272, 474 [Google Scholar]
 Sana, H., de Mink, S. E., de Koter, A., et al. 2012, Science, 337, 444 [Google Scholar]
 Scandariato, G., Robberto, M., Pagano, I., & Hillenbrand, L. A. 2011, A&A, 533, A38 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Spitzer, L., Jr., & Hart, M. H. 1971, ApJ, 164, 399 [NASA ADS] [CrossRef] [Google Scholar]
 Šubr, L., Kroupa, P., & Baumgardt, H. 2012, ApJ, 757, 37 [NASA ADS] [CrossRef] [Google Scholar]
 Thies, I., PflammAltenburg, J., Kroupa, P., & Marks, M. 2015, ApJ, 800, 72 [NASA ADS] [CrossRef] [Google Scholar]
 Wang, L., Kroupa, P., & Jerabkova, T. 2019, MNRAS, 484, 1843 [NASA ADS] [CrossRef] [Google Scholar]
 Weidner, C., & Kroupa, P. 2006, MNRAS, 365, 1333 [NASA ADS] [CrossRef] [Google Scholar]
 Zapartas, E., de Mink, S. E., Izzard, R. G., et al. 2017, A&A, 601, A29 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
Appendix A: Discussion of the initial conditions
In this work, we simulated idealised star clusters, which means that the clusters are isolated, without intracluster gas, and stars were considered point masses without stellar evolution, although an IMF was used. Here, we further discuss how these initial conditions might affect the results presented.
Using a pointmass approximation in star clusters (i.e. within the framework of the Nbody problem) is reasonable since the stellar radii are much smaller than the separation of stars – tidal effects are therefore negligible. Very close random encounters, which could lead to stellar collisions did not occur in the simulations either. The only issue is that a few binaries were initialised as very tight ones (i.e. possible candidates for stellar collisions), but their components would not be able to merge as they are treated as point masses.
Stellar evolution relates to mass loss from stars, hence weakening of the potential well of the cluster and subsequent expansion of the system. Moreover, heavy stars evolving towards black holes (BHs) or neutron stars (NSs) could go through a supernova, receiving a nonnegligible natal kick (e.g. Hansen & Phinney 1997; Jonker & Nelemans 2004) and leave the cluster (cf. Baumgardt & Sollima 2017; Pavlík et al. 2018). Based on the IMF (Kroupa 2001) and the stellar evolutionary algorithm by Hurley et al. (2000), only five stars had the mass required to become a BH, and another ≈30 could become a NS. However, at the time ⟨τ_{v}⟩, which is close to 12 Myr, only the BHs would have been fully evolved, assuming subsolar metallicity (according to the estimate of Pavlík et al. 2018); using the solar metallicity, the evolution would take even longer. Thus, by the time ⟨τ_{v}⟩, the cluster could lose at most 0.1 M_{cl} if all of the BHs escaped. In our models, however, the most massive stars (and all five potential BHs) are put into binaries from the beginning. Those binaries survived for several Myr, which means that they would also evolve a significant part of their life according to binary stellar evolution algorithms (e.g. Hurley et al. 2002) including also the phase of common envelope – redistribution of mass may rejuvenate them and delay the supernova explosion to later times than, for example, the estimated 12 Myr (cf. Podsiadlowski et al. 1992; De Donder & Vanbeveren 2003; Zapartas et al. 2017; Pavlík et al. 2018; Di Carlo et al. 2019). In addition to that, we must also take into account the inner dynamical evolution. As it is visible in Fig. A.1, some of the top 10% most massive stars escape due to three or multibody encounters regardless of their stellar evolution. A similar result in the sense of binary BH ejection was also reported by Di Carlo et al. (2019). It appears, therefore, that in our models, including stellar evolution should not change the initial evolution or the time ⟨τ_{v}⟩. Nevertheless, simulations with full treatment for single and binarystellar evolution (possible with different evolutionary algorithms) would be necessary to fully determine this.
Fig. A.1. Average positions of the 10 most massive stars in each model. For reference, the halfmass radius is plotted with a dashed colour line, and the mean time ⟨τ_{v}⟩ of each model is plotted by a dashed vertical line. 
Modelling an isolated star cluster is also an idealisation. Due to the tidal field, a star cluster would gradually lose its stars through evaporation and shrink. An estimate of this tidal radius is
where M_{ct} is the mass of the cluster within this radius, r_{G} is the distance of the cluster from the centre of the Galaxy, and M_{G} is the mass of the Galaxy comprised in the radius r_{G} (e.g. Binney & Tremaine 1994).
In Fig. A.2, we plot the evolution of our isolated star clusters by the means of the Lagrangian radii. As expected, the cluster goes through the core collapse (at t_{cc} ≈ 0.25 t_{rh}; estimated from the minimum of the innermost radii) and then expands – a measure of the cluster expansion could be the halfmass radius. In Fig. A.3, we plot the evolution of the halfmass radius calculated from the stellar mass contained up to the cutoff at the tidal radius – we assume r_{G} = 5 kpc, which yields M_{G} ≈ 5 × 10^{10} M_{⊙} (Faber & Gallagher 1979; BlandHawthorn & Gerhard 2016). As the tidal radius and the upper boundary for evaluating mass segregation are the same throughout the integration (r_{t} ≈ r_{fbin} = 10 pc), it is unlikely that the resulting time ⟨τ_{v}⟩ would change if the clusters were placed in the background potential of the Galaxy.
Fig. A.2. Evolution of Lagrangian radii in presented models with eigenevolution. The Lagrangian radii are defined as the radii of concentric spheres containing fixed mass fractions (printed to the right of each curve), where the core is identified as the cluster’s density centre (cf. Casertano & Hut 1985; Aarseth 2003). The mean value ⟨τ_{v}⟩ of each model is plotted by a dashed vertical line for reference. 
Fig. A.3. Evolution of the halfmass radius (dashed colour lines) of each model calculated from stars that are bound to the cluster, i.e. up to the tidal radius from Eq. (A.1), which is plotted by a solid line. An artificial halfmass radius, which should represent the cluster evolution after gas expulsion, is plotted by the lighter dashed line for each model. The mean value ⟨τ_{v}⟩ of each model is plotted by a dashed vertical line for reference. 
The most important simplification that we used is that the clusters do not contain gas. Gas expulsion due to radiation pressure from massive stars is responsible for an early cluster expansion. In the case of the ONC, the cluster could inflate on average by a factor of r_{scl} ≈ 3 between approximately 0.5 and 1 Myr (e.g. Šubr et al. 2012). In Fig. A.3, we also plot an artificial halfmass radius represented by lighter colours. Up to 0.75 Myr, it is equal to the halfmass radius calculated from the cluster mass within the tidal radius, and from 0.75 Myr onward, it is scaled up three times to simulate gas expulsion. Since the separations of stars in the cluster become also about three times larger after gas expulsion, the evolution of the whole system slows down by a factor of , which is also included in the modified halfmass radius in Fig. A.3. All halfmass radii in Fig. A.3 (modified or not) reach a similar value at 5 Myr and then evolve almost identically. Nonetheless, we conjure that the time ⟨τ_{v}⟩ would shift to later times in models with gas expulsion due to their slower evolution. Models of embedded clusters with gas should, however, be calculated to clarify this.
All Tables
All Figures
Fig. 1. Initial semimajor axis and eccentricity of the binary populations used in this work, distinguished by the primary component mass – either m_{1} < 5 M_{⊙} (grey) or m_{1} ≥ 5 M_{⊙} (black). Twenty realisations from models with and without eigenevolution are used, while segregated and nonsegregated models are combined in these plots. 

In the text 
Fig. 2. Initial distributions of the semimajor axis, eccentricity and mass ratio of the binary populations used in this work. Two mass groups are distinguished according to the primary component mass: m_{1} < 5 M_{⊙} (grey) or m_{1} ≥ 5 M_{⊙} (black). Models with (solid line) and without (dashed line) eigenevolution are compared. Models with and without primordial mass segregation (ten realisations each) are combined in each curve. 

In the text 
Fig. 3. Evolution of radial distribution of the mean stellar mass in the cluster models with eigenevolution (red lines are for nonsegregated models and blue lines for segregated models). Each row of plots corresponds to one model (top and middle have 50% initial binaries but different pairing methods; the bottom row is with 100% initial binaries). Time increases from left to right in the units of the initial halfmass relaxation time. 

In the text 
Fig. 4. Evolution of the ratio given by Eq. (1) in time of the models used. The plots on the righthand side are with shortperiod binary eigenevolution. The black vertical dashed line and the value ⟨τ_{v}⟩ represent the mean time when the slope of the data points became flat, the thinner vertical lines are the calculated 1σ uncertainties. The horizontal fit in the corresponding time window is indicated by a grey dashed line. 

In the text 
Fig. 5. Comparison of average radial profiles of the projected models of segregated (blue) and nonsegregated clusters (red) with the real ONC (black). Only the models with premain sequence eigenevolution are plotted. Three time frames are shown, and two mass groups are compared in each model (m ≥ 5 M_{⊙} with a solid line and m < 5 M_{⊙} with a dashed line). The pvalues for these groups are plotted in Fig. 6. All three models are modified by artificial extinction and scaling of the lowmass stars along the yaxis, as discussed in Paper I. 

In the text 
Fig. 6. Results of KS test between the ONC data and models with 50% binaries (top and middle panels) and with 100% binaries (bottom panel). Only models with scaling and extinction are shown, thus, the time is also in scaled units. 

In the text 
Fig. A.1. Average positions of the 10 most massive stars in each model. For reference, the halfmass radius is plotted with a dashed colour line, and the mean time ⟨τ_{v}⟩ of each model is plotted by a dashed vertical line. 

In the text 
Fig. A.2. Evolution of Lagrangian radii in presented models with eigenevolution. The Lagrangian radii are defined as the radii of concentric spheres containing fixed mass fractions (printed to the right of each curve), where the core is identified as the cluster’s density centre (cf. Casertano & Hut 1985; Aarseth 2003). The mean value ⟨τ_{v}⟩ of each model is plotted by a dashed vertical line for reference. 

In the text 
Fig. A.3. Evolution of the halfmass radius (dashed colour lines) of each model calculated from stars that are bound to the cluster, i.e. up to the tidal radius from Eq. (A.1), which is plotted by a solid line. An artificial halfmass radius, which should represent the cluster evolution after gas expulsion, is plotted by the lighter dashed line for each model. The mean value ⟨τ_{v}⟩ of each model is plotted by a dashed vertical line for reference. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.