Source author record

James E. Geach

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

30works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

30 published item(s)

preprint2022arXiv

Host Dark Matter Halos of SDSS Red and Blue Quasars: No Significant Difference in Large-scale Environment

The observed optical colors of quasars are generally interpreted in one of two frameworks: unified models which attribute color to random orientation of the accretion disk along the line-of-sight, and evolutionary models which invoke connections between quasar systems and their environments. We test these schema by probing the dark matter halo environments of optically-selected quasars as a function of $g-i$ optical color by measuring the two-point correlation functions of $\sim$ 0.34 million eBOSS quasars as well as the gravitational deflection of cosmic microwave background photons around $\sim$ 0.66 million XDQSO photometric quasar candidates. We do not detect a trend of halo bias with optical color through either analysis, finding that optically-selected quasars at $0.8 < z < 2.2$ occupy halos of characteristic mass $M_{h}\sim 3\times 10^{12} \ h^{-1} M_{\odot}$ regardless of their color. This result implies that a quasar's large-scale halo environment is not strongly connected to its observed optical color. We also confirm findings of fundamental differences in the radio properties of red and blue quasars by stacking 1.4 GHz FIRST images at their positions, suggesting the observed differences cannot be attributed to orientation. Instead, the differences between red and blue quasars likely arise on nuclear-galactic scales, perhaps owing to reddening by a nuclear dusty wind. Finally, we show that optically-selected quasars' halo environments are also independent of their $r-W2$ optical-infrared colors, while previous work has suggested that mid-infrared-selected obscured quasars occupy more massive halos. We discuss implications of this result for models of quasar and galaxy co-evolution.

preprint2022arXiv

Realistic galaxy image simulation via score-based generative models

We show that a Denoising Diffusion Probabalistic Model (DDPM), a class of score-based generative model, can be used to produce realistic mock images that mimic observations of galaxies. Our method is tested with Dark Energy Spectroscopic Instrument (DESI) grz imaging of galaxies from the Photometry and Rotation curve OBservations from Extragalactic Surveys (PROBES) sample and galaxies selected from the Sloan Digital Sky Survey. Subjectively, the generated galaxies are highly realistic when compared with samples from the real dataset. We quantify the similarity by borrowing from the deep generative learning literature, using the `Fréchet Inception Distance' to test for subjective and morphological similarity. We also introduce the `Synthetic Galaxy Distance' metric to compare the emergent physical properties (such as total magnitude, colour and half light radius) of a ground truth parent and synthesised child dataset. We argue that the DDPM approach produces sharper and more realistic images than other generative methods such as Adversarial Networks (with the downside of more costly inference), and could be used to produce large samples of synthetic observations tailored to a specific imaging survey. We demonstrate two potential uses of the DDPM: (1) accurate in-painting of occluded data, such as satellite trails, and (2) domain transfer, where new input images can be processed to mimic the properties of the DDPM training set. Here we `DESI-fy' cartoon images as a proof of concept for domain transfer. Finally, we suggest potential applications for score-based approaches that could motivate further research on this topic within the astronomical community.

preprint2021arXiv

Physical Properties of Massive Compact Starburst Galaxies with Extreme Outflows

We present results on the nature of extreme ejective feedback episodes and the physical conditions of a population of massive ($\rm M_* \sim 10^{11} M_{\odot}$), compact starburst galaxies at z = 0.4-0.7. We use data from Keck/NIRSPEC, SDSS, Gemini/GMOS, MMT, and Magellan/MagE to measure rest-frame optical and near-IR spectra of 14 starburst galaxies with extremely high star formation rate surface densities (mean $\rm Σ_{SFR} \sim 3000 \,M_{\odot} yr^{-1} kpc^{-2}$) and powerful galactic outflows (maximum speeds v$_{98} \sim$ 1000-3000 km s$^{-1}$). Our unique data set includes an ensemble of both emission [OII]$λλ$3726,3729, H$β$, [OIII]$λλ$4959,5007, H$α$, [NII]$λλ$6548,6583, and [SII]$λλ$6716,6731) and absorption MgII$λλ$2796,2803, and FeII$λ$2586) lines that allow us to investigate the kinematics of the cool gas phase (T$\sim$10$^4$ K) in the outflows. Employing a suite of line ratio diagnostic diagrams, we find that the central starbursts are characterized by high electron densities (median n$_e \sim$ 530 cm$^{-3}$), and high metallicity (solar or super-solar). We show that the outflows are most likely driven by stellar feedback emerging from the extreme central starburst, rather than by an AGN. We also present multiple intriguing observational signatures suggesting that these galaxies may have substantial Lyman continuum (LyC) photon leakage, including weak [SII] nebular emission lines. Our results imply that these galaxies may be captured in a short-lived phase of extreme star formation and feedback where much of their gas is violently blown out by powerful outflows that open up channels for LyC photons to escape.

preprint2021arXiv

Rapid sorting of radio galaxy morphology using Haralick features

We demonstrate the use of Haralick features for the automated classification of radio galaxies. The set of thirteen Haralick features represent an extremely compact non-parametric representation of image texture, and are calculated directly from imagery using the Grey Level Co-occurrence Matrix (GLCM). The GLCM is an encoding of the relationship between the intensity of neighbouring pixels in an image. Using 10,000 sources detected in the first data release of the LOFAR Two-metre Sky Survey (LoTSS), we demonstrate that Haralick features are highly efficient, rotationally invariant descriptors of radio galaxy morphology. After calculating Haralick features for LoTSS sources, we employ the fast density-based hierarchical clustering algorithm HDBSCAN to group radio sources into a sequence of morphological classes, illustrating a simple methodology to classify and label new, unseen galaxies in large samples. By adopting a 'soft' clustering approach, we can assign each galaxy a probability of belonging to a given cluster, allowing for more flexibility in the selection of galaxies according to combinations of morphological characteristics and for easily identifying outliers: those objects with a low probability of belonging to any cluster in the Haralick space. Although our demonstration focuses on radio galaxies, Haralick features can be calculated for any image, making this approach also relevant to large optical imaging galaxy surveys.

