Source author record

David M. Alexander

David M. Alexander 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

40works
9topics
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

40 published item(s)

preprint2026arXiv

The Quasar Feedback Survey: Revealing the importance of sensitive radio imaging for AGN identification deeper into the radio-quiet regime

We present new sub-arcsecond ($\sim$0.3-1 arcsec; $\sim$1--3\,kpc) VLA imaging at 1.4\,GHz and 6\,GHz of 29 optically-selected, [O~{\sc iii}] luminous ($L_{\rm [O III]}$ > 10$^{42.1}$\,erg\,s$^{-1}$), $z<0.2$ quasars drawn from the expanded Quasar Feedback Survey (QFeedS; with $L_\mathrm{1.4\,GHz} = 10^{22.6}$--10$^{26.3}$\,W\,Hz$^{-1}$). These 29 new objects occupy the low end of the radio-power distribution ($L_\mathrm{1.4\,GHz}$=$10^{22.63}$--10$^{23.45}$\,W\,Hz$^{-1}$) in the QFeedS sample and are nominally `radio quiet'. Despite this, we find widespread evidence of AGN-driven synchrotron activity. Nearly $\sim 31\,$per\,cent exhibit resolved radio structures on $\sim$0.1--20\,kpc scales consistent with compact jets or wind-driven outflows, and $\sim 90\,$per\,cent display steep spectra ($α\lesssim -1$) indicative of optically thin synchrotron emission. Combining morphology, spectral index and brightness-temperature diagnostics, at least $\sim38\,$per\,cent of the sample show clear AGN signatures that cannot be explained by star formation alone. These constitute the first results from the expanded QFeedS (now 71 quasars spanning $\approx 4$ dex in radio power) and demonstrate that compact, low-power jets and AGN shocks are common deep inside the radio-quiet regime. A thorough understanding of feedback processes from quasars, deep into the `radio-quiet' regime, will be obtained by connecting these high resolution radio observations with multi-wavelength observations.

preprint2023arXiv

The Target-selection Pipeline for the Dark Energy Spectroscopic Instrument

In 2021 May, the Dark Energy Spectroscopic Instrument (DESI) began a 5 yr survey of approximately 50 million total extragalactic and Galactic targets. The primary DESI dark-time targets are emission line galaxies (ELGs), luminous red galaxies (LRGs) and quasars (QSOs). In bright time, DESI will focus on two surveys known as the Bright Galaxy Survey (BGS) and the Milky Way Survey (MWS). DESI also observes a selection of "secondary" targets for bespoke science goals. This paper gives an overview of the publicly available pipeline (desitarget) used to process targets for DESI observations. Highlights include details of the different DESI survey targeting phases, the targeting ID (TARGETID) used to define unique targets, the bitmasks used to indicate a particular type of target, the data model and structure of DESI targeting files, and examples of how to access and use the desitarget code base. This paper will also describe "supporting" DESI target classes, such as standard stars, sky locations, and random catalogs that mimic the angular selection function of DESI targets. The DESI target selection pipeline is complex and sizable; this paper attempts to summarize the most salient information required to understand and work with DESI targeting data.

preprint2022arXiv

Host Dark Matter Halos of SDSS Red and Blue Quasars: No Significant Difference in Large-scale Environment

The observed optical colors of quasars are generally interpreted in one of two frameworks: unified models which attribute color to random orientation of the accretion disk along the line-of-sight, and evolutionary models which invoke connections between quasar systems and their environments. We test these schema by probing the dark matter halo environments of optically-selected quasars as a function of $g-i$ optical color by measuring the two-point correlation functions of $\sim$ 0.34 million eBOSS quasars as well as the gravitational deflection of cosmic microwave background photons around $\sim$ 0.66 million XDQSO photometric quasar candidates. We do not detect a trend of halo bias with optical color through either analysis, finding that optically-selected quasars at $0.8 < z < 2.2$ occupy halos of characteristic mass $M_{h}\sim 3\times 10^{12} \ h^{-1} M_{\odot}$ regardless of their color. This result implies that a quasar's large-scale halo environment is not strongly connected to its observed optical color. We also confirm findings of fundamental differences in the radio properties of red and blue quasars by stacking 1.4 GHz FIRST images at their positions, suggesting the observed differences cannot be attributed to orientation. Instead, the differences between red and blue quasars likely arise on nuclear-galactic scales, perhaps owing to reddening by a nuclear dusty wind. Finally, we show that optically-selected quasars' halo environments are also independent of their $r-W2$ optical-infrared colors, while previous work has suggested that mid-infrared-selected obscured quasars occupy more massive halos. We discuss implications of this result for models of quasar and galaxy co-evolution.

preprint2022arXiv

Localizing narrow Fe K$α$ emission within bright AGN

The 6.4 keV Fe Ka emission line is a ubiquitous feature in X-ray spectra of AGN, and its properties track the interaction between the variable primary X-ray continuum and the surrounding structure from which it arises. We clarify the nature and origin of the narrow Fe Ka emission using X-ray spectral, timing, and imaging constraints, plus possible correlations to AGN and host galaxy properties, for 38 bright nearby AGN ($z<0.5$) from the BAT AGN Spectroscopic Survey. Modeling Chandra and XMM-Newton spectra, we computed line full-width half-maxima (FWHMs) and constructed Fe Ka line and 2-10 keV continuum light curves. The FWHM provides one estimate of the Fe Ka emitting region size, RFeKa, assuming virial motion. A second estimate comes from comparing the degree of correlation between the variability of the continuum and line-only light curves, compared to simulated light curves. Finally, we extracted Chandra radial profiles to place upper limits on RFeKa. We found that for 90% (21/24) of AGN with FWHM measurements, RFeKa is smaller than the fiducial dust sublimation radius, Rsub. Despite a wide range of variability properties, the constraints on the Fe Ka photon reprocessor size independently confirm that RFeKa is smaller than Rsub in 83% of AGN. Finally, the imaging analysis yields loose upper limits for all but two sources; notably, the Circinus Galaxy and NGC 1068 show significant but subdominant extended Fe Ka emission out to $\sim$100 and $\sim$800 pc, respectively. Based on independent constraints, we conclude that the majority of the narrow Fe Ka emission in typical AGN predominantly arises from regions smaller than and presumably inside Rsub, and thus it is associated either with the outer broad line region or outer accretion disk. However, the large diversity of continuum and narrow Fe Ka variability properties are not easily accommodated by a universal scenario.

preprint2022arXiv

Low Rank Approximation of Dual Complex Matrices

Dual complex numbers can represent rigid body motion in 2D spaces. Dual complex matrices are linked with screw theory, and have potential applications in various areas. In this paper, we study low rank approximation of dual complex matrices. We define $2$-norm for dual complex vectors, and Frobenius norm for dual complex matrices. These norms are nonnegative dual numbers. We establish the unitary invariance property of dual complex matrices. We study eigenvalues of square dual complex matrices, and show that an $n \times n$ dual complex Hermitian matrix has exactly $n$ eigenvalues, which are dual numbers. We present a singular value decomposition (SVD) theorem for dual complex matrices, define ranks and appreciable ranks for dual complex matrices, and study their properties. We establish an Eckart-Young like theorem for dual complex matrices, and present an algorithm framework for low rank approximation of dual complex matrices via truncated SVD. The SVD of dual complex matrices also provides a basic tool for Principal Component Analysis (PCA) via these matrices. Numerical experiments are reported.

preprint2022arXiv

The Black Hole-Galaxy Connection: Interplay between Feedback, Obscuration, and Host Galaxy Substructure

