Showing 1–50 of 73 results for author: Sajina, A

Halfway to the Peak: The JWST MIRI 5.6 micron number counts and source population

Authors: Leonid Sajkov, Anna Sajina, Alexandra Pope, Stacey Alberts, Lee Armus, Duncan Farrah, Jamie Lin, Danilo Marchesini, Jed McKinney, Sylvain Veilleux, Lin Yan, Jason Young

Abstract: 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… ▽ More 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\,$μ$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. △ Less

Submitted 6 June, 2024; originally announced June 2024.

Comments: 17 pages, 14 figures, submitted to ApJ

Estimating Galaxy Parameters with Self-Organizing Maps and the Effect of Missing Data

Authors: Valentina La Torre, Anna Sajina, Andy D. Goulding, Danilo Marchesini, Rachel Bezanson, Alan N. Pearl, Laerte Sodré Jr

Abstract: 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-limite… ▽ More 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. △ Less

Submitted 27 March, 2024; originally announced March 2024.

Comments: 22 pages, 12 figures, 1 table, accepted for publication in The Astronomical Journal

Halfway to the peak: Spatially resolved star formation and kinematics in a z=0.54 dusty galaxy with JWST/MIRI

Authors: Jason Young, Alexandra Pope, Anna Sajina, Lin Yan, Thiago S Goncalves, Miriam Eleazer, Stacey Alberts, Lee Armus, Matteo Bonato, Daniel A. Dale, Duncan Farrah, Carl Ferkinhoff, Christopher C. Hayward, Jed McKinney, Eric J. Murphy, Nicole Nesvadba, Patrick Ogle, Leonid Sajkov, Sylvain Veilleux

Abstract: 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 hydroc… ▽ More 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. △ Less

Submitted 12 October, 2023; v1 submitted 10 October, 2023; originally announced October 2023.

Comments: 10 pages, 5 figures, 1 table

The Spitzer Coverage of HSC-Deep with IRAC for Z studies (SHIRAZ) I: IRAC mosaics

Authors: Marianna Annunziatella, Anna Sajina, Mauro Stefanon, Danilo Marchesini, Mark Lacy, Ivo Labbe, Lilianna Houston, Rachel Bezanson, Eiichi Egami, Xiaohui Fan, Duncan Farrah, Jenny Greene, Andy Goulding, Yen-Ting Lin, Xin Liu, Thibaud Moutard, Yoshiaki Ono, Masami Ouchi, Marcin Sawicki, Jason Surace, Katherine Whitaker

Abstract: We present new Spitzer Infrared Array Camera (IRAC) 3.6 and 4.5μ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 improve… ▽ More We present new Spitzer Infrared Array Camera (IRAC) 3.6 and 4.5μ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σ depth of 23.7(23.3) at 3.6(4.5)μ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. △ Less

Submitted 22 May, 2023; originally announced May 2023.

Comments: accepted for publication in AJ

The Past and Future of Mid-Infrared Studies of AGN

Authors: Anna Sajina, Mark Lacy, Alexandra Pope

Abstract: 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 i… ▽ More 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)}. △ Less

Submitted 5 October, 2022; originally announced October 2022.

Comments: 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

Journal ref: Universe 2022, 8(7), 356

Consistent analysis of the AGN LF in X-ray and MIR in the XMM-LSS field

Authors: Jack Runburg, Duncan Farrah, Anna Sajina, Mark Lacy, Jenna Lidua, Evanthia Hatziminaoglou, W. N. Brandt, Chien-Ting J. Chen, Kristina Nyland, Raphael Shirley, D. L. Clements, Lura K. Pitchford

