No Arabic abstract
We have modeled direct collapse of a primordial gas within dark matter halos in the presence of radiative transfer, in high-resolution zoom-in simulations in a cosmological framework, down to the formation of the photosphere and the central object. Radiative transfer has been implemented in the flux-limited diffusion (FLD) approximation. Adiabatic models were run for comparison. We find that (a) the FLD flow forms an irregular central structure and does not exhibit fragmentation, contrary to adiabatic flow which forms a thick disk, driving a pair of spiral shocks, subject to Kelvin-Helmholtz shear instability forming fragments; (b) the growing central core in the FLD flow quickly reaches ~10 Mo and a highly variable luminosity of 10^{38}-10^{39} erg/s, comparable to the Eddington luminosity. It experiences massive recurrent outflows driven by radiation force and thermal pressure gradients, which mix with the accretion flow and transfer the angular momentum outwards; and (c) the interplay between these processes and a massive accretion, results in photosphere at ~10 AU. We conclude that in the FLD model (1) the central object exhibits dynamically insignificant rotation and slower than adiabatic temperature rise with density; (2) does not experience fragmentation leading to star formation, thus promoting the fast track formation of a supermassive black hole (SMBH) seed; (3) inclusion of radiation force leads to outflows, resulting in the mass accumulation within the central 10^{-3} pc, which is ~100 times larger than characteristic scale of star formation. The inclusion of radiative transfer reveals complex early stages of formation and growth of the central structure in the direct collapse scenario of SMBH seed formation.
Direct collapse within dark matter (DM) halos is a promising path to form supermassive black hole (SMBH) seeds at high redshifts. The outer part of this collapse remains optically thin, and has been studied intensively using numerical simulations. However, the innermost region of the collapse is expected to become optically thick and requires us to follow the radiation field in order to understand its subsequent evolution. So far, the adiabatic approximation has been used exclusively for this purpose. We apply radiative transfer in the flux-limited diffusion (FLD) approximation to solve the evolution of coupled gas and radiation, for isolated halos. For direct collapse within isolated DM halos, we find that (1) the photosphere forms at ~10^{-6} pc and rapidly expands outward. (2) A central core forms, with a mass of ~1 Mo, supported by thermal gas pressure gradients and rotation. (3) Growing thermal gas and radiation pressure gradients dissolve it. (4) This process is associated with a strong anisotropic outflow, and another core forms nearby and grows rapidly. (5) Typical radiation luminosity emerging from the photosphere encompassing these cores is ~5 x 10^{37}-5 x 10^{38} erg/s, of order the Eddington luminosity. (6) Two variability timescales are associated with this process: a long one, which is related to the accretion flow within the central ~10^{-4}-10^{-3} pc, and ~0.1 yr, which is related to radiation diffusion. (7) Adiabatic models have been run for comparison and their evolution differs profoundly from that of the FLD models, by forming a central geometrically-thick disk. Overall, an adiabatic equation of state is not a good approximation to the advanced stage of direct collapse, mainly because the radiation in the FLD is capable of escaping due to anisotropy in the optical depth and associated gradients.
Observations of high-redshift quasars imply the presence of supermassive black holes already at z~ 7.5. An appealing and promising pathway to their formation is the direct collapse scenario of a primordial gas in atomic-cooling haloes at z ~ 10 - 20, when the $rm H_2$ formation is inhibited by a strong background radiation field, whose intensity exceeds a critical value, $J_{rm crit}$. To estimate $J_{rm crit}$, typically, studies have assumed idealized spectra, with a fixed ratio of $rm H_{2}$ photo-dissociation rate $k_{rm H_2}$ to the $rm H^-$ photo-detachment rate $k_{rm H^-}$. This assumption, however, could be too narrow in scope as the nature of the background radiation field is not known precisely. In this work we argue that the critical condition for suppressing the $rm H_2$ cooling in the collapsing gas could be described in a more general way by a combination of $k_{rm H_2}$ and $k_{rm H^-}$ parameters. By performing a series of cosmological zoom-in simulations with an encompassing set of $k_{rm H_2}$ and $k_{rm H^-}$, we examine the gas flow by following evolution of basic parameters of the accretion flow. We test under what conditions the gas evolution is dominated by $rm H_{2}$ and/or atomic cooling. We confirm the existence of a critical curve in the $k_{rm H_2}-k_{rm H^-}$ plane, and provide an analytical fit to it. This curve depends on the conditions in the direct collapse, and reveals domains where the atomic cooling dominates over the molecular cooling. Furthermore, we have considered the effect of $rm H_{2}$ self-shielding on the critical curve, by adopting three methods for the effective column density approximation in $rm H_{2}$. We find that the estimate of the characteristic length-scale for shielding can be improved by using $lambda_{rm Jeans25}$, which is 0.25 times that of the local Jeans length.
We study the early stage of the formation of seed supermassive black holes via direct collapse in dark matter (DM) halos, in the cosmological context. We perform high-resolution zoom-in simulations of such collapse at high-$z$. Using the adaptive mesh refinement code ENZO, we resolve the formation and growth of a DM halo, until its virial temperature reaches $sim 10^4$K, atomic cooling turns on, and collapse ensues. We demonstrate that direct collapse proceeds in two stages, although they are not well separated. The first stage is triggered by the onset of atomic cooling, and leads to rapidly increasing accretion rate with radius, from $dot Msim 0.1,M_odot {rm yr^{-1}}$ at the halo virial radius to few $M_odot ,{rm yr^{-1}}$, around the scale radius $R_{rm s}sim 30$pc of the NFW DM density profile. The second stage of the collapse commences when the gas density takes precedence over the DM density. This is associated with the gas decoupling from the DM gravitational potential. The ensuing collapse approximates that of an isothermal sphere with $dot M ( r )sim $const. We confirm that the gas loses its angular momentum through non-axisymmetric perturbations and gravitational torques, to overcome the centrifugal barrier. During the course of the collapse, this angular momentum transfer process happens on nearly all spatial scales, and the angular momentum vector of the gas varies with position and time. Collapsing gas also exhibits supersonic turbulent motions which suppress gas fragmentation, and are characterized by density PDF consisting of a lognormal part and a high-density power law tail.
We study the collapse of rapidly rotating supermassive stars that may have formed in the early Universe. By self-consistently simulating the dynamics from the onset of collapse using three-dimensional general-relativistic hydrodynamics with fully dynamical spacetime evolution, we show that seed perturbations in the progenitor can lead to the formation of a system of two high-spin supermassive black holes, which inspiral and merge under the emission of powerful gravitational radiation that could be observed at redshifts z>10 with the DECIGO or Big Bang Observer gravitational-wave observatories, assuming supermassive stars in the mass range 10^4-10^6 Msol. The remnant is rapidly spinning with dimensionless spin a^*=0.9. The surrounding accretion disk contains ~10% of the initial mass.
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.