Source author record

M. Huynh

M. Huynh 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

40works
4topics
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

40 published item(s)

preprint2021arXiv

Remnant Radio Galaxies Discovered in a Multi-frequency Survey

The remnant phase of a radio galaxy begins when the jets launched from an active galactic nucleus are switched off. To study the fraction of radio galaxies in a remnant phase, we take advantage of a $8.31$\,deg$^2$ sub-region of the GAMA~23~field which comprises of surveys covering the frequency range 0.1--9\,GHz. We present a sample of 104 radio galaxies compiled from observations conducted by the Murchison Wide-field Array (216\,MHz), the Australia Square Kilometer Array Pathfinder (887\,MHz), and the Australia Telescope Compact Array (5.5\,GHz). We adopt an `absent radio core' criterion to identify 10 radio galaxies showing no evidence for an active nucleus. We classify these as new candidate remnant radio galaxies. Seven of these objects still display compact emitting regions within the lobes at 5.5\,GHz; at this frequency the emission is short-lived, implying a recent jet switch-off. On the other hand, only three show evidence of aged lobe plasma by the presence of an ultra-steep spectrum ($α<-1.2$) and a diffuse, low surface-brightness radio morphology. The predominant fraction of young remnants is consistent with a rapid fading during the remnant phase. Within our sample of radio galaxies, our observations constrain the remnant fraction to $4\%\lesssim f_{\mathrm{rem}} \lesssim 10\%$; the lower limit comes from the limiting case in which all remnant candidates with hotspots are simply active radio galaxies with faint, undetected radio cores. Finally, we model the synchrotron spectrum arising from a hotspot to show they can persist for 5--10\,Myr at 5.5\,GHz after the jets switch off -- radio emission arising from such hotspots can therefore be expected in an appreciable fraction of genuine remnants.

preprint2020arXiv

Cataloging the radio-sky with unsupervised machine learning: a new approach for the SKA era

We develop a new analysis approach towards identifying related radio components and their corresponding infrared host galaxy based on unsupervised machine learning methods. By exploiting PINK, a self-organising map algorithm, we are able to associate radio and infrared sources without the a priori requirement of training labels. We present an example of this method using $894,415$ images from the FIRST and WISE surveys centred towards positions described by the FIRST catalogue. We produce a set of catalogues that complement FIRST and describe 802,646 objects, including their radio components and their corresponding AllWISE infrared host galaxy. Using these data products we (i) demonstrate the ability to identify objects with rare and unique radio morphologies (e.g. 'X'-shaped galaxies, hybrid FR-I/FR-II morphologies), (ii) can identify the potentially resolved radio components that are associated with a single infrared host and (iii) introduce a "curliness" statistic to search for bent and disturbed radio morphologies, and (iv) extract a set of 17 giant radio galaxies between 700-1100 kpc. As we require no training labels, our method can be applied to any radio-continuum survey, provided a sufficiently representative SOM can be trained.

preprint2020arXiv

FLASH Early Science -- Discovery of an intervening HI 21-cm absorber from an ASKAP survey of the GAMA 23 field

We present early science results from the First Large Absorption Survey in HI (FLASH), a spectroscopically blind survey for 21-cm absorption lines in cold hydrogen HI gas at cosmological distances using the Australian Square Kilometre Array Pathfinder (ASKAP). We have searched for HI absorption towards 1253 radio sources in the GAMA 23 field, covering redshifts between $z = 0.34$ and $0.79$ over a sky area of approximately 50 deg$^{2}$. In a purely blind search we did not obtain any detections of 21-cm absorbers above our reliability threshold. Assuming a fiducial value for the HI spin temperature of $T_{\rm spin}$ = 100 K and source covering fraction $c_{\rm f} = 1$, the total comoving absorption path length sensitive to all Damped Lyman $α$ Absorbers (DLAs; $N_{\rm HI} \geq 2 \times 10^{20}$ cm$^{-2}$) is $Δ{X} = 6.6 \pm 0.3$ ($Δ{z} = 3.7 \pm 0.2$) and super-DLAs ($N_{\rm HI} \geq 2 \times 10^{21}$ cm$^{-2}$) is $Δ{X} = 111 \pm 6$ ($Δ{z} = 63 \pm 3$). We estimate upper limits on the HI column density frequency distribution function that are consistent with measurements from prior surveys for redshifted optical DLAs, and nearby 21-cm emission and absorption. By cross matching our sample of radio sources with optical spectroscopic identifications of galaxies in the GAMA 23 field, we were able to detect 21-cm absorption at $z = 0.3562$ towards NVSS J224500$-$343030, with a column density of $N_{\rm HI} = (1.2 \pm 0.1) \times 10^{20} (T_{\rm spin}/100~\mathrm{K})$ cm$^{-2}$. The absorber is associated with GAMA J22450.05$-$343031.7, a massive early-type galaxy at an impact parameter of 17 kpc with respect to the radio source and which may contain a massive ($M_{\rm HI} \gtrsim 3 \times 10^{9}$ M$_{\odot}$) gas disc. Such gas-rich early types are rare, but have been detected in the nearby Universe.

preprint2020arXiv

PKS 2250$-$351: A Giant Radio Galaxy in Abell 3936

We present a detailed analysis of the radio galaxy PKS 2250-351, a giant of 1.2 Mpc projected size, its host galaxy, and its environment. We use radio data from the Murchison Widefield Array, the upgraded Giant Metre-wavelength Radio Telescope, the Australian Square Kilometre Array Pathfinder, and the Australia Telescope Compact Array to model the jet power and age. Optical and infra-red data come from the Galaxy And Mass Assembly (GAMA) survey and provide information on the host galaxy and environment. GAMA spectroscopy confirms that PKS 2250-351 lies at z=0.2115 in the irregular, and likely unrelaxed, cluster Abell 3936. We find its host is a massive, `red and dead' elliptical galaxy with negligible star formation but with a highly obscured active galactic nucleus dominating the mid-infrared emission. Assuming it lies on the local M-sigma relation it has an Eddington accretion rate of lambda_EDD~0.014. We find that the lobe-derived jet power (a time-averaged measure) is an order of magnitude greater than the hotspot-derived jet power (an instantaneous measure). We propose that over the lifetime of the observed radio emission (~300 Myr) the accretion has switched from an inefficient advection dominated mode to a thin-disc efficient mode, consistent with the decrease in jet power. We also suggest that the asymmetric radio morphology is due to its environment, with the host of PKS 2250-351 lying to the west of the densest concentration of galaxies in Abell 3936.

preprint2020arXiv

The XXL Survey. XLI. Radio AGN luminosity functions based on the GMRT $610 \ \mathrm{MHz}$ continuum observations