Abstract: 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 (… ▽ More 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-σ$ 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. △ Less

Submitted 30 November, 2021; originally announced December 2021.

Comments: 17 pages, 14 figures, 1 appendix

Dust-Enshrouded AGN can Dominate Host-Galaxy-Scale Cold-Dust Emission

Authors: Jed McKinney, Christopher C. Hayward, Lee J. Rosenthal, Juan Rafael Martinez-Galarza, Alexandra Pope, Anna Sajina, Howard A. Smith

Abstract: 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 $λ\ga100$ \micron. However, these mode… ▽ More 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 $λ\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 $λ\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. △ Less

Submitted 26 July, 2021; v1 submitted 23 March, 2021; originally announced March 2021.

Comments: 11 pages, 7 figures, Accepted to ApJ

The black hole masses of extremely luminous radio-WISE selected galaxies

Authors: E. R. Ferris, A. W. Blain, R. J. Assef, N. A. Hatch, A. Kimball, M. Kim, A. Sajina, A. Silva, D. Stern, T. Diaz-Santos, C-W. Tsai, D. Wylezalek

Abstract: We present near-IR photometry and spectroscopy of 30 extremely luminous radio and mid-IR selected galaxies. With bolometric luminosities exceeding $\sim10^{13}$ $\rm{L_{\odot}}$ and redshifts ranging from $z = 0.880-2.853$, we use VLT instruments X-shooter and ISAAC to investigate this unique population of galaxies. Broad multi-component emission lines are detected in 18 galaxies and we measure th… ▽ More We present near-IR photometry and spectroscopy of 30 extremely luminous radio and mid-IR selected galaxies. With bolometric luminosities exceeding $\sim10^{13}$ $\rm{L_{\odot}}$ and redshifts ranging from $z = 0.880-2.853$, we use VLT instruments X-shooter and ISAAC to investigate this unique population of galaxies. Broad multi-component emission lines are detected in 18 galaxies and we measure the near-IR lines $\rm{H\rmβ}$, $\text{[OIII]}\rmλ\rmλ4959,5007$ and $\rm{H\rmα}$ in six, 15 and 13 galaxies respectively, with 10 $\rm{Lyα}$ and five CIV lines additionally detected in the UVB arm. We use the broad $\text{[OIII]}\rmλ5007$ emission lines as a proxy for the bolometric AGN luminosity, and derive lower limits to supermassive black hole masses of $10^{7.9}$-$10^{9.4}$ $\text{M}_{\odot}$ with expectations of corresponding host masses of $10^{10.4}$-$10^{12.0}$ $\text{M}_{\odot}$. We measure $\rmλ_{Edd}$ > 1 for eight of these sources at a $2σ$ significance. Near-IR photometry and SED fitting are used to compare stellar masses directly. We detect both Balmer lines in five galaxies and use these to infer a mean visual extinction of $A_{V}$ = 2.68 mag. Due to non-detections and uncertainties in our $\rm{H\rmβ}$ emission line measurements, we simulate a broad $\rm{H\rmβ}$ line of FWHM = 1480 $\rm{kms^{-1}}$ to estimate extinction for all sources with measured $\rm{H\rmα}$ emission. We then use this to infer a mean $A_{V}=3.62$ mag, demonstrating the highly-obscured nature of these galaxies, with the consequence of increasing our estimates of black-hole masses by an 0.5 orders of magnitude in the most extreme and obscured cases. △ Less

Submitted 7 January, 2021; v1 submitted 6 January, 2021; originally announced January 2021.

Comments: Accepted for publication in MNRAS. 14 pages (+8 page appendix), 11 figures and 9 tables

A Spitzer survey of Deep Drilling Fields to be targeted by the Vera C. Rubin Observatory Legacy Survey of Space and Time

Authors: M. Lacy, J. A. Surace, D. Farrah, K. Nyland, J. Afonso, W. N. Brandt, D. L. Clements, C. D. P. Lagos, C. Maraston, J. Pforr, A. Sajina, M. Sako, M. Vaccari, G. Wilson, D. R. Ballantyne, W. A. Barkhouse, R. Brunner, R. Cane, T. E. Clarke, M. Cooper, A. Cooray, G. Covone, C. D'Andrea, A. E. Evrard, H. C. Ferguson , et al. (38 additional authors not shown)

Abstract: The Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) will observe several Deep Drilling Fields (DDFs) to a greater depth and with a more rapid cadence than the main survey. In this paper, we describe the ``DeepDrill'' survey, which used the Spitzer Space Telescope Infrared Array Camera (IRAC) to observe three of the four currently defined DDFs in two bands, centered on 3.6 $μ$m and… ▽ More The Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) will observe several Deep Drilling Fields (DDFs) to a greater depth and with a more rapid cadence than the main survey. In this paper, we describe the ``DeepDrill'' survey, which used the Spitzer Space Telescope Infrared Array Camera (IRAC) to observe three of the four currently defined DDFs in two bands, centered on 3.6 $μ$m and 4.5 $μ$m. These observations expand the area which was covered by an earlier set of observations in these three fields by the Spitzer Extragalactic Representative Volume Survey (SERVS). The combined DeepDrill and SERVS data cover the footprints of the LSST DDFs in the Extended Chandra Deep Field-South field (ECDFS), the ELAIS-S1 field (ES1), and the XMM Large-Scale Structure Survey field (XMM-LSS). The observations reach an approximate $5σ$ point-source depth of 2 $μ$Jy (corresponding to an AB magnitude of 23.1; sufficient to detect a 10$^{11} M_{\odot}$ galaxy out to $z\approx 5$) in each of the two bands over a total area of $\approx 29\,$deg$^2$. The dual-band catalogues contain a total of 2.35 million sources. In this paper we describe the observations and data products from the survey, and an overview of the properties of galaxies in the survey. We compare the source counts to predictions from the SHARK semi-analytic model of galaxy formation. We also identify a population of sources with extremely red ([3.6]$-$[4.5] $>1.2$) colours which we show mostly consists of highly-obscured active galactic nuclei. △ Less

