We present the [CII]158$\mu$m line luminosity functions (LFs) at $z\sim 4-6$ using the ALMA observations of 118 sources, which are selected to have UV luminosity $M_{1500A} < -20.2$ and optical spectroscopic redshifts in COSMOS and ECDF-S. Of the 118 targets, 75 have significant [C II] detections and 43 are upper limits. With this UV-selected sample, we set stringent lower bounds to the volume density of [C II] emitters at $z\sim4-6$. Our derived LFs are consistent with the $z\sim0$ [C II] LF at $\leq 10^9L_\odot$, but show a hint of an excess of [C II] emitters at the luminous end ($>10^9L_\odot$) compared to $z\sim0$. This excess is consistent with a potential population of [C II]-bright UV-faint sources that are largely excluded from the ALPINE target selection. To further constrain such sources, in a companion analysis, Loiacono et al. (2020) estimate the [C II] LF based on serendipitous [C II] sources in the ALPINE maps. Their resulting single data point is consistent with the excess we observe at high-L. Combined, these studies suggest possible redshift evolution in the [C II] LFs at high luminosity between $z\sim 5$ and $0$. This conclusion is further supported by the inferred [C II] LFs from the far-IR and CO LFs at $z\sim4-6$. Integrating our estimate on the best-fit [C II] LF, we set a constraint on molecular gas mass density at $z\sim4-6$, with $\rho_{mol} \sim (2-7)\times10^7M_\odot$ Mpc$^{-3}$. This is broadly consistent with previous studies, although the uncertainties span more than an order of magnitude. Finally, we find that available model predictions tend to underestimate the number densities of [C II] emitters at $z\sim5$.

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

We investigate the stellar and dust properties of massive (log$(M_*/M_\odot) \ge 10.5$) and dusty ($A_V \ge 1$) galaxies at $1 \le z \le 4$ by modeling their spectral energy distributions (SEDs) obtained from the combination of UltraVISTA DR3 photometry and \textit{Herschel} PACS-SPIRE data using MAGPHYS. Although the rest-frame U-V vs V-J (UVJ) diagram traces well the star-formation rates (SFR) and dust obscuration (A$_V$) out to $z \sim 3$, $\sim$15-20\% of the sample surprisingly resides in the quiescent region of the UVJ diagram, while $\sim50$\% at $3<z<4$ fall in the unobscured star-forming region. The median SED of massive dusty galaxies exhibits weaker MIR and UV emission, and redder UV slopes with increasing cosmic time. The IR emission for our sample has a significant contribution ($>20\%$) from dust heated by evolved stellar populations rather than star formation, demonstrating the need for panchromatic SED modeling. The local relation between dust mass and SFR is followed only by a sub-sample with cooler dust temperatures, while warmer objects have reduced dust masses at a given SFR. Most star-forming galaxies in our sample do not follow local IRX-$\beta$ relations, though IRX does strongly correlate with A$_V$. Our sample follows local relations, albeit with large scatter, between ISM diagnostics and sSFR. We show that FIR-detected sources represent the extreme of a continuous population of dusty galaxies rather than a fundamentally different population. Finally, using commonly adopted relations to derive SFRs from the combination of the rest-frame UV and the observed 24$\mu$m is found to overestimate the SFR by a factor of 3-5 for the galaxies in our sample.

The bulk of the star-formation rate density peak at cosmic noon was obscured by dust. How accurately we can assess the role of dust obscured star-formation is affected by inherent biases in our empirical methods -- both those that rely on direct dust emission and those that rely on the inferred dust attenuation of starlight. We use a library of hydrodynamic simulations with radiative transfer to explore these biases. We find that for IR luminous galaxies that are in rapidly quenching systems (e.g. post-coalescence) standard luminosity-to-SFR relations can strongly overestimate the true SFRs. We propose using the $L_{IR}/L_{1.6}$ color to both help identify such systems and provide more accurate SFRs. Conversely, we find that the diagnostic UVJ plot misidentifies a subset of dusty star-forming galaxies. This is due to variability in the effective attenuation curves including being much grayer in the optical-to-near-IR regime than the Calzetti starburst law. This is in agreement with recent observations of IR-selected galaxies at cosmic noon. Our results support the view that we need a panchromatic approach from the rest-frame UV through the IR and SED modeling that includes realistic SFHs and allows for variable attenuation curves if we want to fully account for dust obscured star-formation across the epochs of greatest galaxy build-up.

Galaxies grow their supermassive black holes in concert with their stars, although the relationship between these major galactic components is poorly understood. Observations of the cosmic growth of stars and black holes in galaxies suffer from disjoint samples and the strong effects of dust attenuation. The thermal infrared holds incredible potential for simultaneously measuring both the star formation and black hole accretion rates in large samples of galaxies covering a wide range of physical conditions. Spitzer demonstrated this potential at low redshift, and by observing some of the most luminous galaxies at z~2. JWST will apply these methods to normal galaxies at these epochs, but will not be able to generate large spectroscopic samples or access the thermal infrared at high-redshift. An order of magnitude gap in our wavelength coverage will persist between JWST and ALMA. A large, cold infrared telescope can fill this gap to determine when (in cosmic time), and where (within the cosmic web), stars and black holes co-evolve, by measuring these processes simultaneously in statistically complete and unbiased samples of galaxies to z>8. A next-generation radio interferometer will have the resolution and sensitivity to measure star-formation and nuclear accretion in even the dustiest galaxies. Together, the thermal infrared and radio can uniquely determine how stars and supermassive blackholes co-evolve in galaxies over cosmic time.

We have observed the environments of a population of 33 heavily dust obscured, ultra-luminous, high-redshift galaxies, selected using WISE and NVSS at $z>$1.3 with the Infra-Red Array Camera on the $Spitzer$ Space Telescope over $\rm5.12\,'\times5.12\,'$ fields. Colour selections are used to quantify any potential overdensities of companion galaxies in these fields. We find no significant excess of galaxies with the standard colour selection for IRAC colours of $\rm[3.6]-[4.5]>-0.1$ consistent with galaxies at $z>$1.3 across the whole fields with respect to wide-area $Spitzer$ comparison fields, but there is a $\rm>2\sigma$ statistical excess within $\rm0.25\,'$ of the central radio-WISE galaxy. Using a colour selection of $\rm[3.6]-[4.5]>0.4$, 0.5 magnitudes redder than the standard method of selecting galaxies at $z>$1.3, we find a significant overdensity, in which $\rm76\%$ ($\rm33\%$) of the 33 fields have a surface density greater than the $\rm3\sigma$ ($\rm5\sigma$) level. There is a statistical excess of these redder galaxies within $\rm0.5\,'$, rising to a central peak $\rm\sim2$--4 times the average density. This implies that these galaxies are statistically linked to the radio-WISE selected galaxy, indicating similar structures to those traced by red galaxies around radio-loud AGN.

We present images taken using the Gemini South Adaptive Optics Imager (GSAOI) with the Gemini Multiconjugate Adaptive Optics System (GeMS) in three 2 arcmin$^2$ fields in the Spitzer Extragalactic Representative Volume Survey. These GeMS/GSAOI observations are among the first $\approx 0.1^{''}$ resolution data in the near-infrared spanning extragalactic fields exceeding $1.5^{\prime}$ in size. We use these data to estimate galaxy sizes, obtaining results similar to those from studies with the Hubble Space Telescope, though we find a higher fraction of compact star forming galaxies at $z>2$. To disentangle the star-forming galaxies from active galactic nuclei (AGN), we use multiwavelength data from surveys in the optical and infrared, including far-infrared data from Herschel, as well as new radio continuum data from the Australia Telescope Compact Array and Very Large Array. We identify ultraluminous infrared galaxies (ULIRGs) at $z \sim 1-3$, which consist of a combination of pure starburst galaxies and Active Galactic Nuclei (AGN)/starburst composites. The ULIRGs show signs of recent merger activity, such as highly disturbed morphologies and include a rare candidate triple AGN. We find that AGN tend to reside in hosts with smaller scale sizes than purely star-forming galaxies of similar infrared luminosity. Our observations demonstrate the potential for MCAO to complement the deeper galaxy surveys to be made with the James Webb Space Telescope.

The bulk of the stellar growth over cosmic time is dominated by IR luminous galaxies at cosmic noon (z=1-2), many of which harbor a hidden active galactic nucleus (AGN). We use state of the art infrared color diagnostics, combining Spitzer and Herschel observations, to separate dust-obscured AGN from dusty star forming galaxies (SFGs) in the CANDELS and COSMOS surveys. We calculate 24 micron counts of SFGs, AGN/star forming "Composites", and AGN. AGN and Composites dominate the counts above 0.8 mJy at 24 micron, and Composites form at least 25% of an IR sample even to faint detection limits. We develop methods to use the Mid-Infrared Instrument (MIRI) on JWST to identify dust-obscured AGN and Composite galaxies from z~1-2. With the sensitivity and spacing of MIRI filters, we will detect >4 times as many AGN hosts than with Spitzer/IRAC criteria. Any star formation rates based on the 7.7 micron PAH feature (likely to be applied to MIRI photometry) must be corrected for the contribution of the AGN, or the SFR will be overestimated by ~35% for cases where the AGN provides half the IR luminosity and ~50% when the AGN accounts for 90% of the luminosity. Finally, we demonstrate that our MIRI color technique can select AGN with an Eddington ratio of $\lambda_{\rm Edd}\sim0.01$ and will identify AGN hosts with a higher sSFR than X-ray techniques alone. JWST/MIRI will enable critical steps forward in identifying and understanding dust-obscured AGN and the link to their host galaxies.

