Source author record

Seb Oliver

Seb Oliver 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

24works
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

24 published item(s)

preprint2022arXiv

A hyperluminous obscured quasar at a redshift of z ~ 4.3

In this work we report the discovery of the hyperluminous galaxy HELP_J100156.75+022344.7 at the photometric redshift of z ~ 4.3. The galaxy was discovered in the Cosmological Evolution Survey (COSMOS) field, one of the fields studied by the Herschel Extragalactic Legacy Project (HELP). We present the spectral energy distribution (SED) of the galaxy and fit it with the CYprus models for Galaxies and their NUclear Spectra (CYGNUS) multi-component radiative transfer models. We find that its emission is dominated by an obscured quasar with a predicted total 1-1000um luminosity of $3.91^{+1.69}_{-0.55} \times 10^{13} L_\odot$ and an active galactic nucleus (AGN) fraction of ~89%. We also fit HELP_J100156.75+022344.7 with the Code Investigating GALaxy Emission (CIGALE) code and find a similar result. This is only the second z > 4 hyperluminous obscured quasar discovered to date. The discovery of HELP_J100156.75+022344.7 in the ~ 2deg^2 COSMOS field implies that a large number of obscured hyperluminous quasars may lie in the HELP fields which cover ~ 1300deg^2. If this is confirmed, tension between supermassive black hole evolution models and observations will be alleviated. We estimate the space density of objects like HELP_J100156.75+022344.7 at z ~ 4.5 to be $\sim 1.8 \times 10^{-8}$Mpc$^{-3}$. This is slightly higher than the space density of coeval hyperluminous optically selected quasars suggesting that the obscuring torus in z > 4 quasars may have a covering factor $\gtrsim 50\%$.

preprint2022arXiv

Stellar and black hole assembly in z<0.3 infrared-luminous mergers: intermittent starbursts vs. super-Eddington accretion

We study stellar and black hole mass assembly in a sample of 42 infrared-luminous galaxy mergers at z<0.3 by combining results from radiative transfer modelling with archival measures of molecular gas and black hole mass. The ratios of stellar mass, molecular gas mass, and black hole mass to each other are consistent with those of massive gas-rich galaxies at z<0.3. The advanced mergers may show increased black hole mass to stellar mass ratios, consistent with the transition from AGN to ellipticals and implying substantial black hole mass growth over the course of the merger. Star formation rates are enhanced relative to the local main sequence, by factors of ~100 in the starburst and ~1.8 in the host. The starburst star formation rates appear distinct to star formation in the main sequence at all redshifts up to at least z~5. Starbursts may prefer late-stage mergers, but are observed at any merger stage. We do not find evidence that the starbursts in these low-redshift systems substantially increase the total stellar mass, with a soft upper limit on the stellar mass increase from starburst activity of about a factor of two. In contrast, 12 objects show evidence for super-Eddington accretion, associated with late-stage mergers, suggesting that many AGN in infrared-luminous mergers go through a super-Eddington phase. The super-Eddington phase may increase black hole mass by up to an order of magnitude at an accretion efficiency of 42+/-33% over a period of 44+/-22Myr. Our results imply that super-Eddington accretion is an important black hole growth channel in infrared-luminous galaxies at all redshifts.

preprint2020arXiv

An empirical, Bayesian approach to modelling the impact of weather on crop yield: maize in the US

We apply an empirical, data-driven approach for describing crop yield as a function of monthly temperature and precipitation by employing generative probabilistic models with parameters determined through Bayesian inference. Our approach is applied to state-scale maize yield and meteorological data for the US Corn Belt from 1981 to 2014 as an exemplar, but would be readily transferable to other crops, locations and spatial scales. Experimentation with a number of models shows that maize growth rates can be characterised by a two-dimensional Gaussian function of temperature and precipitation with monthly contributions accumulated over the growing period. This approach accounts for non-linear growth responses to the individual meteorological variables, and allows for interactions between them. Our models correctly identify that temperature and precipitation have the largest impact on yield in the six months prior to the harvest, in agreement with the typical growing season for US maize (April to September). Maximal growth rates occur for monthly mean temperature 18-19$^\circ$C, corresponding to a daily maximum temperature of 24-25$^\circ$C (in broad agreement with previous work) and monthly total precipitation 115 mm. Our approach also provides a self-consistent way of investigating climate change impacts on current US maize varieties in the absence of adaptation measures. Keeping precipitation and growing area fixed, a temperature increase of $2^\circ$C, relative to 1981-2014, results in the mean yield decreasing by 8\%, while the yield variance increases by a factor of around 3. We thus provide a flexible, data-driven framework for exploring the impacts of natural climate variability and climate change on globally significant crops based on their observed behaviour. In concert with other approaches, this can help inform the development of adaptation strategies that will ensure food security under a changing climate.