There is growing evidence for physical influence between supermassive black holes and their host galaxies. We present a case study of nearby galaxy NGC 7582, for which we find evidence that galactic substructure plays an important role in affecting the collimation of ionized outflows as well as contributing to the heavy active galactic nucleus (AGN) obscuration. This result contrasts with a simple, small-scale AGN torus model, according to which AGN wind collimation may take place inside the torus itself, at subparsec scale. Using 3D spectroscopy with the MUSE instrument, we probe the kinematics of the stellar and ionized gas components as well as the ionization state of the gas from a combination of emission line ratios. We report for the first time a kinematically distinct core (KDC) in NGC 7582, on a scale of ~600pc. This KDC coincides spatially with dust lanes and starbursting complexes previously observed. We interpret it as a circumnuclear ring of stars and dusty, gas-rich material. We obtain a clear view of the outflowing cones over kpc scales, and demonstrate that they are predominantly photoionized by the central engine. We detect the back cone (behind the galaxy), and confirm previous results of a large nuclear obscuration of both the stellar continuum and HII regions. While we tentatively associate the presence of the KDC to a large-scale bar and/or a minor galaxy merger, we stress the importance of gaining a better understanding of the role of galaxy substructure in controlling the fueling, feedback and obscuration of AGN.

preprint2021arXiv

The Black Hole Mass Function Across Cosmic Times I. Stellar Black Holes and Light Seed Distribution

This is the first paper in a series aimed at modeling the black hole (BH) mass function, from the stellar to the intermediate to the (super)massive regime. In the present work we focus on stellar BHs and provide an ab-initio computation of their mass function across cosmic times. Specifically, we exploit the state-of-the-art stellar and binary evolutionary code \texttt{SEVN}, and couple its outputs with redshift-dependent galaxy statistics and empirical scaling relations involving galaxy metallicity, star-formation rate and stellar mass. The resulting relic mass function ${\rm d}N/{\rm d}V{\rm d}\log m_\bullet$ as a function of the BH mass $m_\bullet$ features a rather flat shape up to $m_\bullet\approx 50\, M_\odot$ and then a log-normal decline for larger masses, while its overall normalization at a given mass increases with decreasing redshift. We highlight the contribution to the local mass function from isolated stars evolving into BHs and from binary stellar systems ending up in single or binary BHs. We also include the distortion on the mass function induced by binary BH mergers, finding that it has a minor effect at the high-mass end. We estimate a local stellar BH relic mass density of $ρ_\bullet\approx 5\times 10^7\, M_\odot$ Mpc$^{-3}$, which exceeds by more than two orders of magnitude that in supermassive BHs; this translates into an energy density parameter $Ω_\bullet\approx 4\times 10^{-4}$, implying that the total mass in stellar BHs amounts to $\lesssim 1\%$ of the local baryonic matter. We show how our mass function for merging BH binaries compares with the recent estimates from gravitational wave observations by LIGO/Virgo, and discuss the possible implications for dynamical formation of BH binaries in dense environments like star clusters. [abridged]

preprint2020arXiv

Galaxy mergers in EAGLE do not induce a significant amount of black hole growth yet do increase the rate of luminous AGN

We investigate the connection between galaxy--galaxy mergers and enhanced black hole (BH) growth using the cosmological hydrodynamical EAGLE simulation. We do this via three methods of analysis, investigating: the merger fraction of AGN, the AGN fraction of merging systems and the AGN fraction of galaxies with close companions. In each case, we find an increased abundance of AGN within merging systems relative to control samples of inactive or isolated galaxies (by up to a factor of $\approx 3$ depending on the analysis method used), confirming that mergers are enhancing BH accretion rates for at least a subset of the galaxy population. The greatest excess of AGN triggered via a merger are found in lower mass ($M_* \sim 10^{10}$~\Msol) gas rich ($f_{\mathrm{gas}} > 0.2$) central galaxies with lower mass BHs ($M_{\mathrm{BH}} \sim 10^{7}$~\Msol) at lower redshifts ($z<1$). We find no enhancement of AGN triggered via mergers in more massive galaxies ($M_* \gtrsim 10^{11}$~\Msol). The enhancement of AGN is not uniform throughout the phases of a merger, and instead peaks within the early \emph{remnants} of merging systems (typically lagging $\approx 300$~Myr post-coalescence of the two galaxies at $z=0.5$). We argue that neither major ($M_{\mathrm{*,1}} / M_{\mathrm{*,2}} \geq \frac{1}{4}$) nor minor mergers ($\frac{1}{10} < M_{\mathrm{*,1}} / M_{\mathrm{*,2}} < \frac{1}{4}$) are statistically relevant for enhancing BH masses globally. Whilst at all redshifts the galaxies experiencing a merger have accretion rates that are on average 2--3 times that of isolated galaxies, the majority of mass that is accreted onto BHs occurs outside the periods of a merger. We compute that on average no more than 15\% of a BHs final day mass comes from the enhanced accretion rates triggered via a merger

preprint2020arXiv

X-ray Binary Luminosity Function Scaling Relations in Elliptical Galaxies: Evidence for Globular Cluster Seeding of Low-Mass X-ray Binaries in Galactic Fields

We investigate X-ray binary (XRB) luminosity function (XLF) scaling relations for Chandra detected populations of low-mass XRBs (LMXBs) within the footprints of 24 early-type galaxies. Our sample includes Chandra and HST observed galaxies at D < 25 Mpc that have estimates of the globular cluster (GC) specific frequency (SN) reported in the literature. As such, we are able to directly classify X-ray-detected sources as being either coincident with unrelated background/foreground objects, GCs, or sources that are within the fields of the galaxy targets. We model the GC and field LMXB population XLFs for all galaxies separately, and then construct global models characterizing how the LMXB XLFs vary with galaxy stellar mass and SN. We find that our field LMXB XLF models require a component that scales with SN, and has a shape consistent with that found for the GC LMXB XLF. We take this to indicate that GCs are "seeding" the galactic field LMXB population, through the ejection of GC-LMXBs and/or the diffusion of the GCs in the galactic fields themselves. However, we also find that an important LMXB XLF component is required for all galaxies that scales with stellar mass, implying that a substantial population of LMXBs are formed "in situ," which dominates the LMXB population emission for galaxies with SN < 2. For the first time, we provide a framework quantifying how directly-associated GC LMXBs, GC-seeded LMXBs, and in-situ LMXBs contribute to LMXB XLFs in the broader early-type galaxy population.

preprint2019arXiv

Planck far-infrared detection of Hyper Suprime-Cam protoclusters at $\bf z\sim4$: hidden AGN and star formation activity

