The Square Kilometer Array is expected to provide the first tomographic observations of the neutral intergalactic medium at redshifts z>6 and pinpoint the locations of individual ionized bubbles during early stages of cosmic reionization. In scenarios where star-forming galaxies provide most of the ionizing photons required for cosmic reionization, one expects the first ionized bubbles to be centered on overdensities of such galaxies. Here, we model the properties of galaxy populations within isolated, ionized bubbles that SKA-1 should be able to resolve at z=7-10, and explore the prospects for galaxy counts within such structures with various upcoming near-infrared telescopes. We find that, for the bubbles that are within reach of SKA-1 tomography, the bubble volume is closely tied to the number of ionizing photons that have escaped from the galaxies within. In the case of galaxy-dominated reionization, galaxies are expected to turn up above the spectroscopic detection threshold of JWST and ELT in even the smallest resolvable bubbles at redshifts z=10 or below. The prospects of detecting galaxies within these structures in purely photometric surveys with Euclid, WFIRST, JWST or ELT are also discussed. While spectroscopy is preferable towards the end of reionization to provide a robust sample of bubble members, multiband imaging may be a competitive option for bubbles at z~10, due to the very small number of line-of-sight interlopers expected at similar redshifts.

Supernovae (SNe) could be powerful probes of the properties of stars and galaxies at high redshifts in future surveys. Wide fields and longer exposure times are required to offset diminishing star formation rates and lower fluxes to detect useful numbers of events at high redshift. In principle, the Large Synoptic Survey Telescope (LSST) could discover large numbers of early SNe because of its wide fields but only at lower redshifts because of its AB mag limit of ~ 24. But gravitational lensing by galaxy clusters and massive galaxies could boost flux from ancient SNe and allow LSST to detect them at earlier times. Here, we calculate detection rates for lensed SNe at z ~ 5 - 7 for LSST. We find that the LSST Wide Deep Fast survey could detect up to 120 lensed Population (Pop) I and II SNe but no lensed Pop III SNe. Deep-drilling programs in a single 10 square degree FoV could detect Pop I and II core-collapse SNe and Pop III pair-instability SNe at AB magnitudes of 27 - 28 and 26, respectively. An alternative deep survey over 80 nights with a one-year cadence could find ~ 8 Pop III SNe.

We built a grid of photoionization models and compiled already available observational emission line intensities ($\rm 1000 \: < \: \lambda(\AA) \: < \: 2000$) of confirmed star formation regions and Active Galactic Nucleus (AGNs) in order to classify five Ly$\alpha$ emitter (LAE) objects at high redshift $(5.7 \: < \: z \: < \:7.2)$. We selected objects for which at least one metal emission-line was measured. The resulting sample is composed by the objects RXCJ2248.7-4431-ID3, HSCJ233408+004403, COSY, A1703-zd6, and CR7 (clump C). The photoionization models were built assuming a Power Law (associated with the presence of an AGN), a Direct Collapse Black Hole (DCBH), and Population II stars for the ionizing source. The resulting models were then compared with observational emission-line ratios in six diagnostic diagrams to produce a spectral classification of the sample. We found that CR7 (clump C), HSCJ233408+004403 and COSY probably have a non thermal ionizing source (AGN or DCBH) while the RXC J2248.7-4431-ID3 and A1703-zd6 seem to host a stellar cluster. Detailed photoionization models were constructed to reproduce observational emission line ratios of the sample of LAEs, and to derive chemical abundances and number of ionizing photons $Q(\rm H)$ of these objects. From these models, we found metallicities in the range $(Z/Z_{\odot})=0.1-0.5$ and $\log Q(\rm H) \: > \: 53$. Values for C/O abundance ratio derived for the LAEs seem to be consistent with those derived for local star forming objects with similar metallicities, while an overabundance of N/O was found for most of the LAEs.

A fundamental quest of modern astronomy is to locate the earliest galaxies and study how they influenced the intergalactic medium a few hundred million years after the Big Bang. The abundance of star-forming galaxies is known to decline from redshifts of about 6 to 10, but a key question is the extent of star formation at even earlier times, corresponding to the period when the first galaxies might have emerged. Here we present spectroscopic observations of MACS1149-JD1, a gravitationally lensed galaxy observed when the Universe was less than four per cent of its present age. We detect an emission line of doubly ionized oxygen at a redshift of $9.1096\pm0.0006$, with an uncertainty of one standard deviation. This precisely determined redshift indicates that the red rest-frame optical colour arises from a dominant stellar component that formed about 250 million years after the Big Bang, corresponding to a redshift of about 15. Our results indicate the it may be possible to detect such early episodes of star formation in similar galaxies with future telescopes.

