Search

We present an analysis of 8 JWST Mid-Infrared Instrument 5.6 micron images with 5sigma depths of ~0.1uJy. We detect 2854 sources within our combined area of 18.4 sq.arcmin -- a >4x increase in source density over earlier IRAC channel 3 data. We compute the MIRI 5.6um number counts including an analysis of the field-to-field variation. Relative to earlier published MIRI 5.6micron counts, our counts have a more pronounced knee, at roughly 2\,$\mu$Jy. The location and amplitude of the counts at the knee are consistent with the Cowley et al. (2018) model predictions, although these models tend to overpredict the counts below the knee. In areas of overlap, 84% of the MIRI sources have a counterpart in the COSMOS2020 catalog. These MIRI sources have redshifts that are mostly in the z~0.5-2, with a tail out to z~5. They are predominantly moderate to low stellar masses $10^8-10^{10}$M$_{\odot}$) main sequence star-forming galaxies suggesting that with $\approx$2hr exposures, MIRI can reach well below $M^*$ at cosmic noon and reach higher mass systems out to z~5. Nearly 70% of the COSMOS2020 sources in areas of overlap now have a data point at 5.6micron (rest-frame near-IR at cosmic noon) which allows for more accurate stellar population parameter estimates. Finally, we discover 31 MIRI-bright sources not in COSMOS2020. A cross-match with IRAC channel 1 suggests that 10-20% of these are likely lower mass (M$_*\approx10^9$M$_{\odot}$), $z\sim1$ dusty galaxies. The rest (80--90%) are consistent with more massive, but still very dusty galaxies at z>3., Comment: 17 pages, 14 figures, submitted to ApJ

The current and upcoming large data volume galaxy surveys require the use of machine learning techniques to maximize their scientific return. This study explores the use of Self-Organizing Maps (SOMs) to estimate galaxy parameters with a focus on handling cases of missing data and providing realistic probability distribution functions for the parameters. We train a SOM with a simulated mass-limited lightcone assuming a ugrizYJHKs+IRAC dataset, mimicking the Hyper Suprime-Cam (HSC) Deep joint dataset. For parameter estimation, we derive SOM likelihood surfaces considering photometric errors to derive total (statistical and systematic) uncertainties. We explore the effects of missing data including which bands are particular critical to the accuracy of the derived parameters. We demonstrate that the parameter recovery is significantly better when the missing bands are "filled-in" rather than if they are completely omitted. We propose a practical method for such recovery of missing data., Comment: 22 pages, 12 figures, 1 table, accepted for publication in The Astronomical Journal

We present JWST/MIRI/MRS observations of an infrared luminous disk galaxy, FLS1, at z=0.54. With a lookback time of 5 Gyr, FLS1 is chronologically at the midpoint between the peak epoch of star formation and the present day. The MRS data provide maps of the atomic fine structure lines [Ar II]6.99 micron, [Ar III]8.99 micron, [Ne II]12.81 micron, and [Ne III]15.55 micron, polycyclic aromatic hydrocarbon (PAH) features at 3.3 micron, 6.2 micron, and 11.3 micron, and the warm molecular gas indicators H2S(5) and H2S(3); all these emission features are spatially resolved. We find that the PAH emission is more extended along the Northern side of the galaxy when compared to the well-studied star-formation tracer [Ne II]. The H2 rotational lines, which are shock indicators, are strongest and most extended on the Southern side of the galaxy. [Ar II] is the second brightest fine structure line detected in FLS1 and we show that it is a useful kinematic probe which can be detected with JWST out to z=3. Velocity maps of [Ar II] show a rotating disk with signs of turbulence. Our results provide an example of how spatially resolved mid-infrared spectroscopy can allow us to better understand the star formation and ISM conditions in a galaxy halfway back to the peak epoch of galaxy evolution., Comment: 10 pages, 5 figures, 1 table

We present new Spitzer Infrared Array Camera (IRAC) 3.6 and 4.5{\mu}m mosaics of three fields, E-COSMOS, DEEP2-F3, and ELAIS-N1. Our mosaics include both new IRAC observations as well as re-processed archival data in these fields. These fields are part of the HSC-Deep grizy survey and have a wealth of additional ancillary data. The addition of these new IRAC mosaics is critical in allowing for improved photometric redshifts and stellar population parameters at cosmic noon and earlier epochs. The total area mapped by this work is {\sim} 17 deg2 with a mean integration time of {\sim}1200s, providing a median 5{\sigma} depth of 23.7(23.3) at 3.6(4.5){\mu}m in AB. We perform SExtractor photometry both on the combined mosaics as well as the single-epoch mosaics taken {\sim}6 months apart. The resultant IRAC number counts show good agreement with previous studies. In combination with the wealth of existing and upcoming spectro-photometric data in these fields, our IRAC mosaics will enable a wide range of galactic evolution and AGN studies. With that goal in mind, we make the combined IRAC mosaics and coverage maps of these three fields publicly available. counts show good agreement with previous studies., Comment: accepted for publication in AJ