preprint2020arXiv

Forecasting vegetation condition for drought early warning systems in pastoral communities in Kenya

Droughts are a recurring hazard in sub-Saharan Africa, that can wreak huge socioeconomic costs.Acting early based on alerts provided by early warning systems (EWS) can potentially provide substantial mitigation, reducing the financial and human cost. However, existing EWS tend only to monitor current, rather than forecast future, environmental and socioeconomic indicators of drought, and hence are not always sufficiently timely to be effective in practice. Here we present a novel method for forecasting satellite-based indicators of vegetation condition. Specifically, we focused on the 3-month Vegetation Condition Index (VCI3M) over pastoral livelihood zones in Kenya, which is the indicator used by the Kenyan National Drought Management Authority(NDMA). Using data from MODIS and Landsat, we apply linear autoregression and Gaussian process modeling methods and demonstrate high forecasting skill several weeks ahead. As a benchmark we predicted the drought alert marker used by NDMA (VCI3M<35). Both of our models were able to predict this alert marker four weeks ahead with a hit rate of around 89% and a false alarm rate of around 4%, or 81% and 6% respectively six weeks ahead. The methods developed here can thus identify a deteriorating vegetation condition well and sufficiently in advance to help disaster risk managers act early to support vulnerable communities and limit the impact of a drought hazard.

preprint2020arXiv

Rise of the Titans: Gas Excitation and Feedback in a Binary Hyper-Luminous Dusty Starburst Galaxy at z~6

We report new observations toward the hyper-luminous dusty starbursting major merger ADFS-27 (z=5.655), using ATCA and ALMA. We detect CO 2-1, 8-7, 9-8, 10-9 and H2O(321-221) emission, and a P-Cygni-shaped OH+(11-01) absorption/emission feature. We also tentatively detect H2O(321-312) and OH+(12-01) emission and CH+(1-0) absorption. We find a total cold molecular mass of M_gas = (2.1+/-0.2) x 10^11 (alpha_CO/1.0) Msun. We also find that the excitation of the star-forming gas is overall moderate for a z>5 dusty starburst, which is consistent with its moderate dust temperature. A high density, high kinetic temperature gas component embedded in the gas reservoir is required to fully explain the CO line ladder. This component is likely associated with the "maximum starburst" nuclei in the two merging galaxies, which are separated by only (140+/-13) km/s along the line of sight and 9.0 kpc in projection. The kinematic structure of both components is consistent with galaxy disks, but this interpretation remains limited by the spatial resolution of the current data. The OH+ features are only detected towards the northern component, which is also the one that is more enshrouded in dust and thus remains undetected up to 1.6 um even in our sensitive new HST/WFC3 imaging. The absorption component of the OH+ line is blueshifted and peaks near the CO and continuum emission peak while the emission is redshifted and peaks offset by 1.7 kpc from the CO and continuum emission peak, suggesting that the gas is associated with a massive molecular outflow from the intensely star-forming nucleus that supplies 125 Msun/yr of enriched gas to its halo.

preprint2019arXiv

The role of environment in galaxy evolution in the SERVS Survey I: density maps and cluster candidates

We use photometric redshifts derived from new $u$-band through 4.5$μ$m Spitzer IRAC photometry in the 4.8\,deg$^2$ of the XMM-LSS field to construct surface density maps in the redshift range 0.1-1.5. Our density maps show evidence for large-scale structure in the form of filaments spanning several tens of Mpc. Using these maps, we identify 339 overdensities that our simulated lightcone analysis suggests are likely associated with dark matter haloes with masses, $M_{\rm halo}$, log($M_{\rm halo}/M_{\odot})>$13.7. From this list of overdensities we recover 43 of 70 known X-ray detected and spectroscopically confirmed clusters. The missing X-ray clusters are largely at lower redshifts and lower masses than our target log($M_{\rm halo}/M_{\odot})>$13.7. The bulk of the overdensities are compact, but a quarter show extended morphologies which include likely projection effects, clusters embedded in apparent filaments as well as at least one potential cluster merger (at $z\sim1.28$). The strongest overdensity in our highest redshift slice (at $z\sim1.5$) shows a compact red galaxy core potentially implying a massive evolved cluster.