We perform a stacking analysis of {\it Planck}, {\it AKARI}, Infrared Astronomical Satellite ($IRAS$), Wide-field Infrared Survey Eplorer ($WISE$), and {\it Herschel} images of the largest number of (candidate) protoclusters at $z\sim3.8$ selected from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). Stacking the images of the $179$ candidate protoclusters, the combined infrared (IR) emission of the protocluster galaxies in the observed $12-850~μ$m wavelength range is successfully detected with $>5σ$ significance (at $Planck$). This is the first time that the average IR spectral energy distribution (SED) of a protocluster has been constrained at $z\sim4$. The observed IR SEDs of the protoclusters exhibit significant excess emission in the mid-IR compared to that expected from typical star-forming galaxies (SFGs). They are reproduced well using SED models of intense starburst galaxies with warm/hot dust heated by young stars, or by a population of active galactic nuclei (AGN)/SFG composites. For the pure star-forming model, a total IR (from 8 to 1000 $μ$m) luminosity of $19.3_{-4.2}^{+0.6}\times10^{13}~L_{\odot}$ and a star formation rate (SFR) of $16.3_{-7.8}^{+1.0}\times10^3~M_{\odot}$ yr$^{-1}$ are found whereas for the AGN/SFG composite model, $5.1_{-2.5}^{+2.5}\times10^{13}~L_{\odot}$ and $2.1^{+6.3}_{-1.7}\times10^3~M_{\odot}$ yr$^{-1}$ are found. Uncertainty remaining in the total SFRs; however, the IR luminosities of the most massive protoclusters are likely to continue increasing up to $z\sim4$. Meanwhile, no significant IR flux excess is observed around optically selected QSOs at similar redshifts, which confirms previous results. Our results suggest that the $z\sim4$ protoclusters trace dense, intensely star-forming environments that may also host obscured AGNs missed by the selection in the optical.

preprint2019arXiv

Probing black hole accretion tracks, scaling relations and radiative efficiencies from stacked X-ray active galactic nuclei

The masses of supermassive black holes at the centres of local galaxies appear to be tightly correlated with the mass and velocity dispersions of their galactic hosts. However, the local Mbh-Mstar relation inferred from dynamically measured inactive black holes is up to an order-of-magnitude higher than some estimates from active black holes, and recent work suggests that this discrepancy arises from selection bias on the sample of dynamical black hole mass measurements. In this work we combine X-ray measurements of the mean black hole accretion luminosity as a function of stellar mass and redshift with empirical models of galaxy stellar mass growth, integrating over time to predict the evolving Mbh-Mstar relation. The implied relation is nearly independent of redshift, indicating that stellar and black hole masses grow, on average, at similar rates. Matching the de-biased local Mbh-Mstar relation requires a mean radiative efficiency ~0.15, in line with theoretical expectations for accretion onto spinning black holes. However, matching the "raw" observed relation for inactive black holes requires a mean radiative efficiency around 0.02, far below theoretical expectations. This result provides independent evidence for selection bias in dynamically estimated black hole masses, a conclusion that is robust to uncertainties in bolometric corrections, obscured active black hole fractions, and kinetic accretion efficiency. For our fiducial assumptions, they favour moderate-to-rapid spins of typical supermassive black holes, to achieve a mean radiative efficiency ~0.12-0.20. Our approach has similarities to the classic Soltan analysis, but by using galaxy-based data instead of integrated quantities we are able to focus on regimes where observational uncertainties are minimized.

preprint2016arXiv

NuSTAR Hard X-ray Survey of the Galactic Center Region II: X-ray Point Sources

We present the first survey results of hard X-ray point sources in the Galactic Center (GC) region by NuSTAR. We have discovered 70 hard (3-79 keV) X-ray point sources in a 0.6 deg^2 region around Sgr A* with a total exposure of 1.7 Ms, and 7 sources in the Sgr B2 field with 300 ks. We identify clear Chandra counterparts for 58 NuSTAR sources and assign candidate counterparts for the remaining 19. The NuSTAR survey reaches X-ray luminosities of ~4 x and ~8 x 10^32 erg s^-1 at the GC (8 kpc) in the 3-10 and 10-40 keV bands, respectively. The source list includes three persistent luminous X-ray binaries and the likely run-away pulsar called the Cannonball. New source-detection significance maps reveal a cluster of hard (>10 keV) X-ray sources near the Sgr A diffuse complex with no clear soft X-ray counterparts. The severe extinction observed in the Chandra spectra indicates that all the NuSTAR sources are in the central bulge or are of extragalactic origin. Spectral analysis of relatively bright NuSTAR sources suggests that magnetic cataclysmic variables constitute a large fraction (>40-60%). Both spectral analysis and logN-logS distributions of the NuSTAR sources indicate that the X-ray spectra of the NuSTAR sources should have kT > 20 keV on average for a single temperature thermal plasma model or an average photon index of Gamma = 1.5 - 2 for a power-law model. These findings suggest that the GC X-ray source population may contain a larger fraction of X-ray binaries with high plasma temperatures than the field population.

preprint2016arXiv

The distribution of radioactive $^{44}$Ti in Cassiopeia A

The distribution of elements produced in the inner-most layers of a supernova explosion is a key diagnostic for studying the collapse of massive stars. Here we present the results of a 2.4 Ms \textit{NuSTAR} observing campaign aimed at studying the supernova remnant Cassiopeia A (Cas A). We perform spatially-resolved spectroscopic analyses of the $^{44}$Ti ejecta which we use to determine the Doppler shift and thus the three-dimensional (3D) velocities of the $^{44}$Ti ejecta. We find an initial $^{44}$Ti mass of 1.54 $\pm$ 0.21 $\times 10^{-4}$ M$_{\odot}$ which has a present day average momentum direction of 340$^{\circ}$ $\pm$ 15$^{\circ}$ projected on to the plane of the sky (measured clockwise from Celestial North) and tilted by 58$^{\circ}$ $\pm$ 20$^{\circ}$ into the plane of the sky away from the observer, roughly opposite to the inferred direction of motion of the central compact object. We find some $^{44}$Ti ejecta that are clearly interior to the reverse shock and some that are clearly exterior to the reverse shock. Where we observe $^{44}$Ti ejecta exterior to the reverse shock we also see shock-heated iron; however, there are regions where we see iron but do not observe $^{44}$Ti. This suggests that the local conditions of the supernova shock during explosive nucleosynthesis varied enough to suppress the production of $^{44}$Ti in some regions by at least a factor of two, even in regions that are assumed to be the result of processes like $α$-rich freezeout that should produce both iron and titanium.

preprint2016arXiv

The Geometry of the Infrared and X-ray Obscurer in a Dusty Hyperluminous Quasar

We study the geometry of the AGN obscurer in IRAS 09104+4109, an IR-luminous, radio-intermediate FR-I source at $z=0.442$, using infrared data from Spitzer and Herschel, X-ray data from Nustar, Swift, Suzaku, and Chandra, and an optical spectrum from Palomar. The infrared data imply a total rest-frame 1-1000$μ$m luminosity of $5.5\times10^{46}$erg s$^{-1}$ and require both an AGN torus and starburst model. The AGN torus has an anisotropy-corrected IR luminosity of $4.9\times10^{46}$erg s$^{-1}$, and a viewing angle and half opening angle both of approximately $36$ degrees from pole-on. The starburst has a star formation rate of $(110\pm34)$M$_{\odot}$ yr$^{-1}$ and an age of $<50$Myr. These results are consistent with two epochs of luminous activity in 09104: one approximately $150$Myr ago, and one ongoing. The X-ray data suggest a photon index of $Γ\simeq 1.8$ and a line-of-sight column of $N_{\rm H} \simeq 5\times10^{23}$cm$^{-2}$. This argues against a reflection-dominated hard X-ray spectrum, which would have implied a much higher $N_{\rm H}$ and luminosity. The X-ray and infrared data are consistent with a bolometric AGN luminosity of $L_{\rm bol}\sim(0.5-2.5)\times10^{47}$erg s$^{-1}$. The X-ray and infrared data are further consistent with coaligned AGN obscurers in which the line of sight "skims" the torus. This is also consistent with the optical spectra, which show both coronal iron lines and broad lines in polarized but not direct light. Combining constraints from the X-ray, optical, and infrared data suggests that the AGN obscurer is within a vertical height of $20$pc, and a radius of $125$pc, of the nucleus.

preprint2015arXiv

