The formation of supermassive stars has generally been studied under the assumption of rapid accretion of pristine metal-free gas. Recently it was found, however, that gas enriched to metallicities up to $Z \sim 10^{-3}$ Z$_{\odot}$ can also facilitate supermassive star formation, as long as the total mass infall rate onto the protostar remains sufficiently high. We extend the analysis further by examining how the abundance of supermassive star candidate haloes would be affected if all haloes with super-critical infall rates, regardless of metallicity were included. We investigate this scenario by identifying all atomic cooling haloes in the Renaissance simulations with central mass infall rates exceeding a fixed threshold. We find that among these haloes with central mass infall rates above 0.1 M$_{\odot}$ yr$^{-1}$ approximately two-thirds of these haloes have metallicities of $Z > 10^{-3}$ Z$_{\odot}$. If metal mixing within these haloes is inefficient early in their assembly and pockets of metal-poor gas can remain then the number of haloes hosting supermassive stars can be increased by at least a factor of four. Additionally the centres of these high infall-rate haloes provide ideal environments in which to grow pre-existing black holes. Further research into the (supermassive) star formation dynamics of rapidly collapsing haloes, with inhomogeneous metal distributions, is required to gain more insight into both supermassive star formation in early galaxies as well as early black hole growth.

Using the Renaissance suite of simulations we examine the emergence of pristine atomic cooling haloes that are both metal-free and star-free in the early Universe. The absence of metals prevents catastrophic cooling, suppresses fragmentation, and may allow for the formation of massive black hole seeds. Here we report on the abundance of pristine atomic cooling haloes found and on the specific physical conditions that allow for the formation of these direct-collapse-black-hole (DCBH) haloes. In total in our simulations we find that 79 DCBH haloes form before a redshift of 11.6. We find that the formation of pristine atomic haloes is driven by the rapid assembly of the atomic cooling haloes with mergers, both minor and/or major, prior to reaching the atomic cooling limit a requirement. However, the ability of assembling haloes to remain free of (external) metal enrichment is equally important and underlines the necessity of following the transport of metals in such simulations. The candidate DCBH hosting haloes we find, have been exposed to mean Lyman-Werner radiation fields of J$_{LW}$ $\sim$ 1 J$_{21}$ and typically lie at least 10 kpc (physical) from the nearest massive galaxy. Growth rates of the haloes reach values of greater than 10$^7$ M$_{\odot}$ per unit redshift, leading to significant dynamical heating and the suppression of efficient cooling until the halo crosses the atomic cooling threshold. Finally, we also find five synchronised halo candidates where pairs of pristine atomic cooling haloes emerge that are both spatially and temporally synchronised.

We propose a space-based interferometer surveying the gravitational wave (GW) sky in the milli-Hz to $\mu$-Hz frequency range. By the 2040s', the $\mu$-Hz frequency band, bracketed in between the Laser Interferometer Space Antenna (LISA) and pulsar timing arrays, will constitute the largest gap in the coverage of the astrophysically relevant GW spectrum. Yet many outstanding questions related to astrophysics and cosmology are best answered by GW observations in this band. We show that a $\mu$-Hz GW detector will be a truly overarching observatory for the scientific community at large, greatly extending the potential of LISA. Conceived to detect massive black hole binaries from their early inspiral with high signal-to-noise ratio, and low-frequency stellar binaries in the Galaxy, this instrument will be a cornerstone for multimessenger astronomy from the solar neighbourhood to the high-redshift Universe.

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."