preprint2016arXiv

The multiplicity of 250-$μ$m Herschel sources in the COSMOS field

We investigate the multiplicity of extragalactic sources detected by the Herschel Space Observatory in the COSMOS field. Using 3.6- and 24-$μ$m catalogues, in conjunction with 250-$μ$m data from Herschel, we seek to determine if a significant fraction of Herschel sources are composed of multiple components emitting at 250 $μ$m. We use the XID+ code, using Bayesian inference methods to produce probability distributions of the possible contributions to the observed 250-$μ$m flux for each potential component. The fraction of Herschel flux assigned to the brightest component is highest for sources with total 250-$μ$m fluxes < 45 mJy; however, the flux in the brightest component is still highest in the brightest Herschel sources. The faintest 250-$μ$m sources (30-45 mJy) have the majority of their flux assigned to a single bright component; the second brightest component is typically significantly weaker, and contains the remainder of the 250-$μ$m source flux. At the highest 250-$μ$m fluxes (45-110 mJy), the brightest and second brightest components are assigned roughly equal fluxes, and together are insufficient to reach 100 per cent of the 250-$μ$m source flux. This indicates that additional components are required, beyond the brightest two components, to reproduce the observed flux. 95 per cent of the sources in our sample have a second component that contains more than 10 per cent of the total source flux. Particularly for the brightest Herschel sources, assigning the total flux to a single source may overestimate the flux contributed by around 150 per cent.

preprint2014arXiv

A Measurement of the Millimeter Emission and the Sunyaev-Zel'dovich Effect Associated with Low-Frequency Radio Sources

We present a statistical analysis of the millimeter-wavelength properties of 1.4 GHz-selected sources and a detection of the Sunyaev-Zel'dovich (SZ) effect associated with the halos that host them. The Atacama Cosmology Telescope (ACT) has conducted a survey at 148 GHz, 218 GHz and 277 GHz along the celestial equator. Using samples of radio sources selected at 1.4 GHz from FIRST and NVSS, we measure the stacked 148, 218 and 277 GHz flux densities for sources with 1.4 GHz flux densities ranging from 5 to 200 mJy. At these flux densities, the radio source population is dominated by active galactic nuclei (AGN), with both steep and flat spectrum populations, which have combined radio-to-millimeter spectral indices ranging from 0.5 to 0.95, reflecting the prevalence of steep spectrum sources at high flux densities and the presence of flat spectrum sources at lower flux densities. The thermal SZ effect associated with the halos that host the AGN is detected at the 5$σ$ level through its spectral signature. When we compare the SZ effect with weak lensing measurements of radio galaxies, we find that the relation between the two is consistent with that measured by Planck for local bright galaxies. We present a detection of the SZ effect in some of the lowest mass halos (average $M_{200}\approx10^{13}$M$_{\odot}h_{70}^{-1}$) studied to date. This detection is particularly important in the context of galaxy evolution models, as it confirms that galaxies with radio AGN also typically support hot gaseous halos. With Herschel observations, we show that the SZ detection is not significantly contaminated by dust. We show that 5 mJy$<S_{1.4}<$200 mJy radio sources contribute $\ell(\ell+1)C_{\ell}/(2π)=0.37\pm0.03μ$K$^2$ to the angular power spectrum at $\ell=3000$ at 148 GHz, after accounting for the SZ effect associated with their host halos.

preprint2014arXiv

Detailed modelling of a large sample of Herschel sources in the Lockman Hole: identification of cold dust and of lensing candidates through their anomalous SEDs