Observational studies of AGN in the mid-infrared regime are crucial to our understanding of AGN and their role in the evolution of galaxies. Mid-IR-based selection of AGN is complementary to more traditional techniques allowing for a more complete census of AGN activity across cosmic time. Mid-IR observations including time variability and spatially resolved imaging have given us unique insights into the nature of the obscuring structures around AGN. The wealth of fine structure, molecular, and dust features in the mid-IR allow us to simultaneously probe multiple components of the ISM allowing us to explore in detail the impact on the host galaxy by the presence of an AGN -- a crucial step toward understanding galaxy-SMBH co-evolution. This review gives a broad overview of this wide range of studies. It also aims to show the evolution of this field starting with its nascency in the 1960s, through major advances thanks to several generations of space-based and ground-based facilities, as well as the promise of upcoming facilities such as the {\sl James Webb Space Telescope (JWST)}., Comment: 37 pages, 20 figures; Invited review (published) for the special issue "Recent Advances in Infrared Galaxies and AGN", edited by Anna Sajina and Asantha R. Cooray, in Universe

The luminosity function (LF) of active galactic nuclei (AGN) probes the history of supermassive black hole assembly and growth across cosmic time. To mitigate selection biases, we present a consistent analysis of the AGN LFs derived for both X-ray and mid-infrared (MIR) selected AGN in the XMM-Large Scale Structure (XMM-LSS) field. There are 4268 AGN used to construct the MIR luminosity function (IRLF) and 3427 AGN used to construct the X-ray luminosity function (XLF), providing the largest census of the AGN population out to $z=4$ in both bands with significant reduction in uncertainties. We are able for the first time to see the knee of the IRLF at $z>2$ and observe a flattening of the faint-end slope as redshift increases. The bolometric luminosity density, a proxy for the cosmic black hole accretion history, computed from our LFs shows a peak at $z\approx2.25$ consistent with recent estimates of the peak in the star formation rate density (SFRD). However, at earlier epochs, the AGN luminosity density is flatter than the SFRD. If confirmed, this result suggests that the build up of black hole mass outpaces the growth of stellar mass in high mass systems at $z\gtrsim 2.5$. This is consistent with observations of redshift $z\sim 6$ quasars which lie above the local $M-\sigma$ relationship. The luminosity density derived from the IRLF is higher than that from the XLF at all redshifts. This is consistent with the dominant role of obscured AGN activity in the cosmic growth of supermassive black holes., Comment: 17 pages, 14 figures, 1 appendix

It is widely assumed that long-wavelength infrared (IR) emission from cold dust (T~20-40K) is a reliable tracer of star formation even in the presence of a bright active galactic nucleus (AGN). Based on radiative transfer (RT) models of clumpy AGN tori, hot dust emission from the torus contributes negligibly to the galaxy spectral energy distribution (SED) at $\lambda\ga100$ \micron. However, these models do not include AGN heating of host-galaxy-scale diffuse dust, which may have far-IR (FIR) colors comparable to cold diffuse dust heated by stars. To quantify the contribution of AGN heating to host-galaxy-scale cold dust emission at $\lambda\ga100$ \micron, we perform dust RT calculations on a simulated galaxy merger both including and excluding the bright AGN that it hosts. By differencing the SEDs yielded by RT calculations with and without AGN that are otherwise identical, we quantify the FIR cold dust emission arising solely from re-processed AGN photons. In extreme cases, AGN-heated host-galaxy-scale dust can increase galaxy-integrated FIR flux densities by factors of 2-4; star formation rates calculated from the FIR luminosity assuming no AGN contribution can overestimate the true value by comparable factors. Because the FIR colors of such systems are similar to those of purely star-forming galaxies and redder than torus models, broadband SED decomposition may be insufficient for disentangling the contributions of stars and heavily dust-enshrouded AGN in the most IR-luminous galaxies. We demonstrate how kpc-scale resolved observations can be used to identify deeply dust-enshrouded AGN with cool FIR colors when spectroscopic and/or X-ray detection methods are unavailable., Comment: 11 pages, 7 figures, Accepted to ApJ

