The leading contenders for the seeds of the first quasars are direct collapse black holes (DCBHs) formed during catastrophic baryon collapse in atomically-cooled halos at $z \sim$ 20. The discovery of the Ly$\alpha$ emitter CR7 at $z =$ 6.6 was initially held to be the first detection of a DCBH, although this interpretation has since been challenged on the grounds of Spitzer IRAC and Very Large Telescope X-Shooter data. Here we determine if radio flux from a DCBH in CR7 could be detected and discriminated from competing sources of radio emission in the halo such as young supernovae and H II regions. We find that a DCBH would emit a flux of 10 - 200 nJy at 1.0 GHz, far greater than the sub-nJy signal expected for young supernovae but on par with continuum emission from star-forming regions. However, radio emission from a DCBH in CR7 could be distinguished from free-free emission from H II regions by its spectral evolution with frequency and could be detected by the Square Kilometer Array in the coming decade.

Earliest quasars at the cosmic dawn are powered by mass accretion onto supermassive black holes of a billion solar masses. Massive black hole seeds forming through the direct collapse mechanism are considered the most promising candidates but how do they grow and co-evolve with their host galaxies at early cosmic times remains unknown. We here present results from a cosmological radiation hydrodynamical simulation including self-consistent modeling of both Pop III and Pop II star formation, their radiative and supernova feedback in the host galaxy along with X-ray feedback from an accreting massive black hole (MBH) of $\rm 10^5 ~M_{\odot}$ in a halo of $\rm 2 \times 10^9~M_{\odot}$ from $z=26$ down to $z=16$. Our results show that energy deposition from X-rays in the proximity of MBH suppresses Pop III star formation for about 12 Myr while at the same time these X-rays catalyze $\rm H_2$ formation which leads to the formation of a Pop III star cluster of 500 $\rm M_{\odot}$ in the close vicinity of the MBH. We find that mode of star formation for Pop III is episodic and bursty due to the clumpy accretion while for Pop II it is continuous. The stellar mass of the host galaxy at $z \sim 16$ is $\rm 2 \times 10^7~M_{\odot}$ with a star formation rate (SFR) of $\rm \sim 0.1-1~M_{\odot}/yr$. In total, the MBH accretes $\rm 1.5 \times 10^6~M_{\odot}$ during 120 Myr with the mean accretion rate of $\rm \sim 0.01~M_{\odot}/yr$ corresponding to an average Eddington fraction of 50\%.

Supermassive primordial stars forming during catastrophic baryon collapse in atomically-cooling halos at $z \sim$ 15 - 20 may be the origin of the first quasars in the universe. However, no simulation to date has followed the evolution of these halos at resolutions that are high enough or for times that are long enough to determine if collapse actually produces SMSs. Here we report new cosmological simulations of baryon collapse in atomically-cooled halos for times that are long enough for SMSs to form and die as direct-collapse black holes (DCBHs). We find that the high infall rates required to build up such stars do persist until the end of their lives and could fuel the rapid growth of their BHs thereafter. Our simulations also demonstrate that binary and even small multiples of SMSs can form in low-spin and high-spin halos, respectively. This discovery raises the exciting prospect of detecting gravitational waves from DCBH mergers with LISA and tidal disruption events in the near infrared with the {\em James Webb Space Telescope} and ground-based telescopes in the coming decade.

Supermassive black holes (BHs) of millions solar masses and above reside in the center of most local galaxies, but they also power active galactic nuclei and quasars, detected up to z=7. This observational evidence puts strong constraints on the BH growth and the mass of the first BH seeds. The scenario of "direct collapse" is very appealing as it leads to the formation of large mass BH seeds in the range 10^4-10^6 Msun, which eases explaining how quasars at z=6-7 are powered by BHs with masses >10^9 Msun. Direct collapse, however, appears to be rare, as the conditions required by the scenario are that gas is metal-free, the presence of a strong photo-dissociating Lyman-Werner flux, and large inflows of gas at the center of the halo, sustained for 10-100 Myr. We performed several cosmological hydrodynamical simulations that cover a large range of box sizes and resolutions, thus allowing us to understand the impact of several physical processes on the distribution of direct collapse BHs. We identify halos where direct collapse can happen, and derive the number density of BHs. We also investigate the discrepancies between hydrodynamical simulations, direct or post-processed, and semi-analytical studies. We find that for direct collapse to account for BHs in normal galaxies, the critical Lyman-Werner flux required for direct collapse must be much less than predicted by 3D simulations that include detailed chemical models. However, when supernova feedback is relatively weak, enough direct collapse BHs to explain z=6-7 quasars can be obtained for more realistic values of the critical Lyman Werner flux.

The first galaxies forming a few hundred million years after the Big Bang are the key drivers of cosmic evolution and ideal laboratories to study theories of galaxy formation. We here study the role of UV radiation in suppressing star formation in primordial galaxies by destroying molecular hydrogen, the main coolant in primordial gas and provide estimates of cold dense gas at the onset of star formation. To accomplish this goal, we perform three dimensional cosmological simulations of minihalos in different environments forming at $ z\sim 25$ by varying strength of background UV flux below the Lyman limit between 0.01-1000 in units of $\rm J_{21}=10^{-21}~erg/cm^2/s/Hz/sr$. Particularly, we include photo-detachment of $\rm H^-$, the self-shielding of $\rm H_2$ which both were neglected in previous studies and use updated reaction rates. Our results show that depending on the background level $\rm H_2$ formation is suppressed, delaying gravitational collapse until halos reach the atomic cooling limit. We find that the formation of cold dense molecular gas and subsequently star formation gets delayed by 100 to 230 Myr depending on the level of the background radiation and the growth history of the dark matter halos. The fraction of dense self-shielded gas is a strong function of the background flux and exponentially declines with the strength of incident UV flux above $\rm J_{21} \geq 1$. We find that taking into account $\rm H_2$ self-shielding is crucial for accurately estimating the amount of cold dense gas available for star formation.