preprint2021arXiv

Reproducing sub-millimetre galaxy number counts with cosmological hydrodynamic simulations

Matching the number counts of high-$z$ sub-millimetre-selected galaxies (SMGs) has been a long standing problem for galaxy formation models. In this paper, we use 3D dust radiative transfer to model the sub-mm emission from galaxies in the SIMBA cosmological hydrodynamic simulations, and compare predictions to the latest single-dish observational constraints on the abundance of 850$\mathrm{μm}$-selected sources. We find good agreement with the shape of the integrated 850$\mathrm{μm}$ luminosity function, and the normalisation is within 0.25 dex at $> 3 \; \mathrm{mJy}$, unprecedented for a fully cosmological hydrodynamic simulation, along with good agreement in the redshift distribution of bright SMGs. The agreement is driven primarily by SIMBA's good match to infrared measures of the star formation rate (SFR) function between $z = 2-4$ at high SFRs. Also important is the self-consistent on-the-fly dust model in SIMBA, which predicts, on average, higher dust masses (by up to a factor of 2.5) compared to using a fixed dust-to-metals ratio of 0.3. We construct a lightcone to investigate the effect of far-field blending, and find that 52% of sources are blends of multiple components, which makes a small contribution to the normalisation of the bright-end of the number counts. We provide new fits to the 850$\mathrm{μm}$ luminosity as a function of SFR and dust mass. Our results demonstrate that exotic solutions to the discrepancy between sub-mm counts in simulations and observations, such as a top-heavy IMF, are unnecessary, and that sub-millimetre-bright phases are a natural consequence of massive galaxy evolution.

preprint2020arXiv

Eigengalaxies: describing galaxy morphology using principal components in image space

We demonstrate how galaxy morphologies can be represented by weighted sums of "eigengalaxies" and how eigengalaxies can be used in a probabilistic framework to enable principled and simplified approaches in a variety of applications. Eigengalaxies can be derived from a Principal Component Analysis (PCA) of sets of single- or multi-band images. They encode the image space equivalent of basis vectors that can be combined to describe the structural properties of large samples of galaxies in a massively reduced manner. As an illustration, we show how a sample of 10,243 galaxies in the Hubble Space Telescope CANDELS survey can be represented by just 12 eigengalaxies. We show in some detail how this image space may be derived and tested. We also describe a probabilistic extension to PCA (PPCA) which enables the eigengalaxy framework to assign probabilities to galaxies. We present four practical applications of the probabilistic eigengalaxy framework that are particularly relevant for the next generation of large imaging surveys: we (i) show how low likelihood galaxies make for natural candidates for outlier detection (ii) demonstrate how missing data can be predicted (iii) show how a similarity search can be performed on exemplars (iv) demonstrate how unsupervised clustering of objects can be implemented.

preprint2020arXiv

Evidence for Infalling Gas in a Lyman-$α$ Blob

Lyman-$α$ blobs (LABs) are spatially extended nebulae of emission in the Ly$α$ line of hydrogen, seen at high redshifts$^{1,2}$, and most commonly found in the dense environment of star-forming galaxies$^{3,4}$. The origin of Ly$α$ emission in the LABs is still unclear and under debate$^{5}$. Proposed powering sources generally fall into two categories: (1) photoionization, galactic super-winds/outflows, resonant scattering of Ly$α$ photons from starbursts or active galactic nuclei (AGNs)$^{6,7,8,9,10}$ and (2) cooling radiation from cold streams of gas accreting onto galaxies$^{12}$. Here we analyze the gas kinematics within a LAB providing rare observational evidence for infalling gas. This is consistent with the release of gravitational accretion energy as cold streams radiate Ly$α$ photons. It also provides direct evidence for possible cold streams feeding the central galaxies. The infalling gas is not important by mass but hints at more than one mechanism to explain the origin of the extended Ly$α$ emission around young galaxies. It is also possible that the infalling gas may represent material falling back to the galaxy from where it originated, forming a galactic fountain.

preprint2019arXiv

Galaxy morphological classification in deep-wide surveys via unsupervised machine learning

Galaxy morphology is a fundamental quantity, that is essential not only for the full spectrum of galaxy-evolution studies, but also for a plethora of science in observational cosmology. While a rich literature exists on morphological-classification techniques, the unprecedented data volumes, coupled, in some cases, with the short cadences of forthcoming 'Big-Data' surveys (e.g. from the LSST), present novel challenges for this field. Large data volumes make such datasets intractable for visual inspection (even via massively-distributed platforms like Galaxy Zoo), while short cadences make it difficult to employ techniques like supervised machine-learning, since it may be impractical to repeatedly produce training sets on short timescales. Unsupervised machine learning, which does not require training sets, is ideally suited to the morphological analysis of new and forthcoming surveys. Here, we employ an algorithm that performs clustering of graph representations, in order to group image patches with similar visual properties and objects constructed from those patches, like galaxies. We implement the algorithm on the Hyper-Suprime-Cam Subaru-Strategic-Program Ultra-Deep survey, to autonomously reduce the galaxy population to a small number (160) of 'morphological clusters', populated by galaxies with similar morphologies, which are then benchmarked using visual inspection. The morphological classifications (which we release publicly) exhibit a high level of purity, and reproduce known trends in key galaxy properties as a function of morphological type at z<1 (e.g. stellar-mass functions, rest-frame colours and the position of galaxies on the star-formation main sequence). Our study demonstrates the power of unsupervised machine learning in performing accurate morphological analysis, which will become indispensable in this new era of deep-wide surveys.

