Search

In order to study transient phenomena in the Universe, existing and forthcoming imaging surveys are covering wide areas of sky repeatedly over time, with a range of cadences, point spread functions, and depths. We describe here a framework that allows an efficient search for different types of time-varying astrophysical phenomena in current and future, large data repositories. We first present a methodology to generate and store key survey parameters that enable researchers to determine if a survey, or a combination of surveys, allows specific time-variable astrophysical phenomena to be discovered. To facilitate further exploration of sources in regions of interest, we then generate a few sample metrics that capture the essential brightness characteristics of a sky pixel at a specific wavelength. Together, we refer to these as "annotated coadds". The techniques presented here for WISE/NEOWISE-R data are sensitive to 10 percent brightness variations at around 12th Vega magnitude at 4.5 microns wavelength. Application of the technique to ZTF data also enabled the detection of 0.5 mag variability at 20 AB mag in the r-band. We demonstrate the capabilities of these metrics for different classes of sources: high proper-motion stars, periodic variable stars, and supernovae, and find that each metric has its advantages depending on the nature of variability. We also present a data structure which will ease the search for temporally varying phenomena in future surveys., Comment: 26 pages, 11 figures, PASP in press

The final explosive fate of massive stars, and the nature of the compact remnants they leave behind (black holes and neutron stars), are major open questions in astrophysics. Many massive stars are stripped of their outer hydrogen envelopes as they evolve. Such Wolf-Rayet (W-R) stars emit strong and rapidly expanding (v_wind>1000 km/s) winds indicating a high escape velocity from the stellar surface. A fraction of this population is also helium depleted, with spectra dominated by highly-ionized emission lines of carbon and oxygen (Types WC/WO). Evidence indicates that the most commonly-observed supernova (SN) explosions that lack hydrogen and helium (Types Ib/Ic) cannot result from massive WC/WO stars, leading some to suggest that most such stars collapse directly into black holes without a visible supernova explosions. Here, we present observations of supernova SN 2019hgp, discovered about a day after explosion. The short rise time and rapid decline place it among an emerging population of rapidly-evolving transients (RETs). Spectroscopy reveals a rich set of emission lines indicating that the explosion occurred within a nebula composed of carbon, oxygen, and neon. Narrow absorption features show that this material is expanding at relatively high velocities (>1500 km/s) requiring a compact progenitor. Our observations are consistent with an explosion of a massive WC/WO star, and suggest that massive W-R stars may be the progenitors of some rapidly evolving transients., Comment: Unedited author version, Nature in press

Early observations of Type Ia supernovae (SNe$\,$Ia) provide essential clues for understanding the progenitor system that gave rise to the terminal thermonuclear explosion. We present exquisite observations of SN$\,$2019yvq, the second observed SN$\,$Ia, after iPTF$\,$14atg, to display an early flash of emission in the ultraviolet (UV) and optical. Our analysis finds that SN$\,$2019yvq was unusual, even when ignoring the initial flash, in that it was moderately underluminous for an SN$\,$Ia ($M_g \approx -18.5\,$mag at peak) yet featured very high absorption velocities ($v \approx 15,000\,\mathrm{km\,s}^{-1}$ for Si II $\lambda$6355 at peak). We find that many of the observational features of SN$\,$2019yvq, aside from the flash, can be explained if the explosive yield of radioactive $^{56}\mathrm{Ni}$ is relatively low (we measure $M_{^{56}\mathrm{Ni}} = 0.31 \pm 0.05\,M_\odot$) and it and other iron-group elements are concentrated in the innermost layers of the ejecta. To explain both the UV/optical flash and peak properties of SN$\,$2019yvq we consider four different models: interaction between the SN ejecta and a nondegenerate companion, extended clumps of $^{56}\mathrm{Ni}$ in the outer ejecta, a double-detonation explosion, and the violent merger of two white dwarfs. Each of these models has shortcomings when compared to the observations; it is clear additional tuning is required to better match SN$\,$2019yvq. In closing, we predict that the nebular spectra of SN$\,$2019yvq will feature either H or He emission, if the ejecta collided with a companion, strong [Ca II] emission, if it was a double detonation, or narrow [O I] emission, if it was due to a violent merger., Comment: 27 pages, 14 figures, accepted in ApJ

We present optical, radio, and X-ray observations of SN2020bvc (=ASASSN20bs; ZTF20aalxlis), a nearby ($z=0.0252$; $d$=114 Mpc) broad-lined (BL) Type Ic supernova (SN). Our observations show that SN2020bvc shares several properties in common with the Ic-BL SN2006aj, which was associated with the low-luminosity gamma-ray burst (LLGRB) 060218. First, the 10 GHz radio light curve is on the faint end of LLGRB-SNe ($L_\mathrm{radio} \approx 10^{37}$erg/s): we model our VLA observations (spanning 13-43 d) as synchrotron emission from a mildly relativistic ($v \gtrsim 0.3c$) forward shock. Second, with Swift and Chandra we detect X-ray emission ($L_X \approx 10^{41}$erg/s) that is not naturally explained as inverse Compton emission or as part of the same synchrotron spectrum as the radio emission. Third, high-cadence ($6\times$/night) data from the Zwicky Transient Facility (ZTF) shows a double-peaked optical light curve, the first peak from shock-cooling emission from extended low-mass material (mass $M<10^{-2} M_\odot$ at radius $R>10^{12}$cm) and the second peak from the radioactive decay of Ni-56. SN2020bvc is the first confirmed double-peaked Ic-BL SN discovered without a GRB trigger, and shows X-ray and radio emission similar to LLGRB-SNe: this is consistent with models in which the same mechanism produces both the LLGRB and the shock-cooling emission. For four of the five other nearby ($z\lesssim0.05$) Ic-BL SNe with ZTF high-cadence data, we rule out a first peak like that seen in SN2006aj and SN2020bvc, i.e. that lasts $\approx 1$d and reaches a peak luminosity $M \approx -18$. X-ray and radio follow-up observations of future such events will establish whether double-peaked optical light curves are indeed predictive of LLGRB-like X-ray and radio emission., Comment: Accepted to ApJ. 39 pages, 17 figures. Comments welcome

We present ZTF18abvkwla (the "Koala"), a fast blue optical transient discovered in the Zwicky Transient Facility (ZTF) One-Day Cadence (1DC) Survey. ZTF18abvkwla has a number of features in common with the groundbreaking transient AT2018cow: blue colors at peak ($g-r\approx-0.5$ mag), a short rise time from half-max of under two days, a decay time to half-max of only three days, a high optical luminosity ($M_{g,\mathrm{peak}}\approx-20.6$mag), a hot ($\gtrsim 40,000$K) featureless spectrum at peak light, and a luminous radio counterpart. At late times ($\Delta t>80$d) the radio luminosity of ZTF18abvkwla ($\nu L_\nu \gtrsim 10^{40}$erg/s at 10 GHz, observer-frame) is most similar to that of long-duration gamma-ray bursts (GRBs). The host galaxy is a dwarf starburst galaxy ($M\approx5\times10^{8}M_\odot$, $\mathrm{SFR}\approx7 M_\odot$/yr) that is moderately metal-enriched ($\log\mathrm{[O/H]} \approx 8.5$), similar to the hosts of GRBs and superluminous supernovae. As in AT2018cow, the radio and optical emission in ZTF18abvkwla likely arise from two separate components: the radio from fast-moving ejecta ($\Gamma \beta c >0.38c$) and the optical from shock-interaction with confined dense material ($<0.07M_\odot$ in $\sim 10^{15}$cm). Compiling transients in the literature with $t_\mathrm{rise} <5$d and $M_\mathrm{peak}<-20$mag, we find that a significant number are engine-powered, and suggest that the high peak optical luminosity is directly related to the presence of this engine. From 18 months of the 1DC survey, we find that transients in this rise-luminosity phase space are at least two to three orders of magnitude less common than CC SNe. Finally, we discuss strategies for identifying such events with future facilities like the Large Synoptic Survey Telescope, and prospects for detecting accompanying X-ray and radio emission., Comment: 27 pages, 11 figures. Resubmitted to ApJ on 13 April following comments by referee. Comments welcome!

