Issue 
A&A
Volume 639, July 2020



Article Number  A92  
Number of page(s)  16  
Section  Galactic structure, stellar clusters and populations  
DOI  https://doi.org/10.1051/00046361/202037843  
Published online  14 July 2020 
The effects of a background potential in star cluster evolution
A delay in the relaxation timescale and runaway collision processes
^{1}
Departamento de Astronomía, Facultad Ciencias Físicas y Matemáticas, Universidad de Concepción, Av. Esteban Iturra s/n Barrio Universitario, Casilla 160C, Concepción, Chile
email: breinoso@udec.cl
^{2}
Universität Heidelberg, Zentrum für Astronomie, Institut für Theoretische Astrophysik, AlbertUeberleStr. 2, 69120 Heidelberg, Germany
^{3}
Department of Astrophysics, American Museum of Natural History, New York, NY 10024, USA
^{4}
Universität Heidelberg, Interdisziplinäres Zentrum für Wissenschaftliches Rechnen, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
Received:
27
February
2020
Accepted:
22
April
2020
Runaway stellar collisions in dense star clusters are invoked to explain the presence of very massive stars or blue stragglers in the center of those systems. This process has also been explored for the first star clusters in the Universe and shown to yield stars that may collapse at some points into an intermediate mass black hole. Although the early evolution of star clusters requires the explicit modeling of the gas out of which the stars form, these calculations would be extremely timeconsuming and often the effects of the gas can be accurately treated by including a background potential to account for the extra gravitational force. We apply this approximation to model the early evolution of the first dense star clusters formed in the Universe by performing Nbody simulations, our goal is to understand how the additional gravitational force affects the growth of a very massive star through stellar mergers in the central parts of the star cluster. Our results show that the background potential increases the velocities of the stars, causing an overall delay in the evolution of the clusters and in the runaway growth of a massive star at the center. The population of binary stars is lower due to the increased kinetic energy of the stars, initially reducing the number of stellar collisions, and we show that relaxation processes are also affected. Despite these effects, the external potential enhances the mass of the merger product by a factor ∼2 if the collisions are maintained for long times.
Key words: galaxies: star clusters: general / galaxies: star formation / binaries: general
© ESO 2020
1. Introduction
The discovery of supermassive black holes (SMBHs) in the first billion years of the Universe (Fan et al. 2006; Mortlock et al. 2011; Wu et al. 2015; Reed et al. 2019) has led to the study of formation channels for massive black hole seeds early on (z ∼ 10−20), following the collapse of molecularcooling halos (∼10^{6} M_{⊙}) or atomiccooling halos (∼10^{8} M_{⊙}). The direct collapse of protogalactic clouds exposed to a moderate LymanWerner flux from neighboring star forming halos produced the most massive seeds (10^{3−5} M_{⊙}) (Wise et al. 2019), and recent numerical simulations suggest that contrary to what was previously thought, the UV flux is not that critical anymore in the case of the efficient merger of fragments (Suazo et al. 2019). Another promising mechanism to form massive seeds is the formation of very massive stars (VMS) in the centers of dense stellar systems through stellar mergers (Fujii & Portegies Zwart 2013; Katz et al. 2015; Sakurai et al. 2017; Reinoso et al. 2018), producing either VMSs with several 10^{2} M_{⊙} or intermediate mass black holes with ∼10^{3} M_{⊙} and potentially up to 10^{4} M_{⊙}. Such massive seeds could be present in the early Universe if accreting population III (Pop. III) protostars can merge before entering the main sequence or if the remaining gas in the cluster can be accreted by the central object (Boekholt et al. 2018). Furthermore, after the formation of a BH in the center of a stellar cluster, additional growth can be expected by tidal disruption events of stars passing close to the black hole (Sakurai et al. 2019; Bonetti et al. 2020).
Although most of these previous studies have focused on mergers in second generation star clusters, that is, clusters that formed from molecular clouds that were polluted by stellar winds or supernovae from the first generation of stars (Katz et al. 2015; Sakurai et al. 2017, 2019), we are mostly interested in the very first star clusters of the Universe given the particular properties of these systems and these stars. Fragmentation occurs at higher densities in primordial gas clouds (10^{9} cm^{−3} or higher) (e.g., Clark et al. 2011a; Greif et al. 2011, 2012; Smith et al. 2011, 2012; Latif et al. 2013a), with clusters potentially having halfmass radii of 0.1 pc, especially if dust grains are present in such clouds which trigger fragmentation at high densities (e.g., Omukai et al. 2005; Klessen et al. 2012; Bovino et al. 2016; Latif et al. 2016).
In those environments, the protostellar radii could also be enhanced due to the rapid accretion expected in primordial or lowmetallicity gas, which in turn are a consequence of much higher gas temperatures than at present day star formation, and so the protostellar radii are enhanced up to 300 R_{⊙} (Stahler et al. 1986; Omukai & Palla 2001, 2003). Moreover, high accretion rates of ∼0.1 M_{⊙} yr^{−1} have been reported in several simulations for these protostars (Hosokawa et al. 2012, 2013; Schleicher et al. 2013; Haemmerlé et al. 2018; Woods et al. 2017), implying ∼500 R_{⊙} for a 10 M_{⊙} star and potentially more than 1000 R_{⊙} for a 100 M_{⊙} star. While some models already included the effects of an external potential in the evolution of dense star clusters (e.g., Leigh et al. 2013a,b, 2014; Boekholt et al. 2018; Sakurai et al. 2019) in order to mimic the effects of gas during the embedded phase of the cluster or to account for the dark matter halo, disentangling the effects of the external potential can be tricky because these models also include mass accretion recipes (Boekholt et al. 2018).
In this paper, we present a systematic investigation on how the formation of very massive stars in dense star clusters and the evolution of the clusters themselves are affected when including a background potential in the evolution of the systems. Our work can be used to constraint the parameter space that future, more sophisticated simulations should focus on. Since each physical collision led to a merger in our runs, we use these terms interchangeably. We describe our setup and initial conditions in Sect. 2 and our results are presented in Sect. 3, including the typical evolution of clusters under the influence of a background potential and clusters without a background potential, the number of collisions, the mass of the resulting object, and our model to estimate the number of mergers at different times. We apply our model to population III (Pop.III) star clusters in Sect. 4. A discussion about neglected effects and considerations for future research on this subject is given in Sect. 5 and finally the conclusions are presented in Sect. 6.
2. Simulation setup
To understand the impact of the gas potential on the number of mergers that occur in star clusters, we first performed simulations in which we did not include stellar mergers, and we compared the evolution of these systems without an external potential and when placed in the center of a background potential. Once we addressed the effects of the extra force on the evolution of the clusters, we performed simulations that include stellar mergers in clusters with and without an external potential. To perform the calculations, we used a modified version of NBODY6^{1} (Aarseth et al. 2000) to treat collisions, where we switched off the stellar evolution package and instead explicitly specify the stellar radii to perform a parameter study. NBODY6 is a fourth order Hermite integrator, and includes a spatial hierarchy to speed up the calculations: This is referred to as the AhmadCohen scheme (Ahmad & Cohen 1973). It also includes routines to treat tidal capture and tidal circularization of binary systems (Mardling & Aarseth 2001) that we activated in our simulations. Another important routine included is the KustaanheimoStiefel regularization (Kustaanheimo & Stiefel 1965), which is an algorithm that can be used to treat binaries and close twobody encounters more accurately. We performed a total of 344 Nbody simulations.
2.1. Simulations without stellar mergers
In order to create a controlled experiment to which we can compare our results including stellar mergers, we first investigated the effects produced by the external potential on the evolution of the star clusters in the absence of collisions. We specifically explored how the mass of the external potential affects the corecollapse, the formation of binary systems, and the evaporation of the clusters. We modeled each cluster as a Plummer (Plummer 1911) distribution of N = 1000 and N = 10 000 equal mass particles with a total cluster mass of M_{stars} = 10^{4} M_{⊙}. Given that we are interested in the overall evolution of the system, and in order to save computational time, we modeled less dense clusters than the ones presented in Sect. 2.2, with a virial radius of R_{v} = 1.0 pc. We then included an analytic background potential that also follows a Plummer density profile with the same virial radius as the stellar distribution R_{v, ext} = 1.0 pc and the mass of the potential was varied as M_{ext} = 0.0, 0.1, 0.5 and 1.0 × M_{stars}. The clusters start in virial equilibrium. We performed a total of eight simulations that do not include stellar mergers and which are listed in Table 1.
Simulations that do not include stellar mergers.
2.2. Simulations that include stellar mergers
Our simulations that include stellar mergers aim to model the first stages after the formation of Pop. III star clusters. We modeled a compact cluster in virial equilibrium by also using a Plummer distribution (Plummer 1911) for the stars that are all equal in mass and radius at the beginning of the simulation. We modeled dense clusters with a virial radius of R_{v} = 0.14 pc and a total mass in stars of M_{stars} = 10^{4} M_{⊙}. We then performed the same set of simulations including an external analytic Plummer potential with M_{ext} = 1.0 × M_{stars} with the same virial radius of R_{v, ext} = 0.14 pc in order to consider at first order the effects of the gas that remains in the cluster after the formation of the stars. Taking the mass of the external potential into account, the crossing time of the clusters were calculated as:
with (see Appendix A); this gives a value of t_{cross} = 0.0078 Myr and t_{cross, ext} = 0.0039 Myr for clusters without and with an external potential, respectively. We investigated how the mass of the final merger product depends on the initial number of stars in the cluster N, which we varied as N = 100, 500, 1000, and 5000 stars. We also varied the initial stellar radius as R_{star} = 20, 50, 100, 200, 500, 1000, and 5000 R_{⊙} using equal radii stars in each run. For each of these configurations, we performed simulations with M_{ext} = 0.0 and 1.0× M_{stars}. Finally, we ran a total of six random simulations per each configuration, which are listed in Table 2. This gives a total of 336 Nbody simulations that include stellar mergers.
Simulations that include stellar mergers.
2.3. Stellar mergers
In our simulations, a merger occurs when two stars are separated by a distance equal to or smaller than the sum of their radii, this is:
with r being the distance between both stars center of mass, and R_{1} and R_{2} are the radii of the stars. Given that the stars are physically in contact at this point, we also call this a stellar collision. Additionally, since each collision in our runs led to a merger, we then consider the terms merger and collision to be synonymous. Once this condition was fulfilled, we merged both stars replacing them by a new star whose total mass is the sum of the masses of the merging stars. Furthermore, the radius was calculated assuming that the star rapidly settles to a stable configuration where the new star has the same density as the previously merging stars, which is consistent with the calculations of Hosokawa et al. (2012) and Haemmerlé et al. (2018). Therefore the new mass and radius of the merger product were calculated as:
3. Results
First we present the results of the simulations that do not include stellar mergers and describe the effects introduced by the addition of an external potential and the variation of the potential mass. Then, we present the results of our simulations that include stellar mergers and describe how their evolution and the formation of massive stars due to the runaway merger process are affected by the external potential.
3.1. Simulations without mergers
Here, we describe the main effects of including an external potential on the evolution of our star cluster models in the absence of stellar mergers. Due to the similarities between the runs, here, we only present the results for simulations number 5 and 8 listed in Table 1 with the setup described in Sect. 2.1. However, a description for the remaining simulations listed in Table 1 is included in Appendix B.
The evolution of clusters with M_{ext} = 0 × M_{stars} is well understood and consistent with the results of Spitzer (1987), where the clusters evolve toward core collapse as a result of twobody relaxation. The core then experiences coreoscillations, contracting and expanding again due to the formation of hard binaries that act as an energy source in the core. The core collapse is expected to occur between 10 and 20 half mass relaxation times t_{rh} that we calculated as (Spitzer 1987):
with γ = 0.4 because we used equal mass stars and t_{cross} was calculated with Eq. (1) setting M_{tot} = M_{stars}. The half mass relaxation time for the cluster with N = 10 000 stars and M_{ext} = 0 × M_{stars} is t_{rh} = 166.38 t_{cross} and core collapse occurs at 2306 t_{cross} which is 14 t_{rh}. We found the core collapse time by visual inspection and looking for the first drop and subsequent rise in the 10% Lagrangian radius as marked by the vertical green line in Fig. 1.
Fig. 1. Evolution of a cluster with N = 10 000 stars with a total stellar mass of M_{stars} = 10^{4} M_{⊙} and virial radius R_{v} = 1.0 pc. Top panel: fraction of ejections from the cluster and bottom panel: Lagrangian radii of the cluster. The vertical green line marks the moment of core collapse. 

