Researcher profile

Stephen Wilkins

Stephen Wilkins contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
14works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

14 published item(s)

preprint2022arXiv

The BlueTides Mock Image Catalogue: Simulated observations of high-redshift galaxies and predictions for JWST imaging surveys

We present a mock image catalogue of ~100,000 MUV=-22.5 to -19.6 mag galaxies at z=7-12 from the BlueTides cosmological simulation. We create mock images of each galaxy with the James Webb (JWST), Hubble, Roman, and Euclid Space Telescopes, as well as Subaru, and VISTA, with a range of near- and mid-infrared filters. We perform photometry on the mock images to estimate the success of these instruments for detecting high-z galaxies. We predict that JWST will have unprecedented power in detecting high-z galaxies, with a 95% completeness limit at least 2.5 magnitudes fainter than VISTA and Subaru, 1.1 magnitudes fainter than Hubble, and 0.9 magnitudes fainter than Roman, for the same wavelength and exposure time. Focusing on JWST, we consider a range of exposure times and filters, and find that the NIRCam F356W and F277W filters will detect the faintest galaxies, with 95% completeness at m=27.4 mag in 10ks exposures. We also predict the number of high-z galaxies that will be discovered by upcoming JWST imaging surveys. We predict that the COSMOS-Web survey will detect ~1000 MUV<-20.1 mag galaxies at 6.5<z<7.5, by virtue of its large survey area. JADES-Medium will detect almost 100% of MUV<-20 mag galaxies at z<8.5 due to its significant depth, however with its smaller survey area it will detect only ~100 of these galaxies at 6.5<z<7.5. Cosmic variance results in a large range in the number of predicted galaxies each survey will detect, which is more evident in smaller surveys such as CEERS and the PEARLS NEP and GOODS-S fields.

preprint2022arXiv

The Impact of Dust on the Sizes of Galaxies in the Epoch of Reionization

We study the sizes of galaxies in the Epoch of Reionization using a sample of ~100,000 galaxies from the BlueTides cosmological hydrodynamical simulation from z=7 to 11. We measure the galaxy sizes from stellar mass and luminosity maps, defining the effective radius as the minimum radius which could enclose the pixels containing 50% of the total mass/light in the image. We find an inverse relationship between stellar mass and effective half-mass radius, suggesting that the most massive galaxies are more compact and dense than lower mass galaxies, which have flatter mass distributions. We find a mildly negative relation between intrinsic far-ultraviolet luminosity and size, while we find a positive size-luminosity relation when measured from dust-attenuated images. This suggests that dust is the predominant cause of the observed positive size-luminosity relation, with dust preferentially attenuating bright sight lines resulting in a flatter emission profile and thus larger measured effective radii. We study the size-luminosity relation across the rest-frame ultraviolet and optical, and find that the slope decreases at longer wavelengths; this is a consequence of the relation being caused by dust, which produces less attenuation at longer wavelengths. We find that the far-ultraviolet size-luminosity relation shows mild evolution from z=7 to 11, and galaxy size evolves with redshift as $R\propto(1+z)^{-m}$, where $m=0.662\pm0.009$. Finally, we investigate the sizes of z=7 quasar host galaxies, and find that while the intrinsic sizes of quasar hosts are small relative to the overall galaxy sample, they have comparable sizes when measured from dust-attenuated images.

preprint2022arXiv

The Star Formation Burstiness and Ionizing Efficiency of Low-mass Galaxies

