Direct N-body simulations of rotating and extremely massive Population III star clusters

TL;DR

This study uses detailed N-body simulations to explore the evolution, black hole formation, and dynamics of extremely massive, rotating Population III star clusters, revealing key processes in early universe black hole seed formation.

Contribution

It introduces comprehensive simulations including rotation, binaries, and relativistic effects, demonstrating new pathways for intermediate-mass black hole formation in Population III clusters.

Findings

Faster rotation leads to earlier core collapse and more black hole mergers.

Clusters form black holes across the pair-instability gap with multi-generation growth.

Initial rotation influences cluster evolution, mass segregation, and IMBH subsystem formation.

Abstract

Aims. We present eight direct N-body simulations with NBODY6++GPU of extremely massive, initially rotating Population III star clusters with 1.01 x 10^5 stars. Methods. Our models include primordial binaries, a continuous initial mass function, differential rotation, tidal mass loss, updated fitting formulae for extremely massive metal-poor Population III stars, and general-relativistic merger recoil kicks. We assess their impact on cluster dynamics. Results. All runs form black holes below, within, and above the pair-instability gap, with multi-generation growth. Faster-rotating clusters core-collapse earlier; post-collapse clusters host a rotating, axisymmetric subsystem of intermediate-mass black holes (IMBHs) at the centre and an expanding halo of lower-mass objects. Pair-instability supernovae and compact-object formation at ~2-3 Myr sharply reduce total mass and a large fraction of the cluster's angular momentum. All Population III clusters in our simulations have the gravothermal-gravogyro catastrophe phase. Conclusions. We confirm two of the hypothesized formation channels of galactic nucleus seed black holes: gravitational runaway mergers of black holes and of Population III stars, which core-collapse into IMBHs thereafter. A higher initial star cluster bulk rotation correlates with earlier core collapse and, in the event counts reported here, with more coalescences and collisions, as well as lower retained (compact) binary abundances. Initial bulk rotation is a primary control parameter of cluster evolution: faster rotation accelerates early angular-momentum transport, gravothermal collapse, mass segregation, and amplifies post-collapse expansion, which also favours the formation of a compact central IMBH subsystem.