Observations reveal the presence of supermassive black holes (SMBH) as early as ~700 million years after the Big Bang. Their formation path is still subject to current debate. We explore the influence of magnetic fields, which are strongly amplified via the turbulent small-scale dynamo, on the formation of SMBH seeds within the direct collapse scenario. In this study, we perform for the first time cosmological magnetohydrodynamic large eddy simulations that employ a model for unresolved, compressible MHD turbulence. In total we perform 36 simulations for 9 haloes each with two different initial magnetic field strengths, and with and without employing the unresolved turbulence model. We make use of the adaptive mesh refinement approach to achieve an effective spatial resolution of less than one proper astronomical unit. We consider a regime where cooling is regulated by atomic hydrogen and the molecular hydrogen gets dissociated by a strong radiation field. Our main finding is that the majority of the gas properties in the haloes at the final output are predominantly determined by the run-away gravitational collapse. Turbulence is supersonic and super-Alfv\'enic in all cases, and magnetic fields are amplified to an approximately dynamically relevant regime. Finally, fragmentation during the collapse is intermittent and mass accretion rates range from 0.2-3 Msun/yr. This suggests that the presence of strongly amplified magnetic fields and turbulence provides additional pressure support on small scales and make the direct collapse a viable scenario for the formation of massive objects under the required ambient conditions.

In recent years, the discovery of massive quasars at z~7 has provided a striking challenge to our understanding of the origin and growth of supermassive black holes in the early Universe. Mounting observational and theoretical evidence indicates the viability of massive seeds, formed by the collapse of supermassive stars, as a progenitor model for such early, massive accreting black holes. Although considerable progress has been made in our theoretical understanding, many questions remain regarding how (and how often) such objects may form, how they live and die, and how next generation observatories may yield new insight into the origin of these primordial titans. This review focusses on our present understanding of this remarkable formation scenario, based on discussions held at the Monash Prato Centre from November 20--24, 2017, during the workshop "Titans of the Early Universe: The Origin of the First Supermassive Black Holes."

Understanding the formation of earliest supermassive black holes is a question of prime astrophysical interest. In this chapter, we focus on the formation of massive black holes via gas dynamical processes. The necessary requirement for this mechanism are large inflow rates of about 0.1 solar mass per year. We discuss how to obtain such inflow rates via an isothermal collapse in the presence of atomic hydrogen cooling, and the outcome of such a collapse from three dimensional cosmological simulations in subsection 2.2. Alternatives to an isothermal direct collapse are discussed in subsection 3 which include trace amounts of metals and/or molecular hydrogen. In the end, we briefly discuss future perspectives and potential detection of massive black hole seeds via upcoming missions.

Black holes with masses of $\rm 10^6-10^9~M_{\odot}$ dwell in the centers of most galaxies, but their formation mechanisms are not well known. A subdominant dissipative component of dark matter with similar properties to the ordinary baryons, known as mirror dark matter, may collapse to form massive black holes during the epoch of first galaxies formation. In this study, we explore the possibility of massive black hole formation via this alternative scenario. We perform three-dimensional cosmological simulations for four distinct halos and compare their thermal, chemical and dynamical evolution in both the ordinary and the mirror sectors. We find that the collapse of halos is significantly delayed in the mirror sector due to the lack of $\rm H_2$ cooling and only halos with masses above $ \rm \geq 10^7~ M_{\odot}$ are formed. Overall, the mass inflow rates are $\rm \geq 10^{-2}~M_{\odot}/yr$ and there is less fragmentation. This suggests that the conditions for the formation of massive objects, including black holes, are more favorable in the mirror sector.

Understanding the growth of high redshift massive black holes (MBHs) is a problem of great astrophysical interest. The most luminous quasars at $z>6$ are frequently observed but they represent only the tip of the iceberg as the majority of the low luminosity AGN population remains undetected. In the present study, we perform a radiation hydrodynamics cosmological simulation to study the growth of "normal" black holes in the high redshift universe. In our simulation we model the formation of Pop III and Pop II stars along with their chemical, mechanical and radiative feedback. We consider both UV and X-ray emission from an accreting BH to simulate its radiative feedback. The selected halo has a mass of $\rm 3 \times 10^{10}~M_{\odot}$ at $z=7.5$ and we turn on radiative feedback from a MBH seed of $\rm 10^5~M_{\odot}$ along with in-situ star formation at $z=12$ when the halo mass reaches well above the atomic cooling limit. We find that the MBH accretes only about 2200 $\rm M_{\odot}$ during 320 Myr and the average mass accretion onto the MBH is a few times $\rm 10^{-6}~M_{\odot}/yr$. Our results suggest that the stunted growth of MBH is a consequence of supernovae in tandem with MBH feedback which drive large outflows and evacuate the gas from MBH vicinity. This may explain why a population of low luminosity AGN has not been detected so-far at $z>6$; the large contrast between the star formation rate and the MBH accretion rate may make then hard to detect even in upcoming deep surveys.

A massive primordial halo near an intensely star-forming galaxy may collapse into a supermassive star (SMS) and leave a massive black hole seed of about $10^5~M_{\odot}$. To investigate the impact of ionizing radiation on the formation of an SMS from a nearby galaxy, we perform three-dimensional radiation hydrodynamical simulations by selecting a pair of massive dark matter halos forming at $z >10$. We find that rich structures such as clumps and filaments around the source galaxy shield the cloud from ionizing radiation. In fact, in some cases, cloud collapse is accelerated under ionizing radiation. This fact suggests that the ionization of the cloud's surroundings helps its collapse. Only strong radiation at the early stage of structure formation can halt the cloud collapse, but this is much stronger than observationally allowed value. We also explored the effect of ionizing radiation on a sample of 68 halos by employing an analytical model and found that increase in the mean density of the gas between the SMS forming cloud and the source galaxy protects the gas cloud from ionizing radiation as they approach each other. Thus, we conclude that ionizing radiation does not prevent the formation of an SMS in most of the cases.

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 detection of quasars at $z>6$ unveils the presence of supermassive black holes (BHs) of a few billion solar masses. The rapid formation process of these extreme objects remains a fascinating and open issue. Such discovery implies that seed black holes must have formed early on, and grown via either rapid accretion or BH/galaxy mergers. In this theoretical review, we discuss in detail various BH seed formation mechanisms and the physical processes at play during their assembly. We discuss the three most popular BH formation scenarios, involving the (i) core-collapse of massive stars, (ii) dynamical evolution of dense nuclear star clusters, (iii) collapse of a protogalactic metal free gas cloud. This article aims at giving a broad introduction and an overview of the most advanced research in the field.