We investigate the burstiness of star formation and the ionizing efficiency of a large sample of galaxies at $0.7 < z < 1.5$ using HST grism spectroscopy and deep ultraviolet (UV) imaging in the GOODS-N and GOODS-S fields. The star formation history (SFH) in these strong emission line low-mass galaxies indicates an elevated star formation rate (SFR) based on the H$α$ emission line at a given stellar mass when compared to the standard main sequence. Moreover, when comparing the H$α$ and UV SFR indicators, we find that an excess in SFR(H$α$) compared to SFR(UV) is preferentially observed in lower-mass galaxies below $10^{9}$ M$\odot$, which are also the highest-EW galaxies. These findings suggest that the burstiness parameters of these strong emission line galaxies may differ from those inferred from hydrodynamical simulations and previous observations. For instance, a larger burstiness duty cycle would explain the observed SFR(H$α$) excess. We also estimate the ionizing photon production efficiency $ξ_{ion}$, finding a median value of Log($ξ_{ion}$/erg$^{-1}$ Hz)$=24.80 \pm 0.26$ when adopting a Galactic dust correction for H$α$ and an SMC one for the stellar component. We observe an increase of $ξ_{ion}$ with redshift, further confirming similar results at higher redshifts. We also find that $ξ_{ion}$ is strongly correlated with EW(H$α$), which provides an approach for deriving $ξ_{ion}$ in early galaxies. Lower-mass, lower-luminosity galaxies have a higher $ξ_{ion}$. Overall, these results provide further support for faint galaxies playing a major role in the reionization of the Universe.

preprint2020arXiv

Cosmic variance of $z>7$ galaxies: Prediction from BlueTides

In the coming decade, a new generation of telescopes, including JWST and WFIRST, will probe the period of the formation of first galaxies and quasars, and open up the last frontier for structure formation. Recent simulations as well as observations have suggested that these galaxies are strongly clustered (with large scale bias $\gtrsim6$), and therefore have significant cosmic variance. In this work, we use \texttt{BlueTides}, the largest volume cosmological simulation of galaxy formation, to directly estimate the cosmic variance for current and upcoming surveys. Given its resolution and volume, \texttt{BlueTides} can probe the bias and cosmic variance of $z>7$ galaxies between magnitude $M_{UV}\sim-16$ to $M_{UV}\sim-22$ over survey areas $\sim0.1\ \mathrm{arcmin}^2$ to $\sim 10~\mathrm{deg}^2$. Within this regime, the cosmic variance decreases with survey area/ volume as a power law with exponents between $\sim-0.25$ to $\sim-0.45$. For the planned $10~\mathrm{deg}^2$ field of WFIRST, the cosmic variance is between $3\%$ to $10\%$. Upcoming JWST medium/ deep surveys with areas up to $A\sim100\ \mathrm{arcmin}^2$ will have cosmic variance ranging from $\sim 20-50\%$. Lensed surveys have the highest cosmic variance $\gtrsim 40\%$; the cosmic variance of $M_{UV}\lesssim-16$ galaxies is $\lesssim100\%$ up to $z\sim11$. At higher redshifts such as $z\sim12~(14)$, effective volumes of $\gtrsim(8~\mathrm{Mpc}/h)^3$ ($\gtrsim(12\ \mathrm{Mpc}/h)^3$) are required to limit the cosmic variance to within $100\%$. Finally, we find that cosmic variance is larger than Poisson variance and forms the dominant component of the overall uncertainty in all current and upcoming surveys. We present our calculations in the form of simple fitting functions and an online cosmic variance calculator (CV_AT_COSMIC_DAWN) which we publicly release.

preprint2016arXiv

Galaxy formation in the Planck cosmology II. Star formation histories and post-processing magnitude reconstruction