preprint2016arXiv

The Impact of the Dusty Torus on Obscured Quasar Halo Mass Measurements

Recent studies have found that obscured quasars cluster more strongly and are thus hosted by dark matter haloes of larger mass than their unobscured counterparts. These results pose a challenge for the simplest unification models, in which obscured objects are intrinsically the same as unobscured sources but seen through a dusty line of sight. There is general consensus that a structure like a "dusty torus" exists, meaning that this intrinsic similarity is likely the case for at least some subset of obscured quasars. However, the larger host halo masses of obscured quasars implies that there is a second obscured population that has an even higher clustering amplitude and typical halo mass. Here, we use simple assumptions about the host halo mass distributions of quasars, along with analytical methods and cosmological $N$-body simulations to isolate the signal from this population. We provide values for the bias and halo mass as a function of the fraction of the "non-torus obscured" population. Adopting a reasonable value for this fraction of $\sim$25% implies a non-torus obscured quasar bias that is much higher than the observed obscured quasar bias, because a large fraction of the obscured population shares the same clustering strength as the unobscured objects. For this non-torus obscured population, we derive a bias of $\sim$3, and typical halo masses of $\sim3\times10^{13}$ M$_{\odot}/h$ at $z=1$. These massive haloes are likely the descendants of high-mass unobscured quasars at high redshift, and will evolve into members of galaxy groups at $z=0$.

preprint2016arXiv

The SCUBA-2 Cosmology Legacy Survey: Multi-wavelengths counterparts to 10$^3$ submillimeter galaxies in the UKIDSS-UDS field

We present multiwavelength identifications for the counterparts of 1088 submillimeter sources detected at 850$μ$m in the SCUBA-2 Cosmology Legacy Survey study of the UKIDSS-UDS field. By utilising an ALMA pilot study on a subset of our bright SCUBA-2 sample as a training set, along with the deep optical-near-infrared data available in this field, we develop a novel technique, Optical-IR Triple Color (OIRTC), using $z-K$, $K-[3.6]$, $[3.6]-[4.5]$ colors to select the candidate submillimeter galaxy (SMG) counterparts. By combining radio identification and the OIRTC technique, we find counterpart candidates for 80% of the Class = 1 $\geq4\,σ$ SCUBA-2 sample, defined as those that are covered by both radio and OIR imaging and the base sample for our scientific analyses. Based on the ALMA training set, we expect the accuracy of these identifications to be $82\pm20$%, with a completeness of $69\pm16$%, essentially as accurate as the traditional $p$-value technique but with higher completeness. We find that the fraction of SCUBA-2 sources having candidate counterparts is lower for fainter 850$μ$m sources, and we argue that for follow-up observations sensitive to SMGs with $S_{850}\gtrsim 1$ mJy across the whole ALMA beam, the fraction with multiple counterparts is likely to be $>40$% for SCUBA-2 sources at $S_{850} \gtrsim 4$ mJy. We find that the photometric redshift distribution for the SMGs is well fit by a lognormal distribution, with a median redshift of $z=2.3\pm0.1$. After accounting for the sources without any radio and/or OIRTC counterpart, we estimate the median redshift to be $z=2.6\pm0.1$ for SMGs with $S_{850} >1$ mJy. We also use this new large sample to study the clustering of SMGs and the the far-infrared properties of the unidentified submillimeter sources by stacking their Herschel SPIRE far-infrared emission.

preprint2015arXiv

A SCUBA-2 survey of FeLoBAL QSOs: Are FeLoBALs in a `transition phase' between ULIRGs and QSOs?

It is thought that a class of broad absorption line (BAL) QSOs, characterised by Fe absorption features in their UV spectra (called `FeLoBALs'), could mark a transition stage between the end of an obscured starburst event and a youthful QSO beginning to shed its dust cocoon, where Fe has been injected into the interstellar medium by the starburst. To test this hypothesis we have undertaken deep SCUBA-2 850 $μ$m observations of a sample of 17 FeLoBAL QSOs with 0.89 $\leq$ z $\leq$ 2.78 and -23.31 $\leq$ M$_{B}$ $\leq$-28.50 to directly detect an excess in the thermal emission of the dust which would probe enhanced star-formation activity. We find that FeLoBALs are not luminous sources in the submillimetre, none of them are individually detected at 850 $μ$m, nor as a population through stacking ($F_{s}=1.14\pm0.58$ mJy). Statistical and survival analyses reveal that FeLoBALs have sub-mm properties consistent with BAL and non-BAL QSOs with matched redshifts and magnitudes. An SED fitting analysis shows that the FIR emission is dominated by AGN activity, and a starburst component is required only in 6/17 sources of our sample; moreover the integrated total luminosity of 16/17 sources is L$\geq$10$^{12}$L$_{\odot}$, high enough to classify FeLoBALs as infrared luminous. In conclusion, we do not find any evidence in support of FeLoBAL QSOs being a transition population between a ULIRG and an unobscured QSO; in particular, FeLoBALs are not characterized by a cold starburst which would support this hypothesis.

preprint2015arXiv

Star formation in quasar hosts and the origin of radio emission in radio-quiet quasars