We study the space density evolution of active galactic nuclei (AGN) using the $610 \ \mathrm{MHz}$ radio survey of the XXL-North field, performed with the Giant Metrewave Radio Telescope (GMRT). The survey covers an area of $30.4 \ \mathrm{deg}^2$, with a beamsize of $6.5 \ \mathrm{arcsec}$. The survey is divided into two parts, one covering an area of $11.9 \ \mathrm{deg}^2$ with $1 σ$ rms noise of $200\ \mathrm{μJy \ beam^{-1}}$ and the other spanning $18.5 \ \mathrm{deg}^2$ with rms noise of $45\ \mathrm{μJy \ beam^{-1}}$. We extracted the catalog of radio components above $7 σ$. The catalog was cross-matched with a multi-wavelength catalog of the XXL-North field (covering about $80 \%$ of the radio XXL-North field) using a likelihood ratio method, which determines the counterparts based on their positions and their optical properties. The multi-component sources were matched visually with the aid of a computer code: Multi-Catalog Visual Cross-Matching (MCVCM). A flux density cut above $1\ \mathrm{mJy}$ selects AGN hosts with a high purity in terms of star formation contamination based on the available source counts. After cross-matching and elimination of observational biases arising from survey incompletenesses, the number of remaining sources was $1150$. We constructed the rest-frame $1.4 \ \mathrm{GHz}$ radio luminosity functions of these sources using the maximum volume method. This survey allows us to probe luminosities of $ 23 \lesssim \log(L_{1.4 \ \mathrm{GHz}}[\mathrm{W/Hz}]) \lesssim 28$ up to redshifts of $z \approx 2.1$. Our results are consistent with the results from the literature in which AGN are comprised of two differently evolving populations, where the high luminosity end of the luminosity functions evolves more strongly than the low-luminosity end.

preprint2019arXiv

Radio Galaxy Zoo: Knowledge Transfer Using Rotationally Invariant Self-Organising Maps

With the advent of large scale surveys the manual analysis and classification of individual radio source morphologies is rendered impossible as existing approaches do not scale. The analysis of complex morphological features in the spatial domain is a particularly important task. Here we discuss the challenges of transferring crowdsourced labels obtained from the Radio Galaxy Zoo project and introduce a proper transfer mechanism via quantile random forest regression. By using parallelized rotation and flipping invariant Kohonen-maps, image cubes of Radio Galaxy Zoo selected galaxies formed from the FIRST radio continuum and WISE infrared all sky surveys are first projected down to a two-dimensional embedding in an unsupervised way. This embedding can be seen as a discretised space of shapes with the coordinates reflecting morphological features as expressed by the automatically derived prototypes. We find that these prototypes have reconstructed physically meaningful processes across two channel images at radio and infrared wavelengths in an unsupervised manner. In the second step, images are compared with those prototypes to create a heat-map, which is the morphological fingerprint of each object and the basis for transferring the user generated labels. These heat-maps have reduced the feature space by a factor of 248 and are able to be used as the basis for subsequent ML methods. Using an ensemble of decision trees we achieve upwards of 85.7% and 80.7% accuracy when predicting the number of components and peaks in an image, respectively, using these heat-maps. We also question the currently used discrete classification schema and introduce a continuous scale that better reflects the uncertainty in transition between two classes, caused by sensitivity and resolution limits.

preprint2015arXiv

A search of the Orion spur for continuous gravitational waves using a "loosely coherent" algorithm on data from LIGO interferometers

We report results of a wideband search for periodic gravitational waves from isolated neutron stars within the Orion spur towards both the inner and outer regions of our Galaxy. As gravitational waves interact very weakly with matter, the search is unimpeded by dust and concentrations of stars. One search disk (A) is $6.87^\circ$ in diameter and centered on $20^\textrm{h}10^\textrm{m}54.71^\textrm{s}+33^\circ33'25.29"$, and the other (B) is $7.45^\circ$ in diameter and centered on $8^\textrm{h}35^\textrm{m}20.61^\textrm{s}-46^\circ49'25.151"$. We explored the frequency range of 50-1500 Hz and frequency derivative from $0$ to $-5\times 10^{-9}$ Hz/s. A multi-stage, loosely coherent search program allowed probing more deeply than before in these two regions, while increasing coherence length with every stage. Rigorous followup parameters have winnowed initial coincidence set to only 70 candidates, to be examined manually. None of those 70 candidates proved to be consistent with an isolated gravitational wave emitter, and 95% confidence level upper limits were placed on continuous-wave strain amplitudes. Near $169$ Hz we achieve our lowest 95% CL upper limit on worst-case linearly polarized strain amplitude $h_0$ of $6.3\times 10^{-25}$, while at the high end of our frequency range we achieve a worst-case upper limit of $3.4\times 10^{-24}$ for all polarizations and sky locations.

preprint2015arXiv

The ASKAP/EMU Source Finding Data Challenge

The Evolutionary Map of the Universe (EMU) is a proposed radio continuum survey of the Southern Hemisphere up to declination +30 deg., with the Australian Square Kilometre Array Pathfinder (ASKAP). EMU will use an automated source identification and measurement approach that is demonstrably optimal, to maximise the reliability, utility and robustness of the resulting radio source catalogues. As part of the process of achieving this aim, a "Data Challenge" has been conducted, providing international teams the opportunity to test a variety of source finders on a set of simulated images. The aim is to quantify the accuracy of existing automated source finding and measurement approaches, and to identify potential limitations. The Challenge attracted nine independent teams, who tested eleven different source finding tools. In addition, the Challenge initiators also tested the current ASKAPsoft source-finding tool to establish how it could benefit from incorporating successful features of the other tools. Here we present the results of the Data Challenge, identifying the successes and limitations for this broad variety of the current generation of radio source finding tools. As expected, most finders demonstrate completeness levels close to 100% at 10sigma dropping to levels around 10% by 5sigma. The reliability is typically close to 100% at 10sigma, with performance to lower sensitivities varying greatly between finders. All finders demonstrate the usual trade-off between completeness and reliability, whereby maintaining a high completeness at low signal-to-noise comes at the expense of reduced reliability, and vice-versa. We conclude with a series of recommendations for improving the performance of the ASKAPsoft source-finding tool.

preprint2015arXiv

The Spitzer Archival Far-InfraRed Extragalactic Survey

We present the Spitzer Archival Far-InfraRed Extragalactic Survey (SAFIRES). This program produces refined mosaics and source lists for all far-infrared extragalactic data taken during the more than six years of the cryogenic operation of the Spitzer Space Telescope. The SAFIRES products consist of far-infrared data in two wavelength bands (70 um and 160 um) across approximately 180 square degrees of sky, with source lists containing far-infrared fluxes for almost 40,000 extragalactic point sources. Thus, SAFIRES provides a large, robust archival far-infrared data set suitable for many scientific goals.