While it is clear that Type Ia supernovae (SNe) are the result of thermonuclear explosions in C/O white dwarfs (WDs), a great deal remains uncertain about the binary companion that facilitates the explosive disruption of the WD. Here, we present a comprehensive analysis of a large, unique data set of 127 SNe$\,$Ia with exquisite coverage by the Zwicky Transient Facility (ZTF). High-cadence (six observations per night) ZTF observations allow us to measure the SN rise time and examine its initial evolution. We develop a Bayesian framework to model the early rise as a power law in time, which enables the inclusion of priors in our model. For a volume-limited subset of normal SNe$\,$Ia, we find that the mean power-law index is consistent with 2 in the $r_\mathrm{ZTF}$-band ($\alpha_r = 2.01\pm0.02$), as expected in the expanding fireball model. There are, however, individual SNe that are clearly inconsistent with $\alpha_r=2$. We estimate a mean rise time of 18.9$\,$d (with a range extending from $\sim$15 to 22$\,$d), though this is subject to the adopted prior. We identify an important, previously unknown, bias whereby the rise times for higher-redshift SNe within a flux-limited survey are systematically underestimated. This effect can be partially alleviated if the power-law index is fixed to $\alpha=2$, in which case we estimate a mean rise time of 21.7$\,$d (with a range from $\sim$18 to 23$\,$d). The sample includes a handful of rare and peculiar SNe$\,$Ia. Finally, we conclude with a discussion of lessons learned from the ZTF sample that can eventually be applied to observations from the Vera C. Rubin Observatory., Comment: 31 pages, 17 figures; ApJ accepted

We present photometry and spectroscopy of four hydrogen-poor luminous supernovae discovered during the two-month science commissioning and early operations of the Zwicky Transient Facility (ZTF) survey. Three of these objects, SN2018bym (ZTF18aapgrxo), SN2018avk (ZTF18aaisyyp) and SN2018bgv (ZTF18aavrmcg) resemble typical SLSN-I spectroscopically, while SN2018don (ZTF18aajqcue) may be an object similar to SN2007bi experiencing considerable host galaxy reddening, or an intrinsically long-lived, luminous and red SN Ic. We analyze the light curves, spectra, and host galaxy properties of these four objects and put them in context of the population of SLSN-I. SN2018bgv stands out as the fastest-rising SLSN-I observed to date, with a rest-frame g-band rise time of just 10 days from explosion to peak -- if it is powered by magnetar spin-down, the implied ejecta mass is only ~1 M$_{\odot}$. SN2018don also displays unusual properties -- in addition to its red colors and comparatively massive host galaxy, the light curve undergoes some of the strongest light curve undulations post-peak seen in a SLSN-I, which we speculate may be due to interaction with circumstellar material. We discuss the promises and challenges of finding SLSNe in large-scale surveys like ZTF given the observed diversity in the population., Comment: Updated to match published version. Expanded discussion following referee comments; conclusions unchanged

The Astropy project supports and fosters the development of open-source and openly-developed Python packages that provide commonly-needed functionality to the astronomical community. A key element of the Astropy project is the core package Astropy, which serves as the foundation for more specialized projects and packages. In this article, we provide an overview of the organization of the Astropy project and summarize key features in the core package as of the recent major release, version 2.0. We then describe the project infrastructure designed to facilitate and support development for a broader ecosystem of inter-operable packages. We conclude with a future outlook of planned new features and directions for the broader Astropy project., Comment: Minor changes to author list and title. Comments and feedback welcome through the paper source repository: https://github.com/astropy/astropy-v2.0-paper For more information about Astropy, see http://www.astropy.org/

We report the discovery and analysis of 36 new eclipsing EL CVn-type binaries, consisting of a core helium-composition pre-white dwarf and an early-type main-sequence companion, more than doubling the known population of these systems. We have used supervised machine learning methods to search 0.8 million lightcurves from the Palomar Transient Factory, combined with SDSS, Pan-STARRS and 2MASS colours. The new systems range in orbital periods from 0.46-3.8 d and in apparent brightness from ~14-16 mag in the PTF $R$ or $g^{\prime}$ filters. For twelve of the systems, we obtained radial velocity curves with the Intermediate Dispersion Spectrograph at the Isaac Newton Telescope. We modelled the lightcurves, radial velocity curves and spectral energy distributions to determine the system parameters. The radii (0.3-0.7 $\mathrm{R_{\odot}}$) and effective temperatures (8000-17000 K) of the pre-He-WDs are consistent with stellar evolution models, but the masses (0.12-0.28 $\mathrm{M_{\odot}}$) show more variance than models predicted. This study shows that using machine learning techniques on large synoptic survey data is a powerful way to discover substantial samples of binary systems in short-lived evolutionary stages.

We used wide area surveys over 39 deg$^2$ by the HerMES collaboration, performed with the Herschel Observatory SPIRE multi-wavelength camera, to estimate the low-redshift, $0.02

We present a five-band Herschel study (100-500um) of three galaxy clusters at z~1.2 from the Spitzer Adaptation of the Red-Sequence Cluster Survey (SpARCS). With a sample of 120 spectroscopically-confirmed cluster members, we investigate the role of environment on galaxy properties utilizing the projected cluster phase space (line-of-sight velocity versus clustercentric radius), which probes the time-averaged galaxy density to which a galaxy has been exposed. We divide cluster galaxies into phase-space bins of (r/r200) x (v/sigma_v), tracing a sequence of accretion histories in phase space. Stacking optically star-forming cluster members on the Herschel maps, we measure average infrared star formation rates, and, for the first time in high-redshift galaxy clusters, dust temperatures for dynamically distinct galaxy populations---namely, recent infalls and those that were accreted onto the cluster at an earlier epoch. Proceeding from the infalling to virialized (central) regions of phase space, we find a steady decrease in the specific star formation rate and increase in the stellar age of star-forming cluster galaxies. We perform a probability analysis to investigate all acceptable infrared spectral energy distributions within the full parameter space and measure a ~4 sigma drop in the average dust temperature of cluster galaxies in an intermediate phase-space bin, compared to an otherwise flat trend with phase space. We suggest one plausible quenching mechanism which may be consistent with these trends, invoking ram-pressure stripping of the warmer dust for galaxies within this intermediate accretion phase., Comment: 14 pages, 7 figures, 2 tables; accepted for publication in ApJ