We adapt the L-Galaxies semi-analytic model to follow the star-formation histories (SFH) of galaxies -- by which we mean a record of the formation time and metallicities of the stars that are present in each galaxy at a given time. We use these to construct stellar spectra in post-processing, which offers large efficiency savings and allows user-defined spectral bands and dust models to be applied to data stored in the Millennium data repository. We contrast model SFHs from the Millennium Simulation with observed ones from the VESPA algorithm as applied to the SDSS-7 catalogue. The overall agreement is good, with both simulated and SDSS galaxies showing a steeper SFH with increased stellar mass. The SFHs of blue and red galaxies, however, show poor agreement between data and simulations, which may indicate that the termination of star formation is too abrupt in the models. The mean star-formation rate (SFR) of model galaxies is well-defined and is accurately modelled by a double power law at all redshifts: SFR proportional to $1/(x^{-1.39}+x^{1.33})$, where $x=(t_a-t)/3.0\,$Gyr, $t$ is the age of the stars and $t_a$ is the loopback time to the onset of galaxy formation; above a redshift of unity, this is well approximated by a gamma function: SFR proportional to $x^{1.5}e^{-x}$, where $x=(t_a-t)/2.0\,$Gyr. Individual galaxies, however, show a wide dispersion about this mean. When split by mass, the SFR peaks earlier for high-mass galaxies than for lower-mass ones, and we interpret this downsizing as a mass-dependence in the evolution of the quenched fraction: the SFHs of star-forming galaxies show only a weak mass dependence.

preprint2016arXiv

Monsters in the Dark: Predictions for Luminous Galaxies in the Early Universe from the BlueTides Simulation

Using deep Hubble and Spitzer observations Oesch et al. (2016) have identified a bright ($M_{\rm UV}\approx -22$) star forming galaxy candidate at $z \approx 11$. The presence of GN-$z11$ implies a number density $\sim 10^{-6}\,{\rm Mpc^{-3}}$, roughly an order of magnitude higher than the expected value based on extrapolations from lower redshift. Using the unprecedented volume and high resolution of the BlueTides cosmological hydrodynamical simulation, we study the population of luminous rare objects at $z > 10$. The luminosity function in BlueTides implies an enhanced number of massive galaxies, consistent with the observation of GN-$z11$. We find about 30 galaxies at $M_{\rm UV}\approx -22$ at $z = 11$ in the BlueTides volume, including a few objects about 1.5 magnitudes brighter. The probability of observing GN-$z11$ in the volume probed by Oesch et al. (2016) is $\sim 13$ per cent. The predicted properties of the rare bright galaxies at $z = 11$ in BlueTides closely match those inferred from the observations of GN-$z11$. BlueTides predicts a negligible contribution from faint AGN in the observed SED. The enormous increase in volume surveyed by WFIRST will provide observations of $\sim1000$ galaxies with $M_{\rm UV} < -22$ beyond $z = 11$ out to $z = 13.5$.

preprint2015arXiv

Galaxy formation in the PLANCK cosmology - III. The high-redshift universe

We present high-redshift predictions of the star-formation-rate distribution function (SFR DF), UV luminosity function (UV LF), galactic stellar mass function (GSMF), and specific star-formation rates (sSFRs) of galaxies from the latest version of the Munich semi-analytic model L-Galaxies. We find a good fit to both the shape and normalisation of the SFR DF at $z=4-7$, apart from a slight under-prediction at the low SFR end at $z=4$. Likewise, we find a good fit to the faint number counts for the observed UV LF; at brighter magnitudes our predictions lie below the observations, increasingly so at higher redshifts. At all redshifts and magnitudes, the raw (unattenuated) number counts for the UV LF lie above the observations. Because of the good agreement with the SFR we interpret our under-prediction as an over-estimate of the amount of dust in the model for the brightest galaxies, especially at high-redshift. While the shape of our GSMF matches that of the observations, we lie between (conflicting) observations at $z=4-5$, and under-predict at $z=6-7$. The sSFRs of our model galaxies show the observed trend of increasing normalisation with redshift, but do not reproduce the observed mass dependence. Overall, we conclude that the latest version of L-Galaxies, which is tuned to match observations at $z\leq3$, does a fair job of reproducing the observed properties of galaxies at $z\geq4$. More work needs to be done on understanding observational bias at high-redshift, and upon the dust model, before strong conclusions can be drawn on how to interpret remaining discrepancies between the model and observations.

preprint2015arXiv

