Gamma-ray burst (GRB) data suggest that the jets from GRBs in the high redshift universe are more narrowly collimated than those at lower redshifts. This implies that we detect relatively fewer long GRB progenitor systems (i.e. massive stars) at high redshifts, because a greater fraction of GRBs have their jets pointed away from us. As a result, estimates of the star formation rate (from the GRB rate) at high redshifts may be diminished if this effect is not taken into account. In this paper, we estimate the star formation rate (SFR) using the observed GRB rate, accounting for an evolving jet opening angle. We find that the SFR in the early universe (z > 3) can be up to an order of magnitude higher than the canonical estimates, depending on the severity of beaming angle evolution and the fraction of stars that make long gamma-ray bursts. Additionally, we find an excess in the SFR at low redshifts, although this lessens when accounting for evolution of the beaming angle. Finally, under the assumption that GRBs do in fact trace canonical forms of the cosmic SFR, we constrain the resulting fraction of stars that must produce GRBs, again accounting for jet beaming-angle evolution. We find this assumption suggests a high fraction of stars in the early universe producing GRBs - a result that may, in fact, support our initial assertion that GRBs do not trace canonical estimates of the SFR.

Jet opening angles of long gamma-ray bursts (lGRBs) appear to evolve in cosmic time, with lGRBs at higher redshifts being on average more narrowly beamed than those at lower redshifts. We examine the nature of this anti-correlation in the context of collimation by the progenitor stellar envelope. First, we show that the data indicate a strong correlation between gamma-ray luminosity and jet opening angle, and suggest this is a natural selection effect - only the most luminous GRBs are able to successfully launch jets with large opening angles. Then, by considering progenitor properties expected to evolve through cosmic time, we show that denser stars lead to more collimated jets; we argue that the apparent anti-correlation between opening angle and redshift can be accounted for if lGRB massive star progenitors at high redshifts have higher average density compared to those at lower redshifts. This may be viable for an evolving IMF - under the assumption that average density scales directly with mass, this relationship is consistent with the form of the IMF mass evolution suggested in the literature. The jet angle-redshift anti-correlation may also be explained if the lGRB progenitor population is dominated by massive stars at high redshift, while lower redshift lGRBs allow for a greater diversity of progenitor systems (that may fail to collimate the jet as acutely). Overall, however, we find both the jet angle-redshift anti-correlation and jet angle-luminosity correlation are consistent with the conditions of jet launch through, and collimation by, the envelope of a massive star progenitor.

The classic Bondi solution remains a common starting point both for studying black hole growth across cosmic time in cosmological simulations and for smaller scale simulations of AGN feedback. In nature, however, there will be inhomogenous distributions of rotational velocity and density along the outer radius ($R_o$) marking the sphere of influence of a black hole. While there have been many studies of how the Bondi solution changes with a prescribed angular momentum boundary condition, they have all assumed a constant density at $R_o$. In this Letter, we show that a non-uniform density at $R_o$ causes a meridional flow and due to conservation of angular momentum, the Bondi solution qualitatively changes into an inflow-outflow solution. Using physical arguments, we analytically identify the critical logarithmic density gradient $|\partial{\ln{\rho}}/\partial{\theta}|$ above which this change of the solution occurs. For realistic $R_o$, this critical gradient is less than 0.01 and tends to 0 as $R_o \rightarrow \infty$. We show using numerical simulations that, unlike for solutions with an imposed rotational velocity, the accretion rate for solutions under an inhomogenous density boundary condition remains constant at nearly the Bondi rate $\dot{M}_B$, while the outflow rate can greatly exceed $\dot{M}_B$.

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.

We examine the relationship between a number of long gamma-ray burst (lGRB) properties (isotropic emitted energy, luminosity, intrinsic duration, jet opening angle) and redshift. We find that even when accounting for conservative detector flux limits, there appears to be a significant correlation between isotropic equivalent energy and redshift, suggesting cosmological evolution of the lGRB progenitor. Analyzing a sub-sample of lGRBs with jet opening angle estimates, we find the beaming-corrected lGRB emitted energy does not correlate with redshift, but jet opening angle does. Additionally, we find a statistically significant anti-correlation between the intrinsic prompt duration and redshift, even when accounting for potential selection effects. We also find that - for a given redshift - isotropic energy is positively correlated with intrinsic prompt duration. None of these GRB properties appear to be correlated with galactic offset. From our selection-effect-corrected redshift distribution, we estimate a co-moving rate density for lGRBs, and compare this to the global cosmic star formation rate (SFR). We find the lGRB rate mildly exceeds the global star formation rate between a redshift of 3 and 5, and declines rapidly at redshifts above this (although we cannot constrain the lGRB rate above a redshift of about 6 due to sample incompleteness). We find the lGRB rate diverges significantly from the SFR at lower redshifts. We discuss both the correlations and lGRB rate density in terms of various lGRB progenitor models and their apparent preference for low-metallicity environments.