preprint2014arXiv

A directed search for gravitational waves from Scorpius X-1 with initial LIGO

We present results of a search for continuously-emitted gravitational radiation, directed at the brightest low-mass X-ray binary, Scorpius X-1. Our semi-coherent analysis covers 10 days of LIGO S5 data ranging from 50-550 Hz, and performs an incoherent sum of coherent $\mathcal{F}$-statistic power distributed amongst frequency-modulated orbital sidebands. All candidates not removed at the veto stage were found to be consistent with noise at a 1% false alarm rate. We present Bayesian 95% confidence upper limits on gravitational-wave strain amplitude using two different prior distributions: a standard one, with no a priori assumptions about the orientation of Scorpius X-1; and an angle-restricted one, using a prior derived from electromagnetic observations. Median strain upper limits of 1.3e-24 and 8e-25 are reported at 150 Hz for the standard and angle-restricted searches respectively. This proof of principle analysis was limited to a short observation time by unknown effects of accretion on the intrinsic spin frequency of the neutron star, but improves upon previous upper limits by factors of ~1.4 for the standard, and 2.3 for the angle-restricted search at the sensitive region of the detector.

preprint2014arXiv

Application of a Hough search for continuous gravitational waves on data from the 5th LIGO science run

We report on an all-sky search for periodic gravitational waves in the frequency range $\mathrm{50-1000 Hz}$ with the first derivative of frequency in the range $-8.9 \times 10^{-10}$ Hz/s to zero in two years of data collected during LIGO's fifth science run. Our results employ a Hough transform technique, introducing a $χ^2$ test and analysis of coincidences between the signal levels in years 1 and 2 of observations that offers a significant improvement in the product of strain sensitivity with compute cycles per data sample compared to previously published searches. Since our search yields no surviving candidates, we present results taking the form of frequency dependent, 95$%$ confidence upper limits on the strain amplitude $h_0$. The most stringent upper limit from year 1 is $1.0\times 10^{-24}$ in the $\mathrm{158.00-158.25 Hz}$ band. In year 2, the most stringent upper limit is $\mathrm{8.9\times10^{-25}}$ in the $\mathrm{146.50-146.75 Hz}$ band. This improved detection pipeline, which is computationally efficient by at least two orders of magnitude better than our flagship Einstein$@$Home search, will be important for "quick-look" searches in the Advanced LIGO and Virgo detector era.

preprint2014arXiv

Characterization of the LIGO detectors during their sixth science run

In 2009-2010, the Laser Interferometer Gravitational-wave Observa- tory (LIGO) operated together with international partners Virgo and GEO600 as a network to search for gravitational waves of astrophysical origin. The sensitiv- ity of these detectors was limited by a combination of noise sources inherent to the instrumental design and its environment, often localized in time or frequency, that couple into the gravitational-wave readout. Here we review the performance of the LIGO instruments during this epoch, the work done to characterize the de- tectors and their data, and the effect that transient and continuous noise artefacts have on the sensitivity of LIGO to a variety of astrophysical sources.

preprint2014arXiv

Constraints on cosmic strings from the LIGO-Virgo gravitational-wave detectors

Cosmic strings can give rise to a large variety of interesting astrophysical phenomena. Among them, powerful bursts of gravitational waves (GWs) produced by cusps are a promising observational signature. In this Letter we present a search for GWs from cosmic string cusps in data collected by the LIGO and Virgo gravitational wave detectors between 2005 and 2010, with over 625 days of live time. We find no evidence of GW signals from cosmic strings. From this result, we derive new constraints on cosmic string parameters, which complement and improve existing limits from previous searches for a stochastic background of GWs from cosmic microwave background measurements and pulsar timing data. In particular, if the size of loops is given by the gravitational backreaction scale, we place upper limits on the string tension $Gμ$ below $10^{-8}$ in some regions of the cosmic string parameter space.

preprint2014arXiv

First all-sky search for continuous gravitational waves from unknown sources in binary systems

We present the first results of an all-sky search for continuous gravitational waves from unknown spinning neutron stars in binary systems using LIGO and Virgo data. Using a specially developed analysis program, the TwoSpect algorithm, the search was carried out on data from the sixth LIGO Science Run and the second and third Virgo Science Runs. The search covers a range of frequencies from 20 Hz to 520 Hz, a range of orbital periods from 2 to ~2,254 h and a frequency- and period-dependent range of frequency modulation depths from 0.277 to 100 mHz. This corresponds to a range of projected semi-major axes of the orbit from ~0.6e-3 ls to ~6,500 ls assuming the orbit of the binary is circular. While no plausible candidate gravitational wave events survive the pipeline, upper limits are set on the analyzed data. The most sensitive 95% confidence upper limit obtained on gravitational wave strain is 2.3e-24 at 217 Hz, assuming the source waves are circularly polarized. Although this search has been optimized for circular binary orbits, the upper limits obtained remain valid for orbital eccentricities as large as 0.9. In addition, upper limits are placed on continuous gravitational wave emission from the low-mass x-ray binary Scorpius X-1 between 20 Hz and 57.25 Hz.

preprint2014arXiv

Gravitational waves from known pulsars: results from the initial detector era

We present the results of searches for gravitational waves from a large selection of pulsars using data from the most recent science runs (S6, VSR2 and VSR4) of the initial generation of interferometric gravitational wave detectors LIGO (Laser Interferometric Gravitational-wave Observatory) and Virgo. We do not see evidence for gravitational wave emission from any of the targeted sources but produce upper limits on the emission amplitude. We highlight the results from seven young pulsars with large spin-down luminosities. We reach within a factor of five of the canonical spin-down limit for all seven of these, whilst for the Crab and Vela pulsars we further surpass their spin-down limits. We present new or updated limits for 172 other pulsars (including both young and millisecond pulsars). Now that the detectors are undergoing major upgrades, and, for completeness, we bring together all of the most up-to-date results from all pulsars searched for during the operations of the first-generation LIGO, Virgo and GEO600 detectors. This gives a total of 195 pulsars including the most recent results described in this paper.

preprint2014arXiv

Implementation of an F-statistic all-sky search for continuous gravitational waves in Virgo VSR1 data