Super-Eddington accretion onto massive black hole seeds may be commonplace in the early Universe, where the conditions exist for rapid accretion. Direct collapse black holes are often invoked as a possible solution to the observation of super massive black holes (SMBHs) in the pre-reionisation Universe. We investigate here how feedback, mainly in the form of bipolar jets, from super-Eddington accreting seed black holes will affect their subsequent growth. We find that, nearly independent of the mass loading of the bipolar jets, the violent outflows generated by the jets evacuate a region of approximately 0.1 pc surrounding the black hole seed. However, the jet outflows are unable to break free of the halo and their impact is limited to the immediate vicinity of the black hole. The outflows suppress any accretion for approximately a dynamical time. The gas then cools, recombines and falls back to the centre where high accretion rates are again observed. The overall effect is to create an effective accretion rate with values of between 0.1 and 0.5 times the Eddington rate. If this episodic accretion rate is maintained for order 500 million years then the black hole will increase in mass by a factor of between 3 and 300 but far short of the factor of $10^4$ required for the seeds to become the SMBHs observed at $z>6$. Therefore, direct collapse black holes born into atomic cooling haloes and which experience strong negative mechanical feedback will require external influences (e.g. rapid major mergers with other haloes) to promote efficient accretion and reach SMBH masses within a few hundred million years.

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 origin of supermassive black holes (SMBHs) that inhabit the centers of massive galaxies is largely unconstrained. Remnants from supermassive stars (SMSs) with masses around 10,000 solar masses provide the ideal seed candidates, known as direct collapse black holes. However, their very existence and formation environment in the early Universe are still under debate, with their supposed rarity further exacerbating the problem of modeling their ab-initio formation. SMS models have shown that rapid collapse, with an infall rate above a critical value, in metal-free haloes is a requirement for the formation of a proto-stellar core which will then form an SMS. Using a radiation hydrodynamics simulation of early galaxy formation, we show the natural emergence of metal-free haloes both massive enough, and with sufficiently high infall rates, to form an SMS. We find that haloes that are exposed to both a Lyman-Werner intensity of J_LW ~ 3 J_21 and that undergo at least one period of rapid growth early in their evolution are ideal cradles for SMS formation. This rapid growth induces substantial dynamical heating, amplifying the existing Lyman-Werner suppression originating from a group of young galaxies 20 kiloparsecs away. Our results strongly indicate that structure formation dynamics, rather than a critical Lyman-Werner (LW) flux, may be the main driver of massive black hole formation in the early Universe. We find that massive black hole seeds may be much more common in overdense regions of the early Universe than previously considered with a comoving number density up to 10^-3 Mpc^-3.

The formation of stellar mass black holes from the remnants of Population III stars provides a source of initial black hole seeds with the potential to grow into intermediate or, in rare cases, possibly supermassive black holes. We use the Renaissance simulation suite to follow the growth of over 15,000 black holes born into mini-haloes in the early Universe. We compute the evolution of the black holes by post-processing individual remnant Population III star particles in the Renaissance simulation snapshots. The black holes populate haloes from 10$^{6}$ M$_{\odot}$ up to 10$^{9}$ M$_{\odot}$. We find that all of the black holes display very inefficient growth. On average the black holes increase their initial mass by a factor 10$^{-5}$, with the most active black holes increasing their mass by approximately 10%. Only a single black hole experiences any period of super-Eddington accretion, but the duration is very short and not repeated. Furthermore, we find no correlation of black hole accretion with halo mass in the mass range sampled. Within most haloes, we identify clumps of cool, dense gas for which accretion rates would be high, but instances of black holes encountering these clumps are rare and short-lived. Star formation competes with black hole growth by consuming available gas and driving down accretion rates through feedback. We conclude that the black holes born from Population III remnants do not form a significant population of intermediate mass black holes in the early Universe and will need to wait until later times to undergo significant accretion, if at all.

