-
PolyCLEAN: When Högbom meets Bayes -- Fast Super-Resolution Imaging with Bayesian MAP Estimation
Authors:
Adrian Jarret,
Sepand Kashani,
Joan Rué-Queralt,
Paul Hurley,
Julien Fageot,
Matthieu Simeoni
Abstract:
Aims: We address two issues for the adoption of convex optimization in radio interferometric imaging. First, a method for a fine resolution setup is proposed which scales naturally in terms of memory usage and reconstruction speed. Second, a new tool to localize a region of uncertainty is developed, paving the way for quantitative imaging in radio interferometry. Methods: The classical $\ell_1$ pe…
▽ More
Aims: We address two issues for the adoption of convex optimization in radio interferometric imaging. First, a method for a fine resolution setup is proposed which scales naturally in terms of memory usage and reconstruction speed. Second, a new tool to localize a region of uncertainty is developed, paving the way for quantitative imaging in radio interferometry. Methods: The classical $\ell_1$ penalty is used to turn the inverse problem into a sparsity-promoting optimization. For efficient implementation, the so-called Frank-Wolfe algorithm is used together with a \textit{polyatomic} refinement. The algorithm naturally produces sparse images at each iteration, leveraged to reduce memory and computational requirements. In that regard, PolyCLEAN reproduces the numerical behavior of CLEAN while guaranteeing that it solves the minimization problem of interest. Additionally, we introduce the dual certificate image, which appears as a numerical byproduct of the Frank-Wolfe algorithm. This image is proposed as a tool for uncertainty quantification on the location of the recovered sources. Results: PolyCLEAN demonstrates good scalability performance, in particular for fine-resolution grids. On simulations, the Python-based implementation is competitive with the fast numerically-optimized CLEAN solver. This acceleration does not affect image reconstruction quality: PolyCLEAN images are consistent with CLEAN-obtained ones for both point sources and diffuse emission recovery. We also highlight PolyCLEAN reconstruction capabilities on observed radio measurements. Conclusions: PolyCLEAN can be considered as an alternative to CLEAN in the radio interferometric imaging pipeline, as it enables the use of Bayesian priors without impacting the scalability and numerical performance of the imaging method.
△ Less
Submitted 3 June, 2024;
originally announced June 2024.
-
Probabilistic and progressive deblended far-infrared and sub-millimetre point source catalogues I. Methodology and first application in the COSMOS field
Authors:
Lingyu Wang,
Antonio La Marca,
Fangyou Gao,
William J. Pearson,
Berta Margalef-Bentabol,
Matthieu Béthermin,
Longji Bing,
James Donnellan,
Peter D. Hurley,
Seb J. Oliver,
Catherine L. Hale,
Matt J. Jarvis,
Lucia Marchetti,
Mattia Vaccari,
Imogen H. Whittam
Abstract:
Single-dish far-infrared (far-IR) and sub-millimetre (sub-mm) point source catalogues and their connections with catalogues at other wavelengths are of paramount importance. However, due to the large mismatch in spatial resolution, cross-matching galaxies at different wavelengths is challenging. This work aims to develop the next-generation deblended far-IR and sub-mm catalogues and present the fi…
▽ More
Single-dish far-infrared (far-IR) and sub-millimetre (sub-mm) point source catalogues and their connections with catalogues at other wavelengths are of paramount importance. However, due to the large mismatch in spatial resolution, cross-matching galaxies at different wavelengths is challenging. This work aims to develop the next-generation deblended far-IR and sub-mm catalogues and present the first application in the COSMOS field. Our progressive deblending used the Bayesian probabilistic framework known as XID+. The deblending started from the Spitzer/MIPS 24 micron data, using an initial prior list composed of sources selected from the COSMOS2020 catalogue and radio catalogues from the VLA and the MeerKAT surveys, based on spectral energy distribution modelling which predicts fluxes of the known sources at the deblending wavelength. To speed up flux prediction, we made use of a neural network-based emulator. After deblending the 24 micron data, we proceeded to the Herschel PACS (100 & 160 micron) and SPIRE wavebands (250, 350 & 500 micron). Each time we constructed a tailor-made prior list based on the predicted fluxes of the known sources. Using simulated far-IR and sub-mm sky, we detailed the performance of our deblending pipeline. After validation with simulations, we then deblended the real observations from 24 to 500 micron and compared with blindly extracted catalogues and previous versions of deblended catalogues. As an additional test, we deblended the SCUBA-2 850 micron map and compared our deblended fluxes with ALMA measurements, which demonstrates a higher level of flux accuracy compared to previous results.We publicly release our XID+ deblended point source catalogues. These deblended long-wavelength data are crucial for studies such as deriving the fraction of dust-obscured star formation and better separation of quiescent galaxies from dusty star-forming galaxies.
△ Less
Submitted 28 May, 2024;
originally announced May 2024.
-
BIPP: An efficient HPC implementation of the Bluebild algorithm for radio astronomy
Authors:
Emma Tolley,
Simon Frasch,
Etienne Orliac,
Shreyam Krishna,
Michele Bianco,
Sepand Kashani,
Paul Hurley,
Matthieu Simeoni,
Jean-Paul Kneib
Abstract:
The Bluebild algorithm is a new technique for image synthesis in radio astronomy which forms a least-squares estimate of the sky intensity function using the theory of sampling and interpolation operators. We present an HPC implementation of the Bluebild algorithm for radio-interferometric imaging: Bluebild Imaging++ (BIPP). BIPP is a spherical imager that leverages functional PCA to decompose the…
▽ More
The Bluebild algorithm is a new technique for image synthesis in radio astronomy which forms a least-squares estimate of the sky intensity function using the theory of sampling and interpolation operators. We present an HPC implementation of the Bluebild algorithm for radio-interferometric imaging: Bluebild Imaging++ (BIPP). BIPP is a spherical imager that leverages functional PCA to decompose the sky into distinct energy levels. The library features interfaces to C++, C and Python and is designed with seamless GPU acceleration in mind. We evaluate the accuracy and performance of BIPP on simulated observations of the upcoming Square Kilometer Array Observatory and real data from the Low-Frequency Array (LOFAR) telescope. We find that BIPP offers accurate wide-field imaging with no need for a w-term approximation and has comparable execution time with respect to the interferometric imaging libraries CASA and WSClean. Futhermore, due to the energy level decomposition, images produced with BIPP can reveal information about faint and diffuse structures before any cleaning iterations. The source code of BIPP is publicly released.
△ Less
Submitted 13 October, 2023;
originally announced October 2023.
-
Getting ready for the LSST data -- estimating the physical properties of $z<2.5$ main sequence galaxies
Authors:
Gabriele Riccio,
Katarzyna Małek,
Ambra Nanni,
Mederic Boquien,
Veronique Buat,
Denis Burgarella,
Darko Donevski,
Mahmoud Hamed,
Peter Hurley,
Raphael Shirley,
Agnieszka Pollo
Abstract:
In this work we study how to employ the upcoming Legacy Survey of Space and Time (LSST) data to constrain physical properties of normal, star forming galaxies. We use simulated LSST data and existing real observations to test the estimations of the physical properties of galaxies, such as star formation rate (SFR), stellar mass ($M_{star}$), and dust luminosity ($L_{dust}$). We focus on normal sta…
▽ More
In this work we study how to employ the upcoming Legacy Survey of Space and Time (LSST) data to constrain physical properties of normal, star forming galaxies. We use simulated LSST data and existing real observations to test the estimations of the physical properties of galaxies, such as star formation rate (SFR), stellar mass ($M_{star}$), and dust luminosity ($L_{dust}$). We focus on normal star-forming galaxies, as they form the majority of the galaxy population in the universe and therefore are more likely to be observed by the LSST. We perform a simulation of LSST observations and uncertainties of 50,385 real galaxies within redshift range $0<z<2.5$. In order to achieve this goal, we used the unique multi-wavelength data from the Herschel Extragalactic Legacy Project (HELP) survey. Our analysis focus on two fields: ELAIS-N1 and COSMOS. To obtain galaxy physical properties we fit their Spectral Energy Distributions (SEDs) using the Code Investigating GALaxy Emission (CIGALE). We compare the main galaxy physical properties obtained from the fit of the observed multi-wavelength photometry of galaxies (from UV to FIR) to the ones obtained from the simulated LSST optical measurements only. The stellar masses estimated based on the LSST measurements are in agreement with the full UV-FIR SED estimations, as they depend mainly on the UV and optical emission, well covered by LSST in the considered redshift range. We obtain a clear overestimation of SFR, $L_{dust}$, $M_{dust}$ estimated with LSST only, highly correlated with redshift. We investigate the cause of this overestimation and we conclude that it is related to an overestimation of the dust attenuation, both UV and NIR. We find that it is necessary to employ auxiliary rest-frame mid-infrared observations, simulated UV observations, or FUV attenuation (AFUV)- Mstar relation, to correct the overestimation.
△ Less
Submitted 23 June, 2021;
originally announced June 2021.
-
Deep Extragalactic VIsible Legacy Survey (DEVILS): Consistent multi-wavelength photometry for the DEVILS regions (COSMOS, XMMLSS & ECDFS)
Authors:
L. J. M. Davies,
J. E. Thorne,
A. S. G. Robotham,
S. Bellstedt,
S. P. Driver,
N. J. Adams,
M. Bilicki,
R. A. A. Bowler,
M. Bravo,
L. Cortese,
C. Foster,
M. W. Grootes,
B. Häußler,
A. Hashemizadeh,
B. W. Holwerda,
P. Hurley,
M. J. Jarvis,
C. Lidman,
N. Maddox,
M. Meyer,
M. Paolillo,
S. Phillipps,
M. Radovich,
M. Siudek,
M. Vaccari
, et al. (1 additional authors not shown)
Abstract:
The Deep Extragalactic VIsible Legacy Survey (DEVILS) is an ongoing high-completeness, deep spectroscopic survey of $\sim$60,000 galaxies to Y$<$21.2 mag, over $\sim$6 deg2 in three well-studied deep extragalactic fields: D10 (COSMOS), D02 (XMM-LSS) and D03 (ECDFS). Numerous DEVILS projects all require consistent, uniformly-derived and state-of-the-art photometric data with which to measure galaxy…
▽ More
The Deep Extragalactic VIsible Legacy Survey (DEVILS) is an ongoing high-completeness, deep spectroscopic survey of $\sim$60,000 galaxies to Y$<$21.2 mag, over $\sim$6 deg2 in three well-studied deep extragalactic fields: D10 (COSMOS), D02 (XMM-LSS) and D03 (ECDFS). Numerous DEVILS projects all require consistent, uniformly-derived and state-of-the-art photometric data with which to measure galaxy properties. Existing photometric catalogues in these regions either use varied photometric measurement techniques for different facilities/wavelengths leading to inconsistencies, older imaging data and/or rely on source detection and photometry techniques with known problems. Here we use the ProFound image analysis package and state-of-the-art imaging datasets (including Subaru-HSC, VST-VOICE, VISTA-VIDEO and UltraVISTA-DR4) to derive matched-source photometry in 22 bands from the FUV to 500μm. This photometry is found to be consistent, or better, in colour-analysis to previous approaches using fixed-size apertures (which are specifically tuned to derive colours), but produces superior total source photometry, essential for the derivation of stellar masses, star-formation rates, star-formation histories, etc. Our photometric catalogue is described in detail and, after internal DEVILS team projects, will be publicly released for use by the broader scientific community.
△ Less
Submitted 11 June, 2021;
originally announced June 2021.
-
HELP: The Herschel Extragalactic Legacy Project
Authors:
R. Shirley,
K. Duncan,
M. C. Campos Varillas,
P. D. Hurley,
K. Malek,
Y. Roehlly,
M. W. L. Smith,
H. Aussel,
T. Bakx,
V. Buat,
D. Burgarella,
N. Christopher,
S. Duivenvoorden,
S. Eales,
A. Efstathiou,
E. A. Gonzalez Solares,
M. Griffin,
M. Jarvis,
B. Lo Faro,
L. Marchetti,
I. McCheyne,
A. Papadopoulos,
K. Penner,
E. Pons,
M. Prescott
, et al. (7 additional authors not shown)
Abstract:
We present the Herschel Extragalactic Legacy Project (HELP). This project collates, curates, homogenises, and creates derived data products for most of the premium multi-wavelength extragalactic data sets. The sky boundaries for the first data release cover 1270 deg2 defined by the Herschel SPIRE extragalactic survey fields; notably the Herschel Multi-tiered Extragalactic Survey (HerMES) and the H…
▽ More
We present the Herschel Extragalactic Legacy Project (HELP). This project collates, curates, homogenises, and creates derived data products for most of the premium multi-wavelength extragalactic data sets. The sky boundaries for the first data release cover 1270 deg2 defined by the Herschel SPIRE extragalactic survey fields; notably the Herschel Multi-tiered Extragalactic Survey (HerMES) and the Herschel Atlas survey (H-ATLAS). Here, we describe the motivation and principal elements in the design of the project. Guiding principles are transparent or "open" methodologies with care for reproducibility and identification of provenance. A key element of the design focuses around the homogenisation of calibration, meta data and the provision of information required to define the selection of the data for statistical analysis. We apply probabilistic methods that extract information directly from the images at long wavelengths, exploiting the prior information available at shorter wavelengths and providing full posterior distributions rather than maximum likelihood estimates and associated uncertainties as in traditional catalogues. With this project definition paper we provide full access to the first data release of HELP; Data Release 1 (DR1), including a monolithic map of the largest SPIRE extragalactic field at 385 deg2 and 18 million measurements of PACS and SPIRE fluxes. We also provide tools to access and analyse the full HELP database. This new data set includes far-infrared photometry, photometric redshifts, and derived physical properties estimated from modelling the spectral energy distributions.
△ Less
Submitted 12 May, 2021;
originally announced May 2021.
-
A hyperluminous obscured quasar at a redshift of z ~ 4.3
Authors:
Andreas Efstathiou,
Katarzyna Malek,
Denis Burgarella,
Peter Hurley,
Seb Oliver,
Veronique Buat,
Raphael Shirley,
Steven Duivenvoorden,
Vicky Papadopoulou Lesta,
Duncan Farrah,
Kenneth J. Duncan,
Maria del Carmen Campos Varillas
Abstract:
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 a…
▽ More
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\%$.
△ Less
Submitted 8 April, 2022; v1 submitted 14 March, 2021;
originally announced March 2021.
-
Have we seen all the galaxies that comprise the cosmic infrared background at 250\,$μ$m $\le λ\le$ 500\,$μ$m?
Authors:
S. Duivenvoorden,
S. Oliver,
M. Bethermin,
D. L. Clements,
G. De Zotti,
A. Efstathiou,
D. Farrah,
P. D. Hurley,
R. J. Ivison,
G. Lagache,
D. Scott,
R. Shirley,
L. Wang,
M. Zemcov
Abstract:
The cosmic infrared background (CIB) provides a fundamental observational constraint on the star-formation history of galaxies over cosmic history. We estimate the contribution to the CIB from catalogued galaxies in the COSMOS field by using a novel map fitting technique on the \textit{Herschel} SPIRE maps. Prior galaxy positions are obtained using detections over a large range in wavelengths in t…
▽ More
The cosmic infrared background (CIB) provides a fundamental observational constraint on the star-formation history of galaxies over cosmic history. We estimate the contribution to the CIB from catalogued galaxies in the COSMOS field by using a novel map fitting technique on the \textit{Herschel} SPIRE maps. Prior galaxy positions are obtained using detections over a large range in wavelengths in the $K_{\rm s}$--3\,GHz range. Our method simultaneously fits the galaxies, the system foreground, and the leakage of flux from galaxies located in masked areas and corrects for an "over-fitting" effect not previously accounted for in stacking methods. We explore the contribution to the CIB as a function of galaxy survey wavelength and depth. We find high contributions to the CIB with the deep $r$ ($m_{\rm AB} \le 26.5$), $K_{\rm s}$ ($m_{\rm AB} \le 24.0$) and 3.6\,$μ$m ($m_{\rm AB} \le 25.5$) catalogues. We combine these three deep catalogues and find a total CIB contributions of 10.5 $\pm$ 1.6, 6.7 $\pm$ 1.5 and 3.1 $\pm$ 0.7\,nWm$^{-2}$sr$^{-1}$ at 250, 350 and 500\,$μ$m, respectively. Our CIB estimates are consistent with recent phenomenological models, prior based SPIRE number counts and with (though more precise than) the diffuse total measured by FIRAS. Our results raise the interesting prospect that the CIB contribution at $λ\le 500\,μ$m from known galaxies has converged. Future large-area surveys like those with the Large Synoptic Survey Telescope are therefore likely to resolve a substantial fraction of the population responsible for the CIB at 250\,$μ$m $\leq λ\leq$ 500\,$μ$m.
△ Less
Submitted 4 November, 2019;
originally announced November 2019.
-
HELP: A catalogue of 170 million objects, selected at 0.36-4.5 $μ$m, from 1270 deg.$^{2}$ of prime extragalactic fields
Authors:
Raphael Shirley,
Yannick Roehlly,
Peter D Hurley,
Veronique Buat,
María del Carmen Campos Varillas,
Steven Duivenvoorden,
Kenneth J Duncan,
Andreas Efstathiou,
Duncan Farrah,
Eduardo González Solares,
Katarzyna Małek,
Lucia Marchetti,
Ian McCheyne,
Andreas Papadopoulos,
Estelle Pons,
Roberto Scipioni,
Mattia Vaccari,
Seb Oliver
Abstract:
We present an optical to near-infrared selected astronomical catalogue covering 1270 deg.$^2$. This is the first attempt to systematically combine data from 23 of the premier extragalactic survey fields - the product of a vast investment of telescope time. The fields are those imaged by the Herschel Space Observatory which form the Herschel Extragalactic Legacy Project (HELP). Our catalogue of 170…
▽ More
We present an optical to near-infrared selected astronomical catalogue covering 1270 deg.$^2$. This is the first attempt to systematically combine data from 23 of the premier extragalactic survey fields - the product of a vast investment of telescope time. The fields are those imaged by the Herschel Space Observatory which form the Herschel Extragalactic Legacy Project (HELP). Our catalogue of 170 million objects is constructed by a positional cross match of 51 public surveys. This high resolution optical, near-infrared, and mid-infrared catalogue is designed for photometric redshift estimation, extraction of fluxes in lower resolution far-infrared maps, and spectral energy distribution modelling. It collates, standardises, and provides value added derived quantities including corrected aperture magnitudes and astrometry correction over the Herschel extragalactic wide fields for the first time. $grizy$ fluxes are available on all fields with $g$ band data reaching $5σ$ point-source depths in a 2 arcsec aperture of 23.5, 24.4, and 24.6 (AB) mag at the 25th, 50th, and 75th percentiles, by area covered, across all HELP fields. It has $K$ or $K_s$ coverage over 1146 deg.$^2$ with depth percentiles of 20.2, 20.4, and 21.0 mag respectively. The IRAC Ch 1 band is available over 273 deg.$^2$ with depth percentiles of 17.7, 21.4, and 22.2 mag respectively. This paper defines the "masterlist" objects for the first data release (DR1) of HELP. This large sample of standardised total and corrected aperture fluxes, uniform quality flags, and completeness measures provides large well understood statistical samples over the full Herschel extragalactic sky.
△ Less
Submitted 22 October, 2019; v1 submitted 9 September, 2019;
originally announced September 2019.
-
A multi-wavelength de-blended Herschel view of the statistical properties of dusty star-forming galaxies across cosmic time
Authors:
L. Wang,
W. J. Pearson,
W. Cowley,
J. W. Trayford,
M. Bethermin,
C. Gruppioni,
P. Hurley,
M. J. Michalowski
Abstract:
We aim to study the statistical properties of dusty star-forming galaxies, such as their number counts, luminosity functions (LF) and dust-obscured star-formation rate density (SFRD). We use state-of-the-art de-blended Herschel catalogue in the COSMOS field, generated by combining the Bayesian source extraction tool XID+ and informative prior on the spectral energy distributions, to measure the nu…
▽ More
We aim to study the statistical properties of dusty star-forming galaxies, such as their number counts, luminosity functions (LF) and dust-obscured star-formation rate density (SFRD). We use state-of-the-art de-blended Herschel catalogue in the COSMOS field, generated by combining the Bayesian source extraction tool XID+ and informative prior on the spectral energy distributions, to measure the number counts and LFs at far-infrared (FIR) and sub-millimetre (sub-mm) wavelengths.
Thanks to our de-confusion technique and deep multi-wavelength photometry, we are able to achieve more accurate measurements while probing ten times below the confusion limit. Our number counts at 250 microns agree well with previous Herschel studies. However, our counts at 350 and 500 microns are considerably below previous Herschel results. This is due to previous studies suffering from source confusion which worsens towards longer wavelength. Our number counts at 450 and 870 microns show excellent agreement with previous determinations derived from single dish and interferometric observations. Our measurements of the monochromatic LF and the total IR LF agree well with previous results. The increased dynamic range of our measurements allows us to better measure the faint-end slope of the LF and measure the dust-obscured SFRD out to z~6. We find that the fraction of dust obscured star-formation activity is at its highest around z~1 which then decreases towards both low and high redshift. We do not find a shift of balance between z~3 and z~4 in the cosmic star-formation history from being dominated by unobscured star formation at higher redshift to obscured star formation at lower redshift. However, we do find 3<z<4 to be an interesting transition period as the fraction of the total SFRD that is obscured by dust is significantly lower at higher redshifts.
△ Less
Submitted 25 February, 2019;
originally announced February 2019.
-
HELP: modelling the spectral energy distributions of Herschel detected galaxies in the ELAIS N1 field
Authors:
K. Malek,
V. Buat,
Y. Roehlly,
D. Burgarella,
P. D. Hurley,
R. Shirley,
K. Duncan,
A. Efstathiou,
A. Papadopoulos,
M. Vaccari,
D. Farrah,
L. Marchetti,
S. Oliver
Abstract:
The HELP project focuses on the data from ESA's Herschel mission, which covered over 1300$deg^2$ and is preparing to publish a multi-wavelength catalogue of millions of objects. Our main goal is to find the best approach to simultaneously fitting SEDs of millions of galaxies across a wide redshift range to obtain homogeneous estimates of the main physical parameters of detected IR galaxies. We per…
▽ More
The HELP project focuses on the data from ESA's Herschel mission, which covered over 1300$deg^2$ and is preparing to publish a multi-wavelength catalogue of millions of objects. Our main goal is to find the best approach to simultaneously fitting SEDs of millions of galaxies across a wide redshift range to obtain homogeneous estimates of the main physical parameters of detected IR galaxies. We perform SED fitting on the UV/NIR to FIR emission of 42 047 galaxies from the pilot HELP field: ELAIS N1. We use the latest release of CIGALE, a galaxy SED fitting code relying on energy balance, to deliver the main physical parameters such as $M_{star}$, SFR, and $L_{dust}$. We implement additional quality criteria to the fits by calculating $χ^2$ values for the stellar and dust part of the spectra independently. These criteria allow us to identify the best fits and to identify peculiar galaxies. We perform the SED fitting by assuming three different dust attenuation laws separately allowing us to test the impact of the assumed law on estimated physical parameters. We implemented two additional quality value checks for the SED fitting method based on $M_{star}$ estimation and energy budget. This method allows us to identify possible objects with incorrect matching in the catalogue and peculiar galaxies; we found 351 possible candidates of lensed galaxies using two complementary $χ^2$s criteria (stellar and IR) and $z_{phot}$ calculated for the IR part of the spectrum only. We find that the attenuation law has an important impact on the $M_{star}$ estimate (on average leading to disparities of a factor of two). We derive the relation between $M_{star}$ estimates obtained by different attenuation laws and we find the best recipe for our sample. We also make independent estimates of the total $L_{dust}$ parameter from stellar emission by fitting the galaxies with and without IR data separately.
△ Less
Submitted 3 September, 2018;
originally announced September 2018.
-
Cosmic happenstance: 24-$μ$m selected, multi-component Herschel sources are line-of-sight projections
Authors:
Jillian M. Scudder,
Seb Oliver,
Peter D. Hurley,
Julie L. Wardlow,
Lingyu Wang,
Duncan Farrah
Abstract:
In this paper, we investigate the physical associations between blended far-infrared (FIR)-emitting galaxies, in order to identify the level of line-of-sight projection contamination in the single-dish Herschel data. Building on previous work, and as part of the Herschel Extragalactic Legacy Project (HELP), we identify a sample of galaxies in the COSMOS field which are found to be both FIR-bright…
▽ More
In this paper, we investigate the physical associations between blended far-infrared (FIR)-emitting galaxies, in order to identify the level of line-of-sight projection contamination in the single-dish Herschel data. Building on previous work, and as part of the Herschel Extragalactic Legacy Project (HELP), we identify a sample of galaxies in the COSMOS field which are found to be both FIR-bright (typically $\sim 15$ mJy) and blended within the Herschel 250 $μ$m beam. We identify a spectroscopic or photometric redshift for each FIR-bright source. We conduct a joint probability distribution analysis on the redshift probability density functions to determine the fraction of the FIR sources with multiple FIR-bright counterparts which are likely to be found at consistent ($Δz$ $< 0.01$) redshifts. We find that only 3 (0.4 per cent) of the pair permutations between counterparts are $>50$ per cent likely to be at consistent redshifts. A majority of counterparts (72 per cent) have no overlap in their redshift probability distributions whatsoever. This is in good agreement with the results of recent simulations, which indicate that single-dish observations of the FIR sky should be strongly contaminated by line of sight projection effects. We conclude that for our sample of 3.6- and 24-$μ$m selected, FIR-bright objects in the COSMOS field, the overwhelming majority of multi-component FIR systems are line of sight projections within the 18.1 arcsec Herschel beam, rather than physical associations.
△ Less
Submitted 25 July, 2018;
originally announced July 2018.
-
Main sequence of star forming galaxies beyond the Herschel confusion limit
Authors:
W. J. Pearson,
L. Wang,
P. D. Hurley,
K. Małek,
V. Buat,
D. Burgarella,
D. Farrah,
S. J. Oliver,
D. J. B. Smith,
F. F. S. van der Tak
Abstract:
Deep far-infrared (FIR) cosmological surveys are known to be affected by source confusion, causing issues when examining the main sequence (MS) of star forming galaxies. This has typically been partially tackled by the use of stacking. However, stacking only provides the average properties of the objects in the stack. This work aims to trace the MS over $0.2\leq z<6.0$ using the latest de-blended…
▽ More
Deep far-infrared (FIR) cosmological surveys are known to be affected by source confusion, causing issues when examining the main sequence (MS) of star forming galaxies. This has typically been partially tackled by the use of stacking. However, stacking only provides the average properties of the objects in the stack. This work aims to trace the MS over $0.2\leq z<6.0$ using the latest de-blended Herschel photometry, which reaches $\approx10$ times deeper than the 5$σ$ confusion limit in SPIRE. This provides more reliable star formation rates (SFRs), especially for the fainter galaxies, and hence a more reliable MS. We built a pipeline that uses the spectral energy distribution (SED) modelling and fitting tool CIGALE to generate flux density priors in the Herschel SPIRE bands. These priors were then fed into the de-blending tool XID+ to extract flux densities from the SPIRE maps. Multi-wavelength data were combined with the extracted SPIRE flux densities to constrain SEDs and provide stellar mass (M$_{\star}$) and SFRs. These M$_{\star}$ and SFRs were then used to populate the SFR-M$_{\star}$ plane over $0.2\leq z<6.0$. No significant evidence of a high-mass turn-over was found; the best fit is thus a simple two-parameter power law of the form log(SFR)$=α[$log(M$_{\star})-10.5]+β$. The normalisation of the power law increases with redshift, rapidly at $z\lesssim1.8$, from $0.58\pm0.09$ at $z\approx0.37$ to $1.31\pm0.08$ at $z\approx1.8$. The slope is also found to increase with redshift, perhaps with an excess around $1.8\leq z<2.9$. The increasing slope indicates that galaxies become more self-similar as redshift increases, implying that the specific SFR of high-mass galaxies increases over $z=0.2$ to $z=6.0$, becoming closer to that of low-mass galaxies. The excess in the slope at $1.8\leq z<2.9$, if present, coincides with the peak of the cosmic star formation history.
△ Less
Submitted 14 May, 2018; v1 submitted 10 April, 2018;
originally announced April 2018.
-
Red, redder, reddest: SCUBA-2 imaging of colour-selected \textit{Herschel} sources
Authors:
S. Duivenvoorden,
S. Oliver,
J. M. Scudder,
J. Greenslade,
D. A. Riechers,
S. M. Wilkins,
V. Buat,
S. C. Chapman,
D. L. Clements,
A. Cooray,
K. E. K. Coppin,
H. Dannerbauer,
G. De Zotti,
J. S. Dunlop,
S. A. Eales,
A. Efstathiou,
D. Farrah,
J. E. Geach,
W. S. Holland,
P. D. Hurley,
R. J. Ivison,
L. Marchetti,
G. Petitpas,
M. T. Sargent,
D. Scott
, et al. (6 additional authors not shown)
Abstract:
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 $μ$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 d…
▽ More
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 $μ$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~$μ$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~$μ$m data provide a better constraint on the photometric redshifts of the candidates, with photometric redshift errors decreasing from $σ_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}$).
△ Less
Submitted 22 January, 2018;
originally announced January 2018.
-
HERUS: The Far-IR/Submm Spectral Energy Distributions of Local ULIRGs & Photometric Atlas
Authors:
D. L. Clements,
C. Pearson,
D. Farrah,
J. Greenslade,
Jeronimo Bernard-Salas,
E. Gonzalez-Alfonso,
J. Afonso,
A. Efstathiou,
D. Rigopoulou,
V. Lebouteiller,
P. D. Hurley,
H. Spoon
Abstract:
We present the Herschel-SPIRE photometric atlas for a complete flux limited sample of 43 local Ultraluminous Infrared Galaxies (ULIRGs), selected at 60$μ$m by IRAS, as part of the HERschel ULIRG Survey (HERUS). Photometry observations were obtained using the SPIRE instrument at 250, 350 and 500$μ$m. We describe these observations, present the results, and combine the new observations with data fro…
▽ More
We present the Herschel-SPIRE photometric atlas for a complete flux limited sample of 43 local Ultraluminous Infrared Galaxies (ULIRGs), selected at 60$μ$m by IRAS, as part of the HERschel ULIRG Survey (HERUS). Photometry observations were obtained using the SPIRE instrument at 250, 350 and 500$μ$m. We describe these observations, present the results, and combine the new observations with data from IRAS to examine the far-IR spectral energy distributions (SEDs) of these sources. We fit the observed SEDs of HERUS objects with a simple parameterised modified black body model where temperature and emissivity $β$ are free parameters. We compare the fitted values to those of non-ULIRG local galaxies, and find, in agreement with earlier results, that HERUS ULIRGs have warmer dust (median temperature T = 37.9+/-4.7 K compared to 21.3+/-3.4 K) but a similar $β$ distribution (median $β$ = 1.7 compared to 1.8) to the Herschel reference sample (HRS, Cortese et al., 2014) galaxies. Dust masses are found to be in the range of 10^7.5 to 10^9 Msun significantly higher than that of Herschel Reference Sample (HRS) sources. We compare our results for local ULIRGs with higher redshift samples selected at 250 and 850$μ$m. These latter sources generally have cooler dust and/or redder 100-to-250 $μ$m colours than our 60$μ$m-selected ULIRGs. We show that this difference may in part be the result of the sources being selected at different wavelengths rather than being a simple indication of rapid evolution in the properties of the population.
△ Less
Submitted 13 December, 2017;
originally announced December 2017.
-
De-blending Deep Herschel Surveys: A Multi-wavelength Approach
Authors:
W. J. Pearson,
L. Wang,
F. F. S. van der Tak,
P. D. Hurley,
D. Burgarella,
S. J. Oliver
Abstract:
Cosmological surveys in the far infrared are known to suffer from confusion. The Bayesian de-blending tool, XID+, currently provides one of the best ways to de-confuse deep Herschel SPIRE images, using a flat flux density prior. This work is to demonstrate that existing multi-wavelength data sets can be exploited to improve XID+ by providing an informed prior, resulting in more accurate and precis…
▽ More
Cosmological surveys in the far infrared are known to suffer from confusion. The Bayesian de-blending tool, XID+, currently provides one of the best ways to de-confuse deep Herschel SPIRE images, using a flat flux density prior. This work is to demonstrate that existing multi-wavelength data sets can be exploited to improve XID+ by providing an informed prior, resulting in more accurate and precise extracted flux densities. Photometric data for galaxies in the COSMOS field were used to constrain spectral energy distributions (SEDs) using the fitting tool CIGALE. These SEDs were used to create Gaussian prior estimates in the SPIRE bands for XID+. The multi-wavelength photometry and the extracted SPIRE flux densities were run through CIGALE again to allow us to compare the performance of the two priors. Inferred ALMA flux densities (F$^i$), at 870$μ$m and 1250$μ$m, from the best fitting SEDs from the second CIGALE run were compared with measured ALMA flux densities (F$^m$) as an independent performance validation. Similar validations were conducted with the SED modelling and fitting tool MAGPHYS and modified black body functions to test for model dependency. We demonstrate a clear improvement in agreement between the flux densities extracted with XID+ and existing data at other wavelengths when using the new informed Gaussian prior over the original uninformed prior. The residuals between F$^m$ and F$^i$ were calculated. For the Gaussian prior, these residuals, expressed as a multiple of the ALMA error ($σ$), have a smaller standard deviation, 7.95$σ$ for the Gaussian prior compared to 12.21$σ$ for the flat prior, reduced mean, 1.83$σ$ compared to 3.44$σ$, and have reduced skew to positive values, 7.97 compared to 11.50. These results were determined to not be significantly model dependent. This results in statistically more reliable SPIRE flux densities.
△ Less
Submitted 7 April, 2017;
originally announced April 2017.
-
XID+ a new prior-based extraction tool for Herschel-SPIRE maps
Authors:
Yannick Roehlly,
Peter D. Hurley
Abstract:
We present XID+ a new generation of software for prior-based photometry extraction in the Herschel SPIRE maps. Based on a Bayesian framework, XID+ allows the inclusion of prior information and gives access to the full posterior probability distribution of fluxes. XID+ is developed within the Herschel Extragalactic Legacy Project (HELP) and is available at https://github.com/H-E-L-P/XID_plus.
We present XID+ a new generation of software for prior-based photometry extraction in the Herschel SPIRE maps. Based on a Bayesian framework, XID+ allows the inclusion of prior information and gives access to the full posterior probability distribution of fluxes. XID+ is developed within the Herschel Extragalactic Legacy Project (HELP) and is available at https://github.com/H-E-L-P/XID_plus.
△ Less
Submitted 18 November, 2016;
originally announced November 2016.
-
HERUS: A CO Atlas from SPIRE Spectroscopy of local ULIRGs
Authors:
Chris Pearson,
Dimitra Rigopoulou,
Peter Hurley,
Duncan Farrah,
Jose Afonso,
Jeronimo Bernard-Salas,
Colin Borys,
David L. Clements,
Diane Cormier,
Andreas Efstathiou,
Eduardo Gonzalez-Alfonso,
Vianney Lebouteiller,
Henrik Spoon
Abstract:
We present the Herschel SPIRE Fourier Transform Spectroscopy (FTS) atlas for a complete flux limited sample of local Ultra-Luminous Infra-Red Galaxies as part of the HERschel ULIRG Survey (HERUS). The data reduction is described in detail and was optimized for faint FTS sources with particular care being taken with the subtraction of the background which dominates the continuum shape of the spectr…
▽ More
We present the Herschel SPIRE Fourier Transform Spectroscopy (FTS) atlas for a complete flux limited sample of local Ultra-Luminous Infra-Red Galaxies as part of the HERschel ULIRG Survey (HERUS). The data reduction is described in detail and was optimized for faint FTS sources with particular care being taken with the subtraction of the background which dominates the continuum shape of the spectra. Special treatment in the data reduction has been given to any observation suffering from artefacts in the data caused by anomalous instrumental effects to improve the final spectra. Complete spectra are shown covering $200 - 671μ$m with photometry in the SPIRE bands at 250$μ$m, 350$μ$m and 500$μ$m. The spectra include near complete CO ladders for over half of our sample, as well as fine structure lines from [CI] 370 $μ$m, [CI] 609 $μ$m, and [NII] 205 $μ$m. We also detect H$_{2}$O lines in several objects. We construct CO Spectral Line Energy Distributions (SLEDs) for the sample, and compare their slopes with the far-infrared colours and luminosities. We show that the CO SLEDs of ULIRGs can be broadly grouped into three classes based on their excitation. We find that the mid-J (5$<$J$<$8) lines are better correlated with the total far-infrared luminosity, suggesting that the warm gas component is closely linked to recent star-formation. The higher J transitions do not linearly correlate with the far-infrared luminosity, consistent with them originating in hotter, denser gas unconnected to the current star-formation. {\bf We conclude that in most cases more than one temperature components are required to model the CO SLEDs.}
△ Less
Submitted 19 October, 2016;
originally announced October 2016.
-
Extreme star formation events in quasar hosts over ${\bf0.5<\textit{z}<4}$
Authors:
Lura K. Pitchford,
Evanthia Hatziminaoglou,
Anna Feltre,
Duncan Farrah,
Charlotte Clarke,
Kathryn Harris,
Peter Hurley,
Sebastian Oliver,
Mathew Page,
Lingyu Wang
Abstract:
We explore the relationship between active galactic nuclei and star formation in a sample of 513 optically luminous type 1 quasars up to redshifts of $\sim$4 hosting extremely high star formation rates (SFRs). The quasars are selected to be individually detected by the \textit{Herschel} SPIRE instrument at $> $3$σ$ at 250 $μ$m, leading to typical SFRs of order of 1000 M$_{\odot}$yr$^{-1}$. We find…
▽ More
We explore the relationship between active galactic nuclei and star formation in a sample of 513 optically luminous type 1 quasars up to redshifts of $\sim$4 hosting extremely high star formation rates (SFRs). The quasars are selected to be individually detected by the \textit{Herschel} SPIRE instrument at $> $3$σ$ at 250 $μ$m, leading to typical SFRs of order of 1000 M$_{\odot}$yr$^{-1}$. We find the average SFRs to increase by almost a factor 10 from $z\sim0.5$ to $z\sim3$, mirroring the rise in the comoving SFR density over the same epoch. However, we find that the SFRs remain approximately constant with increasing accretion luminosity for accretion luminosities above 10$^{12}$ L$_{\odot}$. We also find that the SFRs do not correlate with black hole mass. Both of these results are most plausibly explained by the existence of a self-regulation process by the starburst at high SFRs, which controls SFRs on time-scales comparable to or shorter than the AGN or starburst duty cycles. We additionally find that SFRs do not depend on Eddington ratio at any redshift, consistent with no relation between SFR and black hole growth rate per unit black hole mass. Finally, we find that high-ionisation broad absorption line (HiBAL) quasars have indistinguishable far-infrared properties to those of classical quasars, consistent with HiBAL quasars being normal quasars observed along a particular line of sight, with the outflows in HiBAL quasars not having any measurable effect on the star formation in their hosts.
△ Less
Submitted 21 July, 2016;
originally announced July 2016.
-
HELP: XID+, The Probabilistic De-blender for Herschel SPIRE maps
Authors:
P. D. Hurley,
S. Oliver,
M. Betancourt,
C. Clarke,
W. I. Cowley,
S. Duivenvoorden,
D. Farrah,
M. Griffin,
C. Lacey,
E. Le Floc'h,
A. Papadopoulos,
M. Sargent,
J. M. Scudder,
M. Vaccari,
I. Valtchanov,
L. Wang
Abstract:
We have developed a new prior-based source extraction tool, XID+, to carry out photometry in the Herschel SPIRE maps at the positions of known sources. XID+ is developed using a probabilistic Bayesian framework which provides a natural framework in which to include prior information, and uses the Bayesian inference tool, Stan, to obtain the full posterior probability distribution on flux estimates…
▽ More
We have developed a new prior-based source extraction tool, XID+, to carry out photometry in the Herschel SPIRE maps at the positions of known sources. XID+ is developed using a probabilistic Bayesian framework which provides a natural framework in which to include prior information, and uses the Bayesian inference tool, Stan, to obtain the full posterior probability distribution on flux estimates. In this paper, we discuss the details of XID+ and demonstrate the basic capabilities and performance by running it on simulated SPIRE maps resembling the COSMOS field, and comparing to the current prior-based source extraction tool DESPHOT. We show that not only does XID+ perform better on metrics such as flux accuracy and flux uncertainty accuracy, we illustrate how obtaining the posterior probability distribution can help overcome some of the issues inherent with maximum likelihood based source extraction routines. We run XID+ on the COSMOS SPIRE maps from HerMES, using a 24 $\mathrm{μm}$ catalogue as a prior and show the marginalised SPIRE colour-colour plot and marginalised contribution to the cosmic infrared background at the SPIRE wavelengths. XID+ is a core tool arising from the Herschel Extragalactic Legacy Project (HELP) and we discuss how additional work within HELP providing prior information on fluxes can and will be utilised. The software is available at https://github.com/H-E-L-P/XID_plus. We also provide the data product for COSMOS. We believe this is the first time that the full posterior probability of galaxy photometry has been provided as a data product.
△ Less
Submitted 18 June, 2016;
originally announced June 2016.
-
HELP: star formation as function of galaxy environmentwith Herschel
Authors:
S. Duivenvoorden,
S. Oliver,
V. Buat,
B. Darvish,
A. Efstathiou,
D. Farrah,
M. Griffin,
P. D. Hurley,
E. Ibar,
M. Jarvis,
A. Papadopoulos,
M. T. Sargent,
D. Scott,
J. M. Scudder,
M. Symeonidis,
M. Vaccari,
M. P. Viero,
L. Wang
Abstract:
The Herschel Extragalactic Legacy Project (HELP) brings together a vast range of data from many astronomical observatories. Its main focus is on the Herschel data, which maps dust obscured star formation over 1300 deg$^2$. With this unprecedented combination of data sets, it is possible to investigate how the star formation vs stellar mass relation (main-sequence) of star-forming galaxies depends…
▽ More
The Herschel Extragalactic Legacy Project (HELP) brings together a vast range of data from many astronomical observatories. Its main focus is on the Herschel data, which maps dust obscured star formation over 1300 deg$^2$. With this unprecedented combination of data sets, it is possible to investigate how the star formation vs stellar mass relation (main-sequence) of star-forming galaxies depends on environment. In this pilot study we explore this question between 0.1 < z < 3.2 using data in the COSMOS field. We estimate the local environment from a smoothed galaxy density field using the full photometric redshift probability distribution. We estimate star formation rates by stacking the SPIRE data from the Herschel Multi-tiered Extragalactic Survey (HerMES). Our analysis rules out the hypothesis that the main-sequence for star-forming systems is independent of environment at 1.5 < z < 2, while a simple model in which the mean specific star formation rate declines with increasing environmental density gives a better description. However, we cannot exclude a simple hypothesis in which the main-sequence for star-forming systems is independent of environment at z < 1.5 and z > 2. We also estimate the evolution of the star formation rate density in the COSMOS field and our results are consistent with previous measurements at z < 1.5 and z > 2 but we find a $1.4^{+0.3}_{-0.2}$ times higher peak value of the star formation rate density at $z \sim 1.9$.
△ Less
Submitted 16 June, 2016;
originally announced June 2016.
-
The multiplicity of 250-$μ$m Herschel sources in the COSMOS field
Authors:
Jillian M. Scudder,
Seb Oliver,
Peter D. Hurley,
Matt Griffin,
Mark T. Sargent,
Douglas Scott,
Lingyu Wang,
Julie L. Wardlow
Abstract:
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 prob…
▽ More
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.
△ Less
Submitted 29 April, 2016;
originally announced May 2016.
-
Star Formation in Luminous Quasars at 2<z<3
Authors:
Kathryn Harris,
Duncan Farrah,
Bernhard Schulz,
Evanthia Hatziminaoglou,
Marco Viero,
Nick Anderson,
Matthieu Bethermin,
Scott Chapman,
David L. Clements,
Asantha Cooray,
Andreas Efstathiou,
Anne Feltre,
Peter Hurley,
Eduardo Ibar,
Mark Lacy,
Sebastian Oliver,
Mathew J. Page,
Ismael Perez-Fournon,
Sara M. Petty,
Lura K. Pitchford,
Dimitra Rigopoulou,
Douglas Scott,
Myrto Symeonidis,
Joaquin Vieira,
Lingyu Wang
Abstract:
We investigate the relation between star formation rates ($\dot{M}_{s}$) and AGN properties in optically selected type 1 quasars at $2<z<3$ using data from Herschel and the SDSS. We find that $\dot{\rm{M}}_s$ remains approximately constant with redshift, at $300\pm100~\rm{M}_{\odot}$yr$^{-1}$. Conversely, $\dot{\rm{M}}_s$ increases with AGN luminosity, up to a maximum of $\sim600~\rm{M}_{\odot}$yr…
▽ More
We investigate the relation between star formation rates ($\dot{M}_{s}$) and AGN properties in optically selected type 1 quasars at $2<z<3$ using data from Herschel and the SDSS. We find that $\dot{\rm{M}}_s$ remains approximately constant with redshift, at $300\pm100~\rm{M}_{\odot}$yr$^{-1}$. Conversely, $\dot{\rm{M}}_s$ increases with AGN luminosity, up to a maximum of $\sim600~\rm{M}_{\odot}$yr$^{-1}$, and with CIV FWHM. In context with previous results, this is consistent with a relation between $\dot{\rm{M}}_s$ and black hole accretion rate ($\dot{\rm{M}}_{bh}$) existing in only parts of the $z-\dot{\rm{M}}_{s}-\dot{\rm{M}}_{bh}$ plane, dependent on the free gas fraction, the trigger for activity, and the processes that may quench star formation. The relations between $\dot{\rm{M}}_s$ and both AGN luminosity and CIV FWHM are consistent with star formation rates in quasars scaling with black hole mass, though we cannot rule out a separate relation with black hole accretion rate. Star formation rates are observed to decline with increasing CIV equivalent width. This decline can be partially explained via the Baldwin effect, but may have an additional contribution from one or more of three factors; $M_i$ is not a linear tracer of L$_{2500}$, the Baldwin effect changes form at high AGN luminosities, and high CIV EW values signpost a change in the relation between $\dot{\rm{M}}_s$ and $\dot{\rm{M}}_{bh}$. Finally, there is no strong relation between $\dot{\rm{M}}_s$ and Eddington ratio, or the asymmetry of the CIV line. The former suggests that star formation rates do not scale with how efficiently the black hole is accreting, while the latter is consistent with CIV asymmetries arising from orientation effects.
△ Less
Submitted 8 February, 2016;
originally announced February 2016.
-
Dust properties of Lyman break galaxies at $z\sim3$
Authors:
J. Álvarez-Márquez,
D. Burgarella,
S. Heinis,
V. Buat,
B. Lo Faro,
M. Béthermin,
C. E. López-Fortín,
A. Cooray,
D. Farrah,
P. Hurley,
E. Ibar,
O. Ilbert,
A. M. Koekemoer,
B. C. Lemaux,
I. Pérez-Fournon,
G. Rodighiero,
M. Salvato,
D. Scott,
Y. Taniguchi,
J. D. Vieira,
L. Wang
Abstract:
We explore from a statistical point of view the far-infrared (far-IR) and sub-millimeter (sub-mm) properties of a large sample of LBGs (22,000) at z~3 in the COSMOS field. The large number of galaxies allows us to split it in several bins as a function of UV luminosity, UV slope, and stellar mass to better sample their variety. We perform stacking analysis in PACS (100 and 160 um), SPIRE (250, 350…
▽ More
We explore from a statistical point of view the far-infrared (far-IR) and sub-millimeter (sub-mm) properties of a large sample of LBGs (22,000) at z~3 in the COSMOS field. The large number of galaxies allows us to split it in several bins as a function of UV luminosity, UV slope, and stellar mass to better sample their variety. We perform stacking analysis in PACS (100 and 160 um), SPIRE (250, 350 and 500 um) and AzTEC (1.1 mm) images. Our stacking procedure corrects the biases induced by galaxy clustering and incompleteness of our input catalogue in dense regions. We obtain the full IR spectral energy distributions (SED) of subsamples of LBGs and derive the mean IR luminosity as a function of UV luminosity, UV slope, and stellar mass. The average IRX is roughly constant over the UV luminosity range, with a mean of 7.9 (1.8 mag). However, it is correlated with UV slope, and stellar mass. We investigate using a statistically-controlled stacking analysis as a function of (stellar mass, UV slope) the dispersion of the IRX-UVslope and IRX-M* plane. Our results enable us to study the average relation between star-formation rate (SFR) and stellar mass, and we show that our LBG sample lies on the main sequence of star formation at z~3.
△ Less
Submitted 13 December, 2015;
originally announced December 2015.
-
Collaborative Randomized Beamforming for Phased Array Radio Interferometers
Authors:
Orhan Ocal,
Paul Hurley,
Giovanni Cherubini,
Sanaz Kazemi
Abstract:
The Square Kilometre Array (SKA) will form the largest radio telescope ever built and such a huge instrument in the desert poses enormous engineering and logistic challenges. Algorithmic and architectural breakthroughs are needed.
Data is collected and processed in groups of antennas before transport for central processing. This processing includes beamforming, primarily so as to reduce the amou…
▽ More
The Square Kilometre Array (SKA) will form the largest radio telescope ever built and such a huge instrument in the desert poses enormous engineering and logistic challenges. Algorithmic and architectural breakthroughs are needed.
Data is collected and processed in groups of antennas before transport for central processing. This processing includes beamforming, primarily so as to reduce the amount of data sent. The principal existing technique points to a region of interest independently of the sky model and how the other stations beamform.
We propose a new collaborative beamforming algorithm in order to maximize information captured at the stations (thus reducing the amount of data transported). The method increases the diversity in measurements through randomized beam- forming. We demonstrate through numerical simulation the effectiveness of the method. In particular, we show that randomized beamforming can achieve the same image quality while producing 40% less data when compared to the prevailing method matched beamforming.
△ Less
Submitted 14 November, 2014;
originally announced November 2014.
-
Learning the Fundamental MIR Spectral Components of Galaxies with Non-Negative Matrix Factorisation
Authors:
P. D. Hurley,
S. Oliver,
D. Farrah,
V. Lebouteiller,
H. W. W. Spoon
Abstract:
The mid-infrared (MIR) spectra observed with the \textit{Spitzer} Infrared Spectrograph (IRS) provide a valuable dataset for untangling the physical processes and conditions within galaxies.
This paper presents the first attempt to blindly learn fundamental spectral components of MIR galaxy spectra, using non-negative matrix factorisation (NMF). NMF is a recently developed multivariate technique…
▽ More
The mid-infrared (MIR) spectra observed with the \textit{Spitzer} Infrared Spectrograph (IRS) provide a valuable dataset for untangling the physical processes and conditions within galaxies.
This paper presents the first attempt to blindly learn fundamental spectral components of MIR galaxy spectra, using non-negative matrix factorisation (NMF). NMF is a recently developed multivariate technique shown to be successful in blind source separation problems. Unlike the more popular multivariate analysis technique, principal component analysis, NMF imposes the condition that weights and spectral components are non-negative. This more closely resembles the physical process of emission in the mid-infrared, resulting in physically intuitive components. By applying NMF to galaxy spectra in the Cornell Atlas of Spitzer/IRS sources (CASSIS), we find similar components amongst different NMF sets. These similar components include two for AGN emission and one for star formation. [... ABBREVIATED...] We show an NMF set with seven components can reconstruct the general spectral shape of a wide variety of objects, though struggle to fit the varying strength of emission lines. We also show that the seven components can be used to separate out different types of objects. We model this separation with Gaussian Mixtures modelling and use the result to provide a classification tool.
We also show the NMF components can be used to separate out the emission from AGN and star formation regions and define a new star formation/AGN diagnostic which is consistent with all mid-infrared diagnostics already in use but has the advantage that it can be applied to mid-infrared spectra with low signal to noise or with limited spectral range. The 7 NMF components and code for classification are made public on arxiv and are available at: \url{https://github.com/pdh21/NMF_software/}
△ Less
Submitted 2 October, 2013;
originally announced October 2013.
-
Herschel observations and a model for IRAS 08572+3915: a candidate for the most luminous infrared galaxy in the local (z < 0.2) Universe
Authors:
A. Efstathiou,
C. Pearson,
D. Farrah,
D. Rigopoulou,
J. Gracia-Carpio,
A. Verma,
H. W. W. Spoon,
J. Afonso,
J. Bernard-Salas,
D. L. Clements,
A. Cooray,
D. Cormier,
M. Etxaluze,
J. Fischer,
E. Gonzalez-Alfonso,
P. Hurley,
V. Lebouteiller,
S. J. Oliver,
M. Rowan-Robinson,
E. Sturm
Abstract:
We present Herschel photometry and spectroscopy, carried out as part of the Herschel ULIRG survey (HERUS), and a model for the infrared to submillimetre emission of the ultraluminous infrared galaxy IRAS 08572+3915. This source shows one of the deepest known silicate absorption features and no polycyclic aromatic hydrocarbon (PAH) emission. The model suggests that this object is powered by an acti…
▽ More
We present Herschel photometry and spectroscopy, carried out as part of the Herschel ULIRG survey (HERUS), and a model for the infrared to submillimetre emission of the ultraluminous infrared galaxy IRAS 08572+3915. This source shows one of the deepest known silicate absorption features and no polycyclic aromatic hydrocarbon (PAH) emission. The model suggests that this object is powered by an active galactic nucleus (AGN) with a fairly smooth torus viewed almost edge-on and a very young starburst. According to our model the AGN contributes about 90% of the total luminosity of 1.1 x 10^13 Lo, which is about a factor of five higher than previous estimates. The large correction of the luminosity is due to the anisotropy of the emission of the best fit torus. Similar corrections may be necessary for other local and high-z analogs. This correction implies that IRAS 08572+3915 at a redshift of 0.05835 may be the nearest hyperluminous infrared galaxy and probably the most luminous infrared galaxy in the local (z < 0.2) Universe. IRAS 08572+3915 shows a low ratio of [CII] to IR luminosity (log L_[CII]/L_{IR} < -3.8) and a [OI]63um to [CII]158um line ratio of about 1 that supports the model presented in this letter.
△ Less
Submitted 26 September, 2013;
originally announced September 2013.
-
Mid- to far infrared properties of star-forming galaxies and active galactic nuclei
Authors:
G. E. Magdis,
D. Rigopoulou,
G. Helou,
D. Farrah,
P. Hurley,
A. Alonso-Herrero,
J. Bock,
D. Burgarella,
S. Chapman,
V. Charmandaris,
A. Cooray,
Y. S. Dai,
D. Dale,
D. Elbaz,
A. Feltre,
E. Hatziminaoglou,
J-S Huang,
G. Morrison,
S. Oliver,
M. Page,
D. Scott,
Y. Shi
Abstract:
We study the mid- to far-IR properties of a 24um-selected flux-limited sample (S24 > 5mJy) of 154 intermediate redshift (<z>~0.15), infrared luminous galaxies, drawn from the 5MUSES survey. By combining existing mid-IR spectroscopy and new Herschel SPIRE submm photometry from the HerMES program, we derived robust total infrared luminosity (LIR) and dust mass (Md) estimates and infered the relative…
▽ More
We study the mid- to far-IR properties of a 24um-selected flux-limited sample (S24 > 5mJy) of 154 intermediate redshift (<z>~0.15), infrared luminous galaxies, drawn from the 5MUSES survey. By combining existing mid-IR spectroscopy and new Herschel SPIRE submm photometry from the HerMES program, we derived robust total infrared luminosity (LIR) and dust mass (Md) estimates and infered the relative contribution of the AGN to the infrared energy budget of the sources. We found that the total infrared emission of galaxies with weak 6.2um PAH emission (EW<0.2um) is dominated by AGN activity, while for galaxies with EW>0.2um more than 50% of the LIR arises from star formation. We also found that for galaxies detected in the 250-500um Herschel bands an AGN has a statistically insignificant effect on the temperature of the cold dust and the far-IR colours of the host galaxy, which are primarily shaped by star formation activity. For star-forming galaxies we reveal an anti-correlation between the LIR-to-rest-frame 8um luminosity ratio, IR8 = LIR\L8, and the strength of PAH features. We found that this anti-correlation is primarily driven by variations in the PAHs emission, and not by variations in the 5-15um mid-IR continuum emission. Using the [NeIII]/[NeII] line ratio as a tracer of the hardness of the radiation field, we confirm that galaxies with harder radiation fields tend to exhibit weaker PAH features, and found that they have higher IR8 values and higher dust-mass-weighted luminosities (LIR/Md), the latter being a proxy for the dust temperature (Td). We argue that these trends originate either from variations in the environment of the star-forming regions or are caused by variations in the age of the starburst. Finally, we provide scaling relations that will allow estimating LIR, based on single-band observations with the mid-infrared instrument, on board the upcoming JWST.
△ Less
Submitted 16 September, 2013;
originally announced September 2013.
-
Far-infrared Fine-Structure Line Diagnostics of Ultraluminous Infrared Galaxies
Authors:
Duncan Farrah,
Vianney Lebouteiller,
Henrik Spoon,
Jeronimo Bernard-Salas,
Chris Pearson,
Dimitra Rigopoulou,
Howard Smith,
Eduardo Gonzalez-Alfonso,
David Clements,
Andreas Efstathiou,
Diane Cormier,
Jose Afonso,
Sara Petty,
Kathryn Harris,
Peter Hurley,
Colin Borys,
Aprajita Verma,
Asantha Cooray,
Valentina Salvatelli
Abstract:
We present Herschel observations of six fine-structure lines in 25 Ultraluminous Infrared Galaxies at z<0.27. The lines, [O III]52, [N III]57, [O I]63, [N II]122, [O I]145, and [C II]158, are mostly single gaussians with widths <600 km s-1 and luminosities of 10^7 - 10^9 Solar. There are deficits in the [O I]63/L_IR, [N II]/L_IR, [O I]145/L_IR, and [C II]/L_IR ratios compared to lower luminosity s…
▽ More
We present Herschel observations of six fine-structure lines in 25 Ultraluminous Infrared Galaxies at z<0.27. The lines, [O III]52, [N III]57, [O I]63, [N II]122, [O I]145, and [C II]158, are mostly single gaussians with widths <600 km s-1 and luminosities of 10^7 - 10^9 Solar. There are deficits in the [O I]63/L_IR, [N II]/L_IR, [O I]145/L_IR, and [C II]/L_IR ratios compared to lower luminosity systems. The majority of the line deficits are consistent with dustier H II regions, but part of the [C II] deficit may arise from an additional mechanism, plausibly charged dust grains. This is consistent with some of the [C II] originating from PDRs or the ISM. We derive relations between far-IR line luminosities and both IR luminosity and star formation rate. We find that [N II] and both [O I] lines are good tracers of IR luminosity and star formation rate. In contrast, [C II] is a poor tracer of IR luminosity and star formation rate, and does not improve as a tracer of either quantity if the [C II] deficit is accounted for. The continuum luminosity densities also correlate with IR luminosity and star formation rate. We derive ranges for the gas density and ultraviolet radiation intensity of 10^1 < n < 10^2.5 and 10^2.2 < G_0 < 10^3.6, respectively. These ranges depend on optical type, the importance of star formation, and merger stage. We do not find relationships between far-IR line properties and several other parameters; AGN activity, merger stage, mid-IR excitation, and SMBH mass. We conclude that these far-IR lines arise from gas heated by starlight, and that they are not strongly influenced by AGN activity.
△ Less
Submitted 19 August, 2013;
originally announced August 2013.
-
Herschel-SPIRE-Fourier Transform Spectroscopy of the nearby spiral galaxy IC342
Authors:
D. Rigopoulou,
P. D. Hurley,
B. M. Swinyard,
J. Virdee,
K. V. Croxall,
R. H. B. Hopwood,
T. Lim,
G. E. Magdis,
C. P. Pearson,
E. Pellegrini,
E. Polehampton,
J-D. Smith
Abstract:
We present observations of the nearby spiral galaxy IC342 with the Herschel Spectral and Photometric Imaging Receiver (SPIRE) Fourier Transform Spectrometer. The spectral range afforded by SPIRE, 196-671 microns, allows us to access a number of 12CO lines from J=4--3 to J=13--12 with the highest J transitions observed for the first time. In addition we present measurements of 13CO, [CI] and [NII].…
▽ More
We present observations of the nearby spiral galaxy IC342 with the Herschel Spectral and Photometric Imaging Receiver (SPIRE) Fourier Transform Spectrometer. The spectral range afforded by SPIRE, 196-671 microns, allows us to access a number of 12CO lines from J=4--3 to J=13--12 with the highest J transitions observed for the first time. In addition we present measurements of 13CO, [CI] and [NII]. We use a radiative transfer code coupled with Bayesian likelihood analysis to model and constrain the temperature, density and column density of the gas. We find two 12CO components, one at 35 K and one at 400 K with CO column densities of 6.3x10^{17} cm^{-2} and 0.4x10^{17} cm^{-2} and CO gas masses of 1.26x10^{7} Msolar and 0.15x10^{7} Msolar, for the cold and warm components, respectively. The inclusion of the high-J 12CO line observations, indicate the existence of a much warmer gas component (~400 K) confirming earlier findings from H_{2} rotational line analysis from ISO and Spitzer. The mass of the warm gas is 10% of the cold gas, but it likely dominates the CO luminosity. In addition, we detect strong emission from [NII] 205microns and the {3}P_{1}->{3}P_{0} and {3}P_{2} ->{3}P_{1} [CI] lines at 370 and 608 microns, respectively. The measured 12CO line ratios can be explained by Photon-dominated region (PDR) models although additional heating by e.g. cosmic rays cannot be excluded. The measured [CI] line ratio together with the derived [C] column density of 2.1x10^{17} cm^{-2} and the fact that [CI] is weaker than CO emission in IC342 suggests that [CI] likely arises in a thin layer on the outside of the CO emitting molecular clouds consistent with PDRs playing an important role.
△ Less
Submitted 23 June, 2013;
originally announced June 2013.
-
A Dust-Obscured Massive Maximum-Starburst Galaxy at a Redshift of 6.34
Authors:
Dominik A. Riechers,
C. M. Bradford,
D. L. Clements,
C. D. Dowell,
I. Perez-Fournon,
R. J. Ivison,
C. Bridge,
A. Conley,
Hai Fu,
J. D. Vieira,
J. Wardlow,
J. Calanog,
A. Cooray,
P. Hurley,
R. Neri,
J. Kamenetzky,
J. E. Aguirre,
B. Altieri,
V. Arumugam,
D. J. Benford,
M. Bethermin,
J. Bock,
D. Burgarella,
A. Cabrera-Lavers,
S. C. Chapman
, et al. (39 additional authors not shown)
Abstract:
Massive present-day early-type (elliptical and lenticular) galaxies probably gained the bulk of their stellar mass and heavy elements through intense, dust-enshrouded starbursts - that is, increased rates of star formation - in the most massive dark matter halos at early epochs. However, it remains unknown how soon after the Big Bang such massive starburst progenitors exist. The measured redshift…
▽ More
Massive present-day early-type (elliptical and lenticular) galaxies probably gained the bulk of their stellar mass and heavy elements through intense, dust-enshrouded starbursts - that is, increased rates of star formation - in the most massive dark matter halos at early epochs. However, it remains unknown how soon after the Big Bang such massive starburst progenitors exist. The measured redshift distribution of dusty, massive starbursts has long been suspected to be biased low in redshift owing to selection effects, as confirmed by recent findings of systems out to redshift z~5. Here we report the identification of a massive starburst galaxy at redshift 6.34 through a submillimeter color-selection technique. We unambiguously determined the redshift from a suite of molecular and atomic fine structure cooling lines. These measurements reveal a hundred billion solar masses of highly excited, chemically evolved interstellar medium in this galaxy, which constitutes at least 40% of the baryonic mass. A "maximum starburst" converts the gas into stars at a rate more than 2,000 times that of the Milky Way, a rate among the highest observed at any epoch. Despite the overall downturn of cosmic star formation towards the highest redshifts, it seems that environments mature enough to form the most massive, intense starbursts existed at least as early as 880 million years after the Big Bang.
△ Less
Submitted 15 April, 2013;
originally announced April 2013.
-
The complex physics of dusty star-forming galaxies at high redshifts as revealed by Herschel and Spitzer
Authors:
B. Lo Faro,
A. Franceschini,
M. Vaccari,
L. Silva,
G. Rodighiero,
S. Berta,
J. Bock,
D. Burgarella,
V. Buat,
A. Cava,
D. L. Clements,
A. Cooray,
D. Farrah,
A. Feltre,
E. A. González Solares,
P. Hurley,
D. Lutz,
G. Magdis,
B. Magnelli,
L. Marchetti,
S. J. Oliver,
M. J. Page,
P. Popesso,
F. Pozzi,
D. Rigopoulou
, et al. (7 additional authors not shown)
Abstract:
We combine far-infrared photometry from Herschel (PEP/HERMES) with deep mid-infrared spectroscopy from Spitzer to investigate the nature and the mass assembly history of a sample of 31 Luminous and Ultraluminous Infrared Galaxies at z~1 and 2 selected in GOODS-S with 24 $μ$m fluxes between 0.2 and 0.5 mJy. We model the data with a self-consistent physical model (GRASIL) which includes a state-of-t…
▽ More
We combine far-infrared photometry from Herschel (PEP/HERMES) with deep mid-infrared spectroscopy from Spitzer to investigate the nature and the mass assembly history of a sample of 31 Luminous and Ultraluminous Infrared Galaxies at z~1 and 2 selected in GOODS-S with 24 $μ$m fluxes between 0.2 and 0.5 mJy. We model the data with a self-consistent physical model (GRASIL) which includes a state-of-the-art treatment of dust extinction and reprocessing. We find that all of our galaxies appear to require massive populations of old (>1 Gyr) stars and, at the same time, to host a moderate ongoing activity of SF (SFR < 100 M$_{\odot}$/yr). The bulk of the stars appear to have been formed a few Gyr before the observation in essentially all cases. Only five galaxies of the sample require a recent starburst superimposed on a quiescent star formation history (SFH). We also find discrepancies between our results and those based on optical-only SED fitting for the same objects; by fitting their observed Spectral Energy Distributions with our physical model we find higher extinctions (by $Δ$A_{V} ~ 0.81 and 1.14) and higher stellar masses (by $Δ$Log(M*) ~ 0.16 and 0.36 dex) for z~1 and z~2 (U)LIRGs, respectively. The stellar mass difference is larger for the most dust obscured objects. We also find lower SFRs than those computed from L_{IR} using the Kennicutt relation due to the significant contribution to the dust heating by intermediate-age stellar populations through 'cirrus' emission (~73% and ~66% of total L_{IR} for z~1 and z~2 (U)LIRGs, respectively).
△ Less
Submitted 26 November, 2012;
originally announced November 2012.
-
Galaxy and Mass Assembly (GAMA): Colour and luminosity dependent clustering from calibrated photometric redshifts
Authors:
L. Christodoulou,
C. Eminian,
J. Loveday,
P. Norberg,
I. K. Baldry,
P. D. Hurley,
S. P. Driver,
S. P. Bamford,
A. M. Hopkins,
J. Liske,
J. A. Peacock,
J. Bland-Hawthorn,
S. Brough,
E. Cameron,
C. J. Conselice,
S. M. Croom,
C. S. Frenk,
M. Gunawardhana,
D. H. Jones,
L. S. Kelvin,
K. Kuijken,
R. C. Nichol,
H. Parkinson,
K. A. Pimbblet,
C. C. Popescu
, et al. (9 additional authors not shown)
Abstract:
We measure the two-point angular correlation function of a sample of 4,289,223 galaxies with r < 19.4 mag from the Sloan Digital Sky Survey as a function of photometric redshift, absolute magnitude and colour down to M_r - 5log h = -14 mag. Photometric redshifts are estimated from ugriz model magnitudes and two Petrosian radii using the artificial neural network package ANNz, taking advantage of t…
▽ More
We measure the two-point angular correlation function of a sample of 4,289,223 galaxies with r < 19.4 mag from the Sloan Digital Sky Survey as a function of photometric redshift, absolute magnitude and colour down to M_r - 5log h = -14 mag. Photometric redshifts are estimated from ugriz model magnitudes and two Petrosian radii using the artificial neural network package ANNz, taking advantage of the Galaxy and Mass Assembly (GAMA) spectroscopic sample as our training set. The photometric redshifts are then used to determine absolute magnitudes and colours. For all our samples, we estimate the underlying redshift and absolute magnitude distributions using Monte-Carlo resampling. These redshift distributions are used in Limber's equation to obtain spatial correlation function parameters from power law fits to the angular correlation function. We confirm an increase in clustering strength for sub-L* red galaxies compared with ~L* red galaxies at small scales in all redshift bins, whereas for the blue population the correlation length is almost independent of luminosity for ~L* galaxies and fainter. A linear relation between relative bias and log luminosity is found to hold down to luminosities L~0.03L*. We find that the redshift dependence of the bias of the L* population can be described by the passive evolution model of Tegmark & Peebles (1998). A visual inspection of a random sample of our r < 19.4 sample of SDSS galaxies reveals that about 10 per cent are spurious, with a higher contamination rate towards very faint absolute magnitudes due to over-deblended nearby galaxies. We correct for this contamination in our clustering analysis.
△ Less
Submitted 5 June, 2012;
originally announced June 2012.
-
Principal Component Analysis and Radiative Transfer modelling of Spitzer IRS Spectra of Ultra Luminous Infrared Galaxies
Authors:
Peter D Hurley,
Seb Oliver,
Duncan Farrah,
Lingyu Wang,
Andreas Efstathiou
Abstract:
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…
▽ More
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/
△ Less
Submitted 29 May, 2012;
originally announced May 2012.
-
The Herschel Multi-tiered Extragalactic Survey: HerMES
Authors:
HerMES Collaboration,
S. J. Oliver,
J. Bock,
B. Altieri,
A. Amblard,
V. Arumugam,
H. Aussel,
T. Babbedge,
A. Beelen,
M. Béthermin,
A. Blain,
A. Boselli,
C. Bridge,
D. Brisbin,
V. Buat,
D. Burgarella,
N. Castro-Rodríguez,
A. Cava,
P. Chanial,
M. Cirasuolo,
D. L. Clements,
A. Conley,
L. Conversi,
A. Cooray,
C. D. Dowell
, et al. (80 additional authors not shown)
Abstract:
The Herschel Multi-tiered Extragalactic Survey, HerMES, is a legacy program designed to map a set of nested fields totalling ~380 deg^2. Fields range in size from 0.01 to ~20 deg^2, using Herschel-SPIRE (at 250, 350 and 500 μm), and Herschel-PACS (at 100 and 160 μm), with an additional wider component of 270 deg^2 with SPIRE alone. These bands cover the peak of the redshifted thermal spectral ener…
▽ More
The Herschel Multi-tiered Extragalactic Survey, HerMES, is a legacy program designed to map a set of nested fields totalling ~380 deg^2. Fields range in size from 0.01 to ~20 deg^2, using Herschel-SPIRE (at 250, 350 and 500 μm), and Herschel-PACS (at 100 and 160 μm), with an additional wider component of 270 deg^2 with SPIRE alone. These bands cover the peak of the redshifted thermal spectral energy distribution from interstellar dust and thus capture the re-processed optical and ultra-violet radiation from star formation that has been absorbed by dust, and are critical for forming a complete multi-wavelength understanding of galaxy formation and evolution.
The survey will detect of order 100,000 galaxies at 5σin some of the best studied fields in the sky. Additionally, HerMES is closely coordinated with the PACS Evolutionary Probe survey. Making maximum use of the full spectrum of ancillary data, from radio to X-ray wavelengths, it is designed to: facilitate redshift determination; rapidly identify unusual objects; and understand the relationships between thermal emission from dust and other processes. Scientific questions HerMES will be used to answer include: the total infrared emission of galaxies; the evolution of the luminosity function; the clustering properties of dusty galaxies; and the properties of populations of galaxies which lie below the confusion limit through lensing and statistical techniques.
This paper defines the survey observations and data products, outlines the primary scientific goals of the HerMES team, and reviews some of the early results.
△ Less
Submitted 12 March, 2012;
originally announced March 2012.