Open with DEXTER 
At the moment of core collapse, the 10% Lagrangian radius drops to 0.087 pc, yielding a mean density of 3.6 × 10^{5} M_{⊙} pc^{−3}. After this phase, the entire cluster begins to expand and the onset of ejections takes place with a total number of ejections of about 1500 until a time of 10 000 t_{cross}. A star is considered to be ejected from the cluster if its distance to the center of mass of the system is > 20 R_{v} and its kinetic energy is higher than its potential energy at this location.
When we included an external potential in our simulations, we note first that, since the clusters start virialized, the velocities of the stars are higher compared to clusters without the background potential. This compensates for the extra gravitational force. Then we had to calculate the crossing time of the cluster using Eq. (1) with M_{tot} = M_{stars} + M_{ext}.
We find that all of the clusters behave in a similar way. All of them evolve toward core collapse while there is little expansion of the outer parts (> 30% Lagrangian radii) compared to clusters without an external potential (see Fig. 2).
Fig. 2. Evolution of a cluster with N = 10 000 stars with a total stellar mass of M_{stars} = 10^{4} M_{⊙} and a virial radius of R_{v} = 1.0 pc in the center of an external analytic potential with M_{ext} = 10^{4} M_{⊙}. Top panel: fraction of ejections from the cluster and bottom panel: Lagrangian radii of the cluster. The vertical green line marks the moment of core collapse. 

Open with DEXTER 
For the case when the mass of the external potential is comparable to the total mass (M_{ext} = M_{stars}), the core collapse is delayed significantly, which occurs now at 51 889 t_{cross} (see green vertical line in bottom panel of Fig. 2). This is 23 times (in units of the cluster crossing time) later than the core collapse time for the cluster without a background potential. The 10% Lagrangian radius at this moment drops to 0.073 pc, which yields a mean density of 6.0 × 10^{5} M_{⊙} pc^{−3}. For this simulation, the time is not long enough to see the contraction expansion phase that follows after core collapse. Finally, in this run, only one star is ejected from the cluster due to the higher escape velocity compared to the cluster without the external potential.
In order to understand the delay in the core collapse time when including a background potential, we define the parameter q = M_{ext}/M_{stars} as the ratio between the mass of the external potential and the total mass in stars, then we find that the core collapse time scales as t_{cc} ∝ (1 + q)^{4}t_{cross} as presented in Fig. 3. If we assume, as found by Spitzer (1987) that the core collapse time is proportional to the relaxation time, then according to our data the relaxation time is proportional to t_{rh} ∝ (1 + q)^{4}, which is consistent with Eq. (A.12) for which we derived the relaxation time of a cluster in the center of an analytic external potential (see Appendix A):
Fig. 3. Core collapse time t_{cc} divided by the crossing time t_{cross} for clusters with N = 1000 (blue circles) and N = 10 000 (red circles) as a function of q = M_{ext}/M_{stars}. We also plotted two relations that scale as (1 + q)^{4} which fit our data and therefore suggest that if the core collapse time is proportional to the relaxation time of the cluster, then the relaxation time scales as (1 + q)^{4}, which is the same relation that we derive in Appendix A. 

Open with DEXTER 
where the crossing time must be calculated from Eq. (1) with q = M_{ext}/M_{stars} (see the Appendix in Leigh et al. 2013a and Leigh et al. 2014 for similar derivations but slightly modified to consider different cases).
A second important effect that we note is on the binary population. Our results show that in increasing the mass of the background potential, the number of binary systems that form decreases dramatically as shown in Fig. 4 due to the increased mean kinetic energy of the stars. Consequently, we also find that these binaries tend to be more tightly bound (see Fig. 5) given that they form at the hardsoft boundary, which is smaller for higher velocity dispersion (Leigh et al. 2015). This will certainly have an impact on the growth of a massive star through stellar mergers given that collisions are more likely to occur in binary systems.
Fig. 4. Number of binary systems as a function of time and for different values for the mass of the external potential M_{ext} for simulations with N = 1000 and up to 8000 crossing times. The vertical lines mark the moments of core collapse and coincides with the onset of the formation of binary systems. 

Open with DEXTER 
Fig. 5. Distribution of the semimajor axis of the binary population for the set of simulations listed in Table 1. The data clearly shows that increasing the mass of the external potential the formed binaries tend to be more tightly bound. Despite the fact that a low mass of the external potential (M_{ext} < 1.0 × M_{stars}) may promote the formation of wider binaries, these systems tend to be shortlived due to the higher velocities of all the stars. This figure is only intended to show the distribution of the semimajor axis and not the total number of binaries that formed. 

Open with DEXTER 
3.2. Simulations that include stellar mergers
In the following, we present the results that we obtained when we included stellar mergers in the evolution of star clusters with and without an external Plummer potential. Additionally, we describe the effects introduced by this background potential on the formation of massive merger products. These simulations are listed in Table 2.
The first effect that we note is a delay in the runaway growth of the central object, which is explained by the delay in the core collapse time due to a larger relaxation time when increasing the mass of the potential, as can be seen in Eq. (6). Moreover with this external force, binary systems are harder to form, as shown in Fig. 4 and these binaries tend to be more compact.
Here we describe the time evolution for two clusters with N = 1000 stars, with and without an external potential, which correspond to simulations number 17 and 45 in Table 2 and we also present their evolution in Figs. 6 and 7, respectively. These simulations serve as examples of the typical evolution of the runs that include stellar mergers, given that all the simulated clusters follow very similar patterns.
Fig. 6. Evolution of a cluster with N = 1000 stars, total mass M_{stars} = 10^{4} M_{⊙}, R_{v} = 0.14 pc, and R_{star} = 100 R_{⊙}. Uppermost panel: number of mergers as a function of time, the fraction of binary, hyperbolic, and chain mergers in the second panel along with the fraction of stars that escape from the cluster. Third panel: 10, 50 and 90% Lagrangian radii, and bottom panel: mass growth of the most massive object that formed through mergers in the cluster. 

