Massimo Ricotti

We present the results of a set of radiation magnetohydrodynamic simulations of turbulent molecular clouds in which we vary the initial strength of the magnetic field within a range (1 ≲ μ ≲ 5) consistent with observations of local giant molecular clouds (GMCs). We find that as we increase the strength of the magnetic field, star formation transitions from unimodal (the baseline case, μ = 5, with a single burst of star formation and Salpeter IMF) to bimodal. This effect is clearest in the most strongly magnetized GMC (μ = 1): a first burst of star formation with duration, intensity, and IMF comparable to the baseline case, is followed by a second star formation episode in which only low-mass stars are formed. Overall, due to the second burst of star formation, the strongly magnetized case results in a longer star formation period and a higher efficiency of star formation. The second burst is produced by gas that is not expelled by radiative feedback, instead remains trapped in the GMC by the large-scale B-field, producing a nearly one-dimensional flow of gas along the field lines. The trapped gas has a turbulent and magnetic topology that differs from that of the first phase and strongly suppresses gas accretion onto protostellar cores, reducing their masses. We speculate that this star formation bimodality may be an important ingredient to understand the origin of multiple stellar populations observed in massive globular clusters.

Using radiation-hydrodynamic cosmological simulations, we present a detailed (0.1 pc resolution), physically motivated portrait of a typical-mass dwarf galaxy before the epoch of reionization, resolving the formation, and evolution of star clusters into individual 10 M⊙ star particles. In the rest-frame ultraviolet, the galaxy has an irregular morphology with no bulge or disc, dominated by light emitted from numerous, compact, and gravitationally-bound star clusters. This is especially interesting in light of recent James Webb Space Telescope observations that − aided by the magnifying power of gravitational lenses – have imaged, at parsec-scale resolution, individual young star clusters forming in similar galaxies at z> 6. Because of their low metallicities and high temperatures, star-forming gas clouds in this galaxy have densities ∼100 times higher than typical giant molecular clouds; hence, their expected star formation efficiencies (SFEs) are high enough (around 10 − 70 per cent) to produce a sizeable population of potential globular cluster progenitors, but typically smaller (a few 100 − 2 ×104M⊙, half-mass radii of up to 3 pc) and of lower metallicities (10−3.5– 10−2.5 Z⊙). The initial mass function of the star-forming clouds is log-normal, whereas the bound star cluster mass function is a power-law with a slope that depends mainly on SFE but also on the temporal proximity to a major starburst. We find slopes between −0.5 and −2.5 depending on the assumed sub-grid SFE. Star formation is self-regulated on galactic scales; however, the multimodal metallicity distribution of the star clusters and the fraction of stars locked into surviving bound star clusters depends on SFE.

We model Pop III star formation in different FUV and X-ray backgrounds, including radiation feedback from protostars. We confirm previous results that a moderate X-ray background increases the number of Pop III systems per unit cosmological volume, but masses and multiplicities of the system are reduced. The stellar mass function also agrees with previous results, and we confirm the outward migration of the stars within the protostellar discs. We find that nearly all Pop III star systems are hierarchical, i.e. binaries of binaries. Typically, two equal-mass stars form near the centre of the protostellar disc and migrate outwards. Around these stars, mini-discs fragment forming binaries that also migrate outwards. Stars may also form at Lagrange points L4/L5 of the system. Afterwards, star formation becomes more stochastic due to the large multiplicity, and zero-metallicity low-mass stars can form when rapidly ejected from the disc. Stars in the disc often have eccentric orbits, leading to a periodic modulation of their accretion rates and luminosities. At the pericentre, due to strong accretion, the star can enter a red-supergiant phase reaching nearly Eddington luminosity in the optical bands (mAB ∼ 34 for a 100 M⊙ star at z = 6). During this phase, the star, rather than its nebular lines, can be observed directly by JWST, if sufficiently magnified by a gravitational lens. The ∼10 000 AU separations and high eccentricities of many Pop III star binaries in our simulations are favourable parameters for IMBH mergers – and gravitational waves emission – through orbital excitation by field stars.