We present a method for selecting $z>4$ dusty, star forming galaxies (DSFGs) using Herschel/SPIRE 250/350/500 $\mu m$ flux densities to search for red sources. We apply this method to 21 deg$^2$ of data from the HerMES survey to produce a catalog of 38 high-$z$ candidates. Follow-up of the first 5 of these sources confirms that this method is efficient at selecting high-$z$ DSFGs, with 4/5 at $z=4.3$ to $6.3$ (and the remaining source at $z=3.4$), and that they are some of the most luminous dusty sources known. Comparison with previous DSFG samples, mostly selected at longer wavelengths (e.g., 850 $\mu m$) and in single-band surveys, shows that our method is much more efficient at selecting high-$z$ DSFGs, in the sense that a much larger fraction are at $z>3$. Correcting for the selection completeness and purity, we find that the number of bright ($S_{500\,\mu m} \ge 30$ mJy), red Herschel sources is $3.3 \pm 0.8$ deg$^{-2}$. This is much higher than the number predicted by current models, suggesting that the DSFG population extends to higher redshifts than previously believed. If the shape of the luminosity function for high-$z$ DSFGs is similar to that at $z\sim2$, rest-frame UV based studies may be missing a significant component of the star formation density at $z=4$ to $6$, even after correction for extinction., Comment: Accepted for publication in ApJ

The old, red stars which constitute the bulges of galaxies, and the massive black holes at their centres, are the relics of a period in cosmic history when galaxies formed stars at remarkable rates and active galactic nuclei (AGN) shone brightly from accretion onto black holes. It is widely suspected, but unproven, that the tight correlation in mass of the black hole and stellar components results from the AGN quenching the surrounding star formation as it approaches its peak luminosity. X-rays trace emission from AGN unambiguously, while powerful star-forming galaxies are usually dust-obscured and are brightest at infrared to submillimetre wavelengths. Here we report observations in the submillimetre and X-ray which show that rapid star formation was common in the host galaxies of AGN when the Universe was 2-6 Gyrs old, but that the most vigorous star formation is not observed around black holes above an X-ray luminosity of 10^44 erg/s. This suppression of star formation in the host galaxies of powerful AGN is a key prediction of models in which the AGN drives a powerful outflow, expelling the interstellar medium of its host galaxy and transforming the galaxy's properties in a brief period of cosmic time., Comment: Published in Nature, May 2012. Includes supplementary information

We describe the procedure used to flux calibrate the three-band submillimetre photometer in the Spectral and Photometric Imaging REceiver (SPIRE) instrument on the Herschel Space Observatory. This includes the equations describing the calibration scheme, a justification for using Neptune as the primary calibration source, a description of the observations and data processing procedures used to derive flux calibration parameters (for converting from voltage to flux density) for every bolometer in each array, an analysis of the error budget in the flux calibration for the individual bolometers, and tests of the flux calibration on observations of primary and secondary calibrators. The procedure for deriving the flux calibration parameters is divided into two parts. In the first part, we use observations of astronomical sources in conjunction with the operation of the photometer internal calibration source to derive the unscaled derivatives of the flux calibration curves. To scale the calibration curves in Jy/beam/V, we then use observations of Neptune in which the beam of each bolometer is mapped using Neptune observed in a very fine scan pattern. The total instrumental uncertainties in the flux calibration for the individual bolometers is ~0.5% for most bolometers, although a few bolometers have uncertainties of ~1-5% because of issues with the Neptune observations. Based on application of the flux calibration parameters to Neptune observations performed using typical scan map observing modes, we determined that measurements from each array as a whole have instrumental uncertainties of 1.5%. This is considerably less than the absolute calibration uncertainty associated with the model of Neptune, which is estimated at 4%., Comment: 20 pages, 17 figures, accepted for publication in MNRAS

We have observed four massive galaxy clusters with the SPIRE instrument on the Herschel Space Observatory and measure a deficit of surface brightness within their central region after subtracting sources. We simulate the effects of instrumental sensitivity and resolution, the source population, and the lensing effect of the clusters to estimate the shape and amplitude of the deficit. The amplitude of the central deficit is a strong function of the surface density and flux distribution of the background sources. We find that for the current best fitting faint end number counts, and excellent lensing models, the most likely amplitude of the central deficit is the full intensity of the cosmic infrared background (CIB). Our measurement leads to a lower limit to the integrated total intensity of the CIB of I(250 microns) > 0.69_(-0.03)^(+0.03) (stat.)_(-0.06)^(+0.11) (sys.) MJy/sr, with more CIB possible from both low-redshift sources and from sources within the target clusters. It should be possible to observe this effect in existing high angular resolution data at other wavelengths where the CIB is bright, which would allow tests of models of the faint source component of the CIB., Comment: 5 pages, 3 figures, replacing with version matching published manuscript

We exploit the deep and extended far infrared data sets (at 70, 100 and 160 um) of the Herschel GTO PACS Evolutionary Probe (PEP) Survey, in combination with the HERschel Multi tiered Extragalactic Survey (HerMES) data at 250, 350 and 500 um, to derive the evolution of the restframe 35 um, 60 um, 90 um, and total infrared (IR) luminosity functions (LFs) up to z~4. We detect very strong luminosity evolution for the total IR LF combined with a density evolution. In agreement with previous findings, the IR luminosity density increases steeply to z~1, then flattens between z~1 and z~3 to decrease at z greater than 3. Galaxies with different SEDs, masses and sSFRs evolve in very different ways and this large and deep statistical sample is the first one allowing us to separately study the different evolutionary behaviours of the individual IR populations contributing to the IR luminosity density. Galaxies occupying the well established SFR/stellar mass main sequence (MS) are found to dominate both the total IR LF and luminosity density at all redshifts, with the contribution from off MS sources (0.6 dex above MS) being nearly constant (~20% of the total IR luminosity density) and showing no significant signs of increase with increasing z over the whole 0.8

Using Herschel data from the deepest SPIRE and PACS surveys (HerMES and PEP) in COSMOS and GOODS (N+S), we examine the dust properties of IR-luminous (L_IR>10^10 L_sun) galaxies at 0.145K) SEDs and cold (T<25K), cirrus-dominated SEDs are rare, with most sources being within the range occupied by warm starbursts such as M82 and cool spirals such as M51. We observe a luminosity-temperature (L-T) relation, where the average dust temperature of log[L_IR/L_sun]=12.5 galaxies is about 10K higher than that of their log[L_IR/L_sun]=10.5 counterparts. However, although the increased dust heating in more luminous systems is the driving factor behind the L-T relation, the increase in dust mass and/or starburst size with luminosity plays a dominant role in shaping it. Our results show that the dust conditions in IR-luminous sources evolve with cosmic time: at high redshift, dust temperatures are on average up to 10K lower than what is measured locally. This is manifested as a flattening of the L-T relation, suggesting that (U)LIRGs in the early Universe are typically characterised by a more extended dust distribution and/or higher dust masses than local equivalent sources. Interestingly, the evolution in dust temperature is luminosity dependent, with the fraction of LIRGs with T<35K showing a 2-fold increase from z~0 to z~2, whereas that of ULIRGs with T<35K shows a 6-fold increase., Comment: 27 pages, 23 figures, 2 tables, accepted for publication in MNRAS