In this study we focus on the $z\sim6.6$ Lyman-$\alpha$ CR7 consisting of clump A that is host to a potential direct collapse black hole (DCBH), and two metal enriched star forming clumps B and C. In contrast to claims that signatures of metals rule out the existence of DCBHs, we show that metal pollution of A from star forming clumps clumps B and C is inevitable, and that A can form a DCBH well before its metallicity exceeds the critical threshold of $10^{-5}-10^{-6}\ \rm Z_{\odot}$. Assuming metal mixing happens instantaneously, we derive the metallicity of A based on the star formation history of B and C. We find that treating a final accreting black hole of $10^6-10^7\ \rm M_{\odot}$ in A for nebular emission already pushes its $H_{160}$ - [3.6] and [3.6]-[4.5] colours into the 3$\sigma$ limit of observations. Hence, we show that the presence of metals in DCBH hosts is inevitable, and that it is the coevolution of the LW radiation field and metals originating from neighbouring galaxies that governs DCBH formation in a neighbouring {initially} pristine atomic cooling haloes.

Direct collapse black holes forming in pristine, atomically-cooling haloes at $z \approx 10-20$ may act as the seeds of supermassive black holes (BH) at high redshifts. In order to create a massive BH seed, the host halo needs to be prevented from forming stars. H$_2$ therefore needs to be irradiated by a large flux of Lyman-Werner (LW) UV photons in order to suppress H$_2$ cooling. A key uncertainty in this scenario is the escape fraction of LW radiation from first galaxies, the dominant source of UV photons at this epoch. To better constrain this escape fraction, we have performed radiation-hydrodynamical simulations of the growth of HII regions and their associated photodissociation regions in the first galaxies using the ZEUS-MP code. We find that the LW escape fraction crucially depends on the propagation of the ionisation front (I-front). For an R-type I-front overrunning the halo, the LW escape fraction is always larger than 95%. If the halo recombines later from the outside--in, due to a softened and weakened spectrum, the LW escape fraction in the rest-frame of the halo (the near-field) drops to zero. A detailed and careful analysis is required to analyse slowly moving, D-type I-fronts, where the escape fraction depends on the microphysics and can be as small as 3% in the near-field and 61% in the far-field or as large as 100% in both the near-field and the far-field.

Supernovae are important probes of the properties of stars at high redshifts because they can be detected at early epochs and their masses can be inferred from their light curves. Finding the first cosmic explosions in the universe will only be possible with the James Webb Space Telescope, the Wide-Field Infrared Survey Telescope and the next generation of extremely large telescopes. But strong gravitational lensing by massive clusters, like those in the Cluster Lensing and Supernova Survey with Hubble (CLASH), could reveal such events now by magnifying their flux by factors of 10 or more. We find that CLASH will likely discover at least 2 - 3 core-collapse supernovae at 5 < z < 12 and perhaps as many as ten. Future surveys of cluster lenses similar in scope to CLASH by the James Webb Space Telescope might find hundreds of these events out to z ~ 15 - 17. Besides revealing the masses of early stars, these ancient supernovae will also constrain cosmic star formation rates in the era of first galaxy formation.

While the Lyman-$\alpha$ ($\mathrm{Ly}\alpha$) emission line serves as an important tool in the study of galaxies at $z\lesssim 6$, finding Ly$\alpha$ emitters (LAE) at significantly higher redshifts has been more challenging, probably because of the increasing neutrality of the intergalactic medium above $z\sim6$. Galaxies with extremely high rest-frame Ly$\alpha$ equivalent widths, EW(Ly$\alpha$) $\gtrsim 150$ \AA{}, at $z>6$ are good candidates for Ly$\alpha$ follow-up observations, and can stand out in multiband imaging surveys because of their unusual colors. We have conducted a photometric search for such objects in the Cluster Lensing And Supernova survey with Hubble (CLASH), and report here the identification of three likely gravitationally-lensed images of a single LAE candidate at $z\sim6.3$, behind the galaxy cluster Abell 2261($z = 0.225$). In the process, we also measured with Keck/MOSFIRE the first spectroscopic redshift of a multiply-imaged galaxy behind Abell 2261, at $z = 3.337$. This allows us to calibrate the lensing model, which in turn is used to study the properties of the candidate LAE. Population III galaxy spectral energy distribution (SED) model fits to the CLASH broadband photometry of the possible LAE provide a slightly better fit than Population I/II models. The best fitted model suggests intrinsic EW(Ly$\alpha$) $\approx 160$ \AA{} after absorption in the interstellar and intergalactic medium. Future spectroscopic observations will examine this prediction as well as shed more light on the morphology of this object, which indicates it may be a merger of two smaller galaxies.

