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.

We present Powderday, a flexible, fast, open-source dust radiative transfer package designed to interface with galaxy formation simulations. Powderday builds on FSPS population synthesis models, Hyperion dust radiative transfer, and employs yt to interface between different software packages. We include our stellar population synthesis modeling on the fly, which allows for significant run-time flexibility in the assumed stellar physics. We include a model for nebular line emission that can employ either precomputed Cloudy lookup tables (for efficiency), or direct photoionization calculations for all young stars (for flexibility). The dust content follows either observationally-motivated prescriptions, direct modeling from galaxy formation simulations, or a novel approach that includes the dust content via learning-based algorithms from the SIMBA cosmological galaxy formation simulation. AGN can additionally be included via a range of prescriptions. The output of these models are broadband SEDs, as well as filter-convolved images. Powderday is designed to eliminate last-mile efforts by researchers that employ different hydrodynamic galaxy formation models, and seamlessly interfaces with GIZMO, AREPO, GASOLINE, CHANGA, and ENZO. We demonstrate the capabilities of the code via three applications: a model for the star formation rate (SFR) - infrared luminosity relation in galaxies (including the impact of AGN); the impact of circumstellar dust around AGB stars on the mid-infrared emission from galaxy SEDs; and the impact of galaxy inclination angle on dust attenuation laws.

Hyper-accretion occurs when the gas inflow rate onto a black hole (BH) is so high that the radiative feedback cannot reverse the accretion flow. This extreme process is a promising mechanism for the rapid growth of seed BHs in the early universe, which can explain high-redshift quasars powered by billion solar mass BHs. In theoretical models, spherical symmetry is commonly adopted for hyper-accretion flow, however its sustainability on timescales corresponding to the BH growth has not been addressed yet. Here we show that stochastic interactions between the ionizing radiation from the BH and non-uniform accretion flow can lead to a formation of rotating gas disk around the BH. Once the disk is formed, the supply of gas to the BH preferentially occurs via biconical dominated accretion flow (BDAF) perpendicular to the disk, avoiding the centrifugal barrier of the disk. BDAFs from opposite directions collide in the vicinity of the BH supplying high-density, low angular momentum gas to the BH, whereas most of the gas with non-negligible angular momentum is deflected to the rotationally supported outflowing decretion disk. The disk becomes reinforced progressively as more mass from the biconical flow transfers to the disk and some of the outflowing gas from the disk is redirected to the biconical accretion funnels through a meridional structure. This axisymmetric hydrodynamic structure of a BDAF and decretion disk continues to provide uninterrupted flow of high-density gas to the BH.

Recent observations provide evidence that some cool-core clusters (CCCs) host quasars in their brightest cluster galaxies (BCGs). Motivated by these findings we use 3D radiation-hydrodynamic simulations with the code Enzo to explore the joint role of the kinetic and radiative feedback from supermassive black holes (SMBHs) in BCGs. We implement kinetic feedback as sub-relativistic plasma outflows and model radiative feedback using the ray-tracing radiative transfer or thermal energy injection. In our simulations the central SMBH transitions between the radiatively efficient and radiatively inefficient states on timescales of a few Gyr, as a function of its accretion rate. The timescale for this transition depends primarily on the fraction of power allocated to each feedback mode, and to a lesser degree on the overall feedback luminosity of the active galactic nucleus (AGN). Specifically, we find that (a) kinetic feedback must be present at both low and high accretion rates in order to prevent the cooling catastrophe, and (b) its contribution likely accounts for > 10% of the total AGN feedback power, since below this threshold simulated BCGs tend to host radio-loud quasars most of the time, in apparent contrast with observations. We also find a positive correlation between the AGN feedback power and the mass of the cold gas filaments in the cluster core, indicating that observations of H$\alpha$ filaments can be used as a measure of AGN feedback.

The formation of Population III (Pop III) stars is a critical step in the evolution of the early universe. To understand how these stars affected their metal-enriched descendants, the details of how, why and where Pop III formation takes place needs to be determined. One of the processes that is assumed to greatly affect the formation of Pop III stars is the presence of a Lyman-Werner (LW) radiation background, that destroys H$_2$, a necessary coolant in the creation of Pop III stars. Self-shielding can alleviate the effect the LW background has on the H$_2$ within haloes. In this work, we perform a cosmological simulation to study the birthplaces of Pop III stars, using the adaptive mesh refinement code Enzo. We investigate the distribution of host halo masses and its relationship to the LW background intensity. Compared to previous work, haloes form Pop III stars at much lower masses, up to a factor of a few, due to the inclusion of H$_2$ self-shielding. We see no relationship between the LW intensity and host halo mass. Most haloes form multiple Pop III stars, with a median number of four, up to a maximum of 16, at the instance of Pop III formation. Our results suggest that Pop III star formation may be less affected by LW radiation feedback than previously thought and that Pop III multiple systems are common.

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.

As computational resolution of modern cosmological simulations reach ever so close to resolving individual star-forming clumps in a galaxy, a need for "resolution-appropriate" physics for a galaxy-scale simulation has never been greater. To this end, we introduce a self-consistent numerical framework that includes explicit treatments of feedback from star-forming molecular clouds (SFMCs) and massive black holes (MBHs). In addition to the thermal supernovae feedback from SFMC particles, photoionizing radiation from both SFMCs and MBHs is tracked through full 3-dimensional ray tracing. A mechanical feedback channel from MBHs is also considered. Using our framework, we perform a state-of-the-art cosmological simulation of a quasar-host galaxy at z~7.5 for ~25 Myrs with all relevant galactic components such as dark matter, gas, SFMCs, and an embedded MBH seed of ~> 1e6 Ms. We find that feedback from SFMCs and an accreting MBH suppresses runaway star formation locally in the galactic core region. Newly included radiation feedback from SFMCs, combined with feedback from the MBH, helps the MBH grow faster by retaining gas that eventually accretes on to the MBH. Our experiment demonstrates that previously undiscussed types of interplay between gas, SFMCs, and a MBH may hold important clues about the growth and feedback of quasars and their host galaxies in the high-redshift Universe.

The direct formation of a massive black hole is a potential seeding mechanism of the earliest observed supermassive black holes. We investigate how the existence of a massive black hole seed impacts the ionization and thermal state of its pre-galactic host halo and subsequent star formation. We show that its X-ray radiation ionizes and heats the medium, enhancing $\rm{H}_2$ formation in shielded regions, within the nuclear region in the span of a million years. The enhanced molecular cooling triggers the formation of a $\sim 10^4~{\rm M}_\odot$ metal-free stellar cluster at a star formation efficiency of $\sim 0.1\%$ in a single event. Star formation occurs near the edges of the H II region that is partially ionized by X-rays, thus the initial size depends on the black hole properties and surrounding environment. The simulated metal-free galaxy has an initial half-light radius of $\sim 10$ pc but expands to $\sim 50$ pc after 10 million years because of the outward velocities of their birth clouds. Supernova feedback then quenches any further star formation for tens of millions of years, allowing the massive black hole to dominate the spectrum once the massive metal-free stars die.

The universe goes through several phase transitions during its formative stages. Cosmic reionization is the last of them, where ultraviolet and X-ray radiation escape from the first generations of galaxies heating and ionizing their surroundings and subsequently the entire intergalactic medium. There is strong observational evidence that cosmic reionization ended approximately one billion years after the Big Bang, but there are still uncertainties that will be clarified with upcoming optical, infrared, and radio facilities in the next decade. This article gives an introduction to the theoretical and observational aspects of cosmic reionization and discusses their role in our understanding of early galaxy formation and cosmology.

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.

