Collisions were suggested to potentially play a role in the formation of massive stars in present day clusters, and have likely been relevant during the formation of massive stars and intermediate mass black holes within the first star clusters. In the early Universe, the first stellar clusters were particularly dense, as fragmentation typically only occurred at densities above 109 cm−3, and the radii of the protostars were enhanced as a result of larger accretion rates, suggesting a potentially more relevant role of stellar collisions. We present here a detailed parameter study to assess how the number of collisions and the mass growth of the most massive object depend on the properties of the cluster. We also characterize the time evolution with three effective parameters: the time when most collisions occur, the duration of the collisions period, and the normalization required to obtain the total number of collisions. We apply our results to typical Population III (Pop. III) clusters of about 1000 M⊙, finding that a moderate enhancement of the mass of the most massive star by a factor of a few can be expected. For more massive Pop. III clusters as expected in the first atomic cooling halos, we expect a more significant enhancement by a factor of 15–32. We therefore conclude that collisions in massive Pop. III clusters were likely relevant to form the first intermediate mass black holes.

Models aiming to explain the formation of massive black hole seeds, and in particular the direct collapse scenario, face substantial difficulties. These are rooted in rather ad hoc and fine-tuned initial conditions, such as the simultaneous requirements of extremely low metallicities and strong radiation backgrounds. Here we explore a modification of such scenarios where a massive primordial star cluster is initially produced. Subsequent stellar collisions give rise to the formation of massive (10 4 -10 5 M ) objects. Our calculations demonstrate that the interplay between stellar dynamics, gas accretion and protostellar evolution is particularly relevant. Gas accretion onto the protostars enhances their radii, resulting in an enhanced collisional cross section. We show that the fraction of collisions can increase from 0.1-1% of the initial population to about 10% when compared to gas-free models or models of protostellar clusters in the local Universe. We conclude that very massive objects can form in spite of initial fragmentation, making the first massive protostellar clusters viable candidate birth places for observed supermassive black holes.

In this paper, we compute predictions for the number of stellar collisions derived from analytic models based on the mean free path (MFP) approximation and compare them to the results of N-body simulations. Our goal is to identify the cluster conditions under which the MFP approximation remains valid. Adopting a range of particle numbers (100 ≤ N ≤ 5000) and different combinations of particle masses and radii, we explore three different channels leading to stellar collisions, all of which are expected to occur in realistic stellar environments. At high densities, binaries form from isolated three-body interactions of single stars. Hence, we consider collisions between single stars and collisions involving binary stars, after they form in our simulations. For the latter, we consider two channels for mergers, namely direct stellar collisions during chaotic single–binary interactions and perturbation-driven mergers of binaries due to random walks in eccentricity approaching unity. In the densest systems considered here, a very massive object is formed at the cluster centre, causing local stellar orbits to become increasingly Keplerian and the assumptions going into our analytic model to break down. Before reaching this limit, we obtain excellent agreement between our theoretical predictions and the simulations: The analytic rates are typically accurate to within one standard deviation for the entire parameter space considered here, but the agreement is best for short integration times. Our results have direct implications for blue straggler formation in dense star clusters, and stellar mergers in galactic nuclei hosting massive black holes.

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 time-consuming 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 N-body 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.