Open with DEXTER 
Fig. 7. Evolution of a cluster with N = 1000 stars, stellar mass M_{stars} = 10^{4} M_{⊙}, R_{v} = 0.14 pc, and R_{star} = 100 R_{⊙} in the center of an external potential with M_{ext} = 10^{4} M_{⊙}. Uppermost panel: number of mergers as a function of time, the fraction of binary, hyperbolic, and chain mergers in the second panel, along with the fraction of stars that escape from the cluster. Third panel: 10, 50, and 90% Lagrangian radii, and bottom panel: mass growth of the most massive object that formed through mergers in the cluster. The drop in the 50% Lagrangian radius toward the end of the simulation suggests that the central object may grow up to 5000 M_{⊙}. 

Open with DEXTER 
In our runs, we identify three types of mergers: namely hyperbolic mergers, which are mergers that occur between stars that are not gravitationally bound; binary mergers, which occur between stars that are part of a binary system; and chain mergers, which occur between stars that are part of a higher order system, for example, triples. Before the onset of runaway mergers occurs in all of the clusters, there is a nearly constant merger rate for hyperbolic mergers which do not produce a single massive star, but instead several less massive stars that eventually sink to the cluster center, thus finally contributing to the growth of the most massive star. We note that the hyperbolic merger rate in clusters that do not include an external potential is only maintained for short periods of time (see Fig. 6); whereas in clusters that include the extra force, the hyperbolic merger rate is maintained until the end of the simulation (see Fig. 7). Chain mergers are very rare with only a handful having been identified among all of our simulations.
During the period before the rapid growth of the most massive star, relaxation processes drive an energy flow from the central parts of the cluster (inside the halfmass radius) via high velocity stars that migrate to the cluster halo and cause the expansion of the outer parts. All of the clusters evolve toward core collapse, but not all of them reach that stage because for large stellar radii R_{star} ≥ 100 R_{⊙}, the growth of the central star is fast and the drop in the 10% Lagrangian radius that we see in Fig. 6 occurs because the central object has already gathered 10% of the cluster mass, and thus we see the position of the central star rather than a collapsed core. The onset of growth for the central star is also marked by a rapid increase in the merger rate of stars in binary systems despite the presence of an external potential.
An important difference that we note when including the background potential is on the expansion of the clusters and the fraction of escaping stars. In this sense we found that the external potential prevents the evaporation of the cluster (see panel 2 of Fig. 7) and keeps it compact so the mergers can continue for longer periods of time, potentially gathering up to 50% or more of the cluster mass in the central object, as the third panel of Fig. 7 suggests.
3.3. Formation and growth of the most massive star
One important goal of this research is to address the effect of the external potential on the formation of a massive star by stellar mergers. Thus, we present the results regarding the mass growth of the most massive star in our clusters, along with a description of the effects introduced by the background potential.
The first important effect is a delay in the formation of the most massive star due to the increased stellar velocity, which in turn causes relaxation processes to be slower. This then leads to a reduction in the number of binary systems that form, at a given time, if the mass of the external potential is comparable to the total stellar mass. On the other hand, the external potential can also favor the mass growth of the central star by preventing the evaporation and expansion of the cluster.
In the following, we present a model that we used to estimate the mass of the central object at different times for both a cluster with and without an external potential. For this, we followed the same method described in Reinoso et al. (2018).
In the lower panel of Fig. 8, we first present an example of the mass growth of the most massive star, which undergoes a very rapid growth at around 1000–1500 t_{cross}. This coincides with the peak in the number of mergers with the central object as shown in the top panel of the same figure.
Fig. 8. Example of a Gaussian fit to the number of mergers with the central star for simulation number 45 listed in Table 2. Top panel: number of mergers with the most massive star as a function of time (black solid line) along with the Gaussian function that we used to estimate the total number of mergers at different times. The parameters t_{delay} and t_{duration} are also shown in the figure. Bottom panel: mass growth of the central star whose rapid growth coincides with the peak in the number of mergers. 

Open with DEXTER 
We fit the Gaussian function presented in Eq. (7) to the combined data from the six random realizations for each simulation setup presented in Table 2 in order to get an estimate for the number of mergers with the central object during the rapid growth. Thus, we define the delay time t_{delay} as the time at which the peak in the Gaussian occurs, and the height A of the Gaussian gives us an estimate for the number of mergers at t_{delay}. Additionally, we define the duration time , with FWHM being the full width at half maximum of the Gaussian. An illustration of the fitting formula and two parameters are shown in the top panel of Fig. 8. By doing this, we were able to obtain an estimate for the moment at which there is rapid growth of the most massive star (t_{delay}) and an estimate for the duration of this period (t_{duration}), which we can later combine with A and a normalization factor D to finally obtain the number of mergers experienced by the central star using Eq. (8) and therefore an estimate for its final mass at different times. Furthermore we determined D by comparing our model to the results from our simulations. We also show how well our model reproduces the data of our simulations in Figs. 9 and 10. The total number of mergers at time t can thus be calculated as:
Fig. 9. Mean mass of the most massive star normalized by its initial mass, as a function of N and R_{star} for clusters without an external potential evolved until 2000 t_{cross} = 15.6 Myr. The symbols represent the mean values obtained from six simulations and the solid lines show the results of the calculations with our model presented in Eq. (8) using t = 15.6 Myr and the parameters obtained from Eqs. (9)–(12). 

Open with DEXTER 
Fig. 10. Mean mass of the most massive star normalized by its initial mass, as a function of N and R_{star} for clusters without an external potential evolved until 3000 t_{cross} = 16.5 Myr, except for the point at N = 5000 and R_{star} = 20 R_{⊙}, which are simulations that were evolved until 4000 t_{cross} = 22 Myr. The symbols represent the mean values obtained from six simulations and the solid lines show the results of the calculations with our model presented in Eq. (8) using t = 16.5 Myr (t = 22 Myr for N = 5000 and R_{star} < 50 R_{⊙}) and the parameters obtained from Eqs. (13)–(16). 

Open with DEXTER 
and the total number of mergers until a time t_{end} can be calculated as:
Although the data is not a perfect Gaussian (see Fig. 8), we only need an estimate for the delay time t_{delay}, the duration time t_{duration}, and an estimate for the number of collisions A at t_{delay}. Following this procedure, we first combined the data of the six random realizations per each of the configurations listed in Table 2 and then applied the Gaussian fit of Eq. (7) to find the values of A, t_{delay}, and t_{duration} that we present in Figs. 11 and 12 for clusters with M_{ext}/M_{stars} = 0.0 (simulations 1–28 in Table 2) and in Fig. 13 for clusters with M_{ext}/M_{stars} = 1.0 (simulations 29–56 in Table 2). Then we fit these data using an implementation of the nonlinear leastsquares MarquardtLevenberg algorithm in gnuplot and obtained Eqs. (9)–(11) for clusters without an external potential and Eqs. (13)–(15) for clusters in a background potential. Using these equations, we were able to compute N_{merger, total} assuming first D = 1 and we compared the calculated values to the real values from the simulations to adjust D and reproduce the results correctly. By doing this, we obtain Eqs. (12) and (16):
Fig. 11. Parameters of the Gaussian fit to the number of mergers with the central star for simulations of clusters without an external potential. The data points were obtained after the fit to the combined data of six simulations. Upper panel: values of the normalization factor A for different values of N and R_{star} along with our model from Eq. (9). These data points show that the number of mergers increases with both N and R_{star} as expected. Bottom panel: different values of t_{delay} in addition to the different values of N and R_{star} along with the fitting function presented in Eq. (10). These points show that t_{delay} decreases with R_{star}, but it increases with N. 

Open with DEXTER 
Fig. 12. t_{duration} as a function of the number of stars N for clusters without an external potential. We found no clear correlations between t_{duration} and R_{star}. The solid blue line shows the fitting function from Eq. (11). 

Open with DEXTER 
Fig. 13. Parameters of the Gaussian fit to the number of mergers with the central star for simulations of clusters with an external potential. The data points were obtained after the fit to the combined data of six simulations. Upper panel: values of the normalization factor A_{ext} for different values of N and R_{star} along with our model from Eq. (8). These data points show that the number of mergers increases with both N and R_{star} as expected. Middle panel: different values of t_{delay, ext} in addition to the different values of N and R_{star} along with the fitting function presented in Eq. (14). These points show that t_{delay, ext} decreases with R_{star}, but it increases with N. Additionally, these values are larger than the values of t_{delay} for clusters without an external potential (see Fig. 11). Lower panel: values of t_{duration, ext} along with the fitting function from Eq. (15). We clearly see a positive correlation with N and a negative correlation with R_{star}. 