Recent numerical simulations reveal that the isothermal collapse of pristine gas in atomic cooling haloes may result in stellar binaries of supermassive stars with $M_* \gtrsim 10^4\ \mathrm{M}_{\odot}$. For the first time, we compute the in-situ merger rate for such massive black hole remnants by combining their abundance and multiplicity estimates. For black holes with initial masses in the range $10^{4-6} \ \mathrm{M}_{\odot}$ merging at redshifts $z \gtrsim 15$ our optimistic model predicts that LISA should be able to detect 0.6 mergers per year. This rate of detection can be attributed, without confusion, to the in-situ mergers of seeds from the collapse of very massive stars. Equally, in the case where LISA observes no mergers from heavy seeds at $z \gtrsim 15$ we can constrain the combined number density, multiplicity, and coalesence times of these high-redshift systems. This letter proposes gravitational wave signatures as a means to constrain theoretical models and processes that govern the abundance of massive black hole seeds in the early Universe.

We use high resolution adaptive mesh refinement simulations to model the formation of massive metal-free stars in the early Universe. By applying Lyman-Werner (LW) backgrounds of 100 J$_{21}$ and 1000 J$_{21}$ respectively we construct environments conducive to the formation of massive stars. We find that only in the case of the higher LW backgrounds are super-critical accretion rates realised that are necessary for super-massive star (SMS) formation. Mild fragmentation is observed for both backgrounds. Violent dynamical interactions between the stars that form in the more massive halo formed (1000 J$_{21}$ background) results in the eventual expulsion of the two most massive stars from the halo. In the smaller mass halo (100 J$_{21}$ background) mergers of stars occur before any multibody interactions and a single massive Pop III star is left at the centre of the halo at the end of our simulation. Feedback from the very massive Pop III stars is not effective in generating a large HII region with ionising photons absorbed within a few thousand AU of the star. In all cases a massive black hole seed is the expected final fate of the most massive objects. The seed of the massive Pop III star which remained at the centre of the less massive halo experiences steady accretion rates of almost $10^{-2}$ M$_{\odot}$/yr and if these rates continue could potentially experience super-Eddington accretion rates in the immediate aftermath of collapsing into a black hole.

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 study the influence of a high baryonic streaming velocity on the formation of direct collapse black holes (DCBHs) with the help of cosmological simulations carried out using the moving mesh code {\sc arepo}. We show that a streaming velocity that is as large as three times the root-mean-squared value is effective at suppressing the formation of H$_{2}$-cooled minihaloes, while still allowing larger atomic cooling haloes (ACHs) to form. We find that enough H$_{2}$ forms in the centre of these ACHs to effectively cool the gas, demonstrating that a high streaming velocity by itself cannot produce the conditions required for DCBH formation. However, we argue that high streaming velocity regions do provide an ideal environment for the formation of DCBHs in close pairs of ACHs (the "synchronised halo" model). Due to the absence of star formation in minihaloes, the gas remains chemically pristine until the ACHs form. If two such haloes form with only a small separation in time and space, then the one forming stars earlier can provide enough ultraviolet radiation to suppress H$_{2}$ cooling in the other, allowing it to collapse to form a DCBH. Baryonic streaming may therefore play a crucial role in the formation of the seeds of the highest redshift quasars.

A near pristine atomic cooling halo close to a star forming galaxy offers a natural pathway for forming massive direct collapse black hole (DCBH) seeds which could be the progenitors of the $z>6$ redshift quasars. The close proximity of the haloes enables a sufficient Lyman-Werner flux to effectively dissociate H$_2$ in the core of the atomic cooling halo. A mild background may also be required to delay star formation in the atomic cooling halo, often attributed to distant background galaxies. In this letter we investigate the impact of metal enrichment from both the background galaxies and the close star forming galaxy under extremely unfavourable conditions such as instantaneous metal mixing. We find that within the time window of DCBH formation, the level of enrichment never exceeds the critical threshold (Z$_{cr} \sim 1 \times 10^{-5} \ \rm Z_{\odot})$, and attains a maximum metallicity of Z $\sim 2 \times 10^{-6} \ \rm Z_{\odot}$. As the system evolves, the metallicity eventually exceeds the critical threshold, long after the DCBH has formed.