At $z=1-3$, the formation of new stars is dominated by dusty galaxies whose far-IR emission indicates they contain colder dust than local galaxies of a similar luminosity. We explore the reasons for the evolving IR emission of similar galaxies over cosmic time using: 1) Local galaxies from GOALS $(L_{\rm IR}=10^{11}-10^{12}\,L_\odot)$; 2) Galaxies at $z\sim0.1-0.5$ from the 5MUSES ($L_{\rm IR}=10^{10}-10^{12}\,L_\odot$); 3) IR luminous galaxies spanning $z=0.5-3$ from GOODS and Spitzer xFLS ($L_{\rm IR}>10^{11}\,L_\odot$). All samples have Spitzer mid-IR spectra, and Herschel and ground-based submillimeter imaging covering the full IR spectral energy distribution, allowing us to robustly measure $L_{\rm IR}^{\rm\scriptscriptstyle SF}$, $T_{\rm dust}$, and $M_{\rm dust}$ for every galaxy. Despite similar infrared luminosities, $z>0.5$ dusty star forming galaxies have a factor of 5 higher dust masses and 5K colder temperatures. The increase in dust mass is linked with an increase in the gas fractions with redshift, and we do not observe a similar increase in stellar mass or star formation efficiency. $L_{160}^{\rm\scriptscriptstyle SF}/L_{70}^{\rm\scriptscriptstyle SF}$, a proxy for $T_{\rm dust}$, is strongly correlated with $L_{\rm IR}^{\rm\scriptscriptstyle SF}/M_{\rm dust}$ independently of redshift. We measure merger classification and galaxy size for a subsample, and there is no obvious correlation between these parameters and $L_{\rm IR}^{\rm \scriptscriptstyle SF}/M_{\rm dust}$ or $L_{160}^{\rm\scriptscriptstyle SF}/L_{70}^{\rm\scriptscriptstyle SF}$. In dusty star forming galaxies, the change in $L_{\rm IR}^{\rm\scriptscriptstyle SF}/M_{\rm dust}$ can fully account for the observed colder dust temperatures, suggesting that any change in the spatial extent of the interstellar medium is a second order effect.

We apply The Tractor image modeling code to improve upon existing multi-band photometry for the Spitzer Extragalactic Representative Volume Survey (SERVS). SERVS consists of post-cryogenic Spitzer observations at 3.6 and 4.5 micron over five well-studied deep fields spanning 18 square degrees. In concert with data from ground-based near-infrared (NIR) and optical surveys, SERVS aims to provide a census of the properties of massive galaxies out to z ~ 5. To accomplish this, we are using The Tractor to perform "forced photometry." This technique employs prior measurements of source positions and surface brightness profiles from a high-resolution fiducial band from the VISTA Deep Extragalactic Observations (VIDEO) survey to model and fit the fluxes at lower-resolution bands. We discuss our implementation of The Tractor over a square degree test region within the XMM-LSS field with deep imaging in 12 NIR/optical bands. Our new multi-band source catalogs offer a number of advantages over traditional position-matched catalogs, including 1) consistent source cross-identification between bands, 2) de-blending of sources that are clearly resolved in the fiducial band but blended in the lower-resolution SERVS data, 3) a higher source detection fraction in each band, 4) a larger number of candidate galaxies in the redshift range 5 < z < 6, and 5) a statistically significant improvement in the photometric redshift accuracy as evidenced by the significant decrease in the fraction of outliers compared to spectroscopic redshifts. Thus, forced photometry using The Tractor offers a means of improving the accuracy of multi-band extragalactic surveys designed for galaxy evolution studies. We will extend our application of this technique to the full SERVS footprint in the future.

The assessment of the relationship between radio continuum luminosity and star formation rate (SFR) is of crucial importance to make reliable predictions for the forthcoming ultra-deep radio surveys and to allow a full exploitation of their results to measure the cosmic star formation history. We have addressed this issue by matching recent accurate determinations of the SFR function up to high redshifts with literature estimates of the 1.4 GHz luminosity functions of star forming galaxies (SFGs). This was done considering two options, proposed in the literature, for the relationship between the synchrotron emission ($L_{\rm synch}$), that dominates at 1.4 GHz, and the SFR: a linear relation with a decline of the $L_{\rm synch}$/SFR ratio at low luminosities or a mildly non-linear relation at all luminosities. In both cases we get good agreement with the observed radio luminosity functions but, in the non-linear case, the deviation from linearity must be small. The luminosity function data are consistent with a moderate increase of the $L_{\rm synch}$/SFR ratio with increasing redshift, indicated by other data sets, although a constant ratio cannot be ruled out. A stronger indication of such increase is provided by recent deep 1.4 GHz counts, down to $\mu$Jy levels. This is in contradiction with models predicting a decrease of that ratio due to inverse Compton cooling of relativistic electrons at high redshifts. Synchrotron losses appear to dominate up to $z\simeq 5$. We have also updated the Massardi et al. (2010) evolutionary model for radio loud AGNs.

The James Webb Space Telescope's Medium Resolution Spectrometer (MRS), will offer nearly 2 orders of magnitude improvement in sensitivity and >3X improvement in spectral resolution over our previous space-based mid-IR spectrometer, the Spitzer IRS. In this paper, we make predictions for spectroscopic pointed observations and serendipitous detections with the MRS. Specifically, pointed observations of Herschel sources require only a few minutes on source integration for detections of several star-forming and active galactic nucleus lines, out to z$=$3 and beyond. But the same data will also include tens of serendipitous 0$\lesssim$z$\lesssim$4 galaxies per field with infrared luminosities ranging $\sim10^6-10^{13}$L$_{\odot}$. In particular, for the first time and for free we will be able to explore the $L_{IR}<10^{9}L_{\odot}$ regime out to $z\sim3$. We estimate that with $\sim$100 such fields, statistics of these detections will be sufficient to constrain the evolution of the low-$L$ end of the infrared luminosity function, and hence the star formation rate function. The above conclusions hold for a wide range in potential low-$L$ end of the IR luminosity function, and accounting for the PAH deficit in low-$L$, low-metallicity galaxies.

A key question in extragalactic studies is the determination of the relative roles of stars and AGN in powering dusty galaxies at $z\sim$1-3 where the bulk of star-formation and AGN activity took place. In Paper I, we present a sample of $336$ 24$\mu$m-selected (Ultra)Luminous Infrared Galaxies, (U)LIRGs, at $z \sim 0.3$-$2.8$, where we focus on determining the AGN contribution to the IR luminosity. Here, we use hydrodynamic simulations with dust radiative transfer of isolated and merging galaxies, to investigate how well the simulations reproduce our empirical IR AGN fraction estimates and determine how IR AGN fractions relate to the UV-mm AGN fraction. We find that: 1) IR AGN fraction estimates based on simulations are in qualitative agreement with the empirical values when host reprocessing of the AGN light is considered; 2) for star-forming galaxy-AGN composites our empirical methods may be underestimating the role of AGN, as our simulations imply $>$50% AGN fractions, $\sim$3$\times$ higher than previous estimates; 3) 6% of our empirically classified "SFG" have AGN fractions $\gtrsim$50%. While this is a small percentage of SFGs, if confirmed, would imply the true number density of AGN may be underestimated; 4) this comparison depends on the adopted AGN template -- those that neglect the contribution of warm dust lower the empirical fractions by up to 2$\times$; and 5) the IR AGN fraction is only a good proxy for the intrinsic UV-mm AGN fraction when the extinction is high ($A_V\gtrsim 1$ or up to and including coalescence in a merger).

We characterize infrared spectral energy distributions of 343 (Ultra) Luminous Infrared Galaxies from $z=0.3-2.8$. We diagnose the presence of an AGN by decomposing individual Spitzer mid-IR spectroscopy into emission from star-formation and an AGN-powered continuum; we classify sources as star-forming galaxies (SFGs), AGN, or composites. Composites comprise 30% of our sample and are prevalent at faint and bright $S_{24}$, making them an important source of IR AGN emission. We combine spectroscopy with multiwavelength photometry, including Herschel imaging, to create three libraries of publicly available templates (2-1000 $\mu$m). We fit the far-IR emission using a two temperature modified blackbody to measure cold and warm dust temperatures ($T_c$ and $T_w$). We find that $T_c$ does not depend on mid-IR classification, while $T_w$ shows a notable increase as the AGN grows more luminous. We measure a quadratic relationship between mid-IR AGN emission and total AGN contribution to $L_{\rm IR}$. AGN, composites, and SFGs separate in $S_8/S_{3.6}$ and $S_{250}/S_{24}$, providing a useful diagnostic for estimating relative amounts of these sources. We estimate that >40% of IR selected samples host an AGN, even at faint selection thresholds ($S_{24}$>100 $\mu$Jy). Our decomposition technique and color diagnostics are relevant given upcoming observations with the James Webb Space Telescope.