Open with DEXTER 
and the parameters for clusters in an external potential:
It is important to note that in Eqs. (10) and (11), the delay time t_{delay} and duration time t_{duration} are expressed in units of the crossing time of the cluster defined in Eq. (1) and consequently when calculating the total number of mergers until a time t_{end} using Eq. (8), this time must also be expressed in units of the crossing time t_{cross}. The same principle applies for Eqs. (14) and (15). That is, when calculating the number of mergers until a time t_{end} using Eq. (8), this time must be expressed in units of the cluster crossing time that have been modified by the mass of the external potential t_{cross, ext}, which is defined in Eq. (1).
The delay time for clusters in a background potential is defined in Eq. (14) as t_{delay, ext} = t_{delay} + t_{delay, Mext = Mstars}, which is the delay time for clusters without an external potential t_{delay} plus an additional term t_{delay, Mext = Mstars} = [ − 0.006log_{10}(N)+0.240]log_{10}(R_{star})+0.064log_{10}(N) − 1.170 that, in principle, may depend on the mass of the potential. However we do not examine this potential dependence here.
To help in the reading of the equations, we define a general equation of the form:
and present the parameters α, β, and γ in Table 3. They can be used to recover Eqs. (9)–(16).
Now we can estimate the number of mergers with the central star and thus the mass of this star up to a time t for clusters with different initial conditions, including the presence of an external potential that is comparable to the total stellar mass in the cluster. This can be used as a very simplified model of a star cluster, which still contains gas during or shortly after the process of star formation or a nuclear star cluster. In Fig. 14, we show the expected fraction f_{merger} of the stars that merge with a single central object in a cluster without an external potential. In Fig. 15, this is shown for clusters in an external potential, for a broad combination of N and R_{star} after 1, 2, and 10 Myr assuming a cluster with M_{stars} = 10^{4} M_{⊙}, a virial radius of R_{v} = 0.14 pc, and a mass of the external potential M_{ext} = M_{stars} for the second case. We used Eq. (8) for these calculations along with the parameters from Eqs. (9)–(12) to obtain the results presented in Fig. 14, and Eqs. (13)–(16) to obtain the results presented in Fig. 15. Our model gives the total number of stars N_{merger, total} that merge with the central object in a given time interval, and thus assuming that all stars are equal during the time at which most of the collisions occur, we can estimate the mass of the central object as . Our results indicate that when collisions are maintained for short periods of time only, that is, 1 Myr, then the clusters without an external potential form more massive central stars than clusters with an external potential (see left panels of Figs. 14 and 15). This is due to the fact that the external potential dramatically increases the relaxation time of the cluster, as shown in Eq. (6). In fact the relaxation time depends on the ratio of external potential mass to stellar mass q = M_{ext}/M_{stars} as t_{relax} ∝ q^{4}.
Fig. 14. Fraction of stars f_{merger} that merge with the central star in a cluster with a total stellar mass of M_{stars} = 10^{4} M_{⊙} and a virial radius of R_{v} = 0.14 pc after 1 Myr (left panel), 2 Myr (middle panel), and 10 Myr (right panel) for different combinations of the number of stars N and the stellar radius R_{star} calculated using our model presented in Eq. (8). We can see that when the collisions are maintained for short periods of time, that is, 1 Myr, a massive object with ∼f_{merger}M_{stars} = 0.15 × 10^{4} M_{⊙} = 1500 M_{⊙} is formed in the cluster center only when N is small (given that t_{delay} becomes larger than the integration time for large N) or large R_{star}. If the time limit is long, that is, 10 Myr, most of the clusters will form stars with ∼0.3 M_{stars} = 3000 M_{⊙}. 

Open with DEXTER 
Fig. 15. Fraction of stars f_{merger} that merge with the central star in a cluster with a total stellar mass of M_{stars} = 10^{4} M_{⊙}, a virial radius of R_{v} = 0.14 pc, and an external potential with a mass of M_{ext} = M_{stars} after 1 Myr (left panel), 2 Myr (middle panel), and 10 Myr (right panel) for different combinations of the number of stars N and the stellar radius R_{star} calculated using our model presented in Eq. (8). We can see that when the collisions are maintained for short periods of time, that is, 1 Myr, an object that is not so massive with ∼f_{merger}M_{stars} = 0.1 × 10^{4} M_{⊙} = 1000 M_{⊙} is formed in the cluster center only when N is very small (given that t_{delay} becomes larger than the integration time for large N) or when R_{star} is large. If the time limit is 2 Myr (middle panel), we expect that clusters with a small N experience core collapse and form a massive central star. However, if the number of stars is large and R_{star} is also large, given the higher velocity of the stars in those clusters compared to clusters without the external potential, there is a larger probability for hyperbolic mergers to occur before core collapse, which we see on the righthand side of middle panel. If the time limit is long, that is, 10 Myr, most of the clusters form stars with ∼0.3−0.4 M_{stars} = 3000−4000 M_{⊙} and even 0.7 M_{stars} = 7000 M_{⊙} if N is large enough. 

