On 2019 August 14, the LIGO and Virgo interferometers detected a high-significance event labelled S190814bv. Preliminary analysis of the GW data suggests that the event was likely due to the merger of a compact binary system formed by a BH and a NS. ElectromagNetic counterparts of GRAvitational wave sources at the VEry Large Telescope (ENGRAVE) collaboration members carried out an intensive multi-epoch, multi-instrument observational campaign to identify the possible optical/near infrared counterpart of the event. In addition, the ATLAS, GOTO, GRAWITA-VST, Pan-STARRS and VINROUGE projects also carried out a search on this event. Our observations allow us to place limits on the presence of any counterpart and discuss the implications for the kilonova (KN) possibly generated by this NS-BH merger, and for the strategy of future searches. Altogether, our observations allow us to exclude a KN with large ejecta mass $M\gtrsim 0.1\,\mathrm{M_\odot}$ to a high ($>90\%$) confidence, and we can exclude much smaller masses in a subsample of our observations. This disfavours the tidal disruption of the neutron star during the merger. Despite the sensitive instruments involved in the campaign, given the distance of S190814bv we could not reach sufficiently deep limits to constrain a KN comparable in luminosity to AT 2017gfo on a large fraction of the localisation probability. This suggests that future (likely common) events at a few hundreds Mpc will be detected only by large facilities with both high sensitivity and large field of view. Galaxy-targeted observations can reach the needed depth over a relevant portion of the localisation probability with a smaller investment of resources, but the number of galaxies to be targeted in order to get a fairly complete coverage is large, even in the case of a localisation as good as that of this event.

In this chapter we discuss the population and spectral synthesis of stellar populations. We describe the method required to achieve such synthesis and discuss examples where inclusion of interacting binaries are vital to reproducing the properties of observed stellar systems. These examples include the Hertzsprung-Russell diagram, massive star number counts, core-collapse supernovae and the ionising radiation from stellar populations that power both nearby HII regions and the epoch of reionization. We finally offer some speculations on the future paths of research in spectral synthesis.

The recent identification of a candidate very massive 70 M(Sun) black hole is at odds with our current understanding of stellar winds and pair-instability supernovae. We investigate alternate explanations for this system by searching the BPASS v2.2 stellar and population synthesis models for those that match the observed properties of the system. We find binary evolution models that match the LB-1 system, at the reported Gaia distance, with more moderate black hole masses of 4 to 7 M(Sun). We also examine the suggestion that the binary motion may have led to an incorrect distance determination by Gaia. We find that the Gaia distance is accurate and that the binary system is consistent with the observation at this distance. Consequently it is highly improbable that the black hole in this system has the extreme mass originally suggested. Instead, it is more likely to be representative of the typical black hole binary population expected in our Galaxy.

Comparison with artificial galaxy models is essential for translating the incomplete and low signal-to-noise data we can obtain on astrophysical stellar populations to physical interpretations which describe their composition, physical properties, histories and internal conditions. In particular, this is true for distant galaxies, whose unresolved light embeds clues to their formation and evolution as well as their impact on their wider environs. Stellar population synthesis models are now used as the foundation of analysis at all redshifts, but are not without their problems. Here we review the use of stellar population synthesis models, with a focus on applications in the distant Universe.

Binary stars have been shown to have a substantial impact on the integrated light of stellar populations, particularly at low metallicity and early ages - conditions prevalent in the distant Universe. But the fraction of stars in stellar multiples as a function of mass, their likely initial periods and distribution of mass ratios are all known empirically from observations only in the local Universe. Each has associated uncertainties. We explore the impact of these uncertainties in binary parameters on the properties of integrated stellar populations, considering which properties and timescales are most susceptible to uncertainty introduced by binary fractions and whether observations of the integrated light might be sufficient to determine binary parameters. We conclude that the effects of uncertainty in the empirical binary parameter distributions are likely smaller than those introduced by metallicity and stellar population age uncertainties for observational data. We identify emission in the He II 1640 Angstrom emission line and continuum colour in the ultraviolet-optical as potential indicators of a high mass binary presence, although poorly constrained metallicity, dust extinction and degeneracies in plausible star formation history are likely to swamp any measurable signal.

Accurately determining the age of H\,{ii} regions and the stars they host is as important as it is challenging. Historically the most popular method has been isochrone fitting to Hertzsprung-Russell Diagrams or Colour-Magnitude Diagrams. Here we introduce a different method for age determination using BPASS and hoki. We infer the most likely ages of the regions D118 and D119 NGC 300 to be log(age/years)=6.86$^{0.05}_{0.16}$ and we also deduce stellar mass and number counts by comparison with the BPASS models. We compare how our binary and single star models perform and find that the latter are unable to predict 20 per cent ($\pm$ 10 per cent) of our sample. We also discuss how results obtained from isochrone fitting would differ. We conclude that ages could be underestimated by ~0.2 dex and that the limitations of the isochrone method is not solely due to the lack of binary stars. We propose that the method presented here is more reliable and more widely applicable since it can be used on smaller samples. Alongside this study, we release new hoki features to allow easy implementation of this method.

We report the results of optical follow-up observations of 29 gravitational-wave triggers during the first half of the LIGO-Virgo Collaboration (LVC) O3 run with the Gravitational-wave Optical Transient Observer (GOTO) in its prototype 4-telescope configuration (GOTO-4). While no viable electromagnetic counterpart candidate was identified, we estimate our 3D (volumetric) coverage using test light curves of on- and off-axis gamma-ray bursts and kilonovae. In cases where the source region was observable immediately, GOTO-4 was able to respond to a GW alert in less than a minute. The average time of first observation was $8.79$ hours after receiving an alert ($9.90$ hours after trigger). A mean of $732.3$ square degrees were tiled per event, representing on average $45.3$ per cent of the LVC probability map, or $70.3$ per cent of the observable probability. This coverage will further improve as the facility scales up alongside the localisation performance of the evolving gravitational-wave detector network. Even in its 4-telescope prototype configuration, GOTO is capable of detecting AT2017gfo-like kilonovae beyond 200 Mpc in favourable observing conditions. We cannot currently place meaningful electromagnetic limits on the population of distant ($\hat{D}_L = 1.3$ Gpc) binary black hole mergers because our test models are too faint to recover at this distance. However, as GOTO is upgraded towards its full 32-telescope, 2 node (La Palma & Australia) configuration, it is expected to be sufficiently sensitive to cover the predicted O4 binary neutron star merger volume, and will be able to respond to both northern and southern triggers.

We now know that a large number of stars are born in multiple systems. Additionally, more than 70% of massive stars are found in close binary systems, meaning that they will interact over the course of their lifetime. This has strong implications for their evolution as well as the transients (e.g supernovae) and the potential gravitational wave progenitors they produce. Therefore, in order to understand and correctly interpret astronomical observations of stellar populations, we must use theoretical models able to account for the effects of binary stars. This is the case of the Binary Population and Spectral Synthesis code (BPASS), which has been a staple of the field for over 10 years. As is the case for most other theoretical models, the data products of BPASS are large, varied and complex. As a result, their use requires a level of expertise that is not immediately accessible to a wider community that may hold key observational data. The goal of hoki is to bridge the gap between observation and theory, by providing a set of tools to make BPASS data easily accessible and facilitate analysis. The use of Python is deliberate as it is a ubiquitous language within Astronomy. This allows BPASS results to be used naturally within the pre-existing workflow of most astronomers.

Observations of star-forming galaxies in the distant Universe have confirmed the importance of massive stars in shaping galaxy emission and evolution. Distant stellar populations are unresolved, and the limited data available must be interpreted in the context of stellar population models. Understanding these populations, and their evolution with age and heavy element content is key to interpreting processes such as supernovae, cosmic reionization and the chemical enrichment of the Universe. With the upcoming launch of JWST and observations of galaxies within a billion years of the Big Bang, the uncertainties in modelling massive stars - particularly their interactions with binary companions - are becoming increasingly important to our interpretation of the high redshift Universe. In turn, observations of distant stellar populations provide ever stronger tests against which to gauge the success of, and flaws in, current massive star models. Here we briefly review the current status binary stellar population synthesis.

We compare the impacts of uncertainties in both binary population synthesis models and the cosmic star formation history on the predicted rates of Gravitational Wave compact binary merger (GW) events. These uncertainties cause the predicted rates of GW events to vary by up to an order of magnitude. Varying the volume-averaged star formation rate density history of the Universe causes the weakest change to our predictions, while varying the metallicity evolution has the strongest effect. Double neutron-star merger rates are more sensitive to assumed neutron-star kick velocity than the cosmic star formation history. Varying certain parameters affects merger rates in different ways depending on the mass of the merging compact objects; thus some of the degeneracy may be broken by looking at all the event rates rather than restricting ourselves to one class of mergers.

Long-duration gamma-ray bursts (GRBs) are understood to be the final fate for a subset of massive, stripped envelope, rapidly rotating stars. Beyond this, our knowledge of the progenitor systems is limited. Using the BPASS (Binary Population and Spectral Synthesis) stellar evolution models, we investigate the possibility that some massive stars in binaries can maintain the angular momentum required for jet production, while still loosing their outer envelope through winds or binary interactions. We find that a total hydrogen mass of less than 0.0005 Msun and a helium ejecta mass fraction of less than 0.20 provide the best thresholds for the supernova type II/Ibc and Ib/Ic divisions respectively. Tidal interactions in binaries are accounted for by applying a tidal algorithm to post-process the stellar evolution models output by BPASS. We show that the observed volumetric gamma-ray burst rate evolution can be recreated using two distinct pathways and plausible distributions for burst parameters. In the first pathway, stars are spun up by mass accretion into a quasi-homogeneous state. In the second, tides maintain rotation where otherwise the star would spin down. Both lead to type Ic supernova progenitors, and a metallicity distribution consistent with the GRB host galaxy population. The inferred core angular momentum threshold for jet production is consistent with theoretical requirements for collapsars, given the assumptions made in our model. We can therefore reproduce several aspects of core collapse supernova/GRB observation and theory simultaneously. We discuss the predicted observable properties of GRB progenitors and their surviving companions.

The calibrations linking observed luminosities to the star formation rate depend on the assumed stellar population synthesis model, initial mass function, star formation and metal enrichment history, and whether reprocessing by dust and gas is included. Consequently the shape and normalisation of the inferred cosmic star formation history is sensitive to these assumptions. Using v2.2.1 of the Binary Population and Spectral Synthesis (\bpass) model we determine a new set of calibration coefficients for the ultraviolet, thermal-infrared, and, hydrogen recombination lines. These ultraviolet and thermal infrared coefficients are 0.15-0.2 dex higher than those widely utilised in the literature while the H$\alpha$ coefficient is $\sim 0.35$ dex larger. These differences arise in part due to the inclusion binary evolution pathways but predominantly reflect an extension in the IMF to 300 $M_{\odot}$ and a change in the choice of reference metallicity. We use these new coefficients to recalibrate the cosmic star formation history, and find improved agreement between the integrated cosmic star formation history and the in-situ measured stellar mass density as a function of redshift. However, these coefficients produce new tension between star formation rate densities inferred from the ultraviolet and thermal-infrared and those from H$\alpha$.

We use the results of a supernova light-curve population synthesis to predict the range of possible supernova light curves arising from a population of single-star progenitors that lead to type IIP supernovae. We calculate multiple models varying the initial mass, explosion energy, nickel mass and nickel mixing and then compare these to type IIP supernovae with detailed light curve data and pre-explosion imaging progenitor constraints. Where a good fit is obtained to observations, we are able to achieve initial progenitor and nickel mass estimates from the supernova lightcurve that are comparable in precision to those obtained from progenitor imaging. For two of the eleven IIP supernovae considered our fits are poor, indicating that more progenitor models should be included in our synthesis or that our assumptions, regarding factors such as stellar mass loss rates or the rapid final stages of stellar evolution, may need to be revisited in certain cases. Using the results of our analysis we are able to show that most of the type IIP supernovae have an explosion energy of the order of log(E_exp/ergs)=50.52+/-0.10 and that both the amount of nickel in the supernovae and the amount of mixing may have a dependence on initial progenitor mass.

The number of long gamma-ray bursts (GRBs) known to have occurred in the distant Universe (z greater than 5) is small (approx 15), however these events provide a powerful way of probing star formation at the onset of galaxy evolution. In this paper, we present the case for GRB100205A being a largely overlooked high-redshift event. While initially noted as a high-z candidate, this event and its host galaxy have not been explored in detail. By combining optical and near-infrared Gemini afterglow imaging (at t less than 1.3 days since burst) with deep late-time limits on host emission from the Hubble Space Telescope, we show that the most likely scenario is that GRB100205A arose in the redshift range 4-8. GRB100205A is an example of a burst whose afterglow, even at 1 hour post-burst, could only be identified by 8m class IR observations, and suggests that such observations of all optically dark bursts may be necessary to significantly enhance the number of high-redshift GRBs known.

We study and compare different examples of stellar evolutionary synthesis input parameters used to produce photoionisation model grids using the MAPPINGS V modelling code. The aim of this study is to (a) explore the systematic effects of various stellar evolutionary synthesis model parameters on the interpretation of emission lines in optical strong-line diagnostic diagrams, (b) characterise the combination of parameters able to reproduce the spread of local galaxies located in the star-forming region in the Sloan Digital Sky Survey, and (c) investigate the emission from extremely metal-poor galaxies using photoionisation models. We explore and compare the stellar input ionising spectrum (stellar population synthesis code [Starburst99, SLUG, BPASS], stellar evolutionary tracks, stellar atmospheres, star-formation history, sampling of the initial mass function) as well as parameters intrinsic to the H II region (metallicity, ionisation parameter, pressure, H II region boundedness). We also perform a comparison of the photoionisation codes MAPPINGS and CLOUDY. On the variations in the ionising spectrum model parameters, we find that the differences in strong emission-line ratios between varying models for a given input model parameter are small, on average ~0.1 dex. An average difference of ~0.1 dex in emission-line ratio is also found between models produced with MAPPINGS and CLOUDY. Large differences between the emission-line ratios are found when comparing intrinsic H II region parameters. We find that low-metallicity galaxies are better explained by a density-bounded H II region and higher pressures better encompass the spread of galaxies at high redshift.