Radio emission from radio-quiet quasars may be due to star formation in the quasar host galaxy, to a jet launched by the supermassive black hole, or to relativistic particles accelerated in a wide-angle radiatively-driven outflow. In this paper we examine whether radio emission from radio-quiet quasars is a byproduct of star formation in their hosts. To this end we use infrared spectroscopy and photometry from Spitzer and Herschel to estimate or place upper limits on star formation rates in hosts of ~300 obscured and unobscured quasars at z<1. We find that low-ionization forbidden emission lines such as [NeII] and [NeIII] are likely dominated by quasar ionization and do not provide reliable star formation diagnostics in quasar hosts, while PAH emission features may be suppressed due to the destruction of PAH molecules by the quasar radiation field. While the bolometric luminosities of our sources are dominated by the quasars, the 160 micron fluxes are likely dominated by star formation, but they too should be used with caution. We estimate median star formation rates to be 6-29 Msun/year, with obscured quasars at the high end of this range. This star formation rate is insufficient to explain the observed radio emission from quasars by an order of magnitude, with log(L_radio, observed/L_radio, SF)=0.6-1.3 depending on quasar type and star formation estimator. Although radio-quiet quasars in our sample lie close to the 8-1000 micron infrared / radio correlation characteristic of the star-forming galaxies, both their infrared emission and their radio emission are dominated by the quasar activity, not by the host galaxy.

preprint2015arXiv

Teaching a machine to see: unsupervised image segmentation and categorisation using growing neural gas and hierarchical clustering

We present a novel unsupervised learning approach to automatically segment and label images in astronomical surveys. Automation of this procedure will be essential as next-generation surveys enter the petabyte scale: data volumes will exceed the capability of even large crowd-sourced analyses. We demonstrate how a growing neural gas (GNG) can be used to encode the feature space of imaging data. When coupled with a technique called hierarchical clustering, imaging data can be automatically segmented and labelled by organising nodes in the GNG. The key distinction of unsupervised learning is that these labels need not be known prior to training, rather they are determined by the algorithm itself. Importantly, after training a network can be be presented with images it has never 'seen' before and provide consistent categorisation of features. As a proof-of-concept we demonstrate application on data from the Hubble Space Telescope Frontier Fields: images of clusters of galaxies containing a mixture of galaxy types that would easily be recognised and classified by a human inspector. By training the algorithm using one field (Abell 2744) and applying the result to another (MACS0416.1-2403), we show how the algorithm can cleanly separate image features that a human would associate with early and late type galaxies. We suggest that the algorithm has potential as a tool in the automatic analysis and data mining of next-generation imaging and spectral surveys, and could also find application beyond astronomy.

preprint2014arXiv

A relationship between specific star formation rate and metallicity gradient within z=1 galaxies from KMOS-HiZELS

We have observed a sample of typical z=1 star forming galaxies, selected from the HiZELS survey, with the new KMOS near-infrared, multi-IFU instrument on the VLT, in order to obtain their dynamics and metallicity gradients. The majority of our galaxies have a metallicity gradient consistent with being flat or negative (i.e. higher metallicity cores than outskirts). Intriguingly, we find a trend between metallicity gradient and specific star formation rate (sSFR), such that galaxies with a high sSFR tend to have relatively metal-poor centres, a result which is strengthened when combined with datasets from the literature. This result appears to explain the discrepancies reported between different high redshift studies and varying claims for evolution. From a galaxy evolution perspective, the trend we see would mean that a galaxy's sSFR is governed by the amount of metal poor gas that can be funnelled into its core, triggered either by merging or through efficient accretion. In fact merging may play a significant role as it is the starburst galaxies at all epochs, which have the more positive metallicity gradients. Our results may help to explain the origin of the fundamental metallicity relation, in which galaxies at a fixed mass are observed to have lower metallicities at higher star formation rates, especially if the metallicity is measured in an aperture encompassing only the central regions of the galaxy. Finally, we note that this study demonstrates the power of KMOS as an efficient instrument for large scale resolved galaxy surveys.

preprint2014arXiv

PdBI Cold Dust Imaging of Two Extremely Red H-[4.5]>4 Galaxies Discovered with SEDS and CANDELS

We report Plateau de Bure Interferometer (PdBI) 1.1 mm continuum imaging towards two extremely red H-[4.5]>4 (AB) galaxies at z>3, which we have previously discovered making use of Spitzer SEDS and Hubble Space Telescope CANDELS ultra-deep images of the Ultra Deep Survey field. One of our objects is detected on the PdBI map with a 4.3 sigma significance, corresponding to Snu(1.1mm)=(0.78 +/- 0.18) mJy. By combining this detection with the Spitzer 8 and 24 micron photometry for this source, and SCUBA2 flux density upper limits, we infer that this galaxy is a composite active galactic nucleus/star-forming system. The infrared (IR)-derived star formation rate is SFR~(200 +/- 100) Msun/yr, which implies that this galaxy is a higher-redshift analogue of the ordinary ultra-luminous infrared galaxies (ULIRGs) more commonly found at z~2-3. In the field of the other target, we find a tentative 3.1 sigma detection on the PdBI 1.1 mm map, but 3.7 arcsec away of our target position, so it likely corresponds to a different object. In spite of the lower significance, the PdBI detection is supported by a close SCUBA2 3.3 sigma detection. No counterpart is found on either the deep SEDS or CANDELS maps, so, if real, the PdBI source could be similar in nature to the sub-millimetre source GN10. We conclude that the analysis of ultra-deep near- and mid-IR images offers an efficient, alternative route to discover new sites of powerful star formation activity at high redshifts.

preprint2013arXiv

A fundamental metallicity relation for galaxies at z = 0.84 - 1.47 from HiZELS