We combine a deep 0.5~deg$^2$, 1.4~GHz deep radio survey in the Lockman Hole with infrared and optical data in the same field, including the SERVS and UKIDSS near-infrared surveys, to make the largest study to date of the host galaxies of radio sources with typical radio flux densities $\sim 50 \;\mu$Jy. 87% (1274/1467) of radio sources have identifications in SERVS to $AB\approx 23.1$ at 3.6 or 4.5$\mu$m, and 9% are blended with bright objects (mostly stars), leaving only 4% (59 objects) which are too faint to confidently identify in the near-infrared. We are able to estimate photometric redshifts for 68% of the radio sources. We use mid-infrared diagnostics to show that the source population consists of a mixture of star forming galaxies, rapidly accreting (cold mode) AGN and low accretion rate, hot mode AGN, with neither AGN nor starforming galaxies clearly dominating. We see the breakdown in the $K-z$ relation in faint radio source samples, and show that it is due to radio source populations becoming dominated by sources with radio luminosities $\sim 10^{23}\;{\rm WHz^{-1}}$. At these luminosities, both the star forming galaxies and the cold mode AGN have hosts with stellar luminosities about a factor of two lower than those of hot mode AGN, which continue to reside in only the most massive hosts. We show that out to at least $z\sim 2$, galaxies with stellar masses $>10^{11.5}\, M_{\odot}$ have a radio-loud fraction up to $\sim 30$%. This is consistent with there being a sufficient number of radio sources that radio-mode feedback could play a role in galaxy evolution.

We study the environments of 49 WISE/NVSS-selected dusty, hyper-luminous, z~2 quasars using the Atacama Large Millimeter/Sub-millimeter Array (ALMA) 345GHz images. We find that 17 of the 49 WISE/NVSS sources show additional sub-mm galaxies within the ALMA primary beam, probing scales within ~150 kpc. We find a total of 23 additional sub-mm sources, four of which in the field of a single WISE/NVSS source. The measured 870 um source counts are ~10 times expectations for unbiased regions, suggesting such hyper-luminous dusty quasars are excellent at probing high-density peaks.

We present luminosity functions derived from a spectroscopic survey of AGN selected from Spitzer Space Telescope imaging surveys. Selection in the mid-infrared is significantly less affected by dust obscuration. We can thus compare the luminosity functions of the obscured and unobscured AGN in a more reliable fashion than by using optical or X-ray data alone. We find that the AGN luminosity function can be well described by a broken power-law model in which the break luminosity decreases with redshift. At high redshifts ($z>1.6$), we find significantly more AGN at a given bolometric luminosity than found by either optical quasar surveys or hard X-ray surveys. The fraction of obscured AGN decreases rapidly with increasing AGN luminosity, but, at least at high redshifts, appears to remain at $\approx 50$\% even at bolometric luminosities $\sim 10^{14}L_{\odot}$. The data support a picture in which the obscured and unobscured populations evolve differently, with some evidence that high luminosity obscured quasars peak in space density at a higher redshift than their unobscured counterparts. The amount of accretion energy in the Universe estimated from this work suggests that AGN contribute about 12\% to the total radiation intensity of the Universe, and a high radiative accretion efficiency $\approx 0.18^{+0.12}_{-0.07}$ is required to match current estimates of the local mass density in black holes.

We present a new galaxy survey simulation package, which combines the power of Markov Chain Monte Carlo (MCMC) sampling with a robust and adaptable model of galaxy evolution. The aim of this code is to aid in the characterization and study of new and existing galaxy surveys. In this paper we briefly describe the MCMC implementation and the survey simulation methodology and associated tools. A test case of this full suite was to constrain the evolution of the IR Luminosity Function (LF) based on the HerMES (Herschel SPIRE) survey of the Spitzer First Look Survey field. The initial results are consistent with previous studies, but our more general approach should be of wider benefit to the community.

We present the results of a program of optical and near-infrared spectroscopic follow-up of candidate Active Galactic Nuclei (AGN) selected in the mid-infrared. This survey selects both normal and obscured AGN closely matched in luminosity across a wide range, from Seyfert galaxies with bolometric luminosities L_bol~10^10L_sun, to highly luminous quasars (L_bol~10^14L_sun), and with redshifts from 0-4.3. Samples of candidate AGN were selected through mid-infrared color cuts at several different 24 micron flux density limits to ensure a range of luminosities at a given redshift. The survey consists of 786 candidate AGN and quasars, of which 672 have spectroscopic redshifts and classifications. Of these, 137 (20%) are type-1 AGN with blue continua, 294 (44%) are type-2 objects with extinctions A_V>~5 towards their AGN, 96 (14%) are AGN with lower extinctions (A_V~1) and 145 (22%) have redshifts, but no clear signs of AGN activity in their spectra. 50% of the survey objects have L_bol >10^12L_sun, in the quasar regime. We present composite spectra for type-2 quasars and for objects with no signs of AGN activity in their spectra. We also discuss the mid-infrared - emission-line luminosity correlation and present the results of cross-correlations with serendipitous X-ray and radio sources. The results show that: (1) obscured objects dominate the overall AGN population, (2) there exist mid-infrared selected AGN candidates which lack AGN signatures in their optical spectra, but have AGN-like X-ray or radio counterparts, and (3) X-ray and optical classifications of obscured and unobscured AGN often differ.

We analyze the link between active galactic nuclei (AGN) and mid-infrared flux using dust radiative transfer calculations of starbursts realized in hydrodynamical simulations. Focusing on the effects of galaxy dust, we evaluate diagnostics commonly used to disentangle AGN and star formation in ultraluminous infrared galaxies (ULIRGs). We examine these quantities as a function of time, viewing angle, dust model, AGN spectrum, and AGN strength in merger simulations representing two possible extremes of the ULIRG population: one is a typical gas-rich merger at z ~ 0, and the other is characteristic of extremely obscured starbursts at z ~ 2 to 4. This highly obscured burst begins star-formation-dominated with significant PAH emission, and ends with a ~10^9 yr period of red near-IR colors. At coalescence, when the AGN is most luminous, dust obscures the near-infrared AGN signature, reduces the relative emission from polycyclic aromatic hydrocarbons (PAHs), and enhances the 9.7 micron absorption by silicate grains. Although generally consistent with previous interpretations, our results imply none of these indicators can unambiguously estimate the AGN luminosity fraction in all cases. Motivated by the simulations, we show that a combination of the extinction feature at 9.7 micron, the PAH strength, and a near-infrared slope can simultaneously constrain the AGN fraction and dust grain distribution for a wide range of obscuration. We find that this indicator, accessible to the James Webb Space Telescope, may estimate the AGN power as tightly as the hard X-ray flux alone, thereby providing a valuable future cross-check and constraint for large samples of distant ULIRGs.

(abridged for arXiv) We make use of the Planck all-sky survey to derive number counts and spectral indices of extragalactic sources -- infrared and radio sources -- from the Planck Early Catalogue (ERCSC) at 100 to 857GHz. Our sample contains, after the 80% completeness cut, between 122 and 452 and sources, with flux densities above 0.3 and 1.9Jy at 100 and 857GHz, over about 31 to 40% of the sky. Using Planck HFI, all the sources have been classified as either dust-dominated or synchrotron-dominated on the basis of their spectral energy distributions (SED). Our sample is thus complete, flux-limited and color-selected to differentiate between the two populations. We find an approximately equal number of synchrotron and dusty sources between 217 and 353GHz; at 353GHz or higher (or 217GHz and lower) frequencies, the number is dominated by dusty (synchrotron) sources, as expected. For most of the sources, the spectral indices are also derived. We provide for the first time counts of bright sources from 353 to 857GHz and the contributions from dusty and synchrotron sources at all HFI frequencies in the key spectral range where these spectra are crossing. The observed counts are in the Euclidean regime. The number counts are compared to previously published data (earlier Planck, Herschel, BLAST, SCUBA, LABOCA, SPT, and ACT) and models taking into account both radio or infrared galaxies. We derive the multi-frequency Euclidean level and compare it to WMAP, Spitzer and IRAS results. The submillimetre number counts are not well reproduced by current evolution models of dusty galaxies, whereas the millimetre part appears reasonably well fitted by the most recent model for synchrotron-dominated sources. Finally we provide estimates of the local luminosity density of dusty galaxies, providing the first such measurements at 545 and 857GHz.

We report multiple epoch VLA/JVLA observations of 89 northern hemisphere sources, most with 37\,GHz flux density > 1 Jy, observed at 4.8, 8.5, 33.5, and 43.3 GHz. The high frequency selection leads to a predominantly flat spectrum sample, with 85% of our sources being in the Planck Early Release Compact Source Catalog (ERCSC). These observations allow us to: 1) validate Planck's 30 and 44 GHz flux density scale, 2) extend the radio SEDs of Planck sources to lower frequencies allowing for the full 5-857GHz regime to be studied, and 3) characterize the variability of these sources. At 30 GHz and 44 GHz, the JVLA and Planck flux densities agree to within 3%. On timescales of less than two months the median variability of our sources is 2%. On timescales of about a year the median variability increases to 14%. Using the WMAP 7-year data, the 30 GHz median variability on a 1-6 years timescale is 16%.