CR7 is the brightest Lyman-$\alpha$ emitter observed at $z>6$, which shows very strong Lyman-$\alpha$ and HeII 1640\AA\ line luminosities, but no metal line emission. Previous studies suggest that CR7 hosts either young primordial stars with a total stellar mass of $\sim 10^7\,\mathrm{M}_\odot$ or a black hole of $\gtrsim 10^6\,\mathrm{M}_\odot$. Here, we explore different formation scenarios for CR7 with a semianalytical model, based on the random sampling of dark matter merger trees. We are unable to reproduce the observational constraints with a primordial stellar source, given our model assumptions, due to the short stellar lifetimes and the early metal enrichment. Black holes that are the remnants of the first stars are either not massive enough, or reside in metal-polluted haloes, ruling out this possible explanation of CR7. Our models instead suggest that direct collapse black holes, which form in metal-free haloes exposed to large Lyman-Werner fluxes, are more likely the origin of CR7. However, this result is derived under optimistic assumptions and future observations are necessary to further constrain the nature of CR7.

Observations of quasars at $ z > 6$ suggest the presence of black holes with a few times $\rm 10^9 ~M_{\odot}$. Numerous models have been proposed to explain their existence including the direct collapse which provides massive seeds of $\rm 10^5~M_{\odot}$. The isothermal direct collapse requires a strong Lyman-Werner flux to quench $\rm H_2$ formation in massive primordial halos. In this study, we explore the impact of trace amounts of metals and dust enrichment. We perform three dimensional cosmological simulations for two halos of $\rm > 10^7~M_{\odot}$ with $\rm Z/Z_{\odot}= 10^{-4}-10^{-6}$ illuminated by an intense Lyman Werner flux of $\rm J_{21}=10^5$. Our results show that initially the collapse proceeds isothermally with $\rm T \sim 8000$ K but dust cooling becomes effective at densities of $\rm 10^{8}-10^{12} ~cm^{-3}$ and brings the gas temperature down to a few 100-1000 K for $\rm Z/Z_{\odot} \geq 10^{-6}$. No gravitationally bound clumps are found in $\rm Z/Z_{\odot} \leq 10^{-5}$ cases by the end of our simulations in contrast to the case with $\rm Z/Z_{\odot} = 10^{-4}$. Large inflow rates of $\rm \geq 0.1~M_{\odot}/yr$ are observed for $\rm Z/Z_{\odot} \leq 10^{-5}$ similar to a zero-metallicity case while for $\rm Z/Z_{\odot} = 10^{-4}$ the inflow rate starts to decline earlier due to the dust cooling and fragmentation. For given large inflow rates a central star of $\rm \sim 10^4~M_{\odot}$ may form for $\rm Z/Z_{\odot} \leq 10^{-5}$.

The detection of $\rm z>6$ quasars reveals the existence of supermassive black holes of a few $\rm 10^9~M_{\odot}$. One of the potential pathways to explain their formation in the infant universe is the so-called direct collapse model which provides massive seeds of $\rm 10^5-10^6~M_{\odot}$. An isothermal direct collapse mandates that halos should be of a primordial composition and the formation of molecular hydrogen remains suppressed in the presence of a strong Lyman Werner flux. In this study, we perform high resolution cosmological simulations for two massive primordial halos employing a detailed chemical model which includes $\rm H^-$ cooling as well as realistic opacities for both the bound-free $\rm H^-$ emission and the Rayleigh scattering of hydrogen atoms. We are able to resolve the collapse up to unprecedentedly high densities of $\rm \sim 10^{-3}~g/cm^3$ and to scales of about $\rm 10^{-4}$ AU. Our results show that the gas cools down to $\rm \sim $ 5000 K in the presence of $\rm H^-$ cooling, and induces fragmentation at scales of about 8000 AU in one of the two simulated halos, which may lead to the formation of a binary. In addition, fragmentation also occurs on the AU scale in one of the halos but the clumps are expected to merge on short time scales. Our results confirm that $\rm H^-$ cooling does not prevent the formation of a supermassive star and the trapping of cooling radiation stabilises the collapse on small scales.

Quasars powered by massive black holes (BHs) with mass estimates above a billion solar masses have been identified at redshift 6 and beyond. The existence of such BHs requires almost continuous growth at the Eddington limit for their whole lifetime, of order of one billion years. In this paper, we explore the possibility that positively skewed scale-dependent non-Gaussian primordial fluctuations may ease the assembly of massive BHs. In particular, they produce more low-mass halos at high redshift, thus altering the production of metals and ultra-violet flux, believed to be important factors in BH formation. Additionally, a higher number of progenitors and of nearly equal-mass halo mergers would boost the mass increase provided by BH-BH mergers and merger-driven accretion. We use a set of two cosmological simulations, with either Gaussian or scale-dependent non-Gaussian primordial fluctuations to perform a proof-of-concept experiment to estimate how BH formation and growth are altered. We estimate the BH number density and the fraction of halos where BHs form, for both simulations and for two popular scenarios of BH formation (remnants of the first generation of stars and direct collapse in the absence of metals and molecular hydrogen). We find that the fractions of halos where BHs form are almost identical, but that non-Gaussian primordial perturbations increase the total number density of BHs for the both BH formation scenarios. We also evolve BHs using merger trees extracted from the simulations and find that non-Gaussianities increase both the BH masses and the number of the most massive BHs.

Magnetic fields are considered as a vital ingredient of contemporary star formation, and may have been important during the formation of the first stars in the presence of an efficient amplification mechanism. Initial seed fields are provided via plasma fluctuations, and are subsequently amplified by the small-scale dynamo, leading to a strong tangled magnetic field. Here we explore how the magnetic field provided by the small-scale dynamo is further amplified via the $\alpha-\Omega$ dynamo in a protostellar disk and assess its implications. For this purpose, we consider two characteristic cases, a typical Pop.~III star with $10$~M$_\odot$ and an accretion rate of $10^{-3}$~M$_\odot$~yr$^{-1}$, and a supermassive star with $10^5$~M$_\odot$ and an accretion rate of $10^{-1}$~M$_\odot$~yr$^{-1}$. For the $10$~M$_\odot$ Pop.~III star, we find that coherent magnetic fields can be produced on scales of at least $100$~AU, which are sufficient to drive a jet with a luminosity of $100$~L$_\odot$ and a mass outflow rate of $10^{-3.7}$~M$_\odot$~yr$^{-1}$. For the supermassive star, the dynamical timescales in its environment are even shorter, implying smaller orbital timescales and an efficient magnetization out to at least $1000$~AU. The jet luminosity corresponds to $\sim10^{6.0}$~L$_\odot$, and a mass outflow rate of $10^{-2.1}$~M$_\odot$~yr$^{-1}$. We expect that the feedback from the supermassive star can have a relevant impact on its host galaxy.