Population III galaxies are predicted to exist at high redshifts and may be rendered sufficiently bright for detection with current telescopes when gravitationally lensed by a foreground galaxy cluster. Population III galaxies that exhibit strong Lya emission should furthermore be identifiable from broadband photometry because of their unusual colors. Here, we report on a search for such objects at z > 6 in the imaging data from the Cluster Lensing And Supernova survey with Hubble (CLASH), covering 25 galaxy clusters in 16 filters. Our selection algorithm returns five singly-imaged candidates with Lya-like color signatures, for which ground-based spectroscopy with current 8-10 m class telescopes should be able to test the predicted strength of the Lya line. None of these five objects have been included in previous CLASH compilations of high-redshift galaxy candidates. However, when large grids of spectral synthesis models are applied to the study of these objects, we find that only two of these candidates are significantly better fitted by Population III models than by more mundane, low-metallicity stellar populations.

The first metals in the universe are expected to form in population III stars - primordial stars consisting entirely of hydrogen and helium. However, these stars have so far remained elusive. Simulations indicate that galaxies consisting exclusively, or almost exclusively, of population III stars may form at z>6, and such galaxies may provide one of the best probes of the properties of the population III star formation mode. By fitting Yggdrasil model spectra to multiband photometry data, we have identified four population III galaxy candidates in the Cluster Lensing and Supernova survey with Hubble (CLASH). We rule out alternative mundane galaxies and low redshift interlopers through similar fits to catalogs of spectra from more mundane objects. If confirmed through spectroscopy, this would constitute the first detection of the "missing link" between the early pristine universe and the metal-enriched universe.

Dark halo substructure may reveal itself through secondary, small-scale gravitational lensing effects on light sources that are macrolensed by a foreground galaxy. Here, we explore the prospects of using Very Long Baseline Interferometry (VLBI) observations of multiply-imaged quasar jets to search for submilliarcsecond-scale image distortions produced by various forms of dark substructures in the 1e3-1e8 Msolar mass range. We present lensing simulations relevant for the angular resolutions attainable with the existing European VLBI Network (EVN), the global VLBI array, and an upcoming observing mode in which the Atacama Large Millimeter Array (ALMA) is connected to the global VLBI array. While observations of this type would not be sensitive to standard cold dark matter subhalos, they can be used to detect more compact forms of halo substructure predicted in alternative structure formation scenarios. By mapping ~5 strongly lensed systems, it should be possible to detect or robustly rule out primordial black holes in the 1e3-1e6 Msolar mass range if they constitute >1% percent of the dark matter in these lenses. Ultracompact minihalos are harder to detect using this technique, but 1e6-1e8 Msolar ultracompact minihalos could in principle be detected if they constitute >10% of the dark matter.

The first population III stars are predicted to form in minihalos at a redshift of approximately 10-30. The James Webb Space Telescope (JWST), tentatively scheduled for launch in 2018, will probably be able to detect some of the first galaxies, but whether it will also be able to detect the first stars remains more doubtful. Here, we explore the prospects of detecting an isolated population III star or a small cluster of population III stars down to redshift 2 in either lensed or unlensed fields. Our calculations are based on realistic stellar atmospheres and take into account the potential flux contribution from the surrounding HII region. We find that unlensed population III stars are beyond the reach of JWST, and that even lensed population III stars will be extremely difficult to detect. However, the main problem with the latter approach is not necessarily that the lensed stars are too faint, but that their surface number densities are too low. To detect even one population III star of 60 solar masses when pointing JWST through the galaxy cluster MACS J0717.5+3745, the lensing cluster with the largest Einstein radius detected so far, the cosmic star formation rate of population III stars would need to be approximately an order of magnitude higher than predicted by the most optimistic current models.

Small galaxies consisting entirely of population III (pop III) stars may form at high redshifts, and could constitute one of the best probes of such stars. Here, we explore the prospects of detecting gravitationally lensed pop III galaxies behind the galaxy cluster J0717.5+3745 (J0717) with both the Hubble Space Telescope (HST) and the upcoming James Webb Space Telescope (JWST). By projecting simulated catalogs of pop III galaxies at z~7-15 through the J0717 magnification maps, we estimate the lensed number counts as a function of flux detection threshold. We find that the ongoing HST survey CLASH, targeting a total of 25 galaxy clusters including J0717, potentially could detect a small number of pop III galaxies if ~1% of the baryons in these systems have been converted into pop III stars. Using JWST exposures of J0717, this limit can be pushed to ~0.1% of the baryons. Ultra-deep JWST observations of unlensed fields are predicted to do somewhat worse, but will be able to probe pop III galaxies with luminosities intermediate between those detectable in HST/CLASH and in JWST observations of J0717. We also explain how current measurements of the galaxy luminosity function at z=7-10 can be used to constrain pop III galaxy models with very high star formation efficiencies (~10% of the baryons converted into pop III stars).

