No Arabic abstract
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 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.
The first stars in the Universe, the so-called Population III stars, form in small dark matter minihaloes with virial temperatures $T_{rm vir} < 10^{4}$~K. Cooling in these minihaloes is dominated by molecular hydrogen (H$_{2}$), and so Population III star formation is only possible in those minihaloes that form enough H$_{2}$ to cool on a short timescale. As H$_{2}$ cooling is more effective in more massive minihaloes, there is therefore a critical halo mass scale $M_{rm min}$ above which Population III star formation first becomes possible. Two important processes can alter this minimum mass scale: streaming of baryons relative to the dark matter and the photodissociation of H$_{2}$ by a high redshift Lyman-Werner (LW) background. In this paper, we present results from a set of high resolution cosmological simulations that examine the impact of these processes on $M_{rm min}$ and on $M_{rm ave}$ (the average minihalo mass for star formation), both individually and in combination. We show that streaming has a bigger impact on $M_{rm min}$ than the LW background, but also that both effects are additive. We also provide fitting functions quantifying the dependence of $M_{rm ave}$ and $M_{rm min}$ on the streaming velocity and the strength of the LW background.
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.
We propose that star formation is delayed relative to the inflow rate in rapidly-accreting galaxies at very high redshift (z > 2) because of the energy conveyed by the accreting gas. Accreting gas streams provide fuel for star formation, but they stir the disk and increase turbulence above the usual levels compatible with gravitational instability, reducing the star formation efficiency in the available gas. After the specific inflow rate has sufficiently decreased - typically at z < 3 - galaxies settle in a self-regulated regime with efficient star formation. An analytic model shows that this interaction between infalling gas and young galaxies can significantly delay star formation and maintain high gas fractions (>40%) down to z = 2, in contrast to other galaxy formation models. Idealized hydrodynamic simulations of infalling gas streams onto primordial galaxies confirm the efficient energetic coupling at z > 2, and suggest that this effect is largely under-resolved in existing cosmological simulations.
Radiative transfer plays a major role in the process of star formation. Many simulations of gravitational collapse of a cold gas cloud followed by the formation of a protostellar core use a grey treatment of radiative transfer coupled to the hydrodynamics. However, dust opacities which dominate extinction show large variations as a function of frequency. In this paper, we used frequency-dependent radiative transfer to investigate the influence of the opacity variations on the properties of Larsons first core. We used a multigroup M1 moment model in a 1D radiation hydrodynamics code to simulate the spherically symmetric collapse of a 1 solar mass cloud core. Monochromatic dust opacities for five different temperature ranges were used to compute Planck and Rosseland means inside each frequency group. The results are very consistent with previous studies and only small differences were observed between the grey and multigroup simulations. For a same central density, the multigroup simulations tend to produce first cores with a slightly higher radius and central temperature. We also performed simulations of the collapse of a 10 and 0.1 solar mass cloud, which showed the properties of the first core to be independent of the initial cloud mass, with again no major differences between grey and multigroup models. For Larsons first collapse, where temperatures remain below 2000 K, the vast majority of the radiation energy lies in the IR regime and the system is optically thick. In this regime, the grey approximation does a good job reproducing the correct opacities, as long as there are no large opacity variations on scales much smaller than the width of the Planck function. The multigroup method is however expected to yield more important differences in the later stages of the collapse when high energy (UV and X-ray) radiation is present and matter and radiation are strongly decoupled.