The Direct Collapse Black Hole (DCBH) scenario provides a solution for forming the massive black holes powering bright quasars observed in the early Universe. A prerequisite for forming a DCBH is that the formation of (much less massive) Population III stars be avoided - this can be achieved by destroying H$_2$ via Lyman-Werner (LW) radiation (E$_{\rm{LW}}$ = 12.6 eV). We find that two conditions must be met in the proto-galaxy that will host the DCBH. First, prior star formation must be delayed; this can be achieved with a background LW flux of J$_{\rm BG} \gtrsim 100\ J_{21}$. Second, an intense burst of LW radiation from a neighbouring star-bursting proto-galaxy is required, just before the gas cloud undergoes gravitational collapse, to finally suppress star formation completely. We show here for the first time using high-resolution hydrodynamical simulations, including full radiative transfer, that this low-level background, combined with tight synchronisation and irradiation of a secondary proto-galaxy by a primary proto-galaxy, inevitably moves the secondary proto-galaxy onto the isothermal atomic cooling track, without the deleterious effects of either photo-evaporating the gas or polluting it by heavy elements. These, atomically cooled, massive proto-galaxies are expected to ultimately form a DCBH of mass $10^4 - 10^5 M_{\odot}$.

We present the Grackle chemistry and cooling library for astrophysical simulations and models. Grackle provides a treatment of non-equilibrium primordial chemistry and cooling for H, D, and He species, including H2 formation on dust grains; tabulated primordial and metal cooling; multiple UV background models; and support for radiation transfer and arbitrary heat sources. The library has an easily implementable interface for simulation codes written in C, C++, and Fortran as well as a Python interface with added convenience functions for semi-analytical models. As an open-source project, Grackle provides a community resource for accessing and disseminating astrochemical data and numerical methods. We present the full details of the core functionality, the simulation and Python interfaces, testing infrastructure, performance, and range of applicability. Grackle is a fully open-source project and new contributions are welcome.

We introduce a new set of large scale, high resolution hydrodynamical simulations of the intergalactic medium: the Sherwood simulation suite. These are performed in volumes 10^3-160^3 h^-3 comoving Mpc^{3}, span almost four orders of magnitude in mass resolution with up to 17.2 billion particles, and employ a variety of physics variations including warm dark matter and galactic outflows. We undertake a detailed comparison of the simulations to high resolution, high signal-to-noise observations of the Lyman-alpha forest over the redshift range 2 < z < 5. The simulations are in very good agreement with the observational data, lending further support to the paradigm that the Lyman-alpha forest is a natural consequence of the web-like distribution of matter arising in LCDM cosmological models. Only a small number of minor discrepancies remain with respect to the observational data. Saturated Lyman-alpha absorption lines with column densities N_HI > 10^14.5 cm^-2 at 2 < z < 2.5 are underpredicted in the models. An uncertain correction for continuum placement bias is required to match the distribution and power spectrum of the transmitted flux, particularly at z > 4. Finally, the temperature of intergalactic gas in the simulations may be slightly too low at z=2.7 and a flatter temperature-density relation is required at z=2.4, consistent with the expected effects of non-equilibrium ionisation during He-II reionisation.

We estimate the expected event rate of gravitational wave signals from mergers of supermassive black holes that could be resolved by a space-based interferometer, such as the Evolved Laser Interferometer Space Antenna (eLISA), utilising the reference cosmological hydrodynamical simulation from the EAGLE suite. These simulations assume a $\Lambda$CDM cosmogony with state-of-the-art subgrid models for radiative cooling, star formation, stellar mass loss, and feedback from stars and accreting black holes. They have been shown to reproduce the observed galaxy population with unprecedented fidelity. We combine the merger rates of supermassive black holes in EAGLE with the latest phenomenological waveform models to calculate the gravitational waves signals from the intrinsic parameters of the merging black holes. The EAGLE models predict $\sim2$ detections per year by a gravitational wave detector such as eLISA. We find that these signals are largely dominated by mergers between seed mass black holes merging at redshifts between $z\sim2$ and $z\sim1$. In order to investigate the dependence on the assumed black hole seed mass, we introduce an additional model with a black hole seed mass an order of magnitude smaller than in our reference model. We also consider a variation of the reference model where a prescription for the expected delays in the black hole merger timescale has been included after their host galaxies merge. We find that the merger rate is similar in all models, but that the initial black hole seed mass could be distinguished through their detected gravitational waveforms. Hence, the characteristic gravitational wave signals detected by eLISA will provide profound insight into the origin of supermassive black holes and the initial mass distribution of black hole seeds.

