No Arabic abstract
Pristine, atomically-cooled haloes are leading contenders for the sites of primordial quasar formation because atomic cooling triggers rapid baryon collapse that can create 10$^4$ - 10$^5$ M$_{odot}$ black hole seeds. However, until now no numerical simulations with a wide range of halo spins and assembly histories have followed the collapse for the times required to form a black hole. We have now performed cosmological simulations of baryon collapse in atomically-cooled haloes for times that are sufficient for supermassive stars to form and die as direct-collapse black holes (DCBHs). Our simulations reveal that fragmentation of the accretion disk at the center of the halo after $sim$ 500 kyr is nearly ubiquitous and in most cases leads to the formation of binary or multiple supermassive stellar systems. They also confirm that rapid baryon collapse proceeds for the times required for these stars to form DCBHs. Our discovery raises the exciting possibility of detecting gravitational waves from DCBH mergers with LISA and tidal disruption events in the near infrared with the James Webb Space Telescope and ground-based telescopes in the coming decade.
Population III stars can regulate star formation in the primordial Universe in several ways. They can ionize nearby halos, and even if their ionizing photons are trapped by their own halos, their Lyman-Werner (LW) photons can still escape and destroy H$_2$ in other halos, preventing them from cooling and forming stars. LW escape fractions are thus a key parameter in cosmological simulations of early reionization and star formation but have not yet been parametrized for realistic halos by halo or stellar mass. To do so, we perform radiation hydrodynamical simulations of LW UV escape from 9--120 M$_{odot}$ Pop III stars in $10^5$ to $10^7$ M$_{odot}$ halos with ZEUS-MP. We find that photons in the LW lines (i.e. those responsible for destroying H$_{2}$ in nearby systems) have escape fractions ranging from 0% to 85%. No LW photons escape the most massive halo in our sample, even from the most massive star. Escape fractions for photons elsewhere in the 11.18--13.6~eV energy range, which can be redshifted into the LW lines at cosmological distances, are generally much higher, being above 60% for all but the least massive stars in the most massive halos. We find that shielding of H$_2$ by neutral hydrogen, which has been neglected in most studies to date, produces escape fractions that are up to a factor of three smaller than those predicted by H$_2$ self-shielding alone.
Supermassive stars born in pristine environments in the early Universe hold the promise of being the seeds for the supermassive black holes observed as high redshift quasars shortly after the epoch of reionisation. H$_2$ suppression is thought to be crucial in order to negate normal Population III star formation and allow high accretion rates to drive the formation of supermassive stars. Only in the cases where vigorous fragmentation is avoided will a monolithic collapse be successful giving rise to a single massive central object. We investigate the number of fragmentation sites formed in collapsing atomic cooling haloes subject to various levels of background Lyman-Werner flux. The background Lyman-Werner flux manipulates the chemical properties of the gas in the collapsing halo by destroying H$_2$. We find that only when the collapsing gas cloud shifts from the molecular to the atomic cooling regime is the degree of fragmentation suppressed. In our particular case we find that this occurs above a critical Lyman-Werner background of J $sim 10$ J$_{21}$. The important criterion being the transition to the atomic cooling regime rather than the actual value of J, which will vary locally. Once the temperature of the gas exceeds T $gtrsim$ 10$^4$ K and the gas transitions to atomic line cooling, then vigorous fragmentation is strongly suppressed.
We investigate the process of metal-free star formation in the first galaxies with a high-resolution cosmological simulation. We consider the cosmologically motivated scenario in which a strong molecule-destroying Lyman-Werner (LW) background inhibits effective cooling in low-mass haloes, delaying star formation until the collapse or more massive haloes. Only when molecular hydrogen (H2) can self-shield from LW radiation, which requires a halo capable of cooling by atomic line emission, will star formation be possible. To follow the formation of multiple gravitationally bound objects, at high gas densities we introduce sink particles which accrete gas directly from the computational grid. We find that in a 1 Mpc^3 (comoving) box, runaway collapse first occurs in a 3x10^7 M_sun dark matter halo at z~12 assuming a background intensity of J21=100. Due to a runaway increase in the H2 abundance and cooling rate, a self-shielding, supersonically turbulent core develops abruptly with ~10^4 M_sun in cold gas available for star formation. We analyze the formation of this self-shielding core, the character of turbulence, and the prospects for star formation. Due to a lack of fragmentation on scales we resolve, we argue that LW-delayed metal-free star formation in atomic cooling haloes is very similar to star formation in primordial minihaloes, although in making this conclusion we ignore internal stellar feedback. Finally, we briefly discuss the detectability of metal-free stellar clusters with the James Webb Space Telescope.
Direct collapse models for black hole (BH) formation predict massive ($sim 10^5 M_{odot}$) seeds, which hold great appeal as a means to rapidly grow the observed $sim 10^9 M_{odot}$ quasars by $zgtrsim 7$; however, their formation requires fine-tuned conditions. In this work, we use cosmological zoom simulations to study systematically the impact of requiring: 1) low gas angular momentum, and 2) a minimum incident Lyman Werner (LW) flux radiation in order to form direct-collapse BH seeds. We start with a baseline model (introduced in Bhowmick et. al. 2021) that restricts black hole seed formation (with seed masses of $M_{mathrm{seed}}=1.25times10^{4},1times10^{5} & 8times10^{5}M_{odot}/h$) to occur only in haloes with a minimum total mass ($3000times M_{mathrm{seed}}$) and star forming, metal poor gas mass ($5times M_{mathrm{seed}}$). When seeding is further restricted to halos with low gas spins (i.e. smaller than the minimum value required for the gas disc to be gravitationally stable), the seeding frequency is suppressed by factors of $sim6$ compared to the baseline model regardless of the mass threshold used. In contrast, imposing a minimum LW flux ($>10J_{21}$) disproportionately suppresses seed formation in $lesssim10^9M_{odot}/h$ halos, by factors of $sim100$ compared to the baseline model. Very few BH merger events occur in the models with a LW flux criterion, and because early BH growth is dominated by mergers in our models, this results in only the most massive ($8times10^{5}M_{odot}/h$) seeds being able to grow to the supermassive regime ($gtrsim10^6M_{odot}/h$) by $z=7$. Our results therefore suggest that producing the bulk of the $zgtrsim7$ BH population requires alternate seeding channels, early BH growth dominated by rapid or super-eddington accretion, massive seeding scenarios that do not depend on LW flux, or a combination of these possibilities.
Direct collapse black holes forming in pristine, atomically-cooling haloes at $z approx 10-20$ may act as the seeds of supermassive black holes (BH) at high redshifts. In order to create a massive BH seed, the host halo needs to be prevented from forming stars. H$_2$ therefore needs to be irradiated by a large flux of Lyman-Werner (LW) UV photons in order to suppress H$_2$ cooling. A key uncertainty in this scenario is the escape fraction of LW radiation from first galaxies, the dominant source of UV photons at this epoch. To better constrain this escape fraction, we have performed radiation-hydrodynamical simulations of the growth of HII regions and their associated photodissociation regions in the first galaxies using the ZEUS-MP code. We find that the LW escape fraction crucially depends on the propagation of the ionisation front (I-front). For an R-type I-front overrunning the halo, the LW escape fraction is always larger than 95%. If the halo recombines later from the outside--in, due to a softened and weakened spectrum, the LW escape fraction in the rest-frame of the halo (the near-field) drops to zero. A detailed and careful analysis is required to analyse slowly moving, D-type I-fronts, where the escape fraction depends on the microphysics and can be as small as 3% in the near-field and 61% in the far-field or as large as 100% in both the near-field and the far-field.