We obtained Subaru FMOS observations of Halpha emitting galaxies selected from the HiZELS narrow-band survey, to investigate the relationship between stellar mass, metallicity and star-formation rate at z = 0.84 - 1.47, for comparison with the Fundamental Metallicity Relation seen at low redshift. Our findings demonstrate, for the first time with a homogeneously selected sample, that a relationship exists for typical star-forming galaxies at z = 1 - 1.5 and that it is surprisingly similar to that seen locally. Therefore, star-forming galaxies at z = 1 - 1.5 are no less metal abundant than galaxies of similar mass and star formation rate (SFR) at z = 0.1, contrary to claims from some earlier studies. We conclude that the bulk of the metal enrichment for this star-forming galaxy population takes place in the 4 Gyr before z = 1.5. We fit a new mass-metallicity-SFR plane to our data which is consistent with other high redshift studies. However, there is some evidence that the mass-metallicity component of this high redshift plane is flattened, at all SFR, compared with z = 0.1, suggesting that processes such as star-formation driven winds, thought to remove enriched gas from low mass halos, are yet to have as large an impact at this early epoch. The negative slope of the SFR-metallicity relation from this new plane is consistent with the picture that the elevation in the SFR of typical galaxies at z > 1 is fuelled by the inflow of metal-poor gas and not major merging.

preprint2013arXiv

On the evolution and environmental dependence of the star formation rate versus stellar mass relation since z~2

This paper discusses the evolution of the correlation between galaxy star formation rates (SFRs) and stellar mass (M*) over the last ~10 Gyrs, particularly focusing on its environmental dependence. We first present the mid-infrared (MIR) properties of the H-alpha-selected galaxies in a rich cluster Cl0939+4713 at z=0.4. With a MIR stacking analysis, we find that the median SFR of H-alpha emitters is higher in higher-density environment at z=0.4. We also find that star-forming galaxies in high-density environment tend to have higher specific SFR (SSFR), although the trend is much less significant compared to that of SFR. This increase of SSFR in high-density environment is not visible when we consider the SFR derived from H-alpha alone, suggesting that the dust attenuation in galaxies depends on environment; galaxies in high-density environment tend to be dustier (by up to ~0.5 mag). We then discuss the environmental dependence of the SFR-M* relation for star-forming galaxies since z~2, by compiling our comparable, narrow-band selected, large H-alpha emitter samples in both distant cluster environments (from MAHALO-Subaru) and field environments (from HiZELS). We find that the SSFR of H-alpha-selected galaxies (at the fixed mass of log(M*/Msun)=10) rapidly evolves as (1+z)^3, but the SFR-M* relation is independent of the environment since z~2, as far as we rely on the H-alpha-based SFRs (with M*-dependent extinction correction). Even if we consider the possible environmental variation in the dust attenuation, we conclude that the difference in the SFR-M* relation between cluster and field star-forming galaxies is always small (~<0.2dex level) at any time in the history of the Universe since z~2. [abridged]

preprint2013arXiv

The Evolution of Dusty Star formation in Galaxy Clusters to z = 1: Spitzer IR Observations of the First Red-Sequence Cluster Survey

We present an IR study of high-redshift galaxy clusters with the MIPS camera on the Spitzer Space Telescope. Employing a sample of 42 clusters from the RCS-1 over the redshift range 0.3 < z < 1.0 and spanning an approximate range in mass of 10^{14-15} Msun, we show the number of IR-luminous galaxies in clusters above a fixed IR luminosity of 2x10^{11} Msun per unit cluster mass evolves as (1+z)^{5.1+/-1.9}. These results assume a single star forming galaxy template; the presence of AGN, and an evolution in their relative contribution to the mid-IR galaxy emission, will alter the overall number counts per cluster and their rate of evolution. We infer the total SFR per unit cluster mass and find T_SFR/M_c ~ (1+z)^{5.4+/-1.9}. This evolution can be attributed entirely to the change in the in-falling field galaxy population. The T_SFR/M_c (binned over all redshift) decreases with increasing cluster mass with a slope (T_SFR/M_c ~ M_c^{-1.5+/-0.4}) consistent with the dependence of the stellar-to-total mass per unit cluster mass seen locally. The inferred star formation seen here could produce ~5-10% of the total stellar mass in massive clusters at z = 0. Finally, we show a clear decrease in the number of IR-bright galaxies per unit optical galaxy in the cluster cores, confirming star formation continues to avoid the highest density regions of the universe at z ~ 0.75 (the average redshift of the high-redshift clusters). While several previous studies appear to show enhanced star formation in high-redshift clusters relative to the field we note that these papers have not accounted for the overall increase in galaxy or dark matter density at the location of clusters. Once this is done, clusters at z ~ 0.75 have the same or less star formation per unit mass or galaxy as the field.

preprint2013arXiv

The merger rates and sizes of galaxies across the peak epoch of star formation from the HiZELS survey