Submitted 30 November, 2020; originally announced November 2020.

Comments: 20 pages, 9 figures; MNRAS in press

Active Galactic Nuclei as seen by the Spitzer Space Telescope

Authors: Mark Lacy, Anna Sajina

Abstract: 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-infra… ▽ More 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. △ Less

Submitted 12 August, 2020; originally announced August 2020.

Comments: Review article for the Nature Astronomy Spitzer retrospective

Journal ref: Nature Astronomy, 4, 352-363 (2020)

The ALPINE-ALMA [C II] Survey: [C II]158micron Emission Line Luminosity Functions at $z \sim 4-6$

Authors: Lin Yan, A. Sajina, F. Loiacono, G. Lagache, M. Bèthermin, A. Faisst, M. Ginolfi, O. Le Fèvre, C. Gruppioni, P. L. Capak, P. Cassata, D. Schaerer, J. D. Silverman, S. Bardelli, M. Dessauges-Zavadsky, A. Cimatti, N. P. Hathi, B. C. Lemaux, E. Ibar, G. C. Jones, A. M. Koekemoer, P. A. Oesch, M. Talia, F. Pozzi, D. A. Riechers , et al. (6 additional authors not shown)

Abstract: We present the [CII]158$μ$m line luminosity functions (LFs) at $z\sim4-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 [CII] detections and 43 are upper limits. This is by far the largest sample of [CII] detections which allows us to set constr… ▽ More We present the [CII]158$μ$m line luminosity functions (LFs) at $z\sim4-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 [CII] detections and 43 are upper limits. This is by far the largest sample of [CII] detections which allows us to set constraints to the volume density of [CII] emitters at $z\sim4-6$. But because this is a UV-selected sample, we are missing [CII]-bright but UV-faint sources making our constraints strict lower limits. Our derived LFs are statistically consistent with the $z\sim0$ [CII] LF at $10^{8.25} - 10^{9.75}L_\odot$. We compare our results with the upper limits of the [CII] LF derived from serendipitous sources in the ALPINE maps (Loiacono et al. 2020). We also infer the [CII] LFs based on published far-IR and CO LFs at $z\sim4-6$. Combining our robust lower limits with these additional estimates, we set further constraints to the true number density of [CII] emitters at $z\sim 4 - 6$. These additional LF estimates are largely above our LF at $L_{[CII]}>10^9L_{\odot}$, suggesting that UV-faint but [CII]-bright sources likely make a significant contributions to the [CII] emitter volume density. When we include all the LF estimates, we find that available model predictions underestimate the number densities of [CII] emitters at $z\sim4-6$. Finally, we set a constraint on the molecular gas mass density at $z\sim4-6$, with $ρ_{mol} \sim (2-7)\times10^7M_\odot$\,Mpc$^{-3}$. This is broadly consistent with previous studies. △ Less

Submitted 29 October, 2020; v1 submitted 8 June, 2020; originally announced June 2020.

Comments: 17 pages, 9 Figures, Update to match with the published version. Accepted for the publication in ApJ

The role of environment in galaxy evolution in the SERVS Survey I: density maps and cluster candidates

Authors: Nick Krefting, Anna Sajina, Mark Lacy, Kristina Nyland, Duncan Farrah, Behnam Darvish, Steven Duivenvoorden, Ken Duncan, Violeta Gonzalez-Perez, Claudia del P. Lagos, Seb Oliver, Raphael Shirley, Mattia Vaccari

Abstract: We use photometric redshifts derived from new $u$-band through 4.5$μ$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 s… ▽ More We use photometric redshifts derived from new $u$-band through 4.5$μ$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. △ Less

Submitted 4 December, 2019; originally announced December 2019.

Comments: 26 pages, 12 figures, accepted for publication in ApJ

Stellar and Dust Properties of a Complete Sample of Massive Dusty Galaxies at $1 \le z \le 4$ from MAGPHYS Modeling of UltraVISTA DR3 and Herschel Photometry

Authors: Nicholas S. Martis, Danilo M. Marchesini, Adam Muzzin, Mauro Stefanon, Gabriel Brammer, Elisabete da Cunha, Anna Sajina, Ivo Labbé

