Abstract Recent observations indicate that the progenitors of globular clusters (GCs) at high redshifts had high average stellar surface densities above 105 M⊙pc−2. The internal structure and kinematics of the clusters, however, remain out of reach. Numerical simulations are necessary to decipher the origin of spatio-kinematic features in present-day GCs. Here we study star cluster formation in a star-by-star hydrodynamical simulation of a low-metallicity starburst in a merger of two gas-rich dwarf galaxies. The simulation accounts for the multiphase interstellar medium, stellar radiation, winds and supernovae, and the accurate small-scale gravitational dynamics near massive stars. We also include prescriptions for stellar collisions and tidal disruption events by black holes. Gravitationally bound star clusters up to ∼2 × 105 M⊙ form dense with initial half-mass radii of ∼0.1-1 pc. The most massive cluster approaches the observed high-redshift surface densities throughout its hierarchical and dissipative assembly. The cluster also hosts a collisionally growing very massive star of ∼1000 M⊙ that will eventually collapse, forming an intermediate mass black hole. The assembly leaves an imprint in the spatio-kinematic structure of the cluster. The youngest stars are more centrally concentrated, they show significant bulk rotation and have radially biased velocity components at outer radii. The older population is more round in shape, rotates slowly, its velocity distribution is isotropic and exhibits higher dispersion. If chemically enriched star formation proceeds mainly in the later stages of cluster assembly, these results provide a possible explanation for some of the multiple population features observed in dynamically young GCs.
Lahén et al. (Sat,) studied this question.