Open with DEXTER 
As the time limit increases, we see that in both clusters with and without an external potential a more massive star emerges in the center. However, even more massive stars form in the cluster without an external potential for most values of N. However, there is an important difference in clusters that include a background potential, which is that the potential keeps the cluster compact as previously found by Leigh et al. (2014), and thus most of the cluster mass is able to eventually sink to the center at later times. This effect is visible in the left part of the middle panel in Figs. 14 and 15. We also note that for large N and R_{star} for clusters in a background potential and before core collapse, we may also expect mergers with the central star, not binary mergers but hyperbolic, given that in those cases the cross section for collisions is very large and the probability of hyperbolic collisions increases with stellar velocity, that is, if there is an external potential. Finally when the time limit is very long, that is, 10 Myr, we may expect final masses for the central object in the order of 0.3–0.4 M_{stars} for most clusters without an external potential and 0.4–0.7 for clusters in a background potential because the compactness of the cluster is maintained and collisions still occur at an approximately constant rate after the stage of runaway growth of the central star.
4. Implications for primordial clusters
In the following, we explore the implications of our results with respect to primordial star clusters, including both embedded and gas free. Here we particularly distinguish the case of standard Pop. III clusters as expected in a typical minihalo with about 10^{6} M_{⊙}.
4.1. Standard Pop. III clusters (minihalos)
For a typical Pop. III star cluster, we assume a mass of 1000 M_{⊙}, which is consistent with a baryon fraction of about 10% and a star formation efficiency on the order of 1% in a 10^{6} M_{⊙} minihalo. We adopt a virial radius for the cluster on the order of 0.1 pc, which is consistent with the results from simulations and semianalytic models (Clark et al. 2011a,b; Greif et al. 2011, 2012; Latif et al. 2013b; Latif & Schleicher 2015). We use a stellar radius of 100 R_{⊙}, which is characteristic of primordial protostars with accretion rates on the order of 10^{−3} M_{⊙} yr^{−1} (Hosokawa et al. 2012). The crossing time of the cluster then is 0.015 Myr. The number of stars that can be expected in such a cluster is uncertain, but here we adopt an estimate of about N = 100. Thus, a stellar mass of m_{star} = 10 M_{⊙}. Using the relations that we found, we expect two collisions to occur within 1 Myr which correspond to a final mass of 30 M_{⊙}. If we assume now a lifetime of 10 Myr, we can expect a total of six collisions with the central star which correspond to a final mass of 70 M_{⊙}.
We consider now an embedded and accreting Pop. III star cluster with a stellar mass of 1000 M_{⊙}, a gas mass of 1000 M_{⊙}, a virial radius of 0.1 pc, and a mean stellar radius of 100 R_{⊙}. The crossing time of the cluster, according to Eq. (1), corresponds to 0.007 Myr. We also adopt a number of stars N = 100. Then our model predicts that no collisions occur within 1 Myr, but for a lifetime of 10 Myr we expect a total of 15 collisions and hence a central star with 160 M_{⊙}. We expect the lifetime of a massive primordial star to be in between this range, depending on the precise mass, the amount of rotation, and the effects that the collisions may have on the stellar evolution (Maeder & Meynet 2012). We therefore find that a moderate enhancement can be achieved within a normal cluster. We note that the values given here are the expected mean number of mergers. Individual clusters can deviate from these, both towards lower and higher fractions of mergers, potentially including clusters with zero mergers. This is especially the case for typical Pop. III clusters where the number of collisions can be expected to be comparable to the mean value because it is a chaotic collisional dynamical process.
4.2. Massive primordial clusters (atomic cooling halos)
As a next step, we address now the potential impact of collisions in a more massive atomic cooling halo with a total mass of 10^{8} M_{⊙}. Under the right conditions and in particular if the cooling on larger scales is regulated by atomic hydrogen (Latif et al. 2014), a rather massive cluster of 10^{4} M_{⊙} can form, which is then exposed to larger accretion rates on the order of 10^{−1} M_{⊙} yr^{−1}. We assume that the cluster consists of an initial number N = 1000 of stars, the virial radius is R_{v} = 0.14 pc, and that the stellar radii are somewhat enhanced compared to the standard Pop. III cluster due to the higher accretion rates, with a typical radius of about 300 R_{⊙}. The crossing time of the cluster according to Eq. (1) with M_{tot} = 10^{4} M_{⊙} is then 0.0078 Myr. Using our derived model we expect about 51 collisions in 1 Myr, and about 168 within 10 Myr with a single central star. We again expect the realistic lifetime of the resulting massive star to be in between these extreme cases. In the case of an atomic cooling halo, we thus conclude that a considerable enhancement is possible as a result of stellar mergers. If we take the mean stellar mass, then the expected masses are 520–1690 M_{⊙} after 1 and 10 Myr, respectively.
If we now consider the same cluster configuration but during the embedded phase, assuming a total mass in gas of 10^{4} M_{⊙}, the crossing time of the cluster according to Eq. (1) with q = M_{ext}/M_{stars} = 2 is then 0.0039 Myr, and we expect a total of 16 collisions to occur within 1 Myr and a total of 199 mergers within 10 Myr. Again, by taking the mean stellar, mass this corresponds to a mass of 170–2000 M_{⊙} within 1 and 10 Myr, respectively. Also the reported values here correspond to a mean, and there can be deviations to lower and higher merger fractions. Regardless, we expect the number of collisions and hence the final mass to remain within the same order of magnitude.
5. Discussion
We modeled the runaway growth of stars through mergers in the center of dense star clusters by including an analytic external potential in our Nbody simulations. Our model relies on computing four parameters that depend on the number of stars N and their stellar radii R_{star}. These parameters are the delay time t_{delay}, the duration time for the collision process t_{duration}, and two normalization factors A and D. We then used these four parameters to integrate a Gaussian function and estimate the mass of the central object up to a time t. Although a Gaussian fit is a good tool to estimate t_{delay}, t_{duration}, and A (see Fig. 8), there is a deviation from this function when we include an external potential. In fact, when we compare the mass growth of the central object in clusters with and without an external potential (see, e.g. Figs. 6 and 7), we see a delay for clusters in an external potential and an initially less dramatic growth until the onset of mergers of stars in binary systems. This growth continues after the runaway growth and is seen as tails of a Gaussian in a plot of N_{mergers} versus time (see top panel of Figs. 9 and 10). This later growth is associated to hyperbolic mergers, that is, mergers of stars that are not bound by gravity, and this type of merger becomes important when including an external potential because the cluster remains compact for longer periods of times (Leigh et al. 2014). The onset of collisions is delayed in simulations with an external potential relative to simulations without it.
We also have found a formula for the delay time when we include an external potential t_{delay, ext} = t_{delay} + t_{delay, Mext = Mstars}, which is basically the same expression for clusters without the external potential plus an additional term t_{delay, Mext = Mstars} presented in Eq. (14). We present this parameter in this way so that we can investigate in the future if this can be used in a more general way as a function of the mass of the external field M_{ext} and if it is related to a longer relaxation time by means of its dependence on M_{ext} as suggested by Eq. (6).
The model derived here can be used as a first approximation to obtain the number of stellar mergers in embedded star clusters or to understand the effects of an external potential on the formation of massive merger products, which might also be important for the modelling of nuclear star clusters. Future research in this field will employ hydrodynamic modeling of the gas, which may prevent a large delay in the onset of runaway growth for the central star due to gas accretion and dissipative effects via stargas or binarygas interactions. These simulations also need to include realistic massradius relations for accreting protostars when aiming for a realistic modeling of the first star clusters in the Universe. By performing these suites of Nbody simulations we have covered a large parameter space, and this allows us to provide hints as to which part of this parameter space the future, more sophisticated simulations should be focused on and the subsequent physics that becomes relevant and often even dominant in these regimes.
Our simulations including mergers adopt a virial radius of R_{v} = 0.14 pc and a total stellar mass of M_{stars} = 10^{4} M_{⊙}, but our results can be rescaled for different sizes and masses by means of the crossing time. In this work, we only consider for the initial conditions equal mass and equal radii stars but allow for the mass and radius to vary due to the mergers according to Eqs. (3) and (4), and in this sense our simulations represent an effective model where an average stellar radius is adopted over the period of time considered. This effective stellar radius should correspond to the typical radius when the majority of mergers are expected to occur.
Our model still lacks a more realistic stellar population, namely an initial mass function (IMF), which would naturally lead to mass segregation and subsequent mergers of the most massive stars in the cluster center. Moreover, we have not included the dissipative effects expected for stargas interactions in our simulations that include an external potential, which could lead to the formation of more binary systems and hence more mergers (Leigh et al. 2014). In this sense, we do not expect a significant variation in the final mass of the most massive stars if the total mass in the most massive stars is comparable to the final mass obtained in our models, around 30–2000 M_{⊙}. Additionally, including an IMF, threebody encounters and twobody relaxation would cause the ejection of the smallest stars in the cluster, which could be an important source of stellar relics from the first star clusters in the early Universe. Another important way to proceed is a better modeling of the stellar mergers, in particular the mass loss during the stellar collisions which may be up to 25% of the total mass (Gaburov et al. 2010) and the collision product may even receive a kick velocity > 10 km s^{−1}. More recent work on the effects of mass loss and its effect on runaway growth of a central star in a cluster shows that considering 5% mass loss in every collision produces objects which are 20–40% less massive at the end of the runaway growth compared to models without mass loss (Glebbeek et al. 2009; Alister Seguel et al. 2020). Thus regarding this, the results presented here should be taken as an upper limit on the mass of the merger product.
6. Conclusions
In this study, we have performed a set of 374 Nbody simulations of dense, virialized star clusters and clusters in the center of an external potential including stellar mergers in order to understand the effects of the background force on the formation of massive stars and derive a model to estimate the mass enhancement of these objects in embedded star clusters. We find that the presence of an external potential delays the overall evolution of the star cluster and the formation of a massive central star through runaway collisions due to the increased kinetic energy of the stars, which in turn increases the relaxation time. However, the merger products become more massive (if the collision process is maintained for a long time) given that these clusters expand less than clusters without an external potential and so more stars are able to merge with the central star even after the process of runaway growth.
We also find that the increased velocity dispersion for star clusters in an external potential boosts the number of hyperbolic mergers , that is, mergers of stars that are not part of a binary or triple system, with up to 50–60% of the total number of mergers being hyperbolic mergers. Whereas in clusters without an external potential, this percentage is around 30%. This is due to the fact that in presence of a background potential, the clusters remain more compact (Leigh et al. 2013a, 2014) and hyperbolic mergers still may occur. Stellar ejections are highly suppressed in clusters in the center of an external potential.
We find a set of equations that can be used to estimate the mass of the merger product at different times for both clusters with and without an external potential, and we present an example of such calculations in Figs. 14, 15, and in Sect. 4 proving that if the process is not interrupted for long periods of time (≥10 Myr), the external potential enhances the mass of the central object by a factor of ∼2. However, if the process is interrupted at early times (∼1 Myr), the clusters in a background potential produce objects only half as massive compared to objects formed in clusters without an external potential. When applied to Pop. III star clusters we find, for standard clusters formed in a minihalo, a moderate enhancement for the mass of the most massive star in the range of 10–160 M_{⊙} within 1 and 10 Myr in embedded clusters. Whereas for a massive Pop. III star cluster that formed in an atomic cooling halo, we find that the mass of the most massive star lies in the range of 170–2000 M_{⊙} in embedded clusters within 1 and 10 Myr.
Acknowledgments
We thank the anonymous referee for the useful comments. BR thanks Sverre Aarseth for his help with the code NBODY6. BR also thanks Conicyt for funding through (CONICYTPFCHA/MagísterNacional/201722171385), (CONICYTPFCHA/Doctorado acuerdo bilateral DAAD/62180013), financial support from DAAD (funding program number 57451854). DRGS acknowledges support of Fondecyt Regular 1161247 and through Conicyt PIA ACT172033. BR and DRGS acknowledge financial support through Anillo (project number ACT172033). MF acknowledges funding through Fondecyt Regular No.1180291. NWCL gratefully acknowledges the support of a Fondecyt Iniciacion #11180005.
References
 Aarseth, S. J. 2000, in The Chaotic Universe, eds. V. G. Gurzadyan, & R. Ruffini, 286 [CrossRef] [Google Scholar]
 Ahmad, A., & Cohen, L. 1973, J. Comput. Phys., 12, 389 [NASA ADS] [CrossRef] [Google Scholar]
 Alister Seguel, P. J., Schleicher, D. R. G., Boekholt, T. C. N., Fellhauer, M., & Klessen, R. S. 2020, MNRAS, 493, 2352 [CrossRef] [Google Scholar]
 Binney, J., & Tremaine, S. 1987, Galactic Dynamics (Princeton, N.J.: Princeton University Press) [Google Scholar]
 Boekholt, T. C. N., Schleicher, D. R. G., Fellhauer, M., et al. 2018, MNRAS, 476, 366 [NASA ADS] [CrossRef] [Google Scholar]
 Bonetti, M., Rasskazov, A., Sesana, A., et al. 2020, MNRAS, 493, L114 [CrossRef] [Google Scholar]
 Bovino, S., Grassi, T., Schleicher, D. R. G., & Banerjee, R. 2016, ApJ, 832, 154 [NASA ADS] [CrossRef] [Google Scholar]
 Clark, P. C., Glover, S. C. O., Klessen, R. S., & Bromm, V. 2011a, ApJ, 727, 110 [NASA ADS] [CrossRef] [Google Scholar]
 Clark, P. C., Glover, S. C. O., Smith, R. J., et al. 2011b, Science, 331, 1040 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Fan, X., Strauss, M. A., Richards, G. T., et al. 2006, AJ, 131, 1203 [NASA ADS] [CrossRef] [Google Scholar]
 Fujii, M. S., & Portegies Zwart, S. 2013, MNRAS, 430, 1018 [NASA ADS] [CrossRef] [Google Scholar]
 Gaburov, E., Lombardi, J. C., & Portegies Zwart, S. 2010, MNRAS, 402, 105 [CrossRef] [Google Scholar]
 Glebbeek, E., Gaburov, E., de Mink, S. E., Pols, O. R., & Portegies Zwart, S. F. 2009, A&A, 497, 255 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Greif, T. H., Springel, V., White, S. D. M., et al. 2011, ApJ, 737, 75 [NASA ADS] [CrossRef] [Google Scholar]
 Greif, T. H., Bromm, V., Clark, P. C., et al. 2012, MNRAS, 424, 399 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Haemmerlé, L., Woods, T. E., Klessen, R. S., Heger, A., & Whalen, D. J. 2018, MNRAS, 474, 2757 [NASA ADS] [CrossRef] [Google Scholar]
 Hosokawa, T., Omukai, K., & Yorke, H. W. 2012, ApJ, 756, 93 [NASA ADS] [CrossRef] [Google Scholar]
 Hosokawa, T., Yorke, H. W., Inayoshi, K., Omukai, K., & Yoshida, N. 2013, ApJ, 778, 178 [NASA ADS] [CrossRef] [Google Scholar]
 Katz, H., Sijacki, D., & Haehnelt, M. G. 2015, MNRAS, 451, 2352 [NASA ADS] [CrossRef] [Google Scholar]
 Klessen, R. S., Glover, S. C. O., & Clark, P. C. 2012, MNRAS, 421, 3217 [NASA ADS] [CrossRef] [Google Scholar]
 Kustaanheimo, P., & Stiefel, E. 1965, J. Reine Angew. Math., 218, 204 [Google Scholar]
 Latif, M. A., & Schleicher, D. R. G. 2015, MNRAS, 449, 77 [NASA ADS] [CrossRef] [Google Scholar]
 Latif, M. A., Schleicher, D. R. G., Schmidt, W., & Niemeyer, J. 2013a, ApJ, 772, L3 [NASA ADS] [CrossRef] [Google Scholar]
 Latif, M. A., Schleicher, D. R. G., Schmidt, W., & Niemeyer, J. C. 2013b, MNRAS, 436, 2989 [NASA ADS] [CrossRef] [Google Scholar]
 Latif, M. A., Schleicher, D. R. G., Bovino, S., Grassi, T., & Spaans, M. 2014, ApJ, 792, 78 [NASA ADS] [CrossRef] [Google Scholar]
 Latif, M. A., Omukai, K., Habouzit, M., Schleicher, D. R. G., & Volonteri, M. 2016, ApJ, 823, 40 [NASA ADS] [CrossRef] [Google Scholar]
 Leigh, N., Sills, A., & Böker, T. 2013a, MNRAS, 433, 1958 [CrossRef] [Google Scholar]
 Leigh, N. W. C., Böker, T., Maccarone, T. J., & Perets, H. B. 2013b, MNRAS, 429, 2997 [NASA ADS] [CrossRef] [Google Scholar]
 Leigh, N. W. C., MastrobuonoBattisti, A., Perets, H. B., & Böker, T. 2014, MNRAS, 441, 919 [NASA ADS] [CrossRef] [Google Scholar]
 Leigh, N. W. C., Giersz, M., Marks, M., et al. 2015, MNRAS, 446, 226 [NASA ADS] [CrossRef] [Google Scholar]
 Maeder, A., & Meynet, G. 2012, Rev. Mod. Phys., 84, 25 [Google Scholar]
 Mardling, R. A., & Aarseth, S. J. 2001, MNRAS, 321, 398 [NASA ADS] [CrossRef] [Google Scholar]
 Mortlock, D. J., Warren, S. J., Venemans, B. P., et al. 2011, Nature, 474, 616 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Omukai, K., & Palla, F. 2001, ApJ, 561, L55 [NASA ADS] [CrossRef] [Google Scholar]
 Omukai, K., & Palla, F. 2003, ApJ, 589, 677 [NASA ADS] [CrossRef] [Google Scholar]
 Omukai, K., Tsuribe, T., Schneider, R., & Ferrara, A. 2005, ApJ, 626, 627 [NASA ADS] [CrossRef] [Google Scholar]
 Plummer, H. C. 1911, MNRAS, 71, 460 [NASA ADS] [CrossRef] [Google Scholar]
 Reed, S. L., Banerji, M., Becker, G. D., et al. 2019, MNRAS, 487, 1874 [NASA ADS] [CrossRef] [Google Scholar]
 Reinoso, B., Schleicher, D. R. G., Fellhauer, M., Klessen, R. S., & Boekholt, T. C. N. 2018, A&A, 614, A14 [CrossRef] [EDP Sciences] [Google Scholar]
 Sakurai, Y., Yoshida, N., Fujii, M. S., & Hirano, S. 2017, MNRAS, 472, 1677 [NASA ADS] [CrossRef] [Google Scholar]
 Sakurai, Y., Yoshida, N., & Fujii, M. S. 2019, MNRAS, 484, 4665 [CrossRef] [Google Scholar]
 Schleicher, D. R. G., Palla, F., Ferrara, A., Galli, D., & Latif, M. 2013, A&A, 558, A59 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Smith, R. J., Glover, S. C. O., Clark, P. C., Greif, T., & Klessen, R. S. 2011, MNRAS, 414, 3633 [NASA ADS] [CrossRef] [Google Scholar]
 Smith, R. J., Hosokawa, T., Omukai, K., Glover, S. C. O., & Klessen, R. S. 2012, MNRAS, 424, 457 [NASA ADS] [CrossRef] [Google Scholar]
 Spitzer, L. 1987, Dynamical Evolution of Globular Clusters (Princeton, N.J.: Princeton University Press) [Google Scholar]
 Stahler, S. W., Palla, F., & Salpeter, E. E. 1986, ApJ, 302, 590 [NASA ADS] [CrossRef] [Google Scholar]
 Suazo, M., Prieto, J., Escala, A., & Schleicher, D. R. G. 2019, ApJ, 885, 127 [CrossRef] [Google Scholar]
 Wise, J. H., Regan, J. A., O’Shea, B. W., et al. 2019, Nature, 566, 85 [CrossRef] [Google Scholar]
 Woods, T. E., Heger, A., Whalen, D. J., Haemmerlé, L., & Klessen, R. S. 2017, ApJ, 842, L6 [NASA ADS] [CrossRef] [Google Scholar]
 Wu, X.B., Wang, F., Fan, X., et al. 2015, IAU General Assembly, 29, 2251223 [Google Scholar]