A Connection between Obscuration and Star Formation in Luminous Quasars

We present a measurement of the star formation properties of a uniform sample of mid-IR selected, unobscured and obscured quasars (QSO1s and QSO2s) in the Boötes survey region. We use an spectral energy distribution (SED) analysis for photometric data spanning optical to far-IR wavelengths to decompose AGN and host galaxy components. We find that when compared to a matched sample of QSO1s, the QSO2s have higher far-IR detection fractions, far-IR fluxes and infrared star formation luminosities ($L_{\rm IR}^{\rm SF}$) by a factor of $\sim2$. Correspondingly, we show that the AGN obscured fraction rises from 0.3 to 0.7 between $4-40\times10^{11}L_\odot$. We also find evidence associating the absorption in the X-ray emission with the presence of far-IR emitting dust. Overall, these results are consistent with galaxy evolution models in which quasar obscurations can be associated with a dust-enshrouded starburst galaxies.

preprint2015arXiv

Evolution in the Black Hole - Galaxy Scaling Relations and the Duty Cycle of Nuclear Activity in Star-Forming Galaxies

We measure the location and evolutionary vectors of 69 Herschel-detected broad-line active galactic nuclei (BLAGNs) in the M_BH-M_* plane. BLAGNs are selected from the COSMOS and CDF-S fields, and span the redshift range 0.2< z<2.1. Black-hole masses are calculated using archival spectroscopy and single-epoch virial mass estimators, and galaxy total stellar masses are calculated by fitting the spectral energy distribution (subtracting the BLAGN component). The mass-growth rates of both the black hole and galaxy are calculated using Chandra/XMM-Newton X-ray and Herschel far-infrared data, reliable measures of the BLAGN accretion and galaxy star formation rates, respectively. We use Monte Carlo simulations to account for biases in our sample, due to both selection limits and the steep slope of the massive end of the galaxy stellar-mass distribution. We find our sample is consistent with no evolution in the M_BH-M_* relation from z~2 to z~0. BLAGNs and their host galaxies which lie off the black hole mass galaxy total stellar mass relation tend to have evolutionary vectors anti-correlated with their mass ratios: that is, galaxies with over-massive (under-massive) black holes tend to have a low (high) ratio of the specific accretion rate to the specific star formation rate. We also use the measured growth rates to estimate the preferred AGN duty cycle for our galaxies to evolve most consistently onto the local M_BH-M_Bul relation. Under reasonable assumptions of exponentially declining star formation histories, the data suggest a non-evolving (no more than a factor of a few) BLAGN duty cycle among star-forming galaxies of 10% (1sigma range of 1-42% at z<1 and 2-34% at z>1).

preprint2015arXiv

NuSTAR Hard X-ray Survey of the Galactic Center Region I: Hard X-ray Morphology and Spectroscopy of the Diffuse Emission

We present the first sub-arcminute images of the Galactic Center above 10 keV, obtained with NuSTAR. NuSTAR resolves the hard X-ray source IGR J17456-2901 into non-thermal X-ray filaments, molecular clouds, point sources and a previously unknown central component of hard X-ray emission (CHXE). NuSTAR detects four non-thermal X-ray filaments, extending the detection of their power-law spectra with $Γ\sim1.3$-$2.3$ up to ~50 keV. A morphological and spectral study of the filaments suggests that their origin may be heterogeneous, where previous studies suggested a common origin in young pulsar wind nebulae (PWNe). NuSTAR detects non-thermal X-ray continuum emission spatially correlated with the 6.4 keV Fe K$α$ fluorescence line emission associated with two Sgr A molecular clouds: MC1 and the Bridge. Broad-band X-ray spectral analysis with a Monte-Carlo based X-ray reflection model self-consistently determined their intrinsic column density ($\sim10^{23}$ cm$^{-2}$), primary X-ray spectra (power-laws with $Γ\sim2$) and set a lower limit of the X-ray luminosity of Sgr A* flare illuminating the Sgr A clouds to $L_X \stackrel{>}{\sim} 10^{38}$ erg s$^{-1}$. Above ~20 keV, hard X-ray emission in the central 10 pc region around Sgr A* consists of the candidate PWN G359.95-0.04 and the CHXE, possibly resulting from an unresolved population of massive CVs with white dwarf masses $M_{\rm WD} \sim 0.9 M_{\odot}$. Spectral energy distribution analysis suggests that G359.95-0.04 is likely the hard X-ray counterpart of the ultra-high gamma-ray source HESS J1745-290, strongly favoring a leptonic origin of the GC TeV emission.

preprint2015arXiv

Variable Hard X-ray Emission from the Candidate Accreting Black Hole in Dwarf Galaxy Henize 2-10

We present an analysis of the X-ray spectrum and long-term variability of the nearby dwarf starburst galaxy Henize 2-10. Recent observations suggest that this galaxy hosts an actively accreting black hole with mass ~10^6 M_sun. The presence of an AGN in a low-mass starburst galaxy marks a new environment for active galactic nuclei (AGNs), with implications for the processes by which "seed" black holes may form in the early Universe. In this paper, we analyze four epochs of X-ray observations of Henize 2-10, to characterize the long-term behavior of its hard nuclear emission. We analyze observations with Chandra from 2001 and XMM-Newton from 2004 and 2011, as well as an earlier, less sensitive observation with ASCA from 1997. Based on detailed analysis of the source and background, we find that the hard (2-10 keV) flux of the putative AGN has decreased by approximately an order of magnitude between the 2001 Chandra observation and exposures with XMM-Newton in 2004 and 2011. The observed variability confirms that the emission is due to a single source. It is unlikely that the variable flux is due to a supernova or ultraluminous X-ray source, based on the observed long-term behavior of the X-ray and radio emission, while the observed X-ray variability is consistent with the behavior of well-studied AGNs.

preprint2014arXiv

NuSTAR Spectroscopy of Multi-Component X-ray Reflection from NGC 1068

We report on observations of NGC1068 with NuSTAR, which provide the best constraints to date on its $>10$~keV spectral shape. We find no strong variability over the past two decades, consistent with its Compton-thick AGN classification. The combined NuSTAR, Chandra, XMM-Newton, and Swift-BAT spectral dataset offers new insights into the complex reflected emission. The critical combination of the high signal-to-noise NuSTAR data and a spatial decomposition with Chandra allow us to break several model degeneracies and greatly aid physical interpretation. When modeled as a monolithic (i.e., a single N_H) reflector, none of the common Compton-reflection models are able to match the neutral fluorescence lines and broad spectral shape of the Compton reflection. A multi-component reflector with three distinct column densities (e.g., N_H~1.5e23, 5e24, and 1e25 cm^{-2}) provides a more reasonable fit to the spectral lines and Compton hump, with near-solar Fe abundances. In this model, the higher N_H components provide the bulk of the Compton hump flux while the lower N_H component produces much of the line emission, effectively decoupling two key features of Compton reflection. We note that ~30% of the neutral Fe Kalpha line flux arises from >2" (~140 pc), implying that a significant fraction of the <10 keV reflected component arises from regions well outside of a parsec-scale torus. These results likely have ramifications for the interpretation of poorer signal-to-noise observations and/or more distant objects [Abridged].

preprint2014arXiv

The variable hard X-ray emission of NGC4945 as observed by NuSTAR

