The unique combination of superb spatial resolution, wide-area coverage, and deep depth of the optical imaging from the Hyper Suprime-Cam (HSC) Subaru Strategic Program is utilized to search for dual quasar candidates. Using an automated image analysis routine on 34,476 known SDSS quasars, we identify those with two (or more) distinct optical point sources in HSC images covering 796 deg^2. We find 421 candidates out to a redshift of 4.5 of which one hundred or so are more likely after filtering out contaminating stars. Angular separations of 0.6 - 4.0" correspond to projected separations of 3 - 30 kpc, a range relatively unexplored for population studies of luminous dual quasars. Using Keck-I/LRIS and Gemini-N/NIFS, we spectroscopically confirm three dual quasar systems at z < 1, two of which are previously unknown out of eight observed, based on the presence of characteristic broad emission lines in each component, while highlighting that the continuum of one object in one of the pairs is reddened. In all cases, the [OIII]5007 emission lines have mild velocity offsets, thus the joint [OIII] line profile is not double-peaked. We find a dual quasar fraction of 0.26+/-0.18% and no evidence for evolution. A comparison with the Horizon-AGN simulation seems to support the case of no evolution in the dual quasar fraction when broadly matching the quasar selection. These results may indicate a scenario in which the frequency of the simultaneous triggering of luminous quasars is not as sensitive as expected to the cosmic evolution of the merger rate or gas content of galaxies.

The leading contenders for the seeds of the first quasars are direct collapse black holes (DCBHs) formed during catastrophic baryon collapse in atomically-cooled halos at $z \sim$ 20. The discovery of the Ly$\alpha$ emitter CR7 at $z =$ 6.6 was initially held to be the first detection of a DCBH, although this interpretation has since been challenged on the grounds of Spitzer IRAC and Very Large Telescope X-Shooter data. Here we determine if radio flux from a DCBH in CR7 could be detected and discriminated from competing sources of radio emission in the halo such as young supernovae and H II regions. We find that a DCBH would emit a flux of 10 - 200 nJy at 1.0 GHz, far greater than the sub-nJy signal expected for young supernovae but on par with continuum emission from star-forming regions. However, radio emission from a DCBH in CR7 could be distinguished from free-free emission from H II regions by its spectral evolution with frequency and could be detected by the Square Kilometer Array in the coming decade.

The first supernovae enrich the previously pristine gas with metals, out of which the next generation of stars form. Based on hydrodynamical simulations, we develop a new stochastic model to predict the metallicity of star-forming gas in the first galaxies. On average, in internally enriched galaxies, the metals are well mixed with the pristine gas. However, in externally enriched galaxies, the metals can not easily penetrate into the dense gas, which yields a significant metallicity difference between the star-forming and average gas inside a halo. To study the consequences of this effect, we apply a semi-analytical model to Milky Way-like dark matter merger trees and follow stellar fossils from high redshift until the present day with a novel realistic metal mixing recipe. We calibrate the model to reproduce the metallicity distribution function (MDF) at low metallicities and find that a primordial IMF with slope of $\mathrm{d}N/\mathrm{d}M \propto M^{-0.5}$ from $2 Msun$ to $180 Msun$ best reproduces the MDF. Our improved model for inhomogeneous mixing can have a large impact for individual minihalos, but does not significantly influence the modelled MDF at [Fe/H]$\gtrsim -4$ or the best-fitting Pop~III IMF.

Supermassive primordial stars in hot, atomically-cooling haloes at $z \sim$ 15 - 20 may have given birth to the first quasars in the universe. Most simulations of these rapidly accreting stars suggest that they are red, cool hypergiants, but more recent models indicate that some may have been bluer and hotter, with surface temperatures of 20,000 - 40,000 K. These stars have spectral features that are quite distinct from those of cooler stars and may have different detection limits in the near infrared (NIR) today. Here, we present spectra and AB magnitudes for hot, blue supermassive primordial stars calculated with the TLUSTY and CLOUDY codes. We find that photometric detections of these stars by the James Webb Space Telescope (JWST) will be limited to $z \lesssim$ 10 - 12, lower redshifts than those at which red stars can be found, because of quenching by their accretion envelopes. With moderate gravitational lensing, Euclid and the Wide-Field Infrared Space Telescope (WFIRST) could detect blue supermassive stars out to similar redshifts in wide-field surveys.