We present a study of 21 dark gamma-ray burst (GRB) host galaxies, predominantly using X-ray afterglows obtained with the Chandra X-Ray Observatory (CXO) to precisely locate the burst in deep Hubble Space Telescope (HST) imaging of the burst region. The host galaxies are well-detected in F160W in all but one case and in F606W imaging in approx 60 per cent of cases. We measure magnitudes and perform a morphological analysis of each galaxy. The asymmetry, concentration and ellipticity of the dark burst hosts are compared against the host galaxies of optically bright GRBs. In agreement with other studies, we find that dark GRB hosts are redder and more luminous than the bulk of the GRB host population. The distribution of projected spatial offsets for dark GRBs from their host galaxy centroids is comparable to that of optically-bright bursts. The dark GRB hosts are physically larger, more massive and redder, but are morphologically similar to the hosts of bright GRBs in terms of concentration and asymmetry. Our analysis constrains the fraction of high redshift (z greater than 5) GRBs in the sample to approx 14 per cent, implying an upper limit for the whole long-GRB population of less than 4.4 per cent. If dust is the primary cause of afterglow darkening amongst dark GRBs, the measured extinction may require a clumpy dust component in order to explain the observed offset and ellipticity distributions.

Observations of both galaxies in the distant Universe and local starbursts are showing increasing evidence for very hard ionizing spectra that stellar population synthesis models struggle to reproduce. Here we explore the effects of the assumed stellar initial mass function (IMF) on the ionizing photon output of young populations at wavelengths below key ionization energy thresholds. We use a custom set of binary population and spectral synthesis (BPASS) models to explore the effects of IMF assumptions as a function of metallicity, IMF slope, upper mass limit, IMF power law break mass and sampling. We find that while the flux capable of ionizing hydrogen is only weakly dependent on IMF parameters, the photon flux responsible for the He II and O VI lines is far more sensitive to assumptions. In our current models this flux arises primarily from helium and Wolf-Rayet stars which have partially or fully lost their hydrogen envelopes. The timescales for formation and evolution of both Wolf Rayet stars and helium dwarfs, and hence inferred population age, are affected by choice of model IMF. Even the most extreme IMFs cannot reproduce the He II ionizing flux observed in some high redshift galaxies, suggesting a source other than stellar photospheres. We caution that detailed interpretation of features in an individual galaxy spectrum is inevitably going to be subject to uncertainties in the IMF of its contributing starbursts. We remind the community that the initial mass function is fundamentally a statistical construct, and that stellar population synthesis models are most effective when considering entire galaxy populations rather than individual objects.

We present results of a supernova light-curve population synthesis, predicting the range of possible supernova lightcurves arising from a population of progenitor stars that include interacting binary systems. We show that the known diversity of supernova lightcurves can be interpreted as arising from binary interactions. Given detailed models of the progenitor stars, we are able to the determine what parameters within these stars determine the shape of their supernova lightcurve. The primary factors are the mass of supernova ejecta and the mass of hydrogen in the final progenitor. We find that there is a continuum of light-curve behaviour from type IIP, IIL to IIb supernovae related to the range of hydrogen and ejecta masses. Most type IIb supernovae arise from a relatively narrow range of initial masses from 10 to 15 M$_{\odot}$. We also find a few distinct lightcurves that are the result of stellar mergers.

Gravitational wave transients, resulting from the merger of two stellar remnants, are now detectable. The properties and rates of these directly relates to the stellar population which gave rise to their progenitors, and thus to other, electromagnetic transients which result from stellar death. We aim to estimate simultaneously the event rates and delay time distribution of gravitational wave-driven compact object mergers together with the rates of core collapse and thermonuclear supernovae within a single consistent stellar population synthesis paradigm. We combine event delay-time distributions at different metallicities from the Binary Population and Spectral Synthesis (BPASS) models with an analytic model of the volume-averaged cosmic star formation rate density and chemical evolution to determine the volume-averaged rates of each event rate at the current time. We estimate rates in excellent agreement with extant observational constraints on core-collapse supernovae, thermonuclear supernovae and long GRBs. We predict rates for gravitational wave mergers based on the same stellar populations, and find rates consistent with current LIGO estimates. We note that tighter constraints on the rates of these events will be required before it is possible to determine their redshift evolution, progenitor metallicity dependence or constrain uncertain aspects of stellar evolution.

The binary neutron star merger GW170817 was the first multi-messenger event observed in both gravitational and electromagnetic waves. The electromagnetic signal began approximately 2 seconds post-merger with a weak, short burst of gamma-rays, which was followed over the next hours and days by the ultraviolet, optical and near-infrared emission from a radioactively- powered kilonova. Later, non-thermal rising X-ray and radio emission was observed. The low luminosity of the gamma-rays and the rising non-thermal flux from the source at late times could indicate that we are outside the opening angle of the beamed relativistic jet. Alternatively, the emission could be arising from a cocoon of material formed from the interaction between a jet and the merger ejecta. Here we present late-time optical detections and deep near-infrared limits on the emission from GW170817 at 110 days post-merger. Our new observations are at odds with expectations of late-time emission from kilonova models, being too bright and blue. Instead, the emission arises from the interaction between the relativistic ejecta of GW170817 and the interstellar medium. We show that this emission matches the expectations of a Gaussian structured relativistic jet, which would have launched a high luminosity short GRB to an aligned observer. However, other jet structure or cocoon models can also match current data - the future evolution of the afterglow will directly distinguish the origin of the emission.

Massive stars are the main objects that illuminate H II regions and they evolve quickly to end their lives in core-collapse supernovae (CCSNe). Thus it is important to investigate the association between CCSNe and H II regions. In this paper, we present emission line diagnostics of the stellar populations around nearby CCSNe, that include their host H II regions, from the PMAS/PPAK Integral-field Supernova hosts COmpilation (PISCO). We then use BPASS stellar population models to determine the age, metallicity and gas parameters for H II regions associated with CCSNe, contrasting models that consider either single star evolution alone or incorporate interacting binaries. We find binary-star models, that allow for ionizing photon loss, provide a more realistic fit to the observed CCSN hosts with metallicities that are closer to those derived from the oxygen abundance in O3N2. We also find that type II and type Ibc SNe arise from progenitor stars of similar age, mostly from 7 to 45 Myr, which corresponds to stars with masses < 20 solar mass . However these two types SNe have little preference in their host environment metallicity measured by oxygen abundance or in progenitor initial mass. We note however that at lower metallicities supernovae are more likely to be of type II.

Molecular species, most frequently H_2, are present in a small, but growing, number of gamma-ray burst (GRB) afterglow spectra at redshifts z~2-3, detected through their rest-frame UV absorption lines. In rare cases, lines of vibrationally excited states of H_2 can be detected in the same spectra. The connection between afterglow line-of-sight absorption properties of molecular (and atomic) gas, and the observed behaviour in emission of similar sources at low redshift, is an important test of the suitability of GRB afterglows as general probes of conditions in star formation regions at high redshift. Recently, emission lines of carbon monoxide have been detected in a small sample of GRB host galaxies, at sub-mm wavelengths, but no searches for H_2 in emission have been reported yet. In this paper we perform an exploratory search for rest-frame K band rotation-vibrational transitions of H_2 in emission, observable only in the lowest redshift GRB hosts (z<0.22). Searching the data of four host galaxies, we detect a single significant rotation-vibrational H_2 line candidate, in the host of GRB 031203. Re-analysis of Spitzer mid-infrared spectra of the same GRB host gives a single low significance rotational line candidate. The (limits on) line flux ratios are consistent with those of blue compact dwarf galaxies in the literature. New instrumentation, in particular on the JWST and the ELT, can facilitate a major increase in our understanding of the H_2 properties of nearby GRB hosts, and the relation to H_2 absorption in GRBs at higher redshift.

The Lyman-alpha emission line is the strongest feature in the spectrum of most high redshift galaxies, and is typically observed as being highly asymmetric due to galactic outflows. Quantifying this asymmetry is challenging. Here, we explore how measurements of one parameterisation, Lyman-alpha skewness, are affected by instrumental resolution and detection signal-to-noise and thus whether this can be extended throughout the archive. We model these effects through simulated lines and apply our derived corrections to existing archival data sets (including sources observed with FORS2 and DEIMOS) to reconstruct the intrinsic line emission parameters. We find a large uncertainty in parameter reconstruction at low resolutions (R<3000) and high skew values, as well as substantial random errors resulting from the masking of sky lines. We suggest that interpretations of spectral line asymmetry should be made with caution, while a simpler parametrization, like B/R (blue-red flux asymmetry), is likely to yield more robust results. We see a possible weak trend in velocity width with mass, although there is no evidence in the data for a reliable correlation of skew with galaxy mass, star formation rate or age at z=4-5. Using our results, we investigate the possibilities of recovering emission line asymmetry with current and future instruments.