Abstract: 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) a… ▽ More 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-$β$ 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$μ$m is found to overestimate the SFR by a factor of 3-5 for the galaxies in our sample. △ Less

Submitted 18 July, 2019; originally announced July 2019.

Comments: Accpeted for publication in ApJ

Simulations Find Our Accounting of Dust-Obscured Star Formation May Be Incomplete

Authors: Eric Roebuck, Anna Sajina, Christopher C. Hayward, Nicholas Martis, Danilo Marchesini, Nicholas Krefting, Alexandra Pope

Abstract: 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 e… ▽ More 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. △ Less

Submitted 7 July, 2019; originally announced July 2019.

Comments: 16 pages, 9 figures, accepted for publication in ApJ

Simultaneous Measurements of Star Formation and Supermassive Black Hole Growth in Galaxies

Authors: Alexandra Pope, Lee Armus, Eric Murphy, Susanne Aalto, David Alexander, Philip Appleton, Amy Barger, Matt Bradford, Peter Capak, Caitlin Casey, Vassilis Charmandaris, Ranga Chary, Asantha Cooray, Jim Condon, Tanio Diaz Santos, Mark Dickinson, Duncan Farrah, Carl Ferkinhoff, Norman Grogin, Ryan Hickox, Allison Kirkpatrick, Kohno Kotaro, Allison Matthews, Desika Narayanan, Dominik Riechers , et al. (7 additional authors not shown)

Abstract: 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 for… ▽ More 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. △ Less

Submitted 12 March, 2019; originally announced March 2019.

Comments: Science White paper submitted to Astro2020 Decadal Survey

The Environments of Luminous Radio - WISE Selected Infrared Galaxies

Authors: J. I. Penney, A. W. Blain, D. Wylezalek, N. A. Hatch, C. Lonsdale, A. Kimball, R. J. Assef, J. J. Condon, P. R. M. Eisenhardt, S. F. Jones, M. Kim, M. Lacy, S. I. Muldrew, S. Petty, A. Sajina, A. Silva, D. Stern, T. Diaz-Santos, C-W. Tsai, J. Wu

Abstract: 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 o… ▽ More 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σ$ 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σ$ ($\rm5σ$) 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. △ Less

Submitted 16 November, 2018; originally announced November 2018.

Comments: 17 pages, 16 figures, 2 tables

A subarcsecond near-infrared view of massive galaxies at z > 1 with Gemini Multiconjugate Adaptive Optics

Authors: M. Lacy, K. Nyland, M. Mao, P. Jagannathan, J. Pforr, S. E. Ridgway, J. Afonso, D. Farrah, P. Guarnieri, E. Gonzales-Solares, M. J. Jarvis, C. Maraston, D. M. Nielsen, A. O. Petric, A. Sajina, J. A. Surace, M. Vaccari

Abstract: 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… ▽ More 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. △ Less

Submitted 9 July, 2018; originally announced July 2018.

Comments: 20 pages, AJ, in press

The AGN-Star Formation Connection: Future Prospects with JWST

Authors: Allison Kirkpatrick, Stacey Alberts, Alexandra Pope, Guillermo Barro, Matteo Bonato, Dale D. Kocevski, Pablo Perez-Gonzalez, George H. Rieke, Lucia Rodriguez-Munoz, Anna Sajina, Norman A. Grogin, Kameswara Bharadwaj Mantha, Viraj Pandya, Janine Pforr, Paola Santini