We have studied in detail a sample of 967 SPIRE sources with 5-sigma detections at 350 and 500 micron and associations with Spitzer-SWIRE 24 micron galaxies in the HerMES-Lockman survey area, fitting their mid- and far-infrared, and submillimetre, SEDs in an automatic search with a set of six infrared templates. For almost 300 galaxies we have modelled their SEDs individually to ensure the physicality of the fits. We confirm the need for the new cool and cold cirrus templates, and also of the young starburst template, introduced in earlier work. We also identify 109 lensing candidates via their anomalous SEDs and provide a set of colour-redshift constraints which allow lensing candidates to be identified from combined Herschel and Spitzer data. The picture that emerges of the submillimetre galaxy population is complex, comprising ultraluminous and hyperluminous starbursts, lower luminosity galaxies dominated by interstellar dust emission, lensed galaxies and galaxies with surprisingly cold (10-13K) dust. 11 % of 500$μ$m selected sources are lensing candidates. 70 % of the unlensed sources are ultraluminous infrared galaxies and 26 % are hyperluminous. 34 % are dominated by optically thin interstellar dust ('cirrus') emission, but most of these are due to cooler dust than is characteristic of our Galaxy. At the highest infrared luminosities we see SEDs dominated by M82, Arp220 and young starburst types, in roughly equal proportions.

preprint2014arXiv

Multiwavelength characterization of faint Ultra Steep Spectrum radio sources : A search for high-redshift radio galaxies

Ultra Steep Spectrum (USS) radio sources are one of the efficient tracers of powerful High-z Radio Galaxies (HzRGs). In contrast to searches for powerful HzRGs from radio surveys of moderate depths, fainter USS samples derived from deeper radio surveys can be useful in finding HzRGs at even higher redshifts and in unveiling a population of obscured weaker radio-loud AGN at moderate redshifts. Using our 325 MHz GMRT observations (5-sigma ~ 800 microJy) and 1.4 GHz VLA observations (5-sigma ~ 80 - 100 microJy) available in two subfields (viz., VLA-VIMOS VLT Deep Survey (VLA-VVDS) and Subaru X-ray Deep Field (SXDF)) of the XMM-LSS field, we derive a large sample of 160 faint USS radio sources and characterize their nature. The optical, IR counterparts of our USS sample sources are searched using existing deep surveys, at respective wavelengths. We attempt to unveil the nature of our faint USS sources using diagnostic techniques based on mid-IR colors, flux ratios of radio to mid-IR, and radio luminosities. Redshift estimates are available for 86/116 (~ 74%) USS sources in the VLA-VVDS field and for 39/44 (~ 87%) USS sources in the SXDF fields with median values (z_median) ~ 1.18 and ~ 1.57, which are higher than that for non-USS radio sources (z_median non-USS ~ 0.99 and ~ 0.96), in the two subfields, respectively. The flux ratio of radio to mid-IR (S_1.4 GHz/S_3.6 micron) versus redshift diagnostic plot suggests that more than half of our USS sample sources distributed over z ~ 0.5 to 3.8 are likely to be hosted in obscured environments. A significant fraction (~ 26% in the VLA-VVDS and ~ 13% in the SXDF) of our USS sources without redshift estimates mostly remain unidentified in the existing optical, IR surveys, and exhibit high radio to mid-IR flux ratio limits similar to HzRGs, and thus, can be considered as potential HzRG candidates.

preprint2014arXiv

The co-evolution of black hole growth and star formation from a cross-correlation analysis between quasars and the cosmic infrared background

We present the first cross-correlation measurement between Sloan Digital Sky Survey (SDSS) Type 1 quasars and the cosmic infrared background (CIB) measured by Herschel. The distribution of the quasars at 0.15<z<3.5 covers the redshift range where we expect most of the CIB to originate. We detect the sub-mm emission of the quasars, which dominates on small scales, as well as correlated emission from dusty star-forming galaxies (DSFGs) dominant on larger scales. The mean sub-mm flux densities of the DR7 quasars (median redshift <z>=1.4) is $11.1^{+1.6}_{-1.4}$, $7.1^{+1.6}_{-1.3}$ and $3.6^{+1.4}_{-1.0}$ mJy at 250, 350 and 500 microns, respectively, while the mean sub-mm flux densities of the DR9 quasars (<z>=2.5) is $5.7^{+0.7}_{-0.6}$, $5.0^{+0.8}_{-0.7}$ and $1.8^{+0.5}_{-0.4}$ mJy. We find that the correlated sub-mm emission includes both the emission from satellite DSFGs in the same halo as the central quasar and the emission from DSFGs in separate halos (correlated with the quasar-hosting halo). The amplitude of the one-halo term is ~10 times smaller than the sub-mm emission of the quasars, implying the the satellites have a lower star-formation rate than the quasars. The satellite fraction for the DR7 quasars is $0.008^{+0.008}_{-0.005}$ and the host halo mass scale for the central and satellite quasars is $10^{12.36\pm0.87}$ M$_{\odot}$ and $10^{13.60\pm0.38}$ M$_{\odot}$, respectively. The satellite fraction of the DR9 quasars is $0.065^{+0.021}_{-0.031}$ and the host halo mass scale for the central and satellite quasars is $10^{12.29\pm0.62}$ M$_{\odot}$ and $10^{12.82\pm0.39}$ M$_{\odot}$, respectively. Thus, the typical halo environment of the SDSS Type 1 quasars is found to be similar to that of DSFGs, which supports the generally accepted view that dusty starburst and quasar activity are evolutionarily linked.