A nearby source of Lyman-Werner (LW) photons is thought to be a central component in dissociating H$_2$ and allowing for the formation of a direct collapse black hole seed. Nearby sources are also expected to produce copious amounts of hydrogen ionising photons and X-ray photons. We study here the feedback effects of the X-ray photons by including a spectrum due to high-mass X-ray binaries on top of a galaxy with a stellar spectrum. We explicitly trace photon packages emerging from the nearby source and track the radiative and chemical effects of the multi-frequency source $(E_{\rm photon} = \rm{0.76\ eV \rightarrow 7500\ eV}$). We find that X-rays have a strongly negative feedback effect, compared to a stellar only source, when the radiative source is placed at a separation greater than $\gtrsim 1 \ \rm kpc$. The X-rays heat the low and medium density gas in the envelope surrounding the collapsing halo suppressing the mass inflow. The result is a smaller enclosed mass compared to the stellar only case. However, for separations of $\lesssim 1 \ \rm kpc$, the feedback effects of the X-rays becomes somewhat neutral. The enhanced LW intensity at close separations dissociates more H$_2$ and this gas is heated due to stellar photons alone, the addition of X-rays is then not significant. This distance dependence of X-ray feedback suggests that a Goldilocks zone exists close to a forming galaxy where X-ray photons have a much smaller negative feedback effect and ideal conditions exist for creating massive black hole seeds.

The photo-dissociation of H$_2$ by a nearby anisotropic source of radiation is seen as a critical component in creating an environment in which a direct collapse black hole may form. Employing radiative transfer we model the effect of multi-frequency (0.76 eV - 60 eV) radiation on a collapsing halo at high redshift. We vary both the shape of the spectrum which emits the radiation and the distance to the emitting galaxy. We use blackbody spectra with temperatures of $\rm{T = 10^4\ K}$ and $\rm{T = 10^5\ K}$ and a realistic stellar spectrum. We find that an optimal zone exists between 1 kpc and 4 kpc from the emitting galaxy. If the halo resides too close to the emitting galaxy the photo-ionising radiation creates a large HII region which effectively disrupts the collapsing halo, too far from the source and the radiation flux drops below the level of the expected background and the H$_2$ fraction remains too high. When the emitting galaxy is initially placed between 1 kpc and 2 kpc from the collapsing halo, with a spectral shape consistent with a star-forming high redshift galaxy, then a large central core forms. The mass of the central core is between 5000 and 10000 $\rm{M_{\odot}}$ at a temperature of approximately 1000 K. This core is however surrounded by a reservoir of hotter gas at approximately 8000 K which leads to mass inflow rates of the order of $\sim 0.1$ $\rm{M_{\odot}}$ yr$^{-1}$. This environment has the potential to form a massive primordial star which can then lead to the formation of a direct collapse black hole.

We present the results of two 2.3 micron near-infrared radial velocity surveys to detect exoplanets around 36 nearby and young M dwarfs. We use the CSHELL spectrograph (R ~46,000) at the NASA InfraRed Telescope Facility, combined with an isotopic methane absorption gas cell for common optical path relative wavelength calibration. We have developed a sophisticated RV forward modeling code that accounts for fringing and other instrumental artifacts present in the spectra. With a spectral grasp of only 5 nm, we are able to reach long-term radial velocity dispersions of ~20-30 m/s on our survey targets.