Abstract: 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 count… ▽ More 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 $λ_{\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. △ Less

Submitted 27 June, 2017; originally announced June 2017.

Comments: 14 pages, 10 figures, submitted to ApJ

A controlled study of cold dust content in galaxies from $z=0-2$

Authors: Allison Kirkpatrick, Alexandra Pope, Anna Sajina, Daniel A. Dale, Tanio Diaz-Santos, Christopher C. Hayward, Yong Shi, Rachel S. Somerville, Sabrina Stierwalt, Lee Armus, Jeyhan S. Kartaltepe, Dale D. Kocevski, Daniel H. McIntosh, David B. Sanders, Lin Yan

Abstract: 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 (… ▽ More 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. △ Less

Submitted 1 June, 2017; v1 submitted 30 May, 2017; originally announced May 2017.

Comments: Accepted for publication in ApJ. 21 pages, 11 figures

Does the evolution of the radio luminosity function of star-forming galaxies match that of the star-formation rate function?

Authors: Matteo Bonato, Mattia Negrello, Claudia Mancuso, Gianfranco De Zotti, Paolo Ciliegi, Zhen-Yi Cai, Andrea Lapi, Marcella Massardi, Anna Bonaldi, Anna Sajina, Vernesa Smolcic, Eva Schinnerer

Abstract: 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 re… ▽ More 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 $μ$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. △ Less

Submitted 5 May, 2017; v1 submitted 18 April, 2017; originally announced April 2017.

Comments: 12 pages, 9 figures, accepted for publication in MNRAS

An Application of Multi-band Forced Photometry to One Square Degree of SERVS: Accurate Photometric Redshifts and Implications for Future Science

Authors: Kristina Nyland, Mark Lacy, Anna Sajina, Janine Pforr, Duncan Farrah, Gillian Wilson, Jason Surace, Boris Haeussler, Mattia Vaccari, Matt Jarvis

Abstract: 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 c… ▽ More 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. △ Less

Submitted 24 May, 2017; v1 submitted 5 April, 2017; originally announced April 2017.

Comments: accepted to ApJ, 22 pages, 12 figures

Exploring the evolution of star formation and dwarf galaxy properties with JWST/MIRI serendipitous spectroscopic surveys

Authors: Matteo Bonato, Anna Sajina, Gianfranco De Zotti, Jed McKinney, Ivano Baronchelli, Mattia Negrello, Danilo Marchesini, Eric Roebuck, Heath Shipley, Noah Kurinsky, Alexandra Pope, Alberto Noriega-Crespo, Lin Yan, Allison Kirkpatrick

Abstract: 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… ▽ More 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. △ Less

Submitted 25 January, 2017; originally announced January 2017.

Comments: 31 pages, 12 figures, accepted for publication in ApJ

The Role of Star-Formation and AGN in Dust Heating of z=0.3-2.8 Galaxies - II. Informing IR AGN fraction estimates through simulations

Authors: Eric Roebuck, Anna Sajina, Christopher C. Hayward, Alexandra Pope, Allison Kirkpatrick, Lars Hernquist, Lin Yan

Abstract: 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$μ$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.… ▽ More 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$μ$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). △ Less

Submitted 3 October, 2016; v1 submitted 27 September, 2016; originally announced September 2016.

Comments: 18 pages, 11 figures; ApJ accepted for publication

The role of star-formation and AGN in dust heating of z = 0.3-2.8 galaxies - I. Evolution with redshift and luminosity

Authors: Allison Kirkpatrick, Alexandra Pope, Anna Sajina, Eric Roebuck, Lin Yan, Lee Armus, Tanio Diaz-Santos, Sabrina Stierwalt

Abstract: 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… ▽ More 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 $μ$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 $μ$Jy). Our decomposition technique and color diagnostics are relevant given upcoming observations with the James Webb Space Telescope. △ Less

Submitted 9 October, 2015; originally announced October 2015.

Comments: 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

The host galaxies of micro-Jansky radio sources

Authors: K. M. Luchsinger, M. Lacy, K. M. Jones, J. C. Mauduit, J. Pforr, J. A. Surace, M. Vaccari, D. Farrah, E. Gonzales-Solares, M. J. Jarvis, C. Maraston, L. Marchetti, S. Oliver, J. Afonso, D. Cappozi, A. Sajina

Abstract: 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 \;μ$Jy. 87% (1274/1467) of radio sources have identifications in SERVS to $AB\approx 23.1$ at 3.6 or 4.5… ▽ More 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 \;μ$Jy. 87% (1274/1467) of radio sources have identifications in SERVS to $AB\approx 23.1$ at 3.6 or 4.5$μ$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. △ Less

Submitted 4 July, 2015; originally announced July 2015.

Comments: 13 pages, accepted by AJ

ALMA detected overdensity of sub-mm sources around WISE/NVSS-selected z~2 dusty quasars

Authors: Andrea Silva, Anna Sajina, Carol Lonsdale, Mark Lacy

Abstract: 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… ▽ More 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. △ Less

Submitted 19 May, 2015; originally announced May 2015.

Comments: 6 Pages, 4 Figures, 1 Table. Accepted for publication in ApJL

The Spitzer mid-infrared AGN survey. II-the demographics and cosmic evolution of the AGN population

Authors: Mark Lacy, Susan E. Ridgway, Anna Sajina, Andreea O. Petric, Elinor L. Gates, Tanya Urrutia, Lisa J. Storrie-Lombardi

Abstract: 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 functi… ▽ More 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. △ Less

Submitted 16 January, 2015; originally announced January 2015.

Comments: 14 pages, accepted by ApJ

Robust Constraint of Luminosity Function Evolution Through MCMC Sampling

Authors: Noah Kurinsky, Anna Sajina