Appendix A: Modification of the halfmass relaxation timescale
We have seen that when we include an external potential in a star cluster, the overall evolution seems to be delayed, in particular the core collapse timesscale rapidly increases with the mass of the external potential. Until now, we have used the usual halfmass relaxation timescale defined in Eq. (5) (Spitzer 1987). However, this is probably not adequate when we include an external potential because in that case, considering a virialized cluster, the root mean square (rms) velocity of the stars is higher due to the extra force exerted by the potential and this certainly modifies the timescale on which stellar encounters are going to modify the higher velocity of the stars. In order to account for this effect, here, we derive a modified relaxation timescale with the aid of a new parameter q = M_{ext}/M_{stars}. We begin with the usual derivation of the relaxation time assuming virial equilibrium, a similar derivation can be found in Binney & Tremaine (1987):
We note that is the kinetic energy and is the potential energy. In this case, M_{stars} is the total mass in stars, ⟨v^{2}⟩ is the mean velocity squared of the stars, M = M_{stars} + M_{ext} is the total mass, that is, the mass in stars plus the mass of the external potential M_{ext}, and R_{v} is the virial radius of the cluster, which in this case is the same virial radius for the external potential. Now we write the total mass as a function of q:
We begin now with a derivation for the change of the velocity of a star due to an encounter with another star. The force that one star feels due to the other star is (see Fig. A.1):
Fig. A.1. Trajectory of a star with mass M that passes close to another star of mass m, which causes a small deflection on the velocity vector V. 

Open with DEXTER 
where M is the mass of one star and m is the mass of the other star, and d is the impact parameter or distance of closest approach. Now we assume that the velocity change ΔV of the star with mass M in the direction parallel to V is small, and only the perpendicular component of the velocity is changed due to the perpendicular force:
Next we want to know the time at which the perpendicular velocity changes by an amount ΔV ∼ V, which is the relaxation timescale. Therefore, we integrate Eq. (A.5) so the final perpendicular velocity is:
As the star receives many deflections in different directions, we are interested in the mean value of the squares of the velocity kicks that can be found integrating all the small deflections as:
where dN is the expected number of encounters that occurs in a time t between impact parameters b and b + db for a star with typical velocity V, this is:
where n is the number density of stars, then,
After enough time, the perpendicular velocity of one star grows to its original speed, this time is the relaxation timescale that we can derive using Eq. (A.6):
The term is often written as Λ and is the ratio of the size of the system to the “strong encounter distance” b_{min} = 2Gm/V^{2}, which is the distance at which an encounter with another star would result in a 90 degree deflection. This ratio , for a system of identical N stars is found to be (Spitzer 1987), with γ = 0.4. This definition of the relaxation timescale is useful to estimate the core collapse time, which typically occurs between 15–20 halfmass relaxation times. The problem comes when we include an analytic external potential and evolve the cluster at the center of this potential. In this case, we have an extra force acting on the stars, and therefore a virialized cluster under the influence of this external potential has a larger potential energy and also a larger kinetic energy compared to a cluster in which there is not an external potential. This should increase the relaxation timesale as the velocity of the stars increases and the number and masses of stars is kept constant. We now derive the relaxation time for a cluster under the influence of an external potential, assuming that this potential follows the same mass distribution of the stars and with the same virial radius. In Recalling the condition for virial equilibrium from Eq. (A.1) and the total mass for the cluster from Eq. (A.2) we find that:
where V is the rms velocity of the stars, N is the number of stars, m is the mass of a single star, M = M_{stars} + M_{ext} is the total mass, that is, the mass of the stars plus the mass of the external potential, and and R_{v} is the virial radius of the system. When q = 0, we get the typical velocity for stars in a cluster which is in virial equilibrium:
but for a cluster in virial equilibrium and with a background potential, the typical velocity is modified as:
where V_{ext} is the rms velocity of the stars in the presence of an external potential and V is the rms velocity of stars without a background potential.
The relaxation time is often expressed as a function of the crossing time t_{cross} of the cluster, which is simply the time it takes for a star with the typical velocity V to cross the system:
then, combining Eqs. (A.7) and (A.10) the relaxation time for clusters without an external potential is:
Now, if we include an external potential, q ≠ 0, and in substituting Eqs. (A.9) for (A.10), the crossing time for a cluster in an external potential is:
Subsequently the relaxation time for clusters in an external potential must be:
Now we replace Eq. (A.9) for Eq. (A.13) and we find:
Now in recalling that from the usual definition of the relaxation time we can define the halfmass relaxation time for a cluster of N equal mass stars as (Spitzer 1987):
then, for comparison purposes, we define the halfmass relaxation timescale for clusters in an external potential as:
with and γ = 0.4 for equal mass stars.
Appendix B: Evolution of clusters in a background potential
In this section, we describe in more detail the results of our simulations that do not include stellar collisions and which are listed in Table 1. We show the evolution of the clusters with and without an external potential, and we find a delay in the overall evolution when increasing the mass of the external potential.
All of the clusters evolve toward core collapse, which we found upon visual inspection by looking for the first drop and subsequent rise in the 10% Lagrangian radius. The core collapse occurs for the clusters without a background potential at 456 t_{cross} = 19.78 t_{rh} for the cluster with N = 1000 (see Fig. B.1), and for the cluster with N = 10 000 at 2306 t_{cross} = 13.85 t_{rh} (see Fig. B.1).
Fig. B.1. Evolution of a cluster with N = 1000 stars, total mass M_{stars} = 10^{4} M_{⊙} and R_{v} = 1.0 pc. Top panel: fraction of stars that were ejected from the cluster. Bottom panel: Lagrangian radius and the vertical green line marks the moment of core collapse. The time is presented in units of the crossing time of the cluster. 

