Using cosmological simulations to make useful, scientifically relevant emission line predictions is a relatively new and rapidly evolving field. However, nebular emission lines have been particularly challenging to model because they are extremely sensitive to the local photoionization balance, which can be driven by a spatially dispersed distribution of stars amidst an inhomogeneous absorbing medium of dust and gas. As such, several unmodeled mysteries in observed emission line patterns exist in the literature. For example, there is some question as to why He II $\lambda 4686$/H$\beta$ ratios in observations of lower-metallicity dwarf galaxies tend to be higher than model predictions. Since hydrodynamic cosmological simulations are best suited to this mass and metallicity regime, this question presents a good test case for the development of a robust emission line modeling pipeline. The pipeline described in this work can model a process that produces high He II $\lambda 4686$/H$\beta$ ratios and eliminate some of the modeling discrepancy for ratios below 3% without including AGNs, X-ray binaries, high mass binaries, or a top-heavy stellar initial mass function. These ratios are found to be more sensitive to the presence of 15 Myr or longer gaps in the star formation histories than to extraordinary ionization parameters or specific star formation rates. This work also charts a path forward for the next generation of nebular emission line modeling studies.

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.

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.

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.