We present an implementation of the $\mathcal{F}$-statistic to carry out the first search in data from the Virgo laser interferometric gravitational wave detector for periodic gravitational waves from a priori unknown, isolated rotating neutron stars. We searched a frequency $f_0$ range from 100 Hz to 1 kHz and the frequency dependent spindown $f_1$ range from $-1.6\,(f_0/100\,{\rm Hz}) \times 10^{-9}\,$ Hz/s to zero. A large part of this frequency - spindown space was unexplored by any of the all-sky searches published so far. Our method consisted of a coherent search over two-day periods using the $\mathcal{F}$-statistic, followed by a search for coincidences among the candidates from the two-day segments. We have introduced a number of novel techniques and algorithms that allow the use of the Fast Fourier Transform (FFT) algorithm in the coherent part of the search resulting in a fifty-fold speed-up in computation of the $\mathcal{F}$-statistic with respect to the algorithm used in the other pipelines. No significant gravitational wave signal was found. The sensitivity of the search was estimated by injecting signals into the data. In the most sensitive parts of the detector band more than 90% of signals would have been detected with dimensionless gravitational-wave amplitude greater than $5 \times 10^{-24}$.

preprint2014arXiv

Improved Upper Limits on the Stochastic Gravitational-Wave Background from 2009-2010 LIGO and Virgo Data

Gravitational waves from a variety of sources are predicted to superpose to create a stochastic background. This background is expected to contain unique information from throughout the history of the universe that is unavailable through standard electromagnetic observations, making its study of fundamental importance to understanding the evolution of the universe. We carry out a search for the stochastic background with the latest data from LIGO and Virgo. Consistent with predictions from most stochastic gravitational-wave background models, the data display no evidence of a stochastic gravitational-wave signal. Assuming a gravitational-wave spectrum of Omega_GW(f)=Omega_alpha*(f/f_ref)^alpha, we place 95% confidence level upper limits on the energy density of the background in each of four frequency bands spanning 41.5-1726 Hz. In the frequency band of 41.5-169.25 Hz for a spectral index of alpha=0, we constrain the energy density of the stochastic background to be Omega_GW(f)<5.6x10^-6. For the 600-1000 Hz band, Omega_GW(f)<0.14*(f/900 Hz)^3, a factor of 2.5 lower than the best previously reported upper limits. We find Omega_GW(f)<1.8x10^-4 using a spectral index of zero for 170-600 Hz and Omega_GW(f)<1.0*(f/1300 Hz)^3 for 1000-1726 Hz, bands in which no previous direct limits have been placed. The limits in these four bands are the lowest direct measurements to date on the stochastic background. We discuss the implications of these results in light of the recent claim by the BICEP2 experiment of the possible evidence for inflationary gravitational waves.

preprint2014arXiv

Methods and results of a search for gravitational waves associated with gamma-ray bursts using the GEO600, LIGO, and Virgo detectors

In this paper we report on a search for short-duration gravitational wave bursts in the frequency range 64 Hz-1792 Hz associated with gamma-ray bursts (GRBs), using data from GEO600 and one of the LIGO or Virgo detectors. We introduce the method of a linear search grid to analyse GRB events with large sky localisation uncertainties such as the localisations provided by the Fermi Gamma-ray Burst Monitor (GBM). Coherent searches for gravitational waves (GWs) can be computationally intensive when the GRB sky position is not well-localised, due to the corrections required for the difference in arrival time between detectors. Using a linear search grid we are able to reduce the computational cost of the analysis by a factor of O(10) for GBM events. Furthermore, we demonstrate that our analysis pipeline can improve upon the sky localisation of GRBs detected by the GBM, if a high-frequency GW signal is observed in coincidence. We use the linear search grid method in a search for GWs associated with 129 GRBs observed satellite-based gamma-ray experiments between 2006 and 2011. The GRBs in our sample had not been previously analysed for GW counterparts. A fraction of our GRB events are analysed using data from GEO600 while the detector was using squeezed-light states to improve its sensitivity; this is the first search for GWs using data from a squeezed-light interferometric observatory. We find no evidence for GW signals, either with any individual GRB in this sample or with the population as a whole. For each GRB we place lower bounds on the distance to the progenitor, assuming a fixed GW emission energy of $10^{-2} M_{\odot}c^{2}$, with a median exclusion distance of 0.8 Mpc for emission at 500 Hz and 0.3 Mpc at 1 kHz. The reduced computational cost associated with a linear search grid will enable rapid searches for GWs associated with Fermi GBM events in the Advanced detector era.

preprint2014arXiv

Multimessenger Search for Sources of Gravitational Waves and High-Energy Neutrinos: Results for Initial LIGO-Virgo and IceCube

We report the results of a multimessenger search for coincident signals from the LIGO and Virgo gravitational-wave observatories and the partially completed IceCube high-energy neutrino detector, including periods of joint operation between 2007-2010. These include parts of the 2005-2007 run and the 2009-2010 run for LIGO-Virgo, and IceCube's observation periods with 22, 59 and 79 strings. We find no significant coincident events, and use the search results to derive upper limits on the rate of joint sources for a range of source emission parameters. For the optimistic assumption of gravitational-wave emission energy of $10^{-2}$\,M$_\odot$c$^2$ at $\sim 150$\,Hz with $\sim 60$\,ms duration, and high-energy neutrino emission of $10^{51}$\,erg comparable to the isotropic gamma-ray energy of gamma-ray bursts, we limit the source rate below $1.6 \times 10^{-2}$\,Mpc$^{-3}$yr$^{-1}$. We also examine how combining information from gravitational waves and neutrinos will aid discovery in the advanced gravitational-wave detector era.

preprint2014arXiv

Narrow-band search of continuous gravitational-wave signals from Crab and Vela pulsars in Virgo VSR4 data

In this paper we present the results of a coherent narrow-band search for continuous gravitational-wave signals from the Crab and Vela pulsars conducted on Virgo VSR4 data. In order to take into account a possible small mismatch between the gravitational wave frequency and two times the star rotation frequency, inferred from measurement of the electromagnetic pulse rate, a range of 0.02 Hz around two times the star rotational frequency has been searched for both the pulsars. No evidence for a signal has been found and 95$\%$ confidence level upper limits have been computed both assuming polarization parameters are completely unknown and that they are known with some uncertainty, as derived from X-ray observations of the pulsar wind torii. For Vela the upper limits are comparable to the spin-down limit, computed assuming that all the observed spin-down is due to the emission of gravitational waves. For Crab the upper limits are about a factor of two below the spin-down limit, and represent a significant improvement with respect to past analysis. This is the first time the spin-down limit is significantly overcome in a narrow-band search.

preprint2014arXiv

Search for gravitational wave ringdowns from perturbed intermediate mass black holes in LIGO-Virgo data from 2005-2010