preprint2012arXiv

Principal Component Analysis and Radiative Transfer modelling of Spitzer IRS Spectra of Ultra Luminous Infrared Galaxies

The mid-infrared spectra of ultraluminous infrared galaxies (ULIRGs) contain a variety of spectral features that can be used as diagnostics to characterise the spectra. However, such diagnostics are biased by our prior prejudices on the origin of the features. Moreover, by using only part of the spectrum they do not utilise the full information content of the spectra. Blind statistical techniques such as principal component analysis (PCA) consider the whole spectrum, find correlated features and separate them out into distinct components. We further investigate the principal components (PCs) of ULIRGs derived in Wang et al.(2011). We quantitatively show that five PCs is optimal for describing the IRS spectra. These five components (PC1-PC5) and the mean spectrum provide a template basis set that reproduces spectra of all z<0.35 ULIRGs within the noise. For comparison, the spectra are also modelled with a combination of radiative transfer models of both starbursts and the dusty torus surrounding active galactic nuclei. The five PCs typically provide better fits than the models. We argue that the radiative transfer models require a colder dust component and have difficulty in modelling strong PAH features. Aided by the models we also interpret the physical processes that the principal components represent. The third principal component is shown to indicate the nature of the dominant power source, while PC1 is related to the inclination of the AGN torus. Finally, we use the 5 PCs to define a new classification scheme using 5D Gaussian mixtures modelling and trained on widely used optical classifications. The five PCs, average spectra for the four classifications and the code to classify objects are made available at: http://www.phys.susx.ac.uk/~pdh21/PCA/

preprint2011arXiv

Deep Spitzer observations of infrared-faint radio sources: high-redshift radio-loud AGN?

Infrared-faint radio sources (IFRSs) are a rare class of object which are relatively bright at radio wavelengths but very faint at infrared and optical wavelengths. Here we present sensitive near-infrared observations of a sample of these sources taken as part of the Spitzer Extragalactic Representative Volume Survey (SERVS). Nearly all the IFRSs are undetected at a level of ~ 1 μ$Jy in these new deep observations, and even the detections are consistent with confusion with unrelated galaxies. A stacked image implies that the median flux density is $S_{3.6μm} ~ 0.2$ μ$Jy or less, giving extreme values of the radio-infrared flux density ratio. Comparison of these objects with known classes of object suggests that the majority are probably high-redshift radio-loud galaxies, possibly suffering from significant dust extinction.

preprint2011arXiv

Specific star-formation and the relation to stellar mass from 0<z<2 as seen in the far-infrared at 70 and 160mu

We use the Spitzer Wide-area InfraRed Extragalactic Legacy Survey (SWIRE) to explore the specific star-formation activity of galaxies and their evolution near the peak of the cosmic far-infrared (FIR) background at 70 and 160um. We use a stacking analysis to determine the mean FIR properties of well defined subsets of galaxies at flux levels well below the FIR catalogue detection limits of SWIRE and other Spitzer surveys. We tabulate the contribution of different subsets of galaxies to the FIR background at 70um and 160um. These long wavelengths provide a good constraint on the bolometric, obscured emission. The large area provides good constraints at low z and in finer redshift bins than previous work. At all redshifts we find that the specific FIR Luminosity (sLFIR) decreases with increasing mass, following a trend L_FIR/M* propto M_* ^beta with beta =-0.38\pm0.14. This is a more continuous change than expected from the {Delucia2007} semi-analytic model suggesting modifications to the feedback prescriptions. We see an increase in the sLFIR by about a factor of ~100 from 0<z<2 and find that the sLFIR evolves as (1+z)^alpha with alpha=4.4\pm0.3 for galaxies with 10.5 < log M*/Msun < 12. This is considerably steeper than the {Delucia2007} semi-analytic model (alpha \sim 2.5). When separating galaxies into early and late types on the basis of the optical/IR spectral energy distributions we find that the decrease in sLFIR with stellar mass is stronger in early type galaxies (beta ~ -0.46), while late type galaxies exhibit a flatter trend (beta \sim -0.15). The evolution is strong for both classes but stronger for the early type galaxies. The early types show a trend of decreasing strength of evolution as we move from lower to higher masses while the evolution of the late type galaxies has little dependence on stellar mass. We suggest that in late-type galaxies we are seeing a consistently declining sSFR..