The Planck Early Release Compact Source Catalog (ERCSC) has offered the first opportunity to accurately determine the luminosity function of dusty galaxies in the very local Universe (i.e. distances <~ 100 Mpc), at several (sub-)millimetre wavelengths, using blindly selected samples of low redshift sources, unaffected by cosmological evolution. This project, however, requires careful consideration of a variety of issues including the choice of the appropriate flux density measurement, the separation of dusty galaxies from radio sources and from Galactic sources, the correction for the CO emission, the effect of density inhomogeneities, and more. We present estimates of the local luminosity functions at 857 GHz (350 microns), 545 GHz (550 microns) and 353 GHz (850 microns) extending across the characteristic luminosity L_star, and a preliminary estimate over a limited luminosity range at 217 GHz (1382 microns). At 850 microns and for luminosities L >~ L_star our results agree with previous estimates, derived from the SCUBA Local Universe Galaxy Survey (SLUGS), but are higher than the latter at L <~ L_star. We also find good agreement with estimates at 350 and 500 microns based on preliminary Herschel survey data.

We present the mid-infrared spectra of seven of the most powerful radio-galaxies known to exist at 1.5 < z < 2.6. The radio emission of these sources is dominated by the AGN with 500 MHz luminosities in the range 10^27.8 - 10^29.1 W/Hz. The AGN signature is clearly evident in the mid-infrared spectra, however, we also detect polycyclic aromatic hydrocarbons emission, indicative of prodigious star formation at a rate of up to ~1000 Msun/yr. Interestingly, we observe no significant correlation between AGN power and star formation in the host galaxy. We also find most of these radio galaxies to have weak 9.7 um silicate absorption features (tau_{9.7 um} < 0.8) which implies that their mid-IR obscuration is predominantly due to the dusty torus that surrounds the central engine, rather than the host galaxy. The tori are likely to have an inhomogeneous distribution with the obscuring structure consisting of individual clouds. We estimate that these radio galaxies have already formed the bulk of their stellar mass and appear to lie at a stage in their evolution where the obscured AGN dominates the energy output of the system but star formation is also prevalent.

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.

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.

We have obtained Spitzer IRS spectra and MIPS 24, 70, and 160 micron photometry for a volume-limited sample of 22 SDSS-selected Low-ionization Broad Absorption Line QSOs (LoBALs) at 0.5 < z < 0.6. By comparing their mid-IR spectral properties and far-IR SEDs with those of a control sample of 35 non-LoBALs matched in M_i, we investigate the differences between the two populations in terms of their infrared emission and star formation activity. Twenty five percent of the LoBALs show PAH features and 45% have weak 9.7 micron silicate dust emission. We model the SEDs and decouple the AGN and starburst contributions to the far-infrared luminosity in LoBALs and in non-LoBALs. Their median total, starburst, and AGN infrared luminosities are comparable. Twenty percent (but no more than 60%) of the LoBALs and 26% of the non-LoBALs are ultra-luminous infrared galaxies (ULIRGs; L_IR>10^12*L_sun). We estimate star formation rates (SFRs) corrected for the AGN contribution to the FIR flux and find that LoBALs have comparable levels of star formation activity to non-LoBALs when considering the entire samples. However, the SFRs of the IR-luminous LoBALs are 80% higher than those of their counterparts in the control sample. The median contribution of star formation to the total far-infrared flux in LoBALs and in non-LoBALs is estimated to be 40-50%, in agreement with previous results for PG QSOs. Overall, our results show that there is no strong evidence from the mid- and far-IR properties that LoBALs are drawn from a different parent population than non-LoBALs.

We present simultaneous Planck, Swift, Fermi, and ground-based data for 105 blazars belonging to three samples with flux limits in the soft X-ray, hard X-ray, and gamma-ray bands. Our unique data set has allowed us to demonstrate that the selection method strongly influences the results, producing biases that cannot be ignored. Almost all the BL Lac objects have been detected by Fermi-LAT, whereas ~40% of the flat-spectrum radio quasars (FSRQs) in the radio, soft X-ray, and hard X-ray selected samples are still below the gamma-ray detection limit even after integrating 27 months of Fermi-LAT data. The radio to sub-mm spectral slope of blazars is quite flat up to ~70GHz, above which it steepens to <\alpha>~-0.65. BL Lacs have significantly flatter spectra than FSRQs at higher frequencies. The distribution of the rest-frame synchrotron peak frequency (\nupS) in the SED of FSRQs is the same in all the blazar samples with <\nupS>=10^13.1 Hz, while the mean inverse-Compton peak frequency, <\nupIC>, ranges from 10^21 to 10^22 Hz. The distributions of \nupS and of \nupIC of BL Lacs are much broader and are shifted to higher energies than those of FSRQs and strongly depend on the selection method. The Compton dominance of blazars ranges from ~0.2 to ~100, with only FSRQs reaching values >3. Its distribution is broad and depends strongly on the selection method, with gamma-ray selected blazars peaking at ~7 or more, and radio-selected blazars at values ~1, thus implying that the assumption that the blazar power is dominated by high-energy emission is a selection effect. Simple SSC models cannot explain the SEDs of most of the gamma-ray detected blazars in all samples. The SED of the blazars that were not detected by Fermi-LAT may instead be consistent with SSC emission. Our data challenge the correlation between bolometric luminosity and \nupS predicted by the blazar sequence.

We present preliminary results on fitting of SEDs to 142 z>1 quasars selected in the mid-infrared. Our quasar selection finds objects ranging in extinction from highly obscured, type-2 quasars, through more lightly reddened type-1 quasars and normal type-1s. We find a weak tendency for the objects with the highest far-infrared emission to be obscured quasars, but no bulk systematic offset between the far-infrared properties of dusty and normal quasars as might be expected in the most naive evolutionary schemes. The hosts of the type-2 quasars have stellar masses comparable to those of radio galaxies at similar redshifts. Many of the type-1s, and possibly one of the type-2s require a very hot dust component in addition to the normal torus emission.

We present optical through radio observations of the host galaxies of two dust obscured, luminous quasars selected in the mid-infrared, at z=2.62 and z=2.99, including a search for CO emission. Our limits on the CO luminosities are consistent with these objects having masses of molecular gas <~10^10 solar masses, several times less than those of luminous submillimeter-detected galaxies (SMGs) at comparable redshifts. Their near-infrared spectral energy distributions, however, imply that these galaxies have high stellar masses (~10^11-12 solar masses). The relatively small reservoirs of molecular gas and low dust masses are consistent with them being relatively mature systems at high-z.

(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.

The data reported in Planck's Early Release Compact Source Catalogue (ERCSC) are exploited to measure the number counts (dN/dS) of extragalactic radio sources at 30, 44, 70, 100, 143 and 217 GHz. Due to the full-sky nature of the catalogue, this measurement extends to the rarest and brightest sources in the sky. At lower frequencies (30, 44, and 70 GHz) our counts are in very good agreement with estimates based on WMAP data, being somewhat deeper at 30 and 70 GHz, and somewhat shallower at 44 GHz. Planck's source counts at 143 and 217 GHz join smoothly with the fainter ones provided by the SPT and ACT surveys over small fractions of the sky. An analysis of source spectra, exploiting Planck's uniquely broad spectral coverage, finds clear evidence of a steepening of the mean spectral index above about 70 GHz. This implies that, at these frequencies, the contamination of the CMB power spectrum by radio sources below the detection limit is significantly lower than previously estimated.

We present flux densities and polarization percentages of 159 radio galaxies based on nearly simultaneous VLA observations at four frequencies from 5 to 43GHz. This sample is selected from the high-frequency Australia Telescope 20GHz (AT20G) survey and consists of all sources with S_{20GHz}>40mJy in an equatorial field of the ACT survey. For a subset of 25 of these sources we used the GBT to obtain 90GHz data. The goals of this program are: 1) a characterization of the spectra, polarization and variability of high frequency-selected radio sources; 2) extrapolating from the few GHz regime to the ~150GHz regime of the ACT survey, allowing for more accurate removal of the radio source signal in our particular field; and 3) providing a data set that will allow more accurate modeling of the high-frequency radio source contamination in current and future SZ and CMB experiments. We find that, as expected, this sample consists of flatter spectrum and more compact or point-like sources than low frequency-selected samples. In the K-band, variability is typically <~20%, although there are exceptions. The higher frequency data is well suited to the detection of extreme Giga-Hertz Peak spectrum Sources (GPS). The inclusion of the 43GHz data causes the relative fraction of inverted spectrum sources to go down and of peaked spectrum sources to go up when compared with the AT20G survey results. The trend largely continues with the inclusion of the 90GHz data, although ~10% of the sources with GBT data show a spectral upturn from 43GHz to 90GHz. The measured polarization fractions are typically <5%, although in some cases they are measured to be up to ~20%. For the ~40% of the sample with detected polarized flux in all four bands, the polarization fractions typically increase with frequency. This trend is stronger for steeper spectrum sources as well as for the lower flux density sources.