Abstract: 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… ▽ More 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. △ Less

Submitted 30 September, 2014; v1 submitted 3 July, 2014; originally announced July 2014.

Comments: Accepted to the Proceedings of IAU Symposium 306, "Statistical Challenges in 21st Century Cosmology", September 2014. 3 pages, 1 figure

The Spitzer mid-infrared AGN survey. I - optical and near-infrared spectroscopy of candidate obscured and normal AGN selected in the mid-infrared

Authors: M. Lacy, S. E. Ridgway, E. L. Gates, D. M. Nielsen, A. O. Petric, A. Sajina, T. Urrutia, S. Cox Drews, C. Harrison, N. Seymour, L. J. Storrie-Lombardi

Abstract: 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… ▽ More 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. △ Less

Submitted 19 August, 2013; originally announced August 2013.

Comments: ApJS, in press

VLA/JVLA Monitoring of Bright Northern Radio Sources

Authors: Noah Kurinsky, Anna Sajina, Bruce Partridge, Steve Myers, Xi Chen, Marcos López-Caniego

Abstract: 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… ▽ More 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%. △ Less

Submitted 16 November, 2012; originally announced November 2012.

Comments: 16 Pages, 10 figures, Accepted for Publication at A&A

Journal ref: A&A 549, A133 (2013)

The local luminosity function of star-forming galaxies derived from the Planck Early Release Compact Source Catalogue

Authors: Mattia Negrello, Marcel Clemens, Joaquin Gonzalez-Nuevo, Gianfranco De Zotti, Laura Bonavera, Giorgio Cosco, Gianpaolo Guarese, Luca Boaretto, Stephen Serjeant, Luigi Toffolatti, Andrea Lapi, Matthieu Bethermin, Guillaume Castex, Dave L. Clements, Jacques Delabrouille, Herve' Dole, Alberto Franceschini, Reno Mandolesi, Lucia Marchetti, Bruce Partridge, Anna Sajina

Abstract: 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… ▽ More 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. △ Less

Submitted 16 November, 2012; originally announced November 2012.

Comments: 17 pages, 15 figures, accepted for publication in the Monthly Notices of the Royal Astronomical Society (MNRAS)

PAH Emission in Powerful High-Redshift Radio Galaxies

Authors: Jason Ian Rawlings, Nicholas Seymour, Mathew Page, Carlos De Breuck, Daniel Stern, Myrto Symeonidis, Phil Appleton, Arjun Dey, Mark Dickinson, Minh Huynh, Emeric Le Floc'h, Matt Lehnert, James Mullaney, Nicole Nesvadba, Patrick Ogle, Anna Sajina, Joel Vernet, Andrew Zirm

Abstract: 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… ▽ More 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. △ Less

Submitted 12 November, 2012; originally announced November 2012.

Comments: 16 pages, 7 figures, accepted for publication in MNRAS

Modeling Mid-Infrared Diagnostics of Obscured Quasars and Starbursts

Authors: Gregory F. Snyder, Christopher C. Hayward, Anna Sajina, Patrik Jonsson, Thomas J. Cox, Lars Hernquist, Philip F. Hopkins, Lin Yan

Abstract: 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 angl… ▽ More 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. △ Less

Submitted 11 April, 2013; v1 submitted 23 October, 2012; originally announced October 2012.

Comments: 21 pages, 13 figures, ApJ accepted version

Journal ref: Astrophys.J. 768 (2013) 168

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

Authors: A. Sajina, L. Yan, D. Fadda, K. Dasyra, M. Huynh

Abstract: 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 500… ▽ More 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. △ Less

Submitted 20 July, 2012; originally announced July 2012.

Comments: 50 pages, 15 figures, ApJ in press

Planck intermediate results. VII. Statistical properties of infrared and radio extragalactic sources from the Planck Early Release Compact Source Catalogue at frequencies between 100 and 857 GHz

Authors: Planck Collaboration, P. A. R. Ade, N. Aghanim, F. Argüeso, M. Arnaud, M. Ashdown, F. Atrio-Barandela, J. Aumont, C. Baccigalupi, A. Balbi, A. J. Banday, R. B. Barreiro, E. Battaner, K. Benabed, A. Benoît, J. -P. Bernard, M. Bersanelli, M. Bethermin, R. Bhatia, A. Bonaldi, J. R. Bond, J. Borrill, F. R. Bouchet, C. Burigana, P. Cabella , et al. (166 additional authors not shown)