preprint2010arXiv

A deep radio survey of the AKARI North Ecliptic Pole Field - WSRT 20 cm Radio survey description, observations and data reduction

The Westerbork Radio Synthesis Telescope, WSRT, has been used to make a deep radio survey of an ~ 1.7 sq degree field coinciding with the AKARI North Ecliptic Pole Deep Field. The observations, data reduction and source count analysis are presented, along with a description of the overall scientific objectives. The survey consisted of 10 pointings, mosaiced with enough overlap to maintain a similar sensitivity across the central region that reached as low as 21 microJy per beam at 1.4 GHz. A catalogue containing 462 sources detected with a resolution of 17"x15" is presented. The differential source counts calculated from the WSRT data have been compared with those from the shallow VLA-NEP survey of Kollgaard et al 1994, and show a pronounced excess for sources fainter than ~ 1 mJy, consistent with the presence of a population of star forming galaxies at sub-mJy flux levels. The AKARI North Ecliptic Pole Deep field is the focus of a major observing campaign conducted across the entire spectral region. The combination of these data sets, along with the deep nature of the radio observations will allow unique studies of a large range of topics including the redshift evolution of the luminosity function of radio sources, the clustering environment of radio galaxies, the nature of obscured radio-loud active galactic nuclei, and the radio/far-infrared correlation for distant galaxies. This catalogue provides the basic data set for a future series of paper dealing with source identifications, morphologies, and the associated properties of the identified radio sources.

preprint2010arXiv

A Halo Model of Local IRAS Galaxies Selected at 60 Micron Using Conditional Luminosity Functions

Using conditional luminosity functions (CLFs) which encode the luminosity distribution of galaxies as a function of halo mass, we construct a halo model of IRAS galaxies selected at 60 micron. An abundance matching technique is used to link galaxy luminosity to the host halo mass. The shape of the mass - light relation at 60 micron is different from those derived at r-, K- and B-band. This is because the 60 micron LF can not be fitted by a Schechter function with a sharp exponential cutoff. We then seek the parameters in the CLFs that best fit the LF and power spectrum. We find that the predicted galaxy bias as a function of L60 from the best-fit model agrees well with the clustering measurements. At the faint end of the LF where quiescent star-forming galaxies dominate, most IRAS galaxies are central galaxies in halos of M >~ 10^{10} h^{-1} M_sun but a non-negligible fraction are satellites typically hosted in more massive halos. The majority of IRAS galaxies with L60 >~ 10^{10} h^{-2} L_sun are M82 type starbursts which are central galaxies hosted in halos of M >~ 10^{12.5} h^{-1} M_sun. In comparison, optical galaxies generally reside in much more massive halos. The rate of change in L60 (an indicator of recent star formation) as a function of halo mass at M >~ 10^{12.5} h^{-1} M_sun is much larger that d L_optical / dM or d L_NIR / dM indicating the existence of physical mechanisms which are very efficient in converting cold gas into stars, possibly dynamical effects arising from interactions or mergers. We further calculate the space density of major mergers for halos massive enough to host ultraluminous infrared galaxies (ULIRGs) using the mean merger rate derived from the Millennium simulations. Compared to the space density of local ULIRGs, it implies that either the majority of major mergers at z~0 do not lead to ULIRGs or the ULIRG phase is relatively short.

preprint2010arXiv

Principal component analysis of the Spitzer IRS spectra of ultraluminous infrared galaxies

We present the first principal component analysis (PCA) applied to a sample of 119 Spitzer Infrared Spectrograph (IRS) spectra of local ultraluminous infrared galaxies (ULIRGs) at z<0.35. The purpose of this study is to objectively and uniquely characterise the local ULIRG population using all information contained in the observed spectra. We have derived the first three principal components (PCs) from the covariance matrix of our dataset which account for over 90% of the variance. The first PC is characterised by dust temperatures and the geometry of the mix of source and dust. The second PC is a pure star formation component. The third PC represents an anti-correlation between star formation activity and a rising AGN. Using the first three PCs, we are able to accurately reconstruct most of the spectra in our sample. Our work shows that there are several factors that are important in characterising the ULIRG population, dust temperature, geometry, star formation intensity, AGN contribution, etc. We also make comparison between PCA and other diagnostics such as ratio of the 6.2 microns PAH emission feature to the 9.7 micron silicate absorption depth and other observables such as optical spectral type.

