Source author record

Stephen M. Wilkins

Stephen M. Wilkins 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

30works
3topics
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

30 published item(s)

preprint2022arXiv

First Light And Reionisation Epoch Simulations (FLARES) III: The properties of massive dusty galaxies at cosmic dawn

Using the First Light And Reionisation Epoch Simulations (\textsc{Flares}) we explore the dust driven properties of massive high-redshift galaxies at $z\in[5,10]$. By post-processing the galaxy sample using the radiative transfer code \textsc{skirt} we obtain the full spectral energy distribution. We explore the resultant luminosity functions, IRX-$β$ relations as well as the luminosity-weighted dust temperatures in the Epoch of Reionisation (EoR). We find that most of our results are in agreement with the current set of observations, but under-predict the number densities of bright IR galaxies, which are extremely biased towards the most overdense regions. We see that the \textsc{Flares} IRX-$β$ relation (for $5\le z\le8$) predominantly follows the local starburst relation. The IRX shows an increase with stellar mass, plateauing at the high-mass end ($\sim10^{10}$M$_{\odot}$) and shows no evolution in the median normalisation with redshift. We also look at the dependence of the peak dust temperature ($T_{\mathrm{peak}}$) on various galaxy properties including the stellar mass, IR luminosity and sSFR, finding the correlation to be strongest with sSFR. The luminosity-weighted dust temperatures increase towards higher redshifts, with the slope of the $T_{\mathrm{peak}}$ - redshift relation showing a higher slope than the lower redshift relations obtained from previous observational and theoretical works. The results from \textsc{Flares}, which is able to provide a better statistical sample of high-redshift galaxies compared to other simulations, provides a distinct vantage point for the high-redshift Universe.

preprint2022arXiv

First Light And Reionisation Epoch Simulations (FLARES) IV: The size evolution of galaxies at $z\geq5$

We present the intrinsic and observed sizes of galaxies at $z\geq5$ in the First Light And Reionisation Epoch Simulations (FLARES). We employ the large effective volume of FLARES to produce a sizeable sample of high redshift galaxies with intrinsic and observed luminosities and half light radii in a range of rest frame UV and visual photometric bands. This sample contains a significant number of intrinsically ultra-compact galaxies in the far-UV (1500 angstrom), leading to a negative intrinsic far-UV size-luminosity relation. However, after the inclusion of the effects of dust these same compact galaxies exhibit observed sizes that are as much as 50 times larger than those measured from the intrinsic emission, and broadly agree with a range of observational samples. This increase in size is driven by the concentration of dust in the core of galaxies, heavily attenuating the intrinsically brightest regions. At fixed luminosity we find a galaxy size redshift evolution with a slope of $m=1.21-1.87$ depending on the luminosity sample in question, and we demonstrate the wavelength dependence of the size-luminosity relation which will soon be probed by the Webb Space Telescope.

preprint2022arXiv

On the Stellar Populations of Galaxies at z=9-11: The Growth of Metals and Stellar Mass at Early Times

We present a detailed stellar population analysis of 11 bright ($H<26.6$) galaxies at $z=9-11$ (three spectroscopically confirmed) to constrain the chemical enrichment and growth of stellar mass of early galaxies. We use the flexible Bayesian spectral energy distribution (SED) fitting code Prospector with a range of star-formation histories (SFHs), a flexible dust attenuation law and a self-consistent modeling of emission lines. This approach allows us to assess how different priors affect our results, and how well we can break degeneracies between dust attenuation, stellar ages, metallicity and emission lines using data which probe only the rest-frame ultraviolet to optical wavelengths. We measure a median observed ultraviolet spectral slope $β=-1.87_{-0.43}^{+0.35}$ for relatively massive star-forming galaxies ($9<\log(M_{\star}/M_{\odot})<10$), consistent with no change from $z=4$ to $z=9-10$ at these stellar masses, implying rapid enrichment. Our SED-fitting results are consistent with a star-forming main sequence with sub-linear slope ($0.7\pm0.2$) and specific star-formation rates of $3-10~\mathrm{Gyr}^{-1}$. However, the stellar ages and SFHs are less well constrained. Using different SFH priors, we cannot distinguish between median mass-weighted ages of $\sim50-150$ Myr, which corresponds to 50\% formation redshifts of $z_{50}\sim10-12$ at $z\sim9$ and is of the order of the dynamical timescales of these systems. Importantly, the models with different SFH priors are able to fit the data equally well. We conclude that the current observational data cannot tightly constrain the mass-buildup timescales of these $z=9-11$ galaxies, with our results consistent with SFHs implying both a shallow and steep increase of the cosmic SFR density with time at $z>10$.