We present new results on the physical nature of infrared-luminous sources at 0.5<z<2.8 as revealed by HST/NICMOS imaging and IRS mid-infrared spectroscopy. Our sample consists of 134 galaxies selected at 24\mum with a flux of S(24\mum) > 0.9 mJy. We find many (~60%) of our sources to possess an important bulge and/or central point source component, most of which reveal additional underlying structures after subtraction of a best-fit sersic (or sersic+PSF) profile. Based on visual inspection of the NIC2 images and their residuals, we estimate that ~80% of all our sources are mergers. We calculate lower and upper limits on the merger fraction to be 62% and 91% respectively. At z < 1.5, we observe objects in early (pre-coalescence) merging stages to be mostly disk and star formation dominated, while we find mergers to be mainly bulge-dominated and AGN-starburst composites during coalescence and then AGN-dominated in late stages. This is analogous to what is observed in local ULIRGs. At z \geq 1.5, we find a dramatic rise in the number of objects in pre-coalescence phases of merging, despite an increase in the preponderance of AGN signatures in their mid-IR spectra and luminosities above 10^12.5 L_sun. We further find the majority of mergers at those redshifts to retain a disk-dominated profile during coalescence. We conclude that, albeit still driven by mergers, these high-z ULIRGs are substantially different in nature from their local counterparts and speculate that this is likely due to their higher gas content. Finally, we observe obscured ({\tau}_{9.7\mum} > 3.36) quasars to live in faint and compact hosts and show that these are likely high-redshift analogs of local dense-core mergers. We find late-stage mergers to show predominantly unobscured AGN spectra, but do not observe other morphological classes to occupy any one specific region in the Spoon diagram. [abridged]

Abridged: We present a detailed investigation of the morphological distribution and level of star formation and dust obscuration in the nearby tidally distorted galaxy NGC2442. Spitzer images in the IR at 3.6, 4.5, 5.8, 8.0um, and 24um and GALEX images at 1500\AA{} and 2300\AA{} allow us to resolve the galaxy on scales between 240-600pc. We supplement these with archival data in the B, J, H, and K bands. We use the 8um, 24um and FUV (1500\AA) emission to study the star formation rate (SFR). We find that globally, these tracers of star formation give a range of results of ~6-11\msun/yr, with the dust-corrected FUV giving the highest value of SFR. We can reconcile the UV and IR-based estimates by adopting a steeper UV extinction curve that lies in between the starburst (Calzetti) and SMC extinction curves. However, the regions of highest SFR intensity along the spiral arms are consistent with a starburst-like extinction. Overall, the level of star-formation we find is higher than previously published for this galaxy, by about a factor of two, which, contrary to previous conclusions, implies that the interaction that caused the distorted morphology of NGC2442 likely also triggered increased levels of star-formation activity. Outside of the spiral arms, we discover what appears to be a superbubble, ~1.7kpc across in the IRAC images. Significant H{\alpha}, UV and IR emission in the area also suggest vigorous ongoing star-formation. A known, recent supernova (SN1999ga) is located at the edge of this superbubble. Although speculative at this stage, this area suggests a large star-forming region with a morphology shaped by generations of supernovae. Lastly, we discover an 8um (PAH) circumnuclear ring with an ~0.8kpc radius. The H{\alpha} emission is largely concentrated inside that ring and shows a vague spiral structure in the rest of the galaxy.

We present ultra-deep mid-IR spectra of 48 infrared-luminous galaxies in the GOODS-South field obtained with the InfraRed Spectrograph (IRS) on the Spitzer Space Telescope. These galaxies are selected among faint infrared sources (0.14 - 0.5 mJy at 24 um) in two redshift bins (0.76-1.05 and 1.75-2.4) to sample the major contributors to the cosmic infrared background at the most active epochs. We estimate redshifts for 92% of the sample using PAH and Si absorption features. Only few of these galaxies (5% at z~1 and 12% at z~2) have their total infrared luminosity dominated by emission from AGN. The averaged mid-IR spectra of the z~1 LIRGs and of the z~2 ULIRGs are very similar to the averaged spectrum of local starbursts and HII-like ULIRGs, respectively. We find that 6.2um PAH equivalent widths reach a plateau of ~1 um for L(24 mu) < 1E11 L(sun). At higher luminosities, EW (6.2 mu) anti-correlates with L(24 um). Intriguingly, high-z ULIRGs and SMG lie above the local EW (6.2 um) - L(24 um) relationship suggesting that, at a given luminosity, high-z ULIRGs have AGN contributions to their dust emission lower than those of local counterparts. A quantitative analysis of their morphology shows that most of the luminous IR galaxies have morphologies similar to those of IR-quiet galaxies at the same redshift. All z~2 ULIRGs of our sample are IR-excess BzK galaxies and most of them have L(FIR)/L(1600A) ratios higher than those of starburst galaxies at a given UV slope. The ``IR excess'' (Daddi et al. 2007) is mostly due to strong 7.7 um PAH emission and under-estimation of UV dust extinction. On the basis of the AGN-powered L (6 um) continuum measured directly from the mid-IR spectra, we estimate an average intrinsic X-ray AGN luminosity of L(2-10 keV) = (0.1 +/- 0.6) 1E43 erg/s, a value substantially lower than the prediction by Daddi et al. (2007).

We present CO observations of 9 ULIRGs at z~2 with S(24\mu m)>1mJy, previously confirmed with the mid-IR spectra in the Spitzer First Look Survey. All targets are required to have accurate redshifts from Keck/GEMINI near-IR spectra. Using the Plateau de Bure millimeter-wave Interferometer (PdBI) at IRAM, we detect CO J(3-2) [7 objects] or J(2-1) [1 object] line emission from 8 sources with integrated intensities Ic ~(5-9)sigma. The CO detected sources have a variety of mid-IR spectra, including strong PAH, deep silicate absorption and power-law continuum, implying that these molecular gas rich objects at z~2 could be either starbursts or dust obscured AGNs. The measured line luminosity L'[CO] is (1.28-3.77)e+10[K km/s pc^2]. The averaged molecular gas mass M(H2) is 1.7e+10Msun, assuming CO-to-H2 conversion factor of 0.8Msun/[K km/s pc^2]. Three sources (33%) -- MIPS506, MIPS16144 & MIPS8342 -- have double peak velocity profiles. The CO double peaks in MIPS506 and MIPS16144 show spatial separations of 45kpc and 10.9kpc, allowing the estimates of the dynamical masses of 3.2e+11*sin^(-2)(i)Msun and 5.4e+11*sin^{-2}(i)Msun respectively. The implied gas fraction, M(gas)/M(dyn), is 3% and 4%, assuming an average inclination angle. Finally, the analysis of the HST/NIC2 images, mid-IR spectra and IR SED revealed that most of our sources are mergers, containing dust obscured AGNs dominating the luminosities at (3-6)um. Together, these results provide some evidence suggesting SMGs, bright 24um z~2 ULIRGs and QSOs could represent three different stages of a single evolutionary sequence, however, a complete physical model would require much more data, especially high spatial resolution spectroscopy.

We report X-ray constraints for 20 of 52 high-z ULIRGs identified in the Spitzer xFLS to constrain their obscuration. Notably, decomposition of Spitzer-IRS spectra for the 52 objects already indicates that most are weak-PAH ULIRGs dominated by hot-dust continua, characteristic of AGN. Given their redshifts, they have AGN bolometric luminosities of ~1e45-1e47 erg/s comparable to powerful QSOs. This, coupled with their high IR-to-optical ratios and often significant silicate absorption, strongly argues in favor of these mid-IR objects being heavily obscured QSOs. At X-ray energies, we marginally detect two ULIRGs, while the rest have only upper limits. Using the IRS-derived 5.8um AGN continuum luminosity as a proxy for the expected X-ray luminosities, we find that all of the observed sources must individually be highly obscured, while X-ray stacking limits on the undetected sources suggest that the majority, if not all, are likely to be at least mildly Compton-thick (NH>1e24 cm-2). With a space density of ~1.4e-7 Mpc-3 at z~2, such objects imply a lower limit on the obscured AGN fraction (i.e., the ratio of AGN above and below NH=1e22 cm-2) of >~1.7:1 even among luminous QSOs. Our findings, which are based on extensive multi-wavelength constraints including Spitzer IRS spectra, should aid in the interpretation of similar objects from larger or deeper mid-IR surveys, where considerable uncertainty about the source properties remains and comparable follow-up is not yet feasible.

We model the optical to far-infrared SEDs of a sample of six type-1 and six type-2 quasars selected in the mid-infrared. The objects in our sample are matched in mid-IR luminosity and selected based on their Spitzer IRAC colors. We obtained new targeted Spitzer IRS and MIPS observations and used archival photometry to examine the optical to far-IR SEDs. We investigate whether the observed differences between samples are consistent with orientation-based unification schemes. The type-1 objects show significant emission at 3 micron. They do not show strong PAH emission and have less far-IR emission on average when compared to the type-2 objects. The SEDs of the type-2 objects show a wide assortment of silicate features, ranging from weak emission to deep silicate absorption. Some also show strong PAH features. In comparison, silicate is only seen in emission in the type-1 objects. This is consistent with some of the type-2s being reddened by a foreground screen of cooler dust, perhaps in the host galaxy itself. We investigate the AGN contribution to the far-IR emission and find it to be significant. We also estimate the star formation rate for each of the objects by integrating the modeled far-IR flux and compare this with the SFR found from PAH emission. We find the type-2 quasars have a higher average SFR than the type-1 quasars based on both methods, though this could be due to differences in bolometric luminosities of the objects. While we find pronounced differences between the two types of objects, none of them are inconsistent with orientation-based unification schemes.