We examine the impact of dark matter particle resolution on the formation of a baryonic core in high resolution adaptive mesh refinement simulations. We test the effect that both particle smoothing and particle splitting have on the hydrodynamic properties of a collapsing halo at high redshift (z > 20). Furthermore, we vary the background field intensity, with energy below the Lyman limit (< 13.6 eV), as may be relevant for the case of metal-free star formation and super-massive black hole seed formation. We find that using particle splitting methods greatly increases our particle resolution without introducing any numerical noise and allows us to achieve converged results over a wide range of external background fields. Additionally, we find that for lower values of the background field a lower dark matter particle mass is required. We define the radius of the core as the point at which the enclosed baryonic mass dominates over the enclosed dark matter mass. For our simulations this results in $\rm{R_{core} \sim 5\ pc}$. We find that in order to produce converged results which are not affected by dark matter particles requires that the relationship ${M_{\rm{core}} / M_{\rm{DM}}} > 100.0$ be satisfied, where ${M_{\rm{core}}}$ is the enclosed baryon mass within the core and $M_{\rm{DM}}$ is the minimum dark matter particle mass. This ratio should provide a very useful starting point for conducting convergence tests before any production run simulations. We find that dark matter particle smoothing is a useful adjunct to already highly resolved simulations.

The Spitzer Survey of Stellar Structure in Galaxies (S$^4$G, Sheth et. al. 2010) is a deep 3.6 and 4.5 $\mu$m imaging survey of 2352 nearby ($< 40$ Mpc) galaxies. We describe the S$^4$G data analysis pipeline 4, which is dedicated to 2-dimensional structural surface brightness decompositions of 3.6 $\mu$m images, using GALFIT3.0 \citep{peng2010}. Besides automatic 1-component S\'ersic fits, and 2-component S\'ersic bulge + exponential disk fits, we present human supervised multi-component decompositions, which include, when judged appropriate, a central point source, bulge, disk, and bar components. Comparison of the fitted parameters indicates that multi-component models are needed to obtain reliable estimates for the bulge S\'ersic index and bulge-to-total light ratio ($B/T$), confirming earlier results \citep{laurikainen2007, gadotti2008, weinzirl2009}. In this first paper, we describe the preparations of input data done for decompositions, give examples of our decomposition strategy, and describe the data products released via IRSA and via our web page ({\bf \tt www.oulu.fi/astronomy/S4G\_PIPELINE4/MAIN}). These products include all the input data and decomposition files in electronic form, making it easy to extend the decompositions to suit specific science purposes. We also provide our IDL-based visualization tools (GALFIDL) developed for displaying/running GALFIT-decompositions, as well as our mask editing procedure (MASK\_EDIT) used in data preparation. In the second paper we will present a detailed analysis of the bulge, disk, and bar parameter derived from multi-component decompositions.

The direct collapse model of supermassive black hole seed formation provides an attractive solution to the origin of the quasars now routinely observed at $z \gtrsim 6$. We use the adaptive mesh refinement code Enzo to simulate the collapse of gas at high redshift, including a nine species chemical model of H, He, and H$_2$. The direct collapse model requires that the gas cools predominantly via atomic hydrogen. To this end we simulate the effect of an anisotropic radiation source on the collapse of a halo at high redshift. The radiation source is placed at a distance of 3 kpc (physical) from the collapsing object. The source is set to emit monochromatically in the center of the Lyman-Werner (LW) band only at $12.8 \ \rm{eV}$. The LW radiation emitted from the high redshift source is followed self-consistently using ray tracing techniques. We find that, due to self-shielding, a small amount of H$_2$ is able to form at the very center of the collapsing halo even under very strong LW radiation. Furthermore, we find that a radiation source, emitting $> 10^{54}\ (\sim10^3\ \rm{J_{21}})$ photons per second is required to cause the collapse of a clump of $\rm{M \sim 10^5}$ M$_{\odot}$. The resulting accretion rate onto the collapsing object is $\sim 0.25$ M$_{\odot}$ $\rm{yr^{-1}}$. Our results display significant differences, compared to the isotropic radiation field case, in terms of H$_2$ fraction at an equivalent radius. These differences will significantly effect the dynamics of the collapse. With the inclusion of a strong anisotropic radiation source, the final mass of the collapsing object is found to be $\rm{M \sim 10^5}$ M$_{\odot}$. This is consistent with predictions for the formation of a supermassive star or quasi-star leading to a supermassive black hole.