preprint2022arXiv

The COS Legacy Archive Spectroscopy SurveY (CLASSY) Treasury Atlas

Far-ultraviolet (FUV; ~1200-2000 angstroms) spectra are fundamental to our understanding of star-forming galaxies, providing a unique window on massive stellar populations, chemical evolution, feedback processes, and reionization. The launch of JWST will soon usher in a new era, pushing the UV spectroscopic frontier to higher redshifts than ever before, however, its success hinges on a comprehensive understanding of the massive star populations and gas conditions that power the observed UV spectral features. This requires a level of detail that is only possible with a combination of ample wavelength coverage, signal-to-noise, spectral-resolution, and sample diversity that has not yet been achieved by any FUV spectral database. We present the COS Legacy Spectroscopic SurveY (CLASSY) treasury and its first high level science product, the CLASSY atlas. CLASSY builds on the HST archive to construct the first high-quality (S/N_1500 >~ 5/resel), high-resolution (R~15,000) FUV spectral database of 45 nearby (0.002 < z < 0.182) star-forming galaxies. The CLASSY atlas, available to the public via the CLASSY website, is the result of optimally extracting and coadding 170 archival+new spectra from 312 orbits of HST observations. The CLASSY sample covers a broad range of properties including stellar mass (6.2 < logM_star(M_sol) < 10.1), star formation rate (-2.0 < log SFR (M_sol/yr) < +1.6), direct gas-phase metallicity (7.0 < 12+log(O/H) < 8.8), ionization (0.5 < O_32 < 38.0), reddening (0.02 < E(B-V < 0.67), and nebular density (10 < n_e (cm^-3) < 1120). CLASSY is biased to UV-bright star-forming galaxies, resulting in a sample that is consistent with z~0 mass-metallicity relationship, but is offset to higher SFRs by roughly 2 dex, similar to z >~2 galaxies. This unique set of properties makes the CLASSY atlas the benchmark training set for star-forming galaxies across cosmic time.

preprint2022arXiv

Webb's PEARLS: Bright 1.5--2.0 micron Dropouts in the Spitzer/IRAC Dark Field

Using the first epoch of four-band NIRCam observations obtained by the James Webb Space Telescope (JWST) Prime Extragalactic Areas for Reionization and Lensing Science Program in the Spitzer IRAC Dark Field, we search for F150W and F200W dropouts. In 14.2 arcmin^2, we have found eight F150W dropouts and eight F200W dropouts, all brighter than 27.5 mag (the brightest being ~24 mag) in the band to the red side of the break. As they are detected in multiple bands, these must be real objects. Their nature, however, is unclear, and characterizing their properties is important for realizing the full potential of JWST. If the observed color decrements are due to the Lyman break, these objects should be at z >~ 11.7 and z >~ 15.4, respectively. The color diagnostics show that at least four F150W dropouts are far away from the usual contaminators encountered in dropout searches (red galaxies at much lower redshifts or brown dwarf stars). While the diagnostics of the F200W dropouts are less certain due to the limited number of passbands, at least one of them is likely not a known type of contaminant, and the rest are consistent with either high-redshift galaxies with evolved stellar populations or old galaxies at z ~ 3 to 8. If a significant fraction of our dropouts are indeed at z ~ 12, we have to face the severe problem of explaining their high luminosities and number densities. Spectroscopic identifications of such objects are urgently needed.

preprint2020arXiv

Nebular Line Emission During the Epoch of Reionization

Nebular emission lines associated with galactic HII regions carry information about both physical properties of the ionised gas and the source of ionising photons as well as providing the opportunity of measuring accurate redshifts and thus distances once a cosmological model is assumed. While nebular line emission has been extensively studied at lower redshift there are currently only few constraints within the epoch of reionisation (EoR, $z>6$), chiefly due to the lack of sensitive near-IR spectrographs. However, this will soon change with the arrival of the Webb Telescope providing sensitive near-IR spectroscopy covering the rest-frame UV and optical emission of galaxies in the EoR. In anticipation of Webb we combine the large cosmological hydrodynamical simulation Bluetides with photoionisation modelling to predict the nebular emission line properties of galaxies at $z=8\to 13$. We find good agreement with the, albeit limited, existing direct and indirect observational constraints on equivalent widths though poorer agreement with luminosity function constraints.

preprint2016arXiv

Forecasts for the WFIRST High Latitude Survey using the BlueTides Simulation

We use the BlueTides simulation to predict the properties of the high-$z$ galaxy and active galactic nuclei (AGN) populations for the planned 2200deg$^2$ Wide-Field Infrared Survey Telescope's (WFIRST)-AFTA High Latitude Survey (HLS). BlueTides is a cosmological hydrodynamic simulation, which incorporates a variety of baryon physics in a $(400h^{-1} \mathrm{Mpc})^3$ volume evolved to $z=8$ with 0.7 trillion particles. The galaxy luminosity functions in the simulation show good agreement with all the current observational constraints (up to $z=11$) and predicts an enhanced number of UV bright galaxies. At the proposed depth of the HLS ($m < 26.75$), BlueTides predicts $10^6$ galaxies at $z=8$ with a few up to $z\sim 15$ due to the enhanced bright end of the galaxy luminosity function. At $z=8$, galaxies in the mock HLS have specific star formation rates of $\sim 10 {\rm Gyr}^{-1}$ and ages of $\sim 80 {\rm Myr}$ (both evolving linearly with redshift) and a non-evolving mass-metallicity relation. BlueTides also predicts $\sim 10^4$ AGN in WFIRST HLS from $z=8$ out to $z\sim 14$. These AGN host black holes of $M\sim 10^6-10^8 M_\odot$ accreting close to their Eddington luminosity. Galaxies and AGN have host halo masses of $M_{halo}\sim 10^{11-12} M_\odot$ and a linear bias $b\approx 13-20$. Given the expected galaxy space densities, their high bias and large volume probed we speculate that it may be feasible for WFIRST HLS detect the Baryon Acoustic Oscillation peak in the galaxy power spectrum out to $z=8-9$.

preprint2016arXiv

Inferred Hα Flux as a Star-Formation Rate Indicator at z ~ 4-5: Implications for Dust Properties, Burstiness, and the z = 4-8 Star-Formation-Rate Functions

We derive Hα fluxes for a large spectroscopic and photometric-redshift-selected sample of sources over GOODS-North and South in the redshift range z = 3.8-5.0 with deep HST, Spitzer/IRAC, and ground-based observations. The Hα flux is inferred based on the offset between the IRAC 3.6 μm flux and that predicted from the best-fit SED. We demonstrate that the Hα flux correlates well with dust- corrected UV star-formation rate (SFR) and therefore can serve as an independent SFR indicator. However, we also find a systematic offset in the SFR_Hα/SFR_UV ratios for z ~ 4-5 galaxies relative to local relations (assuming the same dust corrections for nebular regions and stellar light). We show that we can resolve the modest tension in the inferred SFRs by assuming bluer intrinsic UV slopes (increasing the dust correction), a rising star-formation history or assuming a low metallicity stellar population with a hard ionizing spectrum (increasing the L_Hα/SFR ratio). Using Hα as a SFR indicator, we find a normalization of the star formation main sequence in good agreement with recent SED-based determinations and also derive the SFR functions at z ~ 4-8. In addition, we assess for the first time the burstiness of star formation in z ~ 4 galaxies on <100 Myr time scales by comparing UV and Hα-based sSFRs; their one-to-one relationship argues against significantly bursty star-formation histories. Further progress will be made on these results, by incorporating new results from ALMA to constrain the dust-obscured star formation in high-redshift UV-selected samples.

preprint2016arXiv

The Lyman-continuum photon production efficiency in the high-redshift Universe

The Lyman Continuum photon production efficiency ($ξ_{\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}(ξ_{\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

preprint2016arXiv

The Photometric Properties of Galaxies in the Early Universe

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

preprint2015arXiv

Exploring the Dawn of Galaxies

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.

preprint2015arXiv

Galaxy And Mass Assembly (GAMA): Panchromatic Data Release (far-UV --- far-IR) and the low-z energy budget

We present the GAMA Panchromatic Data Release (PDR) constituting over 230deg$^2$ of imaging with photometry in 21 bands extending from the far-UV to the far-IR. These data complement our spectroscopic campaign of over 300k galaxies, and are compiled from observations with a variety of facilities including: GALEX, SDSS, VISTA, WISE, and Herschel, with the GAMA regions currently being surveyed by VST and scheduled for observations by ASKAP. These data are processed to a common astrometric solution, from which photometry is derived for 221,373 galaxies with r<19.8 mag. Online tools are provided to access and download data cutouts, or the full mosaics of the GAMA regions in each band. We focus, in particular, on the reduction and analysis of the VISTA VIKING data, and compare to earlier datasets (i.e., 2MASS and UKIDSS) before combining the data and examining its integrity. Having derived the 21-band photometric catalogue we proceed to fit the data using the energy balance code MAGPHYS. These measurements are then used to obtain the first fully empirical measurement of the 0.1-500$μ$m energy output of the Universe. Exploring the Cosmic Spectral Energy Distribution (CSED) across three time-intervals (0.3-1.1Gyr, 1.1-1.8~Gyr and 1.8---2.4~Gyr), we find that the Universe is currently generating $(1.5 \pm 0.3) \times 10^{35}$ h$_{70}$ W Mpc$^{-3}$, down from $(2.5 \pm 0.2) \times 10^{35}$ h$_{70}$ W Mpc$^{-3}$ 2.3~Gyr ago. More importantly, we identify significant and smooth evolution in the integrated photon escape fraction at all wavelengths, with the UV escape fraction increasing from 27(18)% at z=0.18 in NUV(FUV) to 34(23)% at z=0.06. The GAMA PDR will allow for detailed studies of the energy production and outputs of individual systems, sub-populations, and representative galaxy samples at $z<0.5$. The GAMA PDR can be found at: http://gama-psi.icrar.org/

preprint2015arXiv

Luminosity function of [OII] emission-line galaxies in the MassiveBlack-II simulation

We examine the luminosity function (LF) of [OII] emission-line galaxies in the high-resolution cosmological simulation MassiveBlack-II (MBII). From the spectral energy distribution of each galaxy, we select a sub-sample of star-forming galaxies at $0.06 \le z \le 3.0$ using the [OII] emission line luminosity L([OII]). We confirm that the specific star formation rate matches that in the GAMA survey. We show that the [OII] LF at z=1.0 from the MBII shows a good agreement with the LFs from several surveys below L([OII])=$10^{43.0}$ erg/s while the low redshifts ($z \le 0.3$) show an excess in the prediction of bright [OII] galaxies, but still displaying a good match with observations below L([OII])=$10^{41.6}$ erg/s. Based on the validity in reproducing the properties of [OII] galaxies at low redshift ($z \le 1$), we forecast the evolution of the [OII] LF at high redshift ($z \le 3$), which can be tested by upcoming surveys such as the HETDEX and DESI. The slopes of the LFs at bright and faint ends range from -3 to -2 showing minima at z=2. The slope of the bright end evolves approximately as 1/(z+1) at z=2 while the faint end evolves as ~3/(z+1) at $0.6 \le z \le 2$. In addition, a similar analysis is applied for the evolution of [OIII] LFs, which is to be explored in the forthcoming survey WFIRST-AFTA. Finally, we show that the auto-correlation function of [OII] and [OIII] emitting galaxies shows a rapid evolution from z=2 to 1.

preprint2015arXiv

Quantifying the UV continuum slopes of galaxies to z~10 using deep Hubble and Spitzer/IRAC observations

Measurements of the UV-continuum slopes provide valuable information on the physical properties of galaxies forming in the early universe, probing the dust reddening, age, metal content, and even the escape fraction. While constraints on these slopes generally become more challenging at higher redshifts as the UV continuum shifts out of the Hubble Space Telescope bands (particularly at z>7), such a characterisation actually becomes abruptly easier for galaxies in the redshift window z=9.5-10.5 due to the Spitzer/IRAC 3.6um-band probing the rest-UV continuum and the long wavelength baseline between this Spitzer band and the Hubble H-band. Higher S/N constraints on the UV slope are possible at z~10 than at z=8. Here we take advantage of this opportunity and five recently discovered bright z=9.5-10.5 galaxies to present the first measurements of the mean slope for a multi-object sample of galaxy candidates at z~10. We find the measured observed slopes of these candidates are $-2.1\pm0.3\pm0.2$ (random and systematic), only slightly bluer than the measured slopes at 3.5<z<7.5 for galaxies of similar luminosities. Small increases in the stellar ages, metallicities, and dust content of the galaxy population from z~10 to z~7 could easily explain the apparent evolution in slopes.

preprint2014arXiv

Galaxy And Mass Assembly (GAMA): mass-size relations of z$<$0.1 galaxies subdivided by Sérsic index, colour and morphology

We use data from the Galaxy And Mass Assembly (GAMA) survey in the redshift range 0.01$<$z$<$0.1 (8399 galaxies in $g$ to $K_s$ bands) to derive the stellar mass $-$ half-light radius relations for various divisions of 'early' and 'late'-type samples. We find the choice of division between early and late (i.e., colour, shape, morphology) is not particularly critical, however, the adopted mass limits and sample selections (i.e., the careful rejection of outliers and use of robust fitting methods) are important. In particular we note that for samples extending to low stellar mass limits ($<10^{10}\mathcal{M_{\odot}}$) the Sérsic index bimodality, evident for high mass systems, becomes less distinct and no-longer acts as a reliable separator of early- and late-type systems. The final set of stellar mass $-$ half-light radius relations are reported for a variety of galaxy population subsets in 10 bands ($ugrizZYJHK_s$) and are intended to provide a comprehensive low-z benchmark for the many ongoing high-z studies. Exploring the variation of the stellar mass $-$ half-light radius relations with wavelength we confirm earlier findings that galaxies appear more compact at longer wavelengths albeit at a smaller level than previously noted: at $10^{10}\mathcal{M_{\odot}}$ both spiral systems and ellipticals show a decrease in size of 13% from $g$ to $K_s$ (which is near linear in log wavelength). Finally we note that the sizes used in this work are derived from 2D Sérsic light profile fitting (using GALFIT3), i.e., elliptical semi-major half light radii, improving on earlier low-z benchmarks based on circular apertures.

preprint2014arXiv

High-redshift galaxies and low-mass stars

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.

preprint2014arXiv

Spectroscopy of z ~ 7 candidate galaxies: Using Lyman-alpha to constrain the neutral fraction of hydrogen in the high-redshift universe

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-$α$ in the redshift range spanning 5.06 - 7.80. We target 22 $z$-band dropouts and find no evidence of Lyman-$α$ emission, with the exception of a tentative detection ($<5σ$, which is our adopted criterion for a secure detection) for one object. The upper limits on Lyman-$α$ 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-$α$, with an estimated neutral fraction $χ_{HI}\sim0.5$, in agreement with other estimates.

preprint2014arXiv

The MassiveBlack-II Simulation: The Evolution of Halos and Galaxies to z~0

(Abridged for arXiv)We investigate the properties of halos, galaxies and blackholes to z=0 in the high resolution hydrodynamical simulation MassiveBlack-II (MBII) which evolves a LCDM cosmology in a comoving volume Vbox=100(Mpc/h)^3. MBII is the highest resolution simulation of this size which includes a self-consistent model for star formation, black hole accretion and associated feedback. We provide a simulation browser web application which enables interactive search and tagging of halos, subhalos and their properties and publicly release our galaxy catalogs. Our analysis of the halo mass function (MF) in MBII reveals that baryons have strong effects, with changes in the halo abundance of 20-35% below the knee of the MF (Mhalo < 10^13.2 Msun/h at z=0) when compared to fits based on dark matter only simulations. We provide a fitting function for the halo MF out to redshift z=11 and discuss how the onset of non-universality in the MF limits the accuracy of our fit. We study the halo occupation distribution and clustering of galaxies, in particular the evolution and scale dependence of stochasticity and bias finding reasonable agreement with observational data. The shape of the cosmic spectral energy distribution predicted by MBII is consistent with observations, but lower in amplitude. The Galaxy Stellar Mass Function (GSMF) function is broadly consistent with observations at z>=2. At z<2, the population of passive low mass (for M*<10^9 Msun) galaxies in MBII makes the GSMF too steep compared to observations whereas at the high mass end (M*>10^11 Msun) galaxies hosting bright AGNs make significant contributions to the GSMF. The quasar bolometric luminosity function is also largely consistent with observations. We note however that more efficient AGN feedback (beyond simple thermal coupling used here) is likely necessary for the largest, rarest objects/clusters at low redshifts.

preprint2013arXiv

Confronting predictions of the galaxy stellar mass function with observations at high-redshift

We investigate the evolution of the galaxy stellar mass function at high-redshift ($z\ge 5$) using a pair of large cosmological hydrodynamical simulations: {\em MassiveBlack} and {\em MassiveBlack-II}. By combining these simulations we can study the properties of galaxies with stellar masses greater than $10^{8}\,{\rm M_{\odot}}\,h^{-1}$ and (co-moving) number densities of $\log_{10}(ϕ\, [{\rm Mpc^{-3}\,dex^{-1}}\,h^{3}])>-8$. Observational determinations of the galaxy stellar mass function at very-high redshift typically assume a relation between the observed UV luminosity and stellar mass-to-light ratio which is applied to high-redshift samples in order to estimate stellar masses. This relation can also be measured from the simulations. We do this, finding two significant differences with the usual observational assumption: it evolves strongly with redshift and has a different shape. Using this relation to make a consistent comparison between galaxy stellar mass functions we find that at $z=6$ and above the simulation predictions are in good agreement with observed data over the whole mass range. Without using the correct UV luminosity and stellar mass-to-light ratio, the discrepancy would be up to two orders of magnitude for large galaxies $>10^{10}\,{\rm M_{\odot}}\,h^{-1}$. At $z=5$, however the stellar mass function for low mass $<10^{9}\,{\rm M_{\odot}}\,h^{-1}$ galaxies is overpredicted by factors of a few, consistent with the behaviour of the UV luminosity function, and perhaps a sign that feedback in the simulation is not efficient enough for these galaxies.

preprint2013arXiv

Interpreting the observed UV continuum slopes of high-redshift galaxies

The observed UV continuum slope of star forming galaxies is strongly affected by the presence of dust. Its observation is then a potentially valuable diagnostic of dust attenuation, particularly at high-redshift where other diagnostics are currently inaccesible. Interpreting the observed UV continuum slope in the context of dust attenuation is often achieved assuming the empirically calibrated Meurer et al. (1999) relation. Implicit in this relation is the assumption of an intrinsic UV continuum slope ($β=-2.23$). However, results from numerical simulations suggest that the intrinsic UV continuum slopes of high-redshift star forming galaxies are bluer than this, and moreover vary with redshift. Using values of the intrinsic slope predicted by numerical models of galaxy formation combined with a Calzetti et al. (2000) reddening law we infer UV attenuations ($A_{1500}$) $0.35-0.5\,{\rm mag}$ ($A_{V}$: $0.14-0.2\,{\rm mag}$ assuming Calzetti et al. 2000) greater than simply assuming the Meurer relation. This has significant implications for the inferred amount of dust attenuation at very-high ($z\approx 7$) redshift given current observational constraints on $β$, combined with the Meurer relation, suggest dust attenuation to be virtually zero in all but the most luminous systems.

preprint2013arXiv

Single Colour Diagnostics of the Mass-to-light Ratio: Predictions from Galaxy Formation Models

Accurate galaxy stellar masses are crucial to better understand the physical mechanisms driving the galaxy formation process. We use synthetic star formation and metal enrichment histories predicted by the {\sc galform} galaxy formation model to investigate the precision with which various colours $(m_{a}-m_{b})$ can alone be used as diagnostics of the stellar mass-to-light ratio. As an example, we find that, at $z=0$, the {\em intrinsic} (B$_{f435w}-$V$_{f606w}$) colour can be used to determine the intrinsic rest-frame $V$-band stellar mass-to-light ratio ($\log_{10}Γ_{V}=\log_{10}[(M/M_{\odot})/(L_{V}/L_{V\odot})]$) with a precision of $σ_{lgΓ}\simeq 0.06$ when the initial mass function and redshift are known beforehand. While the presence of dust, assuming a universal attenuation curve, can have a systematic effect on the inferred mass-to-light ratio using a single-colour relation, this is typically small as it is often possible to choose a colour for which the dust reddening vector is approximately aligned with the $(m_{a}-m_{b})-\log_{10}Γ_{V}$ relation. The precision with which the stellar mass-to-light ratio can be recovered using a single colour diagnostic rivals implementations of SED fitting using more information but in which simple parameterisations of the star formation and metal enrichment histories are assumed. To facilitate the wide use of these relations, we provide the optimal observer frame colour to estimate the stellar mass-to-light ratio, along with the associated parameters, as a function of redshift ($0<z<1.5$) for two sets of commonly used filters.

preprint2013arXiv

Theoretical predictions for the effect of nebular emission on the broad band photometry of high-redshift galaxies

By combining optical and near-IR observations from the Hubble Space Telescope with NIR photometry from the Spitzer Space Telescope it is possible to measure the rest-frame UV-optical colours of galaxies at z=4-8. The UV-optical spectral energy distribution of star formation dominated galaxies is the result of several different factors. These include the joint distribution of stellar masses, ages, and metallicities, and the subsequent reprocessing by dust and gas in the ISM. Using a large cosmological hydrodynamical simulation we investigate the predicted spectral energy distributions of galaxies at high-redshift with a particular emphasis on assessing the potential contribution of nebular emission. We find that the average pure stellar UV-optical colour correlates with both luminosity and redshift such that galaxies at lower-redshift and higher-luminosity are typically redder. Assuming the escape fraction of ionising photons is close to zero, the effect of nebular emission is to redden the UV-optical 1500-V_w colour by, on average, 0.4 mag at z=8 declining to 0.25 mag at z=4. Young and low-metallicity stellar populations, which typically have bluer pure stellar UV-optical colours, produce larger ionising luminosities and are thus more strongly affected by the reddening effects of nebular emission. This causes the distribution of 1500-V_w colours to narrow and the trends with luminosity and redshift to weaken. The strong effect of nebular emission leaves observed-frame colours critically sensitive to the source redshift. For example, increasing the redshift by 0.1 can result in observed frame colours changing by up to ~0.6. These predictions reinforce the need to include nebular emission when modelling the spectral energy distributions of galaxies at high-redshift and also highlight the difficultly in interpreting the observed colours of individual galaxies without precise redshifts.

preprint2013arXiv

VLT/XSHOOTER & Subaru/MOIRCS Spectroscopy of HUDF-YD3: No Evidence for Lyman-alpha Emission at z=8.55

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.

preprint2012arXiv

Constraining the Bright-end of the UV Luminosity Function for z ~ 7 - 9 Galaxies: results from CANDELS/GOODS-South

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.

preprint2012arXiv

No Evidence for Lyman-alpha Emission in Spectroscopy of z > 7 Candidate Galaxies

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.

preprint2012arXiv

Predictions for the intrinsic UV continuum properties of star forming galaxies and the implications for inferring dust extinction

The observed ultraviolet continuum (UVC) slope is potentially a powerful diagnostic of dust obscuration in star forming galaxies. However, the intrinsic slope is also sensitive to the form of the stellar initial mass function (IMF) and to the recent star formation and metal enrichment histories of a galaxy. Using the galform semi-analytical model of galaxy formation, we investigate the intrinsic distribution of UVC slopes. For star-forming galaxies, we find that the intrinsic distribution of UVC slopes at z=0, parameterised by the power law index beta, has a standard deviation of sigma_beta=0.30. This suggests an uncertainty on the inferred UV attenuation of A_fuv=0.7$ (assuming a Calzetti attenuation curve) for an individual object, even with perfect photometry. Furthermore, we find that the intrinsic UVC slope correlates with star formation rate, intrinsic UV luminosity, stellar mass and redshift. These correlations have implications for the interpretation of trends in the observed UVC slope with these quantities irrespective of the sample size or quality of the photometry. Our results suggest that in some cases the attenuation by dust has been incorrectly estimated.

preprint2012arXiv

The accuracy of the UV continuum as an indicator of the star formation rate in galaxies

The rest-frame intrinsic UV luminosity is often used as an indicator of the instantaneous star formation rate (SFR) in a galaxy. While it is in general a robust indicator of the ongoing star formation activity, the precise value of the calibration relating the UV luminosity to the SFR ($B_ν$), is sensitive to various physical properties, such as the recent star formation and metal enrichment histories, along with the choice of stellar initial mass function. The distribution of these properties for the star-forming galaxy population then suggests that the adoption of a single calibration is not appropriate unless properly qualified with the uncertainties on the calibration. We investigate, with the aid of the {\sc galform} semi-analytic model of galaxy formation, the distribution of UV-SFR calibrations obtained using realistic star formation and metal enrichment histories. At $z=0$, we find that when the initial mass function is fixed (to the Kennicutt IMF), the median calibration is $B_{\rm fuv}=0.9$ where ${\rm SFR}/[{\rm M_{\odot}\,yr^{-1}}]=B_ν\times 10^{-28}\times L_ν/[{\rm ergs\,s^{-1}\,Hz^{-1}}]$. However, the width of the distribution $B_{\rm fuv}$ suggests that for a single object there is around a 20% {\em intrinsic} uncertainty (at $z=0$, rising to $\simeq 30%$ at $z=6$) on the star formation rate inferred from the FUV luminosity without additional constraints on the star formation history or metallicity. We also find that the median value of the calibration $B_{\rm fuv}$ is correlated with the star formation rate and redshift (at $z>3$) raising implications for the correct determination of the star formation rate from the UV.

preprint2011arXiv

The UV Properties of Star Forming Galaxies I: HST WFC3 Observations of Very-high Redshift Galaxies

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-$α$ emission or the Lyman-$α$ 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.

preprint2010arXiv

New Star Forming Galaxies at z\approx 7 from WFC3 Imaging

The addition of Wide Field Camera 3 (WFC3) on the Hubble Space Telescope (HST) has led to a dramatic increase in our ability to study the z>6 Universe. The increase in the near-infrared (NIR) sensitivity of WFC3 over previous instruments has enabled us to reach apparent magnitudes approaching 29 (AB). This allows us to probe the rest-frame ultraviolet (UV) continuum, redshifted into the NIR at $z>6$. Taking advantage of the large optical depths at this redshift, resulting in the Lyman-alpha break, we use a combination of WFC3 imaging and pre-existing Advanced Camera for Surveys (ACS) imaging to search for z approx 7 over 4 fields. Our analysis reveals 29 new z approx 7 star forming galaxy candidates in addition to 16 pre-existing candidates already discovered in these fields. The improved statistics from our doubling of the robust sample of z-drop candidates confirms the previously observed evolution of the bright end of the luminosity function.

preprint2009arXiv

Probing $\sim L_{*}$ Lyman-break Galaxies at $z\approx 7$ in GOODS-South with WFC3 on HST

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.