In recent years, the discovery of massive quasars at z~7 has provided a striking challenge to our understanding of the origin and growth of supermassive black holes in the early Universe. Mounting observational and theoretical evidence indicates the viability of massive seeds, formed by the collapse of supermassive stars, as a progenitor model for such early, massive accreting black holes. Although considerable progress has been made in our theoretical understanding, many questions remain regarding how (and how often) such objects may form, how they live and die, and how next generation observatories may yield new insight into the origin of these primordial titans. This review focusses on our present understanding of this remarkable formation scenario, based on discussions held at the Monash Prato Centre from November 20--24, 2017, during the workshop "Titans of the Early Universe: The Origin of the First Supermassive Black Holes."

We present a novel scenario for the formation of carbon-enhanced metal-poor (CEMP) stars. Carbon enhancement at low stellar metallicities is usually considered a consequence of faint or other exotic supernovae. An analytical estimate of cooling times in low-metallicity gas demonstrates a natural bias, which favours the formation of CEMP stars as a consequence of inhomogeneous metal mixing: carbon-rich gas has a shorter cooling time and can form stars prior to a potential nearby pocket of carbon-normal gas, in which star formation is then suppressed due to energetic photons from the carbon-enhanced protostars. We demonstrate that this scenario provides a natural formation mechanism for CEMP stars from carbon-normal supernovae, if inhomogeneous metal mixing provides carbonicity differences of at least one order of magnitude separated by >10pc. In our fiducial (optimistic) model, 8% (83%) of observed CEMP-no stars ([Ba/Fe]<0) can be explained by this formation channel. This new scenario may change our understanding of the first supernovae and thereby our concept of the first stars. Future 3D simulations are required to assess the likelihood of this mechanism to occur in typical high-redshift galaxies.

The collapse of supermassive primordial stars in hot, atomically-cooled halos may have given birth to the first quasars at $z \sim$ 15 - 20. Recent numerical simulations of these rapidly accreting stars reveal that they are cool, red hypergiants shrouded by dense envelopes of pristine atomically-cooled gas at 6,000 - 8,000 K, with luminosities $L$ $\gtrsim$ 10$^{10}$ L$_{\odot}$. Could such luminous but cool objects be detected as the first stage of quasar formation in future near infrared (NIR) surveys? We have now calculated the spectra of supermassive primordial stars in their birth envelopes with the Cloudy code. We find that some of these stars will be visible to JWST at $z \lesssim$ 20 and that with modest gravitational lensing Euclid and WFIRST could detect them out to $z \sim$ 10 - 12. Rather than obscuring the star, its accretion envelope enhances its visibility in the NIR today by reprocessing its short-wavelength flux into photons that are just redward of the Lyman limit in the rest frame of the star.

Extremely metal poor (EMP) stars in the Milky Way inherited the chemical composition of the gas out of which they formed. They therefore carry the chemical fingerprint of the first stars in their spectral lines. It is commonly assumed that EMP stars form from gas that was enriched by only one progenitor supernova ('mono-enriched'). However, recent numerical simulations show that the first stars form in small clusters. Consequently, we expect several supernovae to contribute to the abundances of an EMP star ('multi-enriched'). We analyse seven recently observed EMP stars from the TOPoS survey by applying the divergence of the chemical displacement and find that J1035+0641 is mono-enriched ($p_{mono}=53\%$) and J1507+0051 is multi-enriched ($p_{mono}=4\%$). For the remaining five stars we can not make a distinct prediction ($p_{mono} \lesssim 50\%$) due to theoretical and observational uncertainties. Further observations in the near-UV will help to improve our diagnostic and therefore contribute to constrain the nature of the first stars.

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

The recent observations of supermassive black holes (SMBHs) at high redshift challenge our understanding of their formation and growth. There are different proposed pathways to form black hole (BH) seeds, such as the remnants of the first stars (chapter 4), gas-dynamical processes (chapter 5), direct collapse (chapter 6), or stellar collisions in dense nuclear clusters (chapter 7). In this chapter, we discuss the probability of forming supermassive black holes (SMBHs) via these channels and determine the expected number density of the BH seeds. We start with a brief discussion of the observational constraints on SMBHs at low and high redshift that theoretical models have to reproduce (a more detailed account is provided in chapter 12). We further present the most popular formation channels of SMBHs, discuss under which conditions they can reproduce the observations, and compare various estimates in the literature on the expected number density of SMBHs. To account for the density of quasars at $z>6$ requires very efficient gas accretion mechanisms or high BH seeds masses. The bottleneck to obtain sufficiently high number densities of seed BHs with masses $>10^5$M$_\odot$ is the interplay between radiative and chemical feedback, which constrains the conditions for primordial, isothermal gas collapse.