The infrared (IR) emission of M_* galaxies (10^{10.4} < M_{star} < 10^{11.0} M_\sun) in galaxy pairs, derived using data obtained in Herschel (PEP/HerMES) and Spitzer (S-COSMOS) surveys, is compared to that of single disk galaxies in well matched control samples to study the cosmic evolution of the star-formation enhancement induced by galaxy-galaxy interaction. Both the mean IR SED and mean IR luminosity of star-forming galaxies (SFGs) in SFG+SFG (S+S) pairs in the redshift bin of 0.6 < z < 1 are consistent with no star-formation enhancement. SFGs in S+S pairs in a lower redshift bin of 0.2 < z < 0.6 show marginal evidence for a weak star-formation enhancement. Together with the significant and strong sSFR enhancement shown by SFGs in a local sample of S+S pairs (obtained using previously published Spitzer observations), our results reveal a trend for the star-formation enhancement in S+S pairs to decrease with increasing redshift. Between z=0 and z=1, this decline of interaction-induced star-formation enhancement occurs in parallel with the dramatic increase (by a factor of ~10) of the sSFR of single SFGs, both can be explained by the higher gas fraction in higher z disks. SFGs in mixed pairs (S+E pairs) do not show any significant star-formation enhancement at any redshift. The difference between SFGs in S+S pairs and in S+E pairs suggests a modulation of the sSFR by the inter-galactic medium IGM in the dark matter halos (DMH) hosting these pairs., Comment: 7 pages, 4 figures; accepted by ApJ

We present measurements of the auto- and cross-frequency power spectra of the cosmic infrared background (CIB) at 250, 350, and 500um (1200, 860, and 600 GHz) from observations totaling ~ 70 deg^2 made with the SPIRE instrument aboard the Herschel Space Observatory. We measure a fractional anisotropy dI / I = 14 +- 4%, detecting signatures arising from the clustering of dusty star-forming galaxies in both the linear (2-halo) and non-linear (1-halo) regimes; and that the transition from the 2- to 1-halo terms, below which power originates predominantly from multiple galaxies within dark matter halos, occurs at k_theta ~ 0.1 - 0.12 arcmin^-1 (l ~ 2160 - 2380), from 250 to 500um. New to this paper is clear evidence of a dependence of the Poisson and 1-halo power on the flux-cut level of masked sources --- suggesting that some fraction of the more luminous sources occupy more massive halos as satellites, or are possibly close pairs. We measure the cross-correlation power spectra between bands, finding that bands which are farthest apart are the least correlated, as well as hints of a reduction in the correlation between bands when resolved sources are more aggressively masked. In the second part of the paper we attempt to interpret the measurements in the framework of the halo model. With the aim of fitting simultaneously with one model the power spectra, number counts, and absolute CIB level in all bands, we find that this is achievable by invoking a luminosity-mass relationship, such that the luminosity-to-mass ratio peaks at a particular halo mass scale and declines towards lower and higher mass halos. Our best-fit model finds that the halo mass which is most efficient at hosting star formation in the redshift range of peak star-forming activity, z ~ 1-3, is log(M_peak/M_sun) ~ 12.1 +- 0.5, and that the minimum halo mass to host infrared galaxies is log(M_min/M_sun) ~ 10.1 +- 0.6., Comment: 29 pages, 15 figures, 10 tables. Accepted for publication in the Astrophysical Journal. Tabulated data and other software available at http://www.astro.caltech.edu/~viero/viero_homepage/toolbox.html