We report results from a search for gravitational waves produced by perturbed intermediate mass black holes (IMBH) in data collected by LIGO and Virgo between 2005 and 2010. The search was sensitive to astrophysical sources that produced damped sinusoid gravitational wave signals, also known as ringdowns, with frequency $50\le f_{0}/\mathrm{Hz} \le 2000$ and decay timescale $0.0001\lesssim τ/\mathrm{s} \lesssim 0.1$ characteristic of those produced in mergers of IMBH pairs. No significant gravitational wave candidate was detected. We report upper limits on the astrophysical coalescence rates of IMBHs with total binary mass $50 \le M/\mathrm{M}_\odot \le 450$ and component mass ratios of either 1:1 or 4:1. For systems with total mass $100 \le M/\mathrm{M}_\odot \le 150$, we report a 90%-confidence upper limit on the rate of binary IMBH mergers with non-spinning and equal mass components of $6.9\times10^{-8}\,$Mpc$^{-3}$yr$^{-1}$. We also report a rate upper limit for ringdown waveforms from perturbed IMBHs, radiating 1% of their mass as gravitational waves in the fundamental, $\ell=m=2$, oscillation mode, that is nearly three orders of magnitude more stringent than previous results.

preprint2014arXiv

Search for gravitational waves associated with gamma-ray bursts detected by the InterPlanetary Network

We present the results of a search for gravitational waves associated with 223 gamma-ray bursts (GRBs) detected by the InterPlanetary Network (IPN) in 2005-2010 during LIGO's fifth and sixth science runs and Virgo's first, second and third science runs. The IPN satellites provide accurate times of the bursts and sky localizations that vary significantly from degree scale to hundreds of square degrees. We search for both a well-modeled binary coalescence signal, the favored progenitor model for short GRBs, and for generic, unmodeled gravitational wave bursts. Both searches use the event time and sky localization to improve the gravitational-wave search sensitivity as compared to corresponding all-time, all-sky searches. We find no evidence of a gravitational-wave signal associated with any of the IPN GRBs in the sample, nor do we find evidence for a population of weak gravitational-wave signals associated with the GRBs. For all IPN-detected GRBs, for which a sufficient duration of quality gravitational-wave data is available, we place lower bounds on the distance to the source in accordance with an optimistic assumption of gravitational-wave emission energy of $10^{-2}M_{\odot}c^2$ at 150 Hz, and find a median of 13 Mpc. For the 27 short-hard GRBs we place 90% confidence exclusion distances to two source models: a binary neutron star coalescence, with a median distance of 12Mpc, or the coalescence of a neutron star and black hole, with a median distance of 22 Mpc. Finally, we combine this search with previously published results to provide a population statement for GRB searches in first-generation LIGO and Virgo gravitational-wave detectors, and a resulting examination of prospects for the advanced gravitational-wave detectors.

preprint2014arXiv

Searching for stochastic gravitational waves using data from the two co-located LIGO Hanford detectors

Searches for a stochastic gravitational-wave background (SGWB) using terrestrial detectors typically involve cross-correlating data from pairs of detectors. The sensitivity of such cross-correlation analyses depends, among other things, on the separation between the two detectors: the smaller the separation, the better the sensitivity. Hence, a co-located detector pair is more sensitive to a gravitational-wave background than a non-co-located detector pair. However, co-located detectors are also expected to suffer from correlated noise from instrumental and environmental effects that could contaminate the measurement of the background. Hence, methods to identify and mitigate the effects of correlated noise are necessary to achieve the potential increase in sensitivity of co-located detectors. Here we report on the first SGWB analysis using the two LIGO Hanford detectors and address the complications arising from correlated environmental noise. We apply correlated noise identification and mitigation techniques to data taken by the two LIGO Hanford detectors, H1 and H2, during LIGO's fifth science run. At low frequencies, 40 - 460 Hz, we are unable to sufficiently mitigate the correlated noise to a level where we may confidently measure or bound the stochastic gravitational-wave signal. However, at high frequencies, 460-1000 Hz, these techniques are sufficient to set a $95%$ confidence level (C.L.) upper limit on the gravitational-wave energy density of Ω(f)<7.7 x 10^{-4} (f/ 900 Hz)^3, which improves on the previous upper limit by a factor of $\sim 180$. In doing so, we demonstrate techniques that will be useful for future searches using advanced detectors, where correlated noise (e.g., from global magnetic fields) may affect even widely separated detectors.

preprint2014arXiv

The NINJA-2 project: Detecting and characterizing gravitational waveforms modelled using numerical binary black hole simulations

The Numerical INJection Analysis (NINJA) project is a collaborative effort between members of the numerical relativity and gravitational-wave astrophysics communities. The purpose of NINJA is to study the ability to detect gravitational waves emitted from merging binary black holes and recover their parameters with next-generation gravitational-wave observatories. We report here on the results of the second NINJA project, NINJA-2, which employs 60 complete binary black hole hybrid waveforms consisting of a numerical portion modelling the late inspiral, merger, and ringdown stitched to a post-Newtonian portion modelling the early inspiral. In a "blind injection challenge" similar to that conducted in recent LIGO and Virgo science runs, we added 7 hybrid waveforms to two months of data recolored to predictions of Advanced LIGO and Advanced Virgo sensitivity curves during their first observing runs. The resulting data was analyzed by gravitational-wave detection algorithms and 6 of the waveforms were recovered with false alarm rates smaller than 1 in a thousand years. Parameter estimation algorithms were run on each of these waveforms to explore the ability to constrain the masses, component angular momenta and sky position of these waveforms. We also perform a large-scale monte-carlo study to assess the ability to recover each of the 60 hybrid waveforms with early Advanced LIGO and Advanced Virgo sensitivity curves. Our results predict that early Advanced LIGO and Advanced Virgo will have a volume-weighted average sensitive distance of 300Mpc (1Gpc) for $10M_{\odot}+10M_{\odot}$ ($50M_{\odot}+50M_{\odot}$) binary black hole coalescences. We demonstrate that neglecting the component angular momenta in the waveform models used in matched-filtering will result in a reduction in sensitivity for systems with large component angular momenta. [Abstract abridged for ArXiv, full version in PDF]

preprint2013arXiv

A directed search for continuous Gravitational Waves from the Galactic Center

We present the results of a directed search for continuous gravitational waves from unknown, isolated neutron stars in the Galactic Center region, performed on two years of data from LIGO's fifth science run from two LIGO detectors. The search uses a semi-coherent approach, analyzing coherently 630 segments, each spanning 11.5 hours, and then incoherently combining the results of the single segments. It covers gravitational wave frequencies in a range from 78 to 496 Hz and a frequency-dependent range of first order spindown values down to -7.86 x 10^-8 Hz/s at the highest frequency. No gravitational waves were detected. We place 90% confidence upper limits on the gravitational wave amplitude of sources at the Galactic Center. Placing 90% confidence upper limits on the gravitational wave amplitude of sources at the Galactic Center, we reach ~3.35x10^-25 for frequencies near 150 Hz. These upper limits are the most constraining to date for a large-parameter-space search for continuous gravitational wave signals.

