Source author record

E. R. Stanway

E. R. Stanway appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

31works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

31 published item(s)

preprint2022arXiv

On the Simultaneous Modelling of Dust and Stellar Populations for Interpretation of Galaxy Properties

The physical properties of galaxies are encoded within their spectral energy distribution and require comparison with models to be extracted. These models must contain a synthetic stellar population and, where infrared data is to be used, also consider prescriptions for energy reprocessing and re-emission by dust. While many such models have been constructed, there are few analyses of the impact of stellar population model choice on derived dust parameters, or vice versa. Here we apply a simple framework to compare the impact of these choices, combining three commonly-used stellar population synthesis models and three dust emission models. We compare fits to the ultraviolet to far-infrared spectral energy distributions of a validation sample of infrared-luminous galaxies. We find that including different physics, such as binary stellar evolution, in the stellar synthesis model can introduce biases and uncertainties in the derived parameters of the dust and stellar emission models, largely due to differences in the far-ultraviolet emission available for reprocessing. This may help to reconcile the discrepancy between the cosmic star formation rate and stellar mass density histories. Notably the inclusion of a dusty stellar birth cloud component in the dust emission model provides more flexibility in accommodating the stellar population model, as its reemission is highly sensitive to the ultraviolet radiation field spectrum and density. Binary populations favour a longer birth cloud dissipation timescale than is found when assuming only single star population synthesis.

preprint2022arXiv

The Dependence of Theoretical Synthetic Spectra on $α$-enhancement in Young, Binary Stellar Populations

The enhancement of $α$ elements such as oxygen is an important phase in the chemical evolution of the early Universe, with nebular material becoming enriched in these elements sooner than iron. Here we present models which incorporate stellar spectra with $α$-enhanced compositions, focusing on the impact on the integrated light of young stellar populations, including those with large binary star fractions using the Binary Populations and Spectral Synthesis (BPASS) framework, while using Solar-scaled stellar evolution models. We find that broad spectrum outputs such as production of ionising flux, the ultraviolet spectral slope and optical colours are only weakly affected by a change in [$α$/Fe]. A number of features such as ultraviolet line indices (e.g. at 1719 and 1853Å) and optical line indices (such as MgB) are sensitive to such changes in composition for a continuously star-forming population and a single starburst population respectively. We find that at ages of more than 1Gyr, $α$-enhanced stellar populations appear bluer than their Solar-scaled counterparts, and show expected sensitivity of optical line indices to composition, in agreement with previous work. The ultraviolet stellar absorption lines are relatively insensitive to subtleties in the abundances ratios, although with sufficient measurement precision, a combination of UV line indices may enable a simultaneous measurement of total metallicity mass fraction and [$α$/Fe] in young stellar populations. The output models are designated as BPASS v2.3 and made available to the community with the aim of assisting interpretation of observations of high-redshift galaxies with the James Webb Space Telescope.

preprint2022arXiv

Towards an understanding of long gamma-ray burst environments through circumstellar medium population synthesis predictions

The temporal and spectral evolution of gamma-ray burst (GRB) afterglows can be used to infer the density and density profile of the medium through which the shock is propagating. In long-duration (core-collapse) GRBs, the circumstellar medium (CSM) is expected to resemble a wind-blown bubble, with a termination shock separating the stellar wind and the interstellar medium (ISM). A long standing problem is that flat density profiles, indicative of the ISM, are often found at lower radii than expected for a massive star progenitor. Furthermore, the presence of both wind-like environments at high radii and ISM-like environments at low radii remains a mystery. In this paper, we perform a 'CSM population synthesis' with long GRB progenitor stellar evolution models. Analytic results for the evolution of wind blown bubbles are adjusted through comparison with a grid of 2D hydrodynamical simulations. Predictions for the emission radii, ratio of ISM to wind-like environments, wind and ISM densities are compared with the largest sample of afterglow-derived parameters yet compiled, which we make available for the community. We find that high ISM densities around 1000/cm3 best reproduce observations. If long GRBs instead occur in typical ISM densities of approximately 1/cm3, then the discrepancy between theory and observations is shown to persist at a population level. We discuss possible explanations for the origin of variety in long GRB afterglows, and for the overall trend of CSM modelling to over-predict the termination shock radius.

preprint2022arXiv

Where are the magnetar binary companions? Candidates from a comparison with binary population synthesis predictions

It is well established that magnetars are neutron stars with extreme magnetic fields and young ages, but the evolutionary pathways to their creation are still uncertain. Since most massive stars are in binaries, if magnetars are a frequent result of core-collapse supernovae, some fraction are expected to have a bound companion at the time of observation. In this paper, we utilise literature constraints, including deep Hubble Space Telescope imaging, to search for bound stellar companions to magnetars. The magnitude and colour measurements are interpreted in the context of binary population synthesis predictions. We find two candidates for stellar companions associated with CXOU J171405.7-381031 and SGR 0755-2933, based on their J-H colours and H-band absolute magnitudes. Overall, the proportion of the Galactic magnetar population with a plausibly stellar near-infrared counterpart candidate, based on their magnitudes and colours, is between 5 and 10 per cent. This is consistent with a population synthesis prediction of 5 per cent, for the fraction of core-collapse neutron stars arising from primaries which remain bound to their companion after the supernova. These results are therefore consistent with magnetars being drawn in an unbiased way from the natal core-collapse neutron star population, but some contribution from alternative progenitor channels cannot be ruled out.

preprint2021arXiv

Estimating Transient Rates from Cosmological Simulations and BPASS