We present simulations from the new "Figuring Out Gas & Galaxies in Enzo" (FOGGIE) project. In contrast to most extant simulations of galaxy formation, which concentrate computational resources on galactic disks and spheroids with fluid and particle elements of fixed mass, the FOGGIE simulations focus on extreme spatial and mass resolution in the circumgalactic medium (CGM) surrounding galaxies. Using the Enzo code and a new refinement scheme, FOGGIE reaches spatial resolutions of 381 comoving $h^{-1}$ pc and resolves extremely low masses ($\lesssim 1$--$100$ Msun out to 100 comoving $h^{-1}$ kpc from the central halo. At these resolutions, cloud and filament-like structures giving rise to simulated absorption are smaller, and better resolved, than the same structures simulated with standard density-dependent refinement. Most of the simulated absorption arises in identifiable and well-resolved structures with masses $\lesssim 10^4$ Msun, well below the mass resolution of typical zoom simulations. However, integrated quantities such as mass surface density and ionic covering fractions change at only the $\lesssim 30$% level as resolution is varied. This relatively small changes in projected quantities---even when the sizes and distribution of absorbing clouds change dramatically---indicate that commonly used observables provide only weak constraints on the physical structure of the underlying gas. Comparing the simulated absorption features to the KODIAQ (Keck Observatory Database of Ionized Absorption toward Quasars) survey of $z \sim2$--$3.5$ Lyman limit systems, we show that high-resolution FOGGIE runs better resolve the internal kinematic structure of detected absorption, and better match the observed distribution of absorber properties. These results indicate that CGM resolution is key in properly testing simulations of galaxy evolution with circumgalactic observations.

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.

Traditional cosmological hydrodynamics simulations fail to spatially resolve the circumgalatic medium (CGM), the reservoir of tenuous gas surrounding a galaxy and extending to its virial radius. We introduce the technique of Enhanced Halo Resolution (EHR), enabling more realistic physical modeling of the simulated CGM by consistently forcing gas refinement to smaller scales throughout the virial halo of a simulated galaxy. We investigate the effects of EHR in the Tempest simulations, a suite of Enzo-based cosmological zoom simulations following the evolution of an L* galaxy, resolving spatial scales of 500 comoving pc out to 100 comoving kpc in galactocentric radius. Among its many effects, EHR (1) changes the thermal balance of the CGM, increasing its cool gas content and decreasing its warm/hot gas content; (2) preserves cool gas structures for longer periods; and (3) enables these cool clouds to exist at progressively smaller size scales. Observationally, this results in a boost in "low ions" like H I and a drop in "high ions" like O VI throughout the CGM. These effects of EHR do not converge in the Tempest simulations, but extrapolating these trends suggests that the CGM in reality is a mist consisting of ubiquitous, small, long-lived, cool clouds suspended in a hot medium at the virial temperature of the halo. Additionally, we explore the physical mechanisms to explain why EHR produces the above effects, proposing that it works both by (1) better sampling the distribution of CGM phases enabling runaway cooling in the denser, cooler tail of the phase distribution; and (2) preventing cool gas clouds from artificially mixing with the ambient hot halo and evaporating. Evidence is found for both EHR mechanisms occurring in the Tempest simulations.

We investigate the formation of extremely metal-poor (EMP) stars that are observed in the Galactic halo and neighboring ultra-faint dwarf galaxies. Their low metal abundances (${\rm [Fe/H]} < -3$) indicate that their parent clouds were enriched by a single or several supernovae (SNe) from the first (Pop III) stars. In this study, we perform numerical simulations of the entire formation sequence of a EMP star through the feedback effects of photo-ionization and metal-enrichment by a Pop III SN. We for the first time employ a metal/dust properties calculated consistently with the progenitor model, and solve all relevant radiative cooling processes and chemical reactions including metal molecular formation and grain growth until the protostar formation. In a minihalo (MH) with mass $1.77\times 10^{6} \ {\rm M}_{\bigodot}$, a Pop III star with mass $13 \ {\rm M}_{\bigodot}$ forms at redshift $z=12.1$. After its SN explosion, the shocked gas falls back into the central MH internally enriching itself. The metallicity in the recollapsing region is $2.6\times 10^{-4} \ {\rm Z}_{\bigodot}$ (${\rm [Fe/H]} = -3.42$). The recollapsing cloud undergoes cooling by HD, CO, and OH molecules and heating along with H$_2$ formation. Eventually by grain growth and dust cooling, knotty filaments appear in the central 100 au region with the help of turbulence driven by the SN, leading to the formation of low-mass EMP stars surviving until the present day.

Previously we identified a new class of early galaxy that we estimate contributes up to 30\% of the ionizing photons responsible for reionization. These are low mass halos in the range $M_h =10^{6.5}-10^{8} M_{\odot}$ that have been chemically enriched by supernova ejecta from prior Pop III star formation. Despite their low star formation rates, these Metal Cooling halos (MCs) are significant sources of ionizing radiation, especially at the onset of reionization, due to their high number density and ionizing escape fractions. Here we present a fully-coupled radiation hydrodynamic simulation of reionization that includes these MCs as well the more massive hydrogen atomic line cooling halos. Our method is novel: we perform halo finding inline with the radiation hydrodynamical simulation, and assign escaping ionizing fluxes to halos using a probability distribution function (PDF) measured from the galaxy-resolving Renaissance Simulations. The PDF captures the mass dependence of the ionizing escape fraction as well as the probability that a halo is actively forming stars. With MCs, reionization starts earlier than if only halos of $10^8 M_{\odot}$ and above are included, however the redshift when reionization completes is only marginally affected as this is driven by more massive galaxies. Because star formation is intermittent in MCs, the earliest phase of reionization exhibits a stochastic nature, with small H II regions forming and recombining. Only later, once halos of mass $\sim 10^9 M_{\odot}$ and above begin to dominate the ionizing emissivity, does reionization proceed smoothly in the usual manner deduced from previous studies. This occurs at $z\approx 10$ in our simulation.

Space telescope observations of massive black holes during their formation may be key to understanding the origin of supermassive black holes and high-redshift quasars. To create diagnostics for their detection and confirmation, we study a simulation of a nascent massive, so-called direct-collapse, black hole that induces a wave of nearby massive metal-free star formation, unique to this seeding scenario and to very high redshifts. Here we describe a series of distinct colors and emission line strengths, dependent on the relative strength of star formation and black hole accretion. We predict that the forthcoming James Webb Space Telescope might be able to detect and distinguish a young galaxy that hosts a direct-collapse black hole in this configuration at redshift 15 with as little as a 20,000-second total exposure time across four filters, critical for constraining supermassive black hole seeding mechanisms and early growth rates. We also find that a massive seed black hole produces strong, H$_2$-dissociating Lyman-Werner radiation.

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.

Regardless of their initial seed mass, any active galactic nuclei observed at redshifts z > 6 must have grown by several orders of magnitude from their seeds. In this chapter, we will discuss the physical processes and latest research on black hole growth and associated feedback after seed formation. Fueling is initially slowed down by radiative feedback from the black hole itself and supernova explosions from nearby stars. Its growth however accelerates once the host galaxy grows past a critical mass.

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.

We have modeled direct collapse of a primordial gas within dark matter halos in the presence of radiative transfer, in high-resolution zoom-in simulations in a cosmological framework, down to the formation of the photosphere and the central object. Radiative transfer has been implemented in the flux-limited diffusion (FLD) approximation. Adiabatic models were run for comparison. We find that (a) the FLD flow forms an irregular central structure and does not exhibit fragmentation, contrary to adiabatic flow which forms a thick disk, driving a pair of spiral shocks, subject to Kelvin-Helmholtz shear instability forming fragments; (b) the growing central core in the FLD flow quickly reaches ~10 Mo and a highly variable luminosity of 10^{38}-10^{39} erg/s, comparable to the Eddington luminosity. It experiences massive recurrent outflows driven by radiation force and thermal pressure gradients, which mix with the accretion flow and transfer the angular momentum outwards; and (c) the interplay between these processes and a massive accretion, results in photosphere at ~10 AU. We conclude that in the FLD model (1) the central object exhibits dynamically insignificant rotation and slower than adiabatic temperature rise with density; (2) does not experience fragmentation leading to star formation, thus promoting the fast track formation of a supermassive black hole (SMBH) seed; (3) inclusion of radiation force leads to outflows, resulting in the mass accumulation within the central 10^{-3} pc, which is ~100 times larger than characteristic scale of star formation. The inclusion of radiative transfer reveals complex early stages of formation and growth of the central structure in the direct collapse scenario of SMBH seed formation.

Direct collapse within dark matter (DM) halos is a promising path to form supermassive black hole (SMBH) seeds at high redshifts. The outer part of this collapse remains optically thin, and has been studied intensively using numerical simulations. However, the innermost region of the collapse is expected to become optically thick and requires us to follow the radiation field in order to understand its subsequent evolution. So far, the adiabatic approximation has been used exclusively for this purpose. We apply radiative transfer in the flux-limited diffusion (FLD) approximation to solve the evolution of coupled gas and radiation, for isolated halos. For direct collapse within isolated DM halos, we find that (1) the photosphere forms at ~10^{-6} pc and rapidly expands outward. (2) A central core forms, with a mass of ~1 Mo, supported by thermal gas pressure gradients and rotation. (3) Growing thermal gas and radiation pressure gradients dissolve it. (4) This process is associated with a strong anisotropic outflow, and another core forms nearby and grows rapidly. (5) Typical radiation luminosity emerging from the photosphere encompassing these cores is ~5 x 10^{37}-5 x 10^{38} erg/s, of order the Eddington luminosity. (6) Two variability timescales are associated with this process: a long one, which is related to the accretion flow within the central ~10^{-4}-10^{-3} pc, and ~0.1 yr, which is related to radiation diffusion. (7) Adiabatic models have been run for comparison and their evolution differs profoundly from that of the FLD models, by forming a central geometrically-thick disk. Overall, an adiabatic equation of state is not a good approximation to the advanced stage of direct collapse, mainly because the radiation in the FLD is capable of escaping due to anisotropy in the optical depth and associated gradients.

We use a cosmological hydrodynamic simulation calculated with Enzo and the semi-analytic galaxy formation model (SAM) GAMMA to address the chemical evolution of dwarf galaxies in the early universe. The long-term goal of the project is to better understand the origin of metal-poor stars and the formation of dwarf galaxies and the Milky Way halo by cross-validating these theoretical approaches. We combine GAMMA with the merger tree of the most massive galaxy found in the hydrodynamic simulation and compare the star formation rate, the metallicity distribution function (MDF), and the age-metallicity relationship predicted by the two approaches. We found that the SAM can reproduce the global trends of the hydrodynamic simulation. However, there are degeneracies between the model parameters and more constraints (e.g., star formation efficiency, gas flows) need to be extracted from the simulation to isolate the correct semi-analytic solution. Stochastic processes such as bursty star formation histories and star formation triggered by supernova explosions cannot be reproduced by the current version of GAMMA. Non-uniform mixing in the galaxy's interstellar medium, coming primarily from self-enrichment by local supernovae, causes a broadening in the MDF that can be emulated in the SAM by convolving its predicted MDF with a Gaussian function having a standard deviation of ~0.2 dex. We found that the most massive galaxy in the simulation retains nearby 100% of its baryonic mass within its virial radius, which is in agreement with what is needed in GAMMA to reproduce the global trends of the simulation.

In hierarchical structure formation, merging of galaxies is frequent and known to dramatically affect their properties. To comprehend these interactions high-resolution simulations are indispensable because of the nonlinear coupling between pc and Mpc scales. To this end, we present the first adaptive mesh refinement (AMR) simulation of two merging, low mass, initially gas-rich galaxies (1.8e10 Ms each), including star formation and feedback. With galaxies resolved by ~2e7 total computational elements, we achieve unprecedented resolution of the multiphase interstellar medium, finding a widespread starburst in the merging galaxies via shock-induced star formation. The high dynamic range of AMR also allows us to follow the interplay between the galaxies and their embedding medium depicting how galactic outflows and a hot metal-rich halo form. These results demonstrate that AMR provides a powerful tool in understanding interacting galaxies.

Using an isolated Milky Way-mass galaxy simulation, we compare results from 9 state-of-the-art gravito-hydrodynamics codes widely used in the numerical community. We utilize the infrastructure we have built for the AGORA High-resolution Galaxy Simulations Comparison Project. This includes the common disk initial conditions, common physics models (e.g., radiative cooling and UV background by the standardized package Grackle) and common analysis toolkit yt, all of which are publicly available. Subgrid physics models such as Jeans pressure floor, star formation, supernova feedback energy, and metal production are carefully constrained across code platforms. With numerical accuracy that resolves the disk scale height, we find that the codes overall agree well with one another in many dimensions including: gas and stellar surface densities, rotation curves, velocity dispersions, density and temperature distribution functions, disk vertical heights, stellar clumps, star formation rates, and Kennicutt-Schmidt relations. Quantities such as velocity dispersions are very robust (agreement within a few tens of percent at all radii) while measures like newly-formed stellar clump mass functions show more significant variation (difference by up to a factor of ~3). Systematic differences exist, for example, between mesh-based and particle-based codes in the low density region, and between more diffusive and less diffusive schemes in the high density tail of the density distribution. Yet intrinsic code differences are generally small compared to the variations in numerical implementations of the common subgrid physics such as supernova feedback. Our experiment reassures that, if adequately designed in accordance with our proposed common parameters, results of a modern high-resolution galaxy formation simulation are more sensitive to input physics than to intrinsic differences in numerical schemes.

We investigate the spatially-resolved star formation relation using a galactic disk formed in a comprehensive high-resolution (3.8 pc) simulation. Our new implementation of stellar feedback includes ionizing radiation as well as supernova explosions, and we handle ionizing radiation by solving the radiative transfer equation rather than by a subgrid model. Photoheating by stellar radiation stabilizes gas against Jeans fragmentation, reducing the star formation rate. Because we have self-consistently calculated the location of ionized gas, we are able to make spatially-resolved mock observations of star formation tracers, such as H-alpha emission. We can also observe how stellar feedback manifests itself in the correlation between ionized and molecular gas. Applying our techniques to the disk in a galactic halo of 2.3e11 Msun, we find that the correlation between star formation rate density (estimated from mock H-alpha emission) and molecular hydrogen density shows large scatter, especially at high resolutions of <~ 75 pc that are comparable to the size of giant molecular clouds (GMCs). This is because an aperture of GMC size captures only particular stages of GMC evolution, and because H-alpha traces hot gas around star-forming regions and is displaced from the molecular hydrogen peaks themselves. By examining the evolving environment around star clusters, we speculate that the breakdown of the traditional star formation laws of the Kennicutt-Schmidt type at small scales is further aided by a combination of stars drifting from their birthplaces, and molecular clouds being dispersed via stellar feedback.

We describe a new method for simulating ionizing radiation and supernova feedback in the analogues of low-redshift galactic disks. In this method, which we call star-forming molecular cloud (SFMC) particles, we use a ray-tracing technique to solve the radiative transfer equation for ultraviolet photons emitted by thousands of distinct particles on the fly. Joined with high numerical resolution of 3.8 pc, the realistic description of stellar feedback helps to self-regulate star formation. This new feedback scheme also enables us to study the escape of ionizing photons from star-forming clumps and from a galaxy, and to examine the evolving environment of star-forming gas clumps. By simulating a galactic disk in a halo of 2.3e11 Msun, we find that the average escape fraction from all radiating sources on the spiral arms (excluding the central 2.5 kpc) fluctuates between 0.08% and 5.9% during a ~20 Myr period with a mean value of 1.1%. The flux of escaped photons from these sources is not strongly beamed, but manifests a large opening angle of more than 60 degree from the galactic pole. Further, we investigate the escape fraction per SFMC particle, f_esc(i), and how it evolves as the particle ages. We discover that the average escape fraction f_esc is dominated by a small number of SFMC particles with high f_esc(i). On average, the escape fraction from a SFMC particle rises from 0.27% at its birth to 2.1% at the end of a particle lifetime, 6 Myrs. This is because SFMC particles drift away from the dense gas clumps in which they were born, and because the gas around the star-forming clumps is dispersed by ionizing radiation and supernova feedback. The framework established in this study brings deeper insight into the physics of photon escape fraction from an individual star-forming clump, and from a galactic disk.

There is mounting evidence for the coevolution of galaxies and their embedded massive black holes (MBHs) in a hierarchical structure formation paradigm. To tackle the nonlinear processes of galaxy-MBH interaction, we describe a self-consistent numerical framework which incorporates both galaxies and MBHs. The high-resolution adaptive mesh refinement (AMR) code Enzo is modified to model the formation and feedback of molecular clouds at their characteristic scale of 15.2 pc and the accretion of gas onto a MBH. Two major channels of MBH feedback, radiative feedback (X-ray photons followed through full 3D adaptive ray tracing) and mechanical feedback (bipolar jets resolved in high-resolution AMR), are employed. We investigate the coevolution of a 9.2e11 Msun galactic halo and its 1e5 Msun embedded MBH at redshift 3 in a cosmological LCDM simulation. The MBH feedback heats the surrounding ISM up to 1e6 K through photoionization and Compton heating and locally suppresses star formation in the galactic inner core. The feedback considerably changes the stellar distribution there. This new channel of feedback from a slowly growing MBH is particularly interesting because it is only locally dominant, and does not require the heating of gas globally on the disk. The MBH also self-regulates its growth by keeping the surrounding ISM hot for an extended period of time.

Ultrafaint dwarf galaxies (UFDs) are typically assumed to have simple, stellar populations with star formation ending at reionization. Yet as the observations of these galaxies continue to improve, their star formation histories (SFHs) are revealed to be more complicated than previously thought. In this paper, we study how star formation, chemical enrichment, and mixing proceed in small, dark matter halos at early times using a high-resolution, cosmological, hydrodynamical simulation. The goals are to inform the future use of analytic models and to explore observable properties of the simulated halos in the context of UFD data. Specifically, we look at analytic approaches that might inform metal enrichment within and beyond small galaxies in the early Universe. We find that simple assumptions for modeling the extent of supernova-driven winds agree with the simulation on average whereas inhomogeneous mixing and gas flows have a large effect on the spread in simulated stellar metallicities. In the context of the UFDs, this work demonstrates that simulations can form halos with a complex SFH and a large spread in the metallicity distribution function within a few hundred Myr in the early Universe. In particular, bursty and continuous star formation are seen in the simulation and both scenarios have been argued from the data. Spreads in the simulated metallicities, however remain too narrow and too metal-rich when compared to the UFDs. Future work is needed to help reduce these discrepancies and advance our interpretation of the data.

We produce synthetic spectra and observations for metal-free stellar populations and high mass X-ray binaries in the Renaissance Simulations at a redshift of 15. We extend our methodology from the first paper in the series by modelling the production and extinction of emission lines throughout a dusty and metal-enriched interstellar and circum-galactic media extracted from the simulation, using a Monte Carlo calculation. To capture the impact of high-energy photons, we include all frequencies from hard X-ray to far infrared with enough frequency resolution to discern line emission and absorption profiles. The most common lines in our sample in order of their rate of occurrence are Ly$\alpha$, the C IV $\lambda\lambda1548,1551$ doublet, H-$\alpha$, and the Ca II $\lambda\lambda\lambda8498,8542,8662$ triplet. The best scenario for a direct observation of a metal-free stellar population is a merger between two Population III galaxies. In mergers between metal-enriched and metal-free stellar populations, some characteristics may be inferred indirectly. Single Population III galaxies are too dim to be observed photometrically at $z = 15$. Ly$\alpha$ emission is discernible by $JWST$ as an increase in $\rm{J_{200w} - J_{277w}}$ colour off the intrinsic stellar tracks. Observations of metal-free stars will be difficult, though not impossible, with the next generation of space telescopes.

Semi-numeric methods have made it possible to efficiently model the epoch of reionisation (EoR). While most implementations involve a reduction to a simple three-parameter model, we introduce a new mass-dependent ionising efficiency parameter that folds in physical parameters that are constrained by the latest numerical simulations. This new parameterization enables the effective modeling of a broad range of host halo masses containing ionising sources, extending from the smallest Population III host halos with $M \sim 10^6 M_\odot$, which are often ignored, to the rarest cosmic peaks with $M \sim 10^{12} M_\odot$ during EoR. We compare the resulting ionising histories with a typical three-parameter model and also compare with the latest constraints from the Planck mission. Our model results in a optical depth due to Thomson scattering, $\tau_{\mathrm{e}}$ = 0.057, that is consistent with Planck. The largest difference in our model is shown in the resulting bubble size distributions which peak at lower characteristic sizes and are broadened. We also consider the uncertainties of the various physical parameters and comparing the resulting ionising histories broadly disfavors a small contribution from galaxies. As the smallest haloes cease a meaningful contribution to the ionising photon budget after $z = 10$, implying they play a role in determining the start of EoR and little else.

Accretion of gas and interaction of matter and radiation are at the heart of many questions pertaining to black hole (BH) growth and coevolution of massive BHs and their host galaxies. To answer them it is critical to quantify how the ionizing radiation that emanates from the innermost regions of the BH accretion flow couples to the surrounding medium and how it regulates the BH fueling. In this work we use high resolution 3-dimensional (3D) radiation-hydrodynamic simulations with the code Enzo, equipped with adaptive ray tracing module Moray, to investigate radiation-regulated BH accretion of cold gas. Our simulations reproduce findings from an earlier generation of 1D/2D simulations: the accretion powered UV and X-ray radiation forms a highly ionized bubble, which leads to suppression of BH accretion rate characterized by quasi-periodic outbursts. A new feature revealed by the 3D simulations is the highly turbulent nature of the gas flow in vicinity of the ionization front. During quiescent periods between accretion outbursts, the ionized bubble shrinks in size and the gas density that precedes the ionization front increases. Consequently, the 3D simulations show oscillations in the accretion rate of only ~2-3 orders of magnitude, significantly smaller than 1D/2D models. We calculate the energy budget of the gas flow and find that turbulence is the main contributor to the kinetic energy of the gas but corresponds to less than 10% of its thermal energy and thus does not contribute significantly to the pressure support of the gas.

We present synthetic observations for the first generations of galaxies in the Universe and make predictions for future deep field observations for redshifts greater than 6. Due to the strong impact of nebular emission lines and the relatively compact scale of HII regions, high resolution cosmological simulations and a robust suite of analysis tools are required to properly simulate spectra. We created a software pipeline consisting of FSPS, Hyperion, Cloudy and our own tools to generate synthetic IR observations from a fully three-dimensional arrangement of gas, dust, and stars. Our prescription allows us to include emission lines for a complete chemical network and tackle the effect of dust extinction and scattering in the various lines of sight. We provide spectra, 2-D binned photon imagery for both HST and JWST IR filters, luminosity relationships, and emission line strengths for a large sample of high redshift galaxies in the Renaissance Simulations. Our resulting synthetic spectra show high variability between galactic halos with a strong dependence on stellar mass, metallicity, gas mass fraction, and formation history. Halos with the lowest stellar mass have the greatest variability in [OIII]/H$\beta$, [OIII] and CIII] while halos with higher masses are seen to show consistency in their spectra and [OIII] equivalent widths (EW) between 1\AA\ and 10\AA. Viewing angle accounted for three-fold difference in flux due to the presence of ionized gas channels in a halo. Furthermore, JWST color plots show a discernible relationship between redshift, color, and mean stellar age.