preprint2013arXiv

Far-Infrared Properties of Type 1 Quasars

We use the Spitzer Space Telescope Enhanced Imaging Products (SEIP) and the Spitzer Archival Far-InfraRed Extragalactic Survey (SAFIRES) to study the spectral energy distributions of spectroscopically confirmed type 1 quasars selected from the Sloan Digital Sky Survey (SDSS). By combining the Spitzer and SDSS data with the 2-Micron All Sky Survey (2MASS) we are able to construct a statistically robust rest-frame 0.1-100 micron type 1 quasar template. We find the quasar population is well-described by a single power-law SED at wavelengths less than 20 microns, in good agreement with previous work. However, at longer wavelengths we find a significant excess in infrared luminosity above an extrapolated power-law, along with signifiant object-to-object dispersion in the SED. The mean excess reaches a maximum of 0.8 dex at rest-frame wavelengths near 100 microns.

preprint2013arXiv

First Searches for Optical Counterparts to Gravitational-wave Candidate Events

During the LIGO and Virgo joint science runs in 2009-2010, gravitational wave (GW) data from three interferometer detectors were analyzed within minutes to select GW candidate events and infer their apparent sky positions. Target coordinates were transmitted to several telescopes for follow-up observations aimed at the detection of an associated optical transient. Images were obtained for eight such GW candidates. We present the methods used to analyze the image data as well as the transient search results. No optical transient was identified with a convincing association with any of these candidates, and none of the GW triggers showed strong evidence for being astrophysical in nature. We compare the sensitivities of these observations to several model light curves from possible sources of interest, and discuss prospects for future joint GW-optical observations of this type.

preprint2013arXiv

Search for long-lived gravitational-wave transients coincident with long gamma-ray bursts

Long gamma-ray bursts (GRBs) have been linked to extreme core-collapse supernovae from massive stars. Gravitational waves (GW) offer a probe of the physics behind long GRBs. We investigate models of long-lived (~10-1000s) GW emission associated with the accretion disk of a collapsed star or with its protoneutron star remnant. Using data from LIGO's fifth science run, and GRB triggers from the swift experiment, we perform a search for unmodeled long-lived GW transients. Finding no evidence of GW emission, we place 90% confidence level upper limits on the GW fluence at Earth from long GRBs for three waveforms inspired by a model of GWs from accretion disk instabilities. These limits range from F<3.5 ergs cm^-2 to $F<1200 ergs cm^-2, depending on the GRB and on the model, allowing us to probe optimistic scenarios of GW production out to distances as far as ~33 Mpc. Advanced detectors are expected to achieve strain sensitivities 10x better than initial LIGO, potentially allowing us to probe the engines of the nearest long GRBs.

preprint2012arXiv

Spitzer and Herschel-based SEDs of 24um-bright z~0.3-3.0 starbursts and obscured quasars

In this paper, we characterize the infrared spectral energy distributions (SEDs) of mid-IR selected z~0.3-3.0 and L_IR~10^11-10^13Lsun galaxies, and study how their SEDs differ from those of local and high-z analogs. Our mid-IR flux-limited sample of 191 sources is unique in size, and spectral coverage, including Spitzer mid-IR spectroscopy. Here we add Herschel photometry at 250um, 350um, and 500um, which allows us to obtain accurate total IR luminosities, as well as constrain the relative contributions of AGN and starbursts to those luminosities. Our sample constitutes ~23% AGN (i.e. where the AGN contributes >50% of L_IR), ~30% starbursts (where AGN contributes <20% of L_IR and the mid-IR spectra are starburst-like); and ~47% composites (which show both significant AGN and starburst activity). The AGN-dominated sources divide into ones that show a strong silicate 9.7um absorption feature, implying highly obscured systems, and ones that do not. The high-tau_9.7 sources are half of our z>1.2 AGN, but show SEDs that are extremely rare among local AGN. The SEDs of our z~2 starburst-dominated ULIRGs are much closer to those of local LIRGs than ULIRGs. This is consistent with our earlier finding that, unlike local ULIRGs, our high-z starbursts are typically only in the early stages of a merger. The SEDs of the composite sources are most similar to the local archetypal warm ULIRG, Mrk231. In summary, our results show that there is strong evolution in the SEDs between local and z~2 IR-luminous galaxies, as well as that there is a wide range of SEDs among high redshift IR-luminous sources. The publicly-available SED templates we derive from our sample will be particularly useful for infrared population synthesis models, as well as in the interpretation of other mid-IR high-z galaxies in particular those detected by the recent all sky WISE survey.

preprint2012arXiv

The Spitzer Extragalactic Representative Volume Survey (SERVS): survey definition and goals

We present the Spitzer Extragalactic Representative Volume Survey (SERVS), an 18 square degrees medium-deep survey at 3.6 and 4.5 microns with the post-cryogenic Spitzer Space Telescope to ~2 microJy (AB=23.1) depth of five highly observed astronomical fields (ELAIS-N1, ELAIS-S1, Lockman Hole, Chandra Deep Field South and XMM-LSS). SERVS is designed to enable the study of galaxy evolution as a function of environment from z~5 to the present day, and is the first extragalactic survey both large enough and deep enough to put rare objects such as luminous quasars and galaxy clusters at z>1 into their cosmological context. SERVS is designed to overlap with several key surveys at optical, near- through far-infrared, submillimeter and radio wavelengths to provide an unprecedented view of the formation and evolution of massive galaxies. In this paper, we discuss the SERVS survey design, the data processing flow from image reduction and mosaicing to catalogs, as well as coverage of ancillary data from other surveys in the SERVS fields. We also highlight a variety of early science results from the survey.

preprint2011arXiv

GOODS-Herschel: an infrared main sequence for star-forming galaxies