Extremely metal-poor (EMP) stars in the Milky Way (MW) allow us to infer the properties of their progenitors by comparing their chemical composition to the metal yields of the first supernovae. This method is most powerful when applied to mono-enriched stars, i.e. stars that formed from gas that was enriched by only one previous supernova. We present a novel diagnostic to identify this subclass of EMP stars. We model the first generations of star formation semi-analytically, based on dark matter halo merger trees that yield MW-like halos at the present day. Radiative and chemical feedback are included self-consistently and we trace all elements up to zinc. Mono-enriched stars account for only $\sim 1\%$ of second generation stars in our fiducial model and we provide an analytical formula for this probability. We also present a novel analytical diagnostic to identify mono-enriched stars, based on the metal yields of the first supernovae. This new diagnostic allows us to derive our main results independently from the specific assumptions made regarding Pop III star formation, and we apply it to a set of observed EMP stars to demonstrate its strengths and limitations. Our results may provide selection criteria for current and future surveys and therefore contribute to a deeper understanding of EMP stars and their progenitors.

Pair-instability supernovae (PISNe) are very luminous explosions of massive, low metallicity stars. They can potentially be observed out to high redshifts due to their high explosion energies, thus providing a probe of the Universe prior to reionization. The near-infrared camera, NIRCam, on board the James Webb Space Telescope is ideally suited for detecting their redshifted ultraviolet emission. We calculate the photometric signature of high-redshift PISNe and derive the optimal detection strategy for identifying their prompt emission and possible afterglow. We differentiate between PISNe and other sources that could have a similar photometric signature, such as active galactic nuclei or high-redshift galaxies. We demonstrate that the optimal strategy, which maximizes the visibility time of the PISN lightcurve per invested exposure time, consists of the two wide-band filters F200W and F356W with an exposure time of 600s. For such exposures, we expect one PISN at $z \lesssim 7.5$ per at least 50,000 different field of view, which can be accomplished with parallel observations and an extensive archival search. The PISN afterglow, caused by nebular emission and reverberation, is very faint and requires unfeasibly long exposure times to be uniquely identified. However, this afterglow would be visible for several hundred years, about two orders of magnitude longer than the prompt emission, rendering PISNe promising targets for future, even more powerful telescopes.

Galactic outflows, driven by active galactic nuclei (AGN), play a crucial role in galaxy formation and in the self-regulated growth of supermassive black holes (BHs). AGN feedback couples to and affects gas, rather than stars, and in many, if not most, gas-rich galaxies cold gas is rotationally supported and settles in a disc. We present a 2D analytical model for AGN-driven outflows in a gaseous disc and demonstrate the main improvements, compared to existing 1D solutions. We find significant differences for the outflow dynamics and wind efficiency. The outflow is energy-driven due to inefficient cooling up to a certain AGN luminosity ($\sim 10^{43}$erg/s in our fiducial model), above which the outflow remains momentum-driven in the disc up to galactic scales. We reproduce results of 3D simulations that gas is preferentially ejected perpendicular to the disc and find that the fraction of ejected interstellar medium is lower than in 1D models. The recovery time of gas in the disc, defined as the freefall time from the radius to which the AGN pushes the ISM at most, is remarkably short, of the order 1Myr. This indicates that AGN-driven winds cannot suppress BH growth for long. Without the inclusion of supernova feedback, we find a scaling of the black hole mass with the halo velocity dispersion of $M_\mathrm{BH} \propto \sigma ^{4.8}$.

Dwarf galaxy anomalies, such as their abundance and cusp-core problems, remain a prime challenge in our understanding of galaxy formation. The inclusion of baryonic physics could potentially solve these issues, but the efficiency of stellar feedback is still controversial. We analytically explore the possibility of feedback from Active Galactic Nuclei (AGN) in dwarf galaxies and compare AGN and supernova (SN) feedback. We assume the presence of an intermediate mass black hole within low mass galaxies and standard scaling relations between the relevant physical quantities. We model the propagation and properties of the outflow and explore the critical condition for global gas ejection. Performing the same calculation for SNe, we compare the ability of AGN and SNe to drive gas out of galaxies. We find that a critical halo mass exists below which AGN feedback can remove gas from the host halo and that the critical halo mass for AGN is greater than the equivalent for SNe in a significant part of the parameter space, suggesting that AGN could provide an alternative and more successful source of negative feedback than SNe, even in the most massive dwarf galaxies.

