Source author record

Anna Sajina

Anna Sajina appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

21works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

21 published item(s)

preprint2021arXiv

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

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.

preprint2020arXiv

Active Galactic Nuclei as seen by the Spitzer Space Telescope

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.

preprint2019arXiv

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

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.

preprint2016arXiv

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

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

preprint2015arXiv

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

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.

preprint2015arXiv

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

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.

preprint2015arXiv

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

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.

preprint2014arXiv

Robust Constraint of Luminosity Function Evolution Through MCMC Sampling

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.

preprint2013arXiv

Modeling Mid-Infrared Diagnostics of Obscured Quasars and Starbursts

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.

preprint2012arXiv

PAH Emission in Powerful High-Redshift Radio Galaxies

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.

preprint2012arXiv

Spectral energy distributions of quasars selected in the mid-infrared

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.

preprint2012arXiv

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

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.

preprint2012arXiv

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

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.

preprint2012arXiv

VLA/JVLA Monitoring of Bright Northern Radio Sources

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

preprint2011arXiv

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

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.

preprint2011arXiv

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

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]

preprint2011arXiv

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

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.

preprint2010arXiv

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

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.

preprint2009arXiv

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

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.

preprint2008arXiv

Spitzer Mid-Infrared Spectroscopy of Infrared Luminous Galaxies at z~2 III: Far-IR to Radio Properties and Optical Spectral Diagnostics

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

preprint2006arXiv

The Hubble Deep Field North SCUBA Super-map IV - Characterizing submillimetre galaxies using deep Spitzer imaging

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