High redshift quasars emit copious X-ray photons which heat the intergalactic medium to temperatures up to $\sim$ 10$^6$ K. At such high temperatures the primordial gas will not form stars until it is assembled into dark matter haloes with masses of up to $\sim$ 10$^{11}$ M$_{\odot}$, at which point the hot gas collapses and cools under the influence of gravity. Once this occurs, there is a massive reservoir of primordial gas from which stars can form, potentially setting the stage for the brightest Population (Pop) III starbursts in the early Universe. Supporting this scenario, recent observations of quasars at z $\sim$ 6 have revealed a lack of accompanying Lyman $\alpha$ emitting galaxies, consistent with suppression of primordial star formation in haloes with masses below $\sim$ 10$^{10}$ M$_{\odot}$. Here we model the chemical and thermal evolution of the primordial gas as it collapses into such a massive halo irradiated by a nearby quasar in the run-up to a massive Pop III starburst. We find that within $\sim$ 100 kpc of the highest redshift quasars discovered to date the Lyman-Werner flux produced in the quasar host galaxy may be high enough to stimulate the formation of a direct collapse black hole (DCBH). A survey with single pointings of the NIRCam instrument at individual known high-z quasars may be a promising strategy for finding Pop III stars and DCBHs with the James Webb Space Telescope.

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.

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.

Intermediate mass protostarsprovide a bridge between theories of low- and high-mass star formation. Emerging molecular outflows can be used to determine the influence of fragmentation and multiplicity on protostellar evolution through the correlation of outflow forces of intermediate mass protostars with the luminosity. The aim of this paper is to derive outflow forces from outflows of six intermediate mass protostellar regions and validate the apparent correlation between total luminosity and outflow force seen in earlier work, as well as remove uncertainties caused by different methodology. By comparing CO 6--5 observations obtained with APEX with non-LTE radiative transfer model predictions, optical depths, temperatures, densities of the gas of the molecular outflows are derived. Outflow forces, dynamical timescales and kinetic luminosities are subsequently calculated. Outflow parameters, including the forces, were derived for all sources. Temperatures in excess of 50 K were found for all flows, in line with recent low-mass results. However, comparison with other studies could not corroborate conclusions from earlier work on intermediate mass protostars which hypothesized that fragmentation enhances outflow forces in clustered intermediate mass star formation. Any enhancement in comparison with the classical relation between outflow force and luminosity can be attributed the use of a higher excitation line and improvement in methods; They are in line with results from low-mass protostars using similar techniques. The role of fragmentation on outflows is an important ingredient to understand clustered star formation and the link between low and high-mass star formation. However, detailed information on spatial scales of a few 100 AU, covering all individual members is needed to make the necessary progress.

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 large cross-section of the Lyman alpha (Lya) line makes it a sensitive probe of the ionization state of the intergalactic medium (IGM). Here we present the most complete study to date of the IGM Lya opacity, and its application to the redshift evolution of the 'Lya fraction', i.e. the fraction of color-selected galaxies with a detectable Lya emission line. We use a tiered approach, which combines large-scale semi-numeric simulations of reionization with moderate-scale hydrodynamic simulations of the ionized IGM. This allows us to simultaneously account for evolution in both: (i) the opacity from an incomplete (patchy) reionization, parameterized by the filling factor of ionized regions, Q_HII; and (ii) the opacity from self-shielded systems in the ionized IGM, parameterized by the average photo-ionization rate inside HII regions, \Gamma. In contrast to recent empirical models, attenuation from patchy reionization has a unimodal distribution along different sightlines, while attenuation from self-shielded systems is more bimodal. We quantify the average IGM transmission in our (Q_HII, \Gamma) parameter space, which can easily be used to interpret new data sets. Using current observations, we predict that the Lya fraction cannot drop by more than a factor of ~2 with IGM attenuation alone, even for HII filling factors as low as Q_HII>0.1. Larger changes in the Lya fraction could result from a co-evolution with galaxy properties. Marginalizing over \Gamma, we find that current observations constrain Q_HII < 0.6 at z=7 [68% confidence level (C.L.)]. However, all of our parameter space is consistent with observations at 95% C.L., highlighting the need for larger observational samples at z >= 6.

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.