The search for metal-free stars has so far been unsuccessful, proving that if there are surviving stars from the first generation, they are rare, they have been polluted, or we have been looking in the wrong place. To predict the likely location of Population~III (Pop~III) survivors, we semi-analytically model early star formation in progenitors of Milky Way-like galaxies and their environments. We base our model on merger trees from the high-resolution dark matter only simulation suite \textit{Caterpillar}. Radiative and chemical feedback are taken into account self-consistently, based on the spatial distribution of the haloes. Our results are consistent with the non-detection of Pop III survivors in the Milky Way today. We find that possible surviving Population III stars are more common in Milky Way satellites than in the main Galaxy. In particular, low mass Milky Way satellites contain a much larger fraction of Pop~III stars than the Milky Way. Such nearby, low mass Milky Way satellites are promising targets for future attempts to find Pop~III survivors, especially for high-resolution, high signal-to-noise spectroscopic observations. We provide the probabilities for finding a Pop~III survivor in the red giant branch phase for all known Milky Way satellites to guide future observations.

With new observational facilities becoming available soon, discovering and characterising supernovae from the first stars will open up alternative observational windows to the end of the cosmic dark ages. Based on a semi-analytical merger tree model of early star formation we constrain Population III supernova rates. We find that our method reproduces the Population III supernova rates of large-scale cosmological simulations very well. Our computationally efficient model allows us to survey a large parameter space and to explore a wide range of different scenarios for Population III star formation. Our calculations show that observations of the first supernovae can be used to differentiate between cold and warm dark matter models and to constrain the corresponding particle mass of the latter. Our predictions can also be used to optimize survey strategies with the goal to maximize supernova detection rates.

CR7 is the brightest Lyman-$\alpha$ emitter observed at $z>6$, which shows very strong Lyman-$\alpha$ and HeII 1640\AA\ line luminosities, but no metal line emission. Previous studies suggest that CR7 hosts either young primordial stars with a total stellar mass of $\sim 10^7\,\mathrm{M}_\odot$ or a black hole of $\gtrsim 10^6\,\mathrm{M}_\odot$. Here, we explore different formation scenarios for CR7 with a semianalytical model, based on the random sampling of dark matter merger trees. We are unable to reproduce the observational constraints with a primordial stellar source, given our model assumptions, due to the short stellar lifetimes and the early metal enrichment. Black holes that are the remnants of the first stars are either not massive enough, or reside in metal-polluted haloes, ruling out this possible explanation of CR7. Our models instead suggest that direct collapse black holes, which form in metal-free haloes exposed to large Lyman-Werner fluxes, are more likely the origin of CR7. However, this result is derived under optimistic assumptions and future observations are necessary to further constrain the nature of CR7.

Gravitational waves (GWs) provide a revolutionary tool to investigate yet unobserved astrophysical objects. Especially the first stars, which are believed to be more massive than present-day stars, might be indirectly observable via the merger of their compact remnants. We develop a self-consistent, cosmologically representative, semi-analytical model to simulate the formation of the first stars. By extrapolating binary stellar-evolution models at 10% solar metallicity to metal-free stars, we track the individual systems until the coalescence of the compact remnants. We estimate the contribution of primordial stars to the merger rate density and to the detection rate of the Advanced Laser Interferometer Gravitational-Wave Observatory (aLIGO). Owing to their higher masses, the remnants of primordial stars produce strong GW signals, even if their contribution in number is relatively small. We find a probability of $\gtrsim1\%$ that the current detection GW150914 is of primordial origin. We estimate that aLIGO will detect roughly 1 primordial BH-BH merger per year for the final design sensitivity, although this rate depends sensitively on the primordial initial mass function (IMF). Turning this around, the detection of black hole mergers with a total binary mass of $\sim 300\,\mathrm{M}_\odot$ would enable us to constrain the primordial IMF.

The detection of $\rm z>6$ quasars reveals the existence of supermassive black holes of a few $\rm 10^9~M_{\odot}$. One of the potential pathways to explain their formation in the infant universe is the so-called direct collapse model which provides massive seeds of $\rm 10^5-10^6~M_{\odot}$. An isothermal direct collapse mandates that halos should be of a primordial composition and the formation of molecular hydrogen remains suppressed in the presence of a strong Lyman Werner flux. In this study, we perform high resolution cosmological simulations for two massive primordial halos employing a detailed chemical model which includes $\rm H^-$ cooling as well as realistic opacities for both the bound-free $\rm H^-$ emission and the Rayleigh scattering of hydrogen atoms. We are able to resolve the collapse up to unprecedentedly high densities of $\rm \sim 10^{-3}~g/cm^3$ and to scales of about $\rm 10^{-4}$ AU. Our results show that the gas cools down to $\rm \sim $ 5000 K in the presence of $\rm H^-$ cooling, and induces fragmentation at scales of about 8000 AU in one of the two simulated halos, which may lead to the formation of a binary. In addition, fragmentation also occurs on the AU scale in one of the halos but the clumps are expected to merge on short time scales. Our results confirm that $\rm H^-$ cooling does not prevent the formation of a supermassive star and the trapping of cooling radiation stabilises the collapse on small scales.