The BlueTides Simulation: First Galaxies and Reionization

We introduce the BlueTides simulation and report initial results for the luminosity functions of the first galaxies and AGN, and their contribution to reionization. BlueTides was run on the BlueWaters cluster at NCSA from $z=99$ to $z=8.0$ and includes 2$\times$7040$^3$ particles in a $400$Mpc/h per side box, making it the largest hydrodynamic simulation ever performed at high redshift. BlueTides includes a pressure-entropy formulation of smoothed particle hydrodynamics, gas cooling, star formation (including molecular hydrogen), black hole growth and models for stellar and AGN feedback processes. The star formation rate density in the simulation is a good match to current observational data at $z\sim 8-10$. We find good agreement between observations and the predicted galaxy luminosity function in the currently observable range $-18\le M_{\mathrm UV} \le -22.5$ with some dust extinction required to match the abundance of brighter objects. BlueTides implements a patchy reionization model that produces a fluctuating UV background. BlueTides predicts number counts for galaxies fainter than current observational limits which are consistent with extrapolating the faint end slope of the luminosity function with a power law index $α\sim -1.8$ at $z\sim 8$ and redshift dependence of $α\sim (1+z)^{-0.4}$. The AGN population has a luminosity function well fit by a power law with a slope $α\sim -2.4$ that compares favourably with the deepest CANDELS-Goods fields. We investigate how these luminosity functions affect the progress of reionization, and find that a high Lyman-$α$ escape fraction ($f_\mathrm{esc} \sim 0.5$) is required if galaxies dominate the ionising photon budget during reionization. Smaller galaxy escape fractions imply a large contribution from faint AGN (down to $M_\mathrm{UV}=-12$) which results in a rapid reionization, disfavoured by current observations.

preprint2015arXiv

The Case for a James Webb Space Telescope Extragalactic Key Project

The upcoming launch of the James Webb Space Telescope (JWST) in less than three years is certain to bring a revolution in our understanding of many area of astrophysics, with one of the key goals being galaxy evolution. As the first proposals will be due in a little over two years, the time is ripe to take a holistic look at the science goals which the community would wish to accomplish with this observatory. Contrary to our experiences with the Hubble Space Telescope, which has now operated successfully for over two decades due to several timely servicing missions, the lifetime of JWST is finite and relatively short, with a lifetime requirement of five years, and a ten-year goal. Following the discussion session at the &#34;Exploring the Universe with JWST&#34; conference at ESA-ESTEC in October 2015, we highlight in this document the (non-local) extragalactic science goals for JWST. We describe how a concerted community effort could best address these, ensuring that the desired survey can be completed during the JWST mission.

preprint2015arXiv

The formation of Milky Way-mass disk galaxies in the first 500 million years of a cold dark matter universe

Whether among the myriad tiny proto-galaxies there exists a population with similarities to present day galaxies is an open question. We show, using BlueTides, the first hydrodynamic simulation large enough to resolve the relevant scales, that the first massive galaxies to form are %in fact predicted to have extensive rotationally-supported disks. Although their morphology resembles in some ways Milky-way types seen at much lower redshifts, these high-redshift galaxies are smaller, denser, and richer in gas than their low redshift counterparts. From a kinematic analysis of a statistical sample of 216 galaxies at redshift $z=8-10$ we have found that disk galaxies make up 70\% of the population of galaxies with stellar mass $10^{10} M_\odot$ or greater. Cold Dark Matter cosmology therefore makes specific predictions for the population of large galaxies 500 million years after the Big Bang. We argue that wide-field satellite telescopes (e.g. WFIRST) will in the near future discover these first massive disk galaxies. The simplicity of their structure and formation history should make possible new tests of cosmology.

preprint2014arXiv

Asking gender questions: Results from a survey of gender and question asking among UK Astronomers at NAM2014