Open with DEXTER 
As the clusters evolve toward core collapse, the outer parts continually expand which leads, along with the onset of stellar ejections, to the evaporation of the clusters. We also note that the onset of ejections occur just after core collapse when binary systems form in the core and three body interactions lead to the ejection of stars. An ejection occurs if the distance from the center of mass of the cluster to a star is larger than 20 times the virial radius R_{v} and if the total energy of the star is greater than zero, that is, E_{star} = E_{kin} +E_{pot} > 0 with E_{kin} and E_{pot} being the kinetic and potential energy, respectively. These results are shown in Figs. 1 and B.1.
When we include a background potential with a low mass compared to the total mass in stars, M_{ext} = 0.1 × M_{stars} in our simulations with N = 1000, the core collapse time is 460 t_{cross}, which is similar to the core collapse time for the cluster without the external potential in terms of the crossing time of the cluster. However, this cluster show less expansion and only one star is ejected (see Fig. B.2).
Fig. B.2. Evolution of a cluster with N = 1000 stars, total mass M_{stars} = 10^{4} M_{⊙} and R_{v} = 1.0 pc in a background potential with mass M_{ext} = 10^{3} M_{⊙}. Top panel: fraction of stars that were ejected from the cluster. Bottom panel: Lagrangian radius and the vertical green line marks the moment of core collapse. The time is presented in units of the crossing time of the cluster. 

Open with DEXTER 
Fig. B.3. Evolution of a cluster with N = 1000 stars, total mass M_{stars} = 10^{4} M_{⊙}, R_{v} = 1.0 pc in a background potential with mass M_{ext} = 5 × 10^{3} M_{⊙}. Top panel: fraction of stars that were ejected from the cluster. Bottom panel: Lagrangian radius and the vertical green line marks the moment of core collapse. The time is presented in units of the crossing time of the cluster. 

Open with DEXTER 
For the case when the mass of the background potential is half the mass in stars, M_{ext} = 0.5× M_{stars} the cluster evolves toward core collapse, which now occurs at 1820 t_{cross} = 16 t_{rh} (calculated with Eq. (A.14)). The cluster again shows less expansion than the cluster without the external potential (see Fig.B.1) and only one single star is ejected.
For the case when the mass of the external potential is the same as the total mass in stars M_{ext} = M_{stars}, the cluster also evolves toward core collapse. However, this is now delayed now until 6553 t_{cross} as shown by the green vertical line in Fig. B.4.
Fig. B.4. Evolution of a cluster with N = 1000 stars, total mass M_{stars} = 10^{4} M_{⊙}, R_{v} = 1.0 pc in a background potential with mass M_{ext} = 10^{4} M_{⊙}. Top panel: fraction of stars that were ejected from the cluster. Bottom panel: Lagrangian radius and the vertical green line marks the moment of core collapse. The time is presented in units of the crossing time of the cluster. 

Open with DEXTER 
When the number of stars is N = 10 000, the behavior is essentially the same when we include a background potential; the global evolution is delayed due to the increased velocity of the stars. First, we present the case when the mass of the external potential is low M_{ext} = 0.1 × M_{stars} compared to the total mass in stars. The cluster also evolves toward core collapse, which now occurs at 3377 t_{cross} as indicated with a green vertical line in the bottom panel of Fig. B.5, which is ∼14 t_{rh} (we used Eq. (A.14) to calculate the half mass relaxation time). The mean density inside the 10% Lagrangian radius at this moment is of 4.4 × 10^{5} M_{⊙} pc^{−3}, which is higher than the mean density at core collapse for the cluster without a background potential. Moreover the cluster in general shows less expansion when including the external potential.
Fig. B.5. Evolution of a cluster with N = 10 000 stars, total mass M_{stars} = 10^{4} M_{⊙} and R_{v} = 1.0 pc under the influence of a background potential with M_{ext} = 10^{3}M_{ext}. Top panel: fraction of stars that were ejected from the cluster, which is negligible since only one star escaped the cluster. Bottom panel: Lagrangian radius and the vertical green line marks the moment of core collapse. The time is presented in units of the crossing time of the cluster. 

Open with DEXTER 
When the mass of the external potential is half the total mass in stars, that is, M_{ext} = 0.5 × M_{stars}, the core collapse is delayed until 13 519 t_{cross} as marked by the green vertical line in Fig. B.6 which is 16 t_{rh} (we used Eq. (A.14) again to calculate the half mass relaxation time) and during this moment the mean density inside the 10% Lagrangian radius is of 5.5 × 10^{5} M_{⊙} pc^{−3}. Our simulation is not long enough to see the expansion of the outer parts. however we do expect even less expansion than for the cluster with M_{ext} = 0.1 × M_{stars} and again until this time only 1 star has been ejected from the cluster, this was also found in Leigh et al. (2013a); Leigh et al. (2014).
Fig. B.6. Evolution of a cluster with N = 10 000 stars, total mass M_{stars} = 10^{4} M_{⊙} and R_{v} = 1.0 pc under the influence of a background potential with M_{ext} = 5 × 10^{3}M_{ext}. Top panel: fraction of stars that were ejected from the cluster, which is negligible since only one star escaped the cluster. Bottom panel: Lagrangian radius and the vertical green line marks the moment of core collapse. The time is presented in units of the crossing time of the cluster. 

Open with DEXTER 
All Tables
All Figures
Fig. 1. Evolution of a cluster with N = 10 000 stars with a total stellar mass of M_{stars} = 10^{4} M_{⊙} and virial radius R_{v} = 1.0 pc. Top panel: fraction of ejections from the cluster and bottom panel: Lagrangian radii of the cluster. The vertical green line marks the moment of core collapse. 

Open with DEXTER  
In the text 
Fig. 2. Evolution of a cluster with N = 10 000 stars with a total stellar mass of M_{stars} = 10^{4} M_{⊙} and a virial radius of R_{v} = 1.0 pc in the center of an external analytic potential with M_{ext} = 10^{4} M_{⊙}. Top panel: fraction of ejections from the cluster and bottom panel: Lagrangian radii of the cluster. The vertical green line marks the moment of core collapse. 

Open with DEXTER  
In the text 
Fig. 3. Core collapse time t_{cc} divided by the crossing time t_{cross} for clusters with N = 1000 (blue circles) and N = 10 000 (red circles) as a function of q = M_{ext}/M_{stars}. We also plotted two relations that scale as (1 + q)^{4} which fit our data and therefore suggest that if the core collapse time is proportional to the relaxation time of the cluster, then the relaxation time scales as (1 + q)^{4}, which is the same relation that we derive in Appendix A. 

Open with DEXTER  
In the text 
Fig. 4. Number of binary systems as a function of time and for different values for the mass of the external potential M_{ext} for simulations with N = 1000 and up to 8000 crossing times. The vertical lines mark the moments of core collapse and coincides with the onset of the formation of binary systems. 

Open with DEXTER  
In the text 
Fig. 5. Distribution of the semimajor axis of the binary population for the set of simulations listed in Table 1. The data clearly shows that increasing the mass of the external potential the formed binaries tend to be more tightly bound. Despite the fact that a low mass of the external potential (M_{ext} < 1.0 × M_{stars}) may promote the formation of wider binaries, these systems tend to be shortlived due to the higher velocities of all the stars. This figure is only intended to show the distribution of the semimajor axis and not the total number of binaries that formed. 