The detection rate of electromagnetic (EM) and gravitational wave (GW) transients is growing exponentially. As the accuracy of the transient rates will significantly improve over the coming decades, so will our understanding of their evolution through cosmic history. To this end, we present predicted rates for EM and GW transients over the age of the Universe using Binary Population and Spectral Synthesis (BPASS) results combined with four cosmic star formation histories (SFH). These include a widely used empirical SFH of Madau & Dickinson and those from three cosmological simulations: MilliMillennium, EAGLE and IllustrisTNG. We find that the choice of SFH significantly changes our predictions: transients with short delay times are most affected by the star formation rate, while long-delay time events tend to depend on the metallicity evolution of star formation. Importantly we find that the cosmological simulations have very different metallicity evolution that cannot be reproduced by the widely used metallicity model of Langer & Norman, which impacts the binary black hole merger and stripped-envelope supernovae rates in the local Universe most acutely. We recommend against using simple prescriptions for the metallicity evolution of the Universe when predicting the rates of events that can have long delay times and that are sensitive to metallicity evolution.

preprint2021arXiv

The Gravitational-wave Optical Transient Observer (GOTO): prototype performance and prospects for transient science

The Gravitational-wave Optical Transient Observer (GOTO) is an array of wide-field optical telescopes, designed to exploit new discoveries from the next generation of gravitational wave detectors (LIGO, Virgo, KAGRA), study rapidly evolving transients, and exploit multi-messenger opportunities arising from neutrino and very high energy gamma-ray triggers. In addition to a rapid response mode, the array will also perform a sensitive, all-sky transient survey with few day cadence. The facility features a novel, modular design with multiple 40-cm wide-field reflectors on a single mount. In June 2017 the GOTO collaboration deployed the initial project prototype, with 4 telescope units, at the Roque de los Muchachos Observatory (ORM), La Palma, Canary Islands. Here we describe the deployment, commissioning, and performance of the prototype hardware, and discuss the impact of these findings on the final GOTO design. We also offer an initial assessment of the science prospects for the full GOTO facility that employs 32 telescope units across two sites.

preprint2020arXiv

A systematic aging method I: HII regions D118 and D119 in NGC 300

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.06}$ 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.

preprint2020arXiv

Binary Fraction Indicators in Resolved Stellar Populations and Supernova Type Ratios

The binary fraction of a stellar population can have pronounced effects on its properties, and in particular the number counts of different massive star types, and the relative subtype rates of the supernovae which end their lives. Here we use binary population synthesis models with a binary fraction that varies with initial mass to test the effects on resolved stellar populations and supernovae, and ask whether these can constrain the poorly-known binary fraction in different mass and metallicity regimes. We show that Wolf-Rayet star subtype ratios are valuable binary diagnostics, but require large samples to distinguish by models. Uncertainties in which stellar models would be spectroscopically classified as Wolf-Rayet stars are explored. The ratio of thermonuclear, stripped envelope and other core-collapse supernovae may prove a more accessible test and upcoming surveys will be sufficient to constrain both the high mass and low mass binary fraction in the z < 1 galaxy population.

preprint2020arXiv

Evaluating the impact of binary parameter uncertainty on stellar population properties

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.

preprint2020arXiv

Searching for Electromagnetic Counterparts to Gravitational-wave Merger Events with the Prototype Gravitational-wave Optical Transient Observer (GOTO-4)

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.

preprint2020arXiv

Weighing in on black hole binaries with BPASS: LB-1 does not contain a 70M$_{\odot}$ black hole

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.

preprint2019arXiv

Supernova lightCURVE POPulation Synthesis II: Validation against supernovae with an observed progenitor

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.

preprint2016arXiv

A Hubble Space Telescope Survey of the Host Galaxies of Superluminous Supernovae

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.

preprint2016arXiv

A radio pulsing white dwarf binary star

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.

preprint2016arXiv

BPASS predictions for Binary Black-Hole Mergers

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.

preprint2016arXiv

Detection of three Gamma-Ray Burst host galaxies at $z\sim6$

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, $λ_{\rm{obs}}\sim1.4\,μm$) filter. The hosts have magnitudes (corrected for Galactic extinction) of $m_{\rm{λ_{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.

preprint2015arXiv

Dissecting the complex environment of a distant quasar with MUSE

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.

preprint2015arXiv

Late time multi wavelength observations of Swift J1644+5734: A luminous optical/IR bump and quiescent X-ray emission

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.

preprint2015arXiv

Swift J1112.2-8238: A Candidate Relativistic Tidal Disruption Flare

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) $γ$-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/$γ$-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.

preprint2014arXiv

A Detection of Molecular Gas Emission in the Host Galaxy of GRB 080517

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.

preprint2014arXiv

ALMA Observations of the Host Galaxy of GRB090423 at z=8.23: Deep Limits on Obscured Star Formation 630 Million Years After the Big Bang

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.

preprint2014arXiv

GRB 080517: A local, low luminosity GRB in a dusty galaxy at z=0.09

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.

preprint2014arXiv

Interpreting high [O III]/Hbeta ratios with maturing starbursts

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.

preprint2013arXiv

Are z~5 QSOs found in the most massive high redshift over-densities?

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.

preprint2013arXiv

Identifying clustering at high redshift through actively star-forming galaxies

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.

preprint2013arXiv

The detection of FIR emission from high redshift star-forming galaxies in the ECDF-S

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.

preprint2012arXiv

Distant galaxy clusters in the XMM Large Scale Structure survey

(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.

preprint2012arXiv

Limits on dust emission from z~5 LBGs and their local environments

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.

preprint2012arXiv

Star formation in the early universe: beyond the tip of the iceberg

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.

preprint2011arXiv

MESMER: MeerKAT Search for Molecules in the Epoch of Reionization

[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.

preprint2009arXiv

Photometric Selection of z~5 Lyman Break Galaxies in the ESO Remote Galaxy Survey

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.