Numerical simulations show the formation of self-gravitating primordial disks during the assembly of the first structures in the Universe, in particular during the formation of Pop.~III and supermassive stars. Their subsequent evolution is expected to be crucial to determine the mass scale of the first cosmological objects, which depends on the temperature of the gas and the dominant cooling mechanism. Here, we derive a one-zone framework to explore the chemical evolution of such disks and show that viscous heating leads to the collisional dissociation of an initially molecular gas. The effect is relevant on scales of 10 AU (1000 AU) for a central mass of 10 M_sun (10^4 M_sun) at an accretion rate of 10^{-1} M_sun yr^{-1}, and provides a substantial heat input to stabilize the disk. If the gas is initially atomic, it remains atomic during the further evolution, and the effect of viscous heating is less significant. The additional thermal support is particularly relevant for the formation of very massive objects, such as the progenitors of the first supermassive black holes. The stabilizing impact of viscous heating thus alleviates the need for a strong radiation background as a means of keeping the gas atomic.

High redshift quasars at z>6 have masses up to ~$10^9$ M$_\odot$. One of the pathways to their formation includes direct collapse of gas, forming a supermassive star, precursor of the black hole seed. The conditions for direct collapse are more easily achievable in metal-free haloes, where atomic hydrogen cooling operates and molecular hydrogen (H2) formation is inhibited by a strong external UV flux. Above a certain value of UV flux (J_crit), the gas in a halo collapses isothermally at ~$10^4$ K and provides the conditions for supermassive star formation. However, H2 can self-shield, reducing the effect of photodissociation. So far, most numerical studies used the local Jeans length to calculate the column densities for self-shielding. We implement an improved method for the determination of column densities in 3D simulations and analyse its effect on the value of J_crit. This new method captures the gas geometry and velocity field and enables us to properly determine the direction-dependent self-shielding factor of H2 against photodissociating radiation. We find a value of J_crit that is a factor of two smaller than with the Jeans approach (~2000 J_21 vs. ~4000 J_21). The main reason for this difference is the strong directional dependence of the H2 column density. With this lower value of J_crit, the number of haloes exposed to a flux >J_crit is larger by more than an order of magnitude compared to previous studies. This may translate into a similar enhancement in the predicted number density of black hole seeds.

Supermassive black holes are not only common in the present-day galaxies, but billion solar masses black holes also powered $z\geq 6$ quasars. One efficient way to form such black holes is the collapse of a massive primordial gas cloud into a so-called direct collapse black hole. The main requirement for this scenario is the presence of large accretion rates of $\rm \geq 0.1~M_{\odot}/yr$ to form a supermassive star. It is not yet clear how and under what conditions such accretion rates can be obtained. The prime aim of this work is to determine the mass accretion rates under non-isothermal collapse conditions. We perform high resolution cosmological simulations for three primordial halos of a few times $\rm 10^7~M_{\odot}$ illuminated by an external UV flux, $\rm J_{21}=100-1000$. We find that a rotationally supported structure of about parsec size is assembled, with an aspect ratio between $\rm 0.25 - 1$ depending upon the thermodynamical properties. Rotational support, however, does not halt collapse, and mass inflow rates of $\rm \sim 0.1~M_{\odot}/yr$ can be obtained in the presence of even a moderate UV background flux of strength $\rm J_{21} \geq 100$. To assess whether such large accretion rates can be maintained over longer time scales, we employed sink particles, confirming the persistence of accretion rates of $\rm \sim 0.1~M_{\odot}/yr$. We propose that complete isothermal collapse and molecular hydrogen suppression may not always be necessary to form supermassive stars, precursors of black hole seeds. Sufficiently high inflow rates can be obtained for UV flux $\rm J_{21}=500-1000$, at least for some cases. This value brings the estimate of the abundance of direct collapse black hole seeds closer to that high redshift quasars.

Massive primordial halos exposed to moderate UV backgrounds are the potential birthplaces of supermassive black holes. In such a halo, an initially isothermal collapse will occur, leading to high accretion rates of $\sim0.1$~M$_\odot$~yr$^{-1}$. During the collapse, the gas in the interior will turn into a molecular state, and form an accretion disk due to the conservation of angular momentum. We consider here the structure of such an accretion disk and the role of viscous heating in the presence of high accretion rates for a central star of $10$, $100$ and $10^4$~M$_\odot$. Our results show that the temperature in the disk increases considerably due to viscous heating, leading to a transition from the molecular to the atomic cooling phase. We found that the atomic cooling regime may extend out to several $100$~AU for a $10^4$~M$_\odot$ central star and provides substantial support to stabilize the disk. It therefore favors the formation of a massive central object. The comparison of clump migration and contraction time scales shows that stellar feedback from these clumps may occur during the later stages of the evolution. Overall, viscous heating provides an important pathway to obtain an atomic gas phase within the center of the halo, and helps in the formation of very massive objects. The latter may collapse to form a massive black hole of about $\geq 10^4$~M$_\odot$.

Our understanding of population III star formation is still in its infancy. They are formed in dark matter minihalos of $\rm 10^5-10^6 M_{\odot}$ at $z=20-30$. Recent high resolution cosmological simulations show that a protostellar disk forms as a consequence of gravitational collapse and fragments into multiple clumps. However, it is not entirely clear if these clumps will be able to survive to form multiple stars as simulations are unable to follow the disk evolution for longer times. In this study, we employ a simple analytical model to derive the properties of marginally stable steady-state disks. Our results show that the stability of the disk depends on the critical value of the viscous parameter $\alpha$. For $\alpha_{crit} = 1$, the disk is stable for an accretion rate of $\rm \leq 10^{-3} M_{\odot}/yr$ and becomes unstable at radii about $\rm \geq 100 AU$ in the presence of an accretion rate of $\rm 10^{-2} M_{\odot}/yr$. For $0.06 < \alpha_{crit} < 1$, the disk can be unstable for both accretion rates. The comparison of the migration and the Kelvin-Helmholtz time scales shows that clumps are expected to migrate inward before reaching the main sequence. Furthermore, in the presence of a massive central star the clumps within the central 1 AU will be tidally disrupted. We also find that UV feedback from the central star is unable to disrupt the disk, and that photo-evaporation becomes important only once the accretion rate has dropped to $\rm 2 \times 10^{-4} M_{\odot}/yr$. As a result, the central star may reach a mass of 100 $\rm M_{\odot}$ or even higher.