We present a list of 13 candidate gravitationally lensed submillimeter galaxies (SMGs) from 95 square degrees of the Herschel Multi-tiered Extragalactic Survey, a surface density of 0.14\pm0.04deg^{-2}. The selected sources have 500um flux densities (S_500) greater than 100mJy. Gravitational lensing is confirmed by follow-up observations in 9 of the 13 systems (70%), and the lensing status of the four remaining sources is undetermined. We also present a supplementary sample of 29 (0.31\pm0.06deg^{-2}) gravitationally lensed SMG candidates with S_500=80--100mJy, which are expected to contain a higher fraction of interlopers than the primary candidates. The number counts of the candidate lensed galaxies are consistent with a simple statistical model of the lensing rate, which uses a foreground matter distribution, the intrinsic SMG number counts, and an assumed SMG redshift distribution. The model predicts that 32--74% of our S_500>100mJy candidates are strongly gravitationally lensed (mu>2), with the brightest sources being the most robust; this is consistent with the observational data. Our statistical model also predicts that, on average, lensed galaxies with S_500=100mJy are magnified by factors of ~9, with apparently brighter galaxies having progressively higher average magnification, due to the shape of the intrinsic number counts. 65% of the sources are expected to have intrinsic 500micron flux densities less than 30mJy. Thus, samples of strongly gravitationally lensed SMGs, such as those presented here, probe below the nominal Herschel detection limit at 500 micron. They are good targets for the detailed study of the physical conditions in distant dusty, star-forming galaxies, due to the lensing magnification, which can lead to spatial resolutions of ~0.01" in the source plane., Comment: ApJ in press. 31 pages, 16 figures, 5 tables. This version updated to match accepted version

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 \mu m), and Herschel-PACS (at 100 and 160 \mu 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\sigma 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., Comment: 23 pages, 17 figures, 9 Tables, MNRAS accepted

ABRIGED Herschel/SPIRE has provided confusion limited maps of deep fields at 250, 350, and 500um, as part of the HerMES survey. Due to confusion, only a small fraction of the Cosmic Infrared Background can be resolved into individually-detected sources. Our goal is to produce deep galaxy number counts and redshift distributions below the confusion limit, which we then use to place strong constraints on the origins of the cosmic infrared background and on models of galaxy evolution. We individually extracted the bright SPIRE with a method using the positions, the flux densities, and the redshifts of the 24um sources as a prior, and derived the number counts and redshift distributions of the bright SPIRE sources. For fainter SPIRE sources, we reconstructed the number counts and the redshift distribution below the confusion limit using the deep 24um catalogs associated with photometric redshift and information provided by the stacking of these sources into the deep SPIRE maps. Finally, by integrating all these counts, we studied the contribution of the galaxies to the CIB as a function of their flux density and redshift. Through stacking, we managed to reconstruct the source counts per redshift slice down to ~2 mJy in the three SPIRE bands, which lies about a factor 10 below the 5sigma confusion limit. None of the pre-existing population models are able to reproduce our results at better than 3sigma. Finally, we extrapolate our counts to zero flux density in order to derive an estimate of the total contribution of galaxies to the CIB, finding 10.1, 6.5, and 2.8 nW/m2/sr at 250, 350, and 500um, respectively. These values agree well with FIRAS absolute measurements, suggesting our number counts and their extrapolation are sufficient to explain the CIB. Finally, combining our results with other works, we estimate the energy budget contained in the CIB between 8 and 1000um: 26 nW/m2/sr., Comment: 25 pages, 17 figures, 15 tables, accepted by A&A, bug with Fig. 3 fixed

We investigate the potential of submm-mm and submm-mm-radio photometric redshifts using a sample of mm-selected sources as seen at 250, 350 and 500 {\mu}m by the SPIRE instrument on Herschel. From a sample of 63 previously identified mm-sources with reliable radio identifications in the GOODS-N and Lockman Hole North fields 46 (73 per cent) are found to have detections in at least one SPIRE band. We explore the observed submm/mm colour evolution with redshift, finding that the colours of mm-sources are adequately described by a modified blackbody with constant optical depth {\tau} = ({\nu}/{\nu}0)^{\beta} where {\beta} = +1.8 and {\nu}0 = c/100 {\mu}m. We find a tight correlation between dust temperature and IR luminosity. Using a single model of the dust temperature and IR luminosity relation we derive photometric redshift estimates for the 46 SPIRE detected mm-sources. Testing against the 22 sources with known spectroscopic, or good quality optical/near-IR photometric, redshifts we find submm/mm photometric redshifts offer a redshift accuracy of |z|/(1+z) = 0.16 (< |z| >= 0.51). Including constraints from the radio-far IR correlation the accuracy is improved to |z|/(1 + z) = 0.15 (< |z| >= 0.45). We estimate the redshift distribution of mm-selected sources finding a significant excess at z > 3 when compared to ~ 850 {\mu}m selected samples., Comment: 17 pages, 7 figures, accepted for publication in MNRAS

We present the results of a gravitational lensing analysis of the bright $\zs=2.957$ sub-millimeter galaxy (SMG), HERMES J105751.1+573027 found in {\it Herschel}/SPIRE Science Demonstration Phase data from the Herschel Multi-tiered Extragalactic Survey (HerMES) project. The high resolution imaging available in optical and Near-IR channels, along with CO emission obtained with the Plateau de Bure Interferometer, allow us to precisely estimate the intrinsic source extension and hence estimate the total lensing magnification to be $\mu=10.9\pm 0.7$. We measure the half-light radius $R_{\rm eff}$ of the source in the rest-frame Near-UV and $V$ bands that characterize the unobscured light coming from stars and find $R_{\rm eff,*}= [2.0 \pm 0.1]$ kpc, in good agreement with recent studies on the Submillimeter Galaxy population. This lens model is also used to estimate the size of the gas distribution ($R_{\rm eff,gas}= [1.1\pm0.5]$) kpc by mapping back in the source plane the CO (J=5-4) transition line emission. The lens modeling yields a relatively large Einstein radius $R_{\rm Ein}= 4\farcs10 \pm 0\farcs02$, corresponding to a deflector velocity dispersion of [$483\pm 16] \,\kms$. This shows that HERMES J105751.1+573027 is lensed by a {\it galaxy group-size} dark matter halo at redshift $\zl\sim 0.6$. The projected dark matter contribution largely dominates the mass budget within the Einstein radius with $f_{\rm dm}(

We report the discovery of a bright ($f(250\mum) > 400$ mJy), multiply-lensed submillimeter galaxy \obj\ in {\it Herschel}/SPIRE Science Demonstration Phase data from the HerMES project. Interferometric 880\mum\ Submillimeter Array observations resolve at least four images with a large separation of $\sim 9\arcsec$. A high-resolution adaptive optics $K_p$ image with Keck/NIRC2 clearly shows strong lensing arcs. Follow-up spectroscopy gives a redshift of $z=2.9575$, and the lensing model gives a total magnification of $\mu \sim 11 \pm 1$. The large image separation allows us to study the multi-wavelength spectral energy distribution (SED) of the lensed source unobscured by the central lensing mass. The far-IR/millimeter-wave SED is well described by a modified blackbody fit with an unusually warm dust temperature, $88 \pm 3$ K. We derive a lensing-corrected total IR luminosity of $(1.43 \pm 0.09) \times 10^{13}\, \mathrm{L}_{\odot}$, implying a star formation rate of $\sim 2500\, \mathrm{M}_{\odot}\, \mathrm{yr}^{-1}$. However, models primarily developed from brighter galaxies selected at longer wavelengths are a poor fit to the full optical-to-millimeter SED. A number of other strongly lensed systems have already been discovered in early {\it Herschel} data, and many more are expected as additional data are collected., Comment: Accepted for publication in ApJL

The extragalactic background light at far-infrared wavelengths originates from optically-faint, dusty, star-forming galaxies in the universe with star-formation rates at the level of a few hundred solar masses per year. Due to the relatively poor spatial resolution of far-infrared telescopes, the faint sub-millimetre galaxies are challenging to study individually. Instead, their average properties can be studied using statistics such as the angular power spectrum of the background intensity variations. A previous attempt at measuring this power spectrum resulted in the suggestion that the clustering amplitude is below the level computed with a simple ansatz based on a halo model. Here we report a clear detection of the excess clustering over the linear prediction at arcminute angular scales in the power spectrum of brightness fluctuations at 250, 350, and 500 microns. From this excess, we find that sub-millimetre galaxies are located in dark matter halos with a minimum mass of log[M_min/M_sun ]= 11.5^+0.7_-0.2 at 350 microns. This minimum dark matter halo mass corresponds to the most efficient mass scale for star formation in the universe, and is lower than that predicted by semi-analytical models for galaxy formation., Comment: Accepted for publication as a Letter to Editor in Nature; This paper is under press embargo until Nature publication (early February). Supplementary material and electronic tables will be available from http://cosmology.uci.edu at the time of publication. For HerMES team information see http://hermes.sussex.ac.uk/ and for data releases see http://hedam.oamp.fr/HerMES/

We describe the production and verification of sky maps of the five SPIRE fields observed as part of the Herschel Multi-tiered Extragalactic Survey (HerMES) during the Science Demonstration Phase (SDP) of the Herschel mission. We have implemented an iterative map-making algorithm (SHIM; The SPIRE-HerMES Iterative Mapper) to produce high fidelity maps that preserve extended diffuse emission on the sky while exploiting the repeated observations of the same region of the sky with many detectors in multiple scan directions to minimize residual instrument noise. We specify here the SHIM algorithm and outline the various tests that were performed to determine and characterize the quality of the maps and verify that the astrometry, point source flux and power on all relevant angular scales meets the needs of the HerMES science goals. These include multiple jackknife tests, determination of the map transfer function and detailed examination of the power spectra of both sky and jackknife maps. The map transfer function is approximately unity on scales from one arcminute to one degree. Final maps (v1.0), including multiple jackknives, as well as the SHIM pipeline, have been used by the HerMES team for the production of SDP papers., Comment: 9 pages, 9 figure, MNRAS Accepted, 27 Sept, 2010

Dusty, star forming galaxies contribute to a bright, currently unresolved cosmic far-infrared background. Deep Herschel-SPIRE images designed to detect and characterize the galaxies that comprise this background are highly confused, such that the bulk lies below the classical confusion limit. We analyze three fields from the HerMES programme in all three SPIRE bands (250, 350, and 500 microns); parameterized galaxy number count models are derived to a depth of ~2 mJy/beam, approximately 4 times the depth of previous analyses at these wavelengths, using a P(D) (probability of deflection) approach for comparison to theoretical number count models. Our fits account for 64, 60, and 43 per cent of the far-infrared background in the three bands. The number counts are consistent with those based on individually detected SPIRE sources, but generally inconsistent with most galaxy number counts models, which generically overpredict the number of bright galaxies and are not as steep as the P(D)-derived number counts. Clear evidence is found for a break in the slope of the differential number counts at low flux densities. Systematic effects in the P(D) analysis are explored. We find that the effects of clustering have a small impact on the data, and the largest identified systematic error arises from uncertainties in the SPIRE beam., Comment: Accepted for publication in MNRAS

We present the cross-identification and source photometry techniques used to process Herschel SPIRE imaging taken as part of the Herschel Multi-Tiered Extragalactic Survey (HerMES). Cross-identifications are performed in map-space so as to minimise source blending effects. We make use of a combination of linear inversion and model selection techniques to produce reliable cross-identification catalogues based on Spitzer MIPS 24 micron source positions. Testing on simulations and real Herschel observations show that this approach gives robust results for even the faintest sources S250~10 mJy. We apply our new technique to HerMES SPIRE observations taken as part of the science demostration phase of Herschel. For our real SPIRE observations we show that, for bright unconfused sources, our flux density estimates are in good agreement with those produced via more traditional point source detection methods (SussExtractor; Savage & Oliver et al. 2006) by Smith et al. 2010. When compared to the measured number density of sources in the SPIRE bands, we show that our method allows the recovery of a larger fraction of faint sources than these traditional methods. However this completeness is heavily dependant on the relative depth of the existing 24 micron catalogues and SPIRE imaging. Using our deepest multi-wavelength dataset in GOODS-N, we estimate that the use of shallow 24 micron in our other fields introduces an incompleteness at faint levels of between 20-40 per cent at 250 micron., Comment: Accepted for publication in MNRAS, 19 pages, 11 figures

The Spectral and Photometric Imaging Receiver (SPIRE) on Herschel has been carrying out deep extragalactic surveys, one of whose aims is to establish spectral energy distributions (SED)s of individual galaxies spanning the infrared/submillimeter (IR/SMM) wavelength region. We report observations of the (IR/SMM) emission from the Lockman North field (LN) and Great Observatories Origins Deep Survey field North (GOODS-N). Because galaxy images in the wavelength range covered by Herschel generally represent a blend with contributions from neighboring galaxies, we present sets of galaxies in each field especially free of blending at 250, 350, and 500 microns. We identify the cumulative emission of these galaxies and the fraction of the far infrared cosmic background radiation they contribute. Our surveys reveal a number of highly luminous galaxies at redshift z ~< 3 and a novel relationship between infrared and visible emission that shows a dependence on luminosity and redshift., Comment: 9 pages, 3 large tables, 2 tables in text, 5 figures Accepted for publication in MNRAS

We study the dust properties of galaxies in the redshift range 0.10.5 with L_IR>5x10^{10} L_\odot, appears to be 2-5 K colder than that of AKARI-selected local galaxies with similar luminosities; and the dispersion in T_dust for high-z galaxies increases with L_IR due to the existence of cold galaxies that are not seen among local galaxies. We show that this large dispersion of the L_IR-T_dust relation can bridge the gap between local star-forming galaxies and high-z submillimeter galaxies (SMGs). We also find that three SMGs with very low T_dust (<20 K) covered in this study have close neighbouring sources with similar 24-\mum brightness, which could lead to an overestimation of FIR/(sub)millimeter fluxes of the SMGs., Comment: 9 pages, 6 figures, To appear in MNRAS

The reliability of infrared (IR) and ultraviolet (UV) emissions to measure star formation rates in galaxies is investigated for a large sample of galaxies observed with the SPIRE and PACS instruments on Herschel as part of the HerMES project. We build flux-limited 250 micron samples of sources at redshift z<1, cross-matched with the Spitzer/MIPS and GALEX catalogues. About 60 % of the Herschel sources are detected in UV. The total IR luminosities, L_IR, of the sources are estimated using a SED-fitting code that fits to fluxes between 24 and 500 micron. Dust attenuation is discussed on the basis of commonly-used diagnostics: the L_IR/L_UV ratio and the slope, beta, of the UV continuum. A mean dust attenuation A_UV of ~ 3 mag is measured in the samples. L_IR/L_UV is found to correlate with L_IR. Galaxies with L_IR > 10 ^{11} L_sun and 0.5< z<1 exhibit a mean dust attenuation A_UV about 0.7 mag lower than that found for their local counterparts, although with a large dispersion. Our galaxy samples span a large range of beta and L_IR/L_UV values which, for the most part, are distributed between the ranges defined by the relations found locally for starburst and normal star-forming galaxies. As a consequence the recipe commonly applied to local starbursts is found to overestimate the dust attenuation correction in our galaxy sample by a factor ~2-3 ., Comment: accepted for publication in MNRAS, Herschel special issue

The Spectral and Photometric Imaging Receiver (SPIRE), is the Herschel Space Observatory`s submillimetre camera and spectrometer. It contains a three-band imaging photometer operating at 250, 350 and 500 microns, and an imaging Fourier Transform Spectrometer (FTS) which covers simultaneously its whole operating range of 194-671 microns (447-1550 GHz). The SPIRE detectors are arrays of feedhorn-coupled bolometers cooled to 0.3 K. The photometer has a field of view of 4' x 8', observed simultaneously in the three spectral bands. Its main operating mode is scan-mapping, whereby the field of view is scanned across the sky to achieve full spatial sampling and to cover large areas if desired. The spectrometer has an approximately circular field of view with a diameter of 2.6'. The spectral resolution can be adjusted between 1.2 and 25 GHz by changing the stroke length of the FTS scan mirror. Its main operating mode involves a fixed telescope pointing with multiple scans of the FTS mirror to acquire spectral data. For extended source measurements, multiple position offsets are implemented by means of an internal beam steering mirror to achieve the desired spatial sampling and by rastering of the telescope pointing to map areas larger than the field of view. The SPIRE instrument consists of a cold focal plane unit located inside the Herschel cryostat and warm electronics units, located on the spacecraft Service Module, for instrument control and data handling. Science data are transmitted to Earth with no on-board data compression, and processed by automatic pipelines to produce calibrated science products. The in-flight performance of the instrument matches or exceeds predictions based on pre-launch testing and modelling: the photometer sensitivity is comparable to or slightly better than estimated pre-launch, and the spectrometer sensitivity is also better by a factor of 1.5-2., Comment: Accepted for publication in Astronomy & Astrophyics (Herschel first results special issue)

We measure the angular correlation function, w(theta), from 0.5 to 30 arcminutes of detected sources in two wide fields of the Herschel Multi-tiered Extragalactic Survey (HerMES). Our measurements are consistent with the expected clustering shape from a population of sources that trace the dark matter density field, including non-linear clustering at arcminute angular scales arising from multiple sources that occupy the same dark matter halos. By making use of the halo model to connect the spatial clustering of sources to the dark matter halo distribution, we estimate source bias and halo occupation number for dusty sub-mm galaxies at z ~ 2. We find that sub-mm galaxies with 250 micron flux densities above 30 mJy reside in dark matter halos with mass above (5\pm4) x 10^12 M_sun, while (14\pm8)% of such sources appear as satellites in more massive halos., Comment: 5 pages; accepted for publication in A&A Special Issue on Herschel First Science Results (July 2010). HerMES information available at http://hermes.sussex.ac.uk/

The Herschel Space Observatory enables us to accurately measure the bolometric output of starburst galaxies and active galactic nuclei (AGN) by directly sampling the peak of their far-infrared (IR) emission. Here we examine whether the spectral energy distribution (SED) and dust temperature of galaxies have strongly evolved since z~2.5. We use Herschel deep extragalactic surveys from 100 to 500um to compute total IR luminosities in galaxies down to the faintest levels, using PACS and SPIRE in the GOODS-North field (PEP and HerMES key programs). We show that measurements in the SPIRE bands can be used below the statistical confusion limit if information at higher spatial resolution is used to identify isolated galaxies whose flux is not boosted by bright neighbors. Below z~1.5, mid-IR extrapolations are correct for star-forming galaxies with a dispersion of only 40% (0.15dex), therefore similar to z~0 galaxies. This narrow distribution is puzzling when considering the range of physical processes that could have affected the SED of these galaxies. Extrapolations from only one of the 160um, 250um or 350um bands alone tend to overestimate the total IR luminosity. This may be explained by the lack of far-IR constraints around and above ~150um (rest-frame) on local templates. We also note that the dust temperature of luminous IR galaxies around z~1 is mildly colder by 10-15% than their local analogs and up to 20% for ULIRGs at z~1.6. Above z=1.5, distant galaxies are found to exhibit a substantially larger mid- over far-IR ratio, which could either result from stronger broad emission lines or warm dust continuum heated by a hidden AGN. Two thirds of the AGNs identified in the field with a measured redshift exhibit the same behavior as purely star-forming galaxies. Hence a large fraction of AGNs harbor star formation at very high SFR and in conditions similar to purely star-forming galaxies., Comment: Astronomy and Astrophysics, Herschel Special Issue, in press as a Letter; 5 pages

We present colours of sources detected with the Herschel/SPIRE instrument in deep extragalactic surveys of the Lockman Hole, Spitzer-FLS, and GOODS-N fields in three photometric bands at 250, 350 and 500 micrometers. We compare these with expectations from the literature and discuss associated uncertainties and biases in the SPIRE data. We identify a 500 micrometer flux limited selection of sources from the HerMES point source catalogue that appears free from neighbouring/blended sources in all three SPIRE bands. We compare the colours with redshift tracks of various contemporary models. Based on these spectral templates we show that regions corresponding to specific population types and redshifts can be identified better in colour-flux space. The redshift tracks as well as the colour-flux plots imply a majority of detected objects with redshifts at 1

Local luminosity functions are fundamental benchmarks for high-redshift galaxy formation and evolution studies as well as for models describing these processes. Determining the local luminosity function in the submillimeter range can help to better constrain in particular the bolometric luminosity density in the local Universe, and Herschel offers the first opportunity to do so in an unbiased way by imaging large sky areas at several submillimeter wavelengths. We present the first Herschel measurement of the submillimeter 0

We report on the sensitivity of SPIRE photometers on the Herschel Space Observatory. Specifically, we measure the confusion noise from observations taken during the Science Demonstration Phase of the Herschel Multi-tiered Extragalactic Survey. Confusion noise is defined to be the spatial variation of the sky intensity in the limit of infinite integration time, and is found to be consistent among the different fields in our survey at the level of 5.8, 6.3 and 6.8 mJy/beam at 250, 350 and 500 microns, respectively. These results, together with the measured instrument noise, may be used to estimate the integration time required for confusion-limited maps, and provide a noise estimate for maps obtained by SPIRE., Comment: Astronomy and Astrophysics, Herschel Special Issue, in press as a Letter; 5 pages

We have carried out two extremely deep surveys with SPIRE, one of the two cameras on Herschel, at 250 microns, close to the peak of the far-infrared background. We have used the results to investigate the evolution of the rest-frame 250-micron luminosity function out to z=2. We find evidence for strong evolution out to a redshift of around 1 but evidence for at most weak evolution beyond this redshift. Our results suggest that a significant part of the stars and metals in the Universe today were formed at z<1.4 in spiral galaxies., Comment: Astronomy and Astrophysics, Herschel Special Issue, in press as a Letter; 5 pages

Nuclear and starburst activity are known to often occur concomitantly. Herschel-SPIRE provides sampling of the FIR SEDs of type 1 and type 2 AGN, allowing for the separation between the hot dust (torus) and cold dust (starburst) emission. We study large samples of spectroscopically confirmed type 1 and type 2 AGN lying within the Herschel Multi-tiered Extragalactic Survey (HerMES) fields observed during the science demonstration phase, aiming to understand their FIR colour distributions and constrain their starburst contributions. We find that one third of the spectroscopically confirmed AGN in the HerMES fields have 5-sigma detections at 250um, in agreement with previous (sub)mm AGN studies. Their combined Spitzer-MIPS and Herschel-SPIRE colours - specifically S(250)/S(70) vs. S(70)/S(24) - quite clearly separate them from the non-AGN, star-forming galaxy population, as their 24-um flux is dominated by the hot torus emission. However, their SPIRE colours alone do not differ from those of non-AGN galaxies. SED fitting shows that all those AGN need a starburst component to fully account for their FIR emission. For objects at z > 2, we find a correlation between the infrared luminosity attributed to the starburst component, L(SB), and the AGN accretion luminosity, L(acc), with L(SB) propto L(acc)^0.35. Type 2 AGN detected at 250um show on average higher L(SB) than type 1 objects but their number is still too low to establish whether this trend indicates stronger star-formation activity., Comment: Astronomy and Astrophysics, Herschel Special Issue, in press as a Letter; 5 pages

Emission at far-infrared wavelengths makes up a significant fraction of the total light detected from galaxies over the age of Universe. Herschel provides an opportunity for studying galaxies at the peak wavelength of their emission. Our aim is to provide a benchmark for models of galaxy population evolution and to test pre-existing models of galaxies. With the Herschel Multi-tiered Extra-galactic survey, HerMES, we have observed a number of fields of different areas and sensitivity using the SPIRE instrument on Herschel. We have determined the number counts of galaxies down to ~20 mJy. Our constraints from directly counting galaxies are consistent with, though more precise than, estimates from the BLAST fluctuation analysis. We have found a steep rise in the Euclidean normalised counts at <100 mJy. We have directly resolved 15% of the infrared extra-galactic background at the wavelength near where it peaks., Comment: Astronomy and Astrophysics, Herschel Special Issue, in press as a Letter; 5 pages

We set out to determine the ratio, q(IR), of rest-frame 8-1000um flux, S(IR), to monochromatic radio flux, S(1.4GHz), for galaxies selected at far-IR and radio wavelengths, to search for signs that the ratio evolves with redshift, luminosity or dust temperature, and to identify any far-IR-bright outliers - useful laboratories for exploring why the far-IR/radio correlation is generally so tight when the prevailing theory suggests variations are almost inevitable. We use flux-limited 250-um and 1.4-GHz samples, obtained in GOODS-N using Herschel (HerMES; PEP) and the VLA. We determine bolometric IR output using ten bands spanning 24-1250um, exploiting data from PACS and SPIRE, as well as Spitzer, SCUBA, AzTEC and MAMBO. We also explore the properties of an L(IR)-matched sample, designed to reveal evolution of q(IR) with z, spanning log L(IR) = 11-12 L(sun) and z=0-2, by stacking into the radio and far-IR images. For 1.4-GHz-selected galaxies, we see tentative evidence of a break in the flux ratio, q(IR), at L(1.4GHz) ~ 10^22.7 W/Hz, where AGN are starting to dominate the radio power density, and of weaker correlations with z and T(d). From our 250-um-selected sample we identify a small number of far-IR-bright outliers, and see trends of q(IR) with L(1.4GHz), L(IR), T(d) and z, noting that some of these are inter-related. For our L(IR)-matched sample, there is no evidence that q(IR) changes significantly as we move back into the epoch of galaxy formation: we find q(IR) goes as (1+z)^gamma, where gamma = -0.04 +/- 0.03 at z=0-2; however, discounting the least reliable data at z<0.5 we find gamma = -0.26 +/- 0.07, modest evolution which may be related to the radio background seen by ARCADE2, perhaps driven by <10uJy radio activity amongst ordinary star-forming galaxies at z>1., Comment: A&A Herschel Special Issue, in press as a Letter. 5 pages

We use the Spitzer Wide-area InfraRed Extragalactic Legacy Survey (SWIRE) to explore the specific star-formation activity of galaxies and their evolution near the peak of the cosmic far-infrared (FIR) background at 70 and 160um. We use a stacking analysis to determine the mean FIR properties of well defined subsets of galaxies at flux levels well below the FIR catalogue detection limits of SWIRE and other Spitzer surveys. We tabulate the contribution of different subsets of galaxies to the FIR background at 70um and 160um. These long wavelengths provide a good constraint on the bolometric, obscured emission. The large area provides good constraints at low z and in finer redshift bins than previous work. At all redshifts we find that the specific FIR Luminosity (sLFIR) decreases with increasing mass, following a trend L_FIR/M* propto M_* ^beta with beta =-0.38\pm0.14. This is a more continuous change than expected from the {Delucia2007} semi-analytic model suggesting modifications to the feedback prescriptions. We see an increase in the sLFIR by about a factor of ~100 from 0

We present Spitzer 70 and 160 micron observations of the COSMOS Spitzer survey (S-COSMOS). The data processing techniques are discussed for the publicly released products consisting of images and source catalogs. We present accurate 70 and 160 micron source counts of the COSMOS field and find reasonable agreement with measurements in other fields and with model predictions. The previously reported counts for GOODS-North and the extragalactic First Look Survey are updated with the latest calibration, and counts are measured based on the large area SWIRE survey to constrain the bright source counts. We measure an extragalactic confusion noise level of sigma_c = 9.4+/-3.3 mJy (q=5) for the MIPS 160-micron band based on the deep S-COSMOS data and report an updated confusion noise level of sigma_c = 0.35+/-0.15 mJy (q=5) for the MIPS 70-micron band., Comment: Accepted AJ, 15 Aug. 2009. Data available at http://spider.ipac.caltech.edu/staff/frayer/mycosmos/ until released by IRSA

We investigate the clustering of galaxies selected in the 3.6 micron band of the Spitzer Wide-area Infrared Extragalactic (SWIRE) legacy survey. The angular two-point correlation function is calculated for eleven samples with flux limits of S_3.6 > 4-400 mujy, over an 8 square degree field. The angular clustering strength is measured at >5-sigma significance at all flux limits, with amplitudes of A=(0.49-29)\times10^{-3} at one degree, for a power-law model, A\theta^{-0.8}. We estimate the redshift distributions of the samples using phenomological models, simulations and photometric redshifts, and so derive the spatial correlation lengths. We compare our results with the GalICS (Galaxies In Cosmological Simulations) models of galaxy evolution and with parameterized models of clustering evolution. The GalICS simulations are consistent with our angular correlation functions, but fail to match the spatial clustering inferred from the phenomological models or the photometric redshifts. We find that the uncertainties in the redshift distributions of our samples dominate the statistical errors in our estimates of the spatial clustering. At low redshifts (median z<0.5) the comoving correlation length is approximately constant, r_0=6.1\pm0.5h^{-1} Mpc, and then decreases with increasing redshift to a value of 2.9\pm0.3h^{-1} Mpc for the faintest sample, for which the median redshift is z=1. We suggest that this trend can be attributed to a decrease in the average galaxy and halo mass in the fainter flux-limited samples, corresponding to changes in the relative numbers of early- and late-type galaxies. However, we cannot rule out strong evolution of the correlation length over 0.5

We present the stellar calibrator sample and the conversion from instrumental to physical units for the 24 micron channel of the Multiband Imaging Photometer for Spitzer (MIPS). The primary calibrators are A stars, and the calibration factor based on those stars is 4.54*10^{-2} MJy sr^{-1} (DN/s)^{-1}, with a nominal uncertainty of 2%. We discuss the data-reduction procedures required to attain this accuracy; without these procdures, the calibration factor obtained using the automated pipeline at the Spitzer Science Center is 1.6% +/- 0.6% lower. We extend this work to predict 24 micron flux densities for a sample of 238 stars which covers a larger range of flux densities and spectral types. We present a total of 348 measurements of 141 stars at 24 micron. This sample covers a factor of ~460 in 24 micron flux density, from 8.6 mJy up to 4.0 Jy. We show that the calibration is linear over that range with respect to target flux and background level. The calibration is based on observations made using 3-second exposures; a preliminary analysis shows that the calibration factor may be 1% and 2% lower for 10- and 30-second exposures, respectively. We also demonstrate that the calibration is very stable: over the course of the mission, repeated measurements of our routine calibrator, HD 159330, show a root-mean-square scatter of only 0.4%. Finally, we show that the point spread function (PSF) is well measured and allows us to calibrate extended sources accurately; Infrared Astronomy Satellite (IRAS) and MIPS measurements of a sample of nearby galaxies are identical within the uncertainties., Comment: 33 pages, 10 figures, to be published in PASP

Our goal is to probe the populations of obscured and unobscured AGN investigating their optical-IR and X-ray properties as a function of X-ray flux, luminosity and redshift within a hard X-ray selected sample of 136 X-ray sources in the XMM Medium Deep Survey (XMDS) with wide multiwavelength coverage. The XMDS area is covered with optical photometry from the VVDS and CFHTLS surveys and infrared Spitzer data. Based on the X-ray luminosity and X-ray to optical ratio, 132 sources are likely AGN, of which 122 have unambiguous optical - IR identification. The observed optical and IR spectral energy distributions of sources are fitted with AGN/galaxy templates in order to classify them and compute photometric redshifts. 70% of the AGN are fitted by a type 2 AGN or a star forming galaxy template and are grouped together in a single class of ``optically obscured'' AGN. They have ``red'' optical colors and generally show significant X-ray absorption from X-ray spectra or hardness ratios (N$_H > 10^{22}$ cm$^{-2}$). Sources with SEDs typical of type 1 AGN have ``blue'' optical colors and exhibit X-ray absorption in about 30% of cases. We performed a stacking analysis for obscured and type 1 AGN. The stacked X-ray spectrum of obscured AGN is flatter than that of type 1 AGN and has an average spectral slope of Gamma = 1.6. The subsample of objects fitted by a galaxy template has an even harder stacked spectrum, with Gamma = 1.2 - 1.3. The obscured fraction is larger at lower fluxes, lower redshifts and lower luminosities. X-ray absorption is less common than ``optical'' obscuration and its incidence is nearly constant with redshift and luminosity. This implies that X-ray absorption is not necessarily related to optical obscuration., Comment: 33 pages, 21 figures, accepted for publication in A&A

We present the spectral energy distributions (SEDs) of a hard X-ray selected sample. The sample contains 136 sources with F(2-10 keV)>10^-14 erg/cm^2/s and 132 are AGNs. The sources are detected in a 1 square degree area of the XMM-Newton-Medium Deep Survey where optical data from the VVDS, CFHTLS surveys, and infrared data from the SWIRE survey are available. Based on a SED fitting technique we derive photometric redshifts with sigma(1+z)=0.11 and 6% of outliers and identify AGN signatures in 83% of the objects. This fraction is higher than derived when a spectroscopic classification is available. The remaining 17+9-6% of AGNs shows star-forming galaxy SEDs (SF class). The sources with AGN signatures are divided in two classes, AGN1 (33+6-1%) and AGN2 (50+6-11). The AGN1 and AGN2 classes include sources whose SEDs are fitted by type 1 and type 2 AGN templates, respectively. On average, AGN1s show soft X-ray spectra, consistent with being unabsorbed, while AGN2s and SFs show hard X-ray spectra, consistent with being absorbed. The analysis of the average SEDs as a function of X-ray luminosity shows a reddening of the IR SEDs, consistent with a decreasing contribution from the host galaxy at higher luminosities. The AGNs in the SF classes are likely obscured in the mid-infrared, as suggested by their low L(3-20micron)/Lcorr(0.5-10 keV) ratios. We confirm the previously found correlation for AGNs between the radio luminosity and the X-ray and the mid-infrared luminosities. The X-ray-radio correlation can be used to identify heavily absorbed AGNs. However, the estimated radio fluxes for the missing AGN population responsible for the bulk of the background at E>10 keV are too faint to be detected even in the deepest current radio surveys., Comment: Accepted for publication in ApJ

We construct rest-frame luminosity functions at 3.6, 4.5, 5.8, 8 and 24 microns over the redshift range 0