In 1997, an extension of general relativity was proposed that predicts the dark energy density \Lambda\ to vary linearly with the total number of macroscopic black holes in the universe. We explore this prediction and find that \Lambda\ must be roughly constant after the bulk of the stellar mass black holes are in place, so for a redshift z smaller than unity. Conversely, the highest black hole formation rate corresponds to the peak in the cosmic star formation history, earlier than z=1. This yields a fast declining \Lambda, by a factor of about 5, from redshift 1 to 3. At even earlier times, before many stars were formed, the value of \Lambda\ should be much smaller than its current value. These predicted effects are all consistent with current data, and near future observations can definitively confirm or disproof the link between the dark energy density and the total number of black holes in the universe.

The initial mass function (IMF) of the first (Population III) stars and Population II (Pop II) stars is poorly known due to a lack of observations of the period between recombination and reionization. In simulations of the formation of the first stars, it has been shown that, due to the limited ability of metal-free primordial gas to cool, the IMF of the first stars is a few orders of magnitude more massive than the current IMF. The transition from a high-mass IMF of the first stars to a lower-mass current IMF is thus important to understand. To study the underlying physics of this transition, we performed several simulations using the cosmological hydrodynamic adaptive mesh refinement code Enzo for metallicities of 10^{-4}, 10^{-3}, 10^{-2}, and 10^{-1} Z_{\odot}. In our simulations we include a star formation prescription that is derived from a metallicity dependent multi-phase ISM structure, an external UV radiation field, and a mechanical feedback algorithm. We also implement cosmic ray heating, photoelectric heating and gas-dust heating/cooling, and follow the metal enrichment of the ISM. It is found that the interplay between metallicity and UV radiation leads to the co-existence of Pop III and Pop II star formation in non-zero metallicity (Z/Z_{\odot} \geq10^{-2}) gas. A cold (T<100 K) and dense (\rho>10^{-22} g cm^{-3}) gas phase is fragile to ambient UV radiation. In a metal-poor (Z/Z_{\odot} \leq10^{-3}) gas, the cold and dense gas phase does not form in the presence of a radiation field of F_{0}\sim10^{-5}-10^{-4} erg cm^{-2} s^{-1}. Therefore, metallicity by itself is not a good indicator of the Pop III-Pop II transition. Metal-rich (Z/Z_{\odot}\geq10^{-2}) gas dynamically evolves two to three orders of magnitude faster than metal poor gas (Z/Z_{\odot}\leq10^{-3}). The simulations including SNe show that pre-enrichment of the halo does not affect the mixing of metals.

Star formation is regulated through a variety of feedback processes. In this study, we treat feedback by metal injection and a UV background as well as by X-ray irradiation. Our aim is to investigate whether star formation is significantly affected when the ISM of a proto-galaxxy enjoys different metallicities and when a star forming cloud resides in the vicinity of a strong X-ray source. We perform cosmological Enzo simulations with a detailed treatment of non-zero metallicity chemistry and thermal balance. We also perform FLASH simulations with embedded Lagrangian sink particles of a collapsing molecular cloud near a massive, 10^{7} M\odot, black hole that produces X-ray radiation. We find that a multi-phase ISM forms for metallicites as small as 10^{-4} Solar at z = 6, with higher (10^{-2}Z\odot) metallicities supporting a cold (< 100 K) and dense (> 10^{3} cm^{-3}) phase at higher (z = 20) redshift. A star formation recipe based on the presence of a cold dense phase leads to a self-regulating mode in the presence of supernova and radiation feedback. We also find that when there is strong X-ray feedback a collapsing cloud fragments into larger clumps whereby fewer but more massive protostellar cores are formed. This is a consequence of the higher Jeans mass in the warm (50 K, due to ionization heating) molecular gas. Accretion processes dominate the mass function and a near-flat, non-Salpeter IMF results.