Observations of high redshift quasars at $z>6$ indicate that they harbor supermassive black holes (SMBHs) of a billion solar masses. The direct collapse scenario has emerged as the most plausible way to assemble SMBHs. The nurseries for the direct collapse black holes are massive primordial halos illuminated with an intense UV flux emitted by population II (Pop II) stars. In this study, we compute the critical value of such a flux ($J_{21}^{\rm crit}$) for realistic spectra of Pop II stars through three-dimensional cosmological simulations. We derive the dependence of $J_{21}^{\rm crit}$ on the radiation spectra, on variations from halo to halo, and on the impact of X-ray ionization. Our findings show that the value of $J_{21}^{\rm crit}$ is a few times $\rm 10^4$ and only weakly depends on the adopted radiation spectra in the range between $T_{\rm rad}=2 \times 10^4-10^5$ K. For three simulated halos of a few times $\rm 10^{7}$~M$_{\odot}$, $J_{21}^{\rm crit}$ varies from $\rm 2 \times 10^4 - 5 \times 10^4$. The impact of X-ray ionization is almost negligible and within the expected scatter of $J_{21}^{\rm crit}$ for background fluxes of $J_{\rm X,21} \leq 0.1$. The computed estimates of $J_{21}^{\rm crit}$ have profound implications for the quasar abundance at $z=10$ as it lowers the number density of black holes forming through an isothermal direct collapse by a few orders of magnitude below the observed black holes density. However, the sites with moderate amounts of $\rm H_2$ cooling may still form massive objects sufficient to be compatible with observations.

Context. The seeds of the first supermassive black holes may have resulted from the direct collapse of hot primordial gas in $\gtrsim 10^4$ K haloes, forming a supermassive or quasistar as an intermediate stage. Aims. We explore the formation of a protostar resulting from the collapse of primordial gas in the presence of a strong Lyman-Werner radiation background. Particularly, we investigate the impact of turbulence and rotation on the fragmentation behaviour of the gas cloud. We accomplish this goal by varying the initial turbulent and rotational velocities. Methods. We performed 3D adaptive mesh refinement simulations with a resolution of 64 cells per Jeans length using the ENZO code, simulating the formation of a protostar up to unprecedentedly high central densities of $10^{21}$ cm$^{-3}$, and spatial scales of a few solar radii. To achieve this goal, we employed the KROME package to improve modelling of the chemical and thermal processes. Results. We find that the physical properties of the simulated gas clouds become similar on small scales, irrespective of the initial amount of turbulence and rotation. After the highest level of refinement was reached, the simulations have been evolved for an additional ~5 freefall times. A single bound clump with a radius of $2 \times 10^{-2}$ AU and a mass of ~$7 \times 10^{-2}$ M$_{\odot}$ is formed at the end of each simulation, marking the onset of protostar formation. No strong fragmentation is observed by the end of the simulations, regardless of the initial amount of turbulence or rotation, and high accretion rates of a few solar masses per year are found. Conclusions. Given such high accretion rates, a quasistar of $10^5$ M$_{\odot}$ is expected to form within $10^5$ years.

The ability of metal free gas to cool by molecular hydrogen in primordial halos is strongly associated with the strength of ultraviolet (UV) flux produced by the stellar populations in the first galaxies. Depending on the stellar spectrum, these UV photons can either dissociate $\rm H_{2}$ molecules directly or indirectly by photo-detachment of $\rm H^{-}$ as the latter provides the main pathway for $\rm H_{2}$ formation in the early universe. In this study, we aim to determine the critical strength of the UV flux above which the formation of molecular hydrogen remains suppressed for a sample of five distinct halos at $z>10$ by employing a higher order chemical solver and a Jeans resolution of 32 cells. We presume that such flux is emitted by PopII stars implying atmospheric temperatures of $\rm 10^{4}$~K. We performed three-dimensional cosmological simulations and varied the strength of the UV flux below the Lyman limit in units of $\rm J_{21}$. Our findings show that the value of $\rm J_{21}^{crit}$ varies from halo to halo and is sensitive to the local thermal conditions of the gas. For the simulated halos it varies from 400-700 with the exception of one halo where $\rm J_{21}^{crit} \geq 1500$. This has important implications for the formation of direct collapse black holes and their estimated population at z > 6. It reduces the number density of direct collapse black holes by almost three orders of magnitude compared to the previous estimates.

Recent discoveries of carbon-enhanced metal-poor stars like SMSS J031300.36-670839.3 provide increasing observational insights into the formation conditions of the first second-generation stars in the Universe, reflecting the chemical conditions after the first supernova explosion. Here, we present the first cosmological simulations with a detailed chemical network including primordial species as well as C, C$^+$, O, O$^+$, Si, Si$^+$, and Si$^{2+}$ following the formation of carbon-enhanced metal poor stars. The presence of background UV flux delays the collapse from $z=21$ to $z=15$ and cool the gas down to the CMB temperature for a metallicity of Z/Z$_\odot$=10$^{-3}$. This can potentially lead to the formation of lower mass stars. Overall, we find that the metals have a stronger effect on the collapse than the radiation, yielding a comparable thermal structure for large variations in the radiative background. We further find that radiative backgrounds are not able to delay the collapse for Z/Z$_\odot$=10$^{-2}$ or a carbon abundance as in SMSS J031300.36-670839.3.