We present a broadband (~0.5-79 keV) spectral and temporal analysis of multiple NuSTAR observations combined with archival Suzaku and Chandra data of NGC4945, the brightest extragalactic source at 100 keV. We observe hard X-ray (> 10 keV) flux and spectral variability, with flux variations of a factor 2 on timescales of 20 ksec. A variable primary continuum dominates the high energy spectrum (>10 keV) in all the states, while the reflected/scattered flux which dominates at E< 10 keV stays approximately constant. From modelling the complex reflection/transmission spectrum we derive a Compton depth along the line of sight of tau_Thomson ~ 2.9, and a global covering factor for the circumnuclear gas of ~ 0.15. This agrees with the constraints derived from the high energy variability, which implies that most of the high energy flux is transmitted, rather that Compton-scattered. This demonstrates the effectiveness of spectral analysis in constraining the geometric properties of the circumnuclear gas, and validates similar methods used for analyzing the spectra of other bright, Compton-thick AGN. The lower limits on the e-folding energy are between 200-300 keV, consistent with previous BeppoSAX, Suzaku and Swift BAT observations. The accretion rate, estimated from the X-ray luminosity and assuming a bolometric correction typical of type 2 AGN, is in the range ~0.1-0.3 lambda_Edd depending on the flux state. The substantial observed X-ray luminosity variability of NGC4945 implies that large errors can arise from using single-epoch X-ray data to derive L/L_Edd values for obscured AGNs.

preprint2013arXiv

A comparative analysis of virial black-hole mass estimates of moderate-luminosity active galactic nuclei using Subaru/FMOS

We present an analysis of broad emission lines observed in moderate-luminosity active galactic nuclei (AGNs), typical of those found in X-ray surveys of deep fields, with the aim to test the validity of single-epoch virial black hole mass estimates. We have acquired near-infrared (NIR) spectra of AGNs up to z ~ 1.8 in the COSMOS and Extended Chandra Deep Field-South Survey, with the Fiber Multi-Object Spectrograph (FMOS) mounted on the Subaru Telescope. These low-resolution NIR spectra provide a significant detection of the broad Halpha line that has been shown to be a reliable probe of black hole mass at low redshift. Our sample has existing optical spectroscopy which provides a detection of MgII, a broad emission line typically used for black hole mass estimation at z > 1. We carry out a spectral-line fitting procedure using both Halpha and MgII to determine the virial velocity of gas in the broad line region, the monochromatic continuum luminosity at 3000 A, and the total Halpha line luminosity. With a sample of 43 AGNs spanning a range of two decades in luminosity (i.e., L ~ 10^44-46 ergs/s), we find a tight correlation between the continuum and line luminosity with a distribution characterized by <log(L_3000/L_Halpha)> = 1.52 and a dispersion sigma = 0.16. There is also a close one-to-one relationship between the FWHM of Halpha and of MgII up to 10000 km/s with a dispersion of 0.14 in the distribution of the logarithm of their ratios. Both of these then lead to there being very good agreement between Halpha- and MgII-based masses over a wide range in black hole mass (i.e., M_BH ~ 10^7-9 M_sun). We do find a small offset in MgII-based masses, relative to those based on Halpha, of +0.17 dex and a dispersion sigma = 0.32. In general, these results demonstrate that local scaling relations, using MgII or Halpha, are applicable for AGN at moderate luminosities and up to z ~ 2.

preprint2013arXiv

A correlation between star formation rate and average black hole accretion in star forming galaxies

We present a measurement of the average supermassive black hole accretion rate (BHAR) as a function of star formation rate (SFR) for galaxies in the redshift range 0.25<z<0.8. We study a sample of 1,767 far-IR selected star-forming galaxies in the 9 deg^2 Boötes multiwavelength survey field. The SFR is estimated using 250 micron observations from the Herschel Space Observatory, for which the contribution from the AGN is minimal. In this sample, 121 AGNs are directly identified using X-ray or mid-IR selection criteria. We combined these detected AGNs and an X-ray stacking analysis for undetected sources to study the average BHAR for all of the star-forming galaxies in our sample. We find an almost linear relation between the average BHAR (in M_sun/year) and the SFR (in M_sun/year) for galaxies across a wide SFR range 0.85<log SFR<2.56 : log BHAR=(-3.72\pm0.52)+(1.05\pm0.33) log SFR. This global correlation between SFR and average BHAR is consistent with a simple picture in which SFR and AGN activity are tightly linked over galaxy evolution timescales.

preprint2013arXiv

Black hole variability and the star formation-AGN connection: Do all star-forming galaxies host an AGN?

We investigate the effect of active galactic nucleus (AGN) variability on the observed connection between star formation and black hole accretion in extragalactic surveys. Recent studies have reported relatively weak correlations between observed AGN luminosities and the properties of AGN hosts, which has been interpreted to imply that there is no direct connection between AGN activity and star formation. However, AGNs may be expected to vary significantly on a wide range of timescales (from hours to Myr) that are far shorter than the typical timescale for star formation (>~100 Myr). This variability can have important consequences for observed correlations. We present a simple model in which all star-forming galaxies host an AGN when averaged over ~100 Myr timescales, with long-term average AGN accretion rates that are perfectly correlated with the star formation rate (SFR). We show that reasonable prescriptions for AGN variability reproduce the observed weak correlations between SFR and L_AGN in typical AGN host galaxies, as well as the general trends in the observed AGN luminosity functions, merger fractions, and measurements of the average AGN luminosity as a function of SFR. These results imply there may be a tight connection between AGN activity and SFR over galaxy evolution timescales, and that the apparent similarities in rest-frame colors, merger rates, and clustering of AGNs compared to "inactive" galaxies may be due primarily to AGN variability. The results provide motivation for future deep, wide extragalactic surveys that can measure the distribution of AGN accretion rates as a function of SFR.

preprint2013arXiv

Modeling the Nuclear Infrared Spectral Energy Distribution of Type II Active Galactic Nuclei

We present results from model fitting to the Spectral Energy Distribution (SED) of a homogeneous sample of Seyfert II galaxies drawn from the $12μ$m Galaxy Sample. Imaging and nuclear flux measurements are presented in an accompanying paper (Videla et al., 2013). Here we add IRS Spitzer observations to further constrain the SEDs after careful subtraction of a starburst component. We use the library of CLUMPY torus models from Nenkova et al.~(2008ab) and also test the two-phase models recently produced by Stalevski et al.~(2012). We find that photometric and spectroscopic observations in the mid-IR (>5mu) are crucial to properly constrain the best-fit torus models. About half of our sources show clear near-IR excess of their SEDs above the best fit models. This problem can be less severe when using the Stalevski et al.~(2012) models. It is not clear what is the nature of this emission since best fitted black body temperatures are very high (~1700-2500 K) and the Type II classification of our sources would correspond to a small probability to peer directly into the hottest regions of the torus. Crucially, the derived torus parameters when using CLUMPY models are quite robust,, independently of whether the sources require an additional black body component or not. Our findings suggest that tori are characterized by N_0 > 5, sigma > 40, tau < 25, i > 40 degrees, Y < 50 and A_v^los ~ 100-300. From these we can determine that typical torus sizes and masses of 0.1-5.0 pc and 10^{4-6} M_\odot. We find tentative evidence that those nuclei with a detected Hidden Broad Line Regions are characterized by lower levels of extinction than those without one. Finally, we find no correlation between the torus properties and the presence of circumnuclear or more global star-formation.

preprint2013arXiv

Nuclear Infrared Spectral Energy Distribution of Type II Active Galactic Nuclei