We have performed high-resolution numerical simulations with the hydrodynamical AMR code Enzo to investigate the formation of massive seed black holes in a sample of six dark matter haloes above the atomic cooling threshold. The aim of this study is to illustrate the effects of varying the maximum refinement level on the final object formed. The virial temperatures of the simulated haloes range from $\rm{T} \sim 10000\ \rm{K} - 16000\ \rm{K}$ and they have virial masses in the range $\rm{M} \sim 2 \times 10^7 \rm{M_{\odot}}$ to $\rm{M} \sim 7 \times 10^7 \rm{M_{\odot}}$ at $z \sim 15$. The outcome of our six fiducial simulations is both generic and robust. A rotationally supported, marginally gravitationally stable, disk forms with an exponential profile. The mass and scale length of this disk depends strongly on the maximum refinement level used. Varying the maximum refinement level by factors between 1 / 64 to 256 times the fiducial level illustrates the care that must be taken in interpreting the results. The lower resolution simulations show tentative evidence that the gas may become rotationally supported out to 20 pc while the highest resolution simulations show only weak evidence of rotational support due to the shorter dynamical times for which the simulation runs. The higher resolution simulations do, however, point to fragmentation at small scales of the order of $\sim 100$ AU. In the highest resolution simulations a central object of a few times $10^2\ \rm{M_{\odot}}$ forms with multiple strongly bound, Jeans unstable, clumps of $\sim 10\ \rm{M_{\odot}}$ and radii of 10 - 20 AU suggesting the formation of dense star clusters in these haloes.

We use a set of semi-numerical simulations based on Zel'dovich approximation, friends-of-friends algorithm and excursion set formalism to generate reionization maps of high dynamic range with a range of assumptions regarding the distribution and luminosity of ionizing sources and the spatial distribution of sinks for the ionizing radiation. We find that ignoring the inhomogeneous spatial distribution of regions of high gas density where recombinations are important -- as is often done in studies of this kind -- can lead to misleading conclusions regarding the topology of reionization, especially if reionization occurs in the photon-starved regime suggested by Lya forest data. The inhomogeneous spatial distribution of recombinations significantly reduces the mean free path of ionizing photons and the typical size of coherently ionized regions. Reionization proceeds then much more as an outside-in process. Low-density regions far from ionizing sources become ionized before regions of high gas density not hosting sources of ionizing radiation. The spatial distribution of sinks of ionization radiation also significantly affects shape and amplitude the power spectrum of fluctuations of 21cm emission. The slope of the 21cm power spectrum as measured by upcoming 21cm experiments should be able to distinguish to what extent the topology of reionization proceeds outside-in or inside-out while the evolution of the amplitude of the power spectrum with increasing ionized mass fraction should be sensitive to the spatial distribution and the luminosity of ionizing sources.