One viable seeding mechanism for supermassive black holes is the direct gaseous collapse route in pre-galactic dark matter halos, producing objects on the order of $10^4 - 10^6$ solar masses. These events occur when the gas is prevented from cooling below $10^4$ K that requires a metal-free and relatively H$_2$-free medium. The initial collapse cools through atomic hydrogen transitions, but the gas becomes optically thick to the cooling radiation at high densities. We explore the effects ofLyman-$\alpha$ trapping in such a collapsing system with a suite of Monte Carlo radiation transport calculations in uniform density and isotropic cases that are based from a cosmological simulation. Our method includes both non-coherent scattering and two-photon line cooling. We find that Lyman-$\alpha$ radiation is marginally trapped in the parsec-scale gravitationally unstable central cloud, allowing the temperature to increase to 50,000 K at a number density of $3 \times 10^4$ cm$^{-3}$ and increasing the Jeans mass by a factor of five. The effective equation of state changes from isothermal at low densities to have an adiabatic index of 4/3 around the temperature maximum and then slowly retreats back to isothermal at higher densities. Our results suggest that Lyman-$\alpha$ trapping delays the initial collapse by raising the Jeans mass. Afterward the high density core cools back to $10^4$ K that is surrounded by a warm envelope whose inward pressure may alter the fragmentation scales at high densities.

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}$.