We present near and mid--IR observations of a sample of Seyfert II galaxies drawn from the $12μ$m Galaxy Sample. The sample was observed in the J, H, K, L, M and N bands. Galaxy Surface Brightness Profiles are modeled using nuclear, bulge, bar (when necessary) and disk components. To check the reliability of our findings the procedure was tested using {\em Spitzer\/} observations of M\,31. Nuclear Spectral Energy Distributions (SEDs) are determined for 34 objects, and optical spectra are presented for 38, including analysis of their stellar populations using the STARLIGHT spectral synthesis code. Emission line diagnostic-diagrams are used to discriminate between genuine AGN and HII nuclei. Combining our observations with those found in the literature, we have a total of 40 SEDs. It is found that about 40\%\ of the SEDs are characterized by an upturn in the near-IR, which we have quantified as a NIR slope $α< 1$ for an SED characterized as $λf_λ \propto λ^α$. Three objects with an HII nucleus and two Seyfert nuclei with strong contamination from a circumnuclear starburst, also show an upturn. For genuine AGN this component could be explained as emission from the accretion disk, a jet, or from a very hot dust component leaking from the central region through a clumpy obscuring structure. The presence of a very compact nuclear starburst as the origin for this NIR excess emission is not favored by our spectroscopic data for these objects.

preprint2013arXiv

The Nuclear Spectroscopic Telescope Array (NuSTAR) Mission

The Nuclear Spectroscopic Telescope Array (NuSTAR) mission, launched on 13 June 2012, is the first focusing high-energy X-ray telescope in orbit. NuSTAR operates in the band from 3 -- 79 keV, extending the sensitivity of focusing far beyond the ~10 keV high-energy cutoff achieved by all previous X-ray satellites. The inherently low-background associated with concentrating the X-ray light enables NuSTAR to probe the hard X-ray sky with a more than one-hundred-fold improvement in sensitivity over the collimated or coded-mask instruments that have operated in this bandpass. Using its unprecedented combination of sensitivity, spatial and spectral resolution, NuSTAR will pursue five primary scientific objectives, and will also undertake a broad program of targeted observations. The observatory consists of two co-aligned grazing-incidence X-ray telescopes pointed at celestial targets by a three-axis stabilized spacecraft. Deployed into a 600 km, near-circular, 6degree inclination orbit, the Observatory has now completed commissioning, and is performing consistent with pre-launch expectations. NuSTAR is now executing its primary science mission, and with an expected orbit lifetime of ten years, we anticipate proposing a guest investigator program, to begin in Fall 2014.

preprint2013arXiv

Widespread and Hidden Active Galactic Nuclei in Star-Forming Galaxies at Redshift > 0.3

We characterize the incidence of active galactic nuclei (AGNs) is 0.3 < z < 1 star-forming galaxies by applying multi-wavelength AGN diagnostics (X-ray, optical, mid-infrared, radio) to a sample of galaxies selected at 70-micron from the Far-Infrared Deep Extragalactic Legacy survey (FIDEL). Given the depth of FIDEL, we detect "normal" galaxies on the specific star formation rate (sSFR) sequence as well as starbursting systems with elevated sSFR. We find an overall high occurrence of AGN of 37+/-3%, more than twice as high as in previous studies of galaxies with comparable infrared luminosities and redshifts but in good agreement with the AGN fraction of nearby (0.05 < z < 0.1) galaxies of similar infrared luminosities. The more complete census of AGNs comes from using the recently developed Mass-Excitation (MEx) diagnostic diagram. This optical diagnostic is also sensitive to X-ray weak AGNs and X-ray absorbed AGNs, and reveals that absorbed active nuclei reside almost exclusively in infrared-luminous hosts. The fraction of galaxies hosting an AGN appears to be independent of sSFR and remains elevated both on the sSFR sequence and above. In contrast, the fraction of AGNs that are X-ray absorbed increases substantially with increasing sSFR, possibly due to an increased gas fraction and/or gas density in the host galaxies.

preprint2012arXiv

GOODS-Herschel and CANDELS: The Morphologies of Ultraluminous Infrared Galaxies at z~2

Using deep 100-160 micron observations in GOODS-S from the GOODS-H survey, combined with HST/WFC3 NIR imaging from CANDELS, we present the first morphological analysis of a complete, FIR selected sample of 52 ULIRGs at z~2. We also make use of a comparison sample of galaxies without Herschel detections but with the same z and magnitude distribution. Our visual classifications of these two samples indicate that the fraction of objects with disk and spheroid morphologies is roughly the same but that there are significantly more mergers, interactions, and irregular galaxies among the ULIRGs. The combination of disk and irregular/interacting morphologies suggests that early stage interactions and minor mergers could play an important role in ULIRGs at z~2. We compare these fractions with those of a z~1 sample across a wide luminosity range and find that the fraction of disks decreases systematically with L_IR while the fraction of mergers and interactions increases, as has been observed locally. At comparable luminosities, the fraction of ULIRGs with various morphological classifications is similar at z~2 and z~1. We investigate the position of the ULIRGs, along with 70 LIRGs, on the specific star formation rate versus redshift plane, and find 52 systems to be starbursts (lie more than a factor of 3 above the main sequence relation). The morphologies of starbursts are dominated by interacting and merging systems (50%). If irregular disks are included as potential minor mergers, then we find that up to 73% of starbursts are involved in a merger or interaction at some level. Although the final coalescence of a major merger may not be required for the high luminosities of ULIRGs at z~2 as is the case locally, the large fraction of interactions at all stages and potential minor mergers suggest that the high star formation rates of ULIRGs are still largely externally triggered at z~2.

preprint2012arXiv

GOODS-Herschel: Impact of Active Galactic Nuclei and Star Formation Activity on Infrared Spectral Energy Distributions at High Redshift

We explore the effects of active galactic nuclei (AGN) and star formation activity on the infrared (0.3-1000 microns) spectral energy distributions of luminous infrared galaxies from z = 0.5 to 4.0. We have compiled a large sample of 151 galaxies selected at 24 microns (S24 > 100 uJy) in the GOODS-N and ECDFS fields for which we have deep Spitzer IRS spectroscopy, allowing us to decompose the mid-IR spectrum into contributions from star formation and AGN activity. A significant portion (~25%) of our sample is dominated by an AGN in the mid-IR. Based on the mid-IR classification, we divide our full sample into four sub-samples: z~1 star-forming (SF) sources; z~2 SF sources; AGN with clear 9.7 micron silicate absorption; and AGN with featureless mid-IR spectra. From our large spectroscopic sample and wealth of multi-wavelength data, including deep Herschel imaging at 100, 160, 250, 350, and 500 microns, we use 95 galaxies with complete spectral coverage to create a composite spectral energy distribution (SED) for each sub-sample. We then fit a two-temperature component modified blackbody to the SEDs. We find that the IR SEDs have similar cold dust temperatures, regardless of the mid-IR power source, but display a marked difference in the warmer dust temperatures. We calculate the average effective temperature of the dust in each sub-sample and find a significant (~20 K) difference between the SF and AGN systems. We compare our composite SEDs to local templates and find that local templates do not accurately reproduce the mid-IR features and dust temperatures of our high redshift systems. High redshift IR luminous galaxies contain significantly more cool dust than their local counterparts. We find that a full suite of photometry spanning the IR peak is necessary to accurately account for the dominant dust temperature components in high redshift IR luminous galaxies.

preprint2012arXiv

What Drives the Growth of Black Holes?