Radiative feedback from populations II stars played a vital role in early structure formation. Particularly, photons below the Lyman limit can escape the star forming regions and produce a background ultraviolet (UV) flux which consequently may influence the pristine halos far away from the radiation sources. These photons can quench the formation of molecular hydrogen by photo-detachment of $\rm H^{-}$. In this study, we explore the impact of such UV radiation on fragmentation in massive primordial halos of a few times $\rm 10^{7}$~M${_\odot}$. To accomplish this goal, we perform high resolution cosmological simulations for two distinct halos and vary the strength of the impinging background UV field in units of $\rm J_{21}$. We further make use of sink particles to follow the evolution for 10,000 years after reaching the maximum refinement level. No vigorous fragmentation is observed in UV illuminated halos while the accretion rate changes according to the thermal properties. Our findings show that a few 100-10, 000 solar mass protostars are formed when halos are irradiated by $\rm J_{21}=10-500$ at $\rm z>10$ and suggest a strong relation between the strength of UV flux and mass of a protostar. This mode of star formation is quite different from minihalos, as higher accretion rates of about $\rm 0.01-0.1$ M$_{\odot}$/yr are observed by the end of our simulations. The resulting massive stars are the potential cradles for the formation of intermediate mass black holes at earlier cosmic times and contribute to the formation of a global X-ray background.

While Population III stars are typically thought to be massive, pathways towards lower-mass Pop III stars may exist when the cooling of the gas is particularly enhanced. A possible route is enhanced HD cooling during the merging of dark-matter halos. The mergers can lead to a high ionization degree catalysing the formation of HD molecules and may cool the gas down to the cosmic microwave background (CMB) temperature. In this paper, we investigate the merging of mini-halos with masses of a few 10$^5$ M$_\odot$ and explore the feasibility of this scenario. We have performed three-dimensional cosmological hydrodynamics calculations with the ENZO code, solving the thermal and chemical evolution of the gas by employing the astrochemistry package KROME. Our results show that the HD abundance is increased by two orders of magnitude compared to the no-merging case and the halo cools down to $\sim$60 K triggering fragmentation. Based on Jeans estimates the expected stellar masses are about 10 M$_\odot$. Our findings show that the merging scenario is a potential pathway for the formation of low-mass stars.

Baryonic streaming motions produced prior to the epoch of recombination became supersonic during the cosmic dark ages. Various studies suggest that such streaming velocities change the halo statistics and also influence the formation of Population III stars. In this study, we aim to explore the impact of streaming velocities on the formation of supermassive black holes at $z>10$ via the direct collapse scenario. To accomplish this goal, we perform cosmological large eddy simulations for two halos of a few times $\rm 10^{7} M_{\odot}$ with initial streaming velocities of 3, 6 and 9 $\rm km/s$. These massive primordial halos illuminated by the strong Lyman Werner flux are the potential cradles for the formation of direct collapse seed black holes. To study the evolution for longer times, we employ sink particles and track the accretion for 10,000 years. Our findings show that higher streaming velocities increase the circular velocities from about 14 $\rm km/s$ to 16 $\rm km/s$. They also delay the collapse of halos for a few million years, but do not have any significant impact on the halo properties such as turbulent energy, radial velocity, density and accretion rates. Sink particles of about $\rm \sim 10^5 M_{\odot}$ are formed at the end of our simulations and no clear distribution of sink masses is observed in the presence of streaming motions. It is further found that the impact of streaming velocities is less severe in massive halos compared to the minihalos as reported in the previous studies.

Observations of quasars at $\rm z> 6$ report the existence of a billion solar mass black holes. Comprehending their formation in such a short time scale is a matter of ongoing research. One of the most promising scenarios to assemble supermassive black holes is a monolithic collapse of protogalactic gas clouds in atomic cooling halos with $\rm T_{vir} \geq 10^{4} K$. In this article, we study the amplification and impact of magnetic fields during the formation of seed black holes in massive primordial halos. We perform high resolution cosmological magnetohydrodynamics simulations for four distinct halos and follow their collapse for a few free-fall times until the simulations reach a peak density of $\rm 7 \times 10^{-10} g/cm^{3}$. Our findings show that irrespective of the initial seed field, the magnetic field strength reaches a saturated state in the presence of strong accretion shocks. Under such conditions, the growth time becomes very short and amplification occurs rapidly within a small fraction of the free-fall time. We find that the presence of such strong magnetic fields provides additional support against gravity and helps in suppressing fragmentation. Massive clumps of a few hundred solar masses are formed at the end of our simulations and high accretion rates of $\rm 1 M_{\odot}/yr$ are observed. We expect that in the presence of such accretion rates, the clumps will grow to form supermassive stars of $\rm \sim 10^{5} M_{\odot}$. Overall, the role of the magnetic fields seems supportive for the formation of massive black holes.

Black holes of a billion solar masses are observed in the infant universe a few hundred million years after the Big Bang. The direct collapse of protogalactic gas clouds in primordial halos with $\rm T_{vir} \geq 10^{4} K$ provides the most promising way to assemble massive black holes. In this study, we aim to determine the characteristic mass scale of seed black holes and the time evolution of the accretion rates resulting from the direct collapse model. We explore the formation of supermassive black holes via cosmological large eddy simulations (LES) by employing sink particles and following their evolution for twenty thousand years after the formation of the first sink. As the resulting protostars were shown to have cool atmospheres in the presence of strong accretion, we assume here that UV feedback is negligible during this calculation. We confirm this result in a comparison run without sinks. Our findings show that black hole seeds with characteristic mass of $\rm 10^{5} M_{\odot}$ are formed in the presence of strong Lyman Werner flux which leads to an isothermal collapse. The characteristic mass is a about two times higher in LES compared to the implicit large eddy simulations (ILES). The accretion rates increase with time and reach a maximum value of 10 $\rm M_{\odot}/yr$ after $\rm 10^{4}$ years. Our results show that the direct collapse model is clearly feasible as it provides the expected mass of the seed black holes.

Population III stars forming in the infant universe at z=30 heralded the end of the cosmic dark ages. They are presumed to be assembled in so-called minihaloes with virial temperatures of a few thousand K where collapse is triggered by molecular hydrogen cooling. A central question concerns their final masses, and whether fragmentation occurs during their formation. While studies employing Lagrangian codes suggest fragmentation via a self-gravitating disk, recent high resolution simulations indicated that disk formation is suppressed. Here we report the first high-resolution large-eddy simulations performed with the Eulerian grid-based code Enzo following the evolution beyond the formation of the first peak, to investigate the accretion of the central massive clump and potential fragmentation. For a total of 3 halos, we see that a disk forms around the first clump. The central clump reaches $\sim10$ solar masses after 40 years, while subsequent accretion is expected at a rate of $10^{-2}$ solar masses per year. In one of these halos, additional clumps form as a result of fragmentation which proceeds at larger scales. We note that subgrid-scale turbulence yields relevant contributions to the stability of the protostellar disks. We conclude that the first protostar may reach masses up to $\rm 40-100 M_{\odot}$, which are only limited by the effect of radiative feedback.

