We study the thermal evolution of UV-irradiated atomic cooling haloes using high-resolution three-dimensional hydrodynamic simulations. We consider the effect of H^- photodetachment by Ly{\alpha} cooling radiation trapped in the optically-thick cores of three such haloes, a process which has not been included in previous simulations. H^- is a precursor of molecular hydrogen, and therefore, its destruction can diminish the H2 abundance and cooling. Using a simple high-end estimate for the trapped Ly{\alpha} energy density, we find that H^- photodetachment by Ly{\alpha} decreases the critical UV flux for suppressing H2-cooling by up to a factor of \approx 5. With a more conservative estimate of the Ly{\alpha} energy density, we find the critical flux is decreased only by ~15-50 percent. Our results suggest that Ly{\alpha} radiation may have an important effect on the thermal evolution of UV-irradiated haloes, and therefore on the potential for massive black hole formation.

Massive 10^6-10^10 Msun black holes (BHs) are ubiquitous in local galactic nuclei. They were common by the time the Universe is several Gyr old, and many of them were in place within the first 1~Gyr after the Big Bang. Their quick assembly has been attributed to mechanisms such as the rapid collapse of gas into the nuclei of early protogalaxies, accretion and mergers of stellar-mass BHs accompanying structure formation at early times, and the runaway collapse of early, ultra-dense stellar clusters. The origin of the early massive BHs remains an intriguing and long-standing unsolved puzzle in astrophysics. Here we discuss strategies for discerning between BH seeding models using electromagnetic observations. We argue that the most direct answers will be obtained through detection of BHs with masses M<10^5 Msun at redshifts z>10, where we expect them to first form. Reaching out to these redshifts and down to these masses is crucial, because BHs are expected to lose the memory of their initial assembly by the time they grow well above 10^5 Msun and are incorporated into higher-mass galaxies. The best way to detect 10^4-10^5 Msun BHs at high redshifts is by a sensitive X-ray survey. Critical constraining power is augmented by establishing the properties and the environments of their host galaxies in deep optical/IR imaging surveys. Required OIR data can be obtained with the JWST and WFIRST missions. The required X-ray flux limits (down to 10^{-19} erg/s/cm^2) are accessible only with a next-generation X-ray observatory which has both high (sub-1") angular resolution and high throughput. A combination of deep X-ray and OIR surveys will be capable of probing several generic markers of the BH seed scenarios, and resolving the long-stanging puzzle of their origin. These electromagnetic observations are also highly synergistic with the information from LISA on high-z BH mergers.

The abundance of molecular hydrogen (H2), the primary coolant in primordial gas, is critical for the thermodynamic evolution and star-formation histories in early protogalaxies. Determining the photodissociation rate of H2 by an incident Lyman-Werner (LW) flux is thus crucial, but prohibitively expensive to calculate on the fly in simulations. The rate is sensitive to the H2 rovibrational distribution, which in turn depends on the gas density, temperature, and incident LW radiation field. We use the publicly available cloudy package to model primordial gas clouds and compare exact photodissociation rate calculations to commonly-used fitting formulae. We find the fit from Wolcott-Green et al. (2011) is most accurate for moderate densities n~10^3 cm^{-3} and temperatures, T~10^3K, and we provide a new fit, which captures the increase in the rate at higher densities and temperatures, owing to the increased excited rovibrational populations in this regime. Our new fit has typical errors of a few percent percent up to n =<10^7 cm^{-3}, T =< 8000K, and H2 column density NH2 =<10^{17} cm^{-2}, and can be easily utilized in simulations. We also show that pumping of the excited rovibrational states of H2 by a strong LW flux further modifies the level populations when the gas density is low, and noticeably decreases self-shielding for J_21 > 10^3 and n < 10^2 cm^{-3}. This may lower the "critical flux" at which primordial gas remains H2-poor in some protogalaxies, enabling massive black hole seed formation.

Suppression of H2-cooling in early protogalaxies has important implications for the formation of supermassive black holes seeds, the first generation of stars, and the epoch of reionization. This suppression can occur via photodissociation of H2 (by ultraviolet Lyman-Werner [LW] photons) or by photodetachment of H-, a precursor in H2 formation (by infrared [IR] photons). Previous studies have typically adopted idealised spectra, with a blackbody or a power-law shape, in modeling the chemistry of metal-free protogalaxies, and utilised a single parameter, the critical UV flux, or Jcrit, to determine whether H2-cooling is prevented. This can be misleading, and that independent of the spectral shape, there is a critical curve in the (kLW , kH^- ) plane, where kLW and kH^- are the H2-dissocation rates by LW and IR photons, which determines whether a protogalaxy can cool below ~1000 Kelvin. We use a one-zone model to follow the chemical and thermal evolution of gravitationally collapsing protogalactic gas, to compute this critical curve, and provide an accurate analytical fit for it. We improve on previous works by considering a variety of more realistic Pop III or Pop II-type spectra from population synthesis models and perform fully frequency-dependent calculations of the H2-photodissociation rates for each spectrum. We compute the ratio kLW/kH^- for each spectrum, as well as the minimum stellar mass M_star, for various IMFs and metallicities, required to prevent cooling in a neighboring halo a distance d away. We provide critical M_star/d2 values for suppression of H2-cooling, with analytic fits, which can be used in future studies.

It is commonly believed that the earliest stages of star-formation in the Universe were self-regulated by global radiation backgrounds - either by the ultraviolet Lyman-Werner (LW) photons emitted by the first stars (directly photodissociating H_2), or by the X-rays produced by accretion onto the black hole (BH) remnants of these stars (heating the gas but catalyzing H_2 formation). Recent studies have suggested that a significant fraction of the first stars may have had low masses (a few M_sun). Such stars do not leave BH remnants and they have softer spectra, with copious infrared (IR) radiation at photon energies around 1eV. Similar to LW and X-ray photons, these photons have a mean-free path comparable to the Hubble distance, building up an early IR background. Here we show that if soft-spectrum stars, with masses of a few M_sun, contributed more than 1% of the UV background (or their mass fraction exceeded 90%), then their IR radiation dominated radiative feedback in the early Universe. The feedback is different from the UV feedback from high-mass stars, and occurs through the photo-detachment of H^- ions, necessary for efficient H_2 formation. Nevertheless, we find that the baryon fraction which must be incorporated into low-mass stars in order to suppress H_2-cooling is only a factor of few higher than for high-mass stars.

The ability of primordial gas to cool in proto-galactic haloes exposed to Lyman-Werner (LW) radiation is critically dependent on the self-shielding of H_2. We perform radiative transfer calculations of LW line photons, post-processing outputs from three-dimensional adaptive mesh refinement (AMR) simulations of haloes with T_vir > 10^4 K at redshifts around z=10. We calculate the optically thick photodissociation rate numerically, including the effects of density, temperature, and velocity gradients in the gas, as well as line overlap and shielding of H_2 by HI, over a large number of sight-lines. In low-density regions (n<10^4 cm^-3) the dissociation rates exceed those obtained using most previous approximations by more than an order of magnitude; the correction is smaller at higher densities. We trace the origin of the deviations primarily to inaccuracies of (i) the most common fitting formula (Draine & Bertoldi 1996) for the suppression of the dissociation rate and (ii) estimates for the effective shielding column density from local properties of the gas. The combined effects of gas temperature and velocity gradients are comparatively less important, typically altering the spherically averaged rate only by a factor of less than two. We present a simple modification to the DB96 fitting formula for the optically thick rate which improves agreement with our numerical results to within approx. 15 per cent, and can be adopted in future simulations. We find that estimates for the effective shielding column can be improved by using the local Sobolev length. Our correction to the H_2 self-shielding reduces the critical LW flux to suppress H_2-cooling in T_vir>10^4 K haloes by an order of magnitude; this increases the number of such haloes in which supermassive (approx. M=10^5 M_sun) black holes may have formed.

It has been shown that HD molecules can form efficiently in metal-free gas collapsing into massive protogalactic halos at high redshift. The resulting radiative cooling by HD can lower the gas temperature to that of the cosmic microwave background, T_CMB=2.7(1+z)K, significantly below the temperature of a few 100 K achievable via H_2-cooling alone, and thus reduce the masses of the first generation of stars. Here we consider the suppression of HD-cooling by UV irradiation in the Lyman-Werner (LW) bands. We include photo-dissociation of both H_2 and HD, and explicitly compute the self-shielding and shielding of both molecules by neutral hydrogen as well as the shielding of HD by H_2. We use a simplified dynamical collapse model, and follow the chemical and thermal evolution of the gas, in the presence of a UV background. We find that a LW flux of J_crit = 1e-22 erg/cm^2/sr/s/Hz is able to suppress HD cooling and thus prevent collapsing primordial gas from reaching temperatures below 100 K. The main reason for the lack of HD cooling for J>J_crit is the partial photo-dissociation of H_2, which prevents the gas from reaching sufficiently low temperatures (T<150K) for HD to become the dominant coolant; direct HD photo-dissociation is unimportant except for a narrow range of fluxes and column densities. Since the prevention of HD-cooling requires only partial H_2 photo-dissociation, the critical flux J_crit is modest, and is below the UV background required to reionize the universe at redshift z=10-20. We conclude that HD-cooling can reduce the masses of typical stars only in rare halos forming well before the epoch of reionization.