Massive black holes (BHs) are at once exotic and yet ubiquitous, residing in the centers of massive galaxies in the local Universe. Recent years have seen remarkable advances in our understanding of how these BHs form and grow over cosmic time, during which they are revealed as active galactic nuclei (AGN). However, despite decades of research, we still lack a coherent picture of the physical drivers of BH growth, the connection between the growth of BHs and their host galaxies, the role of large-scale environment on the fueling of BHs, and the impact of BH-driven outflows on the growth of galaxies. In this paper we review our progress in addressing these key issues, motivated by the science presented at the "What Drives the Growth of Black Holes?" workshop held at Durham on 26th-29th July 2010, and discuss how these questions may be tackled with current and future facilities.

preprint2011arXiv

A New Diagnostic of Active Galactic Nuclei: Revealing Highly-Absorbed Systems at Redshift>0.3

We introduce the Mass-Excitation (MEx) diagnostic to identify active galactic nuclei (AGN) in galaxies at intermediate redshift. In the absence of near-infrared spectroscopy, necessary to use traditional nebular line diagrams at z>0.4, we demonstrate that combining [OIII]5007/Hbeta and stellar mass successfully distinguishes between star formation and AGN emission. The MEx classification scheme relies on a novel probabilistic approach splitting galaxies into sub-categories with more confidence than alternative high-z diagnostic diagrams. It recognizes that galaxies near empirical boundaries on traditional diagrams have an uncertain classification and thus a non-zero probability of belonging to more than one category. An outcome of this work is a system of statistical weights that can be used to compute global properties of galaxy samples. We apply the MEx diagram to 2,812 galaxies at 0.3<z<1 in the Great Observatories Origins Deep Survey North and Extended Groth Strip fields, and compare it to an independent X-ray classification scheme. We identify Compton-thick AGN candidates with large X-ray absorption, which we infer from the luminosity ratio between hard X-ray emission and [OIII]5007, a nearly isotropic tracer of AGN. X-ray stacking of sources that were not detected individually supports the validity of the MEx diagram and yields a very flat spectral slope for the Compton-thick candidates (Gamma~0.4, unambiguously indicating absorbed AGN). We present evidence that composite galaxies, which are difficult to identify with alternative high-redshift diagrams, host the majority of the highly-absorbed AGN. Our findings suggest that the interstellar medium of the host galaxy provides significant absorption in addition to the torus invoked in AGN unified models.

preprint2011arXiv

CANDELS: The Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey

The Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) is designed to document the first third of galactic evolution, over the approximate redshift (z) range 8--1.5. It will image >250,000 distant galaxies using three separate cameras on the Hubble Space Telescope, from the mid-ultraviolet to the near-infrared, and will find and measure Type Ia supernovae at z>1.5 to test their accuracy as standardizable candles for cosmology. Five premier multi-wavelength sky regions are selected, each with extensive ancillary data. The use of five widely separated fields mitigates cosmic variance and yields statistically robust and complete samples of galaxies down to a stellar mass of 10^9 M_\odot to z \approx 2, reaching the knee of the ultraviolet luminosity function (UVLF) of galaxies to z \approx 8. The survey covers approximately 800 arcmin^2 and is divided into two parts. The CANDELS/Deep survey (5σ point-source limit H=27.7 mag) covers \sim 125 arcmin^2 within GOODS-N and GOODS-S. The CANDELS/Wide survey includes GOODS and three additional fields (EGS, COSMOS, and UDS) and covers the full area to a 5σ point-source limit of H \gtrsim 27.0 mag. Together with the Hubble Ultra Deep Fields, the strategy creates a three-tiered "wedding cake" approach that has proven efficient for extragalactic surveys. Data from the survey are nonproprietary and are useful for a wide variety of science investigations. In this paper, we describe the basic motivations for the survey, the CANDELS team science goals and the resulting observational requirements, the field selection and geometry, and the observing design. The Hubble data processing and products are described in a companion paper.

preprint2011arXiv

Clustering of obscured and unobscured quasars in the Bootes field: Placing rapidly growing black holes in the cosmic web

We present the first measurement of the spatial clustering of mid-infrared selected obscured and unobscured quasars, using a sample in the redshift range 0.7 < z < 1.8 selected from the 9 deg^2 Bootes multiwavelength survey. Recently the Spitzer Space Telescope and X-ray observations have revealed large populations of obscured quasars that have been inferred from models of the X-ray background and supermassive black hole evolution. To date, little is known about obscured quasar clustering, which allows us to measure the masses of their host dark matter halos and explore their role in the cosmic evolution of black holes and galaxies. In this study we use a sample of 806 mid-infrared selected quasars and ~250,000 galaxies to calculate the projected quasar-galaxy cross-correlation function w_p(R). The observed clustering yields characteristic dark matter halo masses of log (M_halo [h^-1 M_sun]) = 12.7^+0.4_-0.6 and 13.3^+0.3_-0.4 for unobscured quasars (QSO-1s) and obscured quasars (Obs-QSOs), respectively. The results for QSO-1s are in excellent agreement with previous measurements for optically-selected quasars, while we conclude that the Obs-QSOs are at least as strongly clustered as the QSO-1s. We test for the effects of photometric redshift errors on the optically-faint Obs-QSOs, and find that our method yields a robust lower limit on the clustering; photo-z errors may cause us to underestimate the clustering amplitude of the Obs-QSOs by at most ~20%. We compare our results to previous studies, and speculate on physical implications of stronger clustering for obscured quasars.

preprint2010arXiv

A multi-wavelength approach to the properties of Extremely Red Galaxy populations: I - Contribution to the Star Formation Rate density and AGN content

We present a multi-wavelength analysis of the properties of Extremely Red Galaxy (ERG) populations, selected in the GOODS-South/Chandra Deep Field South field. By using all the photometric and spectroscopic information available on large deep samples of EROs (645 sources), IEROs (294 sources), and DRGs (350 sources), we derive redshift distributions, identify AGN powered and Star-formation powered galaxies, and, using the radio observations of this field, estimate robust (AGN- and dust-unbiased) Star Formation Rate Densities (SFRD) for these populations. We also investigate the properties of "pure" (galaxies that conform to only one of the three ERG criteria considered) and "combined" (galaxies that verify all three criteria) sub-populations. Overall, a large number of AGN are identified (up to ~30%, based on X-rays, and mid-infrared criteria), the majority of which are type-2 (obscured) objects. Among ERGs with no evidence for AGN activity, we identify sub-populations covering a wide range of average star-formation rates, from below 10 Mo/yr to as high as 200 Mo/yr. Applying a redshift separation (1<z<2 and 2<z<3) we find significant evolution (an increase of a factor of 2 or higher) of SFRD for EROs and DRGs, while none is observed for IEROs. The former populations can contribute more than 20\% to the global SFRD at 2<z<3. The emission from AGN activity is typically not strong in the ERG population, with AGN increasing the average radio luminosity of ERG sub-populations by, nominally, less than 20%. AGN are common, however, and, if no discrimination is attempted, this could significantly increase the SFRD estimate (by over 100% in some cases). Thus, and while the contribution of star forming processes to the radio luminosity in galaxies with AGN remains uncertain, a comprehensive identification of AGN in these populations is necessary to obtain meaningful results.

preprint2010arXiv

Detection of molecular gas in a distant submillimetre galaxy at z=4.76 with ATCA