Supermassive stars and quasi-stars (massive stars with a central black hole) are both considered as potential progenitors for the formation of supermassive black holes. They are expected to form from rapidly accreting protostars in massive primordial halos. We explore how long rapidly accreting protostars remain on the Hayashi track, implying large protostellar radii and weak accretion luminosity feedback. We assess the potential role of energy production in the nuclear core, and determine what regulates the evolution of such protostars into quasi-stars or supermassive stars. We follow the contraction of characteristic mass scales in rapidly accreting protostars, and infer the timescales for them to reach nuclear densities. We compare the characteristic timescales for nuclear burning with those for which the extended protostellar envelope can be maintained. We find that the extended envelope can be maintained up to protostellar masses of 3.6x10^8 \dot{m}^3 solar, where \dot{m} denotes the accretion rate in solar masses per year. We expect the nuclear core to exhaust its hydrogen content in 7x10^6 yrs. If accretion rates \dot{m}>>0.14 can still be maintained at this point, a black hole may form within the accreting envelope, leading to a quasi-star. Alternatively, the accreting object will gravitationally contract to become a main-sequence supermassive star. Due to the limited gas reservoir in dark matter halos with 10^7 solar masses, the accretion rate onto the central object may drop at late times, implying the formation of supermassive stars as the typical outcome of direct collapse. However, if high accretion rates are maintained, a quasi-star with an interior black hole may form.

The formation of the first stars in the Universe is regulated by a sensitive interplay of chemistry and cooling with the dynamics of a self-gravitating system. As the outcome of the collapse and the final stellar masses depend sensitively on the thermal evolution, it is necessary to accurately model the thermal evolution in high resolution simulations. As previous investigations raised doubts regarding the convergence of the temperature at high resolution, we investigate the role of the numerical method employed to model the chemistry and the thermodynamics. Here we compare the standard implementation in the adaptive-mesh refinement code \verb|ENZO|, employing a first order backward differentiation formula (BDF), with the 5th order accurate BDF solver \verb|DLSODES|. While the standard implementation in \verb|ENZO| shows a strong dependence on the employed resolution, the results obtained with \verb|DLSODES| are considerably more robust, both with respect to the chemistry and thermodynamics, but also for dynamical quantities such as density, total energy or the accretion rate. We conclude that an accurate modeling of the chemistry and thermodynamics is central for primordial star formation.

Supermassive black holes with up to a $\rm 10^{9}~M_{\odot}$ dwell in the centers of present-day galaxies, and their presence has been confirmed at z $\geq$ 6. Their formation at such early epochs is still an enigma. Different pathways have been suggested to assemble supermassive black holes in the first billion years after the Big Bang. Direct collapse has emerged as a highly plausible scenario to form black holes as it provides seed masses of $\rm 10^{5}-10^{6}~M_{\odot}$. Gravitational collapse in atomic cooling haloes with virial temperatures T$_{vir} \geq 10^{4}$~K may lead to the formation of massive seed black holes in the presence of an intense background UV flux. Turbulence plays a central role in regulating accretion and transporting angular momentum. We present here the highest resolution cosmological large-eddy simulations to date which track the evolution of high-density regions on scales of $0.25$~AU beyond the formation of the first peak, and study the impact of subgrid-scale turbulence. The peak density reached in these simulations is $\rm 1.2 \times 10^{-8}~g~cm^{-3}$. Our findings show that while fragmentation occasionally occurs, it does not prevent the growth of a central massive object resulting from turbulent accretion and occasional mergers. The central object reaches $\rm \sim 1000~M_{\odot}$ within $4$ free-fall times, and we expect further growth up to $\rm 10^{6}~M_{\odot}$ through accretion in about 1 million years. The direct collapse model thus provides a viable pathway of forming high-mass black holes at early cosmic times.

While present standard model of cosmology yields no clear prediction for the initial magnetic field strength, efficient dynamo action may compensate for initially weak seed fields via rapid amplification. In particular, the small-scale dynamo is expected to exponentially amplify any weak magnetic field in the presence of turbulence. We explore whether this scenario is viable using cosmological magneto-hydrodynamics simulations modeling the formation of the first galaxies, which are expected to form in so-called atomic cooling halos with virial temperatures $\rm T_{vir} \geq 10^{4}$ K. As previous calculations have shown that a high Jeans resolution is needed to resolve turbulent structures and dynamo effects, our calculations employ resolutions of up to 128 cells per Jeans length. The presence of the dynamo can be clearly confirmed for resolutions of at least 64 cells per Jeans length, while saturation occurs at approximate equipartition with turbulent energy. As a result of the large Reynolds numbers in primordial galaxies, we expect saturation to occur at early stages, implying magnetic field strengths of \sim0.1 $\mu$G at densities of 10^4 cm^{-3}.

Atomic cooling haloes with virial temperatures $\rm T_{vir} \geq 10^{4}$ K are the most plausible sites for the formation of the first galaxies and the first intermediate mass black holes. It is therefore important to assess whether one can obtain robust results concerning their main properties from numerical simulations. A major uncertainty is the presence of turbulence, which is barely resolved in cosmological simulations. We explore the latter both by pursuing high-resolution simulations with up to 64 cells per Jeans length and by incorporating a subgrid-scale turbulence model to account for turbulent pressure and viscosity on unresolved scales. We find that the main physical quantities in the halo, in particular the density, temperature and energy density profile, are approximately converged. However, the morphologies in the central 500 AU change significantly with increasing resolution and appear considerably more turbulent. In a systematic comparison of three different haloes, we further found that the turbulence subgrid-scale model gives rise to more compact central structures, and decreases the amount of vorticity. Such compact morphologies may in particular favor the accretion onto the central object.