We present the deepest far-IR observations obtained with Herschel and examine the 3-500um SEDs of galaxies at 0<z<2.5, supplemented by a local reference sample from IRAS, ISO, Spitzer and AKARI data. We find that the ratio of total IR luminosity to rest-frame 8um luminosity, IR8 (=Lir/L8), follows a Gaussian distribution centered on IR8=4 and defines an IR main sequence (MS). A minority population (<20 %) of outliers producing a tail skewed toward higher values of IR8 consist of starbursts (SB) with compact projected star formation densities. IR8 can be used to separate galaxies with normal and extended modes of star formation from compact SBs with high-IR8, high projected IR surface brightness (>3x10^10 Lsun kpc^-2) and a high specific SFR (i.e., SBs). The rest-frame, UV-2700A size of these distant SBs is typically half that of MS galaxies, supporting the correlation between star formation density and SB activity that is measured for the local sample. Locally, (U)LIRGs are systematically in the SB mode, whereas most distant (U)LIRGs form stars in the "normal" MS mode. This confusion between two modes of star formation is the cause of the so-called "mid-IR excess" population of galaxies found at z>1.5 by previous studies. MS galaxies have strong PAH emission line features, a broad far-IR bump resulting from a combination of dust temperatures (Tdust~15-50 K), and an effective Tdust~31 K, as derived from the peak wavelength of their IR SED. Galaxies in the SB regime instead exhibit weak PAH EW and a sharper far-IR bump with an effective Tdust~40 K. Finally, we present evidence that the mid-to-far IR emission of X-ray AGNs is predominantly produced by star formation and that candidate dusty AGNs with a power-law emission in the mid-IR systematically occur in compact, dusty SBs. After correcting for the effect of SBs on IR8, we identify new candidates for extremely obscured AGNs.

preprint2011arXiv

Planck early results XIV: ERCSC validation and extreme radio sources

Planck's all sky surveys at 30-857 GHz provide an unprecedented opportunity to follow the radio spectra of a large sample of extragalactic sources to frequencies 2-20 times higher than allowed by past, large area, ground-based surveys. We combine the results of the Planck Early Release Compact Source Catalog (ERCSC) with quasi-simultaneous ground-based observations, as well as archival data, at frequencies below or overlapping Planck frequency bands, to validate the astrometry and photometry of the ERCSC radio sources and study the spectral features shown in this new frequency window opened by Planck. The ERCSC source positions and flux density scales are found to be consistent with the ground-based observations. We present and discuss the spectral energy distributions (SEDs) of a sample of "extreme" radio sources to illustrate the richness of the ERCSC for the study of extragalactic radio sources. Variability is found to play a role in the unusual spectral features of some of these sources.

preprint2011arXiv

Planck Early Results. V. The Low Frequency Instrument data processing

We describe the processing of data from the Low Frequency Instrument (LFI) used in production of the Planck Early Release Compact Source Catalogue (ERCSC). In particular, we discuss the steps involved in reducing the data from telemetry packets to cleaned, calibrated, time-ordered data (TOD) and frequency maps. Data are continuously calibrated using the modulation of the temperature of the cosmic microwave background radiation induced by the motion of the spacecraft. Noise properties are estimated from TOD from which the sky signal has been removed using a generalized least square map-making algorithm. Measured 1/f noise knee-frequencies range from 100mHz at 30GHz to a few tens of mHz at 70GHz. A destriping code (Madam) is employed to combine radiometric data and pointing information into sky maps, minimizing the variance of correlated noise. Noise covariance matrices required to compute statistical uncertainties on LFI and Planck products are also produced. Main beams are estimated down to the approx -10dB level using Jupiter transits, which are also used for geometrical calibration of the focal plane.

preprint2011arXiv

Planck Early Results. VII. The Early Release Compact Source Catalog

(Abridged) A brief description of the methodology of construction, contents and usage of the Planck Early Release Compact Source Catalogue (ERCSC), including the Early Cold Cores (ECC) and the Early Sunyaev-Zeldovich (ESZ) cluster catalogue is provided. The catalogue is based on data that consist of mapping the entire sky once and 60% of the sky a second time by Planck, thereby comprising the first high sensitivity radio/submillimetre observations of the entire sky. A Monte-Carlo algorithm based on the injection and extraction of artificial sources into the Planck maps was implemented to select reliable sources among all extracted candidates such that the cumulative reliability of the catalogue is >=90%. The 10sigma photometric flux density limit of the catalogue at |b|>30 deg is 0.49, 1.0, 0.67, 0.5, 0.33, 0.28, 0.25, 0.47 and 0.82 Jy at each of the nine frequencies between 30 and 857 GHz. Sources which are up to a factor of ~2 fainter than this limit, and which are present in "clean" regions of the Galaxy where the sky background due to emission from the interstellar medium is low, are included in the ERCSC if they meet the high reliability criterion. The Planck ERCSC sources have known associations to stars with dust shells, stellar cores, radio galaxies, blazars, infrared luminous galaxies and Galactic interstellar medium features. A significant fraction of unclassified sources are also present in the catalogs. In addition, two early release catalogs that contain 915 cold molecular cloud core candidates and 189 SZ cluster candidates that have been generated using multi-frequency algorithms are presented. The entire source list, with more than 15000 unique sources, is ripe for follow-up characterisation with Herschel, ATCA, VLA, SOFIA, ALMA and other ground-based observing facilities.

preprint2010arXiv

The ATLAS Survey of the CDFS and ELAIS-S1 Fields

The first phase of the ATLAS (Australia Telescope Large Area Survey) project surveyed a total 7 square degrees down to 30 micro Jy rms at 1.4 GHz and is the largest sensitive radio survey ever attempted. We report on the scientific achievements of ATLAS to date and plans to extend the project as a path finder for the proposed EMU (Evolutionary map of the Universe) project which has been designed to use ASKAP (Australian Square Kilometre Array Pathfinder).

preprint2009arXiv

Characterising the Far-infrared Properties of Distant X-ray Detected AGNs: Evidence for Evolution in the Infrared--X-ray Luminosity Ratio

(Abridged) We investigate the far-infrared properties of X-ray sources detected in the Chandra Deep Field-South (CDF-S) survey using the ultra-deep 70um and 24um Spitzer observations taken in this field. We rely on stacking analyses of the 70um data to characterise the average 70um properties of the X-ray sources. Using Spitzer-IRS data of the Swift-BAT sample of z~0 active galactic nuclei (hereafter, AGNs), we show that the 70um/24um flux ratio can distinguish between AGN-dominated and starburst-dominated systems out to z~1.5. From stacking analysis we find that both high redshift and z~0 AGNs follow the same tendency toward warmer 70um/24um colours with increasing X-ray luminosity (L_X). We also show that the 70um flux can be used to determine the infrared (8-1000um) luminosities of high redshift AGNs. We use this information to show that L_X=10^{42-43} erg/s AGNs at high redshifts (z=1-2) have infrared to X-ray luminosity ratios (hereafter, L_IR/L_X) that are, on average, 4.7_{-2.0}^{+10.2} and 12.7^{+7.1}_{-2.6} times higher than AGNs with similar X-ray luminosities at z=0.5-1 and z~0, respectively. By comparison, we find that the L_IR/L_X ratios of L_X=10^{43-44} erg/s AGNs remain largely unchanged across this same redshift interval. We explore the consequences that these results may have on the identification of distant, potentially Compton thick AGNs using L_IR/L_X ratios. We discuss possible scenarios for the observed increase in the L_IR/L_X ratio with redshift, including changes in the dust covering factor of AGNs and/or the star formation rates of their host galaxies. Finally, we show how deep observations to be undertaken by the Herschel Space Observatory will enable us to discriminate between these proposed scenarios and also identify Compton-thick AGNs at high redshifts.

