Globular cluster formation from inertial inflows: accreting extremely massive stars as the origin of abundance anomalies

Abstract: We use the inertial-inflow model of massive star formation to describe the formation of globular clusters (GCs) in turbulent molecular clouds. A key aspect of this model is that the maximum stellar mass scales linearly with cloud mass, such that extremely massive stars (EMSs, $10^{3-4}\,{\rm M}\odot$) form in massive GCs ($\gtrsim10^5\,M\odot$). The total wind mass loss is dominated by accreting EMSs (aEMSs), whose wind mass-loss rates become comparable to their accretion rates ($\gtrsim10^{-2}\,{\rm M}\odot\,{\rm yr}^{-1}$). These winds pollute the intra-cluster medium with hot-hydrogen burning yields during GC formation. We propose a parameterised model for the evolution of the stellar mass function during GC formation ($\sim 1-2\,{\rm Myr}$), accounting for gas inflow, wind mass loss and mixing of aEMS yields with pristine gas that has initial proto-GC abundances. Low-mass stars ($\lesssim1\,{\rm M}\odot$) form continuously from this mixed gas and their abundances resemble observed abundance trends with GC mass and metallicity, specifically: (i) the helium spread in a typical GC is small ($\Delta Y \simeq 0.01$) and increases with GC mass; (ii) the fraction of polluted stars increases with GC mass and metallicity; (iii) the extent of the Mg-Al anticorrelations is more pronounced in metal-poor and massive GCs. We conclude that GCs formed with a population of EMSs, that may leave behind black holes with masses above the pair-instability gap ($\gtrsim120\,{\rm M}_\odot$) and that the nitrogen-rich galaxies found by the James Webb Space Telescope ({\it JWST}) are in the first phase of galaxy evolution, when a large fraction of star formation occurs in the form of GCs.