sort results by

Use logical operators AND, OR, NOT and round brackets to construct complex queries. Whitespace-separated words are treated as ANDed.

Show articles per page in mode

Wilkins, Stephen M.

Normalized to: Wilkins, S.

79 article(s) in total. 318 co-authors, from 1 to 19 common article(s). Median position in authors list is 5,0.

[1]  [pdf] - 2070440
Nebular Line Emission During the Epoch of Reionization
Comments: 17 pages, accepted to MNRAS, significant modification from v1.0 data available at
Submitted: 2019-04-16, last modified: 2020-03-26
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.
[2]  [pdf] - 2014471
The host galaxies of z=7 quasars: predictions from the BlueTides simulation
Comments: 15 pages, 18 figures. Submitted to MNRAS
Submitted: 2019-12-06, last modified: 2019-12-12
We examine the properties of the host galaxies of $z=7$ quasars using the large volume, cosmological hydrodynamical simulation BlueTides. We find that the most massive black holes and quasars are hosted by massive galaxies with stellar masses $\log(M_\ast/M_\odot)=10.8\pm0.2$, and $10.2\pm0.4$, which have a wide range of star formation rates, of $50\substack{+119 \\ -35}M_\odot/\rm{yr}$ and $19\substack{+28 \\ -12}M_\odot/\rm{yr}$, respectively. The hosts of the most massive black holes and quasars in BlueTides are generally bulge-dominated, with bulge-to-total mass ratio $B/T\simeq0.85\pm0.1$, however their morphologies are not biased relative to the overall $z=7$ galaxy sample. We find that the hosts of the most massive black holes and quasars are significantly more compact, with half-mass radii $R_{0.5}=0.29\substack{+0.15 \\ -0.10}$ kpc and $0.28\substack{+0.08 \\ -0.06}$ kpc respectively, relative to galaxies with similar masses, which have $R_{0.5}=0.5\substack{+0.3 \\ -0.2}$ kpc. We make mock James Webb Space Telescope (JWST) images of these quasars and their host galaxies. We find that distinguishing the host from the quasar emission will be possible but still challenging with JWST, due to the small sizes of quasar hosts. We find that currently observable quasar samples are biased tracers of the intrinsic black hole--stellar mass relations, following a relation that is 0.27 dex higher than that of the full galaxy sample. Finally, we find that black hole hosts are generally quite isolated. However, the most massive black holes are more likely to be found in denser environments than the typical $M_{\textrm{BH}}>10^{6.5}M_\odot$ black hole, indicating that minor mergers at least play some role in growing black holes in the early Universe.
[3]  [pdf] - 1984015
Recalibrating the Cosmic Star Formation History
Comments: 7 pages, 8 figures, accepted for publication in MNRAS
Submitted: 2019-10-11
The calibrations linking observed luminosities to the star formation rate depend on the assumed stellar population synthesis model, initial mass function, star formation and metal enrichment history, and whether reprocessing by dust and gas is included. Consequently the shape and normalisation of the inferred cosmic star formation history is sensitive to these assumptions. Using v2.2.1 of the Binary Population and Spectral Synthesis (\bpass) model we determine a new set of calibration coefficients for the ultraviolet, thermal-infrared, and, hydrogen recombination lines. These ultraviolet and thermal infrared coefficients are 0.15-0.2 dex higher than those widely utilised in the literature while the H$\alpha$ coefficient is $\sim 0.35$ dex larger. These differences arise in part due to the inclusion binary evolution pathways but predominantly reflect an extension in the IMF to 300 $M_{\odot}$ and a change in the choice of reference metallicity. We use these new coefficients to recalibrate the cosmic star formation history, and find improved agreement between the integrated cosmic star formation history and the in-situ measured stellar mass density as a function of redshift. However, these coefficients produce new tension between star formation rate densities inferred from the ultraviolet and thermal-infrared and those from H$\alpha$.
[4]  [pdf] - 1983830
Learning the Relationship between Galaxies Spectra and their Star Formation Histories using Convolutional Neural Networks and Cosmological Simulations
Comments: 20 pages, 21 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society (MNRAS)
Submitted: 2019-03-25, last modified: 2019-10-09
We present a new method for inferring galaxy star formation histories (SFH) using machine learning methods coupled with two cosmological hydrodynamic simulations. We train Convolutional Neural Networks to learn the relationship between synthetic galaxy spectra and high resolution SFHs from the EAGLE and Illustris models. To evaluate our SFH reconstruction we use Symmetric Mean Absolute Percentage Error (SMAPE), which acts as a true percentage error in the low-error regime. On dust-attenuated spectra we achieve high test accuracy (median SMAPE $= 10.5\%$). Including the effects of simulated observational noise increases the error ($12.5\%$), however this is alleviated by including multiple realisations of the noise, which increases the training set size and reduces overfitting ($10.9\%$). We also make estimates for the observational and modelling errors. To further evaluate the generalisation properties we apply models trained on one simulation to spectra from the other, which leads to only a small increase in the error (median SMAPE $\sim 15\%$). We apply each trained model to SDSS DR7 spectra, and find smoother histories than in the VESPA catalogue. This new approach complements the results of existing SED fitting techniques, providing star formation histories directly motivated by the results of the latest cosmological simulations.
[5]  [pdf] - 1931529
Cosmic variance of $z>7$ galaxies: Prediction from BlueTides
Submitted: 2019-08-07
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 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, BlueTides can probe the bias and cosmic variance of $z>7$ galaxies between magnitude $H\sim30$ to $\sim25$ 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 $6\%$ to $10\%$. Upcoming JWST surveys with areas up to $A\sim100~\mathrm{arcmin}^2$ will have cosmic variance ranging from $\sim 20-40\%$. Lensed surveys have the highest cosmic variance $\gtrsim 60\%$; these surveys require volumes $\gtrsim(6~\mathrm{Mpc}/h)^3$ and $\gtrsim(11~\mathrm{Mpc}/h)^3$ at $z\sim7.5$ and $z\sim9$ respectively to limit the cosmic variance to less than $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.
[6]  [pdf] - 1920876
Detailed dust modelling in the L-Galaxies semi-analytic model of galaxy formation
Comments: 19 pages, 16 figures. Edits: Corrected equation 15 and its parameters. Accepted in MNRAS
Submitted: 2019-04-03, last modified: 2019-07-05
We implement a detailed dust model into the L-Galaxies semi-analytical model which includes: injection of dust by type II and type Ia supernovae (SNe) and AGB stars; grain growth in molecular clouds; and destruction due to supernova-induced shocks, star formation, and reheating. Our grain growth model follows the dust content in molecular clouds and the inter-cloud medium separately, and allows growth only on pre-existing dust grains. At early times, this can make a significant difference to the dust growth rate. Above $z\sim8$, type II SNe are the primary source of dust, whereas below $z\sim8$, grain growth in molecular clouds dominates, with the total dust content being dominated by the latter below $z\sim6$. However, the detailed history of galaxy formation is important for determining the dust content of any individual galaxy. We introduce a fit to the dust-to-metal (DTM) ratio as a function of metallicity and age, which can be used to deduce the DTM ratio of galaxies at any redshift. At $z\lesssim3$, we find a fairly flat mean relation between metallicity and the DTM, and a positive correlation between metallicity and the dust-to-gas (DTG) ratio, in good agreement with the shape and normalisation of the observed relations. We also match the normalisation of the observed stellar mass -- dust mass relation over the redshift range of $0-4$, and to the dust mass function at $z=0$. Our results are important in interpreting observations on the dust content of galaxies across cosmic time, particularly so at high redshift.
[7]  [pdf] - 1850019
On the observability of individual Population III stars and their stellar-mass black hole accretion disks through cluster caustic transits
Comments: 6 pages + references, 3 figures, science white paper submitted to the Astro2020 decadal survey. arXiv admin note: substantial text overlap with arXiv:1801.03584
Submitted: 2019-03-13
Recent near-IR power-spectra and panchromatic Extragalactic Background Light measurements provide upper limits on the near-IR surface brightness (SB>31 mag/arcsec^2) that may come from Pop III stars and accretion disks around resulting stellar-mass black holes (BHs) in the epoch of First Light (z=7-17). Physical parameters for zero metallicity Pop III stars at z>7 can be estimated from MESA stellar evolution models through helium-depletion, and for BH accretion disks from quasar microlensing results and multicolor accretion models. Second-generation stars can form at higher multiplicity, so that BH accretion disks may be fed by Roche-lobe overflow from lower-mass companions in their AGB stage. The near-IR SB constraints can be used to calculate the number of caustic transits behind lensing clusters that JWST and the 25~39 m ground-based telescopes may detect for both Pop III stars and stellar mass BH accretion disks. Because Pop III stars and stellar mass BH accretion disks have sizes of a few x 10^-11 arcsec at z>7, typical caustic magnifications can be mu=10^4~10^5, with rise times of hours and decline times of < 1 year for cluster transverse velocities of v<1000 km/s. Microlensing by intracluster medium objects can modify transit magnifications, and lengthen visibility times. Depending on BH masses, accretion-disk radii and feeding efficiencies, stellar-mass BH accretion-disk caustic transits could outnumber those from Pop III stars. To observe Pop III caustic transits directly may require monitoring 3~30 lensing clusters to AB< 29 mag over a decade or more. Such a program must be started with JWST in Cycle 1, and -- depending on the role of microlensing in the Intra Cluster Light -- should be continued for decades with the GMT and TMT, where JWST and the ground-based telescopes each will play a unique and strongly complementary role.
[8]  [pdf] - 1721160
HDUV: The Hubble Deep UV Legacy Survey
Comments: 11 pages, 10 figures, accepted for publication in ApJS
Submitted: 2018-06-05
We present the Hubble Deep UV Legacy Survey (HDUV), a 132 orbit imaging program with the WFC3/UVIS camera onboard the Hubble Space Telescope (HST). The HDUV extends and builds on the few previous UV imaging surveys in the two GOODS/CANDELS-Deep fields to provide deep images over a total area of ~100 arcmin2 in the two filters F275W and F336W. Our release also includes all the F275W imaging data taken by the CANDELS survey, which were aligned using a novel approach and combined with the HDUV survey data. By reaching depths of 27.5-28.0 mag (5sigma, in 0.4" apertures), these are the deepest high-resolution UV data over such a large area taken to date. Such unique UV imaging enables a wide range of science by the community. Among the main goals of the HDUV survey are: (1) provide a complete sample of faint star-forming galaxies at z~1-3, (2) constrain the ionizing photon escape fraction from galaxies at z~2-3, and (3) track the build-up of bulges and the disappearance of clumpy disk galaxies through reliable internal stellar population properties at sub-kpc resolution out to z~3. The addition of the HDUV data further enhances the legacy value of the two GOODS/CANDELS-Deep fields, which now include deep 11-band HST imaging as well as very deep ancillary data from X-ray to radio, enabling unique multi-wavelength studies. Here, we provide an overview of the survey design, describe the data reduction, and highlight a few basic analyses on the images which are released to the community as high level science products via the Mikulski Archive for Space Telescopes (MAST).
[9]  [pdf] - 1790741
A tiny host galaxy for the first giant black hole: $z= 7.5$ quasar in BlueTides
Comments: 12 pages, 14 figures, submitted to MNRAS
Submitted: 2018-06-01
The most distant known quasar recently discovered by Ba\~nados et al. (2018) is at $z=7.5$ (690 Myr after the Big Bang), at the dawn of galaxy formation. We explore the host galaxy of the brightest quasar in the large volume cosmological hydrodynamic simulation BlueTides, which in Phase II has reached these redshifts. The brightest quasar in BlueTides has a luminosity of a $\sim$ few $10^{13} L_{\odot}$ and a black hole mass of $6.4 \times 10^{8} M_{\odot}$ at $z \sim 7.5$, comparable to the observed quasar (the only one in this large volume). The quasar resides in a rare halo of mass $M_{H} \sim 10^{12} M_{\odot}$ and has a host galaxy of stellar mass of $4 \times 10^{10}M_{\odot}$ with an ongoing (intrinsic) star formation rate of $\sim 80 M_{\odot} yr^{-1}$. The corresponding intrinsic UV magnitude of the galaxy is $-23.1$, which is roughly $2.7$ magnitudes fainter than the quasar's magnitude of $-25.9$. We find that the galaxy is highly metal enriched with a mean metallicity equal to the solar value. We derive quasar and galaxy spectral energy distribution (SED) in the mid and near infrared JWST bands. We predict a significant amount of dust attenuation in the rest-frame UV corresponding to $A_{1500} \sim 1.7$ giving an UV based SFR of $\sim 14 M_{\odot} yr^{-1}$. We present mock JWST images of the galaxy with and without central point source, in different MIRI and NIRCam filters. The host galaxy is detectable in NIRCam filters, but it is extremely compact ($R_{E}=0.35$ kpc). It will require JWST's exquisite sensitivity and resolution to separate the galaxy from the central point source. Finally within the FOV of the quasar in BlueTides there are two more sources that would be detectable by JWST.
[10]  [pdf] - 1795700
Substellar and low-mass dwarf identification with near-infrared imaging space observatories
Comments: Accepted for publication in A&A
Submitted: 2018-05-02
AIMS: We aim to evaluate the near-infrared colors of brown dwarfs as observed with four major infrared imaging space observatories: the Hubble Space Telescope (HST), the James Webb Space Telescope (JWST), the Euclid mission, and the WFIRST telescope. METHODS: We used the SPLAT SPEX/ISPEX spectroscopic library to map out the colors of the M-, L-, and T-type dwarfs. We have identified which color-color combination is optimal for identifying broad type and which single color is optimal to then identify the subtype (e.g., T0-9). We evaluated each observatory separately as well as the narrow-field (HST and JWST) and wide-field (Euclid and WFIRST) combinations. RESULTS: The Euclid filters perform poorly typing brown dwarfs and WFIRST performs only marginally better, despite a wider selection of filters. WFIRST's W146 and F062 combined with Euclid's Y-band discriminates somewhat better between broad brown dwarf categories. However, subtyping with any combination of Euclid and WFIRST observations remains uncertain due to the lack of medium or narrow-band filters. We argue that a medium band added to the WFIRST filter selection would greatly improve its ability to preselect brown dwarfs its imaging surveys. CONCLUSIONS: The HST filters used in high-redshift searches are close to optimal to identify broad stellar type. However, the addition of F127M to the commonly used broad filter sets would allow for unambiguous subtyping. An improvement over HST is one of two broad and medium filter combinations on JWST: pairing F140M with either F150W or F162M discriminates very well between subtypes.
[11]  [pdf] - 1659667
Red, redder, reddest: SCUBA-2 imaging of colour-selected \textit{Herschel} sources
Comments: 22 pages, 14 figures and 3 tables
Submitted: 2018-01-22
High-redshift, luminous, dusty star forming galaxies (DSFGs) constrain the extremity of galaxy formation theories. The most extreme are discovered through follow-up on candidates in large area surveys. Here we present 850 $\mu$m SCUBA-2 follow-up observations of 188 red DSFG candidates from the \textit{Herschel} Multi-tiered Extragalactic Survey (HerMES) Large Mode Survey, covering 274 deg$^2$. We detected 87 per cent with a signal-to-noise ratio $>$ 3 at 850~$\mu$m. We introduce a new method for incorporating the confusion noise in our spectral energy distribution fitting by sampling correlated flux density fluctuations from a confusion limited map. The new 850~$\mu$m data provide a better constraint on the photometric redshifts of the candidates, with photometric redshift errors decreasing from $\sigma_z/(1+z)\approx0.21$ to $0.15$. Comparison spectroscopic redshifts also found little bias ($\langle (z-z_{\rm spec})/(1+z_{\rm spec})\rangle = 0.08 $). The mean photometric redshift is found to be 3.6 with a dispersion of $0.4$ and we identify 21 DSFGs with a high probability of lying at $z > 4$. After simulating our selection effects we find number counts are consistent with phenomenological galaxy evolution models. There is a statistically significant excess of WISE-1 and SDSS sources near our red galaxies, giving a strong indication that lensing may explain some of the apparently extreme objects. Nevertheless, our sample should include examples of galaxies with the highest star formation rates in the Universe ($\gg10^3$ M$_\odot$yr$^{-1}$).
[12]  [pdf] - 1601625
Characterising and Identifying Galaxy Protoclusters
Comments: 19 pages, 16 figures, submitted to MNRAS
Submitted: 2017-10-05, last modified: 2017-12-01
We study the characteristics of galaxy protoclusters using the latest L-galaxies semi-analytic model. Searching for protoclusters on a scale of $\sim 10 \, \mathrm{cMpc}$ gives an excellent compromise between the completeness and purity of their galaxy populations, leads to high distinction from the field in overdensity space, and allows accurate determination of the descendant cluster mass. This scale is valid over a range of redshifts and selection criteria. We present a procedure for estimating, given a measured galaxy overdensity, the protocluster probability and its descendant cluster mass for a range of modelling assumptions, particularly taking into account the shape of the measurement aperture. This procedure produces lower protocluster probabilities compared to previous estimates using fixed size apertures. The relationship between AGN and protoclusters is also investigated, and shows significant evolution with redshift; at $z \sim 2$ the fraction of protoclusters traced by AGN is high, but the fraction of all AGN in protoclusters is low, whereas at $z \geqslant 5$ the fraction of protoclusters containing AGN is low, but most AGN are in protoclusters. We also find indirect evidence for the emergence of a passive sequence in protoclusters at $z \sim 2$, and note that a significant fraction of all galaxies reside in protoclusters at $z \geqslant 2$, particularly the most massive.
[13]  [pdf] - 1613976
Galaxy And Mass Assembly (GAMA): the G02 field, Herschel-ATLAS target selection and Data Release 3
Comments: 14 pages, 15 figures, accepted by MNRAS. GAMA DR3 is available at
Submitted: 2017-11-24
We describe data release 3 (DR3) of the Galaxy And Mass Assembly (GAMA) survey. The GAMA survey is a spectroscopic redshift and multi-wavelength photometric survey in three equatorial regions each of 60.0 deg^2 (G09, G12, G15), and two southern regions of 55.7 deg^2 (G02) and 50.6 deg^2 (G23). DR3 consists of: the first release of data covering the G02 region and of data on H-ATLAS sources in the equatorial regions; and updates to data on sources released in DR2. DR3 includes 154809 sources with secure redshifts across four regions. A subset of the G02 region is 95.5% redshift complete to r<19.8 over an area of 19.5 deg^2, with 20086 galaxy redshifts, that overlaps substantially with the XXL survey (X-ray) and VIPERS (redshift survey). In the equatorial regions, the main survey has even higher completeness (98.5%), and spectra for about 75% of H-ATLAS filler targets were also obtained. This filler sample extends spectroscopic redshifts, for probable optical counterparts to H-ATLAS sub-mm sources, to 0.8 mag deeper (r<20.6) than the GAMA main survey. There are 25814 galaxy redshifts for H-ATLAS sources from the GAMA main or filler surveys. GAMA DR3 is available at the survey website (
[14]  [pdf] - 1637549
GAMA/G10-COSMOS/3D-HST: The 0<z<5 cosmic star-formation history, stellar- and dust-mass densities
Comments: Accepted for publication in MNRAS
Submitted: 2017-10-18, last modified: 2017-10-19
We use the energy-balance code MAGPHYS to determine stellar and dust masses, and dust corrected star-formation rates for over 200,000 GAMA galaxies, 170,000 G10-COSMOS galaxies and 200,000 3D-HST galaxies. Our values agree well with previously reported measurements and constitute a representative and homogeneous dataset spanning a broad range in stellar mass (10^8---10^12 Msol), dust mass (10^6---10^9 Msol), and star-formation rates (0.01---100 Msol per yr), and over a broad redshift range (0.0 < z < 5.0). We combine these data to measure the cosmic star-formation history (CSFH), the stellar-mass density (SMD), and the dust-mass density (DMD) over a 12 Gyr timeline. The data mostly agree with previous estimates, where they exist, and provide a quasi-homogeneous dataset using consistent mass and star-formation estimators with consistent underlying assumptions over the full time range. As a consequence our formal errors are significantly reduced when compared to the historic literature. Integrating our cosmic star-formation history we precisely reproduce the stellar-mass density with an ISM replenishment factor of 0.50 +/- 0.07, consistent with our choice of Chabrier IMF plus some modest amount of stripped stellar mass. Exploring the cosmic dust density evolution, we find a gradual increase in dust density with lookback time. We build a simple phenomenological model from the CSFH to account for the dust mass evolution, and infer two key conclusions: (1) For every unit of stellar mass which is formed 0.0065---0.004 units of dust mass is also formed; (2) Over the history of the Universe approximately 90 to 95 per cent of all dust formed has been destroyed and/or ejected.
[15]  [pdf] - 1604961
Dust Obscured Star Forming Galaxies in the Early Universe
Comments: 8 pages, accepted for publication in MNRAS
Submitted: 2017-10-05
Motivated by recent observational constraints on dust reprocessed emission in star forming galaxies at $z\sim 6$ and above we use the very-large cosmological hydrodynamical simulation \bluetides\ to explore predictions for the amount of dust obscured star formation in the early Universe ($z>8$). \bluetides\ matches current observational constraints on both the UV luminosity function and galaxy stellar mass function and predicts that approximately $90\%$ of the star formation in high-mass ($M_{*}>10^{10}\,{\rm M_{\odot}}$) galaxies at $z=8$ is already obscured by dust. The relationship between dust attenuation and stellar mass predicted by \bluetides\ is consistent with that observed at lower redshift. However, observations of several individual objects at $z>6$ are discrepant with the predictions, though it is possible their uncertainties may have been underestimated. We find that the predicted surface density of $z\ge 8$ sub-mm sources is below that accessible to current {\em Herschel}, SCUBA-2, and ALMA sub-mm surveys. However, as ALMA continues to accrue additional surface area the population of $z>8$ dust-obscured galaxies may become accessible in the near future.
[16]  [pdf] - 1630140
The HDUV Survey: A Revised Assessment of the Relationship between UV Slope and Dust Attenuation for High-Redshift Galaxies
Comments: 22 pages, 3 tables, 11 figures, submitted to the Astrophysical Journal
Submitted: 2017-05-25
We use a newly assembled large sample of 3,545 star-forming galaxies with secure spectroscopic, grism, and photometric redshifts at z=1.5-2.5 to constrain the relationship between UV slope (beta) and dust attenuation (L(IR)/L(UV)=IRX). Our sample benefits from the combination of deep Hubble WFC3/UVIS photometry from the Hubble Deep UV (HDUV) Legacy survey and existing photometric data compiled in the 3D-HST survey, and extends the range of UV luminosity and beta probed in previous UV-selected samples. IRX is measured using stacks of deep Herschel/PACS 100 and 160 micron data, and the results are compared with predictions of the IRX-beta relation for different assumptions of the stellar population model and obscuration curve. We find that z=1.5-2.5 galaxies have an IRX-beta relation that is consistent with the predictions for an SMC extinction curve if we invoke sub-solar metallicity models that are currently favored for high-redshift galaxies, while the commonly assumed starburst attenuation curve over-predicts the IRX at a given beta by a factor of ~3. The IRX of high-mass (M*>10^9.75 Msun) galaxies is a factor of >4 larger than that of low-mass galaxies, lending support for the use of stellar mass as a proxy for attenuation. The commonly observed trend of fainter galaxies having bluer beta may simply reflect bluer intrinsic UV slopes for such galaxies, rather than lower obscurations. The IRX-beta for young/low-mass galaxies implies a dust curve that is steeper than the SMC, suggesting a lower attenuation at a given beta relative to older/more massive galaxies. The lower attenuations and higher ionizing photon output implied by low metallicity stellar population models point to Lyman continuum production efficiencies, xi_ion, that may be elevated by a factor of ~2 relative to the canonical value for L* galaxies, aiding in their ability to keep the universe ionized at z~2. [Abridged]
[17]  [pdf] - 1583685
Galaxy And Mass Assembly: the evolution of the cosmic spectral energy distribution from z = 1 to z = 0
Comments: 19 pages, 14 figures, accepted by MNRAS
Submitted: 2017-05-22
We present the evolution of the Cosmic Spectral Energy Distribution (CSED) from $z = 1 - 0$. Our CSEDs originate from stacking individual spectral energy distribution fits based on panchromatic photometry from the Galaxy and Mass Assembly (GAMA) and COSMOS datasets in ten redshift intervals with completeness corrections applied. Below $z = 0.45$, we have credible SED fits from 100 nm to 1 mm. Due to the relatively low sensitivity of the far-infrared data, our far-infrared CSEDs contain a mix of predicted and measured fluxes above $z = 0.45$. Our results include appropriate errors to highlight the impact of these corrections. We show that the bolometric energy output of the Universe has declined by a factor of roughly four -- from $5.1 \pm 1.0$ at $z \sim 1$ to $1.3 \pm 0.3 \times 10^{35}~h_{70}$~W~Mpc$^{-3}$ at the current epoch. We show that this decrease is robust to cosmic variance, SED modelling and other various types of error. Our CSEDs are also consistent with an increase in the mean age of stellar populations. We also show that dust attenuation has decreased over the same period, with the photon escape fraction at 150~nm increasing from $16 \pm 3$ at $z \sim 1$ to $24 \pm 5$ per cent at the current epoch, equivalent to a decrease in $A_\mathrm{FUV}$ of 0.4~mag. Our CSEDs account for $68 \pm 12$ and $61 \pm 13$ per cent of the cosmic optical and infrared backgrounds respectively as defined from integrated galaxy counts and are consistent with previous estimates of the cosmic infrared background with redshift.
[18]  [pdf] - 1583328
Galaxy And Mass Assembly (GAMA): The galaxy stellar mass function to $z=0.1$ from the r-band selected equatorial regions
Comments: 22 pages; 15 figures; Accepted for publication in MNRAS
Submitted: 2017-05-11
We derive the low redshift galaxy stellar mass function (GSMF), inclusive of dust corrections, for the equatorial Galaxy And Mass Assembly (GAMA) dataset covering 180 deg$^2$. We construct the mass function using a density-corrected maximum volume method, using masses corrected for the impact of optically thick and thin dust. We explore the galactic bivariate brightness plane ($M_\star-\mu$), demonstrating that surface brightness effects do not systematically bias our mass function measurement above 10$^{7.5}$ M$_{\odot}$. The galaxy distribution in the $M-\mu$-plane appears well bounded, indicating that no substantial population of massive but diffuse or highly compact galaxies are systematically missed due to the GAMA selection criteria. The GSMF is {fit with} a double Schechter function, with $\mathcal M^\star=10^{10.78\pm0.01\pm0.20}M_\odot$, $\phi^\star_1=(2.93\pm0.40)\times10^{-3}h_{70}^3$Mpc$^{-3}$, $\alpha_1=-0.62\pm0.03\pm0.15$, $\phi^\star_2=(0.63\pm0.10)\times10^{-3}h_{70}^3$Mpc$^{-3}$, and $\alpha_2=-1.50\pm0.01\pm0.15$. We find the equivalent faint end slope as previously estimated using the GAMA-I sample, although we find a higher value of $\mathcal M^\star$. Using the full GAMA-II sample, we are able to fit the mass function to masses as low as $10^{7.5}$ $M_\odot$, and assess limits to $10^{6.5}$ $M_\odot$. Combining GAMA-II with data from G10-COSMOS we are able to comment qualitatively on the shape of the GSMF down to masses as low as $10^{6}$ $M_\odot$. Beyond the well known upturn seen in the GSMF at $10^{9.5}$ the distribution appears to maintain a single power-law slope from $10^9$ to $10^{6.5}$. We calculate the stellar mass density parameter given our best-estimate GSMF, finding $\Omega_\star= 1.66^{+0.24}_{-0.23}\pm0.97 h^{-1}_{70} \times 10^{-3}$, inclusive of random and systematic uncertainties.
[19]  [pdf] - 1582334
The properties of the first galaxies in the BLUETIDES simulation
Comments: 16 pages, accepted for publication in MNRAS
Submitted: 2017-04-04
We employ the very large cosmological hydrodynamical simulation BLUETIDES to investigate the predicted properties of the galaxy population during the epoch of reionisation ($z>8$). BLUETIDES has a resolution and volume ($(400/h\approx 577)^{3}\,{\rm cMpc^3}$) providing a population of galaxies which is well matched to depth and area of current observational surveys targeting the high-redshift Universe. At $z=8$ BLUETIDES includes almost 160,000 galaxies with stellar masses $>10^{8}\,{\rm M_{\odot}}$. The population of galaxies predicted by BLUETIDES closely matches observational constraints on both the galaxy stellar mass function and far-UV ($150\,{\rm nm}$) luminosity function. Galaxies in BLUETIDES are characterised by rapidly increasing star formation histories. Specific star formation rates decrease with redshift though remain largely insensitive to stellar mass. As a result of the enhanced surface density of metals more massive galaxies are predicted to have higher dust attenuation resulting in a significant steepening of the observed far-UV luminosity function at high luminosities. The contribution of active SMBHs to the UV luminosities of galaxies with stellar masses $10^{9-10}\,{\rm M_{\odot}}$ is around $3\%$ on average. Approximately $25\%$ of galaxies with $M_{*}\approx 10^{10}\,{\rm M_{\odot}}$ are predicted to have active SMBH which contribute $>10\%$ of the total UV luminosity.
[20]  [pdf] - 1534860
Galaxy And Mass Assembly: The 1.4GHz SFR indicator, SFR-M* relation and predictions for ASKAP-GAMA
Comments: 14 pages, 6 figures
Submitted: 2017-01-22
We present a robust calibration of the 1.4GHz radio continuum star formation rate (SFR) using a combination of the Galaxy And Mass Assembly (GAMA) survey and the Faint Images of the Radio Sky at Twenty-cm (FIRST) survey. We identify individually detected 1.4GHz GAMA-FIRST sources and use a late-type, non-AGN, volume-limited sample from GAMA to produce stellar mass-selected samples. The latter are then combined to produce FIRST-stacked images. This extends the robust parametrisation of the 1.4GHz-SFR relation to faint luminosities. For both the individually detected galaxies and our stacked samples, we compare 1.4GHz luminosity to SFRs derived from GAMA to determine a new 1.4GHz luminosity-to-SFR relation with well constrained slope and normalisation. For the first time, we produce the radio SFR-M* relation over 2 decades in stellar mass, and find that our new calibration is robust, and produces a SFR-M* relation which is consistent with all other GAMA SFR methods. Finally, using our new 1.4GHz luminosity-to-SFR calibration we make predictions for the number of star-forming GAMA sources which are likely to be detected in the upcoming ASKAP surveys, EMU and DINGO.
[21]  [pdf] - 1580667
The HDUV Survey: Six Lyman Continuum Emitter Candidates at z~2 Revealed by HST UV Imaging
Comments: 11 pages, 6 figures, submitted to ApJ
Submitted: 2016-11-21
We present six galaxies at z~2 that show evidence of Lyman continuum (LyC) emission based on the newly acquired UV imaging of the Hubble Deep UV legacy survey (HDUV) conducted with the WFC3/UVIS camera on the Hubble Space Telescope (HST). At the redshift of these sources, the HDUV F275W images partially probe the ionizing continuum. By exploiting the HST multi-wavelength data available in the HDUV/GOODS fields, models of the UV spectral energy distributions, and detailed Monte Carlo simulations of the intergalactic medium absorption, we estimate the absolute ionizing photon escape fractions of these galaxies to be very high -- typically >60% (>13% for all sources at 90% likelihood). Our findings are in broad agreement with previous studies that found only a small fraction of galaxies to show high escape fraction. These six galaxies comprise the largest sample yet of LyC leaking candidates at z~2 whose inferred LyC flux has been cleanly observed at HST resolution. While three of our six candidates show evidence of hosting an active galactic nucleus (AGN), two of these are heavily obscured and their LyC emission appears to originate from star-forming regions rather than the central nucleus. This suggests an AGN-aided pathway for LyC escape from these sources. Extensive multi-wavelength data in the GOODS fields, especially the near-IR grism spectra from the 3D-HST survey, enable us to study the candidates in detail and tentatively test some recently proposed indirect methods to probe LyC leakage -- namely, the [OIII]/[OII] line ratio and the H$\beta-$UV slope diagram. High-resolution spectroscopic followup of our candidates will help constrain such indirect methods which are our only hope of studying $f_{esc}$ at z~5-9 in the fast-approaching era of the James Webb Space Telescope.
[22]  [pdf] - 1530379
Inferred H{\alpha} 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
Comments: 22 pages, 13 figures, 9 tables, ApJ, in press
Submitted: 2015-11-27, last modified: 2016-11-01
We derive H{\alpha} 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{\alpha} flux is inferred based on the offset between the IRAC 3.6 {\mu}m flux and that predicted from the best-fit SED. We demonstrate that the H{\alpha} 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{\alpha}/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{\alpha}/SFR ratio). Using H{\alpha} 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{\alpha}-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.
[23]  [pdf] - 1488348
Galaxy formation in the Planck cosmology II. Star formation histories and post-processing magnitude reconstruction
Comments: Accepted version of the paper, to appear in MNRAS. Compared to the original version, contains more detail on the post-processing of magnitudes, including a table of rms magnitude errors. SFHs available on Millennium database
Submitted: 2015-01-22, last modified: 2016-09-28
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.
[24]  [pdf] - 1530849
ALMA Spectroscopic Survey in the Hubble Ultra Deep Field: The Infrared Excess of UV-selected z=2-10 galaxies as a function of UV-continuum Slope and Stellar Mass
Comments: 33 pages, 19 figures, 14 tables, accepted for publication in ApJ, updated to match version in press, paper VI of the ASPECS (ALMA Spectroscopic Survey in the Hubble Ultra Deep Field) series
Submitted: 2016-06-16, last modified: 2016-08-06
We make use of deep 1.2mm-continuum observations (12.7microJy/beam RMS) of a 1 arcmin^2 region in the Hubble Ultra Deep Field to probe dust-enshrouded star formation from 330 Lyman-break galaxies spanning the redshift range z=2-10 (to ~2-3 Msol/yr at 1sigma over the entire range). Given the depth and area of ASPECS, we would expect to tentatively detect 35 galaxies extrapolating the Meurer z~0 IRX-beta relation to z>~2 (assuming T_d~35 K). However, only 6 tentative detections are found at z>~2 in ASPECS, with just three at >3sigma. Subdividing z=2-10 galaxies according to stellar mass, UV luminosity, and UV-continuum slope and stacking the results, we only find a significant detection in the most massive (>10^9.75 Msol) subsample, with an infrared excess (IRX=L_{IR}/L_{UV}) consistent with previous z~2 results. However, the infrared excess we measure from our large selection of sub-L* (<10^9.75 Msol) galaxies is 0.11(-0.42)(+0.32) and 0.14(-0.14)(+0.15) at z=2-3 and z=4-10, respectively, lying below even an SMC IRX-beta relation (95% confidence). These results demonstrate the relevance of stellar mass for predicting the IR luminosity of z>~2 galaxies. We furthermore find that the evolution of the IRX-stellar mass relationship depends on the evolution of the dust temperature. If the dust temperature increases monotonically with redshift (as (1+z)^0.32) such that T_d~44-50 K at z>=4, current results are suggestive of little evolution in this relationship to z~6. We use these results to revisit recent estimates of the z>~3 SFR density. One less obvious implication is in interpreting the high Halpha EWs seen in z~5 galaxies: our results imply that star-forming galaxies produce Lyman-continuum photons at twice the efficiency (per unit UV luminosity) as implied in conventional models. Star-forming galaxies can then reionize the Universe, even if the escape fraction is <10%.
[25]  [pdf] - 1553853
The origin of most massive black holes at high-z: BLUETIDES and the next quasar frontier
Comments: 9 pages, 7 figures, submitted to MNRAS
Submitted: 2016-06-28
The growth of the most massive black holes in the early universe, consistent with the detection of highly luminous quasars at $z> 6$ implies sustained, critical accretion of material to grow and power them. Given a black hole seed scenario, it is still uncertain which conditions in the early Universe allow the fastest black hole growth. Large scale hydrodynamical cosmological simulations of structure formation allow us to explore the conditions conducive to the growth of the earliest supermassive black holes. We use the cosmological hydrodynamic simulation BlueTides, which incorporates a variety of baryon physics in a (400 Mpc/h)^3 volume with 0.7 trillion particles to follow the earliest phases of black hole critical growth. At z=8 the most massive black holes (a handful) approach masses of 10^8 Msun with the most massive (with M_BH = 4 x 10^8 Msun ) being found in an extremely compact spheroid-dominated host galaxy. Examining the large-scale environment of hosts, we find that the initial tidal field is more important than overdensity in setting the conditions for early BH growth. In regions of low tidal fields gas accretes 'cold' onto the black hole and falls along thin, radial filaments straight into the center forming the most compact galaxies and most massive black holes at earliest times. Regions of high tidal fields instead induce larger coherent angular momenta and influence the formation of the first population of massive compact disks. The extreme early growth depends on the early interplay of high gas densities and the tidal field that shapes the mode of accretion. Mergers play a minor role in the formation of the first generation, rare massive BHs.
[26]  [pdf] - 1470666
Forecasts for the WFIRST High Latitude Survey using the BlueTides Simulation
Comments: 12 pages, 17 figures, submitted to MNRAS May 13th, 2016
Submitted: 2016-05-18, last modified: 2016-05-19
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$.
[27]  [pdf] - 1443979
The Photometric Properties of Galaxies in the Early Universe
Comments: 11 pages, accepted in MNRAS
Submitted: 2016-05-17
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$.
[28]  [pdf] - 1396995
Galaxy And Mass Assembly (GAMA): Accurate Panchromatic Photometry from Optical Priors using LAMBDAR
Comments: 38 pages, 27 figures, 3 tables, 3 appendices, accepted for publication in MNRAS
Submitted: 2016-04-07, last modified: 2016-04-08
We present the Lambda Adaptive Multi-Band Deblending Algorithm in R (LAMBDAR), a novel code for calculating matched aperture photometry across images that are neither pixel- nor PSF-matched, using prior aperture definitions derived from high resolution optical imaging. The development of this program is motivated by the desire for consistent photometry and uncertainties across large ranges of photometric imaging, for use in calculating spectral energy distributions. We describe the program, specifically key features required for robust determination of panchromatic photometry: propagation of apertures to images with arbitrary resolution, local background estimation, aperture normalisation, uncertainty determination and propagation, and object deblending. Using simulated images, we demonstrate that the program is able to recover accurate photometric measurements in both high-resolution, low-confusion, and low-resolution, high-confusion, regimes. We apply the program to the 21-band photometric dataset from the Galaxy And Mass Assembly (GAMA) Panchromatic Data Release (PDR; Driver et al. 2016), which contains imaging spanning the far-UV to the far-IR. We compare photometry derived from LAMBDAR with that presented in Driver et al. (2016), finding broad agreement between the datasets. Nonetheless, we demonstrate that the photometry from LAMBDAR is superior to that from the GAMA PDR, as determined by a reduction in the outlier rate and intrinsic scatter of colours in the LAMBDAR dataset. We similarly find a decrease in the outlier rate of stellar masses and star formation rates using LAMBDAR photometry. Finally, we note an exceptional increase in the number of UV and mid-IR sources able to be constrained, which is accompanied by a significant increase in the mid-IR colour-colour parameter-space able to be explored.
[29]  [pdf] - 1436338
Monsters in the Dark: Predictions for Luminous Galaxies in the Early Universe from the BlueTides Simulation
Comments: 5 pages, 7 figures
Submitted: 2016-04-01, last modified: 2016-04-05
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$.
[30]  [pdf] - 1426144
Z > 7 galaxies with red Spitzer/IRAC [3.6]-[4.5] colors in the full CANDELS data set: the brightest-known galaxies at Z ~ 7-9 and a probable spectroscopic confirmation at Z=7.48
Comments: 19 pages, 13 figures, 5 tables, accepted for publication in the Astrophysical Journal. Updated to match the version in press and reflect spectroscopic confirmation of the highest redshift candidate in our sample by Zitrin+2015
Submitted: 2015-06-02, last modified: 2016-04-01
We identify 4 unusually bright (H < 25.5) galaxies from HST and Spitzer CANDELS data with probable redshifts z ~ 7-9. These identifications include the brightest-known galaxies to date at z > 7.5. As Y-band observations are not available over the full CANDELS program to perform a standard Lyman-break selection of z > 7 galaxies, we employ an alternate strategy using deep Spitzer/IRAC data. We identify z ~ 7.1 - 9.1 galaxies by selecting z >~ 6 galaxies from the HST CANDELS data that show quite red IRAC [3.6]-[4.5] colors, indicating strong [OIII]+Hbeta lines in the 4.5 micron band. This selection strategy was validated using a modest sample for which we have deep Y-band coverage, and subsequently used to select the brightest z > 7 sources. Applying the IRAC criteria to all HST-selected optical-dropout galaxies over the full ~900 arcmin**2 of the CANDELS survey revealed four unusually bright z ~ 7.1, 7.6, 7.9 and 8.6 candidates. The median [3.6]-[4.5] color of our selected z ~ 7.1-9.1 sample is consistent with rest-frame [OIII]+Hbeta EWs of ~1500A, in the [4.5] band. Keck/MOSFIRE spectroscopy has been independently reported for two of our selected sources, showing Ly-alpha at redshifts of 7.7302+/-0.0006 and 8.683^+0.001_-0.004, respectively. We present similar Keck/MOSFIRE spectroscopy for a third selected galaxy with a probable 4.7sigma Ly-alpha line at z_spec=7.4770+/-0.0008. All three have H-band magnitudes of ~25 mag and are ~0.5 mag more luminous (M(UV) ~ -22.0) than any previously discovered z ~ 8 galaxy, with important implications for the UV LF. Our 3 brightest, highest redshift z > 7 galaxies all lie within the CANDELS EGS field, providing a dramatic illustration of the potential impact of field-to-field variance.
[31]  [pdf] - 1373170
The Lyman-continuum photon production efficiency in the high-redshift Universe
Comments: 4 pages, 2 figures, accepted to MNRAS. Minor changes reflecting referee comments and community feedback
Submitted: 2015-12-10, last modified: 2016-03-14
The Lyman Continuum photon production efficiency ($\xi_{\rm ion}$) is a critical ingredient for inferring the number of photons available to reionise the intergalactic medium. To estimate the theoretical production efficiency in the high-redshift Universe we couple the BlueTides cosmological hydrodynamical simulation with a range of stellar population synthesis models. We find Lyman Continuum photon production efficiencies of $\log_{10}(\xi_{\rm ion}/{\rm erg^{-1}\, Hz})\approx 25.1-25.5$ depending on the choice of stellar population synthesis model. These results are broadly consistent with recent observational constraints at high-redshift though favour a model incorporating the effects of binary evolution
[32]  [pdf] - 1327624
The Case for a James Webb Space Telescope Extragalactic Key Project
Comments: White paper summarizing and developing the results of a discussion session at the "Exploring the Universe with JWST" meeting at ESA-ESTEC in October 2015. Four pages, one figure, one table
Submitted: 2015-12-14
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 "Exploring the Universe with JWST" 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.
[33]  [pdf] - 1323446
Galaxy And Mass Assembly (GAMA): Panchromatic Data Release (far-UV --- far-IR) and the low-z energy budget
Comments: 31 pages and 30 figures. Accepted in MNRAS. High-resolution copy available from our release site: or directly via
Submitted: 2015-08-09, last modified: 2015-11-20
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$\mu$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:
[34]  [pdf] - 1314288
Two-phase galaxy evolution: the cosmic star-formation histories of spheroids and discs
Comments: Accepted for publication in MNRAS, 12 pages, 6 figures, 3 animated gifs (available via MNRAS online or as indicated in text)
Submitted: 2013-01-06, last modified: 2015-11-20
From two very simple axioms: (1) that AGN activity traces spheroid formation, and (2) that the cosmic star-formation history is dominated by spheroid formation at high redshift, we derive simple expressions for the star-formation histories of spheroids and discs, and their implied metal enrichment histories. Adopting a Baldry-Glazebrook initial mass function we use these relations and apply PEGASE.2 to predict the z=0 cosmic spectral energy distributions (CSEDs) of spheroids and discs. The model predictions compare favourably to the dust-corrected CSED recently reported by the Galaxy And Mass Assembly (GAMA) team from the FUV through to the K band. The model also provides a reasonable fit to the total stellar mass contained within spheroid and disc structures as recently reported by the Millennium Galaxy Catalogue team. Three interesting inferences can be made following our axioms: (1) there is a transition redshift at z ~ 1.7 at which point the Universe switches from what we refer to as "hot mode evolution" (i.e., spheroid formation/growth via mergers and/or collapse) to what we term "cold mode evolution" (i.e., disc formation/growth via gas infall and minor mergers); (2) there is little or no need for any pre-enrichment prior to the main phase of star-formation; (3) in the present Universe mass-loss is fairly evenly balanced with star-formation holding the integrated stellar mass density close to a constant value. The model provides a simple prediction of the energy output from spheroid and disc projenitors, the build-up of spheroid and disc mass, and the mean metallicity enrichment of the Universe.
[35]  [pdf] - 1308351
Quantifying the UV continuum slopes of galaxies to z~10 using deep Hubble and Spitzer/IRAC observations
Comments: 10 pages, 10 figures, accepted in MNRAS
Submitted: 2015-10-06
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.
[36]  [pdf] - 1288607
Exploring the Dawn of Galaxies
Comments: 4 pages, introduction to the "Cosmic Dawn" series of articles published in Astronomy & Astrophysics
Submitted: 2015-10-06
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.
[37]  [pdf] - 1288344
Reionization after Planck: The Derived Growth of the Cosmic Ionizing Emissivity now matches the Growth of the Galaxy UV Luminosity Density
Comments: 20 pages, 12 figures, 5 tables, Astrophysical Journal, updated to match version in press, Figure 6 shows the main result of the paper
Submitted: 2015-03-27, last modified: 2015-09-05
Thomson optical depth tau measurements from Planck provide new insights into the reionization of the universe. In pursuit of model-independent constraints on the properties of the ionising sources, we determine the empirical evolution of the cosmic ionizing emissivity. We use a simple two-parameter model to map out the evolution in the emissivity at z>~6 from the new Planck optical depth tau measurements, from the constraints provided by quasar absorption spectra and from the prevalence of Ly-alpha emission in z~7-8 galaxies. We find the redshift evolution in the emissivity dot{N}_{ion}(z) required by the observations to be d(log Nion)/dz=-0.15(-0.11)(+0.08), largely independent of the assumed clumping factor C_{HII} and entirely independent of the nature of the ionising sources. The trend in dot{N}_{ion}(z) is well-matched by the evolution of the galaxy UV-luminosity density (dlog_{10} rho_UV/dz=-0.11+/-0.04) to a magnitude limit >~-13 mag, suggesting that galaxies are the sources that drive the reionization of the universe. The role of galaxies is further strengthened by the conversion from the UV luminosity density rho_UV to dot(N)_{ion}(z) being possible for physically-plausible values of the escape fraction f_{esc}, the Lyman-continuum photon production efficiency xi_{ion}, and faint-end cut-off $M_{lim}$ to the luminosity function. Quasars/AGN appear to match neither the redshift evolution nor normalization of the ionizing emissivity. Based on the inferred evolution in the ionizing emissivity, we estimate that the z~10 UV-luminosity density is 8(-4)(+15)x lower than at $z~6, consistent with the observations. The present approach of contrasting the inferred evolution of the ionizing emissivity with that of the galaxy UV luminosity density adds to the growing observational evidence that faint, star-forming galaxies drive the reionization of the universe.
[38]  [pdf] - 1284963
Luminosity function of [OII] emission-line galaxies in the MassiveBlack-II simulation
Comments: 12 pages, 13 figures, 3 tables, accepted for publication in MNRAS
Submitted: 2015-08-20
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.
[39]  [pdf] - 1247613
Galaxy And Mass Assembly (GAMA): the effect of close interactions on star formation in galaxies
Comments: 22 pages, 17 figure, Accepted MNRAS
Submitted: 2015-07-16
The modification of star formation (SF) in galaxy interactions is a complex process, with SF observed to be both enhanced in major mergers and suppressed in minor pair interactions. Such changes likely to arise on short timescales and be directly related to the galaxy-galaxy interaction time. Here we investigate the link between dynamical phase and direct measures of SF on different timescales for pair galaxies, targeting numerous star-formation rate (SFR) indicators and comparing to pair separation, individual galaxy mass and pair mass ratio. We split our sample into the higher (primary) and lower (secondary) mass galaxies in each pair and find that SF is indeed enhanced in all primary galaxies but suppressed in secondaries of minor mergers. We find that changes in SF of primaries is consistent in both major and minor mergers, suggesting that SF in the more massive galaxy is agnostic to pair mass ratio. We also find that SF is enhanced/suppressed more strongly for short-time duration SFR indicators (e.g. H-alpha), highlighting recent changes to SF in these galaxies, which are likely to be induced by the interaction. We propose a scenario where the lower mass galaxy has its SF suppressed by gas heating or stripping, while the higher mass galaxy has its SF enhanced, potentially by tidal gas turbulence and shocks. This is consistent with the seemingly contradictory observations for both SF suppression and enhancement in close pairs.
[40]  [pdf] - 1238888
Galaxy formation in the PLANCK cosmology - III. The high-redshift universe
Comments: Updated to reflect the published version in MNRAS (12 pages, 7 figures). Note that this was previously paper IV, but was renamed to III by the journal. Some references have been fixed, and the Bowler2014b observational errors in Figure 3 have been corrected to their proper values
Submitted: 2015-04-13, last modified: 2015-06-29
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.
[41]  [pdf] - 1259091
Galaxy And Mass Assembly (GAMA): end of survey report and data release 2
Comments: Accepted for publication in MMRAS, 40 pages, 33 figures
Submitted: 2015-06-26
The Galaxy And Mass Assembly (GAMA) survey is one of the largest contemporary spectroscopic surveys of low-redshift galaxies. Covering an area of ~286 deg^2 (split among five survey regions) down to a limiting magnitude of r < 19.8 mag, we have collected spectra and reliable redshifts for 238,000 objects using the AAOmega spectrograph on the Anglo-Australian Telescope. In addition, we have assembled imaging data from a number of independent surveys in order to generate photometry spanning the wavelength range 1 nm - 1 m. Here we report on the recently completed spectroscopic survey and present a series of diagnostics to assess its final state and the quality of the redshift data. We also describe a number of survey aspects and procedures, or updates thereof, including changes to the input catalogue, redshifting and re-redshifting, and the derivation of ultraviolet, optical and near-infrared photometry. Finally, we present the second public release of GAMA data. In this release we provide input catalogue and targeting information, spectra, redshifts, ultraviolet, optical and near-infrared photometry, single-component S\'ersic fits, stellar masses, H$\alpha$-derived star formation rates, environment information, and group properties for all galaxies with r < 19.0 mag in two of our survey regions, and for all galaxies with r < 19.4 mag in a third region (72,225 objects in total). The database serving these data is available at
[42]  [pdf] - 1258844
The formation of Milky Way-mass disk galaxies in the first 500 million years of a cold dark matter universe
Comments: Updated Figure 3, fixing bug in V/sigma calculation
Submitted: 2015-04-24, last modified: 2015-06-15
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.
[43]  [pdf] - 1319542
The BlueTides Simulation: First Galaxies and Reionization
Submitted: 2015-04-24
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 $\alpha\sim -1.8$ at $z\sim 8$ and redshift dependence of $\alpha\sim (1+z)^{-0.4}$. The AGN population has a luminosity function well fit by a power law with a slope $\alpha\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-$\alpha$ 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.
[44]  [pdf] - 908100
Asking gender questions: Results from a survey of gender and question asking among UK Astronomers at NAM2014
Comments: 15 pages, 7 figures, comments welcome
Submitted: 2014-12-15
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.
[45]  [pdf] - 1216899
Galaxy And Mass Assembly (GAMA): Curation and reanalysis of 16.6k redshifts in the G10/COSMOS region
Comments: 16 pages, 15 figures, Accepted to MNRAS
Submitted: 2014-09-11, last modified: 2014-12-08
We discuss the construction of the Galaxy And Mass Assembly (GAMA) 10h region (G10) using publicly available data in the Cosmic Evolution Survey region (COSMOS) in order to extend the GAMA survey to z~1 in a single ~1 deg^2 field. In order to obtain the maximum number of high precision spectroscopic redshifts we re-reduce all archival zCOSMOS-bright data and use the GAMA automatic cross-correlation redshift fitting code autoz. We use all available redshift information (autoz, zCOSMOS-bright 10k, PRIMUS, VVDS, SDSS and photometric redshifts) to calculate robust best-fit redshifts for all galaxies and visually inspect all 1D and 2D spectra to obtain 16,583 robust redshifts in the full COSMOS region. We then define the G10 region to be the central ~1deg^2 of COSMOS, which has relatively high spectroscopic completeness, and encompasses the CHILES VLA region. We define a combined r < 23.0mag & i < 22.0mag G10 sample (selected to have the highest bijective overlap) with which to perform future analysis, containing 9,861 sources with reliable high precision VLT-VIMOS spectra. All tables, spectra and imaging are available at: .
[46]  [pdf] - 1451820
Galaxy And Mass Assembly (GAMA): mass-size relations of z$<$0.1 galaxies subdivided by S\'ersic index, colour and morphology
Submitted: 2014-11-24
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\'ersic 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\'ersic light profile fitting (using GALFIT3), i.e., elliptical semi-major half light radii, improving on earlier low-z benchmarks based on circular apertures.
[47]  [pdf] - 1223085
Galaxy And Mass Assembly (GAMA): Bivariate functions of H$\alpha$ star forming galaxies
Comments: 27 pages, 21 figures, accepted for publication in MNRAS
Submitted: 2014-11-10
We present bivariate luminosity and stellar mass functions of H$\alpha$ star forming galaxies drawn from the Galaxy And Mass Assembly (GAMA) survey. While optically deep spectroscopic observations of GAMA over a wide sky area enable the detection of a large number of $0.001<{SFR}_{H\alpha}$ (M$_{\odot}$ yr$^{-1}$)$<100$ galaxies, the requirement for an H$\alpha$ detection in targets selected from an $r$-band magnitude limited survey leads to an incompleteness due to missing optically faint star forming galaxies. Using $z<0.1$ bivariate distributions as a reference we model the higher-$z$ distributions, thereby approximating a correction for the missing optically faint star forming galaxies to the local SFR and stellar mass densities. Furthermore, we obtain the $r$-band LFs and stellar mass functions of H$\alpha$ star forming galaxies from the bivariate LFs. As our sample is selected on the basis of detected H$\alpha$ emission, a direct tracer of on-going star formation, this sample represents a true star forming galaxy sample, and is drawn from both photometrically classified blue and red sub-populations, though mostly from the blue population. On average 20-30% of red galaxies at all stellar masses are star forming, implying that these galaxies may be dusty star forming systems.
[48]  [pdf] - 1216367
Interpreting high [O III]/Hbeta ratios with maturing starbursts
Comments: 8 pages. MNRAS accepted
Submitted: 2014-08-18
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.
[49]  [pdf] - 1216158
Galaxy And Mass Assembly (GAMA): Galaxy close-pairs, mergers, and the future fate of stellar mass
Comments: Accepted for publication in MNRAS: 25 Pages, 21 Figures, 3 Tables
Submitted: 2014-08-07
We use a highly complete subset of the GAMA-II redshift sample to fully describe the stellar mass dependence of close-pairs and mergers between 10^8 Msun and 10^12 Msun. Using the analytic form of this fit we investigate the total stellar mass accreting onto more massive galaxies across all mass ratios. Depending on how conservatively we select our robust merging systems, the fraction of mass merging onto more massive companions is 2.0%-5.6%. Using the GAMA-II data we see no significant evidence for a change in the close-pair fraction between redshift $z = 0.05-0.2$. However, we find a systematically higher fraction of galaxies in similar mass close-pairs compared to published results over a similar redshift baseline. Using a compendium of data and the function $\gamma_M =A(1+z)m$ to predict the major close-pair fraction, we find fitting parameters of $A = 0.021 \pm 0.001$ and $m = 1.53 \pm 0.08$, which represents a higher low-redshift normalisation and shallower power-law slope than recent literature values. We find that the relative importance of in-situ star-formation versus galaxy merging is inversely correlated, with star-formation dominating the addition of stellar material below Mstar and merger accretion events dominating beyond Mstar. We find mergers have a measurable impact on the whole extent of the GSMF, manifest as a deepening of the dip in the GSMF over the next Gyr and an increase in Mstar by as much as 0.01-0.05 dex.
[50]  [pdf] - 1180409
Spectroscopy of z ~ 7 candidate galaxies: Using Lyman-alpha to constrain the neutral fraction of hydrogen in the high-redshift universe
Comments: accepted for publication in MNRAS
Submitted: 2013-10-31, last modified: 2014-07-03
Following our previous spectroscopic observations of $z>7$ galaxies with Gemini/GNIRS and VLT/XSHOOTER, which targeted a total of 8 objects, we present here our results from a deeper and larger VLT/FORS2 spectroscopic sample of Wide Field Camera 3 selected $z>7$ candidate galaxies. With our FORS2 setup we cover the 737-1070nm wavelength range, enabling a search for Lyman-$\alpha$ in the redshift range spanning 5.06 - 7.80. We target 22 $z$-band dropouts and find no evidence of Lyman-$\alpha$ emission, with the exception of a tentative detection ($<5\sigma$, which is our adopted criterion for a secure detection) for one object. The upper limits on Lyman-$\alpha$ flux and the broad-band magnitudes are used to constrain the rest-frame Equivalent Widths for this line emission. We analyse our FORS2 observations in combination with our previous GNIRS and XSHOOTER observations, and suggest that a simple model where the fraction of high rest-frame Equivalent Width emitters follows the trend seen at $z=3-6.5$ is inconsistent with our non-detections at $z\sim7.8$ at the 96\% confidence level. This may indicate that a significant neutral HI fraction in the intergalactic medium suppresses Lyman-$\alpha$, with an estimated neutral fraction $\chi_{HI}\sim0.5$, in agreement with other estimates.
[51]  [pdf] - 1203177
The MassiveBlack-II Simulation: The Evolution of Halos and Galaxies to z~0
Comments: 26 pages, 25 figures. Submitted to MNRAS. High-resolution version and MBII galaxy catalogs can be found at
Submitted: 2014-02-04
(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.
[52]  [pdf] - 791931
High-redshift galaxies and low-mass stars
Comments: 14 pages, 14 figures, accepted in MNRAS
Submitted: 2014-01-27
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.
[53]  [pdf] - 1202547
Galaxy And Mass Assembly (GAMA): Mid-Infrared Properties and Empirical Relations from $WISE$
Comments: 20 pages. Accepted to ApJ
Submitted: 2014-01-04
The Galaxy And Mass Assembly (GAMA) survey furnishes a deep redshift catalog that, when combined with the Wide-field Infrared Explorer ($WISE$), allows us to explore for the first time the mid-infrared properties of $> 110, 000$ galaxies over 120 deg$^2$ to $z\simeq 0.5$. In this paper we detail the procedure for producing the matched GAMA-$WISE$ catalog for the G12 and G15 fields, in particular characterising and measuring resolved sources; the complete catalogs for all three GAMA equatorial fields will be made available through the GAMA public releases. The wealth of multiwavelength photometry and optical spectroscopy allows us to explore empirical relations between optically determined stellar mass (derived from synthetic stellar population models) and 3.4micron and 4.6micron WISE measurements. Similarly dust-corrected Halpha-derived star formation rates can be compared to 12micron and 22micron luminosities to quantify correlations that can be applied to large samples to $z<0.5$. To illustrate the applications of these relations, we use the 12micron star formation prescription to investigate the behavior of specific star formation within the GAMA-WISE sample and underscore the ability of WISE to detect star-forming systems at $z\sim0.5$. Within galaxy groups (determined by a sophisticated friends-of-friends scheme), results suggest that galaxies with a neighbor within 100$\,h^{-1} $kpc have, on average, lower specific star formation rates than typical GAMA galaxies with the same stellar mass.
[54]  [pdf] - 1178836
Theoretical predictions for the effect of nebular emission on the broad band photometry of high-redshift galaxies
Comments: 13 pages, 17 figures, accepted in MNRAS
Submitted: 2013-08-28
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.
[55]  [pdf] - 1166007
Single Colour Diagnostics of the Mass-to-light Ratio: Predictions from Galaxy Formation Models
Comments: 11 pages, 7 figures, 2 tables, accepted in MNRAS
Submitted: 2013-04-16
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}\Gamma_{V}=\log_{10}[(M/M_{\odot})/(L_{V}/L_{V\odot})]$) with a precision of $\sigma_{lg\Gamma}\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}\Gamma_{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.
[56]  [pdf] - 1159545
Interpreting the observed UV continuum slopes of high-redshift galaxies
Comments: 6 pages, 4 figures, accepted in MNRAS
Submitted: 2013-02-06
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 ($\beta=-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 $\beta$, combined with the Meurer relation, suggest dust attenuation to be virtually zero in all but the most luminous systems.
[57]  [pdf] - 1159368
Galaxy And Mass Assembly (GAMA): Spectroscopic analysis
Comments: 19 pages, 19 figures, accepted for publication in MNRAS
Submitted: 2013-01-29
The Galaxy And Mass Assembly (GAMA) survey is a multiwavelength photometric and spectroscopic survey, using the AAOmega spectrograph on the Anglo-Australian Telescope to obtain spectra for up to ~300000 galaxies over 280 square degrees, to a limiting magnitude of r_pet < 19.8 mag. The target galaxies are distributed over 0<z<0.5 with a median redshift of z~0.2, although the redshift distribution includes a small number of systems, primarily quasars, at higher redshifts, up to and beyond z=1. The redshift accuracy ranges from sigma_v~50km/s to sigma_v~100km/s depending on the signal-to-noise of the spectrum. Here we describe the GAMA spectroscopic reduction and analysis pipeline. We present the steps involved in taking the raw two-dimensional spectroscopic images through to flux-calibrated one-dimensional spectra. The resulting GAMA spectra cover an observed wavelength range of 3750<lambda<8850 A at a resolution of R~1300. The final flux calibration is typically accurate to 10-20%, although the reliability is worse at the extreme wavelength ends, and poorer in the blue than the red. We present details of the measurement of emission and absorption features in the GAMA spectra. These measurements are characterised through a variety of quality control analyses detailing the robustness and reliability of the measurements. We illustrate the quality of the measurements with a brief exploration of elementary emission line properties of the galaxies in the GAMA sample. We demonstrate the luminosity dependence of the Balmer decrement, consistent with previously published results, and explore further how Balmer decrement varies with galaxy mass and redshift. We also investigate the mass and redshift dependencies of the [NII]/Halpha vs [OIII]/Hbeta spectral diagnostic diagram, commonly used to discriminate between star forming and nuclear activity in galaxies.
[58]  [pdf] - 1159143
VLT/XSHOOTER & Subaru/MOIRCS Spectroscopy of HUDF-YD3: No Evidence for Lyman-alpha Emission at z=8.55
Comments: Accepted for publication in MNRAS
Submitted: 2013-01-18
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.
[59]  [pdf] - 1158988
Confronting predictions of the galaxy stellar mass function with observations at high-redshift
Comments: 7 pages, 4 figures, accepted in MNRAS
Submitted: 2013-01-12
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}(\phi\, [{\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.
[60]  [pdf] - 1151291
The UV colours and dust attenuation of Lyman-break galaxies
Comments: (20 pages, 14 figures, MNRAS)
Submitted: 2012-09-10, last modified: 2012-11-20
Using GALFORM, a semi-analytical model of galaxy formation in the Lambda cold dark matter cosmology, we study the rest-frame ultraviolet (UV) colours of Lyman-break galaxies (LBGs) in the redshift range 2.5 < z < 10. As the impact of dust on UV luminosity can be dramatic, our model includes a self-consistent computation of dust attenuation based on a radiative transfer model. We find that intrinsically brighter galaxies suffer stronger dust attenuation than fainter ones, though the relation has a large scatter. The model predicts galaxies with UV colours consistent with the colour selection regions designed to select LBGs in observational surveys. We find that the drop-out technique that selects LBGs based on two rest-frame UV colours is robust and effective, selecting more than 70 per cent of UV bright galaxies at a given redshift. We investigate the impact on the predicted UV colours of varying selected model parameters. We find that the UV colours are most sensitive to the modelling of dust attenuation and in particular, to the extinction curve used in the radiative transfer calculation. If we assume a Milky Way dust extinction curve, the predicted UV continuum slopes are, in general, bluer than observed. However, we find that the opposite is true when using the Small Magellanic Cloud dust extinction curve. This demonstrates the strong dependence of UV colours on dust properties and highlights the inadequacy of using the UV continuum slope as a tracer of dust attenuation without any further knowledge of the galaxy inclination or dust characteristics in high redshift galaxies.
[61]  [pdf] - 1152617
Constraining the Bright-end of the UV Luminosity Function for z ~ 7 - 9 Galaxies: results from CANDELS/GOODS-South
Comments: Accepted for publication in MNRAS. Originally submitted on 2012 June 18. arXiv admin note: text overlap with arXiv:1006.3545
Submitted: 2012-10-31
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.
[62]  [pdf] - 1151426
The accuracy of the UV continuum as an indicator of the star formation rate in galaxies
Comments: 7 pages, 9 figures, accepted in MNRAS
Submitted: 2012-09-16
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_{\nu}$), 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_{\nu}\times 10^{-28}\times L_{\nu}/[{\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.
[63]  [pdf] - 1151050
Galaxy And Mass Assembly (GAMA): The 0.013 < z < 0.1 cosmic spectral energy distribution from 0.1 micron to 1mm
Comments: Accepted for publication in MNRAS. 23 pages, 17 figures. Figs. 1, 2 & 4 degraded. Based on GAMA data see:
Submitted: 2012-09-03
We use the GAMA I dataset combined with GALEX, SDSS and UKIDSS imaging to construct the low-redshift (z<0.1) galaxy luminosity functions in FUV, NUV, ugriz, and YJHK bands from within a single well constrained volume of 3.4 x 10^5 (Mpc/h)^{3}. The derived luminosity distributions are normalised to the SDSS DR7 main survey to reduce the estimated cosmic variance to the 5 per cent level. The data are used to construct the cosmic spectral energy distribution (CSED) from 0.1 to 2.1 \mum free from any wavelength dependent cosmic variance for both the elliptical and non-elliptical populations. The two populations exhibit dramatically different CSEDs as expected for a predominantly old and young population respectively. Using the Driver et al. (2008) prescription for the azimuthally averaged photon escape fraction, the non-ellipticals are corrected for the impact of dust attenuation and the combined CSED constructed. The final results show that the Universe is currently generating (1.8 +/- 0.3) x 10^{35} h W Mpc^{-3} of which (1.2 +/- 0.1) x 10^{35} h W Mpc^{-3} is directly released into the inter-galactic medium and (0.6 +/- 0.1) x 10^{35} h W Mpc^{-3} is reprocessed and reradiated by dust in the far-IR. Using the GAMA data and our dust model we predict the mid and far-IR emission which agrees remarkably well with available data. We therefore provide a robust description of the pre- and post dust attenuated energy output of the nearby Universe from 0.1micron to 0.6mm. The largest uncertainty in this measurement lies in the mid and far-IR bands stemming from the dust attenuation correction and its currently poorly constrained dependence on environment, stellar mass, and morphology.
[64]  [pdf] - 1150973
No Evidence for Lyman-alpha Emission in Spectroscopy of z > 7 Candidate Galaxies
Comments: Accepted for publication in MNRAS
Submitted: 2012-08-29
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.
[65]  [pdf] - 1124085
Predictions for the intrinsic UV continuum properties of star forming galaxies and the implications for inferring dust extinction
Comments: 9 pages, Accepted to MNRAS
Submitted: 2012-06-13
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.
[66]  [pdf] - 1077645
The UV Properties of Star Forming Galaxies I: HST WFC3 Observations of Very-high Redshift Galaxies
Comments: 14 pages, accepted to MNRAS (27/6/2011)
Submitted: 2011-06-29
The acquisition of deep Near-IR imaging with Wide Field Camera 3 on the Hubble Space Telescope has provided the opportunity to study the very-high redshift Universe. For galaxies up to $z\approx 7.7$ sufficient wavelength coverage exists to probe the rest-frame ultraviolet (UV) continuum without contamination from either Lyman-$\alpha$ emission or the Lyman-$\alpha$ break. In this work we use Near-IR imaging to measure the rest-frame UV continuum colours of galaxies at $4.7<z<7.7$. We are have carefully defined a colour-colour selection to minimise any inherent bias in the measured UV continuum slope for the drop-out samples. For the highest-redshift sample ($6.7<z<7.7$), selected as $z_{f850lp}$-band dropouts, we find mean UV continuum colours approximately equal to zero (AB), consistent with a dust-free, solar metallicity, star forming population (or a moderately dusty population of low metallicity). At lower-redshift we find that the mean UV continuum colours of galaxies (over the same luminosity range) are redder, and that galaxies with higher luminosities are also slightly redder on average. One interpretation of this is that lower-redshift and more luminous galaxies are dustier, however this interpretation is complicated by the effects of the star formation history and metallicity and potentially the initial mass function on the UV continuum colours.
[67]  [pdf] - 1053326
Galaxy And Mass Assembly (GAMA): The star formation rate dependence of the stellar initial mass function
Comments: 18 pages, 15 figures, accepted for publication in MNRAS. Image resolution has been reduced. For the full resolution version, see
Submitted: 2011-04-13
The stellar initial mass function (IMF) describes the distribution in stellar masses produced from a burst of star formation. For more than fifty years, the implicit assumption underpinning most areas of research involving the IMF has been that it is universal, regardless of time and environment. We measure the high-mass IMF slope for a sample of low-to-moderate redshift galaxies from the Galaxy And Mass Assembly survey. The large range in luminosities and galaxy masses of the sample permits the exploration of underlying IMF dependencies. A strong IMF-star formation rate dependency is discovered, which shows that highly star forming galaxies form proportionally more massive stars (they have IMFs with flatter power-law slopes) than galaxies with low star formation rates. This has a significant impact on a wide variety of galaxy evolution studies, all of which rely on assumptions about the slope of the IMF. Our result is supported by, and provides an explanation for, the results of numerous recent explorations suggesting a variation of or evolution in the IMF.
[68]  [pdf] - 1033166
Candidate z~8-9 Galaxies from WFC3 Imaging
Comments: Accepted for publication in MNRAS
Submitted: 2010-06-17, last modified: 2011-02-04
We present a search for galaxies at 7.6<z<9.8 using the latest HST WFC3 near-infrared data, based on the Lyman-break technique. We search for galaxies which have large (Y-J) colours (the "Y-drops") on account of the Lyman-alpha forest absorption, and with (J-H) colours inconsistent with being low-redshift contaminants. We identify 24 candidates at redshift z~8-9 (15 are robust and a further 9 more marginal but consistent with being high redshift) over an area of ~50 square arcminutes. Previous searches for Y-drops with WFC3 have focussed only on the Hubble Ultra Deep Field (HUDF), and our larger survey (involving two other nearby deep fields and a wider area survey) has trebelled the number of robust Y-drop candidates. For the first time, we have sufficient z~8-9 galaxies to fit both phi^* and M^* of the UV Schechter luminosity function. There is evidence for evolution in this luminosity function from z=6-7 to z=8-9, in the sense that there are fewer UV-bright galaxies at z~8-9, consistent with an evolution mainly in M^*. The candidate z~8-9 galaxies we detect have insufficient ionizing flux to reionize the Universe, and it is probable that galaxies below our detection limit provide a significant UV contribution. The faint-end slope, alpha, is not well constrained. However, adopting a similiar faint-end slope to that determined at z=3-6 (alpha=-1.7) and a Salpeter initial mass function, then the ionizing photon budget still falls short if f_escape<0.5, even integrating down to M(UV)=-8. A steeper faint end slope or a low-metallicity population (or a top-heavy IMF) might still provide sufficient photons for star-forming galaxies to reionize the Universe, but confirmation of this might have to await the James Webb Space Telescope.
[69]  [pdf] - 460992
New Star Forming Galaxies at z\approx 7 from WFC3 Imaging
Comments: 15 pages, accepted in MNRAS
Submitted: 2010-02-25, last modified: 2010-10-01
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.
[70]  [pdf] - 1026469
Extending the Canada-France brown Dwarfs Survey to the near-infrared: first ultracool brown dwarfs from CFBDSIR
Comments: A&A, Accepted
Submitted: 2010-04-22, last modified: 2010-08-20
We present the first results of the ongoing Canada-France Brown Dwarfs Survey-InfraRed, hereafter CFBDSIR, a Near InfraRed extension to the optical wide-field survey CFBDS. Our final objectives are to constrain ultracool atmosphere physics by finding a statistically significant sample of objects cooler than 650K and to explore the ultracool brown dwarf mass function building on a well defined sample of such objects. Candidates are identified in CFHT/WIRCam J and CFHT/MegaCam z' images using optimised psf-fitting, and we follow them up with pointed near infrared imaging with SOFI at NTT. We finally obtain low resolution spectroscopy of the coolest candidates to characterise their atmospheric physics. We have so far analysed and followed up all candidates on the first 66 square degrees of the 335 square degrees survey. We identified 55 T-dwarfs candidates with z'-J > 3:5 and have confirmed six of them as T-dwarfs, including 3 that are strong later-than-T8 candidates, based on their far-red and NIR colours. We also present here the NIR spectra of one of these ultracool dwarfs, CFBDSIR1458+1013 which confirms it as one of the coolest brown dwarf known, possibly in the 550-600K temperature range. From the completed survey we expect to discover 10 to 15 dwarfs later than T8, more than doubling the known number of such objects. This will enable detailed studies of their extreme atmospheric properties and provide a stronger statistical base for studies of their luminosity function.
[71]  [pdf] - 31465
A Matched Catalogue of z> 5.9 Galaxies in the WFC3 Hubble Ultra Deep Field
Comments: This is a matched catalogue giving the various different names for the same Lyman-break sources in the HUDF from several independent studies. This is intended for general use to reduce confusion over naming
Submitted: 2009-12-08
There have been several independent analyses of the recent Wide Field Camera 3 images of the Hubble Deep Field, selecting galaxies at z>6 through the Lyman break technique. Presented here is a matched catalogue of objects in common between the analyses posted to this preprint server, listing the different catalogue names associated with the same sources.
[72]  [pdf] - 1017500
The Contribution of High Redshift Galaxies to Cosmic Reionization: New Results from Deep WFC3 Imaging of the Hubble Ultra Deep Field
Comments: To appear in MNRAS. Updated in response to referee report
Submitted: 2009-09-14, last modified: 2009-12-07
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'-band (0.85micron), we identify objects with red colours, (z'-Y)_AB>1.3), consistent with the Ly-alpha forest absorption at z~6.7-8.8. We identify 12 of these z'-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'-drops is a probably a T-dwarf star. The z~7 z'-drops have much bluer spectral slopes than Lyman-break galaxies at lower redshift. Our brightest z'-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.
[73]  [pdf] - 460987
Probing $\sim L_{*}$ Lyman-break Galaxies at $z\approx 7$ in GOODS-South with WFC3 on HST
Comments: Accepted by MNRAS; minor changes to text to reflect referee's comments
Submitted: 2009-10-07, last modified: 2009-12-07
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.
[74]  [pdf] - 900691
A public catalogue of stellar masses, star formation and metallicity histories and dust content from the Sloan Digital Sky Survey using VESPA
Comments: Accepted for publication in ApJS. New and extended sections, references added, some figures changed. Conclusions and catalogue unchanged
Submitted: 2009-04-06, last modified: 2009-09-08
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 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.
[75]  [pdf] - 18141
Detecting Mass Substructure in Galaxy Clusters: An Aperture Mass Statistic for Gravitational Flexion
Comments: 11 pages, 10 figures, MNRAS in press
Submitted: 2008-11-03, last modified: 2009-03-02
Gravitational flexion has recently been introduced as a technique by which one can map out and study substructure in clusters of galaxies. Previous analyses involving flexion have measured the individual galaxy-galaxy flexion signal, or used either parametric techniques or a KSB-type inversion to reconstruct the mass distribution in Abell 1689. In this paper, we present an aperture mass statistic for flexion, and apply it to the lensed images of background galaxies obtained by ray-tracing simulations through a simple analytic mass distribution and through a galaxy cluster from the Millennium simulation. We show that this method is effective at detecting and accurately tracing structure within clusters of galaxies on sub-arcminute scales with high signal-to-noise even using a moderate background source number density and image resolution. In addition, the method provides much more information about both the overall shape and the small-scale structure of a cluster of galaxies than can be achieved through a weak lensing mass reconstruction using gravitational shear data. Lastly, we discuss how the zero-points of the aperture mass might be used to infer the masses of structures identified using this method.
[76]  [pdf] - 16334
Extragalactic Constraints on the Initial Mass Function
Comments: 6 pages, 5 figures, accepted to MNRAS
Submitted: 2008-09-15
The local stellar mass density is observed to be significantly lower than the value obtained from integrating the cosmic star formation history (SFH), assuming that all the stars formed with a Salpeter initial mass function (IMF). Even other favoured IMFs, more successful in reconciling the observed $z=0$ stellar mass density with that inferred from the SFH, have difficulties in reproducing the stellar mass density observed at higher redshift. In this study we investigate to what extent this discrepancy can be alleviated for any universal power-law IMF. We find that an IMF with a high-mass slope shallower (2.15) than the Salpeter slope (2.35) reconciles the observed stellar mass density with the cosmic star formation history, but only at low redshifts. At higher redshifts $z>0.5$ we find that observed stellar mass densities are systematically lower than predicted from the cosmic star formation history, for any universal power-law IMF.
[77]  [pdf] - 11280
The Evolution of Stellar Mass Density and its Implied Star Formation History
Comments: 4 pages, in "Panoramic Views of Galaxy Formation and Evolution" the 1st Subaru international conference, Hayama, Japan, Dec. 2007
Submitted: 2008-03-27
Using a compilation of measurements of the stellar mass density as a function of redshift we can infer the cosmic star formation history. For z < 0.7 there is good agreement between the two star formation histories. At higher redshifts the instantaneous indicators suggest star formation rates larger than that implied by the evolution of the stellar mass density. This discrepancy peaks at z = 3 where instantaneous indicators suggest a star formation rate around 0.6 dex higher than those of the best fit to the stellar mass history. We discuss a variety of explanations for this inconsistency, such as inaccurate dust extinction corrections, incorrect measurements of stellar masses and a possible evolution of the stellar initial mass function.
[78]  [pdf] - 8867
The Evolution of Stellar Mass and the Implied Star Formation History
Comments: 10 pages, 3 figures and 1 table, MNRAS accepted
Submitted: 2008-01-10
We present a compilation of measurements of the stellar mass density as a function of redshift. Using this stellar mass history we obtain a star formation history and compare it to the instantaneous star formation history. For z<0.7 there is good agreement between the two star formation histories. At higher redshifts the instantaneous indicators suggest star formation rates larger than that implied by the evolution of the stellar mass density. This discrepancy peaks at z=3 where instantaneous indicators suggest a star formation rate around 0.6 dex higher than those of the best fit to the stellar mass history. We discuss a variety of explanations for this inconsistency, such as inaccurate dust extinction corrections, incorrect measurements of stellar masses and a possible evolution of the stellar initial mass function.
[79]  [pdf] - 86564
The Galaxy Power Spectrum: 2dFGRS-SDSS tension?
Comments: To appear in Astronomical Society of The Pacific conference series: Cosmic Frontiers
Submitted: 2006-11-06
Published galaxy power spectra from the 2dFGRS and SDSS are not in good agreement. We revisit this issue by analyzing both the 2dFGRS and SDSS DR5 catalogues using essentially identical technqiues. We confirm that the 2dFGRS exhibits relatively more large scale power than the SDSS, or, equivalently, SDSS has more small scale power. We demonstrate that this difference is due the r-band selected SDSS catalogue being dominated by more strongly clustered red galaxies, due to these galaxies having a stronger scale dependent bias. The power spectra of galaxies of the same rest frame colours from the two surveys match well. It is therefore important to accurately model scale dependent bias to get accurate estimates of cosmological parameters from these power spectra.