We use the moving-mesh code arepo to investigate the effects of different levels of rotation and turbulence on the fragmentation of primordial gas and the formation of Population III stars. We consider 9 different combinations of turbulence and rotation and carry out 5 different realizations of each setup, yielding one of the largest sets of simulations of Population III star formation ever performed. We find that fragmentation in Population III star-forming systems is a highly chaotic process and show that the outcomes of individual realizations of the same initial conditions often vary significantly. However, some general trends are apparent. Increasing the turbulent energy promotes fragmentation, while increasing the rotational energy inhibits fragmentation. Within the ∼ 1000 yr period that we simulate, runs including turbulence yield flat protostellar mass functions while purely rotational runs show a more top-heavy distribution. The masses of the individual protostars are distributed over a wide range from a few 10 −3 M to several tens of M . The total mass growth rate of the stellar systems remains high throughout the simulations and depends only weakly on the degree of rotation and turbulence. Mergers between protostars are common, but predictions of the merger fraction are highly sensitive to the criterion used to decide whether two protostars should merge. Previous studies of Population III star formation have often considered only one realization per set of initial conditions. However, our results demonstrate that robust trends can only be reliably identified by considering averages over a larger sample of runs.

Metals from Population III (Pop III) supernovae led to the formation of less massive Pop II stars in the early universe, altering the course of evolution of primeval galaxies and cosmological reionization. There are a variety of scenarios in which heavy elements from the first supernovae were taken up into second-generation stars, but cosmological simulations only model them on the largest scales. We present small-scale, high-resolution simulations of the chemical enrichment of a primordial halo by a nearby supernova after partial evaporation by the progenitor star. We find that ejecta from the explosion crash into and mix violently with ablative flows driven off the halo by the star, creating dense, enriched clumps capable of collapsing into Pop II stars. Metals may mix less efficiently with the partially exposed core of the halo, so it might form either Pop III or Pop II stars. Both Pop II and III stars may thus form after the collision if the ejecta do not strip all the gas from the halo. The partial evaporation of the halo prior to the explosion is crucial to its later enrichment by the supernova.

Radiative feedback from massive Population III (Pop. III) stars in the form of ionising and photodissociating photons is widely believed to play a central role in shutting off accretion on to these stars. Understanding whether and how this occurs is vital for predicting the final masses reached by these stars and the form of the Pop. III stellar initial mass function. To help us better understand the impact of UV radiation from massive Pop. III stars on the gas surrounding them, we carry out high resolution simulations of the formation and early evolution of these stars, using the AREPO moving-mesh code coupled with the innovative radiative transfer module SPRAI. Contrary to most previous results, we find that the ionising radiation from these stars is trapped in the dense accretion disk surrounding them. Consequently, the inclusion of radiative feedback has no significant impact on either the number or the total mass of protostars formed during the 20 kyr period that we simulate. We show that the reason that we obtain qualitatively different results from previous studies of Pop. III stellar feedback lies in how the radiation is injected into the simulation. HII region trapping only occurs if the photons are injected on scales smaller than the local scale height of the accretion disk, a criterion not fulfilled in previous 3D simulations of this process. Finally, we speculate as to whether outflows driven by the magnetic field or by Lyman-α radiation pressure may be able to clear enough gas away from the star to allow the HII region to escape from the disk.