High redshift quasars at z>6 have masses up to ~$10^9$ M$_\odot$. One of the pathways to their formation includes direct collapse of gas, forming a supermassive star, precursor of the black hole seed. The conditions for direct collapse are more easily achievable in metal-free haloes, where atomic hydrogen cooling operates and molecular hydrogen (H2) formation is inhibited by a strong external UV flux. Above a certain value of UV flux (J_crit), the gas in a halo collapses isothermally at ~$10^4$ K and provides the conditions for supermassive star formation. However, H2 can self-shield, reducing the effect of photodissociation. So far, most numerical studies used the local Jeans length to calculate the column densities for self-shielding. We implement an improved method for the determination of column densities in 3D simulations and analyse its effect on the value of J_crit. This new method captures the gas geometry and velocity field and enables us to properly determine the direction-dependent self-shielding factor of H2 against photodissociating radiation. We find a value of J_crit that is a factor of two smaller than with the Jeans approach (~2000 J_21 vs. ~4000 J_21). The main reason for this difference is the strong directional dependence of the H2 column density. With this lower value of J_crit, the number of haloes exposed to a flux >J_crit is larger by more than an order of magnitude compared to previous studies. This may translate into a similar enhancement in the predicted number density of black hole seeds.

We present a new near-field cosmological probe of the initial mass function (IMF) of the first stars. Specifically, we constrain the lower-mass limit of the Population III (Pop III) IMF with the total number of stars in large, unbiased surveys of the Milky Way. We model the early star formation history in a Milky Way-like halo with a semi-analytic approach, based on Monte-Carlo sampling of dark matter merger trees, combined with a treatment of the most important feedback mechanisms. Assuming a logarithmically flat Pop III IMF and varying its low mass limit, we derive the number of expected survivors of these first stars, using them to estimate the probability to detect any such Pop III fossil in stellar archaeological surveys. Following our analysis, the most promising region to find possible Pop III survivors is the stellar halo of the Milky Way, which is the best target for future surveys. We find that if no genuine Pop III survivor is detected in a sample size of $4 \times 10^6$ ($2 \times 10^7$) halo stars with well-controlled selection effects, then we can exclude the hypothesis that the primordial IMF extended down below $0.8 M_\odot$ at a confidence level of 68% (99%). With the sample size of the Hamburg/ESO survey, we can tentatively exclude Pop III stars with masses below $0.65 M_\odot$ with a confidence level of 95%, although this is subject to significant uncertainties. To fully harness the potential of our approach, future large surveys are needed that employ uniform, unbiased selection strategies for high-resolution spectroscopic follow-up.

We present a new method for estimating the H2 cooling rate in the optically thick regime in simulations of primordial star formation. Our new approach is based on the TreeCol algorithm, which projects matter distributions onto a spherical grid to create maps of column densities for each fluid element in the computational domain. We have improved this algorithm by using the relative gas velocities, to weight the individual matter contributions with the relative spectral line overlaps, in order to properly account for the Doppler effect. We compare our new method to the widely used Sobolev approximation, which yields an estimate for the column density based on the local velocity gradient and the thermal velocity. This approach generally underestimates the photon escape probability, because it neglects the density gradient and the actual shape of the cloud. We present a correction factor for the true line overlap in the Sobolev approximation and a new method based on local quantities, which fits the exact results reasonably well during the collapse of the cloud, with the error in the cooling rates always being less than 10%. Analytical fitting formulae fail at determining the photon escape probability after formation of the first protostar (error of 40%) because they are based on the assumption of spherical symmetry and therefore break down once a protostellar accretion disc has formed. Our method yields lower temperatures and hence promotes fragmentation for densities above 10^{10}/ccm at a distance of 200AU from the first protostar. Since the overall accretion rates are hardly affected by the cooling implementation, we expect Pop III stars to have lower masses in our simulations, compared to the results of previous simulations that used the Sobolev approximation.