The Spitzer Space Telescope revolutionized studies of Active Galactic Nuclei (AGNs). Its combined sensitivity and mapping speed at mid-infrared wavelengths revealed a substantial population of highly-obscured AGNs. This population implies a higher radiative accretion efficiency, and thus possibly a higher spin for black holes than indicated by surveys in the optical and X-ray. The unique mid-infrared spectrographic capability of Spitzer gave important insights into the distribution and nature of the dust surrounding AGNs, enabling the separation of AGN and starburst components, the detection of silicate features in emission from hot dust, and the identification of shocked gas associated with AGN activity. The sensitivity of Spitzer allowed almost complete identification of the host galaxies of samples of AGNs selected in the X-ray and radio. As we look forward to the James Webb Space Telescope, the lessons learned from Spitzer studies will inform observational programs with new and upcoming infrared facilities., Comment: Review article for the Nature Astronomy Spitzer retrospective

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., Comment: 26 pages, 12 figures, accepted for publication in ApJ

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 $320\%$) 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., Comment: Accpeted for publication in ApJ

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., Comment: 16 pages, 9 figures, accepted for publication in ApJ

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., Comment: Science White paper submitted to Astro2020 Decadal Survey

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., Comment: 14 pages, 10 figures, submitted to ApJ

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., Comment: Accepted for publication in ApJ. 21 pages, 11 figures

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., Comment: 12 pages, 9 figures, accepted for publication in MNRAS

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., Comment: accepted to ApJ, 22 pages, 12 figures

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., Comment: 31 pages, 12 figures, accepted for publication in ApJ

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)., Comment: 18 pages, 11 figures; ApJ accepted for publication

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., Comment: Accepted for publication in ApJ. Publicly available templates are at http://www.astro.umass.edu/~pope/Kirkpatrick2015. Due to size limitations, Appendix B could not be included in the arXiv submission

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., Comment: 6 Pages, 4 Figures, 1 Table. Accepted for publication in ApJL

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., Comment: 14 pages, accepted by ApJ

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., Comment: Accepted to the Proceedings of IAU Symposium 306, "Statistical Challenges in 21st Century Cosmology", September 2014. 3 pages, 1 figure

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., Comment: 17 pages, 15 figures, accepted for publication in the Monthly Notices of the Royal Astronomical Society (MNRAS)

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%., Comment: 16 Pages, 10 figures, Accepted for Publication at A&A

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., Comment: 16 pages, 7 figures, accepted for publication in MNRAS

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., Comment: 21 pages, 13 figures, ApJ accepted version

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., Comment: 28 pages, 8 figures, 7 tables; Accepted for publication in ApJ

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., Comment: 4 pages, 2 figures, to appear in the proceedings of The Spectral Energy Distribution of Galaxies, Preston, September 2011, eds R.J. Tuffs & C.C. Popescu

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., Comment: AJ, in press

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., Comment: 48 pages, 13 figures, ApJ, in press

We present new results on the physical nature of infrared-luminous sources at 0.5 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], Comment: 31 pages, 19 figures. Accepted for publication in ApJ

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., Comment: 16 pages, 13 figures, accepted for publication in ApJ

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., Comment: Accepted for publication in ApJ

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]., Comment: accepted for publication in the Astrophysical Journal

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 ~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 =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., Comment: 22 pages, 14 figures, accepted for publication in ApJ

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., Comment: ApJ, in press. Document revised to match the journal version

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., Comment: Accepted by ApJ. 25 pages, 5 figures and 2 Tables

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., Comment: 13 pages, 3 figures, accepted for publication in ApJL

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 =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., Comment: 60 pages, 15 figures, accepted for publication in ApJ

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 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)., Comment: Accepted by the ApJ. In vol. 658, March, 2007

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., Comment: 24 pages, 11 eps figures. Accepted for publication in MNRAS 2006. Full resolution paper available from http://www.phas.ubc.ca/~pope/MF1289rv2.pdf

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., Comment: 6 pages, 4 figures, Proceedings for the conference "The Spitzer Science Center 2005 Conference: Infrared Diagnostics of Galaxy Evolution", held in Pasadena (USA), in November 2005

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., Comment: 5 pages, 3 figures. Proceeding for the conference "Starbursts: From 30 Doradus to Lyman Break Galaxies", held in Cambridge (UK) in September, 2004

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., Comment: 8 pages, 4 figures, UMass/INAOE conference proceedings on Deep Millimeter Surveys, eds. J. Lowenthal and D. Hughes, references included for extragalactic backgrounds figure

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., Comment: 4 pages, including 1 figure, uses aa.cls, to appear in Astronomy & Astrophysics Letters, slight changes to some numbers and to match published version