preprint2010arXiv

The Herschel-SPIRE Legacy Survey (HSLS): the scientific goals of a shallow and wide submillimeter imaging survey with SPIRE

A large sub-mm survey with Herschel will enable many exciting science opportunities, especially in an era of wide-field optical and radio surveys and high resolution cosmic microwave background experiments. The Herschel-SPIRE Legacy Survey (HSLS), will lead to imaging data over 4000 sq. degrees at 250, 350, and 500 micron. Major Goals of HSLS are: (a) produce a catalog of 2.5 to 3 million galaxies down to 26, 27 and 33 mJy (50% completeness; 5 sigma confusion noise) at 250, 350 and 500 micron, respectively, in the southern hemisphere (3000 sq. degrees) and in an equatorial strip (1000 sq. degrees), areas which have extensive multi-wavelength coverage and are easily accessible from ALMA. Two thirds of the of the sources are expected to be at z > 1, one third at z > 2 and about a 1000 at z > 5. (b) Remove point source confusion in secondary anisotropy studies with Planck and ground-based CMB data. (c) Find at least 1200 strongly lensed bright sub-mm sources leading to a 2% test of general relativity. (d) Identify 200 proto-cluster regions at z of 2 and perform an unbiased study of the environmental dependence of star formation. (e) Perform an unbiased survey for star formation and dust at high Galactic latitude and make a census of debris disks and dust around AGB stars and white dwarfs.

preprint2009arXiv

A New Approach to Multi-wavelength Associations of Astronomical Sources

One of the biggest problems faced by current and next-generation astronomical surveys is trying to produce large numbers of accurate cross identifications across a range of wavelength regimes with varying data quality and positional uncertainty. Until recently simple spatial "nearest neighbour" associations have been sufficient for most applications. However as advances in instrumentation allow more sensitive images to be made the rapid increase in the source density has meant that source confusion across multiple wavelengths is a serious problem. The field of far-IR and sub-mm astronomy has been particularly hampered by such problems. The poor angular resolution of current sub-mm and far-IR instruments is such that in a lot of cases there are multiple plausible counterparts for each source at other wavelengths. Here we present a new automated method of producing associations between sources at different wavelengths using a combination of spatial and SED information set in the Bayesian framework presented by Budavari & Szalay (2008). Testing of the technique is performed on both simulated catalogues of sources from GaLICS and real data from multi-wavelength observations of the SXDF. It is found that a single figure of merit, the Bayes factor, can be effectively used to describe the confidence in the match. Further applications of this technique to future Herschel datasets are discussed.

preprint2009arXiv

GAMA: towards a physical understanding of galaxy formation