We report the detection of the radio afterglow of a long gamma-ray burst (GRB) 111005A at 5-345 GHz, including the very long baseline interferometry observations with the positional error of 0.2 mas. The afterglow position is coincident with the disk of a galaxy ESO 580-49 at z= 0.01326 (~1" from its center), which makes GRB 111005A the second closest GRB known to date, after GRB 980425. The radio afterglow of GRB 111005A was an order of magnitude less luminous than those of local low-luminosity GRBs, and obviously than those of cosmological GRBs. The radio flux was approximately constant and then experienced an unusually rapid decay a month after the GRB explosion. Similarly to only two other GRBs, we did not find the associated supernovae (SN), despite deep near- and mid-infrared observations 1-9 days after the GRB explosion, reaching ~20 times fainter than other SNe associated with GRBs. Moreover, we measured twice solar metallicity for the GRB location. The low gamma-ray and radio luminosities, rapid decay, lack of a SN, and super-solar metallicity suggest that GRB 111005A represents a different rare class of GRBs than typical core-collapse events. We modelled the spectral energy distribution of the GRB 111005A host finding that it is a dwarf, moderately star-forming galaxy, similar to the host of GRB 980425. The existence of two local GRBs in such galaxies is still consistent with the hypothesis that the GRB rate is proportional to the cosmic star formation rate (SFR) density, but suggests that the GRB rate is biased towards low SFRs. Using the far-infrared detection of ESO 580-49, we conclude that the hosts of both GRBs 111005A and 980425 exhibit lower dust content than what would be expected from their stellar masses and optical colours.

Determining the properties of old stellar populations (those with age >1 Gyr) has long involved the comparison of their integrated light, either in the form of photometry or spectroscopic indexes, with empirical or synthetic templates. Here we reevaluate the properties of old stellar populations using a new set of stellar population synthesis models, designed to incorporate the effects of binary stellar evolution pathways as a function of stellar mass and age. We find that single-aged stellar population models incorporating binary stars, as well as new stellar evolution and atmosphere models, can reproduce the colours and spectral indices observed in both globular clusters and quiescent galaxies. The best fitting model populations are often younger than those derived from older spectral synthesis models, and may also lie at slightly higher metallicities.

The part played by stars in the ionization of the intergalactic medium remains an open question. A key issue is the proportion of the stellar ionizing radiation that escapes the galaxies in which it is produced. Spectroscopy of gamma-ray burst afterglows can be used to determine the neutral hydrogen column-density in their host galaxies and hence the opacity to extreme ultra-violet radiation along the lines-of-sight to the bursts. Thus, making the reasonable assumption that long-duration GRB locations are representative of the sites of massive stars that dominate EUV production, one can calculate an average escape fraction of ionizing radiation in a way that is independent of galaxy size, luminosity or underlying spectrum. Here we present a sample of NH measures for 138 GRBs in the range 1.6<z<6.7 and use it to establish an average escape fraction at the Lyman limit of <fesc>~0.005, with a 98% confidence upper limit of ~0.015. This analysis suggests that stars provide a small contribution to the ionizing radiation budget of the IGM at z<5, where the bulk of the bursts lie. At higher redshifts, z>5, firm conclusions are limited by the small size of the GRB sample, but any decline in average HI column-density seems to be modest. We also find no indication of a significant correlation of NH with galaxy UV luminosity or host stellar mass, for the subset of events for which these are available. We discuss in some detail a number of selection effects and potential biases. Drawing on a range of evidence we argue that such effects, while not negligible, are unlikely to produce systematic errors of more than a factor ~2, and so would not affect the primary conclusions. Given that many GRB hosts are low metallicity, high specific star-formation rate, dwarf galaxies, these results present a particular problem for the hypothesis that such galaxies dominated the reionization of the universe.

We identify and explore the properties of an infrared-bright gamma-ray burst (GRB) host population. Candidate hosts are selected by coincidence with sources in WISE, with matching to random coordinates and a false alarm probability analysis showing that the contamination fraction is approx 0.5. This methodology has already identified the host galaxy of GRB 080517. We combine survey photometry from Pan-STARRS, SDSS, APASS, 2MASS, GALEX and WISE with our own WHT/ACAM and VLT/X- shooter observations to classify the candidates and identify interlopers. Galaxy SED fitting is performed using MAGPHYS, in addition to stellar template fitting, yielding 13 possible IR-bright hosts. A further 7 candidates are identified from previously published work. We report a candidate host for GRB 061002, previously unidentified as such. The remainder of the galaxies have already been noted as potential hosts. Comparing the IR-bright population properties including redshift z, stellar mass M*, star formation rate SFR and V-band attenuation Av to GRB host catalogues in the literature, we find that the infrared-bright population is biased toward low z, high M* and high Av. This naturally arises from their initial selection - local and dusty galaxies are more likely to have the required IR flux to be detected in WISE. We conclude that while IR-bright GRB hosts are not a physically distinct class, they are useful for constraining existing GRB host populations, particularly for long GRBs.

THESEUS is a space mission concept aimed at exploiting Gamma-Ray Bursts for investigating the early Universe and at providing a substantial advancement of multi-messenger and time-domain astrophysics. These goals will be achieved through a unique combination of instruments allowing GRB and X-ray transient detection over a broad field of view (more than 1sr) with 0.5-1 arcmin localization, an energy band extending from several MeV down to 0.3 keV and high sensitivity to transient sources in the soft X-ray domain, as well as on-board prompt (few minutes) follow-up with a 0.7 m class IR telescope with both imaging and spectroscopic capabilities. THESEUS will be perfectly suited for addressing the main open issues in cosmology such as, e.g., star formation rate and metallicity evolution of the inter-stellar and intra-galactic medium up to redshift $\sim$10, signatures of Pop III stars, sources and physics of re-ionization, and the faint end of the galaxy luminosity function. In addition, it will provide unprecedented capability to monitor the X-ray variable sky, thus detecting, localizing, and identifying the electromagnetic counterparts to sources of gravitational radiation, which may be routinely detected in the late '20s / early '30s by next generation facilities like aLIGO/ aVirgo, eLISA, KAGRA, and Einstein Telescope. THESEUS will also provide powerful synergies with the next generation of multi-wavelength observatories (e.g., LSST, ELT, SKA, CTA, ATHENA).

We present numerical models of the nebular emission from H II regions around young stellar populations over a range of compositions and ages. The synthetic stellar pop- ulations include both single stars and interacting binary stars. We compare these models to the observed emission lines of 254 HII regions of 13 nearby spiral galax- ies and 21 dwarf galaxies drawn from archival data. The models are created using the combination of the Binary Population and Spectral Synthesis (bpass) code with the photoionization code cloudy to study the differences caused by the inclusion of interacting binary stars in the stellar population. We obtain agreement with the observed emission line ratios from the nearby star-forming regions and discuss the effect of binary star evolution pathways on the nebular ionization of H II regions. We find that at population ages above 10 Myr, single-star models rapidly decrease in flux and ionization strength, while binary-star model still produce strong flux and high O III/H\b{eta} ratios. Our models can reproduce the metallicity of H II regions from spiral galaxies but we find higher metallicities than previously estimated for the H II regions from dwarf galaxies. Comparing the equivalent width of H\b{eta} emission between mod- els and observations, we find that accounting for ionizing photon leakage can affect age estimates for HII regions. When it is included, the typical age derived for HII regions is 5Myr from single-star models, and up to 10Myr with binary-star models. This is due to the existence of binary-star evolution pathways which produce more hot WR and helium stars at older ages. For future reference, we calculate new BPASS binary maximal starburst lines as a function of metallicity, and for the total model population, and present these in Appendix A.

AR Scorpii is unique amongst known white dwarf binaries in showing powerful pulsations extending to radio frequencies. Here we aim to investigate the multi-frequency radio emission of AR Sco in detail, in order to constrain its origin and emission mechanisms. We present interferometric radio frequency imaging of AR Sco at 1.5, 5 and 9 GHz, analysing the total flux and polarization behaviour of this source at high time resolution (10, 3 and 3 s), across a full 3.6 hr orbital period in each band. We find strong modulation of the radio flux on the orbital period and the orbital sideband of the white dwarf's spin period (also known as the "beat" period). This indicates that, like the optical flux, the radio flux arises predominantly from on or near the inner surface of the M-dwarf companion star. The beat-phase pulsations of AR Sco decrease in strength with decreasing frequency. They are strongest at 9 GHz and at an orbital phase ~0.5. Unlike the optical emission from this source, radio emission from AR Sco shows weak linear polarization but very strong circular polarization, reaching ~30% at an orbital phase ~0.8. We infer the probable existence of a non-relativistic cyclotron emission component, which dominates at low radio frequencies. Given the required magnetic fields, this also likely arises from on or near the M-dwarf.

We combine inferred galaxy properties from a semi-analytic galaxy evolution model incorporating dark matter halo merger trees with new estimates of supernova and gamma ray burst rates as a function of metallicity from stellar population synthesis models incorporating binary interactions. We use these to explore the stellar mass fraction of galaxies irradiated by energetic astrophysical transients and its evolution over cosmic time, and thus the fraction which is potentially habitable by life like our own. We find that 18 per cent of the stellar mass in the Universe is likely to have been irradiated within the last 260 Myr, with GRBs dominating that fraction. We do not see a strong dependence of irradiated stellar mass fraction on stellar mass or richness of the galaxy environment. We consider a representative merger tree as a Local Group analogue, and find that there are galaxies at all masses which have retained a high habitable fraction (>40 per cent) over the last 6 Gyr, but also that there are galaxies at all masses where the merger history and associated star formation have rendered galaxies effectively uninhabitable. This illustrates the need to consider detailed merger trees when evaluating the cosmic evolution of habitability.

We present Hubble Space Telescope and Chandra imaging, combined with Very Large Telescope MUSE integral field spectroscopy of the counterpart and host galaxy of the first binary neutron star merger detected via gravitational wave emission by LIGO & Virgo, GW170817. The host galaxy, NGC 4993, is an S0 galaxy at z=0.009783. There is evidence for large, face-on spiral shells in continuum imaging, and edge-on spiral features visible in nebular emission lines. This suggests that NGC 4993 has undergone a relatively recent (<1 Gyr) ``dry'' merger. This merger may provide the fuel for a weak active nucleus seen in Chandra imaging. At the location of the counterpart, HST imaging implies there is no globular or young stellar cluster, with a limit of a few thousand solar masses for any young system. The population in the vicinity is predominantly old with <1% of any light arising from a population with ages <500 Myr. Both the host galaxy properties and those of the transient location are consistent with the distributions seen for short-duration gamma-ray bursts, although the source position lies well within the effective radius (r_e ~ 3 kpc), providing an r_e-normalized offset that is closer than ~90% of short GRBs. For the long delay time implied by the stellar population, this suggests that the kick velocity was significantly less than the galaxy escape velocity. We do not see any narrow host galaxy interstellar medium features within the counterpart spectrum, implying low extinction, and that the binary may lie in front of the bulk of the host galaxy.

The Binary Population and Spectral Synthesis (BPASS) suite of binary stellar evolution models and synthetic stellar populations provides a framework for the physically motivated analysis of both the integrated light from distant stellar populations and the detailed properties of those nearby. We present a new version 2.1 data release of these models, detailing the methodology by which BPASS incorporates binary mass transfer and its effect on stellar evolution pathways, as well as the construction of simple stellar populations. We demonstrate key tests of the latest BPASS model suite demonstrating its ability to reproduce the colours and derived properties of resolved stellar populations, including well- constrained eclipsing binaries. We consider observational constraints on the ratio of massive star types and the distribution of stellar remnant masses. We describe the identification of supernova progenitors in our models, and demonstrate a good agreement to the properties of observed progenitors. We also test our models against photometric and spectroscopic observations of unresolved stellar populations, both in the local and distant Universe, finding that binary models provide a self-consistent explanation for observed galaxy properties across a broad redshift range. Finally, we carefully describe the limitations of our models, and areas where we expect to see significant improvement in future versions.

We present late-time follow-up of the relativistic tidal disruption flare candidate Swift J1112.2-8238. We confirm the previously determined redshift of $z=0.8900\pm0.0005$ based on multiple emission line detections. {\em HST} imaging of the host galaxy indicates a complex and distorted morphology with at least two spatially distinct components. These are offset in velocity space by less than 350\,km\,s$^{-1}$ in VLT/X-Shooter observations, suggesting that the host is undergoing interaction with another galaxy. The transient position is consistent to 2.2$\sigma$ with the centre of a bulge-like component at a distance of 1.1$\pm$0.5\,kpc from its centre. Luminous, likely variable radio emission has also been observed, strengthening the similarities between Swift J1112.2-8238 and other previously identified relativistic tidal disruption flares. While the transient location is $\sim2\sigma$ from the host centroid, the disrupted nature of the host may provide an explanation for this. The tidal disruption model remains a good description for these events.

SN 2002cx-like Type Ia supernovae (also known as SNe Iax) represent one of the most numerous peculiar SN classes. They differ from normal SNe Ia by having fainter peak magnitudes, faster decline rates and lower photospheric velocities, displaying a wide diversity in these properties. We present both integral-field and long-slit visual-wavelength spectroscopy of the host galaxies and explosion sites of SNe Iax to provide constraints on their progenitor formation scenarios. The SN Iax explosion site metallicity distribution is similar to that of core-collapse (CC) SNe and metal-poor compared to normal SNe Ia. Fainter members, speculated to form distinctly from brighter SN Iax, are found at a range of metallicities, extending to very metal-poor environments. Although the SN Iax explosion sites' ages and star-formation rates are comparatively older and less intense than the distribution of star forming regions across their host galaxies, we confirm the presence of young stellar populations (SP) at explosion environments for most SNe Iax, expanded here to a larger sample. Ages of the young SP (several $\times 10^{7}$ to $10^8$~yrs) are consistent with predictions for young thermonuclear and electron-capture SN progenitors. The lack of extremely young SP at the explosion sites disfavours very massive progenitors such as Wolf-Rayet explosions with significant fall-back. We find weak ionised gas in the only SN Iax host without obvious signs of star-formation. The source of the ionisation remains ambiguous but appears unlikely to be mainly due to young, massive stars.

We present radio observations at 1.5 GHz of 32 local objects selected to reproduce the physical properties of $z\sim5$ star-forming galaxies. We also report non-detections of five such sources in the sub-millimetre. We find a radio-derived star formation rate which is typically half that derived from H$\alpha$ emission for the same objects. These observations support previous indications that we are observing galaxies with a young dominant stellar population, which has not yet established a strong supernova-driven synchrotron continuum. We stress caution when applying star formation rate calibrations to stellar populations younger than 100 Myr. We calibrate the conversions for younger galaxies, which are dominated by a thermal radio emission component. We improve the size constraints for these sources, compared to previous unresolved ground-based optical observations. Their physical size limits indicate very high star formation rate surface densities, several orders of magnitude higher than the local galaxy population. In typical nearby galaxies, this would imply the presence of galaxy-wide winds. Given the young stellar populations, it is unclear whether a mechanism exists in our sources that can deposit sufficient kinetic energy into the interstellar medium to drive such outflows.

We present a new model of the optical nebular emission from HII regions by combin- ing the results of the Binary Population and Spectral Synthesis (bpass) code with the photoion- ization code cloudy (Ferland et al. 1998). We explore a variety of emission-line diagnostics of these star-forming HII regions and examine the effects of metallicity and interacting binary evo- lution on the nebula emission-line production. We compare the line emission properties of HII regions with model stellar populations, and provide new constraints on their stellar populations and supernova progenitors. We find that models including massive binary stars can successfully match all the observational constraints and provide reasonable age and mass estimation of the HII regions and supernova progenitors.

The binary system AR Sco contains an M star and the only known radio-pulsing white dwarf. The system shows emission from radio to X-rays, likely dominated by synchrotron radiation. The mechanism that produces most of this emission remains unclear. Two competing scenarios have been proposed: Collimated outflows, and direct interaction between the magnetospheres of the white dwarf and the M star. The two proposed scenarios can be tested via very long baseline interferometric radio observations. We conducted a radio observation with the Australian Long Baseline Array (LBA) on 20 Oct 2016 at 8.5 GHz to study the compactness of the radio emission. Simultaneous data with the Australian Telescope Compact Array (ATCA) were also recorded for a direct comparison of the obtained flux densities. AR Sco shows radio emission compact on milliarcsecond angular scales ($\lesssim 0.02\ \mathrm{AU}$, or $4\ \mathrm{R_{\odot}}$). The emission is orbitally modulated, with an average flux density of $\approx 6.5\ \mathrm{mJy}$. A comparison with the simultaneous ATCA data shows that no flux is resolved out on mas scales, implying that the radio emission is produced in this compact region. Additionally, the obtained radio light curves on hour timescales are consistent with the optical light curve. The radio emission in AR Sco is likely produced in the magnetosphere of the M star or the white dwarf, and we see no evidence for a radio outflow or collimated jets significantly contributing to the radio emission.

The z = 6.6 Lyman-$\alpha$ emitter 'CR7' has been claimed to have a Population III-like stellar population, or alternatively, be a candidate Direct Collapse Black Hole (DCBH). In this paper we investigate the evidence for these exotic scenarios using recently available, deeper, optical, near-infrared and mid-infrared imaging. We find strong Spitzer/IRAC detections for the main component of CR7 at 3.6 and 4.5 microns, and show that it has a blue colour ([3.6] - [4.5] $= -1.2\pm 0.3$). This colour cannot be reproduced by current Pop. III or pristine DCBH models. Instead, the results suggest that the [3.6] band is contaminated by the [OIII]4959,5007 emission line with an implied rest-frame equivalent width of EW_0 (H$\beta$ + [OIII]) $\gtrsim 2000$\AA. Furthermore, we find that new near-infrared data from the UltraVISTA survey supports a weaker He II 1640 emission line than previously measured, with EW_0 $= 40 \pm 30$\AA. For the fainter components of CR7 visible in Hubble Space Telescope imaging, we find no evidence that they are particularly red as previously claimed, and show that the derived masses and ages are considerably uncertain. In light of the likely detection of strong [OIII] emission in CR7 we discuss other more standard interpretations of the system that are consistent with the data. We find that a low-mass, narrow-line AGN can reproduce the observed features of CR7, including the lack of radio and X-ray detections. Alternatively, a young, low-metallicity (~1/200 solar) starburst, modelled including binary stellar pathways, can reproduce the inferred strength of the He II line and simultaneously the strength of the observed [OIII] emission, but only if the gas shows super-solar $\alpha$-element abundances (O/Fe ~ 5 O/Fe solar).

Gamma-ray bursts (GRBs) are powerful probes of early stars and galaxies, during and potentially even before the era of reionization. Although the number of GRBs identified at z>6 remains small, they provide a unique window on typical star-forming galaxies at that time, and thus are complementary to deep field observations. We report the identification of the optical drop-out afterglow of Swift GRB 120923A in near-infrared Gemini-North imaging, and derive a redshift of z=7.84_{-0.12}^{+0.06} from VLT/X-shooter spectroscopy. At this redshift the peak 15-150 keV luminosity of the burst was 3.2x10^52 erg/s, and in fact the burst was close to the Swift/BAT detection threshold. The X-ray and near-infrared afterglow were also faint, and in this sense it was a rather typical long-duration GRB in terms of rest-frame luminosity. We present ground- and space-based follow-up observations spanning from X-ray to radio, and find that a standard external shock model with a constant-density circumburst environment with density, n~4x10^-2 cm^-3 gives a good fit to the data. The near-infrared light curve exhibits a sharp break at t~3.4 days in the observer frame, which if interpreted as being due to a jet corresponds to an opening angle of ~5 degrees. The beaming corrected gamma-ray energy is then E_gamma~2x10^50 erg, while the beaming-corrected kinetic energy is lower, E_K~10^49 erg, suggesting that GRB 120923A was a comparatively low kinetic energy event. We discuss the implications of this event for our understanding of the high-redshift population of GRBs and their identification.

Observations of star-forming galaxies in the distant Universe (z > 2) are starting to confirm the importance of massive stars in shaping galaxy emission and evolution. Inevitably, these distant stellar populations are unresolved, and the limited data available must be interpreted in the context of stellar population synthesis models. With the imminent launch of JWST and the prospect of spectral observations of galaxies within a gigayear of the Big Bang, the uncertainties in modelling of massive stars are becoming increasingly important to our interpretation of the high redshift Universe. In turn, these observations of distant stellar populations will provide ever stronger tests against which to gauge the success of, and flaws in, current massive star models.

White dwarfs are compact stars, similar in size to Earth but ~200,000 times more massive. Isolated white dwarfs emit most of their power from ultraviolet to near-infrared wavelengths, but when in close orbits with less dense stars, white dwarfs can strip material from their companions, and the resulting mass transfer can generate atomic line and X-ray emission, as well as near- and mid-infrared radiation if the white dwarf is magnetic. However, even in binaries, white dwarfs are rarely detected at far-infrared or radio frequencies. Here we report the discovery of a white dwarf / cool star binary that emits from X-ray to radio wavelengths. The star, AR Scorpii (henceforth AR Sco), was classified in the early 1970s as a delta-Scuti star, a common variety of periodic variable star. Our observations reveal instead a 3.56 hr period close binary, pulsing in brightness on a period of 1.97 min. The pulses are so intense that AR Sco's optical flux can increase by a factor of four within 30 s, and they are detectable at radio frequencies, the first such detection for any white dwarf system. They reflect the spin of a magnetic white dwarf which we find to be slowing down on a 10^7 yr timescale. The spin-down power is an order of magnitude larger than that seen in electromagnetic radiation, which, together with an absence of obvious signs of accretion, suggests that AR Sco is primarily spin-powered. Although the pulsations are driven by the white dwarf's spin, they originate in large part from the cool star. AR Sco's broad-band spectrum is characteristic of synchrotron radiation, requiring relativistic electrons. These must either originate from near the white dwarf or be generated in situ at the M star through direct interaction with the white dwarf's magnetosphere.

Using the Binary Population and Spectral Synthesis code BPASS, we have calculated the rates, timescales and mass distributions for binary black hole mergers as a function of metallicity. We consider these in the context of the recently reported 1st LIGO event detection. We find that the event has a very low probability of arising from a stellar population with initial metallicity mass fraction above Z=0.010 (Z>0.5Zsun). Binary black hole merger events with the reported masses are most likely in populations below 0.008 (Z<0.4Zsun). Events of this kind can occur at all stellar population ages from ~3 Myr up to the age of the universe, but constitute only 0.1 to 0.4 per cent of binary BH mergers between metallicities of Z=0.001 to 0.008. However at metallicity Z=0.0001, 26 per cent of binary BH mergers would be expected to have the reported masses. At this metallicity the progenitor merger times can be close to ~10Gyr and rotationally-mixed stars evolving through quasi-homogeneous evolution, due to mass transfer in a binary, dominate the rate. The masses inferred for the black holes in the binary progenitor of GW\,150914 are amongst the most massive expected at anything but the lowest metallicities in our models. We discuss the implications of our analysis for the electromagnetic follow-up of future LIGO event detections.

Long-duration Gamma-Ray Bursts (GRBs) allow us to pinpoint and study star-forming galaxies in the early universe, thanks to their orders of magnitude brighter peak luminosities compared to other astrophysical sources, and their association with deaths of massive stars. We present Hubble Space Telescope Wide Field Camera 3 detections of three Swift GRB host galaxies lying at redshifts $z = 5.913$ (GRB 130606A), $z = 6.295$ (GRB 050904), and $z = 6.327$ (GRB 140515A) in the F140W (wide-$JH$ band, $\lambda_{\rm{obs}}\sim1.4\,\mu m$) filter. The hosts have magnitudes (corrected for Galactic extinction) of $m_{\rm{\lambda_{obs},AB}}= 26.34^{+0.14}_{-0.16}, 27.56^{+0.18}_{-0.22},$ and $28.30^{+0.25}_{-0.33}$ respectively. In all three cases the probability of chance coincidence of lower redshift galaxies is $\lesssim2\,\%$, indicating that the detected galaxies are most likely the GRB hosts. These are the first detections of high redshift ($z > 5$) GRB host galaxies in emission. The galaxies have luminosities in the range $0.1-0.6\,L^{*}_{z=6}$ (with $M_{1600}^{*}=-20.95\pm0.12$), and half-light radii in the range $0.6-0.9\,\rm{kpc}$. Both their half-light radii and luminosities are consistent with existing samples of Lyman-break galaxies at $z\sim6$. Spectroscopic analysis of the GRB afterglows indicate low metallicities ($[\rm{M/H}]\lesssim-1$) and low dust extinction ($A_{\rm{V}}\lesssim0.1$) along the line of sight. Using stellar population synthesis models, we explore the implications of each galaxy's luminosity for its possible star formation history, and consider the potential for emission-line metallicity determination with the upcoming James Webb Space Telescope.

We use the large cosmological hydro-dynamic simulation BlueTides to predict the photometric properties of galaxies during the epoch of reionisation ($z=8-15$). These properties include the rest-frame UV to near-IR broadband spectral energy distributions, the Lyman continuum photon production, the UV star formation rate calibration, and intrinsic UV continuum slope. In particular we focus on exploring the effect of various modelling assumptions, including the assumed choice of stellar population synthesis model, initial mass function, and the escape fraction of Lyman continuum photons, upon these quantities. We find that these modelling assumptions can have a dramatic effect on photometric properties leading to consequences for the accurate determination of physical properties from observations. For example, at $z=8$ we predict that nebular emission can account for up-to $50\%$ of the rest-frame $R$-band luminosity, while the choice of stellar population synthesis model can change the Lyman continuum production rate up to a factor of $\times 2$.

Intense, compact, star-forming galaxies are rare in the local Universe but ubiquitous at high redshift. We interpret the 0.1-22 um spectral energy distributions (SED) of a sample of 180 galaxies at 0.05<z<0.25 selected for extremely high surface densities of inferred star formation in the ultraviolet. By comparison with well-established stellar population synthesis models we find that our sample comprises young (~ 60 - 400 Myrs), moderate mass (~ $6\times 10^9$ Msun) star-forming galaxies with little dust extinction (mean stellar continuum extinction $E_\mathrm{cont}$(B-V) ~ 0.1) and find star formation rates of a few tens of Solar masses per year. We use our inferred masses to determine a mean specific star formation rate for this sample of ~ $10^{-9}$ yr$^{-1}$, and compare this to the specific star formation rates in distant Lyman break galaxies (LBGs), and in other low redshift populations. We conclude that our sample's characteristics overlap significantly with those of the z~5 LBG population, making ours the first local analogue population well tuned to match those high redshift galaxies. We consider implications for the origin and evolution of early galaxies.

The Lyman Continuum photon production efficiency ($\xi_{\rm ion}$) is a critical ingredient for inferring the number of photons available to reionise the intergalactic medium. To estimate the theoretical production efficiency in the high-redshift Universe we couple the BlueTides cosmological hydrodynamical simulation with a range of stellar population synthesis models. We find Lyman Continuum photon production efficiencies of $\log_{10}(\xi_{\rm ion}/{\rm erg^{-1}\, Hz})\approx 25.1-25.5$ depending on the choice of stellar population synthesis model. These results are broadly consistent with recent observational constraints at high-redshift though favour a model incorporating the effects of binary evolution

We present Hubble Space Telescope (HST) WFC3 UV and near-IR (nIR) imaging of 21 Superluminous Supernovae (SLSNe) host galaxies, providing a sensitive probe of star formation and stellar mass with the hosts. Comparing the photometric and morphological properties of these host galaxies with those of core collapse supernovae (CCSNe) and long-duration gamma-ray bursts (LGRBs), we find SLSN hosts are fainter and more compact at both UV and nIR wavelengths, in some cases we barely recover hosts with absolute magnitude around MV ~ -14. With the addition of ground based optical observations and archival results, we produce spectral energy distribution (SED) fits to these hosts, and show that SLSN hosts possess lower stellar mass and star formation rates. This is most pronounced for the hydrogen deficient Type-I SLSN hosts, although Type-II H-rich SLSN host galaxies remain distinct from the bulk of CCSNe, spanning a remarkably broad range of absolute magnitudes, with ~30% of SLSNe-II arising from galaxies fainter than Mn I R ~ -14. The detection of our faintest SLSN hosts increases the confidence that SLSNe-I hosts are distinct from those of LGRBs in star formation rate and stellar mass, and suggests that apparent similarities in metallicity may be due to the limited fraction of hosts for which emission line metallicity measurements are feasible. The broad range of luminosities of SLSN-II hosts is difficult to describe by metallicity cuts, and does not match the expectations of any reasonable UV-weighted luminosity function, suggesting additional environmental constraints are likely necessary to yield hydrogen rich SLSNe.

The relationship between stellar populations and the ionizing flux with which they irradiate their surroundings has profound implications for the evolution of the intergalactic medium. We quantify the ionizing flux arising from synthetic stellar populations which incorporate the evolution of interacting binary stars. We determine that these show ionizing flux boosted by 60 per cent at 0.05 < Z < 0.3 Z_sun and a more modest 10-20 per cent at near-Solar metallicities relative to star-forming populations in which stars evolve in isolation. The relation of ionizing flux to observables such as 1500A continuum and ultraviolet spectral slope is sensitive to attributes of the stellar population including age, star formation history and initial mass function. For a galaxy forming 1 M_sun yr^{-1}, observed at > 100 Myr after the onset of star formation, we predict a production rate of photons capable of ionizing hydrogen, N_ion = 1.4 x 10^{53} s^{-1} at Z = Z_sun and 3.5 x 10^{53} s^{-1} at 0.1 Z_sun, assuming a Salpeter-like initial mass function. We evaluate the impact of these issues on the ionization of the intergalactic medium, finding that the known galaxy populations can maintain the ionization state of the Universe back to z ~ 9, assuming that their luminosity functions continue to M_UV = -10, and that constraints on the intergalactic medium at z ~ 2 - 5 can be satisfied with modest Lyman continuum photon escape fractions of 4 - 24 per cent depending on assumed metallicity.

Some few hundred million years after the big bang the Universe was illuminated by the first stars and galaxies thereby bringing an end to the cosmological dark ages. Since the installation of WFC3 on the Hubble Space Telescope our ability to probe this critical period of the Universe's history has dramatically changed with thousands of objects now identified within the first billion years of the Universe's history. Our understanding of this period of the Universe's history will further grow thanks to both the Atacama Large millimetre/sub-millimetre Array and the James Webb Space Telescope.

We present late-time multi-wavelength observations of Swift J1644+57, suggested to be a relativistic tidal disruption flare (TDF). Our observations extend to >4 years from discovery, and show that 1.4 years after outburst the relativistic jet switched-off on a timescale less than tens of days, corresponding to a power-law decay faster than $t^{-70}$. Beyond this point weak X-rays continue to be detected at an approximately constant luminosity of $L_X \sim 5 \times 10^{42}$ erg s$^{-1}$, and are marginally inconsistent with a continuing decay of $t^{-5/3}$, similar to that seen prior to the switch-off. Host photometry enables us to infer a black hole mass of $M_{BH}=3 \times 10^6$ M$_{\odot}$, consistent with the late time X-ray luminosity arising from sub-Eddington accretion onto the black hole in the form of either an unusually optically faint AGN or a slowly varying phase of the transient. Optical/IR observations show a clear bump in the light curve at timescales of 30-50 days, with a peak magnitude (corrected for host galaxy extinction) of $M_R \sim -22-23$. The luminosity of the bump is significantly higher than seen in other, non-relativistic TDFs and does not match any re-brightening seen at X-ray or radio wavelengths. Its luminosity, light curve shape and spectrum are broadly similar to those seen in superluminous SNe, although subject to large uncertainties in the correction of the significant host extinction. We discuss these observations in the context of both TDF and massive star origins for Swift J1644+5734 and other candidate relativistic tidal flares.

High redshift quasars can be used to trace the early growth of massive galaxies and may be triggered by galaxy-galaxy interactions. We present MUSE science verification data on one such interacting system consisting of the well-studied z=3.2 PKS1614+051 quasar, its AGN companion galaxy and bridge of material radiating in Lyalpha between the quasar and its companion. We find a total of four companion galaxies (at least two galaxies are new discoveries), three of which reside within the likely virial radius of the quasar host, suggesting that the system will evolve into a massive elliptical galaxy by the present day. The MUSE data are of sufficient quality to split the extended Lyalpha emission line into narrow velocity channels. In these the gas can be seen extending towards each of the three neighbouring galaxies suggesting that the emission-line gas originates in a gravitational interaction between the galaxies and the quasar host. The photoionization source of this gas is less clear but is probably dominated by the two AGN. The quasar's Lyalpha emission spectrum is double-peaked, likely due to absorbing neutral material at the quasar's systemic redshift with a low column density as no damping wings are present. The spectral profiles of the AGN and bridge's Lyalpha emission are also consistent with absorption at the same redshift indicating this neutral material may extend over > 50 kpc. The fact that the neutral material is seen in the line of sight to the quasar and transverse to it, and the fact that we see the quasar and it also illuminates the emission-line bridge, suggests the quasar radiates isotropically and any obscuring torus is small. These results demonstrate the power of MUSE for investigating the dynamics of interacting systems at high redshift.

We present observations of Swift J1112.2-8238, and identify it as a candidate relativistic tidal disruption flare (rTDF). The outburst was first detected by Swift/BAT in June 2011 as an unknown, long-lived (order of days) $\gamma$-ray transient source. We show that its position is consistent with the nucleus of a faint galaxy for which we establish a likely redshift of $z=0.89$ based on a single emission line that we interpret as the blended [OII]$\lambda3727$ doublet. At this redshift, the peak X/$\gamma$-ray luminosity exceeded $10^{47}$ ergs s$^{-1}$, while a spatially coincident optical transient source had $i^{\prime} \sim 22$ (M$_g \sim -21.4$ at $z=0.89$) during early observations, $\sim 20$ days after the Swift trigger. These properties place Swift J1112.2-8238 in a very similar region of parameter space to the two previously identified members of this class, Swift J1644+57 and Swift J2058+0516. As with those events the high-energy emission shows evidence for variability over the first few days, while late time observations, almost 3 years post-outburst, demonstrate that it has now switched off. Swift J1112.2-8238 brings the total number of such events observed by Swift to three, interestingly all detected by Swift over a $\sim$3 month period ($<3\%$ of its total lifetime as of March 2015). While this suggests the possibility that further examples may be uncovered by detailed searches of the BAT archives, the lack of any prime candidates in the years since 2011 means these events are undoubtedly rare.

We have observed the host galaxy of the low redshift, low luminosity GRB 080517 at 105.8 GHz using the IRAM Plateau de Bure interferometer. We detect an emission line with integrated flux S.delta{nu} = 0.39 +/- 0.05 Jy km/s - consistent both spatially and in velocity with identification as the J=1-0 rotational transition of carbon monoxide (CO) at the host galaxy redshift. This represents only the third long GRB host galaxy with molecular gas detected in emission. The inferred molecular gas mass, M_H2 ~ 6.3 x 10^8 M_sun, implies a gas consumption timescale of ~40 Myr if star formation continues at its current rate. Similar short timescales appear characteristic of the long GRB population with CO observations to date, suggesting that the gamma-ray burst in these sources occurs towards the end of their star formation episode.

We present an analysis of the photometry and spectroscopy of the host galaxy of Swift-detected GRB 080517. From our optical spectroscopy, we identify a redshift of z = 0.089 +/- 0.003, based on strong emission lines, making this a rare example of a very local, low luminosity, long gamma ray burst. The galaxy is detected in the radio with a flux density of S(4.8GHz) =0.22 +/- 0.04mJy - one of relatively few known GRB hosts with a securely measured radio flux. Both optical emission lines and a strong detection at 22 um suggest that the host galaxy is forming stars rapidly, with an inferred star formation rate ~16 Msun/yr and a high dust obscuration (E(B-V )>1, based on sight-lines to the nebular emission regions). The presence of a companion galaxy within a projected distance of 25 kpc, and almost identical in redshift, suggests that star formation may have been triggered by galaxy-galaxy interaction. However, fitting of the remarkably flat spectral energy distribution from the ultraviolet through to the infrared suggests that an older, 500Myr post-starburst stellar population is present along with the ongoing star formation. We suggest that that the host galaxy of GRB 080517 is a valuable addition to the still very small sample of well-studied local gamma-ray burst hosts.

Star forming galaxies at high redshift show ubiquitously high ionization parameters, as measured by the ratio of optical emission lines. We demonstrate that local (z < 0.2) sources selected as Lyman break analogues also manifest high line ratios with a typical [O III]/H beta = 3.36(+0.14,-0.04) - comparable to all but the highest ratios seen in star forming galaxies at z ~ 2-4. We argue that the stellar population synthesis code BPASS can explain the high ionization parameters required through the ageing of rapidly formed star populations, without invoking any AGN contribution. Binary stellar evolution pathways prolong the age interval over which a starburst is likely to show elevated line ratios, relative to those predicted by single stellar evolution codes. As a result, model galaxies at near-Solar metallicities and with ages of up to ~100 Myr after a starburst typically have a line ratio [O III]/H beta~3, consistent with those seen in Lyman break galaxies and local sources with similar star formation densities. This emphasises the importance of including binary evolution pathways when simulating the nebular line emission of young or bursty stellar populations.

We present rest-frame far-infrared (FIR) and optical observations of the host galaxy of GRB090423 at z=8.23 from the Atacama Large Millimeter Array (ALMA) and the Spitzer Space Telescope, respectively. The host remains undetected to 3-sigma limits of Fnu(222 GHz)<33 microJy and Fnu(3.6 micron)<81 nJy. The FIR limit is about 20 times fainter than the luminosity of the local ULIRG Arp220, and comparable to the local starburst M82. Comparing to model spectral energy distributions we place a limit on the IR luminosity of L_IR(8-1000 micron)<3e10 Lsun, corresponding to a limit on the obscured star formation rate of SFR_IR<5 Msun/yr; for comparison, the limit on the unobscured star formation rate from Hubble Space Telescope rest-frame UV observations is SFR_UV<1 Msun/yr. We also place a limit on the host galaxy stellar mass of <5e7 Msun (for a stellar population age of 100 Myr and constant star formation rate). Finally, we compare our millimeter observations to those of field galaxies at z>4 (Lyman break galaxies, Ly-alpha emitters, and submillimeter galaxies), and find that our limit on the FIR luminosity is the most constraining to date, although the field galaxies have much larger rest-frame UV/optical luminosities than the host of GRB090423 by virtue of their selection techniques. We conclude that GRB host galaxies at z>4, especially those with measured interstellar medium metallicities from afterglow spectroscopy, are an attractive sample for future ALMA studies of high redshift obscured star formation.

We present 5.5 and 9.0 GHz observations of a sample of seventeen GRB host galaxies at 0.5<z<1.4, using the radio continuum to explore their star formation properties in the context of the small but growing sample of galaxies with similar observations. Four sources are detected, one of those (GRB 100418A) likely due to lingering afterglow emission. We suggest that the previously-reported radio afterglow of GRB 100621A may instead be due to host galaxy flux. We see no strong evidence for redshift evolution in the typical star formation rate of GRB hosts, but note that the fraction of `dark' bursts with detections is higher than would be expected given constraints on the more typical long GRB population. We also determine the average radio-derived star formation rates of core collapse supernovae at comparable redshift, and show that these are still well below the limits obtained for GRB hosts, and show evidence for a rise in typical star formation rate with redshift in supernova hosts.

Following our previous spectroscopic observations of $z>7$ galaxies with Gemini/GNIRS and VLT/XSHOOTER, which targeted a total of 8 objects, we present here our results from a deeper and larger VLT/FORS2 spectroscopic sample of Wide Field Camera 3 selected $z>7$ candidate galaxies. With our FORS2 setup we cover the 737-1070nm wavelength range, enabling a search for Lyman-$\alpha$ in the redshift range spanning 5.06 - 7.80. We target 22 $z$-band dropouts and find no evidence of Lyman-$\alpha$ emission, with the exception of a tentative detection ($<5\sigma$, which is our adopted criterion for a secure detection) for one object. The upper limits on Lyman-$\alpha$ flux and the broad-band magnitudes are used to constrain the rest-frame Equivalent Widths for this line emission. We analyse our FORS2 observations in combination with our previous GNIRS and XSHOOTER observations, and suggest that a simple model where the fraction of high rest-frame Equivalent Width emitters follows the trend seen at $z=3-6.5$ is inconsistent with our non-detections at $z\sim7.8$ at the 96\% confidence level. This may indicate that a significant neutral HI fraction in the intergalactic medium suppresses Lyman-$\alpha$, with an estimated neutral fraction $\chi_{HI}\sim0.5$, in agreement with other estimates.

The sensitivity available to near-infrared surveys has recently allowed us to probe the galaxy population at $z\approx 7$ and beyond. The existing {\em Hubble} Wide Field Camera 3 (WFC3) and the Visible and Infrared Survey Telescope for Astronomy (VISTA) Infrared Camera (VIRCam) instruments allow deep surveys to be undertaken well beyond one micron - a capability that will be further extended with the launch and commissioning of the {\em James Webb Space Telescope (JWST)}. As new regions of parameter space in both colour and depth are probed new challenges for distant galaxy surveys are identified. In this paper we present an analysis of the colours of L and T dwarf stars in widely used photometric systems. We also consider the implications of the newly-identified Y dwarf population - stars that are still cooler and less massive than T dwarfs for both the photometric selection and spectroscopic follow-up of faint and distant galaxies. We highlight the dangers of working in the low-signal-to-noise regime, and the potential contamination of existing and future samples. We find that {\em Hubble}/WFC3 and VISTA/VIRCam $Y$-drop selections targeting galaxies at $z\sim7.5$ are vulnerable to contamination from T and Y class stars. Future observations using {\em JWST}, targeting the $z\sim7$ galaxy population, are also likely to prove difficult without deep medium-band observations. We demonstrate that single emission line detections in typical low signal-to-noise spectroscopic observations may also be suspect, due to the unusual spectral characteristics of the cool dwarf star population.

Lyman break analogues (LBAs) are local galaxies selected to match a more distant (usually z~3) galaxy population in luminosity, UV-spectral slope and physical characteristics, and so provide an accessible laboratory for exploring their properties. However, as the Lyman break technique is extended to higher redshifts, it has become clear that the Lyman break galaxies (LBGs) at z~3 are more massive, luminous, redder, more extended and at higher metallicities than their z~5 counterparts. Thus extrapolations from the existing LBA samples (which match z=3 properties) have limited value for characterising z>5 galaxies, or inferring properties unobservable at high redshift. We present a new pilot sample of twenty-one compact star forming galaxies in the local (0.05<z<0.25) Universe, which are tuned to match the luminosities and star formation volume densities observed in z>~5 LBGs. Analysis of optical emission line indices suggests that these sources have typical metallicities of a few tenths Solar (again, consistent with the distant population). We also present radio continuum observations of a subset of this sample (13 sources) and determine that their radio fluxes are consistent with those inferred from the ultraviolet, precluding the presence of a heavily obscured AGN or significant dusty star formation.

Identifying galaxy clustering at high redshift (i.e. z > 1) is essential to our understanding of the current cosmological model. However, at increasing redshift, clusters evolve considerably in star-formation activity and so are less likely to be identified using the widely-used red sequence method. Here we assess the viability of instead identifying high redshift clustering using actively star-forming galaxies (SMGs associated with over-densities of BzKs/LBGs). We perform both a 2- and 3-D clustering analysis to determine whether or not true (3D) clustering can be identified where only 2D data are available. As expected, we find that 2D clustering signals are weak at best and inferred results are method dependant. In our 3D analysis, we identify 12 SMGs associated with an over-density of galaxies coincident both spatially and in redshift - just 8% of SMGs with known redshifts in our sample. Where an SMG in our target fields lacks a known redshift, their sightline is no more likely to display clustering than blank sky fields; prior redshift information for the SMG is required to identify a true clustering signal. We find that the strength of clustering in the volume around typical SMGs, while identifiable, is not exceptional. However, we identify a small number of highly clustered regions, all associated with an SMG. The most notable of these, surrounding LESSJ033336.8-274401, potentially contains an SMG, a QSO and 36 star-forming galaxies (a > 20sig over-density) all at z~1.8. This region is highly likely to represent an actively star-forming cluster and illustrates the success of using star-forming galaxies to select sites of early clustering. Given the increasing number of deep fields with large volumes of spectroscopy, or high quality and reliable photometric redshifts, this opens a new avenue for cluster identification in the young Universe.

ABRIDGED: We have used the LABOCA Survey of the ECDF-S (LESS) to investigate rest-frame FIR emission from typical SF systems (LBGs) at redshift 3, 4, and 5. We initially concentrate on LBGs at z~3 and select three subsamples on stellar mass, extinction corrected SF and rest-frame UV-magnitude. We produce composite 870micron images of the typical source in our subsamples, obtaining ~4sigma detections and suggesting a correlation between FIR luminosity and stellar mass. We apply a similar procedure to our full samples at z~3, 4, 4.5 and 5 and do not obtain detections - consistent with a simple scaling between FIR luminosity and stellar mass. In order to constrain the FIR SED of these systems we explore their emission at multiple wavelengths spanning the peak of dust emission at z~3 using the Herschel SPIRE observations of the field. We obtain detections at multiple wavelengths for both our stellar mass and UV-magnitude selected samples, and find a best-fit SED with T_dust in the ~33-41K range. We calculate L_FIR, obscured SFRs and M_dust, and find that a significant fraction of SF in these systems is obscured. Interestingly, our extinction corrected SFR sample does not display the large FIR fluxes predicted from its red UV-spectral slope. This suggests that the method of assuming an intrinsic UV-slope and correcting for dust attenuation may be invalid for this sample - and that these are not in fact the most actively SF systems. All of our z~3 samples fall on the `main sequence' of SF galaxies at z~3 and our detected subsamples are likely to represent the high obscuration end of LBGs at their epoch. We compare the FIR properties of our subsamples with various other populations, finding that our stellar mass selected sample shows similar FIR characteristics to SMGs at the same epoch and therefore potentially represents the low L_FIR end of the high redshift FIR luminosity function.

Luminous high-redshift QSOs are thought to exist within the most massive dark matter haloes in the young Universe. As a consequence they are likely to be markers for biased, over-dense regions where early galaxies cluster, regions that eventually grow into the groups and clusters seen in the lower redshift universe. In this paper we explore the clustering of galaxies around z ~ 5 QSOs as traced by Lyman Break Galaxies (LBGs). We target the fields of three QSOs using the same optical imaging and spectroscopy techniques used in the ESO Remote Galaxy Survey (ERGS, Douglas et al. 2009, 2010), which was successful in identifying individual clustered structures of LBGs. We use the statistics of the redshift clustering in ERGS to show that two of the three fields show significant clustering of LBGs at the QSO redshifts. Neither of these fields is obviously over-dense in LBGs from the imaging alone; a possible reason why previous imaging-only studies of high-redshift QSO environments have given ambiguous results. This result shows that luminous QSOs at z ~ 5 are typically found in over-dense regions. The richest QSO field contains at least nine spectroscopically confirmed objects at the same redshift including the QSO itself, seven LBGs and a second fainter QSO. While this is a very strong observational signal of clustering at z ~ 5, it is of similar strength to that seen in two structures identified in the 'blank sky' ERGS fields. This indicates that, while over-dense, the QSO environments are not more extreme than other structures that can be identified at these redshifts. The three richest structures discovered in this work and in ERGS have properties consistent with that expected for proto-clusters and likely represent the early stages in the build-up of massive current-day groups and clusters.

We present spectroscopic observations with VLT/XSHOOTER and Subaru/MOIRCS of a relatively bright Y-band drop-out galaxy in the Hubble Ultra Deep Field, first selected by Bunker et al. (2010), McLure et al. (2010) and Bouwens et al. (2010) to be a likely z~8-9 galaxy on the basis of its colours in the HST ACS and WFC3 images. This galaxy, HUDF.YD3 (also known as UDFy-38135539) has been targetted for VLT/SINFONI integral field spectroscopy by Lehnert et al. (2010), who published a candidate Lyman-alpha emission line at z=8.55 from this source. In our independent spectroscopy using two different infrared spectrographs (5 hours with VLT/XSHOOTER and 11 hours with Subaru/MOIRCS) we are unable to reproduce this line. We do not detect any emission line at the spectral and spatial location reported in Lehnert et al. (2010), despite the expected signal in our combined MOIRCS & XSHOOTER data being 5-sigma. The line emission also seems to be ruled out by the faintness of this object in recently extremely deep F105W (Y-band) HST/WFC3 imaging from HUDF12; the line would fall within this filter and such a galaxy should have been detected at Y(AB)=28.6mag (~20 sigma) rather than the marginal Y(AB)~30mag observed in the Y-band image, >3 times fainter than would be expected if the emission lie was real. Hence it appears highly unlikely that the reported Lyman-alpha line emission at z>8 is real, meaning that the highest-redshift sources for which Lyman-alpha emission has been seen are at z=6.9-7.2. It is conceivable that Lyman-alpha does not escape galaxies at higher redshifts, where the Gunn-Peterson absorption renders the Universe optically thick to this line. However, deeper spectroscopy on a larger sample of candidate z>7 galaxies will be needed to test this.

(Abridged) Distant galaxy clusters provide important tests of the growth of large scale structure in addition to highlighting the process of galaxy evolution in a consistently defined environment at large look back time. We present a sample of 22 distant (z>0.8) galaxy clusters and cluster candidates selected from the 9 deg2 footprint of the overlapping X-ray Multi Mirror (XMM) Large Scale Structure (LSS), CFHTLS Wide and Spitzer SWIRE surveys. Clusters are selected as extended X-ray sources with an accompanying overdensity of galaxies displaying optical to mid-infrared photometry consistent with z>0.8. Nine clusters have confirmed spectroscopic redshifts in the interval 0.8<z<1.2, four of which are presented here for the first time. A further 11 candidate clusters have between 8 and 10 band photometric redshifts in the interval 0.8<z<2.2, while the remaining two candidates do not have information in sufficient wavebands to generate a reliable photometric redshift. All of the candidate clusters reported in this paper are presented for the first time. Those confirmed and candidate clusters with available near infrared photometry display evidence for a red sequence galaxy population, determined either individually or via a stacking analysis, whose colour is consistent with the expectation of an old, coeval stellar population observed at the cluster redshift. We further note that the sample displays a large range of red fraction values indicating that the clusters may be at different stages of red sequence assembly. We compare the observed X-ray emission to the flux expected from a suite of model clusters and find that the sample displays an effective mass limit M200 ~ 1e14 Msolar with all clusters displaying masses consistent with M200 < 5e14 Msolar. This XMM distant cluster study represents a complete sample of X-ray selected z>0.8 clusters.

The recent Hubble Space Telescope near-infrared imaging with the Wide-Field Camera #3 (WFC3) of the GOODS-South field in the CANDELS program covering nearly 100arcmin^2, along with already existing Advanced Camera for Surveys optical data, makes possible the search for bright galaxy candidates at redshift z ~ 7 - 9 using the Lyman-break technique. We present the first analysis of z'-drop z ~ 7 candidate galaxies in this area, finding 19 objects. We also analyse Y-drops at z ~ 8, trebling the number of bright (H_AB < 27 mag) Y-drops from our previous work, and compare our results with those of other groups based on the same data. The bright high redshift galaxy candidates we find serve to better constrain the bright end of the luminosity function at those redshift, and may also be more amenable to spectroscopic confirmation than the fainter ones presented in various previous work on the smaller fields (the Hubble Ultra Deep Field and the WFC 3 Early Release Science observations). We also look at the agreement with previous luminosity functions derived from WFC3 drop-out counts, finding a generally good agreement, except for the luminosity function of Yan et al. (2010) at z ~ 8, which is strongly ruled out.

We present Gemini/GNIRS spectroscopic observations of 4 z-band (z~7) dropout galaxies and VLT/XSHOOTER observations of one z-band dropout and 3 Y-band (z~8-9) dropout galaxies in the Hubble Ultra Deep Field, which were selected with Wide Field Camera 3 imaging on the Hubble Space Telescope. We find no evidence of Lyman-alpha emission with a typical 5-sigma sensitivity of 5X10^-18erg/cm^2/s, and we use the upper limits on Lyman-alpha flux and the broad-band magnitudes to constrain the rest-frame equivalent widths for this line emission. Accounting for incomplete spectral coverage, we survey 3.0 z-band dropouts and 2.9 Y-band dropouts to a Lyman-alpha rest-frame equivalent width limit > 120Ang (for an unresolved emission line); for an equivalent width limit of 50Ang the effective numbers of drop-outs surveyed fall to 1.2 z-band drop-outs and 1.5 Y-band drop-outs. A simple model where the fraction of high rest-frame equivalent width emitters follows the trend seen at z=3-6.5 is inconsistent with our non-detections at z=7-9 at the ~ 1-sigma level for spectrally unresolved lines, which may indicate that a significant neutral HI fraction in the intergalactic medium suppresses the Lyman-alpha line in z-drop and Y-drop galaxies at z > 7.

We present 1.2mm MAMBO-2 observations of a field which is over-dense in Lyman Break Galaxies (LBGs) at z~5. The field includes seven spectroscopically-confirmed LBGs contained within a narrow (z=4.95+/-0.08) redshift range and an eighth at z=5.2. We do not detect any individual source to a limit of 1.6 mJy/beam (2*rms). When stacking the flux from the positions of all eight galaxies, we obtain a limit to the average 1.2 mm flux of these sources of 0.6mJy/beam. This limit is consistent with FIR imaging in other fields which are over-dense in UV-bright galaxies at z~5. Independently and combined, these limits constrain the FIR luminosity (8-1000 micron) to a typical z~5 LBG of LFIR<~3x10^11 Lsun, implying a dust mass of Mdust<~10^8 Msun (both assuming a grey body at 30K). This LFIR limit is an order of magnitude fainter than the LFIR of lower redshift sub-mm sources (z~1-3). We see no emission from any other sources within the field at the above level. While this is not unexpected given millimetre source counts, the clustered LBGs trace significantly over-dense large scale structure in the field at z = 4.95. The lack of any such detection in either this or the previous work, implies that massive, obscured star-forming galaxies may not always trace the same structures as over-densities of LBGs, at least on the length scale probed here. We briefly discuss the implications of these results for future observations with ALMA.

We present late-time Hubble Space Telescope imaging of the fields of six Swift GRBs lying at 5.0<z<9.5. Our data includes very deep observations of the field of the most distant spectroscopically confirmed burst, GRB 090423, at z=8.2. Using the precise positions afforded by their afterglows we can place stringent limits on the luminosities of their host galaxies. In one case, that of GRB 060522 at z=5.11, there is a marginal excess of flux close to the GRB position which may be a detection of a host at a magnitude J(AB)=28.5. None of the others are significantly detected meaning that all the hosts lie below L\star at their respective redshifts, with star formation rates SFR<4Mo/yr in all cases. Indeed, stacking the five fields with WFC3-IR data we conclude a mean SFR<0.17Mo/yr per galaxy. These results support the proposition that the bulk of star formation, and hence integrated UV luminosity, at high redshifts arises in galaxies below the detection limits of deep-field observations. Making the reasonable assumption that GRB rate is proportional to UV luminosity at early times allows us to compare our limits with expectations based on galaxy luminosity functions derived from the Hubble Ultra-Deep Field (HUDF) and other deep fields. We infer that a luminosity function which is evolving rapidly towards steeper faint-end slope (alpha) and decreasing characteristic luminosity (L\star), as suggested by some other studies, is consistent with our observations, whereas a non-evolving LF shape is ruled out at >90% confidence. Although it is not yet possible to make stronger statements, in the future, with larger samples and a fuller understanding of the conditions required for GRB production, studies like this hold great potential for probing the nature of star formation, the shape of the galaxy luminosity function, and the supply of ionizing photons in the early universe.

Young, massive stars dominate the rest-frame ultraviolet spectra of star-forming galaxies. At high redshifts (z>2), these rest-UV features are shifted into the observed-frame optical and a combination of gravitational lensing, deep spectroscopy and spectral stacking analysis allows the stellar population characteristics of these sources to be investigated. We use our stellar population synthesis code BPASS to fit two strong rest-UV spectral features in published Lyman-break galaxy spectra, taking into account the effects of binary evolution on the stellar spectrum. In particular, we consider the effects of quasi-homogeneous evolution (arising from the rotational mixing of rapidly-rotating stars), metallicity and the relative abundance of carbon and oxygen on the observed strengths of HeII (1640 Angstroms) and CIV (1548,1551 Angstroms) spectral lines. We find that Lyman-break galaxy spectra at z=2-3 are best fit with moderately sub-solar metallicities, and with a depleted carbon-to-oxygen ratio. We also find that the spectra of the lowest metallicity sources are best fit with model spectra in which the HeII emission line is boosted by the inclusion of the effect of massive stars being spun-up during binary mass-transfer so these rapidly-rotating stars experiencing quasi-homogeneous evolution.

The acquisition of deep Near-IR imaging with Wide Field Camera 3 on the Hubble Space Telescope has provided the opportunity to study the very-high redshift Universe. For galaxies up to $z\approx 7.7$ sufficient wavelength coverage exists to probe the rest-frame ultraviolet (UV) continuum without contamination from either Lyman-$\alpha$ emission or the Lyman-$\alpha$ break. In this work we use Near-IR imaging to measure the rest-frame UV continuum colours of galaxies at $4.7<z<7.7$. We are have carefully defined a colour-colour selection to minimise any inherent bias in the measured UV continuum slope for the drop-out samples. For the highest-redshift sample ($6.7<z<7.7$), selected as $z_{f850lp}$-band dropouts, we find mean UV continuum colours approximately equal to zero (AB), consistent with a dust-free, solar metallicity, star forming population (or a moderately dusty population of low metallicity). At lower-redshift we find that the mean UV continuum colours of galaxies (over the same luminosity range) are redder, and that galaxies with higher luminosities are also slightly redder on average. One interpretation of this is that lower-redshift and more luminous galaxies are dustier, however this interpretation is complicated by the effects of the star formation history and metallicity and potentially the initial mass function on the UV continuum colours.

[Abridged] Observations of molecular gas at all redshifts are critical for measuring the cosmic evolution in molecular gas density and understanding the star-formation history of the Universe. The 12CO molecule (J=1-0 transition = 115.27 GHz) is the best proxy for extragalactic H2, which is the gas reservoir from which star formation occurs, and has been detected out to z~6. Typically, redshifted high-J lines are observed at mm-wavelengths, the most commonly targeted systems exhibiting high SFRs (e.g. submm galaxies), and far-IR-bright QSOs. While the most luminous objects are the most readily observed, detections of more typical galaxies with modest SFRs are essential for completing the picture. ALMA will be revolutionary in terms of increasing the detection rate and pushing the sensitivity limit down to include such galaxies, however the limited FoV when observing at such high frequencies makes it difficult to use ALMA for studies of the large-scale structure traced out by molecular gas in galaxies. This article introduces a strategy for a systematic search for molecular gas during the EoR (z~7 and above), capitalizing on the fact that the J=1-0 transition of 12CO enters the upper bands of cm-wave instruments at high-z. The FoV advantage gained by observing at such frequencies, coupled with modern broadband correlators allows significant cosmological volumes to be probed on reasonable timescales. In this article we present an overview of our future observing programme which has been awarded 6,500 hours as one of the Large Survey Projects for MeerKAT, the forthcoming South African SKA pathfinder instrument. Its large FoV and correlator bandwidth, and high-sensitivity provide unprecedented survey speed for such work. An existing astrophysical simulation is coupled with instrumental considerations to demonstrate the feasibility of such observations and predict detection rates.

We present a search for galaxies at 7.6<z<9.8 using the latest HST WFC3 near-infrared data, based on the Lyman-break technique. We search for galaxies which have large (Y-J) colours (the "Y-drops") on account of the Lyman-alpha forest absorption, and with (J-H) colours inconsistent with being low-redshift contaminants. We identify 24 candidates at redshift z~8-9 (15 are robust and a further 9 more marginal but consistent with being high redshift) over an area of ~50 square arcminutes. Previous searches for Y-drops with WFC3 have focussed only on the Hubble Ultra Deep Field (HUDF), and our larger survey (involving two other nearby deep fields and a wider area survey) has trebelled the number of robust Y-drop candidates. For the first time, we have sufficient z~8-9 galaxies to fit both phi^* and M^* of the UV Schechter luminosity function. There is evidence for evolution in this luminosity function from z=6-7 to z=8-9, in the sense that there are fewer UV-bright galaxies at z~8-9, consistent with an evolution mainly in M^*. The candidate z~8-9 galaxies we detect have insufficient ionizing flux to reionize the Universe, and it is probable that galaxies below our detection limit provide a significant UV contribution. The faint-end slope, alpha, is not well constrained. However, adopting a similiar faint-end slope to that determined at z=3-6 (alpha=-1.7) and a Salpeter initial mass function, then the ionizing photon budget still falls short if f_escape<0.5, even integrating down to M(UV)=-8. A steeper faint end slope or a low-metallicity population (or a top-heavy IMF) might still provide sufficient photons for star-forming galaxies to reionize the Universe, but confirmation of this might have to await the James Webb Space Telescope.

GRB 090423 is the most distant spectroscopically-confirmed source observed in the universe. Using observations at 37.5 GHz, we place constraints on molecular gas emission in the CO(3-2) line from its host galaxy and immediate environs. The source was not detected either in line emission or in the rest-frame 850 micron continuum, yielding an upper limit of S_{8mm}=9.3 milli-Jy and M(H_2)<4.3x10^9 M_sun (3 sigma), applying standard conversions. This implies that the host galaxy of GRB 090423 did not possess a large reservoir of warm molecular gas but was rather modest either in star formation rate or in mass. It suggests that this was not an extreme starburst, and hence that gamma ray bursts at high redshift trace relatively modest star formation rates, in keeping with the behaviour seen at lower redshifts. We do, however, identify a millimetre emission line source in the field of GRB 090423. Plausible interpretations include a CO(1-0) emitting galaxy at z=2.1, CO(2-1) at z=5.2 and CO(3-2) at z=8.3. Efforts to identify a counterpart for the molecular line emitter and to further characterise this source are continuing.

We present limits on the molecular gas content of Lyman Break Galaxies (LBGs) at z~5 from observations targetting redshifted CO(1-0) and CO(2-1) line emission. We observed a single field containing eight spectroscopically-confirmed z~5 LBGs, seven of which are contained within a narrow (z=4.95\pm 0.08) redshift range and the eighth is at z=5.2. No source was individually detected. Assuming the CO to H2 conversion factor for vigorous starbursts, we place upper limits on the molecular gas content of individual z~5 LBGs of M(H2)<~10^10 Msun. From a stacking analysis combining all of the non-detections, the typical z~5 LBG has an H2 mass limit comparable to their stellar mass, < 3.1 x 10^9 Msun. This limit implies that, given the star formation rates of these systems (measured from their UV emission), star formation could be sustained for at most ~100Myr, similar to the typical ages of their stellar populations. The lack of a substantially larger reservoir of cold gas argues against the LBGs being UV luminous super starbursts embedded in much larger UV-dark systems and as a result increases the likelihood that at least those LBGs with multiple components are starbursts triggered by mergers. The sources responsible for reionization are expected to be starbursts similar to these systems, but with lower luminosities, masses and consequently with star formation timescales far shorter than the recombination timescale. If so, the ionized bubbles expected in the IGM during the reionization era will infrequently have UV-luminous sources at their centres.

(ABRIDGED) We present the global results of a large spectroscopic survey carried out in order to identify z~ 5 Lyman break galaxies (LBGs) across ten widely-separated fields to I_{AB}=26.3. The redshifts of seventy 4.6<z<5.6 LBGs were identified. We find no significant difference in the frequency of high equivalent-width line emitters between z~3 and our z~5 samples. The rest-frame UV continuum slope of a typical z~5 line-emitting galaxy is bluer than that of a typical break-only galaxy, a difference that is difficult to explain purely by differences in the ages of their stellar populations. Variation in metallicity and/or dust extinction can more straightforwardly account for this difference. Given their UV-continuum slopes, the typical z~5 LBGs have metallicities a factor of three lower than those of LBGs at z~3. HST imaging indicates that a large majority of the spectroscopically-confirmed LBGs in our sample are members of multiple systems and/or show disturbed morphology. Using local LBG analogues as a model, this multiplicity could be explained either by super-starburst regions within a larger unseen structure, or by a high incidence of merging events at this epoch. The current data cannot distinguish between these two possibilities. The surface density of z~5 LBGs in two fields is considerably higher than in the rest. Both show clear spikes in their redshift distributions indicating strong three-dimensional clustering. Neither structure can be bound given their depth in redshift and probably extend beyond the observed fields. The three-dimensional distances between LBGs in the structures are too large for them to have triggered their starbursts through mutual gravitational interaction, and so it is likely that the short-lived LBGs represent only a small fraction of the baryons in the structures.

We have carried out 870 micron observations in the J1040.7-1155 field, known to host an overdensity of Lyman break galaxies at z=5.16 +/- 0.05. We do not detect any individual source at the S(870)=3.0 mJy/beam (2 sigma) level. A stack of nine spectroscopically confirmed z>5 galaxies also yields a non-detection, constraining the submillimeter flux from a typical galaxy at this redshift to S(870)<0.85 mJy, which corresponds to a mass limit M(dust)<1.2x10^8 M_sun (2 sigma). This constrains the mass of thermal dust in distant Lyman break galaxies to less than one tenth of their typical stellar mass. We see no evidence for strong submillimeter galaxies associated with the ultraviolet-selected galaxy overdensity, but cannot rule out the presence of fainter, less massive sources.

We analyse recently acquired near-infrared Hubble Space Telescope imaging of the GOODS-South field to search for star forming galaxies at z~7.0. By comparing WFC 3 0.98 micron Y-band images with ACS z-band (0.85 micron) images, we identify objects with colours consistent with Lyman break galaxies at z~6.4-7.4. This new data covers an area five times larger than that previously reported in the WFC3 imaging of the Hubble Ultra Deep Field, and affords a valuable constraint on the bright end of the luminosity function. Using additional imaging of the region in the ACS B, V and i-bands from GOODS v2.0 and the WFC3 J-band we attempt to remove any low-redshift interlopers. Our selection criteria yields 6 candidates brighter than Y_AB = 27.0, of which all except one are detected in the ACS z-band imaging and are thus unlikely to be transients. Assuming all 6 candidates are at z~7 this implies a surface density of objects brighter than Y_AB = 27.0 of 0.30\pm0.12 arcmin-2, a value significantly smaller than the prediction from the z ~ 6 luminosity function. This suggests continued evolution of the bright end of the luminosity function between z = 6 to 7, with number densities lower at higher redshift.

The quest to discover the most distant galaxies has developed rapidly in the last decade. We are now exploring redshifts of 6 and beyond, when the Universe was less than a billion years old, an epoch when the previously-neutral intergalactic medium was reionized. The continuing discovery of galaxies at progressively higher and higher redshifts has been driven by the availability of large telescopes on the ground and in space, improvements in detector technology, and new search strategies. Over the past 4 years, the Lyman break technique has been shown to be effective in isolating z~6 star-forming i'-drop galaxies through spectroscopic confirmation with large ground-based telescopes (Keck, Gemini and the ESO VLTs). Narrow-band imaging, notably with the wide field of the Subaru telescope, has also produced samples of Lyman-alpha emitters at these redshifts. A The discovery of this i'-drop galaxy population has been used to infer the global star formation rate density at this epoch (z~6), and we are now beginning to constrain the contribution to reionization of the UV flux from these galaxies. Infrared data from the Spitzer Space Telescope has been used to determine the spectral energy distributions (SEDs) from the rest-frame UV to the optical, and constrain the previous star formation histories, masses and ages. The indications are that much of the stellar mass of these galaxies might have formed in vigorous bursts at z>6. The next big advances would be to test the population synthesis modelling of these z~6 galaxies through spectroscopy of the rest-frame optical (rather than crude broad-band SEDs), and also to push the observational horizon for galaxies further to directly explore star formation during the reionization epoch. JWST is likely to have a profound impact on realising these goals.

We describe the selection of a sample of photometrically-defined Lyman break galaxies (LBGs) at z~5 using the multi-wavelength imaging data of the ESO Remote Galaxy Survey (ERGS). The data is drawn from ten widely-separated fields covering a total sky area of 275 arcmin squared. Starting with a simple colour (R-I>1.3) and magnitude (I<26.3) cut to isolate the Lyman break and then refining the sample by applying further optical and near-infrared photometric criteria we identify a sample of 253 LBG candidates. We carefully model the completeness of this sample and the factors that affect its reliability. There is considerable overlap between this sample and a spectroscopically-confirmed sample drawn from the same survey and this allows us to determine the reliability of the optical photometric selection (~60 per cent) and to show that the reliability can be significantly improved (to ~80 per cent) by applying near-infrared waveband criteria to exclude very red contaminants. Even this high level of reliability may compromise some statistical studies of LBG properties. We show that over 30 per cent of the highest reliability candidates have multiple UV-luminous components and/or disturbed morphology in HST imaging, though it is unclear whether this represents multiple interacting/merging sources or individual large sources with multiple UV bright regions. Using this sample we confirm that the normalisation of the bright end of the z=5 UV luminosity function (down to M*) is lower than the same at z=4 by a factor of 3. Using a Schechter fit we determine M*uv=-20.9+/-0.2. We discuss whether it is reasonable to expect the UV luminosity function to follow a Schechter function, given the UV emission is short-lived and stochastic, and does not necessarily trace the underlying mass of the galaxy.

We have constructed a new code to produce synthetic spectra of stellar populations that includes massive binaries. We have tested this code against the broadband colours of unresolved young massive stellar clusters in nearby galaxies, the equivalent widths of the Red and Blue Wolf-Rayet bumps in star-forming SDSS galaxies and the UV and optical spectra of the star forming regions Tol-A and B in NGC5398. In each case we find a good agreement between our models and observations. We find that in general binary populations are bluer and have fewer red supergiants, and thus significantly less flux in the I-band and at longer wavelengths, than single star populations. Also we find that Wolf-Rayet stars occur over a wider range of ages up to 10^7 years in a stellar population including binaries, increasing the UV flux and Wolf-Rayet spectral features at later times. In addition we find that nebula emission contributes significantly to these observed properties and must be considered when comparing stellar models with observations of unresolved stellar populations. We conclude that incorporation of massive stellar binaries can improve the agreement between observations and synthetic spectral synthesis codes, particularly for systems with young stellar populations.

We present observations of redshifted CO(1-0) and CO(2-1) in a field containing an overdensity of Lyman break galaxies (LBGs) at z=5.12. Our Australia Telescope Compact Array observations were centered between two spectroscopically-confirmed z=5.12 galaxies. We place upper limits on the molecular gas masses in these two galaxies of M(H_2) <1.7 x 10^10 M_sun and <2.9 x 10^9 M_sun (2 sigma), comparable to their stellar masses. We detect an optically-faint line emitter situated between the two LBGs which we identify as warm molecular gas at z=5.1245 +/- 0.0001. This source, detected in the CO(2-1) transition but undetected in CO(1-0), has an integrated line flux of 0.106 +/- 0.012 Jy km/s, yielding an inferred gas mass M(H_2)=(1.9 +/- 0.2) x 10^10 M_sun. Molecular line emitters without detectable counterparts at optical and infrared wavelengths may be crucial tracers of structure and mass at high redshift.

We present a survey of bright optical dropout sources in two deep, multiwavelength surveys comprising eleven widely-separated fields, aimed at constraining the galaxy luminosity function at z~7 for sources at 5-10L*(z=6). Our combined survey area is 225 arcmin^2 to a depth of J=24.2 (AB, 3 sigma) and 135 arcmin^2 to J=25.3 (AB, 4 sigma). We find that infrared data longwards of 2 microns is essential for classifying optical dropout sources, and in particular for identifying cool Galactic star contaminants. Our limits on the number density of high redshift sources are consistent with current estimates of the Lyman break galaxy luminosity function at z=6.

A large population of z>5 Lyman break galaxies has been identified in recent years. However, the high redshift galaxies selected by different surveys are subject to a variety of selection effects - some overt, others more subtle. We present an analysis of sample completeness and contamination issues in high redshift surveys, focusing on surveys at z=5 and using a spectroscopically-confirmed low redshift sample from the DEEP2 survey in order to characterise contaminant galaxies. We find that most surveys underestimate their contamination from highly clustered galaxies at z=1 and stars. We consider the consequences of this for both the rest-frame ultraviolet luminosity function and the clustering signal from z=5 galaxies. We also find that sources with moderate strength Lyman-alpha emission lines can be omitted from dropout surveys due to their blue colours, again effecting the derived luminosity functions. We discuss the points of comparison between different samples, and the applicability of survey-specific results to the population at z>5 in general.

We present an analysis of the faint M star population seen as foreground contaminants in deep extragalactic surveys. We use space-based data to separate such stars from high redshift galaxies in a publically-available dataset, and consider the photometric properties of the resulting sample in the optical and infrared. The inferred distances place these stars well beyond the scale height of the thick disk. We find strong similarities between this faint sample (reaching i'_{AB}=25) and the brighter disk M dwarf population studied by other authors. The optical-infrared properties of the bulk of our sources spanning 6000A-4.5microns are consistent with those 5-10 magnitudes brighter. We also present deep spectroscopy of faint M dwarf stars reaching continuum limits of i'_{AB}~26, and measure absorption line strengths in the CaH2 and TiO5 bands. Both photometrically and spectroscopically, our sources are consistent with metallicities as low as a tenth solar: metal-rich compared with halo stars at similar heliocentric distances. We comment on the possible MACHO identification of M stars at faint magnitudes.

As part of a large spectroscopic survey of z>5 Lyman break galaxies (LBGs), we have identified a single source which is clearly hosting an AGN. Out of a sample of more than fifty spectroscopically-confirmed R-band dropout galaxies at z~5 and above, only J104048.6-115550.2 at z=5.44 shows evidence for a high ionisation-potential emission-line indicating the presence of a hard ionising continuum from an AGN. Like most objects in our sample the rest-frame-UV spectrum shows the UV continuum breaking across a Ly_alpha line. Uniquely within this sample of LBGs, emission from NV is also detected, a clear signature of AGN photo-ionisation. The object is spatially resolved in HST imaging. This, and the comparatively high Ly_alpha/NV flux ratio indicates that the majority of the Ly_alpha (and the UV continuum longward of it) originates from stellar photo-ionisation, a product of the ongoing starburst in the Lyman break galaxy. Even without the AGN emission, this object would have been photometrically-selected and spectroscopically-confirmed as a Lyman break in our survey. The measured optical flux (I(AB)=26.1) is therefore an upper limit to that from the AGN and is of order 100 times fainter than the majority of known quasars at these redshifts. The detection of a single object in our survey volume is consistent with the best current models of high redshift AGN luminosity function, providing a substantial fraction of such AGN are found within luminous starbursting galaxies. We discuss the cosmological implications of this discovery.

The galaxy population at z~6 has been the subject of intense study in recent years, culminating in the Hubble Ultra Deep Field (HUDF) -- the deepest imaging survey yet. A large number of high redshift galaxy candidates have been identified within the HUDF, but until now analysis of their properties has been hampered by the difficulty of obtaining spectroscopic redshifts for these faint galaxies. Our ''Gemini Lyman-Alpha at Reionisation Era'' (GLARE) project has been designed to undertake spectroscopic follow up of faint (z'<28.5) i'-drop galaxies at z~6 in the HUDF. In a previous paper we presented preliminary results from the first 7.5 hours of data from GLARE. In this paper we detail the complete survey. We have now obtained 36 hours of spectroscopy on a single GMOS slitmask from Gemini-South, with a spectral resolution of lambda/Delta(lambda) ~ 1000. We identify five strong Lyman-alpha emitters at z>5.5, and a further nine possible line emitters with detections at lower significance. We also place tight constraints on the equivalent width of Lyman-alpha emission for a further ten i'-drop galaxies and examine the equivalent width distribution of this faint spectroscopic sample of z~6 galaxies. We find that the fraction of galaxies with little or no emission is similar to that at z~3, but that the z~6 population has a tail of sources with high rest frame equivalent widths. Possible explanations for this effect include a tendency towards stronger line emission in faint sources, which may arise from extreme youth or low metallicity in the Lyman-break population at high redshift, or possibly a top-heavy initial mass function.

We measure the ages, stellar masses, and star formation histories of z~6 galaxies, observed within 1Gyr of the Big Bang. We use imaging from HST and Spitzer from the public "Great Observatories Origins Deep Survey", coupled with ground-based near-infrared imaging, to measure the spectral energy distributions from 0.8-5microns, spanning the rest-frame UV and optical. From our sample of ~50 i-drop Lyman-break star-forming galaxies in GOODS-South with z'(AB)<27mag, we focus on ~30 with reliable photometric or spectroscopic redshifts. Half of these are confused with foreground sources at Spitzer resolution, but from the 16 with clean photometry we find that a surprisingly large fraction (40%) have evidence for substantial Balmer-breaks. This indicates the presence of old underlying stellar populations that dominate the stellar masses. For these objects, we find ages of 200-700Myr, implying formation redshifts of 7<z<14, and large stellar masses in the range 1-3x10^10M_sun. Analysis of 7 i-drops that are undetected at 3.6microns indicates that these are younger, considerably less massive systems. We calculate that line contamination should not severely affect our photometry or derived results. Using data out to 8microns, we find little evidence for substantial intrinsic dust reddening. Correcting for incompleteness in our sample, we find a lower limit on the comoving stellar mass density at z~6 to be 2.5x10^6M_sun/Mpc^3. We are able to explore the star formation histories of our selected galaxies, and we infer that the past global star formation rate may have been much higher than that observed at z~6. The associated UV flux we infer at z>7 could have played a major role in reionizing the universe.

We determine the abundance of i'-band drop-outs in the HST/ACS GOODS surveys and the Hubble Ultra Deep Field (UDF). The majority of these sources are likely to be z~6 galaxies whose flux decrement arises from Lyman-alpha absorption. We have shown with Keck/DEIMOS spectroscopy that this technique does indeed select high redshift galaxies, and we discovered Lyman-alpha emission for about a third of the galaxies with z'_AB<25.6 The increased depth of UDF enables us to reach a ~10sigma limiting magnitude of z'_AB=28.5 (equivalent to 1.5M_sun/yr at z=6.1, or 0.1L*_UV for the z~3 U-drop population). The star formation rate at z~6 was approximately x6 LESS than at z~3. This declining comoving star formation rate poses an interesting challenge for models which suggest that L_UV>0.1L* star forming galaxies at z~6 reionized the universe. The short-fall in ionizing photons might be alleviated by galaxies fainter than our limit, or a radically different IMF. Alternatively, the bulk of reionization might have occurred at z>>6. We have recently discovered evidence of an early epoch of star formation in some of the i'-drops at z~6. Spitzer images with IRAC at 3.6-4.5microns show evidence of the age-sensitive Balmer/4000Ang break, dominated by stars older than 100Myr (and most probably 400Myr old). This pushes the formation epoch for these galaxies to z_form=7.5-13.5. There are at least some galaxies already assembled with stellar masses ~3x10^10M_sun (equivalent to 0.2M* today) within the first billion years. The early formation of such systems may have played a key role in reionizing the Universe at z~10.

We present new evidence for mature stellar populations with ages >100Myr in massive galaxies (M_stellar>10^10M_sun) seen at a time when the Universe was less than 1Gyr old. We analyse the prominent detections of two z~6 star-forming galaxies (SBM03#1 & #3) made at wavelengths corresponding to the rest-frame optical using the IRAC camera onboard the Spitzer Space Telescope. We had previously identified these galaxies in HST/ACS GOODS images of Chandra Deep Field South through the "i-drop" Lyman break technique, and subsequently confirmed spectroscopically with the Keck telescope. The new Spitzer photometry reveals significant Balmer/4000Ang discontinuities, indicative of dominant stellar populations with ages >100Myr. Fitting a range of population synthesis models (for normal initial mass functions) to the HST/Spitzer photometry yields ages of 250-650Myr and implied formation redshifts z~7.5-13.5 in presently-accepted world models. Remarkably, our sources have best-fit stellar masses of 1.3-3.8x10^10M_sun (95% confidence) assuming a Salpeter initial mass function. This indicates that at least some galaxies with stellar masses >20% of those of a present-day L* galaxy had already assembled within the first Gyr after the Big Bang. We also deduce that the past average star formation rate must be comparable to the current observed rate (SFR_UV~5-30M_sun/yr), suggesting that there may have been more vigorous episodes of star formation in such systems at higher redshifts. Although a small sample, limited primarily by Spitzer's detection efficiency, our result lends support to the hypothesis advocated in our earlier analyses of the Ultra Deep Field and GOODS HST/ACS data. The presence of established systems at z~6 suggests long-lived sources at earlier epochs (z>7) played a key role in reionizing the Universe.

We analyse near-infrared HST/NICMOS F110W (J) and F160W (H) band photometry of a sample of 27 i'-drop candidate z~6 galaxies in the central region of the HST/ACS Ultra Deep Field (HUDF). The infrared colours of the 20 objects not affected by near neighbours are consistent with a high redshift interpretation. This suggests that the low redshift contamination of this i'-drop sample is smaller than that observed at brighter magnitudes where values of 10-40% have been reported. The J-H colours are consistent with a slope flat in f_nu (f_lam \propto lam^{-2}), as would be expected for an unreddened starburst. There is, however, evidence for a marginally bluer spectral slope (f_lam propto lam^{-2.2) which is perhaps indicative of an extremely young starburst (~10 Myr old) or a top heavy initial mass function and little dust. The low levels of contamination, median photometric redshift of z~6.0 and blue spectral slope, inferred using the near-infrared data, supports the validity of the assumptions in our earlier work in estimating the star formation rates and, that the majority of the i-drop candidates galaxies lie at z~6.

We determine the abundance of i'-band drop-outs in the recently-released HST/ACS Hubble Ultra Deep Field (UDF). Since the majority of these sources are likely to be z~6 galaxies whose flux decrement between the F775W i'-band and F850LP z'-band arises from Lyman-alpha absorption, the number of detected candidates provides a valuable upper limit to the unextincted star formation rate at this redshift. We demonstrate that the increased depth of UDF enables us to reach an 8-sigma limiting magnitude of z'(AB)=28.5 (equivalent to 1.5/h{70}^2 M_sun/yr at z=6, or 0.1 L*(UV) for the z~3 U-drop population), permitting us to address earlier ambiguities arising from the unobserved form of the luminosity function. We identify 54 galaxies (and only one star) at z'(AB)<28.5 with (i'-z')>1.3 over the deepest 11arcmin^2 portion of the UDF field. The characteristic luminosity (L*) is consistent with values observed at z~3. The faint end slope (alpha) is less well constrained, but is consistent with only modest evolution. The main change appears to be in the number density (Phi*). Specifically, and regardless of possible contamination from cool stars and lower redshift sources, the UDF data support our previous result that the star formation rate at z~6 was at least x6 LESS than at z~3 (Stanway, Bunker & McMahon 2003). This declining comoving star formation rate (0.005 h{70}M_sun/yr/Mpc^3 at z~6 for a Salpeter IMF) poses an interesting challenge for models which suggest that L>0.1L* star forming galaxies at z~6 reionized the universe. The short-fall in ionizing photons might be alleviated by galaxies fainter than our limit, or a radically different IMF. Alternatively, the bulk of reionization might have occurred at z>>6.

A recent preprint by Yan & Windhorst (astro-ph/0407493) independently repeats our selection of candidate z~6 galaxies in the Hubble Ultra Deep Field (astro-ph/0403223). The agreement of this independent study with our original work is excellent, and most of the i-band drop-out galaxies are reproduced. To avoid confusion with various ID numbers for sources, we present the community with a matched list of i-band drop-outs in the Hubble Ultra Deep Field.

We report spectroscopic detection of three z~6 Lyman-alpha emitting galaxies, in the vicinity of the Hubble Ultra Deep Field, from the early data of the Gemini Lyman-$\alpha$ at Reionisation Era (GLARE) project. Two objects, GLARE#3001 (z =5.79) and GLARE#3011 (z =5.94), are new detections and are fainter in $z'$ (z'_{AB} =26.37 and 27.15) than any Lyman break galaxy previously detected in Lyman-alpha. A third object, GLARE#1042 (z =5.83) has previously been detected in line emission from the ground; we report here a new spectroscopic continuum detection. Gemini/GMOS-S spectra of these objects, obtained using nod & shuffle, are presented together with a discussion of their photometric properties. All three objects were selected for spectroscopy via the i-drop Lyman Break technique, the two new detections from the GOODS v1.0 imaging data. The red i'-z' colors and high equivalent widths of these objects suggest a high-confidence z>5 Lyman-alpha identification of the emission lines. This brings the total number of known z>5 galaxies within 9 arcmin of the Hubble Ultra Deep Field to four, of which three are at the same redshift (z=5.8 within 2000 km/s suggesting the existence of a large-scale structure at this redshift.

We measure the surface density of i'-band dropout galaxies at z~6 through wide field HST/ACS imaging and ultra-deep Keck/DEIMOS spectroscopy. Using deep HST/ACS SDSS-i' (F775W) and SDSS-z' (F850LP) imaging from GOODS-N (200 arcmin^2), we identify 9 i'-drops satisfying an (i'-z')_AB>1.5 selection criterion to a depth of z'_AB=25.6 (corresponding to L*_UV at z~3-4). We use HK' imaging data to improve the fidelity of our sample, discriminating against lower redshift red galaxies and cool Galactic stars. Three i'-drops are consistent with M/L/T dwarf stars. We present ultra-deep Keck/DEIMOS spectroscopy of 10 objects from our combined GOODS-N and GOODS-S i'-drop sample. We detect Lyman-alpha emission at z=5.83 from one object in the GOODS-S field, which lies only 8arcmin away (i.e. 3Mpc/h_70) from the z=5.78 object already confirmed by Bunker et al. (2003). One possible Lyman-alpha emitter at z=6.24 is found in the GOODS-N field (although identification of this spatially-offset emission line is ambiguous). Using the rest-frame UV continuum from our 6 candidate z~6 galaxies from the GOODS-N field, we determine a lower limit to the unobscured volume-averaged global star formation rate at z~6 of (5.4+/-2.2)x10^-4 h_70 M_sun/yr/Mpc^3. We find that the cosmic star formation density in galaxies with unobscured star formation rates 15M_sun/yr/h_70^2 falls by a factor of 8 between z~3 and z~6. Hence the luminosity function of LBGs must evolve in this redshift interval: a constant integrated star formation density at $z>3$ requires a much steeper faint-end slope, or a brighter characteristic luminosity. This result is in agreement with our previous measurement from the Chandra Deep Field South (Stanway et al. 2003), indicating that cosmic variance is not a dominant source of uncertainty.

We report the discovery of a luminous z=5.78 star-forming galaxy in the Chandra Deep Field South. This galaxy was selected as an `i-drop' from the GOODS public survey imaging with HST/ACS (object 3 in Stanway, Bunker & McMahon 2003, astro-ph/0302212). The large colour of (i'-z')_AB=1.6 indicated a spectral break consistent with the Lyman-alpha forest absorption short-ward of Lyman-alpha at z~6. The galaxy is very compact (marginally resolved with ACS with a half-light radius of 0.08arcsec, so r_hl<0.5kpc/h_70). We have obtained a deep (5.5-hour) spectrum of this z'(AB)=24.7 galaxy with the DEIMOS optical spectrograph on Keck, and here we report the discovery of a single emission line centred on 8245Ang detected at 20sigma with a flux of f~2E-17 ergs/s/cm^2. The line is clearly resolved with detectable structure at our resolution of better than 55km/s, and the only plausible interpretation consistent with the ACS photometry is that we are seeing Lyman-alpha emission from a z=5.78 galaxy. This is the highest redshift galaxy to be discovered and studied using HST data. The velocity width Delta(v)_FWHM=260km/s and rest-frame equivalent width (W=20Ang) indicate that this line is most probably powered by star formation, as an AGN would typically have larger values. The starburst interpretation is supported by our non-detection of the high-ionization NV1240Ang emission line, and the absence of this source from the deep Chandra X-ray images. The star formation rate inferred from the rest-frame UV continuum is 34M_sun/yr/h^2_70 (Omega_M=0.3, Omega_Lambda=0.7). This is the most luminous starburst known at z>5. Our spectroscopic redshift for this object confirms the validity of the i'-drop technique of Stanway, Bunker & McMahon (2003) to select star-forming galaxies at z~6.

We determine the space density of UV-luminous star-burst galaxies at z~6 using deep HST ACS SDSS-i' (F775W) and SDSS-z' (F850LP) and VLT ISAAC J and K_s band imaging of the Chandra Deep Field South. We find 8 galaxies and one star with (i'-z')>1.5 to a depth of z'{AB}= 25.6 (an 8-sigma detection in each of the 3 available ACS epochs). This corresponds to an unobscured star formation rate of ~15M_sun/yr/h{70}^2 at z=5.8, equivalent to L* for the Lyman break population at z =3-4 (Omega_{Lambda}=0.7, Omega_M=0.3). We are sensitive to star forming galaxies at 5.6<z<7.0 with an effective comoving volume of \~1.8x10^5 Mpc^3/h{70}^3 after accounting for incompleteness at the higher redshifts due to luminosity bias. This volume should encompass the primeval sub-galactic scale fragments of the progenitors of about a thousand L* galaxies at the current epoch. We determine a volume-averaged global star formation rate of (6.7 +/- 2.7)x10^{-4}h{70}M_sun/yr/Mpc^3 at z~6 from rest-frame UV selected star-bursts at the bright end of the luminosity function: this is a lower limit because of dust obscuration and galaxies below our sensitivity limit. This measurement shows that at z~6 the star formation density at the bright end is a factor of ~6 times less than that determined by Steidel et al. (1999) for a comparable sample of UV selected galaxies at z=3-4, and so extends our knowledge of the star formation history of the Universe to earlier times than previous work and into the epoch where reionization may have occurred.

Photometric redshifts have proven a powerful tool in identifying galaxies over a large range of lookback times. We have been generalising this technique to incorporate the selection of candidate high redshift QSOs. We have applied this to a large optical/near-infrared imaging survey in 6 wavebands aiming to push farther in redshift (and fainter in luminosity) than previous studies. We believe that study of these very faint and distant objects provides valuable insights into galaxy formation and evolution. Here we present work in progress and preliminary results for a catalogue of objects detected as part of the Las Campanas Infrared Survey. This is a stepping stone to the type of survey data that will become available in the next few years from projects such as UKIDSS and VISTA.