Large dynamic range numerical simulations of atomic cooling driven collapse of gas in pre-galactic DM haloes with T_vir ~ 10000 K show that the gas loses 90% and more of its angular momentum before rotational support sets in. In a fraction of these haloes where the metallicity is low and UV radiation suppresses H_2 cooling, conditions are thus very favourable for the rapid build-up of massive black holes. Depending on the progression of metal enrichment, the continued suppression of H_2 cooling by external and internal UV radiation and the ability to trap the entropy produced by the release of gravitational energy, the gas at the centre of the halo is expected to form a supermassive star, a stellar-mass black hole accreting at super-Eddington accretion rates or a compact star-cluster undergoing collisional run-away of massive stars at its centre. In all three cases a massive black hole of initially modest mass finds itself at the center of a rapid inflow of gas with inflow rates of ~ 1 M_solar\yr. The massive black hole will thus grow quickly to a mass of 10^5 to 10^6 M_solar until further inflow is halted either by consumption of gas by star formation or by the increasing energy and momentum feedback from the growing massive black hole. Conditions for the formation of massive seed black holes in this way are most favourable in haloes with T_vir ~ 15000 K and V_vir ~ 20 km\s with less massive haloes not allowing collapse of gas by atomic cooling and more massive haloes being more prone to fragmentation. This should imprint a characteristic mass on the mass spectrum of an early population of massive black hole seeds in pre-galactic haloes which will later grow into the observed population of supermassive black holes in galactic bulges.

We have used the hydrodynamical AMR code ENZO to investigate the dynamical evolution of the gas at the centre of dark matter haloes with virial velocities of ~ 20 - 30 kms and virial temperatures of ~ 13000-30000 K at z ~ 15 in a cosmological context. The virial temperature of the dark matter haloes is above the threshold where atomic cooling by hydrogen allows the gas to cool and collapse. We neglect cooling by molecular hydrogen and metals, as may be plausible if H_2 cooling is suppressed by a meta-galactic Lyman-Werner background or an internal source of Lyman-Werner photons, and metal enrichment has not progressed very far. The gas in the haloes becomes gravitationally unstable and develops turbulent velocities comparable to the virial velocities of the dark matter haloes. Within a few dynamical times it settles into a nearly isothermal density profile over many decades in radius losing most of its angular momentum in the process. About 0.1 - 1 % of the baryons, at the centre of the dark matter haloes, collapse into a self-gravitating, fat, ellipsoidal, centrifugally supported exponential disc with scale-length of ~ 0.075-0.27 pc and rotation velocities of 25-60 kms. We are able to follow the settling of the gas into centrifugal support and the dynamical evolution of the compact disc in each dark matter halo for a few dynamical times. The dynamical evolution of the gas at the centre of the haloes is complex. In one of the haloes the gas at the centre fragments into a triple system leading to strong tidal perturbations and eventually to the in-fall of a secondary smaller clump into the most massive primary clump. The formation of centrifugally supported self-gravitating massive discs is likely to be an important intermediary stage en route to the formation of a massive black hole seed.

We compare simulations of the Lyman-alpha forest performed with two different hydrodynamical codes, Gadget-2 and Enzo. A comparison of the dark matter power spectrum for simulations run with identical initial conditions show differences of 1-3% at the scales relevant for quantitative studies of the Lyman-alpha forest. This allows a meaningful comparison of the effect of the different implementations of the hydrodynamic part of the two codes. Using the same cooling and heating algorithm in both codes the differences in the temperature and the density probability distribution function are of the order of 10%. These differences can be mainly attributed to a slight mismatch in the resolution. The differences are smaller than or equal to the effects of boxsize and resolution on these statistics. Once resolution effects are taken into account the differences in the flux power spectrum - the statistics most widely used for estimating the matter power spectrum and cosmological parameters from Lyman-alpha forest data - are about 5%. This is again smaller than or equal to the effects of boxsize and resolution. Numerical uncertainties due to a particular implementation of solving the hydrodynamic or gravitational equations appear therefore to contribute only moderately to the error budget in estimates of the flux power spectrum from numerical simulations. We further find that the differences in the flux power spectrum for Enzo simulations run with and without adaptive mesh refinement are also of order 5% or smaller. The latter require 10 times less CPU time making the CPU time requirement similar to that of a version of Gadget-2 that is optimised for Lyman-alpha forest simulations.