Abstract: (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.… ▽ More (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. △ Less

Submitted 22 January, 2013; v1 submitted 19 July, 2012; originally announced July 2012.

Comments: in press, A&A, 2013 -- matches published version (updated text for clarity after referee's report; conclusions unchanged) 20 pages (2 column), 14 figs (+3 in appendix), 9 tables

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

Authors: J. -C. Mauduit, M. Lacy, D. Farrah, J. A. Surace, M. Jarvis, S. Oliver, C. Maraston, M. Vaccari, L. Marchetti, G. Zeimann, E. A. Gonzalez-Solares, J. Pforr, A. O. Petric, B. Henriques, P. A. Thomas, J. Afonso, A. Rettura, G. Wilson, J. T. Falder, J. E. Geach, M. Huynh, R. P. Norris, N. Seymour, G. T. Richards, S. A. Stanford , et al. (59 additional authors not shown)

Abstract: 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 fun… ▽ More 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. △ Less

Submitted 17 September, 2012; v1 submitted 18 June, 2012; originally announced June 2012.

Comments: 21 pages, 23 figures, published in PASP. This version differs slightly from PASP, mainly due to formatting issues. Figure 2 was also corrected from the previous version

Journal ref: Publications of the Astronomical Society of the Pacific, 2012, Volume 124, issue 917, pp.714-736

The Nature of LoBAL QSOs: I. SEDs and mid-infrared spectral properties

Authors: Mariana S. Lazarova, Gabriela Canalizo, Mark Lacy, Anna Sajina

Abstract: 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… ▽ More 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. △ Less

Submitted 20 June, 2012; v1 submitted 8 June, 2012; originally announced June 2012.

Comments: 28 pages, 8 figures, 7 tables; Accepted for publication in ApJ

Spectral energy distributions of quasars selected in the mid-infrared

Authors: Mark Lacy, Anna Sajina, Susan E. Ridgway, Danielle M. Nielsen, Tanya Urrutia, Duncan Farrah, Elinor L. Gates

Abstract: 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… ▽ More 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. △ Less

Submitted 20 March, 2012; originally announced March 2012.

Comments: 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

The stellar, molecular gas and dust content of the host galaxies of two z~2.8 dust obscured quasars

Authors: Mark Lacy, Andreea O. Petric, Alejo Martinez-Sansigre, Susan E. Ridgway, Anna Sajina, Tanya Urrutia, Duncan Farrah

Abstract: 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 c… ▽ More 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. △ Less

Submitted 11 October, 2011; originally announced October 2011.

Comments: AJ, in press

Simultaneous Planck, Swift, and Fermi observations of X-ray and gamma-ray selected blazars

Authors: P. Giommi, G. Polenta, A. Lahteenmaki, D. J. Thompson, M. Capalbi, S. Cutini, D. Gasparrini, J. Gonzalez-Nuevo, J. Leon-Tavares, M. Lopez-Caniego, M. N. Mazziotta, C. Monte, M. Perri, S. Raino, G. Tosti, A. Tramacere, F. Verrecchia, H. D. Aller, M. F. Aller, E. Angelakis, D. Bastieri, A. Berdyugin, A. Bonaldi, L. Bonavera, C. Burigana , et al. (70 additional authors not shown)

Abstract: 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… ▽ More 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 <α>~-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. △ Less

Submitted 23 May, 2012; v1 submitted 4 August, 2011; originally announced August 2011.

Comments: Version accepted by A&A. Joint Planck, Swift, and Fermi collaborations paper

Journal ref: A&A 541, A160, 2012

High frequency radio SEDs and polarization fractions of sources in an ACT survey field

Authors: Anna Sajina, Bruce Partridge, Tyler Evans, Shannon Stefl, Nicholas Vechik, Steve Myers, Simon Dicker, Phillip Korngut

Abstract: 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 obtai… ▽ More 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. △ Less

Submitted 24 January, 2011; originally announced January 2011.

Comments: 48 pages, 13 figures, ApJ, in press

HST/NICMOS Imaging of Bright High-Redshift 24μm-selected Galaxies: Merging Properties

Authors: Michel Zamojski, Lin Yan, Kalliopi Dasyra, Anna Sajina, Jason Surace, Tim Heckman, George Helou

Abstract: 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 struc… ▽ More 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 (τ_{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] △ Less

Submitted 12 January, 2011; originally announced January 2011.

Comments: 31 pages, 19 figures. Accepted for publication in ApJ

Planck Early Results: Statistical properties of extragalactic radio sources in the Planck Early Release Compact Source Catalogue

Authors: Planck Collaboration, P. A. R. Ade, N. Aghanim, F. Argüeso, M. Arnaud, M. Ashdown, J. Aumont, C. Baccigalupi, A. Balbi, A. J. Banday, R. B. Barreiro, J. G. Bartlett, E. Battaner, K. Benabed, J. -P. Bernard, M. Bersanelli, R. Bhatia, A. Bonaldi, L. Bonavera, J. R. Bond, J. Borrill, F. R. Bouchet, M. Bucher, C. Burigana, P. Cabella , et al. (179 additional authors not shown)

Abstract: 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 e… ▽ More 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. △ Less

Submitted 28 February, 2011; v1 submitted 11 January, 2011; originally announced January 2011.

Comments: Planck Collaboration 2011i, A&A accepted, 11 pages, 7 figures

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

Authors: Planck Collaboration, P. A. R. Ade, N. Aghanim, M. Arnaud, M. Ashdown, J. Aumont, C. Baccigalupi, A. Balbi, A. J. Banday, R. B. Barreiro, J. G. Bartlett, E. Battaner, K. Benabed, A. Benoît, J. -P. Bernard, M. Bersanelli, R. Bhatia, A. Bonaldi, L. Bonavera, J. R. Bond, J. Borrill, F. R. Bouchet, M. Bucher, C. Burigana, R. C. Butler , et al. (206 additional authors not shown)

Abstract: (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 firs… ▽ More (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. △ Less

Submitted 8 July, 2011; v1 submitted 11 January, 2011; originally announced January 2011.

Comments: Updated to match A&A press version

Star formation and dust obscuration in the tidally distorted galaxy NGC 2442

Authors: Anna Pancoast, Anna Sajina, Mark Lacy, Alberto Noriega-Crespo, Jeonghee Rho

Abstract: 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Å and 2300Å 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.… ▽ More 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Å and 2300Å 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Å) 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α, 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α emission is largely concentrated inside that ring and shows a vague spiral structure in the rest of the galaxy. △ Less

Submitted 9 September, 2010; originally announced September 2010.

Comments: 16 pages, 13 figures, accepted for publication in ApJ

Ultra-deep Spitzer Mid-Infrared Spectroscopy of LIRGs and ULIRGs at z ~ 1-2

Authors: D. Fadda, L. Yan, G. Lagache, A. Sajina, D. Lutz, S. Wuyts, D. T. Frayer, D. Marcillac, E. Le Floc'h, K. Caputi, H. W. W. Spoon, S. Veilleux, A. Blain, G. Helou

Abstract: 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 e… ▽ More 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). △ Less

Submitted 14 June, 2010; originally announced June 2010.

Comments: 45 pages, 36 figures, 2 tables, accepted for publication on ApJ

Detections of CO Molecular Gas in 24um-Bright ULIRGs at z~2 in the Spitzer First Look Survey

Authors: Lin Yan, L. ~J. Tacconi, N. Fiolet, A. Sajina, A. Omont, D. Lutz, M. Zamojski, R. Neri, P. Cox, K. ~M. Dasyra

Abstract: We present CO observations of 9 ULIRGs at z~2 with S(24μ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 integrat… ▽ More We present CO observations of 9 ULIRGs at z~2 with S(24μ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. △ Less

Submitted 11 March, 2010; originally announced March 2010.

Comments: 15 pages, 8 figures, accepted for publication in ApJ.

X-ray Constraints on the AGN Properties in Spitzer-IRS identified z~2 Ultraluminous Infrared Galaxies

Authors: F. E. Bauer, Lin Yan, A. Sajina, D. M. Alexander

Abstract: 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. Th… ▽ More 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. △ Less

Submitted 21 December, 2009; originally announced December 2009.

Comments: Accepted for publication in ApJ. 16 pages, 4 color figs, and 4 tables in latex using emuluateapj style

Journal ref: Astrophys.J.710:212-226,2010

An Infrared Comparison of Type-1 and Type-2 Quasars

Authors: Kyle D. Hiner, Gabriela Canalizo, Mark Lacy, Anna Sajina, Lee Armus, Susan Ridgway, Lisa Storrie-Lombardi

Abstract: 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 differ… ▽ More 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. △ Less

Submitted 7 October, 2009; originally announced October 2009.

Comments: Accepted for publication in ApJ

Journal ref: Astrophys.J.706:508-515,2009

Detections of water ice, hydrocarbons, and 3.3um PAH in z~2 ULIRGs

Authors: A. Sajina, H. Spoon, L. Yan, M. Imanishi, D. Fadda, M. Elitzur

Abstract: 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' sugge… ▽ More 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. △ Less

Submitted 30 July, 2009; originally announced July 2009.

Comments: 39 pages, 14 figures, accepted for publication in ApJ