By heating the intergalactic medium (IGM) before reionization, X-rays are expected to play a prominent role in the early Universe. The cosmic 21-cm signal from this "Epoch of Heating" (EoH) could serve as a clean probe of high-energy processes inside the first galaxies. Here we improve on prior estimates of this signal by using high-resolution hydrodynamic simulations to calculate the X-ray absorption due to the interstellar medium (ISM) of the host galaxy. X-rays absorbed inside the host galaxy are unable to escape into the IGM and contribute to the EoH. We find that the X-ray opacity through these galaxies can be approximated by a metal-free ISM with a typical column density of log[N / cm^-2] = 21.4 +0.40-0.65. We compute the resulting 21-cm signal by combining these ISM opacities with public spectra of high-mass X-ray binaries (thought to be important X-ray sources in the early Universe). Our results support "standard scenarios" in which the X-ray heating of the IGM is inhomogeneous, and occurs before the bulk of reionization. The large-scale (k ~ 0.1/Mpc) 21-cm power reaches a peak of ~100 mK^2 at z = 10 - 15, with the redshift depending on the cosmic star formation history. This is in contrast to some recent work, motivated by the much larger X-ray absorption towards local HMXBs inside the Milky Way. Our main results can be reproduced by approximating the X-ray emission from HMXBs with a power-law spectrum with energy index alpha = 1, truncated at energies below 0.5 keV.

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.

Cosmic reionization is thought to be primarily fueled by the first generations of galaxies. We examine their stellar and gaseous properties, focusing on the star formation rates and the escape of ionizing photons, as a function of halo mass, redshift, and environment using the full suite of the {\it Renaissance Simulations} with an eye to provide better inputs to global reionization simulations. This suite, carried out with the adaptive mesh refinement code Enzo, is unprecedented in terms of their size and physical ingredients. The simulations probe overdense, average, and underdense regions of the universe of several hundred comoving Mpc$^3$, each yielding a sample of over 3,000 halos in the mass range $10^7 - 10^{9.5}~\Ms$ at their final redshifts of 15, 12.5, and 8, respectively. In the process, we simulate the effects of radiative and supernova feedback from 5,000 to 10,000 metal-free (Population III) stars in each simulation. We find that halos as small as $10^7~\Ms$ are able to form stars due to metal-line cooling from earlier enrichment by massive Population III stars. However, we find such halos do not form stars continuously. Using our large sample, we find that the galaxy-halo occupation fraction drops from unity at virial masses above $10^{8.5}~\Ms$ to $\sim$50\% at $10^8 ~\Ms$ and $\sim$10\% at $10^7~\Ms$, quite independent of redshift and region. Their average ionizing escape fraction is $\sim$5\% in the mass range $10^8 - 10^9~\Ms$ and increases with decreasing halo mass below this range, reaching 40--60\% at $10^7~\Ms$. Interestingly, we find that the escape fraction varies between 10--20\% in halos with virial masses $\sim 3 \times 10^9~\Ms$. Taken together, our results confirm the importance of the smallest galaxies as sources of ionizing radiation contributing to the reionization of the universe.

Due to their long mean free paths, X-rays are expected to have many significant impacts globally on the properties of the intergalactic medium (IGM) by their heating and ionizing processes on large scales. At high redshifts, X-rays from Population (Pop) III binaries might have important effects on cosmic reionization and the Lyman alpha forest. As a continuation of our previous work on Pop III binary X-rays (Xu et al. 2014), we use the Pop III distribution and evolution from the Renaissance Simulations, a suite of self-consistent cosmological radiation hydrodynamics simulations of the formation of the first galaxies, to calculate the X-ray luminosity density and background over the redshift range 20 > z > 7.6. As we find that Pop III star formation continues at a low, nearly constant rate to the end of reionization, X-rays are being continuously produced at significant rates compared to other possible X-ray sources, such as AGNs and normal X-ray binaries during the same period of time. We estimate that Pop III binaries produce approximately 6 eV of energy in the X-rays per hydrogen atom. We calculate the X-ray background for different monochromatic photon energies. KeV X-rays redshift and accumulate to produce a strong X-ray background spectrum extending to roughly 500 eV. The X-ray background is strong enough to heat the IGM to ~ 1000 K and to ionize a few percent of the neutral hydrogen. These effects are important for an understanding of the neutral hydrogen hyperfine transition 21-cm line signatures, the Ly alpha forest, and optical depth of the CMB to Thomson scattering.

Magnetic fields permeate the Universe on all scales and play a key role during star formation. We study the evolution of magnetic fields around a massive metal-free (Population III) star at $z \sim 15$ during the growth of its HII region and subsequent supernova explosion by conducting three cosmological magnetohydrodynamic simulations with radiation transport. Given the theoretical uncertainty and weak observational constraints of magnetic fields in the early universe, we initialize the simulations with identical initial conditions only varying the seed field strength. We find that magnetic fields grow as $\rho^{2/3}$ during the gravitational collapse preceding star formation, as expected from ideal spherical collapse models. Massive Population III stars can expel a majority of the gas from the host halo through radiative feedback, and we find that the magnetic fields are not amplified above the spherical collapse scaling relation during this phase. However, afterwards when its supernova remnant can radiatively cool and fragment, the turbulent velocity field in and around the shell causes the magnetic field to be significantly amplified on average by $\sim$100 in the shell and up to 6 orders of magnitude behind the reverse shock. Within the shell, field strengths are on the order of a few nG at a number density of 1 cm$^{-3}$. We show that this growth is primarily caused by small-scale dynamo action in the remnant. These strengthened fields will propagate into the first generations of galaxies, possibly affecting the nature of their star formation.

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 results on the formation of Pop III stars at redshift 7.6 from the Renaissance Simulations, a suite of extremely high-resolution and physics-rich radiation transport hydrodynamics cosmological adaptive-mesh refinement simulations of high redshift galaxy formation performed on the Blue Waters supercomputer. In a survey volume of about 220 comoving Mpc$^3$, we found 14 Pop III galaxies with recent star formation. The surprisingly late formation of Pop III stars is possible due to two factors: (i) the metal enrichment process is local and slow, leaving plenty of pristine gas to exist in the vast volume; and (ii) strong Lyman-Werner radiation from vigorous metal-enriched star formation in early galaxies suppresses Pop III formation in ("not so") small primordial halos with mass less than $\sim$ 3 $\times$ 10$^7$ M$_\odot$. We quantify the properties of these Pop III galaxies and their Pop III star formation environments. We look for analogues to the recently discovered luminous Ly $\alpha$ emitter CR7 (Sobral et al. 2015), which has been interpreted as a Pop III star cluster within or near a metal-enriched star forming galaxy. We find and discuss a system similar to this in some respects, however the Pop III star cluster is far less massive and luminous than CR7 is inferred to be.

We examine radiation-regulated accretion onto intermediate-mass and massive black holes (BHs) embedded in a bulge component. Using spherically symmetric one-dimensional radiation-hydrodynamics simulations, we track the growth of BHs accreting from a cold, neutral gas reservoir with temperature T=10^4 K. We find that the accretion rate of BHs embedded in bulges is proportional to r_{B,eff}/r_B, where r_{B,eff} is the increased effective Bondi radius that includes the gravitational potential of the bulge, and r_B is the Bondi radius of the BH. The radiative feedback from the BH suppresses the cold accretion rate to ~1 percent of the Bondi rate when a bulge is not considered. However, we find that the BH fueling rate increases rapidly when the bulge mass M_bulge is greater than the critical value of 10^6 M_sun and is proportional to M_bulge. Since the critical bulge mass is independent of the central BH mass M_{BH}, the growth rate of BHs with masses of 10^2, 10^4, and 10^6 M_sun exhibits distinct dependencies on the bulge-to-BH mass ratio. Our results imply that light seed BHs (<= 10^2 M_sun) which might be the remnants of the Pop III stars, cannot grow through accretion coevally with the early assembly of the bulge of the host galaxies until the bulge reaches the critical mass. However, massive BH seeds (>= 10^5 M_sun) that may form via direct collapse, are more likely to be embedded in a supercritical bulge and thus can grow efficiently coupling to the host galaxies and driving the early evolution of the M_{BH}-$\sigma$ relationship.

We present a simulation of the formation of the earliest Population II stars, starting from cosmological initial conditions and ending when metals created in the first supernovae are incorporated into a collapsing gas-cloud. This occurs after a supernova blast-wave collides with a nearby mini-halo, inducing further turbulence that efficiently mixes metals into the dense gas in the center of the halo. The gas that first collapses has been enriched to a metallicity of Z ~ 2e-5 Zsun. Due to the extremely low metallicity, collapse proceeds similarly to metal-free gas until dust cooling becomes efficient at high densities, causing the cloud to fragment into a large number of low mass objects. This external enrichment mechanism provides a plausible origin for the most metal-poor stars observed, such as SMSS J031300.36-670839.3, that appear to have formed out of gas enriched by a single supernova. This mechanism operates on shorter timescales than the time for low-mass mini-halos (M < 5e5 Msun) to recover their gas after experiencing a supernova. As such, metal-enriched stars will likely form first via this channel if the conditions are right for it to occur. We identify a number of other externally enriched halos that may form stars in this manner. These halos have metallicities as high as 0.01 Zsun, suggesting that some members of the first generation of metal-enriched stars may be hiding in plain sight in current stellar surveys.