We use the HiZELS narrow-band H-alpha survey in combination with CANDELS, UKIDSS and WIRDS near-infrared imaging, to investigate the morphologies, merger rates and sizes of a sample of H-alpha emitting galaxies in the redshift range z=0.40 - 2.23, an epoch encompassing the rise to the peak of the star formation rate density. Merger rates are estimated from space- and ground-based imaging using the M20 coefficient. To account for the increase in the specific star-formation rate (sSFR) of the star forming `main-sequence' with redshift, we normalise the star-formation rate of galaxies at each epoch to the typical value derived from the H-alpha luminosity function. Once this trend in sSFR is removed we see no evidence for an increase in the number density of star-forming galaxies or the merger rate with redshift. We thus conclude that neither is the main driver of the enhanced star-formation rate density at z=1-2, with secular processes such as instabilities within efficiently fuelled, gas-rich discs or multiple minor mergers the most likely alternatives. However, we find that 40-50% of starburst galaxies, those with enhanced specific star formation at their epoch, are major mergers and this fraction is redshift independent. Finally, we find the surprising result that the typical size of a star-forming galaxy of a given mass does not evolve across the redshift range considered, suggesting a universal size-mass relation. Taken in combination, these results indicate a star-forming galaxy population that is statistically similar in physical size, merger rate and mass over the 6 Gyr covered in this study, despite the increase in typical sSFR.

preprint2012arXiv

A large Hα survey at z=2.23, 1.47, 0.84 & 0.40: the 11 Gyr evolution of star-forming galaxies from HiZELS

This paper presents new deep and wide narrow-band surveys undertaken with UKIRT, Subaru and the VLT; a combined effort to select large, robust samples of H-alpha (Ha) emitters at z=0.40, 0.84, 1.47 and 2.23 (corresponding to look-back times of 4.2, 7.0, 9.2 and 10.6 Gyrs) in a uniform manner over ~2 deg^2 in the COSMOS and UDS fields. The deep Ha surveys reach ~3M_sun/yr out to z=2.2 for the first time, while the wide area and the coverage over two independent fields allow to greatly overcome cosmic variance. A total of 1742, 637, 515 and 807 Ha emitters are homogeneously selected at z=0.40, 0.84, 1.47 and 2.23, respectively, and used to determine the Ha luminosity function and its evolution. The faint-end slope is found to be -1.60+-0.08 over z=0-2.23, showing no evolution. The characteristic luminosity of SF galaxies, L*, evolves significantly as log[L*(z)]=0.45z+log[L*(z=0)]. This is the first time Ha has been used to trace SF activity with a single homogeneous survey at z=0.4-2.23. Overall, the evolution seen in Ha is in good agreement with the evolution seen using inhomogeneous compilations of other tracers of star formation, such as FIR and UV, jointly pointing towards the bulk of the evolution in the last 11 Gyrs being driven by a strong luminosity increase from z~0 to z~2.2. Our uniform analysis reveals an Ha star formation history of the Universe which can be simply parameterised by log(SFRD)=-2.1/(1+z) for z<2.23. Both the shape and normalisation of the Ha star formation history are consistent with the measurements of the stellar mass density growth, confirming that our Ha analysis traces the bulk of the formation of stars in the Universe up to z~2.2. Star formation activity since z~2.2 is responsible for ~95% of the total stellar mass density observed locally today (with half being assembled between z=1.2-2.2, and the rest since z<1.2).

preprint2012arXiv

Effects of cosmological model assumptions on galaxy redshift survey measurements

The clustering of galaxies observed in future redshift surveys will provide a wealth of cosmological information. Matching the signal at different redshifts constrains the dark energy driving the acceleration of the expansion of the Universe. In tandem with these geometrical constraints, redshift-space distortions (RSD) depend on the build up of large-scale structure. As pointed out by many authors measurements of these effects are intrinsically coupled. We investigate this link, and argue that it strongly depends on the cosmological assumptions adopted when analysing data. Using representative assumptions for the parameters of the Euclid survey in order to provide a baseline future experiment, we show how the derived constraints change due to different model assumptions. We argue that even the assumption of a Friedman-Robertson-Walker (FRW) space-time is sufficient to reduce the importance of the coupling to a significant degree. Taking this idea further, we consider how the data would actually be analysed and argue that we should not expect to be able to simultaneously constrain multiple deviations from the standard $Λ$CDM model. We therefore consider different possible ways in which the Universe could deviate from the $Λ$CDM model, and show how the coupling between geometrical constraints and structure growth affects the measurement of such deviations.

preprint2012arXiv

Herschel-PACS observations of [OI]63um towards submillimetre galaxies at z~1

We present Herschel-PACS spectroscopy of the [OI]63um far-infrared cooling line from a sample of six unlensed and spectroscopically-confirmed 870um-selected submillimetre (submm) galaxies (SMGs) at 1.1<z<1.6 from the LABOCA Extended Chandra Deep Field South (ECDFS) Submm Survey (LESS). This is the first survey of [OI]63um, one of the main photodissociation region (PDR) cooling lines, in SMGs. New high-resolution ALMA interferometric 870um continuum imaging confirms that these six Herschel-targeted SMG counterparts are bona fide sources of submm emission. We detect [OI]63um in two SMGs with a SNR >3, tentatively detect [OI]63um in one SMG, and constrain the line flux for the non-detections. We also exploit the combination of submm continuum photometry from 250-870um and our new PACS continuum measurements to constrain the far-infrared (FIR) luminosity, L_FIR, in these SMGs to < 30%. We find that SMGs do not show a deficit in their [OI]63um-to-far-infrared continuum luminosity ratios (with ratios ranging from ~0.5-1.5%), similar to what was seen previously for the [CII]158um-to-FIR ratios in SMGs. These observed ratios are about an order of magnitude higher than what is seen typically for local ultra luminous infrared galaxies (ULIRGs), which adds to the growing body of evidence that SMGs are not simply `scaled up' versions of local ULIRGs. Rather, the PDR line-to-L_FIR ratios suggest that the star formation modes of SMGs are likely more akin to that of local normal (lower-luminosity) star-forming galaxies, with the bulk of the star formation occurring in extended regions, galaxy-scale (~kpc) in size. These observations represent the first step towards a census of the major PDR cooling lines in typical SMGs that will be attainable with ALMA, enabling detailed modelling to probe the global properties of the star formation and the evolutionary status of SMGs.

preprint2012arXiv

Molecular and atomic line surveys of galaxies I: the dense, star-forming phase as a beacon