We present the first detections of the 3um water ice and 3.4um amorphous hydrocarbon (HAC) absorption features in z~2 ULIRGs. These are based on deep rest-frame 2-8um Spitzer IRS spectra of 11 sources selected for their appreciable silicate absorption. The HAC-to-silicate ratio for our z~2 sources is typically higher by a factor of 2-5 than that observed in the Milky Way. This HAC `excess' suggests compact nuclei with steep temperature gradients as opposed to predominantly host obscuration. Beside the above molecular absorption features, we detect the 3.3um PAH emission feature in one of our sources with three more individual spectra showing evidence for it. Stacking analysis suggests that water ice, hydrocarbons, and PAH are likely present in the bulk of this sample even when not individually detected. The most unexpected result of our study is the lack of clear detections of the 4.67um CO gas absorption feature. Only three of the sources show tentative signs of this feature and at significantly lower levels than has been observed in local ULIRGs. Overall, we find that the closest local analogs to our sources, in terms of 3-4um color, HAC-to-silicate and ice-to-silicate ratios, as well as low PAH equivalent widths are sources dominated by deeply obscured nuclei. Such sources form only a small fraction of ULIRGs locally and are commonly believed to be dominated by buried AGN. Our sample suggests that, in absolute number, such buried AGN are at least an order of magnitude more common at z~2 than today. The presence of PAH suggests that significant levels of star-formation are present even if the obscured AGN typically dominate the power budget.

We present a catalog of mid-infrared (MIR) spectra of 150 infrared (IR) luminous galaxies in the Spitzer extragalactic first look survey obtained with IRS on board Spitzer. The sample is selected to be brighter than ~0.9 mJy at 24 micron and it has a z distribution in the range [0.3,3.5] with a peak at z=1. It primarily comprises ultraluminous IR galaxies at z>1 and luminous IR galaxies at z<1, as estimated from their monochromatic 14 micron luminosities. The number of sources with spectra that are dominated by an active galactic nucleus (AGN) continuum is 49, while 39 sources have strong, star-formation related features. For this classification, we used the equivalent width (EW) of the 11.3 micron polycyclic aromatic hydrocarbon (PAH) feature. Several intermediate/high z starbursts have higher PAH EW than local ULIRGs. An increase in the AGN activity is observed with increasing z and luminosity, based on the decreasing EW of PAHs and the increasing [NeIII]/[NeII] ratio. Spectral stacking leads to the detection of the 3.3 micron PAH, the H2 0-0 S(1) and S(3) lines, and the [NeV] line. We observe differences in the flux ratios of PAHs in the stacked spectra of IR-luminous galaxies with z or luminosity, which are not due to extinction effects. When placing the observed galaxies on IR color-color diagrams, we find that the wedge defining AGN comprises most sources with continuum-dominated spectra, but also contains many starbursts. The comparison of the 11.3 micron PAH EW and the H-band effective radius, measured from HST data, indicates that sources with EW>2 micron, are typically more extended than ~3 kpc. However, there is no strong correlation between the MIR spectral type and the near-IR extent of the sources. [Abridged].

We present preliminary results from the highest available signal-to-noise rest-frame 2-8um spectra of z~2 ULIRGs. Our 10 targets are selected for their deep silicate absorption features based on previous shallower IRS spectra. The goal of this follow-up program is: 1) allow for a more accurate analysis of inner/hot dust continuum, 2) detecting the 3.3um and 6.2um PAH features, and 3) detecting molecular absorption features such as due to water ice and hydrocarbons (HACs). We find that the 3.4um HAC absorption feature is observed in four sources, while the 3.05um water ice feature is observed in three of the sources. The HAC detectability is higher and ice detectability lower than expected from local ULIRGs, but consistent with a more AGN-dominated sample such as this one. Where ice is detected, the ice-to-silicate ratio is somewhat lower than many local ULIRGs implying on average thinner ice mantles. One source shows the, to our knowledge, highest redshift reported detection of the 3.3um PAH feature (along with a previously detected 6.2um feature). The strength of the 3.3um feature is as expected for a starburst-dominated ULIRG.

We present Hubble Space Telescope NICMOS H-band imaging of 33 Ultraluminous Infrared Galaxies (ULIRGs) at z~2 that were selected from the 24 micron catalog of the Spitzer Extragalactic First Look Survey. The images reveal that at least 17 of the 33 objects are associated with interactions. Up to one fifth of the sources in our sample could be minor mergers whereas only 2 systems are merging binaries with luminosity ratio <=3:1, which is characteristic of local ULIRGs. The rest-frame optical luminosities of the sources are of the order 10^10-10^11 L_sun and their effective radii range from 1.4 to 4.9 kpc. The most compact sources are either those with a strong active nucleus continuum or those with a heavy obscuration in the mid-infrared regime, as determined from Spitzer Infra-Red Spectrograph data. The luminosity of the 7.7 micron feature produced by polycyclic aromatic hydrocarbon molecules varies significantly among compact systems whereas it is typically large for extended systems. A bulge-to-disk decomposition performed for the 6 brightest (m_H<20) sources in our sample indicates that they are best fit by disk-like profiles with small or negligible bulges, unlike the bulge-dominated remnants of local ULIRGs. Our results provide evidence that the interactions associated with ultraluminous infrared activity at z~2 can differ from those at z~0.

We present the first mid-infrared Spitzer/Infrared Spectrograph (IRS) observations of powerful radio galaxies at z>2. These radio galaxies, 4C +23.56 (z=2.48) and 6C J1908+7220 (z=3.53), both show strong mid-infrared continua, but with 6C J1908+7220 also showing strong PAH emission at rest-frame 6.2 and 7.7um. In 4C+23.56 we see no obvious PAH features above the continuum. The PAH emission in 6C J1908+7220 is the amongst the most distant observed to date and implies that there is a large instantaneous star formation rate (SFR). This is consistent with the strong detection of 6C J1908+7220 at far-IR and sub-mm wavelengths, indicative of large amounts of cold dust, ~10^9Msun. Powerful radio galaxies at lower redshifts tend to have weak or undetectable PAH features and typically have lower far-IR luminosities. In addition, 4C 23.56 shows moderate silicate absorption as seen in less luminous radio galaxies, indicating tau_{9.7um}=0.3+/-0.05. This feature is shifted out of the observed wavelength range for 6C J1908+7220. The correlation of strong PAH features with large amounts of cold dust, despite the presence of a powerful AGN, is in agreement with other recent results and implies that star formation at high redshift is, in some cases at least, associated with powerful, obscured AGN.

We present the far-IR, millimeter, and radio photometry as well as optical and near-IR spectroscopy of a sample of 48 z~1-3 Spitzer-selected ULIRGs with IRS mid-IR spectra. Our goals are to compute their bolometric emission, and to determine both the presence and relative strength of their AGN and starburst components. We find that strong-PAH sources tend to have higher 160um and 1.2mm fluxes than weak-PAH sources. The depth of the 9.7um silicate feature does not affect MAMBO detectability. We fit the far-IR SEDs of our sample and find an average <L_{IR}>~7x10^{12}Lsun for our z>1.5 sources. Spectral decomposition suggests that strong-PAH sources typically have ~20-30% AGN fractions. Weak-PAH sources by contrast tend to have >~70% AGN fractions, with a few sources having comparable contributions of AGN and starbursts. The optical line diagnostics support the presence of AGN in the bulk of the weak-PAH sources. With one exception, our sources are narrow-line sources, show no obvious correspondence between the optical extinction and the silicate feature depth, and, in two cases, show evidence for outflows. Radio AGN are present in both strong-PAH and weak-PAH sources. This is supported by our sample's far-IR-to-radio ratios (q) being consistently below the average value of 2.34 for local star-forming galaxies. We use survival analysis to include the lower-limits given by the radio-undetected sources, arriving at <q>=2.07+/-0.01 for our z>1.5 sample. In total, radio and, where available, optical line diagnostics support the presence of AGN in 57% of the z>1.5 sources, independent of IR-based diagnostics. For higher-z sources, the AGN luminosities alone are estimated to be >10^{12}Lsun, which, supported by the [OIII] luminosities, implies that the bulk of our sources host obscured quasars.

We present mid-infrared observations of 18 sources from a sample of 21 z~2 radio-intermediate obscured quasars. The mid-infrared spectra of the sources are continuum dominated, and 12 sources show deep silicate absorption with tau_9.7~1-2. Combining mid-infrared and optical spectra, we achieve 86% spectroscopic completeness which allows us to confirm that most (63 +14 -22 %) z~2 radio-intermediate quasars are obscured. The new spectra also prove that many high-redshift type-2 quasars do not show any rest-frame ultraviolet emission lines. From the 18 individual mid-infrared spectra, we classify most of the sources into three subsamples: those with hints of the 7.7 and 6.2 micron polyaromatic hydrocarbons (3/18 sources show PAHs, subsample A), those with an excess of emission around 8 micron but no hint of the 6.2 micron PAH (7/18 cases, subsample B) and pure-continuum sources with no visible excess (4/18 sources, subsample C). The remaining 4/18 sources have spectra that are featureless or too noisy for any features to be visible. In subsample A, averaging the spectra leads to a statistical detection of both 6.2 and 7.7 micron PAHs over the continuum, with the strength of the 7.7 micron PAH comparable to that of submillimetre-selected galaxies (SMGs) at similar redshifts. These sources are in a phase of coeval growth of a supermassive black hole and a host galaxy.

We present Hubble Space Telescope images, and spectral energy distributions from optical to infrared wavelengths for a sample of six 0.3<z<0.8 type-2 quasars selected in the mid-infrared using data from the Spitzer Space Telescope. All the host galaxies show some signs of disturbance. Most seem to possess dusty, star-forming disks. The disk inclination, estimated from the axial ratio of the hosts, correlates with the depth of the silicate feature in the mid-infrared spectra, implying that at least some of the reddening towards the AGN arises in the host galaxy. The star formation rates in these objects, as inferred from the strengths of the PAH features and far-infrared continuum, range from 3-90 Msun/yr, but are mostly much larger than those inferred from the [OII]3727 emission line luminosity, due to obscuration. Taken together with studies of type-2 quasar hosts from samples selected in the optical and X-ray, this is consistent with previous suggestions that two types of extinction processes operate within the type-2 quasar population, namely a component due to the dusty torus in the immediate environment of the AGN, and a more extended component due to a dusty, star forming disk.

Recent Spitzer observations have revealed a substantial population of z~2 ULIRGs with deep silicate absorption (\tau_{9.7}>1). This paper reports a 20cm radio study of such a sample to elucidate their physical nature. We discover that a substantial fraction (40%) of deep silicate absorption ULIRGs at z~2 are moderately radio-loud with L_{1.4GHz}=10^{25}--10^{26}WHz^{-1}. This is in strong contrast with z<1 radio galaxies and radio-loud quasars where none of the sources with available IRS spectra have \tau_{9.7}>1. In addition, we observe radio jets in two of our sources where one has a double lobe structure ~200kpc in extent, and another shows a one-sided jet extending ~90kpc from the nucleus. The likely high inclination of the latter, coupled with its deep silicate absorption, implies the mid-IR obscuration does not share an axis with the radio jets. These sources are highly obscured quasars, observed in the transition stage after the birth of the radio source, but before feedback effects dispel the ISM and halt the black hole accretion and starburst activity.

We present deep polarimetric observations at 1420 MHz of the European Large Area ISO Survey North 1 region (ELAIS N1) as part of the Dominion Radio Astrophysical Observatory Planck Deep Fields project. By combining closely spaced aperture synthesis fields, we image a region of 7.43 square degrees to a maximum sensitivity in Stokes Q and U of 78 microJy/beam, and detect 786 compact sources in Stokes I. Of these, 83 exhibit polarized emission. We find that the differential source counts (log N - log p) for polarized sources are nearly constant down to p > 500 microJy, and that these faint polarized radio sources are more highly polarized than the strong source population. The median fractional polarization is (4.8 +/- 0.7)% for polarized sources with Stokes I flux density between 1 and 30 mJy; approximately three times larger than sources with I > 100 mJy. The majority of the polarized sources have been identified with galaxies in the Spitzer Wide Area Infrared Extragalactic Survey (SWIRE) image of ELAIS N1. Most of the galaxies occupy regions in the IRAC 5.8/3.6 micron vs. 8.0/4.5 micron color-color diagram associated with dusty AGNs, or with ellipticals with an aging stellar population. A few host galaxies have colors that suggests significant PAH emission in the near-infrared. A small fraction, 12%, of the polarized sources are not detected in the SWIRE data. None of the polarized sources in our sample appears to be associated with an actively star-forming galaxy.

We present mid-IR spectral decomposition of a sample of 48 Spitzer-selected ULIRGs spanning z~1-3 and likely L_IR~10^12-10^13Lsun. Our study aims at quantifying the star-formation and AGN processes in these sources which recent results suggest have evolved strongly between the observed epoch and today. To do this, we study the mid-IR contribution of PAH emission, continuum, and extinction. About 3/4 of our sample are continuum- (i.e. AGN) dominated sources, but ~60% of these show PAH emission, suggesting the presence of star-formation activity. These sources have redder mid-IR colors than typical optically-selected quasars. About 25% of our sample have strong PAH emission, but none are likely to be pure starbursts as reflected in their relatively high 5um hot dust continua. However, their steep 30um-to-14um slopes suggest that star-formation might dominate the total infrared luminosity. Six of our z~2 sources have EW6.2>~0.3um and L_14um>~10^12Lsun (implying L_IR>~10^13Lsun). At these luminosities, such high EW6.2 ULIRGs do not exist in the local Universe. We find a median optical depth at 9.7um of <tau_9.7>=1.4. This is consistent with local IRAS-selected ULIRGs, but differs from early results on SCUBA-selected z~2 ULIRGs. Similar to local ULIRGs about 25% of our sample show extreme obscuration (tau_9.7>~3) suggesting buried nuclei. In general, we find that our sources are similar to local ULIRGs, but are an order of magnitude more luminous. It is not clear whether our z~2 ULIRGs are simply scaled-up versions of local ULIRGs, or subject to fundamentally different physical processes.

We present the mid-infrared (MIR) spectra obtained with the Spitzer InfraRed Spectrograph (IRS) for a sample of 52 sources, selected as infrared luminous, z>1 candidates in the Extragalactic First Look Survey (XFLS). The sample selection criteria are f(24um) > 0.9mJy, nu fnu(24um)/nu fnu(8um) > 3.16 and nu fnu(24um)/nu fnu(0.7um) > 10. Of the 52 spectra, 47 (90%) produced measurable redshifts based solely on the mid-IR spectral features, with the majority (35/47=74%) at 1.5<z<3.2. Keck spectroscopy of a sub-sample (17/47) agrees with the mid-IR redshift measurements. The observed spectra fall crudely into three categories -- (1) 33% (17/52) have strong PAH emission, and are probably powered by star formation with total IR luminosity roughly a factor of 5 higher than the local starburst ULIRGs. (2) 33% (17/52) have only deep silicate absorption at 9.8um, indicative of deeply embedded dusty systems. The energetic nature of the heating sources in these systems can not be determined by these data alone. (3) The remainder 34% are mid-IR continuum dominated systems with either weak PAH emission and/or silicate absorption. This third of the sample are probably AGNs. From the silicate absorption feature, we estimate that roughly two-thirds of the sample have optical depth tau(9.8um) > 1. Their $L_{1600\AA}$ and $L_{\rm IR}$ suggest that our sample is among the most luminous and most dust enshrouded systems of its epoch . Our study has revealed a significant population of dust enshrouded galaxies at z~2, whose enormous energy output, comparable to that of quasars, is generated by AGN as well as starburst. This IR luminous population has very little overlap with sub-mm and UV-selected populations (Abridged).

We present optical spectroscopy of a sample of 77 luminous AGN and quasars selected on the basis of their mid-infrared colors. Our objects are selected from the Spitzer Extragalactic First Look Survey and SWIRE XMM-LSS fields, with a typical 24mu flux density of 5mJy. The median redshift is 0.6, with a range of ~0.05-4. Only 33% (25/77) of these objects are normal type-1 quasars, with no obscuration. 44% (34/77) are type-2 objects, with high-ionization, narrow emission lines, 14% (11/77) are dust-reddened type-1 quasars, showing broad lines but a dust-reddened or unusually weak quasar continuum. 9% (7/77) show no sign of an AGN in the optical spectrum, having either starburst spectra or spectra which could be of either starburst or LINER types. These latter objects are analogous to the X-ray detected population of AGN with weak or non-existent optical AGN emission (the ``XBONGs''). 21 of our objects from the SWIRE field fall within moderately-deep XMM exposures. All the unobscured quasars, and about half the obscured quasars are detected in these exposures. This sample, when taken together with other samples of Spitzer selected AGN and quasars, and results from X-ray studies, confirms that obscured AGN dominate the AGN and quasar number counts of all rapidly-accreting supermassive black hole systems, at least for z~<4. This implies a high radiative efficiency for the black hole accretion process.

We present SEDs, Spitzer colours, and IR luminosities for 850 micron selected galaxies in the GOODS-N field. Using the deep Spitzer Legacy images and new data and reductions of the VLA-HDF radio data, we find statistically secure counterparts for 60 per cent (21/35) of our submm sample, and identify tentative counterparts for another 12 objects. This is the largest sample of submm galaxies with statistically secure counterparts detected in the radio and with Spitzer. We find that in most cases the 850 micron emission is dominated by a single 24 micron source. A composite rest-frame SED shows that the submm sources peak at longer wavelengths than those of local ULIRGs of the same luminosity and therefore appear to be cooler. The SEDs of submm galaxies are also different from those of their high redshift neighbours, the near-IR selected BzK galaxies, whose mid-IR to radio SEDs are more like those of local ULIRGs. Using 24 micron. 850 micron and 1.4 GHz observations, we fit templates that span the mid-IR through radio to derived the integrated IR luminosity of the submm galaxies. By themselves, 24 micron and radio fluxes are able to predict LIR reasonably well because they are relatively insensitive to temperature. However, the submm flux by itself consistently overpredicts LIR when using spectral templates which obey the local ULIRG temperature-luminosity relation. The shorter Spitzer wavelengths sample the stellar bump at the redshifts of the submm sources, and we find that the Spitzer photometry alone provides a model independent estimate of the redshift. Using X-ray and mid-IR data, only 5 per cent of our secure counterparts show strong evidence for an active galactic nucleus dominating the IR luminosity.

We present a preliminary investigation of the spectral energy distributions (SEDs), and star-formation properties of a sample of Mid-IR selected Quasars. The mid-infrared SEDs of our objects are consistent with that expected from clumpy torus models. At longer wavelengths, the radio to infrared ratios of several objects are consistent with those of star-forming galaxies.

Galaxies selected at 170um by the ISO FIRBACK survey represent the brightest \~10% of the Cosmic Infrared Background. Examining their nature in detail is therefore crucial for constraining models of galaxy evolution. Here we combine Spitzer archival data with previous near-IR, far-IR, and sub-mm observations of a representative sample of 22 FIRBACK galaxies spanning three orders of magnitude in infrared luminosity. We fit a flexible, multi-component, empirical SED model of star-forming galaxies designed to model the entire ~1-1000um wavelength range. The fits are performed with a Markov Chain Monte Carlo (MCMC) approach, allowing for meaningful uncertainties to be derived. This approach also highlights degeneracies such as between Td and beta, which we discuss in detail. From these fits and standard relations we derive: L_IR, L_PAH, SFR, tau_V, M_star, M_dust, Td, and beta. We look at a variety of correlations between these and combinations thereof in order to examine the physical nature of these galaxies. Our conclusions are supplemented by morphological examination of the sources, and comparison with local samples. We find the bulk of our sample to be consistent with fairly standard size and mass disk galaxies with somewhat enhanced star-formation relative to local spirals, but likely not bona fide starbursts. A few higher-z LIGs and ULIGs are also present, but contrary to expectation, they are weak mid-IR emitters and overall are consistent with star-formation over an extended cold region rather than concentrated in the nuclear regions. We discuss the implications of this study for understanding populations detected at other wavelengths, such as the bright 850um SCUBA sources or the faint Spitzer 24um sources.

How does the submm galaxy population detected by SCUBA fit into galaxy evolution? How do these rare starbursting systems, which contribute significantly to high redshift star-formation, relate to other galaxy populations? Deep radio observations have been most useful for studying these systems, but still leave a significant fraction of the population unidentified. Now with the deep IRAC and MIPS images from the GOODS Spitzer Legacy program and a re-analysis of the deep radio data, we are able to identify counterparts for a large fraction of SCUBA galaxies in GOODS-N. All of these counterparts are detected by Spitzer. Given the vast multi-wavelength data in this field, we can study the spectral energy distributions (SEDs) of these systems and determine what is fueling their intense infrared luminosities. A rest-frame composite optical-to-radio SED for all spectroscopically identified submm sources shows that the average SCUBA galaxy is consistent with models of ultraluminous starburst galaxies, although cooler than those observed locally. Because of this, the submm flux alone consistently overestimates IR when using spectral templates which obey the local ULIRG temperature-luminosity relation. The wide range of 24/850 micron flux ratios as a function of redshift indicates the presence of strong mid-IR features, to be confirmed with deep IRS spectroscopy. The IRAC colours of the submm systems provide useful redshift constraints, since, at these redshifts, IRAC samples the stellar bump. The Spitzer photometry of this large sample of submm galaxies has allowed us to put constraints on many of the outstanding issues in submm astronomy.

We present preliminary results on a study of the 2--850 micron SEDs of a sample of 30 FIRBACK galaxies selected at 170 micron. These sources are representative of the brightest ~10% of the Cosmic Infrared Background. They are a mixture of mostly local (z<~0.3) starforming galaxies, and a tail of ULIGs that extend up to z~1, and are likely to be a similar population to faint SCUBA sources. We use archival Spitzer IRAC and MIPS data to extend the spectral coverage to the mid-IR regime, resulting in an unprecended (for this redshift range) census of their infrared SEDs. This allows us to study in far greater detail this important population linking the near-IR stellar emission with PAH and thermal dust emission. We do this using a Markov Chain Monte Carlo method, which easily allows for the inclusion of ~6 free parameters, as well as an estimate of parameter uncertainties and correlations.

We present early results from the spectroscopic follow-up of a sample of candidate obscured AGN selected in the mid-infrared from the Spitzer First Look Survey. Our selection allows a direct comparison of the numbers of obscured and unobscured AGN at a given luminosity for the first time, and shows that the ratio of obscured to unobscured AGN at infrared luminosities corresponding to low luminosity quasars is ~1:1 at z~0.5. Most of our optically-faint candidate obscured AGN have the high-ionization, narrow-line spectra expected from type-2 AGN. A composite spectrum shows evidence for Balmer absorption lines, indicating recent star-formation activity in the host galaxies. There is tentative evidence for a decrease in the obscured AGN fraction with increasing AGN luminosity.

We use a simple parameterization of the mid-IR spectra of a wide range of galaxy types in order to predict their distribution in the Infrared Array Camera (IRAC) 3.6um, 4.5um, 5.8um, and 8.0um and Multiband Photometer for Spitzer (MIPS) 24um color-color diagrams. We distinguish three basic spectral types by the energetically dominant component in the 3--12 micron regime: stellar-dominated; PAH-dominated; and continuum-dominated. We have used a Markov chain Monte Carlo (MCMC) approach to arrive at a more systematic and robust representation of the mid-IR spectra of galaxies than more traditional approaches. We find that IRAC color-color plots are well suited to distinguishing the above spectral types, while the addition of 24 micron data allows us to suggest practical 3-color cuts which preferentially select higher redshift sources of specific type. We compare our simulations with the color-color plot obtained by the Spitzer First Look Survey (FLS) and find reasonable agreement. Lastly, we discuss other applications as well as future directions for this work.

Selection of active galactic nuclei (AGN) in the infrared allows the discovery of AGN whose optical emission is extinguished by dust. In this paper, we use the Spitzer Space Telescope First Look Survey (FLS) to assess what fraction of AGN with mid-infrared luminosities comparable to quasars are missed in optical quasar surveys due to dust obscuration. We begin by using the Sloan Digital Sky Survey (SDSS) database to identify 54 quasars within the 4 deg^2 extragalactic FLS. These quasars occupy a distinct region in mid-infrared color space by virtue of their strong, red, continua. This has allowed us to define a mid-infrared color criterion for selecting AGN candidates. About 2000 FLS objects have colors consistent with them being AGN, but most are much fainter in the mid-infrared than the SDSS quasars, which typically have 8 micron flux densities, S(8.0), ~1 mJy. We have investigated the properties of the 43 objects with S(8.0) >= 1 mJy satisfying our AGN color selection. This sample should contain both unobscured quasars, and AGN which are absent from the SDSS survey due to extinction in the optical. After removing 16 known quasars, three probable normal quasars, and eight spurious or confused objects from the initial sample of 43, we are left with 16 objects which are likely to be obscured quasars or luminous Seyfert-2 galaxies. This suggests the numbers of obscured and unobscured AGN are similar in samples selected in the mid-infrared at S(8.0)~1 mJy.

We present results from JCMT sub-mm observations of sources selected from the {\sl ISO} FIRBACK (Far-IR BACKground) survey, along with UKIRT near-IR imaging of a sub-sample. This gives valuable insight into the brightest $\sim$10% of galaxies which contribute to the Cosmic Infrared Background (CIB). We estimate the photometric redshifts and luminosities of these sources by fitting their Spectral Energy Distributions (SEDs). The data appear to show a bimodal galaxy distribution, with normal star-forming galaxies at $z\simeq0$, and a much more luminous population at $z\sim0.4$--0.9. These are similar to the ultraluminous infrared galaxies which are found to evolve rapidly with redshift in other surveys. The detectability threshold of FIRBACK biases the sample away from much higher redshift ($z\stackrel{>}{_{\sim}}1.5$) objects. Nevertheless, the handful of $z\sim0.5$ sources which we identify are likely to be the low-$z$ counterparts of the typically higher-$z$ sources found in blank field sub-mm observations. This sub-sample, being much more nearby than the average SCUBA galaxies, has the virtue of being relatively easy to study in the optical. Hence their detailed investigation could help elucidate the nature of the sub-mm bright galaxies.

Summary and discussion of some projects to use SCUBA to target sources selected at other wavebands, as well as to find new sub-mm galaxies in `blank fields': FIRBACK galaxies; Lyman break galaxies and `the Blob'; HDF flanking fields and the Groth Strip; survey of lensing cluster fields.

The FIRBACK (Far InfraRed BACKground) survey represents the deepest extensive 170 micron images obtained by the ISO satellite. The sources detected comprise about 10% of the Cosmic IR Background (CIB) seen by COBE, and, importantly, were observed at a wavelength near the peak of the CIB. Detailed follow-up of these sources should help pin down the redshifts, Hubble types, and other properties of the galaxies which constitute this background. We have used the Submillimetre Common-User Bolometer Array (SCUBA) instrument on the James Clerk Maxwell Telescope (JCMT) to search for sub-mm emission from a sample of 10 galaxies which are bright at 170 microns and for which we had accurate radio positions. Statistically we detect this sample at both 450 microns (~4 sigma) and 850 microns (~ 7 sigma); individual objects are convincingly detected in four cases at 850 microns, and one case at 450 microns. Fits to simple spectral energy distributions suggest a range of low to moderate redshifts, perhaps z=0 to 1.5 for this FIRBACK sub-sample.