We report on a survey of astronomers asking questions at the most recent National Astronomy Meeting (NAM2014). The gender balance of both speakers and session chairs at NAM (31% and 29% women respectively) closely matched that of attendees (28% female). However, we find that women were under-represented among question askers (just 18% female). Women were especially under-represented in asking the first question (only 14% of first questions asked by women), but when the Q&A session reached four or more questions, women and men were observed to ask roughly equal numbers of questions. We found a small, but statistically insignificant, increase in the fraction of questions from women in sessions where the chair was also female. On average $2.2\pm0.1$ questions were asked per talk, with no detectable difference in the number of questions asked of female and male speakers, but on average female chairs solicited slightly fewer questions than male chairs. We compare these results to a similar study by Davenport et al. (2014) for the AAS, who also found under-representation of women amongst question askers, but saw more pronounced gender effects when a session chair or speaker were female. We conclude with suggestions for improving the balance of questions at future astronomy meetings.

preprint2011arXiv

Candidate z~8-9 Galaxies from WFC3 Imaging

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

preprint2009arXiv

A public catalogue of stellar masses, star formation and metallicity histories and dust content from the Sloan Digital Sky Survey using VESPA

We applied the VESPA algorithm to the Sloan Digital Sky Survey final data release of the Main Galaxies and Luminous Red Galaxies samples. The result is a catalogue of stellar masses, detailed star formation and metallicity histories and dust content of nearly 800,000 galaxies. We make the catalogue public via a T-SQL database, which is described in detail in this paper. We present the results using a range of stellar population and dust models, and will continue to update the catalogue as new and improved models are made public. The data and documentation are currently online, and can be found at http://www-wfau.roe.ac.uk/vespa/. We also present a brief exploration of the catalogue, and show that the quantities derived are robust: luminous red galaxies can be described by one to three populations, whereas a main galaxy sample galaxy needs on average two to five; red galaxies are older and less dusty; the dust values we recover are well correlated with measured Balmer decrements and star formation rates are also in agreement with previous measurements.

preprint2009arXiv

The Contribution of High Redshift Galaxies to Cosmic Reionization: New Results from Deep WFC3 Imaging of the Hubble Ultra Deep Field

We have searched for star-forming galaxies at z~7 by applying the Lyman-break technique to newly-released 1.1micron Y-band images from WFC3 on HST. By comparing these images of the Hubble Ultra Deep Field with the ACS z&#39;-band (0.85micron), we identify objects with red colours, (z&#39;-Y)_AB>1.3), consistent with the Ly-alpha forest absorption at z~6.7-8.8. We identify 12 of these z&#39;-drops down to a limiting magnitude Y_AB<28.5 (equivalent to a star formation rate of 1.3M_sun/yr at z=7.1), which are undetected in the other ACS filters. We use the WFC3 J-band image to eliminate contaminant low mass Galactic stars, which typically have redder colours than z~7 galaxies. One of our z&#39;-drops is a probably a T-dwarf star. The z~7 z&#39;-drops have much bluer spectral slopes than Lyman-break galaxies at lower redshift. Our brightest z&#39;-drop is not present in the NICMOS J-band image of the same field taken 5 years before, and is a possible transient object. From the 10 remaining z~7 candidates we determine a lower limit on the star formation rate density of 0.0017M_sun/yr/Mpc^3 for a Salpeter initial mass function, which rises to 0.0025-0.0034M_sun/yr/Mpc^3 after correction for luminosity bias. The star formation rate density is a factor of ~10 less than that at z=3-4, and is about half the value at z~6. While based on a single deep field, our results suggest that this star formation rate density would produce insufficient Lyman continuum photons to reionize the Universe unless the escape fraction of these photons is extremely high (f_esc>0.5), and the clumping factor of the Universe is low. Even then, we need to invoke a large contribution from galaxies below our detection limit. The apparent shortfall in ionizing photons might be alleviated if stellar populations at high redshift are low metallicity or have a top-heavy IMF.