We predict the space density of molecular gas reservoirs in the Universe, and place a lower limit on the number counts of carbon monoxide (CO), hydrogen cyanide (HCN) molecular and [CII] atomic emission lines in blind redshift surveys in the submillimeter-centimeter spectral regime. Our model uses: (a) recently available HCN Spectral Line Energy Distributions (SLEDs) of local Luminous Infrared Galaxies (LIRGs, L_IR>10^11 L_sun), (b) a value for epsilon=SFR/M_dense(H_2) provided by new developments in the study of star formation feedback on the interstellar medium and (c) a model for the evolution of the infrared luminosity density. Minimal 'emergent' CO SLEDs from the dense gas reservoirs expected in all star-forming systems in the Universe are then computed from the HCN SLEDs since warm, HCN-bright gas will necessarily be CO-bright, with the dense star-forming gas phase setting an obvious minimum to the total molecular gas mass of any star-forming galaxy. We include [CII] as the most important of the far-infrared cooling lines. Optimal blind surveys with the Atacama Large Millimeter Array (ALMA) could potentially detect very distant (z~10-12) [CII] emitters in the >ULIRG galaxy class at a rate of ~0.1-1 per hour (although this prediction is strongly dependent on the star formation and enrichment history at this early epoch), whereas the (high-frequency) Square Kilometer Array (SKA) will be capable of blindly detecting z>3 low-J CO emitters at a rate of ~40-70 per hour. The [CII] line holds special promise for the detection of metal-poor systems with extensive reservoirs of CO-dark molecular gas where detection rates with ALMA can reach up to 2-7 per hour in Bands 4-6.

preprint2012arXiv

Molecular and atomic line surveys of galaxies II: unbiased estimates of their star formation mode

We make use of our 'minimal' cold interstellar medium (ISM) emission line model that predicts the molecular and atomic line emission per unit dense, star-forming gas mass (Geach & Papadopoulos 2012; Paper I) to examine the utility of key line ratios in surveys of the so-called star formation 'mode' as traced by xi_SF = M_dense(H_2)/M_total(H_2). We argue that xi_SF and its proxies provide very sensitive, extinction-free discriminators of rapid starburst/merger-driven versus secular quiescent/disk-like stellar mass assembly, with the most promising diagnostic to be applied in the near-future being CO(4-3)/[CI](1-0). These lines are accessible across nearly the full range 0<z<2 (thus covering the bulk of galaxy evolution) with the Atacama Large Millimeter Array. In addition to their diagnostic power, another advantage of this combination is the similar observed frequencies (Delta nu_0 ~ 30 GHz) of the lines, resulting in nearly spatially-matched beams for a fixed aperture, thus mitigating the effects of resolution/morphology bias in the interpretation of galaxy-averaged line ratios. Finally we discuss the capability of deep blind redshift surveys with the high frequency component of the Square Kilometer Arrray (SKA) in discovering H_2-rich galaxies with very low xi_SF values. These could be the progenitors of starburst galaxies seen prior to the onset of star formation; such galaxies could be a class of extreme outliers from local (gas surface density)-(star formation rate) scaling laws, which would exclude them from current star formatation or stellar mass selected samples. Our conservative model suggests that SKA could detect such systems residing at z~3 at a rate of 20-200 per hour.

preprint2011arXiv

4098 galaxy clusters to z~0.6 in the Sloan Digital Sky Survey equatorial Stripe 82

We present a catalogue of 4098 photometrically selected galaxy clusters with a median redshift <z> = 0.32 in the 270 square degree 'Stripe 82' region of the Sloan Digital Sky Survey (SDSS), covering the celestial equator in the Southern Galactic Cap (-50 < RA < 59 deg, |Dec| < 1.25 deg). Owing to the multi-epoch SDSS coverage of this region, the ugriz photometry is ~2 magnitudes deeper than single scans within the main SDSS footprint. We exploit this to detect clusters of galaxies using an algorithm that searches for statistically significant overdensities of galaxies in a Voronoi tessellation of the projected sky. 32% of the clusters have at least one member with a spectroscopic redshift from existing public data (SDSS Data Release 7, 2SLAQ & WiggleZ), and the remainder have a robust photometric redshift (accurate to ~5-9% at the median redshift of the sample). The weighted average of the member galaxies' redshifts provides a reasonably accurate estimate of the cluster redshift. The cluster catalogue is publicly available for exploitation by the community to pursue a range of science objectives. In addition to the cluster catalogue, we provide a linked catalogue of 18,295 V<21 mag quasar sight-lines with impact parameters within <3 Mpc of the cluster cores selected from the catalogue of Veron et al. (2010). The background quasars cover 0.25 < z < 2, where MgII absorption-line systems associated with the clusters are detectable in optical spectra.

preprint2011arXiv

On the evolution of the molecular gas fraction of star forming galaxies

We present IRAM Plateau de Bure interferometric detections of CO(1-0) emission from a 24um-selected sample of star-forming galaxies at z=0.4. The galaxies have PAH 7.7um-derived star formation rates of SFR~30-60 M_Sun/yr and stellar masses M*~10^{11} M_Sun. The CO(1-0) luminosities of the galaxies imply that the disks still contain a large reservoir of molecular gas, contributing ~20% of the baryonic mass, but have star-formation 'efficiencies' similar to local quiescent disks and gas-dominated disks at z~1.5-2. We reveal evidence that the average molecular gas fraction has undergone strong evolution since z~2, with f_gas ~ (1+z)^{2 +/- 0.5}. The evolution of f_gas encodes fundamental information about the relative depletion/replenishment of molecular fuel in galaxies, and is expected to be a strong function of halo mass. We show that the latest predictions for the evolution of the molecular gas fraction in semi-analytic models of galaxy formation within a LCDM Universe are supported by these new observations.

preprint2011arXiv

Unsupervised self-organised mapping: a versatile empirical tool for object selection, classification and redshift estimation in large surveys