In this paper, we present the first results from the Renaissance Simulations, a suite of extremely high-resolution and physics-rich AMR calculations of high redshift galaxy formation performed on the Blue Waters supercomputer. These simulations contain hundreds of well-resolved galaxies at $z \sim 25-8$, and make several novel, testable predictions. Most critically, we show that the ultraviolet luminosity function of our simulated galaxies is consistent with observations of high-z galaxy populations at the bright end of the luminosity function (M$_{1600} \leq -17$), but at lower luminosities is essentially flat rather than rising steeply, as has been inferred by Schechter function fits to high-z observations, and has a clearly-defined lower limit in UV luminosity. This behavior of the luminosity function is due to two factors: (i) the strong dependence of the star formation rate on halo virial mass in our simulated galaxy population, with lower-mass halos having systematically lower star formation rates and thus lower UV luminosities; and (ii) the fact that halos with virial masses below $\simeq 2 \times 10^8$ M$_\odot$ do not universally contain stars, with the fraction of halos containing stars dropping to zero at $\simeq 7 \times 10^6$ M$_\odot$. Finally, we show that the brightest of our simulated galaxies may be visible to current and future ultra-deep space-based surveys, particularly if lensed regions are chosen for observation.

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.

We describe Rabacus, a Python package for calculating the transfer of hydrogen ionizing radiation in simplified geometries relevant to astronomy and cosmology. We present example solutions for three specific cases: 1) a semi-infinite slab gas distribution in a homogeneous isotropic background, 2) a spherically symmetric gas distribution with a point source at the center, and 3) a spherically symmetric gas distribution in a homogeneous isotropic background. All problems can accommodate arbitrary spectra and density profiles as input. The solutions include a treatment of both hydrogen and helium, a self-consistent calculation of equilibrium temperatures, and the transfer of recombination radiation. The core routines are written in Fortran 90 and then wrapped in Python leading to execution speeds thousands of times faster than equivalent routines written in pure Python. In addition, all variables have associated units for ease of analysis. The software is part of the Python Package Index and the source code is available on Bitbucket at https://bitbucket.org/galtay/rabacus . In addition, installation instructions and a detailed users guide are available at http://pythonhosted.org//rabacus .