We explore the amplification of magnetic fields in the high-redshift Universe. For this purpose, we perform high-resolution cosmological simulations following the formation of primordial halos with \sim10^7 M_solar, revealing the presence of turbulent structures and complex morphologies at resolutions of at least 32 cells per Jeans length. Employing a turbulence subgrid-scale model, we quantify the amount of unresolved turbulence and show that the resulting turbulent viscosity has a significant impact on the gas morphology, suppressing the formation of low-mass clumps. We further demonstrate that such turbulence implies the efficient amplification of magnetic fields via the small-scale dynamo. We discuss the properties of the dynamo in the kinematic and non-linear regime, and explore the resulting magnetic field amplification during primordial star formation. We show that field strengths of \sim10^{-5} G can be expected at number densities of \sim5 cm^{-3}.

Numerical simulations suggest that the first galaxies are formed in protogalactic halos with virial temperatures >= 10^4 K. It is likely that such halos are polluted with trace amounts of metals produced by the first generation of stars. The presence of dust can significantly change the chemistry and dynamics of early galaxies. In this article, we aim to assess the role of dust on the thermal and dynamical evolution of the first galaxies in the presence of a background UV flux, and its implications for the observability of Lyman alpha emitters and sub-mm sources. We have performed high resolution cosmological simulations using the adaptive mesh refinement code FLASH to accomplish this goal. We have developed a chemical network appropriate for these conditions and coupled it with the FLASH code. The main ingredients of our chemical model include the formation of molecules, a multi-level treatment of atomic hydrogen, line trapping of Lyman alpha photons and, photoionization and photodissociation processes in a UV background. We found that the formation of molecules (H_{2} and HD) is significantly enhanced in the presence of dust grains. The presence of a background UV flux strongly influences the formation of molecules by photodissociating them. We explore the evolution after a major merger, leading to the formation of a binary disk. These disks have gas masses of ~10^{7} M_sun at a redshift of 5.4. Each disk lies in a separate subhalo as a result of the merger event. The disks are supported by turbulent pressure due to the highly supersonic turbulence present in the halo. The presence of dust does not significantly reduce the Lyman alpha emission. The emission of Lyman alpha is extended and originates from the envelope of the halo due to line trapping effects. We also find that dust masses of 10^8 M_sun are required to observe the dust continuum emission from z 5 galaxies with ALMA.

The Lyman alpha line is a robust tracer of high redshift galaxies. We present estimates of Lyman alpha emission from a protogalactic halo illuminated by UV background radiation fields with various intensities. For this purpose, we performed cosmological hydrodynamics simulations with the adaptive mesh refinement code FLASH, including a detailed network for primordial chemistry,comprising the formation of primordial molecules, a multi-level model for the hydrogen atom as well as the photo-ionization and photo-dissociation processes in a UV background. We find that the presence of a background radiation field J_21 excites the emission of Lyman alpha photons, increasing the Lyman alpha luminosity up to two orders of magnitude. For a halo of \sim 10^10 M_sun, we find that a maximum flux of 5 \times 10^-15 erg cm^-2 s^-1 is obtained for J21 \times f_esc = 0.1, where f_esc is the escape fraction of the ionizing radiation. Depending on the environmental conditions, the flux may vary by three orders of magnitude. For J_21 \times f_esc > 0.1 the Lyman alpha luminosity decreases as the atomic hydrogen abundance becomes rather small. The fluxes derived here can be probed using Subaru and the upcoming James Webb Space Telescope. The emission of Lyman alpha photons is extended and comes from the envelope of the halo rather than its core. In the center of the halo, line trapping becomes effective above columns of 10^22 cm^-2 and suppresses the emission of Lyman alpha. In addition, cooling by primordial molecules may decrease the gas temperature in the central region, which further reduces Lyman alpha emission. In the central core, H_2 is photo-dissociated for a background flux of J_21 \geq 1000. For weaker radiation fields, i.e. J_21 < 0.1, H_2 and HD cooling are particularly strong in the center of the halo, leading to gas temperatures as low as \sim 100 K.

The formation of the first galaxies is accompanied by large accretion flows and virialization shocks, during which the gas is shock-heated to temperatures of $\sim10^4$ K, leading to potentially strong fluxes in the Lyman alpha line. Indeed, a number of Lyman alpha blobs has been detected at high redshift. In this letter, we explore the origin of such Lyman alpha emission using cosmological hydrodynamical simulations that include a detailed model of atomic hydrogen as a multi-level atom and the effects of line trapping with the adaptive mesh refinement code FLASH. We see that baryons fall into the center of a halo through cold streams of gas, giving rise to a Lyman alpha luminosity of at least $\rm 10^{44} erg s^{-1}$ at $\rm z=4.7$, similar to observed Lyman alpha blobs. We find that a Lyman alpha flux of $\rm 5.0\times 10^{-17} erg cm^{-2} s^{-1}$ emerges from the envelope of the halo rather than its center, where the photons are efficiently trapped. Such emission can be probed in detail with the upcoming James Webb Space Telescope (JWST) and will constitute an important probe of gas infall and accretion.

Numerous cosmological simulations have been performed to study the formation of the first objects. We present the results of high resolution 3-D cosmological simulations of primordial objects formation using the adaptive mesh refinement code FLASH by including in an approximate manner the radiative transfer effects of Lyman alpha photons. We compare the results of a Lyman alpha trapping case inside gas clouds with atomic and molecular hydrogen cooling cases.The principal objective of this research is to follow the collapse of a zero metallicity halo with an effective equation of state (that accounts for the trapping) and to explore the fate of a halo in each of the three cases, specifically, the impact of thermodynamics on fragmentation of halos.Our results show that in the case of Lyman alpha trapping, fragmentation is halted and a massive object is formed at the center of a halo. The temperature of the gas remains well above $10^{4}$ K and the halo is not able to fragment to stellar masses. In the atomic cooling case, gas collapses into one or two massive clumps in contrast to the Lyman alpha trapping case. For the molecular hydrogen cooling case, gas cools efficiently and fragments.The formation of massive primordial objects is thus strongly dependent on the thermodynamics of the gas. A salient feature of our results is that for the formation of massive objects, e.g. intermediate mass black holes, feedback effects are not required to suppress $H_{2}$ cooling, as molecular hydrogen is collisionally dissociated at temperatures higher than $10^{4}$ K as a consequence of Lyman alpha trapping.