We present an application of unsupervised machine learning - the self-organised map (SOM) - as a tool for visualising, exploring and mining the catalogues of large astronomical surveys. Self-organisation culminates in a low-resolution representation of the 'topology' of a parameter volume, and this can be exploited in various ways pertinent to astronomy. Using data from the Cosmological Evolution Survey (COSMOS), we demonstrate two key astronomical applications of the SOM: (i) object classification and selection, using the example of galaxies with active galactic nuclei as a demonstration, and (ii) photometric redshift estimation, illustrating how SOMs can be used as totally empirical predictive tools. With a training set of ~3800 galaxies with z_spec<1, we achieve photometric redshift accuracies competitive with other (mainly template fitting) techniques that use a similar number of photometric bands (sigma(Dz)=0.03 with a ~2% outlier rate when using u*-band to 8um photometry). We also test the SOM as a photo-z tool using the PHoto-z Accuracy Testing (PHAT) synthetic catalogue of Hildebrandt et al. (2010), which compares several different photo-z codes using a common input/training set. We find that the SOM can deliver accuracies that are competitive with many of the established template-fitting and empirical methods. This technique is not without clear limitations, which are discussed, but we suggest it could be a powerful tool in the era of extremely large - 'petabyte' - databases where efficient data-mining is a paramount concern.

preprint2010arXiv

Designing a space-based galaxy redshift survey to probe dark energy

A space-based galaxy redshift survey would have enormous power in constraining dark energy and testing general relativity, provided that its parameters are suitably optimized. We study viable space-based galaxy redshift surveys, exploring the dependence of the Dark Energy Task Force (DETF) figure-of-merit (FoM) on redshift accuracy, redshift range, survey area, target selection, and forecast method. Fitting formulae are provided for convenience. We also consider the dependence on the information used: the full galaxy power spectrum P(k, P(k) marginalized over its shape, or just the Baryon Acoustic Oscillations (BAO). We find that the inclusion of growth rate information (extracted using redshift space distortion and galaxy clustering amplitude measurements) leads to a factor of ~ 3 improvement in the FoM, assuming general relativity is not modified. This inclusion partially compensates for the loss of information when only the BAO are used to give geometrical constraints, rather than using the full P(k) as a standard ruler. We find that a space-based galaxy redshift survey covering ~20,000 deg^2 over 0.5< z < 2 with σ_z/(1+z)<= 0.001 exploits a redshift range that is only easily accessible from space, extends to sufficiently low redshifts to allow both a vast 3-D map of the universe using a single tracer population, and overlaps with ground-based surveys to enable robust modeling of systematic effects. We argue that these parameters are close to their optimal values given current instrumental and practical constraints.

preprint2010arXiv

The clustering and evolution of H-alpha emitters at z~1 from HiZELS

(Abridged) The clustering properties of a well-defined sample of 734 H-alpha emitters at z=0.84 obtained as part of the Hi-z Emission Line Survey (HiZELS) are investigated. The spatial correlation function is very well-described by (r/r_0)^-1.8, with r_0=2.7+-0.3Mpc/h. The correlation length r_0 increases strongly with H-alpha luminosity, L_H-alpha, from r_0~2Mpc/h for the most quiescent galaxies (star-formation rates of ~4M_sun/yr), up to r_0>5Mpc/h for the brightest galaxies in H-alpha. The correlation length also increases with increasing rest-frame K-band luminosity (M_K), but the r_0-L_H-alpha correlation maintains its full statistical significance at fixed M_K. At z=0.84, star-forming galaxies classified as irregulars or mergers are much more clustered than discs and non-mergers, but once the samples are matched in L_H-alpha and M_K, the differences vanish, implying that the clustering is independent of morphological type at z~1. The typical H-alpha emitters found at z=0.84 reside in dark-matter haloes of ~10^12M_sun, but those with the highest SFRs reside in more massive haloes of ~10^13M_sun. Comparing the results with those of H-alpha surveys at different redshifts, it is seen that although the break of the H-alpha luminosity function, L*, evolves by a factor of ~30 from z=0.24 to z=2.23, galaxies with the same L_H-alpha/L*(z) are found in dark matter haloes of similar masses, independently of cosmic time. This not only confirms that star-formation is more efficient at higher redshift, but also suggests a fundamental connection between the strong decrease of L* since z=2.23 and the quenching of star-formation in galaxies residing within dark-matter haloes significantly more massive than 10^12M_sun at any given epoch.

preprint2009arXiv

Obscured star formation at z = 0.84 with HiZELS: the relationship between star formation rate and H-alpha or ultra-violet dust extinction

[Abridged] We compare H-alpha, ultraviolet (UV) and infrared (IR) indicators of star formation rate (SFR) for a sample of z = 0.84 galaxies from the High-z Emission Line Survey (HiZELS). Using multi-wavelength diagnostics, we estimate that 5 - 11 per cent of H-alpha emitters at this redshift are active galactic nuclei. We detect 35 per cent of the H-alpha emitters individually at 24 microns, and stack the star-forming emitters on deep 24-micron images in order to calculate the typical SFRs of our galaxies. These are compared to the observed H-alpha line fluxes in order to estimate the extinction at z = 0.84, and we find a significant increase in dust extinction for galaxies with higher SFRs. We demonstrate that the local relationship between SFR and extinction is suitable for our sample, and attribute the overall increase in dust extinction for z = 0.84 galaxies to an increase in the average SFR, rather than to a change in dust properties. We calculate the UV extinction, and find a similar dependence on SFR, but no evidence for a 2175 angstrom UV bump in the dust attenuation law. We calculate the conversion between the dust attenuation of nebular and stellar radiation, gamma, and show that gamma = 0.50 +/- 0.14. The extinction / SFR relationship is shown to be applicable to galaxies with a range of morphologies and bulge-to-disk ratios, to both merging and non-merging galaxies, and to galaxies within high- and low-density environments, implying that it is a fundamental property of star-forming regions. In order to allow future studies to correct for a SFR-dependent amount of dust extinction, we present an equation to predict the extinction of a galaxy, based solely upon its observed H-alpha luminosity, and use this to recalculate the H-alpha luminosity function and star formation rate density at z = 0.84.