The Galaxy And Mass Assembly (GAMA) project is the latest in a tradition of large galaxy redshift surveys, and is now underway on the 3.9m Anglo-Australian Telescope at Siding Spring Observatory. GAMA is designed to map extragalactic structures on scales of 1kpc - 1Mpc in complete detail to a redshift of z~0.2, and to trace the distribution of luminous galaxies out to z~0.5. The principal science aim is to test the standard hierarchical structure formation paradigm of Cold Dark Matter (CDM) on scales of galaxy groups, pairs, discs, bulges and bars. We will measure (1) the Dark Matter Halo Mass Function (as inferred from galaxy group velocity dispersions); (2) baryonic processes, such as star formation and galaxy formation efficiency (as derived from Galaxy Stellar Mass Functions); and (3) the evolution of galaxy merger rates (via galaxy close pairs and galaxy asymmetries). Additionally, GAMA will form the central part of a new galaxy database, which aims to contain 275,000 galaxies with multi-wavelength coverage from coordinated observations with the latest international ground- and space-based facilities: GALEX, VST, VISTA, WISE, HERSCHEL, GMRT and ASKAP. Together, these data will provide increased depth (over 2 magnitudes), doubled spatial resolution (0.7"), and significantly extended wavelength coverage (UV through Far-IR to radio) over the main SDSS spectroscopic survey for five regions, each of around 50 deg^2. This database will permit detailed investigations of the structural, chemical, and dynamical properties of all galaxy types, across all environments, and over a 5 billion year timeline.

preprint2007arXiv

Bayesian methods of astronomical source extraction

We present two new source extraction methods, based on Bayesian model selection and using the Bayesian Information Criterion (BIC). The first is a source detection filter, able to simultaneously detect point sources and estimate the image background. The second is an advanced photometry technique, which measures the flux, position (to sub-pixel accuracy), local background and point spread function. We apply the source detection filter to simulated Herschel-SPIRE data and show the filter's ability to both detect point sources and also simultaneously estimate the image background. We use the photometry method to analyse a simple simulated image containing a source of unknown flux, position and point spread function; we not only accurately measure these parameters, but also determine their uncertainties (using Markov-Chain Monte Carlo sampling). The method also characterises the nature of the source (distinguishing between a point source and extended source). We demonstrate the effect of including additional prior knowledge. Prior knowledge of the point spread function increase the precision of the flux measurement, while prior knowledge of the background has onlya small impact. In the presence of higher noise levels, we show that prior positional knowledge (such as might arise from a strong detection in another waveband) allows us to accurately measure the source flux even when the source is too faint to be detected directly. These methods are incorporated in SUSSEXtractor, the source extraction pipeline for the forthcoming Akari FIS far-infrared all-sky survey. They are also implemented in a stand-alone, beta-version public tool that can be obtained at http://astronomy.sussex.ac.uk/$\sim$rss23/sourceMiner\_v0.1.2.0.tar.gz

preprint2006arXiv

The spatial clustering of ultraluminous infrared galaxies over 1.5<z<3

We present measurements of the spatial clustering of galaxies with stellar masses >10^11Msun, infrared luminosities >10^12 Lsun, and star formation rates >200Msun per year in two redshift intervals; 1.5<z<2.0 and 2<z<3. Both samples cluster very strongly, with spatial correlation lengths of r_0=14.40+/-1.99 h^-1Mpc for the 2<z<3 sample, and r_0=9.40+/-2.24 h^-1Mpc for the 1.5<z<2.0 sample. These clustering amplitudes are consistent with both populations residing in dark matter haloes with masses of ~6x10^13Msun, making them among the most biased galaxies at these epochs. We infer, from this and previous results, that a minimum dark matter halo mass is an important factor for all forms of luminous, obscured activity in galaxies at z>1, both starbursts and AGN. Adopting plausible models for the growth of DM haloes with redshift, then the haloes hosting the 2<z<3 sample will likely host the richest clusters of galaxies at z=0, whereas the haloes hosting the 1.5<z<2.0 sample will likely host poor to rich clusters at z=0. We conclude that ULIRGs at z>1 signpost stellar buildup in galaxies that will reside in clusters at z=0, with ULIRGs at increasing redshifts signposting the buildup of stars in galaxies that will reside in increasingly rich clusters.

preprint1997arXiv

Observations of the Hubble Deep Field with the Infrared Space Observatory. III. Source counts and P(D) analysis

We present source counts at 6.7 micron and 15 micron from our maps of the Hubble Deep Field region, reaching 38.6 microJy at 6.7 micron and 255 microJy at 15 micron. These are the first ever extra-galactic number counts to be presented at 6.7 micron and are 3 decades fainter than IRAS at 12 micron. Both source counts and a P(D) analysis suggest we have reached the ISO confusion limit at 15 micron: this will have important implications for future space missions. These data provide an excellent reference point for other ongoing ISO surveys. A no-evolution model at 15 micron is ruled out at >3 sigma, while two models which fit the steep IRAS 60 micron counts are acceptable. This provides important confirmation of the strong evolution seen in IRAS surveys. One of these models can then be ruled out from the 6.7 micron data.

preprint1996arXiv

The European Large Area ISO Survey: ELAIS

I describe a European collaborative project to survey \sim 20 square degrees of the sky at 15\micron and 90\micron with ISO. This is the largest open time project being undertaken by ISO. The depth and areal coverage were designed to complement the various Guaranteed Time surveys. The main science thrust is to explore star formation in galaxies to a much higher redshift than was probed by IRAS. We expect to detect around 8000 extra-galactic objects and a similar number of Galactic sources. The maps and source catalogues will represent a major legacy from ISO, inspiring follow up work for many years to come.