Open with DEXTER  
In the text 
Fig. 6. Evolution of a cluster with N = 1000 stars, total mass M_{stars} = 10^{4} M_{⊙}, R_{v} = 0.14 pc, and R_{star} = 100 R_{⊙}. Uppermost panel: number of mergers as a function of time, the fraction of binary, hyperbolic, and chain mergers in the second panel along with the fraction of stars that escape from the cluster. Third panel: 10, 50 and 90% Lagrangian radii, and bottom panel: mass growth of the most massive object that formed through mergers in the cluster. 

Open with DEXTER  
In the text 
Fig. 7. Evolution of a cluster with N = 1000 stars, stellar mass M_{stars} = 10^{4} M_{⊙}, R_{v} = 0.14 pc, and R_{star} = 100 R_{⊙} in the center of an external potential with M_{ext} = 10^{4} M_{⊙}. Uppermost panel: number of mergers as a function of time, the fraction of binary, hyperbolic, and chain mergers in the second panel, along with the fraction of stars that escape from the cluster. Third panel: 10, 50, and 90% Lagrangian radii, and bottom panel: mass growth of the most massive object that formed through mergers in the cluster. The drop in the 50% Lagrangian radius toward the end of the simulation suggests that the central object may grow up to 5000 M_{⊙}. 

Open with DEXTER  
In the text 
Fig. 8. Example of a Gaussian fit to the number of mergers with the central star for simulation number 45 listed in Table 2. Top panel: number of mergers with the most massive star as a function of time (black solid line) along with the Gaussian function that we used to estimate the total number of mergers at different times. The parameters t_{delay} and t_{duration} are also shown in the figure. Bottom panel: mass growth of the central star whose rapid growth coincides with the peak in the number of mergers. 

Open with DEXTER  
In the text 
Fig. 9. Mean mass of the most massive star normalized by its initial mass, as a function of N and R_{star} for clusters without an external potential evolved until 2000 t_{cross} = 15.6 Myr. The symbols represent the mean values obtained from six simulations and the solid lines show the results of the calculations with our model presented in Eq. (8) using t = 15.6 Myr and the parameters obtained from Eqs. (9)–(12). 

Open with DEXTER  
In the text 
Fig. 10. Mean mass of the most massive star normalized by its initial mass, as a function of N and R_{star} for clusters without an external potential evolved until 3000 t_{cross} = 16.5 Myr, except for the point at N = 5000 and R_{star} = 20 R_{⊙}, which are simulations that were evolved until 4000 t_{cross} = 22 Myr. The symbols represent the mean values obtained from six simulations and the solid lines show the results of the calculations with our model presented in Eq. (8) using t = 16.5 Myr (t = 22 Myr for N = 5000 and R_{star} < 50 R_{⊙}) and the parameters obtained from Eqs. (13)–(16). 

Open with DEXTER  
In the text 
Fig. 11. Parameters of the Gaussian fit to the number of mergers with the central star for simulations of clusters without an external potential. The data points were obtained after the fit to the combined data of six simulations. Upper panel: values of the normalization factor A for different values of N and R_{star} along with our model from Eq. (9). These data points show that the number of mergers increases with both N and R_{star} as expected. Bottom panel: different values of t_{delay} in addition to the different values of N and R_{star} along with the fitting function presented in Eq. (10). These points show that t_{delay} decreases with R_{star}, but it increases with N. 

Open with DEXTER  
In the text 
Fig. 12. t_{duration} as a function of the number of stars N for clusters without an external potential. We found no clear correlations between t_{duration} and R_{star}. The solid blue line shows the fitting function from Eq. (11). 

Open with DEXTER  
In the text 
Fig. 13. Parameters of the Gaussian fit to the number of mergers with the central star for simulations of clusters with an external potential. The data points were obtained after the fit to the combined data of six simulations. Upper panel: values of the normalization factor A_{ext} for different values of N and R_{star} along with our model from Eq. (8). These data points show that the number of mergers increases with both N and R_{star} as expected. Middle panel: different values of t_{delay, ext} in addition to the different values of N and R_{star} along with the fitting function presented in Eq. (14). These points show that t_{delay, ext} decreases with R_{star}, but it increases with N. Additionally, these values are larger than the values of t_{delay} for clusters without an external potential (see Fig. 11). Lower panel: values of t_{duration, ext} along with the fitting function from Eq. (15). We clearly see a positive correlation with N and a negative correlation with R_{star}. 

Open with DEXTER  
In the text 
Fig. 14. Fraction of stars f_{merger} that merge with the central star in a cluster with a total stellar mass of M_{stars} = 10^{4} M_{⊙} and a virial radius of R_{v} = 0.14 pc after 1 Myr (left panel), 2 Myr (middle panel), and 10 Myr (right panel) for different combinations of the number of stars N and the stellar radius R_{star} calculated using our model presented in Eq. (8). We can see that when the collisions are maintained for short periods of time, that is, 1 Myr, a massive object with ∼f_{merger}M_{stars} = 0.15 × 10^{4} M_{⊙} = 1500 M_{⊙} is formed in the cluster center only when N is small (given that t_{delay} becomes larger than the integration time for large N) or large R_{star}. If the time limit is long, that is, 10 Myr, most of the clusters will form stars with ∼0.3 M_{stars} = 3000 M_{⊙}. 

Open with DEXTER  
In the text 
Fig. 15. Fraction of stars f_{merger} that merge with the central star in a cluster with a total stellar mass of M_{stars} = 10^{4} M_{⊙}, a virial radius of R_{v} = 0.14 pc, and an external potential with a mass of M_{ext} = M_{stars} after 1 Myr (left panel), 2 Myr (middle panel), and 10 Myr (right panel) for different combinations of the number of stars N and the stellar radius R_{star} calculated using our model presented in Eq. (8). We can see that when the collisions are maintained for short periods of time, that is, 1 Myr, an object that is not so massive with ∼f_{merger}M_{stars} = 0.1 × 10^{4} M_{⊙} = 1000 M_{⊙} is formed in the cluster center only when N is very small (given that t_{delay} becomes larger than the integration time for large N) or when R_{star} is large. If the time limit is 2 Myr (middle panel), we expect that clusters with a small N experience core collapse and form a massive central star. However, if the number of stars is large and R_{star} is also large, given the higher velocity of the stars in those clusters compared to clusters without the external potential, there is a larger probability for hyperbolic mergers to occur before core collapse, which we see on the righthand side of middle panel. If the time limit is long, that is, 10 Myr, most of the clusters form stars with ∼0.3−0.4 M_{stars} = 3000−4000 M_{⊙} and even 0.7 M_{stars} = 7000 M_{⊙} if N is large enough. 

Open with DEXTER  
In the text 
Fig. A.1. Trajectory of a star with mass M that passes close to another star of mass m, which causes a small deflection on the velocity vector V. 

Open with DEXTER  
In the text 
Fig. B.1. Evolution of a cluster with N = 1000 stars, total mass M_{stars} = 10^{4} M_{⊙} and R_{v} = 1.0 pc. Top panel: fraction of stars that were ejected from the cluster. Bottom panel: Lagrangian radius and the vertical green line marks the moment of core collapse. The time is presented in units of the crossing time of the cluster. 

Open with DEXTER  
In the text 
Fig. B.2. Evolution of a cluster with N = 1000 stars, total mass M_{stars} = 10^{4} M_{⊙} and R_{v} = 1.0 pc in a background potential with mass M_{ext} = 10^{3} M_{⊙}. Top panel: fraction of stars that were ejected from the cluster. Bottom panel: Lagrangian radius and the vertical green line marks the moment of core collapse. The time is presented in units of the crossing time of the cluster. 

Open with DEXTER  
In the text 
Fig. B.3. Evolution of a cluster with N = 1000 stars, total mass M_{stars} = 10^{4} M_{⊙}, R_{v} = 1.0 pc in a background potential with mass M_{ext} = 5 × 10^{3} M_{⊙}. Top panel: fraction of stars that were ejected from the cluster. Bottom panel: Lagrangian radius and the vertical green line marks the moment of core collapse. The time is presented in units of the crossing time of the cluster. 

Open with DEXTER  
In the text 
Fig. B.4. Evolution of a cluster with N = 1000 stars, total mass M_{stars} = 10^{4} M_{⊙}, R_{v} = 1.0 pc in a background potential with mass M_{ext} = 10^{4} M_{⊙}. Top panel: fraction of stars that were ejected from the cluster. Bottom panel: Lagrangian radius and the vertical green line marks the moment of core collapse. The time is presented in units of the crossing time of the cluster. 

Open with DEXTER  
In the text 
Fig. B.5. Evolution of a cluster with N = 10 000 stars, total mass M_{stars} = 10^{4} M_{⊙} and R_{v} = 1.0 pc under the influence of a background potential with M_{ext} = 10^{3}M_{ext}. Top panel: fraction of stars that were ejected from the cluster, which is negligible since only one star escaped the cluster. Bottom panel: Lagrangian radius and the vertical green line marks the moment of core collapse. The time is presented in units of the crossing time of the cluster. 

Open with DEXTER  
In the text 
Fig. B.6. Evolution of a cluster with N = 10 000 stars, total mass M_{stars} = 10^{4} M_{⊙} and R_{v} = 1.0 pc under the influence of a background potential with M_{ext} = 5 × 10^{3}M_{ext}. Top panel: fraction of stars that were ejected from the cluster, which is negligible since only one star escaped the cluster. Bottom panel: Lagrangian radius and the vertical green line marks the moment of core collapse. The time is presented in units of the crossing time of the cluster. 

Open with DEXTER  
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.