Population III galaxies, made partly or exclusively of metal-free stars, are predicted to exist at high redshifts and may produce very strong Lya emission. A substantial fraction of these Lya photons are likely absorbed in the intergalactic medium at z>6, but recent simulations suggest that significant Lya emission may be detectable up to z~8.5, i.e. well into the reionization epoch. Here, we argue that high-redshift population III galaxies with strong Lya emission can be identified in Hubble Space Telescope imaging data because of their unusual colours. We quantify this effect in some of the filters used in Y-band dropout searches for galaxies at z~8 and find that population III galaxies with high Lya fluxes may exhibit much bluer J-H colours at z=8-10 than any normal type of galaxy at these redshifts. This colour signature can arise even if pop III stars account for as little as ~1e-3 to ~1e-2 of the stellar mass in these galaxies. Some of the anomalously blue objects reported in current Y-band dropout samples do in fact meet the colour criteria for Lya-emitting population III galaxies.

The James Webb Space Telescope (JWST) is expected to revolutionize our understanding of the high-redshift Universe, and may be able to test the prediction that the first, chemically pristine (population III) stars formed with very high characteristic masses. Since isolated population III stars are likely to be beyond the reach of JWST, small population III galaxies may offer the best prospects of directly probing the properties of metal-free stars. Here, we present Yggdrasil, a new spectral synthesis code geared towards the first galaxies. Using this model, we explore the JWST imaging detection limits for population III galaxies and investigate to what extent such objects may be identified based on their JWST colours. We predict that JWST should be able to detect population III galaxies with stellar population masses as low as ~10^5 Msolar at z=10 in ultra deep exposures. Over limited redshift intervals, it may also be possible to use colour criteria to select population III galaxy candidates for follow-up spectroscopy. The colours of young population III galaxies dominated by direct star light can be used to probe the stellar initial mass function (IMF), but this requires almost complete leakage of ionizing photons into the intergalactic medium. The colours of objects dominated by nebular emission show no corresponding IMF sensitivity. We also note that a clean selection of population III galaxies at z~7-8 can be achieved by adding two JWST/MIRI filters to the JWST/NIRCam filter sets usually discussed in the context of JWST ultra deep fields.

Isolated population III stars are postulated to exist at approximately z=10-30 and may attain masses up to a few hundred solar masses. The James Webb Space telescope (JWST) is the next large space based infrared telescope and is scheduled for launch in 2014. Using a 6.5 meter primary mirror, it will probably be able to detect some of the first galaxies forming in the early Universe. A natural question is whether it will also be able to see any isolated population III stars. Here, we calculate the apparent broadband AB-magnitudes for 300 solar masses population III stars in JWST filters at z=10-20. Our calculations are based on realistic stellar atmospheres and take into account the potential flux contribution from the surrounding HII region. The gravitational magnification boost achieved when pointing JWST through a foreground galaxy cluster is also considered. Using this machinery, we derive the conditions required for JWST to be able to detect population III stars in isolation. We find that a detection of individual population III stars with JWST is unlikely at these redshifts. However, the main problem is not necessarily that these stars are too faint, once gravitational lensing is taken into account, but that their surface number densities are too low.

Some of the first stars could be cooler and more massive than standard stellar models would suggest, due to the effects of dark matter annihilation in their cores. It has recently been argued that such objects may attain masses in the 10^4--10^7 solar mass range, and that such supermassive dark stars should be within reach of the upcoming James Webb Space Telescope. Notwithstanding theoretical difficulties with this proposal, we argue here that some of these objects should also be readily detectable with both the Hubble Space Telescope and ground-based 8--10 m class telescopes. Existing survey data already place strong constraints on 10^7 solar mass dark stars at z~10. We show that such objects must be exceedingly rare or short-lived to have avoided detection.

The first stars in the history of the Universe are likely to form in the dense central regions of 10^5-10^6 Msolar cold dark matter halos at z=10-50. The annihilation of dark matter particles in these environments may lead to the formation of so-called dark stars, which are predicted to be cooler, larger, more massive and potentially more long-lived than conventional population III stars. Here, we investigate the prospects of detecting high-redshift dark stars with the upcoming James Webb Space Telescope (JWST). We find that dark stars at z>6 are intrinsically too faint to be detected by JWST. However, by exploiting foreground galaxy clusters as gravitational telescopes, certain varieties of cool (Teff < 30000 K) dark stars should be within reach at redshifts up to z=10. If the lifetimes of dark stars are sufficiently long, many such objects may also congregate inside the first galaxies. We demonstrate that this could give rise to peculiar features in the integrated spectra of galaxies at high redshifts, provided that dark stars make up at least 1 percent of the total stellar mass in such objects.