We present our prediction for the local 21 cm differential brightness temperature ($\delta T_{b}$) from a set of strongly clustered sources of Population III (Pop III) and II (Pop II) objects in the early Universe, by a numerical simulation of their formation and radiative feedback. These objects are located inside a highly biased environment, which is a rare, high-density peak ("Rarepeak") extending to $\sim7$ comoving Mpc. We study the impact of ultraviolet (UV) and X-ray photons on the intergalactic medium (IGM) and the resulting $\delta T_{b}$, when Pop III stars are assumed to emit X-ray photons by forming X-ray binaries very efficiently. We parameterize the rest-frame spectral energy distribution (SED) of X-ray photons, which regulates X-ray photon-trapping, IGM-heating, secondary Lyman-alpha pumping and the resulting morphology of $\delta T_{b}$. A combination of emission ($\delta T_{b}>0$) and absorption ($\delta T_{b}<0$) regions appears in varying amplitudes and angular scales. The boost of the signal by the high-density environment ($\delta\sim0.64$) and on a relatively large scale combine to make Rarepeak a discernible, spatially-extended ($\theta\sim10'$) object for 21 cm observation at $13\lesssim z\lesssim17$, which is found to be detectable as a single object by SKA with integration time of $\sim1000$ hours. Power spectrum analysis by some of the SKA precursors (LOFAR, MWA, PAPER) of such rare peaks is found difficult due to the rarity of these peaks, and the contribution only by these rare peaks to the total power spectrum remains subdominant compared to that by all astrophysical sources.

We describe an extension of the Enzo code to enable fully-coupled radiation hydrodynamical simulation of inhomogeneous reionization in large $\sim (100 Mpc)^3$ cosmological volumes with thousands to millions of point sources. We solve all dynamical, radiative transfer, thermal, and ionization processes self-consistently on the same mesh, as opposed to a postprocessing approach which coarse-grains the radiative transfer. We do, however, employ a simple subgrid model for star formation which we calibrate to observations. Radiation transport is done in the grey flux-limited diffusion (FLD) approximation, which is solved by implicit time integration split off from the gas energy and ionization equations, which are solved separately. This results in a faster and more robust scheme for cosmological applications compared to the earlier method. The FLD equation is solved using the hypre optimally scalable geometric multigrid solver from LLNL. By treating the ionizing radiation as a grid field as opposed to rays, our method is scalable with respect to the number of ionizing sources, limited only by the parallel scaling properties of the radiation solver. We test the speed and accuracy of our approach on a number of standard verification and validation tests. We show by direct comparison with Enzo's adaptive ray tracing method Moray that the well-known inability of FLD to cast a shadow behind opaque clouds has a minor effect on the evolution of ionized volume and mass fractions in a reionization simulation validation test. We illustrate an application of our method to the problem of inhomogeneous reionization in a 80 Mpc comoving box resolved with $3200^3$ Eulerian grid cells and dark matter particles.

In the last decade, the growth of supermassive black holes (SMBHs) has been intricately linked to galaxy formation and evolution and is a key ingredient in the assembly of galaxies. To investigate the origin of SMBHs, we perform cosmological simulations that target the direct collapse black hole (DCBH) seed formation scenario in the presence of two different strong Lyman-Werner (LW) background fields. These simulations include the X-ray irradiation from a central massive black hole (MBH), $\rm{H}_2$ self-shielding and stellar feedback from metal-free and metal-enriched stars. We find in both simulations that local X-ray feedback induces metal-free star formation $\sim 0.5$ Myr after the MBH forms. The MBH accretion rate reaches a maximum of $10^{-3}$ $M_{\odot}$ yr$^{-1}$ in both simulations. However, the duty cycle differs which is derived to be $6\%$ and $50\%$ for high and low LW cases, respectively. The MBH in the high LW case grows only $\sim 6\%$ in 100 Myr compared to $16\%$ in the low LW case. We find that the maximum accretion rate is determined by the local gas thermodynamics whereas the duty cycle is determined by the large scale gas dynamics and gas reservoir. We conclude that radiative feedback from the central MBH plays an important role in star formation in the nuclear regions and stifling initial MBH growth, relative to the typical Eddington rate argument, and that initial MBH growth might be affected by the local LW radiation field.

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.

The first galaxies in the Universe are the building blocks of all observed galaxies. We present scaling relations for galaxies forming at redshifts $z \ge 15$ when reionization is just beginning. We utilize the ``Rarepeak' cosmological radiation hydrodynamics simulation that captures the complete star formation history in over 3,300 galaxies, starting with massive Population III stars that form in dark matter halos as small as ~$10^6 M_\odot$. We make various correlations between the bulk halo quantities, such as virial, gas, and stellar masses and metallicities and their respective accretion rates, quantifying a variety of properties of the first galaxies up to halo masses of $10^9 M_\odot$. Galaxy formation is not solely relegated to atomic cooling halos with virial temperatures greater than $10^4$ K, where we find a dichotomy in galaxy properties between halos above and below this critical mass scale. Halos below the atomic cooling limit have a stellar mass -- halo mass relationship $\log M_\star \simeq 3.5 + 1.3\log(M_{\rm vir} / 10^7 M_\odot)$. We find a non-monotonic relationship between metallicity and halo mass for the smallest galaxies. Their initial star formation events enrich the interstellar medium and subsequent star formation to a median of $10^{-2} Z_\odot$ and $10^{-1.5} Z_\odot$, respectively, in halos of total mass $10^7 M_\odot$ that is then diluted by metal-poor inflows, well beyond Population III pre-enrichment levels of $10^{-3.5} Z_\odot$. The scaling relations presented here can be employed in models of reionization, galaxy formation and chemical evolution in order to consider these galaxies forming prior to reionization.

Due to their long mean free path, X-rays are expected to have an important impact on cosmic reionization by heating and ionizing the IGM on large scales, especially after simulations have suggested that Population III stars may form in pairs at redshifts as high as 20 - 30. We use the Pop III distribution and evolution from a self-consistent cosmological radiation hydrodynamics simulation of the formation of the first galaxies and a simple Pop III X-ray binaries model to estimate their X-ray output in a high density region larger than 100 comoving (Mpc)$^3$. We then combine three different methods --- ray tracing, a one-zone model, and X-ray background modeling --- to investigate the X-ray propagation, intensity distribution, and long term effects on the IGM thermal and ionization state. The efficiency and morphology of photo-heating and photo-ionization are dependent on the photon energies. The sub-keV X-rays only impact the IGM near the sources while the keV photons contribute significantly to the X-ray background and heat and ionize the IGM smoothly. The X-rays just below 1 keV are most efficient in heating and ionizing the IGM. We find that the IGM might be heated to over 100 K by $z=10$ and the high density source region might reach 10$^4$ K, limited by atomic hydrogen cooling. This may be important for predicting the 21-cm neutral hydrogen signals. But, on the other hand, the free electrons from X-ray ionizations are not enough to contribute significantly to the optical depth of CMB to the Thomson scattering.

Starlight from galaxies plays a pivotal role throughout the process of cosmic reionisation. We present the statistics of dwarf galaxy properties at z > 7 in haloes with masses up to 10^9 solar masses, using a cosmological radiation hydrodynamics simulation that follows their buildup starting with their Population III progenitors. We find that metal-enriched star formation is not restricted to atomic cooling ($T_{\rm vir} \ge 10^4$ K) haloes, but can occur in haloes down to masses ~10^6 solar masses, especially in neutral regions. Even though these smallest galaxies only host up to 10^4 solar masses of stars, they provide nearly 30 per cent of the ionising photon budget. We find that the galaxy luminosity function flattens above M_UV ~ -12 with a number density that is unchanged at z < 10. The fraction of ionising radiation escaping into the intergalactic medium is inversely dependent on halo mass, decreasing from 50 to 5 per cent in the mass range $\log M/M_\odot = 7.0-8.5$. Using our galaxy statistics in a semi-analytic reionisation model, we find a Thomson scattering optical depth consistent with the latest Planck results, while still being consistent with the UV emissivity constraints provided by Ly$\alpha$ forest observations at z = 4-6.

We introduce the AGORA project, a comprehensive numerical study of well-resolved galaxies within the LCDM cosmology. Cosmological hydrodynamic simulations with force resolutions of ~100 proper pc or better will be run with a variety of code platforms to follow the hierarchical growth, star formation history, morphological transformation, and the cycle of baryons in and out of 8 galaxies with halo masses M_vir ~= 1e10, 1e11, 1e12, and 1e13 Msun at z=0 and two different ("violent" and "quiescent") assembly histories. The numerical techniques and implementations used in this project include the smoothed particle hydrodynamics codes GADGET and GASOLINE, and the adaptive mesh refinement codes ART, ENZO, and RAMSES. The codes will share common initial conditions and common astrophysics packages including UV background, metal-dependent radiative cooling, metal and energy yields of supernovae, and stellar initial mass function. These are described in detail in the present paper. Subgrid star formation and feedback prescriptions will be tuned to provide a realistic interstellar and circumgalactic medium using a non-cosmological disk galaxy simulation. Cosmological runs will be systematically compared with each other using a common analysis toolkit, and validated against observations to verify that the solutions are robust - i.e., that the astrophysical assumptions are responsible for any success, rather than artifacts of particular implementations. The goals of the AGORA project are, broadly speaking, to raise the realism and predictive power of galaxy simulations and the understanding of the feedback processes that regulate galaxy "metabolism." The proof-of-concept dark matter-only test of the formation of a galactic halo with a z=0 mass of M_vir ~= 1.7e11 Msun by 9 different versions of the participating codes is also presented to validate the infrastructure of the project.

The electron scattering optical depth constraints on reionization suggest that there may be other sources that contribute to the ionization of hydrogen aside from observable star forming galaxies. Often the calculated value of the electron scattering optical depth, \tau_{es}, falls below the measurements derived from observations of the CMB or an assumption about non-observable sources must be made in order to reach agreement. Here, we calculate the hydrogen ionization fraction as a function of redshift and the electron scattering optical depth from both galaxies and active galactic nuclei (AGN) factoring in the secondary collisional ionizations from the AGN X-ray emission. In this paper we use the most current determination of the evolving hard X-ray luminosity function and extrapolate its evolution beyond z = 6. The AGN spectral energy distributions (SEDs) include both UV and X-ray ionizing photons. To search for the largest possible effect, all AGN are assumed to have \lambda_{Edd} = 1.0 and be completely unobscured. The results show that AGNs produce a perturbative effect on the reionization of hydrogen and remains in agreement with current constraints. Our calculations find the epoch of reionization still ends at z \approx 6 and only increases the electron scattering optical depth by ~2.3% under the most optimal conditions. This can only be moderately increased by assuming a constant black hole mass of M_{BH} = 10^(5) solar masses. As a result, we conclude that there is a need for other sources beyond observable galaxies and AGNs that contribute to the reionization of hydrogen at z > 6.

This paper describes the open-source code Enzo, which uses block-structured adaptive mesh refinement to provide high spatial and temporal resolution for modeling astrophysical fluid flows. The code is Cartesian, can be run in 1, 2, and 3 dimensions, and supports a wide variety of physics including hydrodynamics, ideal and non-ideal magnetohydrodynamics, N-body dynamics (and, more broadly, self-gravity of fluids and particles), primordial gas chemistry, optically-thin radiative cooling of primordial and metal-enriched plasmas (as well as some optically-thick cooling models), radiation transport, cosmological expansion, and models for star formation and feedback in a cosmological context. In addition to explaining the algorithms implemented, we present solutions for a wide range of test problems, demonstrate the code's parallel performance, and discuss the Enzo collaboration's code development methodology.

Recent simulations of Population III star formation have suggested that some fraction form in binary systems, in addition to having a characteristic mass of tens of solar masses. The deaths of metal-free stars result in the initial chemical enrichment of the universe and the production of the first stellar-mass black holes. Here we present a cosmological adaptive mesh refinement simulation of an overdense region that forms a few 10^9 Msun dark matter halos and over 13,000 Population III stars by redshift 15. We find that most halos do not form Population III stars until they reach Mvir ~ 10^7 Msun because this biased region is quickly enriched from both Population III and galaxies, which also produce high levels of ultraviolet radiation that suppress H2 formation. Nevertheless, Population III stars continue to form, albeit in more massive halos, at a rate of ~ 10^{-4} Msun yr^{-1} Mpc^{-3} at redshift 15. The most massive starless halo has a mass of 7 x 10^7 Msun, which could host massive black hole formation through the direct gaseous collapse scenario. We show that the multiplicity of the Population III remnants grows with halo mass above 10^8 Msun, culminating in 50 remnants located in 10^9 Msun halos on average. This has implications that high mass X-ray binaries and intermediate mass black holes that originate from metal-free stars may be abundant in high-redshift galaxies.

Observational studies show that there is a strong link between the formation and evolution of galaxies and the growth of supermassive black holes (SMBH) at their centers. However, the underlying physics of this observed relation is poorly understood. In order to study the effects of X-ray radiation on the surroundings of the black hole, we implement X-ray Dominated Region (XDR) physics into Enzo and use the radiation transport module Moray to calculate the radiative transfer for a polychromatic spectrum. In this work, we investigate the effects of X-ray irradiation, produced by a central massive black hole (MBH) with a mass of M = 5x10^4 M_(solar), on ambient gas with solar and zero metallicity. We find that in the solar metallicity case, due to high opacity of the metals, the energy deposition rate in the central region (< 20 pc) is high and hence the temperatures in the center are on the order of 10^(5-7) K. Moreover, due to the cooling ability and high intrinsic opacity of solar metallicity gas, column densities of 10^(24) cm^(-2) are reached at a radius of 20 pc from the MBH. These column densities are about 3 orders of magnitudes higher than in the zero metallicity case. Furthermore, in the zero metallicity case an X-ray induced H II region is formed already after 5.8 Myr. This causes a significant outflow of gas (~8x10^6 M_(solar) from the central region, with the gas reaching outflow velocities up to ~100 km s^(-1). At later times, ~23 Myr after we insert the MBH, we find that the solar metallicity case also develops an X-ray induced H II region, but delayed by ~17 Myr.

DECal is a new calibration system for the CTIO 4 m Blanco telescope. It is currently being installed as part of the Dark Energy Survey and will provide both broadband flat fields and narrowband (about 1 nm bandwidth) spectrophotometric calibration for the new Dark Energy Camera (DECam). Both of these systems share a new Lambertian flat field screen. The broadband flat field system uses LEDs to illuminate each photometric filter. The spectrophotometric calibration system consists of a monochromator-based tunable light source that is projected onto the flat field screen using a custom line-to-spot fiber bundle and an engineered diffuser. Several calibrated photodiodes positioned along the beam monitor the telescope throughput as a function of wavelength. This system will measure the wavelength-dependent instrumental response function of the total telescope+instrument system in the range 300 < lambda < 1100nm. The spectrophotometric calibration will be performed regularly (roughly once per month) to determine the spectral response of the DECam system and to monitor changes in instrumental throughput during the five year Dark Energy Survey project.

Traditional color and airmass corrections can typically achieve ~0.02 mag precision in photometric observing conditions. A major limiting factor is the variability in atmospheric throughput, which changes on timescales of less than a night. We present preliminary results for a system to monitor the throughput of the atmosphere, which should enable photometric precision when coupled to more traditional techniques of less than 1% in photometric conditions. The system, aTmCam, consists of a set of imagers each with a narrow-band filter that monitors the brightness of suitable standard stars. Each narrowband filter is selected to monitor a different wavelength region of the atmospheric transmission, including regions dominated by the precipitable water absorption and aerosol scattering. We have built a prototype system to test the notion that an atmospheric model derived from a few color indices measurements can be an accurate representation of the true atmospheric transmission. We have measured the atmospheric transmission with both narrowband photometric measurements and spec- troscopic measurements; we show that the narrowband imaging approach can predict the changes in the throughput of the atmosphere to better than ~10% across a broad wavelength range, so as to achieve photometric precision less than 0.01 mag.

We explore structure formation in the dark ages ($z\sim 30-6$) using two well-known methods for initializing cosmological $N$-body simulations. Overall, both the Zel'dovich approximation (\za) and second order Lagrangian perturbation theory (\lpt) are known to produce accurate present-day dark matter halo mass functions. However, since the \lpt method drives more rapid evolution of dense regions, it increases the occurrence of rare massive objects -- an effect that is most pronounced at high redshift. We find that \lpt produces more halos that could harbor Population III stars and their black hole remnants, and they produce them earlier. Although the differences between the \lpt and \za mass functions are nearly erased by $z=6$, this small boost to the number and mass of black holes more than doubles the reionized volume of the early Universe. We discuss the implications for reionization and massive black hole growth.

Massive stars provide feedback that shapes the interstellar medium of galaxies at all redshifts and their resulting stellar populations. Here we present three adaptive mesh refinement radiation hydrodynamics simulations that illustrate the impact of momentum transfer from ionising radiation to the absorbing gas on star formation in high-redshift dwarf galaxies. Momentum transfer is calculated by solving the radiative transfer equation with a ray tracing algorithm that is adaptive in spatial and angular coordinates. We find that momentum input partially affects star formation by increasing the turbulent support to a three-dimensional rms velocity equal to the circular velocity of early haloes. Compared to a calculation that neglects radiation pressure, the star formation rate is decreased by a factor of five to 1.8 x 10^{-2} Msun/yr in a dwarf galaxy with a dark matter and stellar mass of 2.0 x 10^8 and 4.5 x 10^5 solar masses, respectively, when radiation pressure is included. Its mean metallicity of 10^{-2.1} Z_sun is consistent with the observed dwarf galaxy luminosity-metallicity relation. However, what one may naively expect from the calculation without radiation pressure, the central region of the galaxy overcools and produces a compact, metal-rich stellar population with an average metallicity of 0.3 Z_sun, indicative of an incorrect physical recipe. In addition to photo-heating in HII regions, radiation pressure further drives dense gas from star forming regions, so supernovae feedback occurs in a warmer and more diffuse medium, launching metal-rich outflows. Capturing this aspect and a temporal separation between the start of radiative and supernova feedback are numerically important in the modeling of galaxies to avoid the "overcooling problem". We estimate that dust in early low-mass galaxies is unlikely to aid in momentum transfer from radiation to the gas.

The recent implementation of radiative transfer algorithms in numerous hydrodynamics codes has led to a dramatic improvement in studies of feedback in various astrophysical environments. However, because of methodological limitations and computational expense, the spectra of radiation sources are generally sampled at only a few evenly-spaced discrete emission frequencies. Using one-dimensional radiative transfer calculations, we investigate the discrepancies in gas properties surrounding model stars and accreting black holes that arise solely due to spectral discretization. We find that even in the idealized case of a static and uniform density field, commonly used discretization schemes induce errors in the neutral fraction and temperature by factors of two to three on average, and by over an order of magnitude in certain column density regimes. The consequences are most severe for radiative feedback operating on large scales, dense clumps of gas, and media consisting of multiple chemical species. We have developed a method for optimally constructing discrete spectra, and show that for two test cases of interest, carefully chosen four-bin spectra can eliminate errors associated with frequency resolution to high precision. Applying these findings to a fully three-dimensional radiation-hydrodynamic simulation of the early universe, we find that the HII region around a primordial star is substantially altered in both size and morphology, corroborating the one-dimensional prediction that discrete spectral energy distributions can lead to sizable inaccuracies in the physical properties of a medium, and as a result, the subsequent evolution and observable signatures of objects embedded within it.

The first stars in the universe are thought to be massive, forming in dark matter halos with masses around 10^6 solar masses. Recent simulations suggest that these metal-free (Population III) stars may form in binary or multiple systems. Because of their high stellar masses and small host halos, their feedback ionizes the surrounding 3 kpc of intergalactic medium and drives the majority of the gas from the potential well. The next generation of stars then must form in this gas-poor environment, creating the first galaxies that produce the majority of ionizing radiation during cosmic reionization. I will review the latest developments in the field of Population III star formation and feedback and its impact on galaxy formation prior to reionization. In particular, I will focus on the numerical simulations that have demonstrated this sequence of events, ultimately leading to cosmic reionization.

We describe cosmological galaxy formation simulations with the adaptive mesh refinement code Enzo that incorporate a star formation prescription regulated by the local abundance of molecular hydrogen. We show that this H2-regulated prescription leads to a suppression of star formation in low mass halos (M_h < ~10^10 M_sun) at z>4, alleviating some of the dwarf galaxy problems faced by theoretical galaxy formation models. H2 regulation modifies the efficiency of star formation of cold gas directly, rather than indirectly reducing the cold gas content with "supernova feedback". We determine the local H2 abundance in our most refined grid cells (76 proper parsec in size at z=4) by applying the model of Krumholz, McKee, & Tumlinson, which is based on idealized 1D radiative transfer calculations of H2 formation-dissociation balance in ~100 pc atomic--molecular complexes. Our H2-regulated simulations are able to reproduce the empirical (albeit lower z) Kennicutt-Schmidt relation, including the low Sigma_gas cutoff due to the transition from atomic to molecular phase and the metallicity dependence thereof, without the use of an explicit density threshold in our star formation prescription. We compare the evolution of the luminosity function, stellar mass density, and star formation rate density from our simulations to recent observational determinations of the same at z=4-8 and find reasonable agreement between the two.

By definition, Population III stars are metal-free, and their protostellar collapse is driven by molecular hydrogen cooling in the gas-phase, leading to large characteristic masses. Population II stars with lower characteristic masses form when the star-forming gas reaches a critical metallicity of 10^{-6} - 10^{-3.5} Z_\odot. We present an adaptive mesh refinement radiation hydrodynamics simulation that follows the transition from Population III to II star formation. The maximum spatial resolution of 1 comoving parsec allows for individual molecular clouds to be well-resolved and their stellar associations to be studied in detail. We model stellar radiative feedback with adaptive ray tracing. A top-heavy initial mass function for the Population III stars is considered, resulting in a plausible distribution of pair-instability supernovae and associated metal enrichment. We find that the gas fraction recovers from 5 percent to nearly the cosmic fraction in halos with merger histories rich in halos above 10^7 solar masses. A single pair-instability supernova is sufficient to enrich the host halo to a metallicity floor of 10^{-3} Z_\odot and to transition to Population II star formation. This provides a natural explanation for the observed floor on damped Lyman alpha (DLA) systems metallicities reported in the literature, which is of this order. We find that stellar metallicities do not necessarily trace stellar ages, as mergers of halos with established stellar populations can create superpositions of t-Z evolutionary tracks. A bimodal metallicity distribution is created after a starburst occurs when the halo can cool efficiently through atomic line cooling.

We describe a photon-conserving radiative transfer algorithm, using a spatially-adaptive ray tracing scheme, and its parallel implementation into the adaptive mesh refinement (AMR) cosmological hydrodynamics code, Enzo. By coupling the solver with the energy equation and non-equilibrium chemistry network, our radiation hydrodynamics framework can be utilised to study a broad range of astrophysical problems, such as stellar and black hole (BH) feedback. Inaccuracies can arise from large timesteps and poor sampling, therefore we devised an adaptive time-stepping scheme and a fast approximation of the optically-thin radiation field with multiple sources. We test the method with several radiative transfer and radiation hydrodynamics tests that are given in Iliev et al. (2006, 2009). We further test our method with more dynamical situations, for example, the propagation of an ionisation front through a Rayleigh-Taylor instability, time-varying luminosities, and collimated radiation. The test suite also includes an expanding H II region in a magnetised medium, utilising the newly implemented magnetohydrodynamics module in Enzo. This method linearly scales with the number of point sources and number of grid cells. Our implementation is scalable to 512 processors on distributed memory machines and can include radiation pressure and secondary ionisations from X-ray radiation. It is included in the newest public release of Enzo.

Coalescing massive black hole binaries are produced by the mergers of galaxies. The final stages of the black hole coalescence produce strong gravitational radiation that can be detected by the space-borne LISA. In cases where the black hole merger takes place in the presence of gas and magnetic fields, various types of electromagnetic signals may also be produced. Modeling such electromagnetic counterparts of the final merger requires evolving the behavior of both gas and fields in the strong-field regions around the black holes. We have taken a step towards solving this problem by mapping the flow of pressureless matter in the dynamic, 3-D general relativistic spacetime around the merging black holes. We find qualitative differences in collision and outflow speeds, including a signature of the merger when the net angular momentum of the matter is low, between the results from single and binary black holes, and between nonrotating and rotating holes in binaries. If future magnetohydrodynamic results confirm these differences, it may allow assessment of the properties of the binaries as well as yielding an identifiable electromagnetic counterpart to the attendant gravitational wave signal.

It has been argued that low-luminosity dwarf galaxies are the dominant source of ionizing radiation during cosmological reionization. The fraction of ionizing radiation that escapes into the intergalactic medium from dwarf galaxies with masses less than ~10^9.5 solar masses plays a critical role during this epoch. Using an extensive suite of very high resolution (0.1 pc), adaptive mesh refinement, radiation hydrodynamical simulations of idealized and cosmological dwarf galaxies, we characterize the behavior of the escape fraction in galaxies between 3 x 10^6 and 3 x 10^9 solar masses with different spin parameters, amounts of turbulence, and baryon mass fractions. For a given halo mass, escape fractions can vary up to a factor of two, depending on the initial setup of the idealized halo. In a cosmological setting, we find that the time-averaged photon escape fraction always exceeds 25% and reaches up to 80% in halos with masses above 10^8 solar masses with a top-heavy IMF. The instantaneous escape fraction can vary up to an order of magnitude in a few million years and tend to be positively correlated with star formation rate. We find that the mean of the star formation efficiency times ionizing photon escape fraction, averaged over all atomic cooling (T_vir > 8000 K) galaxies, ranges from 0.02 for a normal IMF to 0.03 for a top-heavy IMF, whereas smaller, molecular cooling galaxies in minihalos do not make a significant contribution to reionizing the universe due to a much lower star formation efficiency. These results provide the physical basis for cosmological reionization by stellar sources, predominately atomic cooling dwarf galaxies.

The first stars in the universe, forming at redshifts z>15 in minihalos with masses of order 10^6 Msun, may leave behind black holes as their remnants. These objects could conceivably serve as "seeds" for much larger black holes observed at redshifts z~6. We study the growth of the remnant black holes through accretion including for the first time the emitted accretion radiation with adaptive mesh refinement cosmological radiation-hydrodynamical simulations. The effects of photo-ionization and heating dramatically affect the accretion flow from large scales, resulting in negligible mass growth of the black hole. We compare cases with the accretion luminosity included and neglected to show that the accretion radiation drastically changes the environment within 100 pc of the black hole, where gas temperatures are increased by an order of magnitude. The gas densities are reduced and further star formation in the same minihalo prevented for the two hundred million years of evolution we followed. These calculations show that even without the radiative feedback included most seed black holes do not gain mass as efficiently as has been hoped for in previous theories, implying that black hole remnants of Pop III stars that formed in minihalos are not likely to be the origin of miniquasars. Most importantly, however, these calculations demonstrate that if early stellar mass black holes are indeed accreting close to the Bondi-Hoyle rate with ten percent efficiency they have a dramatic local effect in regulating star formation in the first galaxies.

(Abridged) Using ab initio cosmological Eulerian adaptive mesh refinement radiation hydrodynamical calculations, we discuss how very massive stars start the process of cosmological reionization. The models include non-equilibrium primordial gas chemistry and cooling processes and accurate radiation transport in the Case B approximation using adaptively ray traced photon packages, retaining the time derivative in the transport equation. Supernova feedback is modeled by thermal explosions triggered at parsec scales. All calculations resolve the local Jeans length by at least 16 grid cells at all times and as such cover a spatial dynamic range of ~10^6. These first sources of reionization are highly intermittent and anisotropic and first photoionize the small scales voids surrounding the halos they form in, rather than the dense filaments they are embedded in. As the merging objects form larger, dwarf sized galaxies, the escape fraction of UV radiation decreases and the HII regions only break out on some sides of the galaxies making them even more anisotropic. In three cases, SN blast waves induce star formation in overdense regions that were formed earlier from ionization front instabilities. These stars form tens of parsecs away from the center of their parent DM halo. Approximately 5 ionizing photons are needed per sustained ionization when star formation in 10^6 M_sun halos are dominant in the calculation. As the halos become larger than ~10^7 M_sun, the ionizing photon escape fraction decreases, which in turn increases the number of photons per ionization to 15-50, in calculations with stellar feedback only. Supernova feedback in these more massive halos creates a more diffuse medium, allowing the stellar radiation to escape more easily and maintaining the ratio of 5 ionizing photons per sustained ionization.

The first stars form in dark matter halos of masses ~10^6 M_sun as suggested by an increasing number of numerical simulations. Radiation feedback from these stars expels most of the gas from their shallow potential well of their surrounding dark matter halos. We use cosmological adaptive mesh refinement simulations that include self-consistent Population III star formation and feedback to examine the properties of assembling early dwarf galaxies. Accurate radiative transport is modeled with adaptive ray tracing. We include supernova explosions and follow the metal enrichment of the intergalactic medium. The calculations focus on the formation of several dwarf galaxies and their progenitors. In these halos, baryon fractions in 10^8 solar mass halos decrease by a factor of 2 with stellar feedback and by a factor of 3 with supernova explosions. We find that radiation feedback and supernova explosions increase gaseous spin parameters up to a factor of 4 and vary with time. Stellar feedback, supernova explosions, and H_2 cooling create a complex, multi-phase interstellar medium whose densities and temperatures can span up to 6 orders of magnitude at a given radius. The pair-instability supernovae of Population III stars alone enrich the halos with virial temperatures of 10^4 K to approximately 10^{-3} of solar metallicity. We find that 40% of the heavy elements resides in the intergalactic medium (IGM) at the end of our calculations. The highest metallicity gas exists in supernova remnants and very dilute regions of the IGM.

Cosmic structure originated from minute density perturbations in an almost homogeneous universe. The first stars are believed to be very massive and luminous, providing the first ionizing radiation and heavy elements to the universe and forming 100 million years after the Big Bang. The impact from primordial stellar radiation is far reaching and affects subsequent star and galaxy formation. In this thesis, we present results from adaptive mesh refinement calculations of the formation of the first galaxies. We gradually introduce important physical processes, such as molecular hydrogen cooling and stellar feedback, to base models that only consider atomic hydrogen and helium cooling. In these base models, we find that gas in dark matter halos with masses ~10^8 solar masses centrally collapse before multiple fragmentation occurs in a global disc. We then investigate the importance of molecular hydrogen cooling in early structure formation in the presence of a soft ultraviolet radiation background. We find that molecular hydrogen plays an important role in star formation in halos well below a virial temperature of 10,000 K even in the most extreme assumptions of negative radiative feedback. We also present results from the first radiation hydrodynamics calculations of early dwarf galaxy formation. We develop a novel technique, adaptive ray tracing, to accurately transport radiation from primordial stars. We find primordial stellar feedback alters the landscape of early galaxy formation in that its angular momentum is increased and baryon fractions are decreased. We also describe the metal enrichment of the intergalactic medium and early dwarf galaxies. Finally we explore cosmological reionization by these massive, metal-free stars and its effects on star formation in early galaxies.

Numerous cosmological hydrodynamic studies have addressed the formation of galaxies. Here we choose to study the first stages of galaxy formation, including non-equilibrium atomic primordial gas cooling, gravity and hydrodynamics. Using initial conditions appropriate for the concordance cosmological model of structure formation, we perform two adaptive mesh refinement simulations of ~10^8 M_sun galaxies at high redshift. The calculations resolve the Jeans length at all times with more than 16 cells and capture over 14 orders of magnitude in length scales. In both cases, the dense, 10^5 solar mass, one parsec central regions are found to contract rapidly and have turbulent Mach numbers up to 4. Despite the ever decreasing Jeans length of the isothermal gas, we only find one site of fragmentation during the collapse. However, rotational secular bar instabilities transport angular momentum outwards in the central parsec as the gas continues to collapse and lead to multiple nested unstable fragments with decreasing masses down to sub-Jupiter mass scales. Although these numerical experiments neglect star formation and feedback, they clearly highlight the physics of turbulence in gravitationally collapsing gas. The angular momentum segregation seen in our calculations plays an important role in theories that form supermassive black holes from gaseous collapse.

The first luminous objects in the concordance cosmology form by molecular hydrogen cooling in dark matter dominated halos of masses ~10^6 M_sun. We use Eulerian adaptive mesh refinement simulations to demonstrate that in the presence of a large soft ultraviolet radiation background, molecular hydrogen is the dominant coolant. Even for very large radiation backgrounds, the halo masses that cool and collapse are up to two orders of magnitude smaller than the halos that cool via atomic hydrogen line cooling. The abundance of cooling halos and the cosmic mass fraction contained within them depends exponentially on this critical mass scale. Consequently, the majority of current models of cosmological reionization, chemical evolution, supermassive black hole formation, and galaxy formation underestimate the number of star forming progenitors of a given system by orders of magnitude. At the highest redshifts, this disagreement is largest. We also show that even in the absence of residual electrons, collisional ionization in central shocks create a sufficient amount of electrons to form molecular hydrogen and cool the gas in halos of virial temperatures far below the atomic cooling limit.

(Abridged) Galaxies form in hierarchically assembling dark matter halos. With cosmological three dimensional adaptive mesh refinement simulations, we explore in detail the virialization of baryons in the concordance cosmology, including optically thin primordial gas cooling. We focus on early protogalaxies with virial temperatures of 10^4 K and their progenitors. Without cooling, virial heating occurs in shocks close to the virial radius for material falling in from voids. Material in dense filaments penetrates deeper to about half that radius. With cooling the virial shock position shrinks and also the filaments reach scales as small as a third the virial radius. The temperatures in protogalaxies found in adiabatic simulations decrease by a factor of two from the center and show flat entropy cores. In cooling halos the gas reaches virial equilibrium with the dark matter potential through its turbulent velocities. We observe turbulent Mach numbers ranging from one to three in the cooling cases. This turbulence is driven by the large scale merging and interestingly remains supersonic in the centers of these early galaxies even in the absence of any feedback processes. The virial theorem is shown to approximately hold over 3 orders of magnitude in length scale with the turbulent pressure prevailing over the thermal energy. The turbulent velocity distributions are Maxwellian and by far dominate the small rotation velocities associated with the total angular momentum of the galaxies. Decomposing the velocity field using the Cauchy-Stokes theorem, we show that ample amounts of vorticity are present around shocks even at the very centers of these objects.

The concordance model of cosmology and structure formation predicts the formation of isolated very massive stars at high redshifts in dark matter dominated halos of 10^5 to 10^6 Msun. These stars photo-ionize their host primordial molecular clouds, expelling all the baryons from their halos. When the stars die, a relic HII region is formed within which large amounts of molecular hydrogen form which will allow the gas to cool efficiently when gravity assembles it into larger dark matter halos. The filaments surrounding the first star hosting halo are largely shielded and provide the pathway for gas to stream into the halo when the star has died. We present the first fully three dimensional cosmological radiation hydrodynamical simulations that follow all these effects. A novel adaptive ray casting technique incorporates the time dependent radiative transfer around point sources. This approach is fast enough so that radiation transport, kinetic rate equations, and hydrodynamics are solved self-consistently. It retains the time derivative of the transfer equation and is explicitly photon conserving. This method is integrated with the cosmological adaptive mesh refinement code enzo, and runs on distributed and shared memory parallel architectures. Where applicable the three dimensional calculation not only confirm expectations from earlier one dimensional results but also illustrate the multi-fold hydrodynamic complexities of HII regions. In the absence of stellar winds the circumstellar environments of the first supernovae and putative early gamma-ray bursts will be of low density ~1 cm^-3. Albeit marginally resolved, ionization front instabilities lead to cometary and elephant trunk like small scale structures reminiscent of nearby star forming regions.

Recent simulations of the formation of the first luminous objects in the universe predict isolated very massive stars to form in dark matter halos with virial temperatures large enough to allow significant amounts of molecular hydrogen to form. We construct a semi-analytic model based on the Press-Schechter formalism and calibrate the minimum halos mass that may form a primordial star with the results from extensive adaptive mesh refinement simulations. The model also includes star formation in objects with virial temperatures in excess of ten thousand Kelvin. The free parameters are tuned to match the optical depth measurements by the WMAP satellite. The models explicitly includes the negative feedback of the destruction of molecular hydrogen by a soft UV background which is computed self-consistently. We predict high redshift supernova rates as one of the most promising tools to test the current scenario of primordial star formation. The supernova rate from primordial stars peaks at redshifts ~20. Using an analytic model for the luminosities of pair-instability supernovae we predict observable magnitudes and discuss possible observational strategies. Such supernovae would release enough metals corresponding to a uniform enrichment to a few hundred thousands of solar metalicity. If some of these stars produce gamma ray bursts our rates will be directly applicable to understanding the anticipated results from the SWIFT satellite. This study highlights the great potential for the James Webb space telescope in probing cosmic structure at redshifts greater than 20.

We present the results of a search for variability in the equivalent widths (EWs) of narrow, associated (Delta v < 5,000 km/s) absorption lines found in the UV spectra of z < 1.5 quasars. The goal of this search was to use variability as a means of identifying absorption lines arising in gas that is intrinsic to the quasar central engine. We have compared archival HST/FOS spectra of quasars with recent spectra obtained as part of our own snapshot survey of the same objects with STIS. The intervals between observations are 4-10 years. We primarily focused on the C IV absorption lines, although we also studied other lines available in the same spectra (e.g., Ly-alpha, N V, O VI). Our main result is that 4 out of 15 quasars, or 4 out of 19 associated absorption systems, contained variable narrow absorption lines, which are indicative of intrinsic absorption. We conclude that a minimum of 21% of the associated absorption-line systems are variable. Because not all systems will have necessarily varied, this is a lower limit on this fraction and it is consistent with previous estimates based on variability, partial coverage analysis, or statistical arguments. If we interpret the upper limits on the variability time scale as upper limits on the recombination time of the absorber, we constrain the density of the absorber to be n_e > 3000 cm^{-3} and its distance from the ionizing source to be R < 100pc. Moreover, we are now able to pick out specific intrinsic absorption-line systems to be followed up with high-dispersion spectroscopy in order to constrain the geometry, location, and physical conditions of the absorber. We briefly illustrate how followup studies can yield such constraints by means of a simulation.

We present the first orbital elements for the massive close binary, HD 115071, a double-lined spectroscopic binary in a circular orbit with a period of 2.73135 +/- 0.00003 days. The orbital semiamplitudes indicate a mass ratio of M_2/M_1 = 0.58 +/- 0.02 and yet the stars have similar luminosities. We used a Doppler tomography algorithm to reconstruct the individual component optical spectra, and we applied well known criteria to arrive at classifications of O9.5 V and B0.2 III for the primary and secondary, respectively. We present models of the Hipparcos light curve of the ellipsoidal variations caused by the tidal distortion of the secondary, and the best fit model for a Roche-filling secondary occurs for an inclination of i = 48.7 +/- 2.1 degrees. The resulting masses are 11.6 +/- 1.1 and 6.7 +/- 0.7 solar masses for the primary and secondary, respectively, so that both stars are very overluminous for their mass. The system is one of only a few known semi-detached, Algol-type binaries that contain O-stars. We suggest that the binary has recently emerged from extensive mass transfer (possibly through a delayed contact and common envelope process).