We have detected the CO(2-1) transition from the submillimetre galaxy (SMG) LESSJ033229.4-275619 at z=4.755 using the new Compact Array Broadband Backend system on the Australian Telescope Compact Array. These data have identified a massive gas reservoir available for star formation for the first time in an SMG at z~5. We use the luminosity and velocity width (FWHM of 160 km/s) of the CO(2--1) line emission to constrain the gas and dynamical mass of Mgas~1.6x10^10 Msun and Mdyn(<2kpc)~5x10^10 (0.25/sin^2(i)) Msun, respectively, similar to that observed for SMGs at lower redshifts of z~2-4, although we note that our observed CO FWHM is a factor of ~3 narrower than typically seen in SMGs. Together with the stellar mass we estimate a total baryonic mass of Mbary~1x10^11 Msun, consistent with the dynamical mass for this young galaxy within the uncertainties. Dynamical and baryonic mass limits of high-redshift galaxies are useful tests of galaxy formation models: using the known z~4-5 SMGs as examples of massive baryonic systems, we find that their space density is consistent with that predicted by current galaxy formation models. In addition, these observations have helped to confirm that z~4-5 SMGs possess the baryonic masses and gas consumption timescales necessary to be the progenitors of the luminous old red galaxies seen at z~3. Our results provide a preview of the science that ALMA will enable on the formation and evolution of the earliest massive galaxies in the Universe.

preprint2010arXiv

MOIRCS Deep Survey. VI. Near-Infrared Spectroscopy of K-Selected Star-Forming Galaxies at z~2

We present near-infrared (NIR) spectroscopy for 37 BzK-color-selected star-forming galaxies with MOIRCS on Subaru. The sample is drawn from the Ks-band selected catalog of the MOIRCS Deep Survey (MODS) in the GOODS-N region. About half of our samples are selected from the MIPS 24um-source catalog. Ha emission lines are detected from 23 galaxies, of which the median redshift is 2.12. We derived the star formation rates (SFRs) from extinction-corrected Ha luminosities. The extinction correction is estimated from the SED fitting of multi-band photometric data covering UV to NIR wavelengths. The Balmer decrement of the stacked emission lines shows that the amount of extinction for the ionized gas is larger than that for the stellar continuum. From a comparison of the extinction corrected Ha luminosity and other SFR indicators we found that the relation between the dust properties of stellar continuum and ionized gas is different depending on the intrinsic SFR. The comparison between the Ha SFR and stellar mass estimated from SED fitting shows no correlation between them. Some galaxies with stellar mass smaller than ~10^10 Msun show SFRs higher than ~100 Msun/yr. The specific SFRs (SSFRs) of these galaxies are remarkably high. The average stellar-population age of these high-SSFR galaxies is younger than 100 Myr, which is consistent with the implied high SSFR. The large SFR implies the possibility that the high SSFR galaxies significantly contribute to the cosmic SFR density of the universe at z~2. The total SFR density estimated from the Ha emission line galaxies is 0.089-0.136 Msun/yr/Mpc^3, which is consistent with the total SFR densities in the literature. The metallicity of the high-SSFR galaxies is larger than that expected from the mass-metallicity relation of UV-selected galaxies at z~2 by Erb et al. (2006a). [Abridged]

preprint2010arXiv

Starburst or AGN Dominance in Submillimetre-Luminous Candidate AGN?

It is widely believed that ultraluminous infrared (IR) galaxies and active galactic nuclei (AGN) activity are triggered by galaxy interactions and merging, with the peak of activity occurring at z~2, where submillimetre galaxies are thousands of times more numerous than local ULIRGs. In this evolutionary picture, submillimetre galaxies (SMGs) would host an AGN, which would eventually grow a black hole (BH) strong enough to blow off all of the gas and dust leaving an optically luminous QSO. To probe this evolutionary sequence we have focussed on the 'missing link' sources, which demonstrate both strong starburst (SB) and AGN signatures, in order to determine if the SB is the main power source even in SMGs when we have evidence that an AGN is present from their IRAC colours. The best way to determine if a dominant AGN is present is to look for their signatures in the mid-infrared with the Spitzer IRS, since often even deep X-ray observations miss identifying the presence of AGN in heavily dust-obscured SMGs. We present the results of our audit of the energy balance between star-formation and AGN within this special sub-population of SMGs -- where the BH has grown appreciably to begin heating the dust emission.

preprint2010arXiv

Towards a Complete Census of AGNs in Nearby Galaxies: The Incidence of Growing Black Holes

We investigate the local supermassive black hole (SMBH) density function and relative mass accretion rates of all active galactic nuclei (AGNs) identified in a volume-limited sample of infrared (IR) bright galaxies (L_IR > 3 x 10^9 L_sun) to D<15 Mpc (Goulding & Alexander 2009). A database of accurate SMBH mass (M_BH) estimates is compiled from literature sources using physically motivated AGN modeling techniques (reverberation mapping, maser mapping and gas kinematics) and well-established indirect M_BH estimation methods (the M-sigma and M_BH-L_(K,bul) relations). For the three sources without previously published M_BH estimates, we use 2MASS K-band imaging and GALFIT to constrain the bulge luminosities, and hence SMBH masses. In general, we find the AGNs in the sample host SMBHs which are spread over a wide mass range (M_BH ~ (0.1-30) x 10^7 M_sun), but with the majority in the poorly studied M_BH ~ 10^6-10^7 M_sun region. Using sensitive hard X-ray (2-10 keV) and mid-IR constraints we calculate the bolometric luminosities of the AGNs (L_(Bol,AGN)) and use them to estimate relative mass accretion rates. We use these data to calculate the volume-average SMBH growth rate of galaxies in the local Universe and find that the AGNs hosting SMBHs in the mass range M_BH ~ 10^6-10^7 M_sun are dominated by optically unidentified AGNs. These relatively small SMBHs are acquiring a significant proportion of their mass in the present-day, and are amongst the most rapidly growing in the local Universe (SMBH mass doubling times of ~6 Gyrs). Additionally, we find tentative evidence for an increasing volume-weighted AGN fraction with decreasing SMBH mass in the M_BH ~ 10^6-10^8 M_sun range. Overall, we conclude that significant mass accretion onto small SMBHs may be missed in even the most sensitive optical surveys due to absent or weak optical AGN signatures.

preprint2009arXiv

BLAST: the far-infrared/radio correlation in distant galaxies

We investigate the correlation between far-infrared (FIR) and radio luminosities in distant galaxies, a lynchpin of modern astronomy. We use data from the Balloon-borne Large Aperture Submillimetre Telescope (BLAST), Spitzer, the Large Apex BOlometer CamerA (LABOCA), the Very Large Array (VLA) and the Giant Metre-wave Radio Telescope (GMRT) in the Extended Chandra Deep Field South (ECDFS). For a catalogue of BLAST 250-micron-selected galaxies, we re-measure the 70--870-micron flux densities at the positions of their most likely 24-micron counterparts, which have a median [interquartile] redshift of 0.74 [0.25, 1.57]. From these, we determine the monochromatic flux density ratio, q_250 = log_10 (S_250micron / S_1400MHz), and the bolometric equivalent, q_IR. At z ~= 0.6, where our 250-micron filter probes rest-frame 160-micron emission, we find no evolution relative to q_160 for local galaxies. We also stack the FIR and submm images at the positions of 24-micron- and radio-selected galaxies. The difference between q_IR seen for 250-micron- and radio-selected galaxies suggests star formation provides most of the IR luminosity in ~< 100-uJy radio galaxies, but rather less for those in the mJy regime. For the 24-micron sample, the radio spectral index is constant across 0 < z < 3, but q_IR exhibits tentative evidence of a steady decline such that q_IR is proportional to (1+z)^(-0.15 +/- 0.03) - significant evolution, spanning the epoch of galaxy formation, with major implications for techniques that rely on the FIR/radio correlation. We compare with model predictions and speculate that we may be seeing the increase in radio activity that gives rise to the radio background.

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.