preprint2009arXiv

Investigating the Far-IR/Radio Correlation of Star Forming Galaxies to z=3

In order to examine the far-IR/radio correlation at high redshift we have studied the Spitzer 70um/far-infrared (far-IR) properties of sub-mJy radio sources from the 13^H XMM-Newton/Chandra Deep Field by redshift and galaxy type: active galactic nucleus (AGN) or star forming galaxy (SFG). We directly detect 70um counterparts (at >3sigma significance) for 22.5% (92/408) of the radio sources, while for the rest we perform stacking analysis by redshift and galaxy type. For the sources detected at 70um we find that the median and scatter of the observed flux density ratio, q70, are similar to previous results in the literature, but with a slight decrease in q70 towards higher redshifts. Of the radio sources detected at 70um 8/92 were already classified as AGN, but two of which maybe SFGs. For the stacked sources we obtain a detection for the SFGs at every redshift bin which implies they have mean flux densities a factor ~5 below the original 70um detection limit. For the stacked AGN we obtain a detection only in our highest redshift bin (1<~z<~5) where we may be sampling hot dust associated with the AGN at rest-frame 12-35um. The combined observed mean value of q70 for the SFGs (detected and non-detected at 70um) decreases gradually with redshift, consistent with tracks derived from empirical spectral energy distributions (SEDs) of local SFGs. Upon closer inspection and when comparing with tracks of appropriate luminosity, the values of q70 broadly agree at low redshift. However, at z~1, the observed q70 (for ULIRGs) is 2sigma below the value seen for local ULIRGs tracks, implying a difference in the SED between local and z~1 ULIRGs. At higher redshifts, the convergence of the tracks and the higher uncertainties in q70 do not allow us to determine if such a difference persists.

preprint2009arXiv

The Comoving Infrared Luminosity Density: Domination of Cold Galaxies across 0<z<1

In this paper we examine the contribution of galaxies with different infrared (IR) spectral energy distributions (SEDs) to the comoving infrared luminosity density, a proxy for the comoving star formation rate (SFR) density. We characterise galaxies as having either a cold or hot IR SED depending upon whether the rest-frame wavelength of their peak IR energy output is above or below 90um. Our work is based on a far-IR selected sample both in the local Universe and at high redshift, the former consisting of IRAS 60um-selected galaxies at z<0.07 and the latter of Spitzer 70um selected galaxies across 0.1<z<1. We find that the total IR luminosity densities for each redshift/luminosity bin agree well with results derived from other deep mid/far-IR surveys. At z<0.07 we observe the previously known results: that moderate luminosity galaxies (L_IR<10^11 Lsun) dominate the total luminosity density and that the fraction of cold galaxies decreases with increasing luminosity, becoming negligible at the highest luminosities. Conversely, above z=0.1 we find that luminous IR galaxies (L_IR>10^11 Lsun), the majority of which are cold, dominate the IR luminosity density. We therefore infer that cold galaxies dominate the IR luminosity density across the whole 0<z<1 range, hence appear to be the main driver behind the increase in SFR density up to z~1 whereas local luminous galaxies are not, on the whole, representative of the high redshift population.

preprint2009arXiv

The Link Between SCUBA and Spitzer: Cold Galaxies at z<1

We show that the far-IR properties of distant Luminous and Ultraluminous InfraRed Galaxies (LIRGs and ULIRGs) are on average divergent from analogous sources in the local Universe. Our analysis is based on Spitzer MIPS and IRAC data of L_IR>10^10 L_solar, 70um-selected objects in the 0.1<z<2 redshift range and supported by a comparison with the IRAS Bright Galaxy Sample. The majority of the objects in our sample are described by Spectral Energy Distributions (SEDs) which peak at longer wavelengths than local sources of equivalent total infrared luminosity. This shift in SED peak wavelength implies a noticeable change in the dust and/or star-forming properties from z~0 to the early Universe, tending towards lower dust temperatures, indicative of strong evolution in the cold dust, `cirrus', component. We show that these objects are potentially the missing link between the well-studied local IR-luminous galaxies, Spitzer IR populations and SCUBA sources -- the z<1 counterparts of the cold z>1 SubMillimetre Galaxies (SMGs) discovered in blank-field submillimetre surveys. The Herschel Space Observatory is well placed to fully characterise the nature of these objects, as its coverage extends over a major part of the far-IR/submm SED for a wide redshift range.

preprint2007arXiv

Multiwavelength study of massive galaxies at z~2. I. Star formation and galaxy growth

Examining a sample of massive galaxies at 1.4<z<2.5 with K_{Vega}<22 from the Great Observatories Origins Deep Survey, we compare photometry from Spitzer at mid- and far-IR, to submillimeter, radio and rest-frame ultraviolet wavelengths, to test the agreement between different tracers of star formation rates (SFRs) and to explore the implications for galaxy assembly. For z~2 galaxies with moderate luminosities(L_{8um}<10^{11}L_sun), we find that the SFR can be estimated consistently from the multiwavelength data based on local luminosity correlations. However,20--30% of massive galaxies, and nearly all those with L_{8um}>10^{11}L_sun, show a mid-IR excess which is likely due to the presence of obscured active nuclei, as shown in a companion paper. There is a tight and roughly linear correlation between stellar mass and SFR for 24um-detected galaxies. For a given mass, the SFR at z=2 was larger by a factor of ~4 and ~30 relative to that in star forming galaxies at z=1 and z=0, respectively. Typical ultraluminous infrared galaxies (ULIRGs) at z=2 are relatively 'transparent' to ultraviolet light, and their activity is long lived (~400 Myr), unlike that in local ULIRGs and high redshift submillimeter-selected galaxies. ULIRGs are the common mode of star formation in massive galaxies at z=2, and the high duty cycle suggests that major mergers are not the dominant trigger for this activity.Current galaxy formation models underpredict the normalization of the mass-SFR correlation by about a factor of 4, and the space density of ULIRGs by an orderof magnitude, but give better agreement for z>1.4 quiescent galaxies.