Source author record

D. M. Alexander

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

ResearcherUnclaimed source record

astro-ph.CO astro-ph.GA astro-ph.HE astro-ph

Catalog footprint

What is connected

140works

4topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2022arXiv

GOODS-ALMA 2.0: Starbursts in the main sequence reveal compact star formation regulating galaxy evolution prequenching

Compact star formation appears to be generally common in dusty star-forming galaxies (SFGs). However, its role in the framework set by the scaling relations in galaxy evolution remains to be understood. In this work we follow up on the galaxy sample from the GOODS-ALMA 2.0 survey, an ALMA blind survey at 1.1mm covering a continuous area of 72.42arcmin$^2$ using two array configurations. We derived physical properties, such as star formation rates, gas fractions, depletion timescales, and dust temperatures for the galaxy sample built from the survey. There exists a subset of galaxies that exhibit starburst-like short depletion timescales, but they are located within the scatter of the so-called main sequence of SFGs. These are dubbed starbursts in the main sequence and display the most compact star formation and they are characterized by the shortest depletion timescales, lowest gas fractions, and highest dust temperatures of the galaxy sample, compared to typical SFGs at the same stellar mass and redshift. They are also very massive, accounting for $\sim 60\%$ of the most massive galaxies in the sample ($\log (M_{\rm{*}}/M_{\odot}) > 11.0$). We find trends between the areas of the ongoing star formation regions and the derived physical properties for the sample, unveiling the role of compact star formation as a physical driver of these properties. Starbursts in the main sequence appear to be the extreme cases of these trends. We discuss possible scenarios of galaxy evolution to explain the results drawn from our galaxy sample. Our findings suggest that the star formation rate is sustained in SFGs by gas and star formation compression, keeping them within the main sequence even when their gas fractions are low and they are presumably on the way to quiescence.

preprint2022arXiv

NuSTAR Observations of Intrinsically X-ray Weak Quasar Candidates: An Obscuration-Only Scenario

We utilize recent NuSTAR observations (co-added depth $\approx55$-120 ks) of PG $1001+054$, PG $1254+047$, and PHL 1811 to constrain their hard X-ray ($\gtrsim5$ keV) weakness and spectral shapes, and thus to investigate the nature of their extreme X-ray weakness. These quasars showed very weak soft X-ray emission, and they were proposed to be intrinsically X-ray weak, with the X-ray coronae producing weak continuum emission relative to their optical/UV emission. However, the new observations suggest an alternative explanation. The NuSTAR 3-24 keV spectral shapes for PG $1001+054$ and PHL 1811 are likely flat (effective power-law photon indices $Γ_{\rm eff}=1.0^{+0.5}_{-0.6}$ and $Γ_{\rm eff}=1.4^{+0.8}_{-0.7}$, respectively), while the shape is nominal for PG $1254+047$ ($Γ_{\rm eff}=1.8\pm0.3$). PG $1001+054$ and PHL 1811 are significantly weak at hard X-ray energies (by factors of $\approx26$-74 at rest-frame 8 keV) compared to the expectations from their optical/UV emission, while PG $1254+047$ is only hard X-ray weak by a factor of $\approx3$. We suggest that X-ray obscuration is present in all three quasars. We propose that, as an alternative to the intrinsic X-ray weakness + X-ray obscuration scenario, the soft and hard X-ray weakness of these quasars can be uniformly explained under an obscuration-only scenario. This model provides adequate descriptions of the multi-epoch soft and hard X-ray data of these quasars, with variable column density and leaked fraction of the partial-covering absorber. We suggest that the absorber is the clumpy dust-free wind launched from the accretion disk. These quasars probably have super-Eddington accretion rates that drive powerful and high-density winds.

preprint2022arXiv

Quasar Feedback Survey: Multi-phase outflows, turbulence and evidence for feedback caused by low power radio jets inclined into the galaxy disk

We present a study of a luminous, z=0.15, type-2 quasar (log [L([OIII])/(erg/s)]=42.8) from the Quasar Feedback Survey. It is classified as 'radio-quiet' (log [L(1.4 GHz)/(W/Hz)]=23.8); however, radio imaging reveals ~1 kpc low-power jets (log [Pjet/(erg/s)]=44) inclined into the plane of the galaxy disk. We combine MUSE and ALMA observations to map stellar kinematics and ionised and molecular gas properties. The jets are seen to drive galaxy-wide bi-conical turbulent outflows, reaching W80 = 1000-1300 km/s, in the ionised phase (traced via optical emission-lines), which also have increased electron densities compared to the quiescent gas. The turbulent gas is driven perpendicular to the jet axis and is escaping along the galaxy minor axis, reaching 7.5 kpc on both sides. Traced via CO(3-2) emission, the turbulent material in molecular gas phase is one-third as spatially extended and has 3 times lower velocity-dispersion as compared to ionised gas. The jets are seen to be strongly interacting with the interstellar medium (ISM) through enhanced ionised emission and disturbed/depleted molecular gas at the jet termini. We see further evidence for jet-induced feedback through significantly higher stellar velocity-dispersion aligned, and co-spatial with, the jet axis (<5 deg). We discuss possible negative and positive feedback scenarios arising due to the interaction of the low-power jets with the ISM in the context of recent jet-ISM interaction simulations, which qualitatively agree with our observations. We discuss how jet-induced feedback could be an important feedback mechanism even in bolometrically luminous 'radio-quiet' quasars.

preprint2022arXiv

The Black Hole Mass Function Across Cosmic Times II. Heavy Seeds and (Super)Massive Black Holes

This is the second paper in a series aimed at modeling the black hole (BH) mass function, from the stellar to the (super)massive regime. In the present work we focus on (super)massive BHs and provide an ab-initio computation of their mass function across cosmic times. We consider two main mechanisms to grow the central BH, that are expected to cooperate in the high-redshift star-forming progenitors of local massive galaxies. The first is the gaseous dynamical friction process, that can cause the migration toward the nuclear regions of stellar-mass BHs originated during the intense bursts of star formation in the gas-rich host progenitor galaxy, and the buildup of a central heavy BH seed $M_\bullet\sim 10^{3-5}\, M_\odot$ within short timescales $\lesssim$ some $10^7$ yr. The second mechanism is the standard Eddington-type gas disk accretion onto the heavy BH seed, through which the central BH can become (super)massive $M_\bullet\sim 10^{6-10}\, M_\odot$ within the typical star-formation duration $\lesssim 1$ Gyr of the host. We validate our semi-empirical approach by reproducing the observed redshift-dependent bolometric AGN luminosity functions and Eddington ratio distributions, and the relationship between the star-formation and the bolometric luminosity of the accreting central BH. We then derive the relic (super)massive BH mass function at different redshifts via a generalized continuity equation approach, and compare it with present observational estimates. Finally, we reconstruct the overall BH mass function from the stellar to the (super)massive regime, over more than ten orders of magnitudes in BH mass.

preprint2022arXiv

Warm molecular and ionized gas kinematics in the Type-2 quasar J0945+1737

We analyze Near-Infrared Integral Field Spectrograph (NIFS) observations of the type-2 quasar (QSO2) SDSS J094521.33+173753.2 to investigate its warm molecular and ionized gas kinematics. This QSO2 has a bolometric luminosity of 10$^{45.7}$ erg s$^{-1}$ and a redshift of z = 0.128. The K-band spectra provided by NIFS cover a range of 1.99-2.40 $μ$m where low-ionization (Pa$α$ and Br$δ$), high ionization ([S XI]$λ$1.920 $μ$m and [Si~VI]$λ$1.963 $μ$m) and warm molecular lines (from H$_2$ 1-0S(5) to 1-0S(1)) are detected, allowing us to study the multi-phase gas kinematics. Our analysis reveals gas in ordinary rotation in all the emission lines detected and also outflowing gas in the case of the low- and high-ionization emission lines. In the case of the nuclear spectrum, which corresponds to a circular aperture of 0.3\arcsec~(686 pc) in diameter, the warm molecular lines can be characterized using a single Gaussian component of full width at half maximum (FWHM)= 350-400 km s$^{-1}$, while Pa$α$, Br$δ$, and [Si~VI] are best fitted with two blue-shifted Gaussian components of FWHM$\sim$800 and 1700 km s$^{-1}$, in addition to a narrow component of $\sim$300 km s$^{-1}$. We interpret the blue-shifted broad components as outflowing gas, which reaches the highest velocities, of up to $-$840 km s$^{-1}$, in the south-east direction (PA$\sim$125$^{\circ}$), extending up to a distance of $\sim$3.4 kpc from the nucleus. The ionized outflow has a maximum mass outflow rate of $\dot{\text{M}}_{\text{out, max}}$=42-51 M$_\odot$ yr$^{-1}$, and its kinetic power represents 0.1$\%$ of the quasar bolometric luminosity.

preprint2021arXiv

GOODS-ALMA 2.0: Source catalog, number counts, and prevailing compact sizes in 1.1 mm galaxies

Submillimeter/millimeter observations of dusty star-forming galaxies with the Atacama Large Millimeter/submillimeter Array (ALMA) have shown that dust continuum emission generally occurs in compact regions smaller than the stellar distribution. However, it remains to be understood how systematic these findings are. Studies often lack homogeneity in the sample selection, target discontinuous areas with inhomogeneous sensitivities, and suffer from modest $uv$ coverage coming from single array configurations. GOODS-ALMA is a 1.1mm galaxy survey over a continuous area of 72.42arcmin$^2$ at a homogeneous sensitivity. In this version 2.0, we present a new low resolution dataset and its combination with the previous high resolution dataset from the survey, improving the $uv$ coverage and sensitivity reaching an average of $σ= 68.4μ$Jy beam$^{-1}$. A total of 88 galaxies are detected in a blind search (compared to 35 in the high resolution dataset alone), 50% at $S/N_{peak} \geq 5$ and 50% at $3.5 \leq S/N_{peak} \leq 5$ aided by priors. Among them, 13 out of the 88 are optically dark or faint sources ($H$- or $K$-band dropouts). The sample dust continuum sizes at 1.1mm are generally compact, with a median effective radius of $R_{e} = 0"10 \pm 0"05$ (a physical size of $R_{e} = 0.73 \pm 0.29$kpc at the redshift of each source). Dust continuum sizes evolve with redshift and stellar mass resembling the trends of the stellar sizes measured at optical wavelengths, albeit a lower normalization compared to those of late-type galaxies. We conclude that for sources with flux densities $S_{1.1mm} > 1$mJy, compact dust continuum emission at 1.1mm prevails, and sizes as extended as typical star-forming stellar disks are rare. The $S_{1.1mm} < 1$mJy sources appear slightly more extended at 1.1mm, although they are still generally compact below the sizes of typical star-forming stellar disks.

preprint2021arXiv

The post-Herschel view of intrinsic AGN emission: constructing templates for galaxy and AGN emission at IR wavelengths

Measuring the star-forming properties of AGN hosts is key to our understanding of galaxy formation and evolution. However, this topic remains debated, partly due to the difficulties in separating the infrared (i.e. 1--1000 $μ$m) emission into AGN and star-forming components. Taking advantage of archival far-infrared data from Herschel, we present a new set of AGN and galaxy infrared templates, and introduce the spectral energy distribution fitting code IRAGNSEP. Both can be used to measure infrared host galaxy properties, free of AGN contamination. To build these, we used a sample of 100 local ($z$ < 0.3), low-to-high luminosity AGNs (i.e. $L_{\rm bol}~\sim~10^{42--46}~\rm erg~s^{-1}$), selected from the 105-month Swift - BAT X-ray survey, which have archival Spitzer - IRS spectra and Herschel photometry. We first built a set of seven galaxy templates using a sample of 55 star-forming galaxies selected via infrared diagnostics. Using these templates, combined with a flexible model for the AGN contribution, we extracted the intrinsic infrared emission of our AGN sample. We further demonstrate that we can reduce the diversity in the intrinsic shapes of AGN spectral energy distributions down to a set of three AGN templates, of which two represent AGN continuum, and one represents silicate emission. Our results indicate that, on average, the contribution of AGNs to the far-infrared ($λ~\gtrsim$ 50 $μ$m) is not as high as suggested by some recent work. We further show that the need for two infrared AGN continuum templates could be related to nuclear obscuration, where one of our templates appears dominated by the emission of the extended polar dust.

preprint2020arXiv

$NuSTAR$ Observations of Four Nearby X-ray Faint AGN: Low Luminosity or Heavy Obscuration?

We present $NuSTAR$ observations of four active galactic nuclei (AGN) located within 15 Mpc. These AGN, namely ESO 121-G6, NGC 660, NGC 3486 and NGC 5195, have observed X-ray luminosities of $L_{\rm 2-10\ keV, obs} \lesssim$ 10$^{39}$ erg s$^{-1}$, classifying them as low luminosity AGN (LLAGN). We perform broadband X-ray spectral analysis for the AGN by combining our $NuSTAR$ data with $Chandra$ or $XMM-Newton$ observations to directly measure their column densities ($N_{\rm H}$) and infer their intrinsic power. We complement our X-ray data with archival and new high angular resolution mid-infrared (mid-IR) data for all objects, except NGC 5195. Based on our X-ray spectral analysis, we found that both ESO 121-G6 and NGC 660 are heavily obscured ($N_{\rm H}$ > 10$^{23}$ cm$^{-2}$; $L_{\rm 2-10\ keV,\ int} \sim$ 10$^{41}$ erg s$^{-1}$), and NGC 660 may be Compton-thick. We also note that the X-ray flux and spectral slope for ESO 121-G6 have significantly changed over the last decade, indicating significant changes in the obscuration and potentially accretion rate. On the other hand, NGC 3486 and NGC 5195 appear to be unobscured and just mildly obscured, respectively, with $L_{\rm 2-10\ keV,\ int} <$ 10$^{39}$ erg s$^{-1}$; i.e., genuine LLAGN. Both of the heavily obscured AGN have $L_{\rm bol} >$ 10$^{41}$ erg s$^{-1}$ and $λ_{\rm Edd} \gtrsim$ 10$^{-3}$, and are detected in high angular resolution mid-IR imaging, indicating the presence of obscuring dust on nuclear scale. NGC 3486 however, is undetected in high-resolution mid-IR imaging, and the current data do not provide stringent constraints on the presence or absence of obscuring nuclear dust in the AGN.

preprint2020arXiv

A binning-free method reveals a continuous relationship between galaxies' AGN power and offset from main sequence

Studies investigating the relationship between AGN power and the star formation rates (SFRs) of their host galaxies often rely on averaging techniques -- such as stacking -- to incorporate information from non-detections. However, averages, and especially means, can be strongly affected by outliers and can therefore give a misleading indication of the "typical" case. Recently, a number of studies have taken a step further by binning their sample in terms of AGN power (approximated by the 2-10keV luminosity of the AGN), and investigating how the SFR distribution differs between these bins. These bin thresholds are often weakly motivated, and binning implicitly assumes that sources within the same bin have similar (or even identical) properties. In this paper, we investigate whether the distribution of host SFRs -- relative to the locus of the star-forming main sequence (i.e., $R_{\rm MS}$) -- changes continuously as a function of AGN power. We achieve this by using a hierarchical Bayesian model that completely removes the need to bin in AGN power. In doing so, we find strong evidence that the $R_{\rm MS}$ distribution changes with 2-10keV X-ray luminosity. The results suggest that higher X_ray luminosity AGNs have a tighter physical connection to the star-forming process than lower X-ray luminosity AGNs, at least within the $0.8 < z < 1.2$ redshift range considered here.

preprint2020arXiv

Fundamental differences in the radio properties of red and blue quasars: enhanced compact AGN emission in red quasars

We have recently used the Faint Images of the Radio Sky at Twenty-centimeters (FIRST) survey to show that red quasars have fundamentally different radio properties to typical blue quasars: a significant (factor $\sim3$) enhancement in the radio-detection fraction, which arises from systems around the radio-quiet threshold with compact ($<5''$) radio morphologies. To gain greater insight into these physical differences, here we use the DR14 Sloan Digital Sky Survey (SDSS) and more sensitive, higher resolution radio data from the Very Large Array (VLA) Stripe 82 (S82) and VLA-COSMOS 3 GHz (C3GHz) surveys. With the S82 data, we perform morphological analyses at a resolution and depth three times that of the FIRST radio survey, and confirm an enhancement in radio-faint and compact red quasars over typical quasars; we now also find tentative evidence for an enhancement in red quasars with slightly extended radio structures ($16-43$ kpc at $z=1.5$). These analyses are complemented by C3GHz, which is deep enough to detect radio emission from star-formation processes. From our data we find that the radio enhancement from red quasars is due to AGN activity on compact scales ($< 43$ kpc) for radio-intermediate-radio-quiet sources ($-5<R<-3.4$, where $R=L_{1.4GHz}/L_{6 μm}$), which decreases at $R<-5$ as the radio emission from star-formation starts to dilute the AGN component. Overall our results argue against a simple orientation scenario and are consistent with red quasars representing a younger, earlier phase in the overall evolution of quasars.

preprint2020arXiv

Fundamental differences in the radio properties of red and blue quasars: Insight from the LOFAR Two-metre Sky Survey (LoTSS)

Red quasi-stellar objects (QSOs) are a subset of the luminous end of the cosmic population of active galactic nuclei (AGN), most of which are reddened by intervening dust along the line-of-sight towards their central engines. In recent work from our team, we developed a systematic technique to select red QSOs from the Sloan Digital Sky Survey (SDSS), and demonstrated that they have distinctive radio properties using the Faint Images of the Radio Sky at Twenty centimeters (FIRST) radio survey. Here we expand our study using low-frequency radio data from the LOFAR Two-metre Sky Survey (LoTSS). With the improvement in depth that LoTSS offers, we confirm key results: compared to a control sample of normal "blue" QSOs matched in redshift and accretion power, red QSOs have a higher radio detection rate and a higher incidence of compact radio morphologies. For the first time, we also demonstrate that these differences arise primarily in sources of intermediate radio-loudness: radio-intermediate red QSOs are $\times 3$ more common than typical QSOs, but the excess diminishes among the most radio-loud and the most radio-quiet systems in our study. We develop Monte-Carlo simulations to explore whether differences in star formation could explain these results, and conclude that, while star formation is an important source of low-frequency emission among radio-quiet QSOs, a population of AGN-driven compact radio sources is the most likely cause for the distinct low-frequency radio properties of red QSOs. Our study substantiates the conclusion that fundamental differences must exist between the red and normal blue QSO populations.

preprint2020arXiv

KASHz: No evidence for ionised outflows instantaneously suppressing star formation in moderate luminosity AGN at $z$$\sim$$1.4$-$2.6$

As part of our KMOS AGN Survey at High-redshift (KASHz), we present spatially-resolved VLT/KMOS and VLT/SINFONI spectroscopic data and ALMA 870$μ$m continuum imaging of eight $z$=1.4--2.6 moderate AGN ($L_{\rm 2-10 \rm kev}$ = $10^{42} - 10^{45}$ ergs s$^{-1}$). We map [OIII], H$α$ and rest-frame FIR emission to search for any spatial anti-correlation between ionised outflows (traced by the [OIII] line) and star formation (SF; traced by H$α$ and FIR), that has previously been claimed for some high-z AGN and used as evidence for negative and/or positive AGN feedback. Firstly, we conclude that H$α$ is unreliable to map SF inside our AGN host galaxies based on: (i) SF rates inferred from attenuation-corrected H$α$ can lie below those inferred from FIR; (ii) the FIR continuum is more compact than the H$α$ emission by a factor of $\sim 2$ on average; (iii) in half of our sample, we observe significant spatial offsets between the FIR and H$α$ emission, with an average offset of $1.4\pm0.6$ kpc. Secondly, for the five targets with outflows we find no evidence for a spatial anti-correlation between outflows and SF using either H$α$ or FIR as a tracer. This holds for our re-analysis of a famous $z$=1.6 X-ray AGN (`XID 2028') where positive and negative feedback has been previously claimed. Based on our results, any impact on SF by ionised outflows must be subtle, either occurring on scales below our resolution, or on long timescales.

preprint2020arXiv

Local AGN Survey (LASr): I. Galaxy sample, infrared colour selection and predictions for AGN within 100 Mpc

In order to answer some of the major open questions in the fields of supermassive black hole (SMBH) and galaxy evolution, a complete census of SMBH growth, i.e., active galactic nuclei (AGN), is required. Thanks to deep all-sky surveys, such as those by the Wide-field Infrared Survey Explorer (WISE) and the Spectrum-Roentgen-Gamma (SRG) missions, this task is now becoming feasible in the nearby Universe. Here, we present a new survey, the Local AGN Survey (LASr), with the goal of identifying AGN unbiased against obscuration and determining the intrinsic Compton-thick (CT) fraction. First, we construct the most complete all-sky sample of galaxies within 100 Mpc from astronomical databases (90% completeness for log M_star/M_sun ~9.4), four times deeper than the current local galaxy reference, the Two Micron All-Sky Survey Redshift Survey (2MRS), which turns out to miss 20% of known luminous AGN. These 49k galaxies serve as parent sample for LASr, called LASr-GPS. It contains 4.3k already known AGN, >=82% of these are estimated to have L(12um) < 10^42.3 erg/s, i.e., are low-luminosity AGN. As a first method for identifying Seyfert-like AGN, we use WISE-based infrared colours, finding 221 galaxies at L(12um) >= 10^42.3 erg/s to host an AGN at 90% reliability, This includes 61 new AGN candidates and implies and optical type 2 fraction of 50 to 71%. We quantify the efficiency of this technique and estimate the total number of AGN with L_int(2-10 keV) >= 10^42 erg/s in the volume to be 362^+145_-116 (8.6^+3.5_-2.8 10^-5 Mpc^-3). X-ray brightness estimates indicate the CT fraction to be 40-55% to explain the Swift non-detections of the infrared selected objects. One third of the AGN within 100 Mpc remain to be identified and we discuss the prospects for the eROSITA all-sky survey.

preprint2020arXiv

Searching for obscured AGN in z $\sim$ 2 submillimetre galaxies

Submillimetre-selected galaxies (SMGs) at high redshift ($z$ $\sim$ 2) are potential host galaxies of active galactic nuclei (AGN). If the local Universe is a good guide, $\sim$ 50$\%$ of the obscured AGN amongst the SMG population could be missed even in the deepest X-ray surveys. Radio observations are insensitive to obscuration; therefore, very long baseline interferometry (VLBI) can be used as a tool to identify AGN in obscured systems. A well-established upper limit to the brightness temperature of 10$^5$ K exists in star-forming systems, thus VLBI observations can distinguish AGN from star-forming systems via brightness temperature measurements. We present 1.6 GHz European VLBI Network (EVN) observations of four SMGs (with measured redshifts) to search for evidence of compact radio components associated with AGN cores. For two of the sources, e-MERLIN images are also presented. Out of the four SMGs observed, we detect one source, J123555.14, that has an integrated EVN flux density of 201 $\pm$ 15.2 $μ$Jy, corresponding to a brightness temperature of 5.2 $\pm$ 0.7 $\times$ 10$^5$ K. We therefore identify that the radio emission from J123555.14 is associated with an AGN. We do not detect compact radio emission from a possible AGN in the remaining sources (J123600.10, J131225.73, and J163650.43). In the case of J131225.73, this is particularly surprising, and the data suggest that this may be an extended, jet-dominated AGN that is resolved by VLBI. Since the morphology of the faint radio source population is still largely unknown at these scales, it is possible that with a $\sim$ 10 mas resolution, VLBI misses (or resolves) many radio AGN extended on kiloparsec scales.

preprint2020arXiv

The Bolometric Quasar Luminosity Function at z = 0-7

In this paper, we provide updated constraints on the bolometric quasar luminosity function (QLF) from $z=0$ to $z=7$. The constraints are based on an observational compilation that includes observations in the rest-frame IR, B band, UV, soft and hard X-ray in past decades. Our method follows Hopkins et al. 2007 with an updated quasar SED model and bolometric and extinction corrections. The new best-fit bolometric quasar luminosity function behaves qualitatively different from the Hopkins et al. 2007 model at high redshift. Compared with the old model, the number density normalization decreases towards higher redshift and the bright-end slope is steeper at $z\gtrsim 2$. Due to the paucity of measurements at the faint end, the faint end slope at $z\gtrsim 5$ is quite uncertain. We present two models, one featuring a progressively steeper faint-end slope at higher redshift and the other featuring a shallow faint-end slope at $z\gtrsim 5$. Further multi-band observations of the faint-end QLF are needed to distinguish between these models. The evolutionary pattern of the bolometric QLF can be interpreted as an early phase likely dominated by the hierarchical assembly of structures and a late phase likely dominated by the quenching of galaxies. We explore the implications of this model on the ionizing photon production by quasars, the CXB spectrum, the SMBH mass density and mass functions. The predicted hydrogen photoionization rate contributed by quasars is subdominant during the epoch of reionization and only becomes important at $z\lesssim 3$. The predicted CXB spectrum, cosmic SMBH mass density and SMBH mass function are generally consistent with existing observations.

preprint2020arXiv

The Star-Formation Properties of the Observed and Simulated AGN Universe: BAT vs EAGLE

In this paper we present data from 72 low redshift, hard X-ray selected AGN taken from the {\it Swift}-BAT 58 month catalogue. We utilise spectral energy distribution fitting to the optical to IR photometry in order to estimate host galaxy properties. We compare this observational sample to a volume and flux matched sample of AGN from the EAGLE hydrodynamical simulations in order to verify how accurately the simulations can reproduce observed AGN host galaxy properties. After correcting for the known +0.2 dex offset in the SFRs between EAGLE and previous observations, we find agreement in the SFR and X-ray luminosity distributions; however we find that the stellar masses in EAGLE are $0.2 - 0.4$ dex greater than the observational sample, which consequently leads to lower sSFRs. We compare these results to our previous study at high redshift, finding agreement in both the observations and simulations, whereby the widths of sSFR distributions are similar ($\sim0.4-0.6$ dex) and the median of the SFR distributions lie below the star forming main sequence by $\sim0.3-0.5$ dex across all samples. We also use EAGLE to select a sample of AGN host galaxies at high and low redshift and follow their characteristic evolution from $z=8$ to $z=0$. We find similar behaviour between these two samples, whereby star formation is quenched when the black hole goes through its phase of most rapid growth. Utilising EAGLE we find that 23\% of AGN selected at $z\sim0$ are also AGN at high redshift, and that their host galaxies are among the most massive objects in the simulation. Overall we find EAGLE reproduces the observations well, with some minor inconsistencies ($\sim$ 0.2 dex in stellar masses and $\sim$ 0.4 dex in sSFRs).

preprint2019arXiv

Do AGN really suppress star formation?

Active galactic nuclei (AGN) are believed to regulate star formation inside their host galaxies through "AGN feedback". We summarise our on-going study of luminous AGN (z~0.2-3; L_(AGN,bol)>~10^43 erg/s), which is designed to search for observational signatures of feedback by combining observed star-formation rate (SFR) measurements from statistical samples with cosmological model predictions. Using the EAGLE hydrodynamical cosmological simulations, in combination with our Herschel+ALMA surveys, we show that - even in the presence of AGN feedback - we do not necessarily expect to see any relationships between average galaxy-wide SFRs and instantaneous AGN luminosities. We caution that the correlation with stellar mass for both SFR and AGN luminosity can contribute to apparent observed positive trends between these two quantities. On the other hand, the EAGLE simulations, which reproduce our observations, predict that a signature of AGN feedback can be seen in the wide specific SFR distributions of $all$ massive galaxies (not just AGN hosts). Overall, whilst we can not rule out that AGN have an immediate small-scale impact on in-situ star-formation, all of our results are consistent with a feedback model where galaxy-wide in-situ star formation is not rapidly suppressed by AGN, but where the feedback likely acts over a longer timescale than a single AGN episode.

preprint2016arXiv

A New Population of Compton-Thick AGN Identified Using the Spectral Curvature Above 10 keV

We present a new metric that uses the spectral curvature (SC) above 10 keV to identify Compton-thick AGN in low-quality Swift BAT X-ray data. Using NuSTAR, we observe nine high SC-selected AGN. We find that high-sensitivity spectra show the majority are Compton-thick (78% or 7/9) and the remaining two are nearly Compton-thick (NH~5-8x10^23 cm^-2). We find the SC_bat and SC_nustar measurements are consistent, suggesting this technique can be applied to future telescopes. We tested the SC method on well-known Compton-thick AGN and find it is much more effective than broad band ratios (e.g. 100% using SC vs. 20% using 8-24/3-8 keV). Our results suggest that using the >10 keV emission may be the only way to identify this population since only two sources show Compton-thick levels of excess in the OIII to X-ray emission ratio (F_OIII/F_2-10 keV>1) and WISE colors do not identify most of them as AGN. Based on this small sample, we find that a higher fraction of these AGN are in the final merger stage than typical BAT AGN. Additionally, these nine obscured AGN have, on average, four times higher accretion rates than other BAT-detected AGN (Edd ratio=0.068 compared to 0.016). The robustness of SC at identifying Compton-thick AGN implies a higher fraction of nearby AGN may be Compton-thick (~22%) and the sum of black hole growth in Compton-thick AGN (Eddington ratio times population percentage), is nearly as large as mildly obscured and unobscured AGN.

preprint2016arXiv

ALMA observations of a z~3.1 Protocluster: Star Formation from Active Galactic Nuclei and Lyman-Alpha Blobs in an Overdense Environment

We exploit ALMA 870um observations to measure the star-formation rates (SFRs) of eight X-ray detected Active Galactic Nuclei (AGNs) in a z~3.1 protocluster, four of which reside in extended Ly-alpha haloes (often termed Ly-alpha blobs: LABs). Three of the AGNs are detected by ALMA and have implied SFRs of ~220-410~M_sun/yr; the non detection of the other five AGNs places SFR upper limits of <210 M_sun/yr. The mean SFR of the protocluster AGNs (~110-210 M_sun/yr) is consistent (within a factor of ~0.7-2.3) with that found for co-eval AGNs in the field, implying that galaxy growth is not significantly accelerated in these systems. However, when also considering ALMA data from the literature, we find evidence for elevated mean SFRs (up-to a factor of ~5.9 over the field) for AGNs at the protocluster core, indicating that galaxy growth is significantly accelerated in the central regions of the protocluster. We also show that all of the four protocluster LABs are associated with an ALMA counterpart within the extent of their Ly-alpha emission. The SFRs of the ALMA sources within the LABs (~150-410 M_sun/yr) are consistent with those expected for co-eval massive star-forming galaxies in the field. Furthermore, the two giant LABs (with physical extents of >100 kpc) do not host more luminous star formation than the smaller LABs, despite being an order of magnitude brighter in Ly-alpha emission. We use these results to discuss star formation as the power source of LABs.

preprint2016arXiv

ALMA resolves extended star formation in high-z AGN host galaxies

We present high resolution (0.3") Atacama Large Millimeter Array (ALMA) 870um imaging of five z~1.5-4.5 X-ray detected AGN with luminosities of L(2-8keV)>10^42 erg/s. These data provide a >~20x improvement in spatial resolution over single-dish rest-frame FIR measurements. The sub-millimetre emission is extended on scales of FWHM~0.2"-0.5", corresponding to physical sizes of 1-3 kpc (median value of 1.8 kpc). These sizes are comparable to the majority of z=1-5 sub-millimetre galaxies (SMGs) with equivalent ALMA measurements. In combination with spectral energy distribution analyses, we attribute this rest-frame far-infrared (FIR) emission to dust heated by star formation. The implied star-formation rate surface densities are ~20-200 Mo/yr/kpc^2, which are consistent with SMGs of comparable FIR luminosities (i.e., L(IR)~[1-5]x10^(12)Lo). Although limited by a small sample of AGN, which all have high FIR luminosities, our study suggests that the kpc-scale spatial distribution and surface density of star formation in high-redshift star-forming galaxies is the same irrespective of the presence of X-ray detected AGN.

preprint2016arXiv

First NuSTAR observations of the BL Lac - type blazar PKS~2155-304: constraints on the jet content and distribution of radiating particles

We report the first hard X-ray observations with NuSTAR of the BL Lac type blazar PKS 2155-304, augmented with soft X-ray data from XMM-Newton and gamma-ray data from the Fermi Large Area Telescope, obtained in April 2013 when the source was in a very low flux state. A joint NuSTAR and XMM spectrum, covering the energy range 0.5 - 60 keV, is best described by a model consisting of a log-parabola component with curvature beta = 0.3(+0.2,-0.1) and a (local) photon index 3.04 +/- 0.15 at photon energy of 2 keV, and a hard power-law tail with photon index 2.2 +/- 0.4. The hard X-ray tail can be smoothly joined to the quasi-simultaneous gamma-ray spectrum by a synchrotron self-Compton component produced by an electron distribution with index p = 2.2. Assuming that the power-law electron distribution extends down to the minimum electron Lorentz factor gamma_min = 1 and that there is one proton per electron, an unrealistically high total jet power L_p of roughly 10^47 erg/s is inferred. This can be reduced by two orders of magnitude either by considering a significant presence of electron-positron pairs with lepton-to-proton ratio of at least 30, or by introducing an additional, low-energy break in the electron energy distribution at the electron Lorentz factor gamma_br1 of roughly 100. In either case, the jet composition is expected to be strongly matter-dominated.

preprint2016arXiv

Herschel protocluster survey: A search for dusty star-forming galaxies in protoclusters at z=2-3

We present a Herschel/SPIRE survey of three protoclusters at z=2-3 (2QZCluster, HS1700, SSA22). Based on the SPIRE colours (S350/S250 and S500/S350) of 250 $μ$m sources, we selected high redshift dusty star-forming galaxies potentially associated with the protoclusters. In the 2QZCluster field, we found a 4-sigma overdensity of six SPIRE sources around 4.5' (~2.2 Mpc) from a density peak of H$α$ emitters at z=2.2. In the HS1700 field, we found a 5-sigma overdensity of eight SPIRE sources around 2.1' (~1.0 Mpc) from a density peak of LBGs at z=2.3. We did not find any significant overdensities in SSA22 field, but we found three 500 $μ$m sources are concentrated 3' (~1.4 Mpc) east to the LAEs overdensity. If all the SPIRE sources in these three overdensities are associated with protoclusters, the inferred star-formation rate densities are 10$^3$-10$^4$ times higher than the average value at the same redshifts. This suggests that dusty star-formation activity could be very strongly enhanced in z~2-3 protoclusters. Further observations are needed to confirm the redshifts of the SPIRE sources and to investigate what processes enhance the dusty star-formation activity in z~2-3 protoclusters.

preprint2016arXiv

Kiloparsec-scale dust disks in high-redshift luminous submillimeter galaxies

We present high-resolution (0.16$"$) 870um Atacama Large Millimeter/submillimeter Array (ALMA) imaging of 16 luminous (L_IR ~ 4 x 10^12 L_sun) submillimeter galaxies (SMGs) from the ALESS survey of the Extended Chandra Deep Field South. This dust imaging traces the dust-obscured star formation in these z~2.5 galaxies on ~1.3 kpc scales. The emission has a median effective radius of $R_e=0.24" \pm 0.02"$, corresponding to a typical physical size of $R_{e}=1.8\pm$0.2 kpc. We derive a median Sérsic index of $n=0.9\pm0.2$, implying that the dust emission is remarkably disk-like at the current resolution and sensitivity. We use different weighting schemes with the visibilities to search for clumps on 0.12$"$ (~1.0 kpc) scales, but we find no significant evidence for clumping in the majority of cases. Indeed, we demonstrate using simulations that the observed morphologies are generally consistent with smooth exponential disks, suggesting that caution should be exercised when identifying candidate clumps in even moderate S/N interferometric data. We compare our maps to comparable-resolution HST H$_{160}$-band images, finding that the stellar morphologies appear significantly more extended and disturbed, and suggesting that major mergers may be responsible for driving the formation of the compact dust disks we observe. The stark contrast between the obscured and unobscured morphologies may also have implications for SED fitting routines that assume the dust is co-located with the optical/near-IR continuum emission. Finally, we discuss the potential of the current bursts of star formation to transform the observed galaxy sizes and light profiles, showing that the $z\sim0$ descendants of these SMGs are expected to have stellar masses, effective radii, and gas surface densities consistent with the most compact massive (M_* ~ 1-2 x 10^11 M_sun) early-type galaxies observed locally.

preprint2016arXiv

The 2 Ms Chandra Deep Field-North Survey and the 250 ks Extended Chandra Deep Field-South Survey: Improved Point-Source Catalogs

We present improved point-source catalogs for the 2 Ms Chandra Deep Field-North (CDF-N) and the 250 ks Extended Chandra Deep Field-South (E-CDF-S), implementing a number of recent improvements in Chandra source-cataloging methodology. For the CDF-N/E-CDF-S, we provide a main catalog that contains 683/1003 X-ray sources detected with wavdetect at a false-positive probability threshold of $10^{-5}$ that also satisfy a binomial-probability source-selection criterion of $P<0.004$/$P<0.002$. Such an approach maximizes the number of reliable sources detected: a total of 196/275 main-catalog sources are new compared to the Alexander et al. (2003) CDF-N/Lehmer et al. (2005) E-CDF-S main catalogs. We also provide CDF-N/E-CDF-S supplementary catalogs that consist of 72/56 sources detected at the same wavdetect threshold and having $P$ of $0.004-0.1$/$0.002-0.1$ and $K_s\le22.9/K_s\le22.3$ mag counterparts. For all $\approx1800$ CDF-N and E-CDF-S sources, including the $\approx500$ newly detected ones (these being generally fainter and more obscured), we determine X-ray source positions utilizing centroid and matched-filter techniques; we also provide multiwavelength identifications, apparent magnitudes of counterparts, spectroscopic and/or photometric redshifts, basic source classifications, and estimates of observed AGN and galaxy source densities around respective field centers. Simulations show that both the CDF-N and E-CDF-S main catalogs are highly reliable and reasonably complete. Background and sensitivity analyses indicate that the on-axis mean flux limits reached represent a factor of $\approx1.5-2.0$ improvement over the previous CDF-N and E-CDF-S limits. We make our data products publicly available.

preprint2016arXiv

The average submillimetre properties of Lyman-alpha Blobs at z=3

Ly-alpha blobs (LABs) offer insight into the complex interface between galaxies and their circumgalactic medium. Whilst some LABs have been found to contain luminous star-forming galaxies and active galactic nuclei that could potentially power the Ly-alpha emission, others appear not to be associated with obvious luminous galaxy counterparts. It has been speculated that LABs may be powered by cold gas streaming on to a central galaxy, providing an opportunity to directly observe the `cold accretion' mode of galaxy growth. Star-forming galaxies in LABs could be dust obscured and therefore detectable only at longer wavelengths. We stack deep SCUBA-2 observations of the SSA22 field to determine the average 850um flux density of 34 LABs. We measure S_850 = 0.6 +/- 0.2mJy for all LABs, but stacking the LABs by size indicates that only the largest third (area > 1794 kpc^2) have a mean detection, at 4.5 sigma, with S_850 = 1.4 +/- 0.3mJy. Only two LABs (1 and 18) have individual SCUBA-2 > 3.5 sigma detections at a depth of 1.1mJy/beam. We consider two possible mechanisms for powering the LABs and find that central star formation is likely to dominate the emission of Ly-alpha, with cold accretion playing a secondary role.

preprint2016arXiv

The Chandra COSMOS Legacy survey: overview and point source catalog

The COSMOS-Legacy survey is a 4.6 Ms Chandra program that has imaged 2.2 deg$^2$ of the COSMOS field with an effective exposure of $\simeq$160 ks over the central 1.5 deg$^2$ and of $\simeq$80 ks in the remaining area. The survey is the combination of 56 new observations, obtained as an X-ray Visionary Project, with the previous C-COSMOS survey. We describe the reduction and analysis of the new observations and the properties of 2273 point sources detected above a spurious probability of 2$\times 10^{-5}$. We also present the updated properties of the C-COSMOS sources detected in the new data. The whole survey includes 4016 point sources (3814, 2920 and 2440 in the full, soft and hard band). The limiting depths are 2.2 $\times$ 10$^{-16}$, 1.5 $\times$ 10$^{-15}$ and 8.9$\times$ 10$^{-16}$ ${\rm erg~cm}^{-2}~{\rm s}^{-1}$ in the 0.5-2, 2-10 and 0.5-10 keV bands, respectively. The observed fraction of obscured AGN with column density $> 10^{22}$ cm$^{-2}$ from the hardness ratio (HR) is $\sim$50$^{+17}_{-16}$%. Given the large sample, we compute source number counts in the hard and soft bands, significantly reducing the uncertainties of 5-10%. For the first time, we compute number counts for obscured (HR$>$-0.2) and unobscured (HR$<$-0.2) sources and find significant differences between the two populations in the soft band. Due to the un-precedent large exposure, COSMOS-Legacy area is 3 times larger than surveys at similar depth and its depth is 3 times fainter than surveys covering similar area. The area-flux region occupied by COSMOS-Legacy is likely to remain unsurpassed for years to come.

preprint2016arXiv

The Chandra Deep Field-South Survey: 7 Ms Source Catalogs

We present X-ray source catalogs for the $\approx7$ Ms exposure of the Chandra Deep Field-South (CDF-S), which covers a total area of 484.2 arcmin$^2$. Utilizing WAVDETECT for initial source detection and ACIS Extract for photometric extraction and significance assessment, we create a main source catalog containing 1008 sources that are detected in up to three X-ray bands: 0.5-7.0 keV, 0.5-2.0 keV, and 2-7 keV. A supplementary source catalog is also provided including 47 lower-significance sources that have bright ($K_s\le23$) near-infrared counterparts. We identify multiwavelength counterparts for 992 (98.4%) of the main-catalog sources, and we collect redshifts for 986 of these sources, including 653 spectroscopic redshifts and 333 photometric redshifts. Based on the X-ray and multiwavelength properties, we identify 711 active galactic nuclei (AGNs) from the main-catalog sources. Compared to the previous $\approx4$ Ms CDF-S catalogs, 291 of the main-catalog sources are new detections. We have achieved unprecedented X-ray sensitivity with average flux limits over the central $\approx1$ arcmin$^2$ region of $\approx1.9\times10^{-17}$, $6.4\times10^{-18}$, and $2.7\times10^{-17}$ erg cm$^{-2}$ s$^{-1}$ in the three X-ray bands, respectively. We provide cumulative number-count measurements observing, for the first time, that normal galaxies start to dominate the X-ray source population at the faintest 0.5-2.0 keV flux levels. The highest X-ray source density reaches $\approx50\,500$ deg$^{-2}$, and $47\%\pm4\%$ of these sources are AGNs ($\approx23\,900$ deg$^{-2}$).

preprint2016arXiv

The NuSTAR Extragalactic Surveys: The Number Counts of Active Galactic Nuclei and the Resolved Fraction of the Cosmic X-ray Background

We present the 3-8 keV and 8-24 keV number counts of active galactic nuclei (AGN) identified in the NuSTAR extragalactic surveys. NuSTAR has now resolved 33-39% of the X-ray background in the 8-24 keV band, directly identifying AGN with obscuring columns up to approximately 1e25 /cm2. In the softer 3-8 keV band the number counts are in general agreement with those measured by XMM-Newton and Chandra over the flux range 5e-15 < S(3 - 8 keV)/(erg/cm2/s) < 1e-12 probed by NuSTAR. In the hard 8-24 keV band NuSTAR probes fluxes over the range 2e-14 < S(8-24 keV)/(erg/cm2/s) < 1e-12, a factor of approximately 100 fainter than previous measurements. The 8-24 keV number counts match predictions from AGN population synthesis models, directly confirming the existence of a population of obscured and/or hard X-ray sources inferred from the shape of the integrated cosmic X-ray background. The measured NuSTAR counts lie significantly above simple extrapolation with a Euclidian slope to low flux of the Swift/BAT 15-55 keV number counts measured at higher fluxes S(15-55 keV) > 1e-11 erg/cm2/s, reflecting the evolution of the AGN population between the Swift/BAT local (z<0.1) sample and NuSTAR's z~1 sample. CXB synthesis models, which account for AGN evolution, lie above the Swift/BAT measurements, suggesting that they do not fully capture the evolution of obscured AGN at low redshifts.

preprint2016arXiv

The SCUBA-2 Cosmology Legacy Survey: 850um maps, catalogues and number counts

We present a catalogue of nearly 3,000 submillimetre sources detected at 850um over ~5 square degrees surveyed as part of the James Clerk Maxwell Telescope (JCMT) SCUBA-2 Cosmology Legacy Survey (S2CLS). This is the largest survey of its kind at 850um, probing a meaningful cosmic volume at the peak of star formation activity and increasing the sample size of submillimetre galaxies selected at 850um by an order of magnitude. We describe the wide 850um survey component of S2CLS, which covers the key extragalactic survey fields: UKIDSS-UDS, COSMOS, Akari-NEP, Extended Groth Strip, Lockman Hole North, SSA22 and GOODS-North. The average 1-sigma depth of S2CLS is 1.2 mJy/beam, approaching the SCUBA-2 850um confusion limit, which we determine to be ~0.8 mJy/beam. We measure the single dish 850um number counts to unprecedented accuracy, reducing the Poisson errors on the differential counts to approximately 4% at S_850~3mJy. With several independent fields, we investigate field-to-field variance, finding that the number counts on 0.5-1 degree scales are generally within 50% of the S2CLS mean for S_850>3mJy, with scatter consistent with the Poisson and estimated cosmic variance uncertainties, although there is a marginal (2-sigma) density enhancement in the GOODS-North field. The observed number counts are in reasonable agreement with recent phenomenological and semi-analytic models. Finally, the large solid angle of S2CLS allows us to measure the bright-end counts: at S_850>10mJy there are approximately ten sources per square degree, and we detect the distinctive up-turn in the number counts indicative of the detection of local sources of 850um emission and strongly lensed high-redshift galaxies. Here we describe the data collection and reduction procedures and present calibrated maps and a catalogue of sources; these are made publicly available.

preprint2016arXiv

VLA and ALMA Imaging of Intense, Galaxy-Wide Star Formation in z ~ 2 Galaxies

We present $\simeq$0$.\!\!^{\prime\prime}4$-resolution extinction-independent distributions of star formation and dust in 11 star-forming galaxies (SFGs) at $z = 1.3-3.0$. These galaxies are selected from sensitive, blank-field surveys of the $2' \times 2'$ Hubble Ultra-Deep Field at $λ= 5$ cm and 1.3 mm using the Karl G. Jansky Very Large Array (VLA) and Atacama Large Millimeter/submillimeter Array (ALMA). They have star-formation rates (SFRs), stellar masses, and dust properties representative of massive main-sequence SFGs at $z \sim 2$. Morphological classification performed on spatially-resolved stellar mass maps indicates a mixture of disk and morphologically disturbed systems; half of the sample harbor X-ray active galactic nuclei (AGN), thereby representing a diversity of $z \sim 2$ SFGs undergoing vigorous mass assembly. We find that their intense star formation most frequently occurs at the location of stellar-mass concentration and extends over an area comparable to their stellar-mass distribution, with a median diameter of $4.2 \pm 1.8$ kpc. This provides direct evidence for galaxy-wide star formation in distant, blank-field-selected main-sequence SFGs. The typical galactic-average SFR surface density is 2.5 M$_{\odot}$yr$^{-1}$kpc$^{-2}$, sufficiently high to drive outflows. In X-ray-selected AGN where radio emission is enhanced over the level associated with star formation, the radio excess pinpoints the AGN, which are found to be co-spatial with star formation. The median extinction-independent size of main-sequence SFGs is two times larger than those of bright submillimeter galaxies whose SFRs are $3-8$ times larger, providing a constraint on the characteristic SFR ($\sim300$ M$_{\odot}$yr$^{-1}$) above which a significant population of more compact star-forming galaxies appears to emerge.

preprint2015arXiv

A NuSTAR Survey of Nearby Ultraluminous Infrared Galaxies

We present a NuSTAR, Chandra, and XMM--Newton survey of nine of the nearest ultraluminous infrared galaxies (ULIRGs). The unprecedented sensitivity of NuSTAR at energies above 10 keV enables spectral modeling with far better precision than was previously possible. Six of the nine sources observed were detected sufficiently well by NuSTAR to model in detail their broadband X-ray spectra, and recover the levels of obscuration and intrinsic X-ray luminosities. Only one source (IRAS 13120--5453) has a spectrum consistent with a Compton--thick AGN, but we cannot rule out that a second source (Arp 220) harbors an extremely highly obscured AGN as well. Variability in column density (reduction by a factor of a few compared to older observations) is seen in IRAS 05189--2524 and Mrk 273, altering the classification of these border-line sources from Compton-thick to Compton-thin. The ULIRGs in our sample have surprisingly low observed fluxes in high energy (>10 keV) X-rays, especially compared to their bolometric luminosities. They have lower ratios of unabsorbed 2--10 keV to bolometric luminosity, and unabsorbed 2--10 keV to mid-IR [O IV] line luminosity than do Seyfert 1 galaxies. We identify IRAS 08572+3915 as another candidate intrinsically X-ray weak source, similar to Mrk 231. We speculate that the X-ray weakness of IRAS 08572+3915 is related to its powerful outflow observed at other wavelengths.

preprint2015arXiv

A remarkably flat relationship between the average star formation rate and AGN luminosity for distant X-ray AGN

In this study we investigate the relationship between the star formation rate, SFR, and AGN luminosity, L(AGN), for ~2000 X-ray detected AGN. The AGN span over three orders of magnitude in X-ray luminosity (10^(42) < L(2-8keV) < 10^(45.5) erg/s) and are in the redshift range z = 0.2 - 2.5. Using infrared (IR) photometry (8 - 500um), including deblended Spitzer and Herschel images and taking into account photometric upper limits, we decompose the IR spectral energy distributions into AGN and star formation components. Using the IR luminosities due to star formation, we investigate the average SFRs as a function of redshift and AGN luminosity. In agreement with previous studies, we find a strong evolution of the average SFR with redshift, tracking the observed evolution of the overall star forming galaxy population. However, we find that the relationship between the average SFR and AGN luminosity is flat at all redshifts and across all the AGN luminosities investigated; in comparison to previous studies, we find less scatter amongst the average SFRs across the wide range of AGN luminosities investigated. By comparing to empirical models, we argue that the observed flat relationship is due to short timescale variations in AGN luminosity, driven by changes in the mass accretion rate, which wash out any underlying correlations between SFR and L(AGN). Furthermore, we show that the exact form of the predicted relationship between SFR and AGN luminosity (and it's normalisation) is highly sensitive to the assumed intrinsic Eddington ratio distribution.

preprint2015arXiv

ALMA and Herschel reveal that X-ray selected AGN and main-sequence galaxies have different star formation rate distributions

Using deep Herschel and ALMA observations, we investigate the star formation rate (SFR) distributions of X-ray selected AGN host galaxies at 0.5<z<1.5 and 1.5<z<4, comparing them to that of normal, star-forming (i.e., "main-sequence", or MS) galaxies. We find 34--55 per cent of AGNs in our sample have SFRs at least a factor of two below that of the average MS galaxy, compared to ~15 per cent of all MS galaxies, suggesting significantly different SFR distributions. Indeed, when both are modelled as log-normal distributions, the mass and redshift-normalised SFR distributions of X-ray AGNs are roughly twice as broad, and peak ~0.4 dex lower, than that of MS galaxies. However, like MS galaxies, the normalised SFR distribution of AGNs in our sample appears not to evolve with redshift. Despite X-ray AGNs and MS galaxies having different SFR distributions, the linear-mean SFR of AGNs derived from our distributions is remarkably consistent with that of MS galaxies, and thus with previous results derived from stacked Herschel data. This apparent contradiction is due to the linear-mean SFR being biased by bright outliers, and thus does not necessarily represent a true characterisation of the typical SFR of X-ray AGNs.

preprint2015arXiv

An enhanced merger fraction within the galaxy population of the SSA22 protocluster at z ~ 3.1

The over-dense environments of protoclusters of galaxies in the early Universe (z>2) are expected to accelerate the evolution of galaxies, with an increased rate of stellar mass assembly and black hole accretion compared to co-eval galaxies in the average density `field'. These galaxies are destined to form the passive population of massive systems that dominate the cores of rich clusters today. While signatures of accelerated growth of galaxies in the SSA22 protocluster (z=3.1) have been observed, the mechanism driving this remains unclear. Here we show an enhanced rate of galaxy-galaxy mergers could be responsible. We morphologically classify Lyman-break Galaxies (LBGs) in the SSA22 protocluster and compare these to those of galaxies in the field at z=3.1 as either active mergers or non-merging using Hubble Space Telescope ACS/F814W imaging, probing the rest-frame ultraviolet stellar light. We measure a merger fraction of 48+/-10 per cent for LBGs in the protocluster compared to 30+/-6 per cent for the field. Although the excess is marginal, an enhanced rate of mergers in SSA22 hints that galaxy-galaxy mergers are one of the key channels driving accelerated star formation and AGN growth in protocluster environments.

preprint2015arXiv

An Extragalactic Spectroscopic Survey of the SSA22 Field

We present VLT VIMOS, Keck DEIMOS and Keck LRIS multi-object spectra of 367 sources in the field of the z ~ 3.09 protocluster SSA22. Sources are spectroscopically classified via template matching, allowing new identifications for 206 extragalactic sources, including 36 z > 2 Lyman-break galaxies (LBGs) and Lyman-α emitters (LAEs), 8 protocluster members, and 94 X-ray sources from the ~ 400 ks Chandra deep survey of SSA22. Additionally, in the area covered by our study, we have increased by ~ 4, 13, and 6 times the number of reliable redshifts of sources at 1.0 < z < 2.0, at z > 3.4, and with X-Ray emission, respectively. We compare our results with past spectroscopic surveys of SSA22 to investigate the completeness of the LBGs and the X-Ray properties of the new spectroscopically-classified sources in the SSA22 field.

preprint2015arXiv

Cold Dust Emission from X-ray AGN in the SCUBA-2 Cosmology Legacy Survey: Dependence on Luminosity, Obscuration & AGN Activity

We study the 850um emission in X-ray selected AGN in the 2 sq-deg COSMOS field using new data from the SCUBA-2 Cosmology Legacy Survey. We find 19 850um bright X-ray AGN in a high-sensitivity region covering 0.89 sq-deg with flux densities of S850=4-10 mJy. The 19 AGN span the full range in redshift and hard X-ray luminosity covered by the sample - 0.7<z<3.5 and 43.2<log10(LX) <45. We report a highly significant stacked 850um detection of a hard X-ray flux-limited population of 699 z>1 X-ray AGN - S850=0.71+/-0.08mJy. We explore trends in the stacked 850um flux densities with redshift, finding no evolution in the average cold dust emission over the redshift range probed. For Type 1 AGN, there is no significant correlation between the stacked 850um flux and hard X-ray luminosity. However, in Type 2 AGN the stacked submm flux is a factor of 2 higher at high luminosities. When averaging over all X-ray luminosities, no significant differences are found in the stacked submm fluxes of Type 1 and Type 2 AGN as well as AGN separated on the basis of X-ray hardness ratios and optical-to-infrared colours. However, at log10(LX) >44.4, dependences in average submm flux on the optical-to-infrared colours become more pronounced. We argue that these high luminosity AGN represent a transition from a secular to a merger-driven evolutionary phase where the star formation rates and accretion luminosities are more tightly coupled. Stacked AGN 850um fluxes are compared to the stacked fluxes of a mass-matched sample of K-band selected non-AGN galaxies. We find that at 10.5<log10(M*/M0)<11.5, the non-AGN 850um fluxes are 1.5-2x higher than in Type 2 AGN of equivalent mass. We suggest these differences are due to the presence of massive dusty, red starburst galaxies in the K-band selected non-AGN sample, which are not present in optically selected catalogues covering a smaller area.

preprint2015arXiv

Cosmic X-ray Surveys of Distant Active Galaxies: The Demographics, Physics, and Ecology of Growing Supermassive Black Holes

We review results from cosmic X-ray surveys of active galactic nuclei (AGNs) over the past ~ 15 yr that have dramatically improved our understanding of growing supermassive black holes in the distant universe. First, we discuss the utility of such surveys for AGN investigations and the capabilities of the missions making these surveys, emphasizing Chandra, XMM-Newton, and NuSTAR. Second, we briefly describe the main cosmic X-ray surveys, the essential roles of complementary multiwavelength data, and how AGNs are selected from these surveys. We then review key results from these surveys on the AGN population and its evolution ("demographics"), the physical processes operating in AGNs ("physics"), and the interactions between AGNs and their environments ("ecology"). We conclude by describing some significant unresolved questions and prospects for advancing the field.

preprint2015arXiv

Mid-infrared luminous quasars in the GOODS-Herschel fields: a large population of heavily-obscured, Compton-thick quasars at z~2

We present the infrared (IR) and X-ray properties of a sample of 33 mid-IR luminous quasars ($ν$L(6 micron)>6x10$^{44}$ erg/s) at redshift z~1-3, identified through detailed spectral energy distribution analyses of distant star-forming galaxies, using the deepest IR data from Spitzer and Herschel in the GOODS-Herschel fields. The aim is to constrain the fraction of obscured, and Compton-thick (CT, N$_H$>1.5x10$^{24}$ cm$^{-2}$) quasars at the peak era of nuclear and star-formation activities. Despite being very bright in the mid-IR band, ~30% of these quasars are not detected in the extremely deep 2 Ms and 4 Ms Chandra X-ray data available in these fields. X-ray spectral analysis of the detected sources reveals that the majority (~67%) are obscured by column densities N$_H$>10$^{22}$ cm$^{-2}$; this fraction reaches ~80% when including the X-ray undetected sources (9 out of 33), which are likely to be the most heavily-obscured, CT quasars. We constrain the fraction of CT quasars in our sample to be ~24-48%, and their space density to be $Φ$=(6.7$\pm$2.2)x10$^{-6}$ Mpc$^{-3}$. From the investigation of the quasar host galaxies in terms of star-formation rates (SFRs) and morphological distortions, as a sign of galaxy mergers/interactions, we do not find any direct relation between SFRs and quasar luminosity or X-ray obscuration. On the other hand, there is tentative evidence that the most heavily-obscured quasars have, on average, more disturbed morphologies than the unobscured/moderately-obscured quasar hosts, which preferentially live in undisturbed systems. However, the fraction of quasars with disturbed morphology amongst the whole sample is ~40%, suggesting that galaxy mergers are not the main fuelling mechanism of quasars at z~2.

preprint2015arXiv

NuSTAR Observations of the Compton-thick Active Galactic Nucleus and Ultraluminous X-ray Source Candidate in NGC 5643

We present two NuSTAR observations of the local Seyfert 2 active galactic nucleus (AGN) and an ultraluminous X-ray source (ULX) candidate in NGC 5643. Together with archival data from Chandra, XMM-Newton and Swift-BAT, we perform a high-quality broadband spectral analysis of the AGN over two decades in energy ($\sim$0.5-100 keV). Previous X-ray observations suggested that the AGN is obscured by a Compton-thick (CT) column of obscuring gas along our line-of-sight. However, the lack of high-quality $\gtrsim$ 10 keV observations, together with the presence of a nearby X-ray luminous source, NGC 5643 X-1, had left significant uncertainties in the characterization of the nuclear spectrum. NuSTAR now enables the AGN and NGC 5643 X-1 to be separately resolved above 10 keV for the first time and allows a direct measurement of the absorbing column density toward the nucleus. The new data show that the nucleus is indeed obscured by a CT column of $N_{\rm{H}}$ $\gtrsim$ 5 $\times$ 10$^{24}$ cm$^{-2}$. The range of 2-10 keV absorption-corrected luminosity inferred from the best fitting models is $L_{2-10,\rm{int}} =$ (0.8-1.7) $\times$ 10$^{42}$ erg s$^{-1}$, consistent with that predicted from multiwavelength intrinsic luminosity indicators. We also study the NuSTAR data for NGC 5643 X-1, and show that it exhibits evidence for a spectral cut-off at energy, $E$ $\sim$ 10 keV, similar to that seen in other ULXs observed by NuSTAR. Along with the evidence for significant X-ray luminosity variations in the 3-8 keV band from 2003-2014, our results further strengthen the ULX classification of NGC 5643 X-1.

preprint2015arXiv

NuSTAR Reveals Extreme Absorption in z < 0.5 Type 2 Quasars

The intrinsic column density (NH) distribution of quasars is poorly known. At the high obscuration end of the quasar population and for redshifts z<1, the X-ray spectra can only be reliably characterized using broad-band measurements which extend to energies above 10 keV. Using the hard X-ray observatory NuSTAR, along with archival Chandra and XMM-Newton data, we study the broad-band X-ray spectra of nine optically selected (from the SDSS), candidate Compton-thick (NH > 1.5e24 cm^-2) type 2 quasars (CTQSO2s); five new NuSTAR observations are reported herein, and four have been previously published. The candidate CTQSO2s lie at z<0.5, have observed [OIII] luminosities in the range 8.4 < log (L_[OIII]/L_solar) < 9.6, and show evidence for extreme, Compton-thick absorption when indirect absorption diagnostics are considered. Amongst the nine candidate CTQSO2s, five are detected by NuSTAR in the high energy (8-24 keV) band: two are weakly detected at the ~ 3 sigma confidence level and three are strongly detected with sufficient counts for spectral modeling (>~ 90 net source counts at 8-24 keV). For these NuSTAR-detected sources direct (i.e., X-ray spectral) constraints on the intrinsic AGN properties are feasible, and we measure column densities ~2.5-1600 times higher and intrinsic (unabsorbed) X-ray luminosities ~10-70 times higher than pre-NuSTAR constraints from Chandra and XMM-Newton. Assuming the NuSTAR-detected type 2 quasars are representative of other Compton-thick candidates, we make a correction to the NH distribution for optically selected type 2 quasars as measured by Chandra and XMM-Newton for 39 objects. With this approach, we predict a Compton-thick fraction of f_CT = 36^{+14}_{-12} %, although higher fractions (up to 76%) are possible if indirect absorption diagnostics are assumed to be reliable.

preprint2015arXiv

Rapid Variability of Blazar 3C 279 during Flaring States in 2013-2014 with Joint Fermi-LAT, NuSTAR, Swift, and Ground-Based Multi-wavelength Observations

We report the results of a multi-band observing campaign on the famous blazar 3C 279 conducted during a phase of increased activity from 2013 December to 2014 April, including first observations of it with NuSTAR. The $γ$-ray emission of the source measured by Fermi-LAT showed multiple distinct flares reaching the highest flux level measured in this object since the beginning of the Fermi mission, with $F(E > 100\,{\rm MeV})$ of $10^{-5}$ photons cm$^{-2}$ s$^{-1}$, and with a flux doubling time scale as short as 2 hours. The $γ$-ray spectrum during one of the flares was very hard, with an index of $Γ_γ= 1.7 \pm 0.1$, which is rarely seen in flat spectrum radio quasars. The lack of concurrent optical variability implies a very high Compton dominance parameter $L_γ/L_{\rm syn} > 300$. Two 1-day NuSTAR observations with accompanying Swift pointings were separated by 2 weeks, probing different levels of source activity. While the 0.5$-$70 keV X-ray spectrum obtained during the first pointing, and fitted jointly with Swift-XRT is well-described by a simple power law, the second joint observation showed an unusual spectral structure: the spectrum softens by $ΔΓ_{\rm X} \simeq 0.4$ at $\sim$4 keV. Modeling the broad-band SED during this flare with the standard synchrotron plus inverse Compton model requires: (1) the location of the $γ$-ray emitting region is comparable with the broad line region radius, (2) a very hard electron energy distribution index $p \simeq 1$, (3) total jet power significantly exceeding the accretion disk luminosity $L_{\rm j}/L_{\rm d} \gtrsim 10$, and (4) extremely low jet magnetization with $L_{\rm B}/L_{\rm j} \lesssim 10^{-4}$. We also find that single-zone models that match the observed $γ$-ray and optical spectra cannot satisfactorily explain the production of X-ray emission.

preprint2015arXiv

Storm in a "Teacup": a radio-quiet quasar with ~10kpc radio-emitting bubbles and extreme gas kinematics

We present multi-frequency (1-8 GHz) VLA data, combined with VIMOS IFU data and HST imaging, of a z=0.085 radio-quiet type 2 quasar (with L(1.4GHz)~5e23 W/Hz and L(AGN)~2e45 erg/s). Due to the morphology of its emission-line region, the target (J1430+1339) has been referred to as the Teacup AGN in the literature. We identify "bubbles" of radio emission that are extended ~10-12 kpc to both the east and west of the nucleus. The edge of the brighter eastern bubble is co-spatial with an arc of luminous ionized gas. We also show that the Teacup AGN hosts a compact radio structure, located ~0.8 kpc from the core position, at the base of the eastern bubble. This radio structure is co-spatial with an ionized outflow with an observed velocity of v=-740 km/s. This is likely to correspond to a jet, or possibly a quasar wind, interacting with the interstellar medium at this position. The large-scale radio bubbles appear to be inflated by the central AGN, which indicates that the AGN can also interact with the gas on >~10 kpc scales. Our study highlights that even when a quasar is formally "radio-quiet" the radio emission can be extremely effective for observing the effects of AGN feedback.

preprint2015arXiv

The hard X-ray emission of the luminous infrared galaxy NGC 6240 as observed by NuSTAR

We present a broad-band (~0.3-70 keV) spectral and temporal analysis of NuSTAR observations of the luminous infrared galaxy NGC 6240, combined with archival Chandra, XMM-Newton and BeppoSAX data. NGC 6240 is a galaxy in a relatively early merger state with two distinct nuclei separated by ~1."5. Previous Chandra observations have resolved the two nuclei, showing that they are both active and obscured by Compton-thick material. Although they cannot be resolved by NuSTAR, thanks to the unprecedented quality of the NuSTAR data at energies >10 keV, we clearly detect, for the first time, both the primary and the reflection continuum components. The NuSTAR hard X-ray spectrum is dominated by the primary continuum piercing through an absorbing column density which is mildly optically thick to Compton scattering (tau ~ 1.2, N_H ~ 1.5 x 10^(24) cm^-2). We detect moderate hard X-ray (> 10 keV) flux variability up to 20% on short (15-20 ksec) timescales. The amplitude of the variability is maximum at ~30 keV and is likely to originate from the primary continuum of the southern nucleus. Nevertheless, the mean hard X-ray flux on longer timescales (years) is relatively constant. Moreover, the two nuclei remain Compton-thick, although we find evidence of variability of the material along the line of sight with column densities N_H <~ 2 x 10^(23) cm-2 over long (~3-15 years) timescales. The observed X-ray emission in the NuSTAR energy range is fully consistent with the sum of the best-fit models of the spatially resolved Chandra spectra of the two nuclei.

preprint2015arXiv

The NuSTAR Extragalactic Surveys: First Direct Measurements of the >10 keV X-Ray Luminosity Function for Active Galactic Nuclei at z>0.1

We present the first direct measurements of the rest-frame 10-40 keV X-ray luminosity function (XLF) of Active Galactic Nuclei (AGNs) based on a sample of 94 sources at 0.1 < z <3, selected at 8-24 keV energies from sources in the NuSTAR extragalactic survey program. Our results are consistent with the strong evolution of the AGN population seen in prior, lower-energy studies of the XLF. However, different models of the intrinsic distribution of absorption, which are used to correct for selection biases, give significantly different predictions for the total number of sources in our sample, leading to small, systematic differences in our binned estimates of the XLF. Adopting a model with a lower intrinsic fraction of Compton-thick sources and a larger population of sources with column densities N_H ~ 10^{23-24} /cm2 or a model with stronger Compton reflection component (with a relative normalization of R ~ 2 at all luminosities) can bring extrapolations of the XLF from 2-10 keV into agreement with our NuSTAR sample. Ultimately, X-ray spectral analysis of the NuSTAR sources is required to break this degeneracy between the distribution of absorbing column densities and the strength of the Compton reflection component and thus refine our measurements of the XLF. Furthermore, the models that successfully describe the high-redshift population seen by NuSTAR tend to over-predict previous, high-energy measurements of the local XLF, indicating that there is evolution of the AGN population that is not fully captured by the current models.

preprint2015arXiv

The NuSTAR Extragalactic Surveys: Initial Results and Catalog from the Extended Chandra Deep Field South

We present initial results and the source catalog from the NuSTAR survey of the Extended Chandra Deep Field South (hereafter, ECDFS) - currently the deepest contiguous component of the NuSTAR extragalactic survey program. The survey covers the full ~30 arcmin x 30 arcmin area of this field to a maximum depth of ~360 ks (~220 ks when corrected for vignetting at 3-24 keV), reaching sensitivity limits of ~1.3 x 10^-14 erg/cm2/s (3-8 keV), ~3.4 x 10^-14 erg/cm2/s (8-24 keV) and ~3.0 x 10^-14 erg/cm2/s (3-24 keV). Fifty four (54) sources are detected over the full field, although five of these are found to lie below our significance threshold once contaminating flux from neighboring (i.e., blended) sources is taken into account. Of the remaining 49 that are significant, 19 are detected in the 8-24 keV band. The 8-24 keV to 3-8 keV band ratios of the twelve sources that are detected in both bands span the range 0.39-1.7, corresponding to a photon index range of Gamma ~ 0.5-2.3, with a median photon index of 1.70 +/- 0.52. The redshifts of the 49 sources in our main sample span the range z = 0.21-2.7, and their rest-frame 10-40 keV luminosities (derived from the observed 8-24 keV fluxes) span the range L(10-40 keV) ~ (0.7-300) x 10^43 erg/s, sampling below the "knee" of the X-ray luminosity function out to z ~ 0.8-1. Finally, we identify one NuSTAR source that has neither a Chandra nor an XMM-Newton counterpart, but that shows evidence of nuclear activity at infrared wavelengths, and thus may represent a genuine, new X-ray source detected by NuSTAR in the ECDFS.

preprint2015arXiv

The NuSTAR Extragalactic Surveys: Overview and Catalog from the COSMOS Field

To provide the census of the sources contributing to the X-ray background peak above 10 keV, NuSTAR is performing extragalactic surveys using a three-tier "wedding cake" approach. We present the NuSTAR survey of the COSMOS field, the medium sensitivity and medium area tier, covering 1.7 deg2 and overlapping with both Chandra and XMM-Newton data. This survey consists of 121 observations for a total exposure of ~3 Ms. To fully exploit these data, we developed a new detection strategy, carefully tested through extensive simulations. The survey sensitivity at 20% completeness is 5.9, 2.9 and 6.4 x 10^-14 erg/cm2/s in the 3-24 keV, 3-8 keV and 8-24 keV bands, respectively. By combining detections in 3 bands, we have a sample of 91 NuSTAR sources with 10^42 -10^45.5 erg/s luminosities and redshift z=0.04-2.5. Thirty two sources are detected in the 8-24 keV band with fluxes ~100 times fainter than sources detected by Swift-BAT. Of the 91 detections, all but four are associated with a Chandra and/or XMM-Newton point-like counterpart. One source is associated with an extended lower energy X-ray source. We present the X-ray (hardness ratio and luminosity) and optical-to-X-ray properties. The observed fraction of candidate Compton-thick AGN measured from the hardness ratio is between 13-20%. We discuss the spectral properties of NuSTAR J100259+0220.6 (ID 330) at z=0.044, with the highest hardness ratio in the entire sample. The measured column density exceeds 10^24 /cm2, implying the source is Compton-thick. This source was not previously recognized as such without the >10 keV data.

preprint2014arXiv

{\it NuSTAR} Reveals an Intrinsically X-ray Weak Broad Absorption Line Quasar in the Ultraluminous Infrared Galaxy Markarian 231

We present high-energy (3--30 keV) {\it NuSTAR} observations of the nearest quasar, the ultraluminous infrared galaxy (ULIRG) Markarian 231 (Mrk 231), supplemented with new and simultaneous low-energy (0.5--8 keV) data from {\it Chandra}. The source was detected, though at much fainter levels than previously reported, likely due to contamination in the large apertures of previous non-focusing hard X-ray telescopes. The full band (0.5--30 keV) X-ray spectrum suggests the active galactic nucleus (AGN) in Mrk 231 is absorbed by a patchy and Compton-thin (N$_{\rm H} \sim1.2^{+0.3}_{-0.3}\times10^{23}$ cm$^{-2}$) column. The intrinsic X-ray luminosity (L$_{\rm 0.5-30 keV}\sim1.0\times10^{43}$ erg s$^{1}$) is extremely weak relative to the bolometric luminosity where the 2--10 keV to bolometric luminosity ratio is $\sim$0.03% compared to the typical values of 2--15%. Additionally, Mrk 231 has a low X-ray-to-optical power law slope ($α_{\rm OX}\sim-1.7$). It is a local example of a low-ionization broad absorption line (LoBAL) quasar that is intrinsically X-ray weak. The weak ionizing continuum may explain the lack of mid-infrared [O IV], [Ne V], and [Ne VI] fine-structure emission lines which are present in sources with otherwise similar AGN properties. We argue that the intrinsic X-ray weakness may be a result of the super-Eddington accretion occurring in the nucleus of this ULIRG, and may also be naturally related to the powerful wind event seen in Mrk 231, a merger remnant escaping from its dusty cocoon.

preprint2014arXiv

A deficit of ultraluminous X-ray sources in luminous infrared galaxies

We present results from a Chandra study of ultraluminous X-ray sources (ULXs) in a sample of 17 nearby (D_L<60 Mpc) luminous infrared galaxies (LIRGs), selected to have star formation rates (SFRs) in excess of 7 M_sun yr^-1 and low foreground Galactic column densities (N_H < 5*10^20 cm^-2). A total of 53 ULXs were detected and we confirm that this is a complete catalogue of ULXs for the LIRG sample. We examine the evolution of ULX spectra with luminosity by stacking the spectra of individual objects in three luminosity bins, finding a distinct change in spectral index at luminosity ~2 *10^39 erg s^-1. This may be a change in spectrum as 10 M_sun black holes transit from a ~Eddington to a super-Eddington accretion regime, and is supported by a plausible detection of partially-ionised absorption imprinted on the spectrum of the luminous ULX (L_X ~5*10^39 erg s^-1) CXOU J024238.9-000055 in NGC 1068, consistent with the highly ionised massive wind that we would expect to see driven by a super-Eddington accretion flow. This sample shows a large deficit in the number of ULXs detected per unit SFR (0.2 ULXs M_sun^-1 yr^-1). This deficit also manifests itself as a lower differential X-ray luminosity function normalisation for the LIRG sample than for samples of other star forming galaxies. We show that it is unlikely that this deficit is a purely observational effect. Part of this deficit might be attributable to the high metallicity of the LIRGs impeding the production efficiency of ULXs and/or a lag between the star formation starting and the production of ULXs; however, we argue that the evidence -- including very low N_ULX/L_FIR, and an even lower ULX incidence in the central regions of the LIRGs -- shows that the main culprit for this deficit is likely to be the high column of gas and dust in these galaxies, that fuels the high SFR but also acts to obscure many ULXs from our view.

preprint2014arXiv

A submillimeter galaxy illuminating its circumgalactic medium: Ly-alpha scattering in a cold, clumpy outflow

We report the detection at 850um of the central source in SSA22-LAB1, the archetypal Lyman-alpha Blob (LAB), a 100kpc-scale radio-quiet emission-line nebula at z=3.1. The flux density of the source, $S_{850}=4.6\pm1.1$mJy implies the presence of a galaxy, or group of galaxies, with a total luminosity of $L_{\rm IR}\approx10^{12}L_\odot$. The position of an active source at the center of a ~50kpc-radius ring of linearly polarized Ly-alpha emission detected by Hayes et al. (2011) suggests that the central source is leaking Ly-alpha photons preferentially in the plane of the sky, which undergo scattering in HI clouds at large galactocentric radius. The Ly-alpha morphology around the submillimeter detection is reminiscent of biconical outflow, and the average Ly-alpha line profiles of the two `lobes' are dominated by a red peak, expected for a resonant line emerging from a medium with a bulk velocity gradient that is outflowing relative to the line center. Taken together, these observations provide compelling evidence that the central active galaxy (or galaxies) is responsible for a large fraction of the extended Ly-alpha emission and morphology. Less clear is the history of the cold gas in the circumgalactic medium being traced by Ly-alpha: is it mainly pristine material accreting into the halo that has not yet been processed through an interstellar medium (ISM), now being blown back as it encounters an outflow, or does it mainly comprise gas that has been swept-up within the ISM and expelled from the galaxy?

preprint2014arXiv

An ALMA Survey of Submillimeter Galaxies in the Extended Chandra Deep Field South: Near-infrared morphologies and stellar sizes

We analyse HST WFC3/$H_{160}$-band observations of a sample of 48 ALMA-detected submillimeter galaxies (SMGs) in the Extended Chandra Deep Field South field, to study their stellar morphologies and sizes. We detect 79$\pm$17% of the SMGs in the $H_{160}$-band imaging with a median sensitivity of 27.8 mag, and most (80%) of the non-detections are SMGs with 870$μ$m fluxes of $S_{870} < $3 mJy. With a surface brightness limit of $μ_H \sim $26 mag arcsec$^{-2}$, we find that 82$\pm$9% of the $H_{160}$-band detected SMGs at $z =$ 1-3 appear to have disturbed morphologies, meaning they are visually classified as either irregulars or interacting systems, or both. By determining a Sérsic fit to the $H_{160}$ surface-brightness profiles we derive a median Sérsic index of $n = $1.2$\pm$0.3 and a median half-light radius of $r_e = $4.4$^{+1.1}_{-0.5}$ kpc for our SMGs at $z = $1-3. We also find significant displacements between the positions of the $H_{160}$-component and 870$μ$m emission in these systems, suggesting that the dusty star-burst regions and less-obscured stellar distribution are not co-located. We find significant differences in the sizes and the Sérsic index between our $z = $2-3 SMGs and $z \sim $2 quiescent galaxies, suggesting a major transformation of the stellar light profile is needed in the quenching processes if SMGs are progenitors of the red-and-dead $z\sim$2 galaxies. Given the short-lived nature of SMGs, we postulate that the majority of the $z = $2-3 SMGs with $S_{870} \gtrsim $2 mJy are early/mid-stage major mergers.

preprint2014arXiv

An ALMA survey of submillimetre galaxies in the Extended Chandra Deep Field South: radio properties and the far-infrared/radio correlation

We present a study of the radio properties of 870$μ$m-selected submillimetre galaxies (SMGs), observed at high resolution with ALMA in the Extended Chandra Deep Field South. From our initial sample of 76 ALMA SMGs, we detect 52 SMGs at $>3σ$ significance in VLA 1400MHz imaging, of which 35 are also detected at $>3σ$ in new 610MHz GMRT imaging. Within this sample of radio-detected SMGs, we measure a median radio spectral index $α_{610}^{1400} = -0.79 \pm 0.06$, (with inter-quartile range $α=[-1.16,-0.56]$) and investigate the far-infrared/radio correlation via the parameter $q_{\rm IR}$, the logarithmic ratio of the rest-frame 8-1000$μ$m flux and monochromatic radio flux. Our median $q_{\rm IR} = 2.56 \pm 0.05$ (inter-quartile range $q_{\rm IR}=[2.42,2.78]$) is higher than that typically seen in single-dish 870$μ$m-selected sources ($q_{\rm IR} \sim 2.4$), which may reflect the fact that our ALMA-based study is not biased to radio-bright counterparts, as previous samples were. Finally, we search for evidence that $q_{\rm IR}$ and $α$ evolve with age in a co-dependent manner, as predicted by starburst models: the data populate the predicted region of parameter space, with the stellar mass tending to increase along tracks of $q_{\rm IR}$ versus $α$ in the direction expected, providing the first observational evidence in support of these models.

preprint2014arXiv

Asymmetries in core-collapse supernovae from maps of radioactive $^{44}$Ti in CassiopeiaA

Asymmetry is required by most numerical simulations of stellar core-collapse explosions, but the form it takes differs significantly among models. The spatial distribution of radioactive 44Ti, synthesized in an exploding star near the boundary between material falling back onto the collapsing core and that ejected into the surrounding medium, directly probes the explosion asymmetries. CassiopeiaA is a young, nearby, core-collapse remnant from which 44Ti emission has previously been detected but not imaged. Asymmetries in the explosion have been indirectly inferred from a high ratio of observed 44Ti emission to estimated 56Ni emission, from optical light echoes, and from jet-like features seen in the X-ray and optical ejecta. Here we report spatial maps and spectral properties of the 44Ti in Cassiopeia A. This may explain the unexpected lack of correlation between the 44Ti and iron X-ray emission, the latter being visible only in shock-heated material. The observed spatial distribution rules out symmetric explosions even with a high level of convective mixing, as well as highly asymmetric bipolar explosions resulting from a fast-rotating progenitor. Instead, these observations provide strong evidence for the development of low-mode convective instabilities in core-collapse supernovae.

preprint2014arXiv

Kiloparsec-scale outflows are prevalent among luminous AGN: outflows and feedback in the context of the overall AGN population

We present integral field unit (IFU) observations covering the [O III]4959,5007 and H-Beta emission lines of sixteen z<0.2 type 2 active galactic nuclei (AGN). Our targets are selected from a well-constrained parent sample of 24,000 AGN so that we can place our observations into the context of the overall AGN population. Our targets are radio-quiet with star formation rates (<~[10-100] Msol/yr) that are consistent with normal star-forming galaxies. We decouple the kinematics of galaxy dynamics and mergers from outflows. We find high-velocity ionised gas (velocity widths of 600-1500 km/s and maximum velocities of <=1700 km/s) with observed spatial extents of >~(6-16) kpc in all targets and observe signatures of spherical outflows and bi-polar superbubbles. We show that our targets are representative of z<0.2, luminous (i.e., L([O III]) > 5x10^41 erg/s) type 2 AGN and that ionised outflows are not only common but also in >=70% (3 sigma confidence) of cases, they are extended over kiloparsec scales. Our study demonstrates that galaxy-wide energetic outflows are not confined to the most extreme star-forming galaxies or radio-luminous AGN; however, there may be a higher incidence of the most extreme outflow velocities in quasars hosted in ultra-luminous infrared galaxies. Both star formation and AGN activity appear to be energetically viable to drive the outflows and we find no definitive evidence that favours one process over the other. Although highly uncertain, we derive mass outflow rates (typically ~10x the SFRs), kinetic energies (~0.5-10% of L[AGN]) and momentum rates (typically >~10-20x L[AGN]/c) consistent with theoretical models that predict AGN-driven outflows play a significant role in shaping the evolution of galaxies.

preprint2014arXiv

NuSTAR and XMM-Newton Observations of Luminous, Heavily Obscured, WISE-Selected Quasars at z ~ 2

We report on a NuSTAR and XMM-Newton program that has observed a sample of three extremely luminous, heavily obscured WISE-selected AGN at z~2 in a broad X-ray band (0.1 - 79 keV). The parent sample, selected to be faint or undetected in the WISE 3.4um (W1) and 4.6um (W2) bands but bright at 12um (W3) and 22um (W4), are extremely rare, with only ~1000 so-called W1W2-dropouts across the extragalactic sky. Optical spectroscopy reveals typical redshifts of z~2 for this population, implying rest-frame mid-IR luminosities of L(6um)~6e46 erg/s and bolometric luminosities that can exceed L(bol)~1e14 L(sun). The corresponding intrinsic, unobscured hard X-ray luminosities are L(2-10)~4e45 erg/s for typical quasar templates. These are amongst the most luminous AGN known, though the optical spectra rarely show evidence of a broad-line region and the selection criteria imply heavy obscuration even at rest-frame 1.5um. We designed our X-ray observations to obtain robust detections for gas column densities N(H)<1e24 /cm2. In fact, the sources prove to be fainter than these predictions. Two of the sources were observed by both NuSTAR and XMM-Newton, with neither being detected by NuSTAR and one being faintly detected by XMM-Newton. A third source was observed only with XMM-Newton, yielding a faint detection. The X-ray data require gas column densities N(H)>1e24 /cm2, implying the sources are extremely obscured, consistent with Compton-thick, luminous quasars. The discovery of a significant population of heavily obscured, extremely luminous AGN does not conform to the standard paradigm of a receding torus, in which more luminous quasars are less likely to be obscured. If a larger sample conforms with this finding, then this suggests an additional source of obscuration for these extreme sources.

preprint2014arXiv

NuSTAR J033202-2746.8: direct constraints on the Compton reflection in a heavily obscured quasar at z~2

We report NuSTAR observations of NuSTAR J033202-2746.8, a heavily obscured, radio-loud quasar detected in the Extended Chandra Deep Field-South, the deepest layer of the NuSTAR extragalactic survey (~400 ks, at its deepest). NuSTAR J033202-2746.8 is reliably detected by NuSTAR only at E>8 keV and has a very flat spectral slope in the NuSTAR energy band (Gamma=0.55^{+0.62}_{-0.64}; 3-30 keV). Combining the NuSTAR data with extremely deep observations by Chandra and XMM-Newton (4 Ms and 3 Ms, respectively), we constrain the broad-band X-ray spectrum of NuSTAR J033202-2746.8, indicating that this source is a heavily obscured quasar (N_H=5.6^{+0.9}_{-0.8}x10^23 cm^-2) with luminosity L_{10-40 keV}~6.4x10^44 erg s^-1. Although existing optical and near-infrared (near-IR) data, as well as follow-up spectroscopy with the Keck and VLT telescopes, failed to provide a secure redshift identification for NuSTAR J033202-2746.8, we reliably constrain the redshift z=2.00+/-0.04 from the X-ray spectral features (primarily from the iron K edge). The NuSTAR spectrum shows a significant reflection component (R=0.55^{+0.44}_{-0.37}), which was not constrained by previous analyses of Chandra and XMM-Newton data alone. The measured reflection fraction is higher than the R~0 typically observed in bright radio-loud quasars such as NuSTAR J033202-2746.8, which has L_{1.4 GHz}~10^27 W Hz^-1. Constraining the spectral shape of AGN, including bright quasars, is very important for understanding the AGN population, and can have a strong impact on the modeling of the X-ray background. Our results show the importance of NuSTAR in investigating the broad-band spectral properties of quasars out to high redshift.

preprint2014arXiv

NuSTAR Observations of Heavily Obscured Quasars at z ~ 0.5

We present NuSTAR hard X-ray (3-79 keV) observations of three Type 2 quasars at z ~ 0.4-0.5, optically selected from the Sloan Digital Sky Survey (SDSS). Although the quasars show evidence for being heavily obscured Compton-thick systems on the basis of the 2-10 keV to [OIII] luminosity ratio and multiwavelength diagnostics, their X-ray absorbing column densities (N_H) are poorly known. In this analysis: (1) we study X-ray emission at >10 keV, where X-rays from the central black hole are relatively unabsorbed, in order to better constrain N_H; (2) we further characterize the physical properties of the sources through broad-band near-UV to mid-IR spectral energy distribution (SED) analyses. One of the quasars is detected with NuSTAR at >8 keV with a no-source probability of <0.1%, and its X-ray band ratio suggests near Compton-thick absorption with N_H \gtrsim 5 x 10^23 cm^-2. The other two quasars are undetected, and have low X-ray to mid-IR luminosity ratios in both the low energy (2-10 keV) and high energy (10-40 keV) X-ray regimes that are consistent with extreme, Compton-thick absorption (N_H \gtrsim 10^24 cm^-2). We find that for quasars at z ~ 0.5, NuSTAR provides a significant improvement compared to lower energy (<10 keV) Chandra and XMM-Newton observations alone, as higher column densities can now be directly constrained.

preprint2014arXiv

NuSTAR unveils a Compton-thick Type 2 quasar in Mrk 34

We present Nustar 3-40 keV observations of the optically selected Type 2 quasar (QSO2) SDSS J1034+6001 or Mrk 34. The high-quality hard X-ray spectrum and archival XMM-Newton data can be fitted self-consistently with a reflection-dominated continuum and strong Fe Kalpha fluorescence line with equivalent-width >1 keV. Prior X-ray spectral fitting below 10 keV showed the source to be consistent with being obscured by Compton-thin column densities of gas along the line-of-sight, despite evidence for much higher columns from multiwavelength data. NuSTAR now enables a direct measurement of this column, and shows that Nh lies in the Compton-thick (CT) regime. The new data also show a high intrinsic 2-10 keV luminosity of L_{2-10}~10^{44} erg/s, in contrast to previous low-energy X-ray measurements for which L_{2-10}<~10^{43} erg/s (i.e. X-ray selection below 10 keV does not pick up this source as an intrinsically luminous obscured quasar). Both the obscuring column and the intrinsic power are about an order of magnitude (or more) larger than inferred from pre-NuSTAR X-ray spectral fitting. Mrk34 is thus a 'gold standard' CT QSO2 and is the nearest non-merging system in this class, in contrast to the other local CT quasar NGC6240 which is currently undergoing a major merger coupled with strong star-formation. For typical X-ray bolometric correction factors, the accretion luminosity of Mrk34 is high enough to potentially power the total infrared luminosity. X-ray spectral fitting also shows that thermal emission related to star-formation is unlikely to drive the observed bright soft component below ~3 keV, favoring photionization instead.

preprint2014arXiv

Photometric Redshifts in the Hawaii-Hubble Deep Field-North (H-HDF-N)

We derive photometric redshifts (\zp) for sources in the entire ($\sim0.4$ deg$^2$) Hawaii-Hubble Deep Field-North (\hdfn) field with the EAzY code, based on point spread function-matched photometry of 15 broad bands from the ultraviolet (\bandu~band) to mid-infrared (IRAC 4.5 $μ$m). Our catalog consists of a total of 131,678 sources. We evaluate the \zp~quality by comparing \zp~with spectroscopic redshifts (\zs) when available, and find a value of normalized median absolute deviation \sigm$=$0.029 and an outlier fraction of 5.5\% (outliers are defined as sources having $\rm |\zp - \zs|/(1+\zs) > 0.15$) for non-X-ray sources. More specifically, we obtain \sigm$=0.024$ with 2.7\% outliers for sources brighter than $R=23$~mag, \sigm$=0.035$ with 7.4\% outliers for sources fainter than $R=23$~mag, \sigm$=$0.026 with 3.9\% outliers for sources having $z<1$, and \sigm$=$0.034 with 9.0\% outliers for sources having $z>1$. Our \zp\ quality shows an overall improvement over an earlier \zp\ work that focused only on the central \hdfn\ area. We also classify each object as star or galaxy through template spectral energy distribution fitting and complementary morphological parametrization, resulting in 4959 stars and 126,719 galaxies. Furthermore, we match our catalog with the 2~Ms {\it Chandra} Deep Field-North main \xray~catalog. For the 462 matched non-stellar \xray~sources (281 having \zs), we improve their \zp~quality by adding three additional AGN templates, achieving \sigm$=0.035$ and an outlier fraction of 12.5\%. We make our catalog publicly available presenting both photometry and \zp, and provide guidance on how to make use of our catalog.

preprint2014arXiv

The host galaxies of X-ray selected Active Galactic Nuclei to z=2.5: Structure, star-formation and their relationships from CANDELS and Herschel/PACS

We study the relationship between the structure and star-formation rate (SFR) of X-ray selected low and moderate luminosity active galactic nuclei (AGNs) in the two Chandra Deep Fields, using Hubble Space Telescope imaging from the Cosmic Assembly Near Infrared Extragalactic Legacy Survey (CANDELS) and deep far-infrared maps from the PEP+GOODS-Herschel survey. We derive detailed distributions of structural parameters and FIR luminosities from carefully constructed control samples of galaxies, which we then compare to those of the AGNs. At z~1, AGNs show slightly diskier light profiles than massive inactive (non-AGN) galaxies, as well as modestly higher levels of gross galaxy disturbance (as measured by visual signatures of interactions and clumpy structure). In contrast, at z~2, AGNs show similar levels of galaxy disturbance as inactive galaxies, but display a red central light enhancement, which may arise due to a more pronounced bulge in AGN hosts or due to extinguished nuclear light. We undertake a number of tests of these alternatives, but our results do not strongly favour one interpretation over the other. The mean SFR and its distribution among AGNs and inactive galaxies are similar at z>1.5. At z<1, however, clear and significant enhancements are seen in the SFRs of AGNs with bulge-dominated light profiles. These trends suggest an evolution in the relation between nuclear activity and host properties with redshift, towards a minor role for mergers and interactions at z>1.5.

preprint2014arXiv

The NuSTAR View of Nearby Compton-thick AGN: The Cases of NGC 424, NGC 1320 and IC 2560

We present X-ray spectral analyses for three Seyfert 2 active galactic nuclei, NGC 424, NGC 1320, and IC 2560, observed by NuSTAR in the 3-79 keV band. The high quality hard X-ray spectra allow detailed modeling of the Compton reflection component for the first time in these sources. Using quasi-simultaneous NuSTAR and Swift/XRT data, as well as archival XMM-Newton data, we find that all three nuclei are obscured by Compton-thick material with column densities in excess of ~5 x $10^{24}$ cm$^{-2}$, and that their X-ray spectra above 3 keV are dominated by reflection of the intrinsic continuum on Compton-thick material. Due to the very high obscuration, absorbed intrinsic continuum components are not formally required by the data in any of the sources. We constrain the intrinsic photon indices and the column density of the reflecting medium through the shape of the reflection spectra. Using archival multi-wavelength data we recover the intrinsic X-ray luminosities consistent with the broadband spectral energy distributions. Our results are consistent with the reflecting medium being an edge-on clumpy torus with a relatively large global covering factor and overall reflection efficiency of the order of 1%. Given the unambiguous confirmation of the Compton-thick nature of the sources, we investigate whether similar sources are likely to be missed by commonly used selection criteria for Compton-thick AGN, and explore the possibility of finding their high-redshift counterparts.

preprint2014arXiv

Weak Hard X-ray Emission from Broad Absorption Line Quasars: Evidence for Intrinsic X-ray Weakness

We report NuSTAR observations of a sample of six X-ray weak broad absorption line (BAL) quasars. These targets, at z=0.148-1.223, are among the optically brightest and most luminous BAL quasars known at z<1.3. However, their rest-frame 2 keV luminosities are 14 to >330 times weaker than expected for typical quasars. Our results from a pilot NuSTAR study of two low-redshift BAL quasars, a Chandra stacking analysis of a sample of high-redshift BAL quasars, and a NuSTAR spectral analysis of the local BAL quasar Mrk 231 have already suggested the existence of intrinsically X-ray weak BAL quasars, i.e., quasars not emitting X-rays at the level expected from their optical/UV emission. The aim of the current program is to extend the search for such extraordinary objects. Three of the six new targets are weakly detected by NuSTAR with <45 counts in the 3-24 keV band, and the other three are not detected. The hard X-ray (8-24 keV) weakness observed by NuSTAR requires Compton-thick absorption if these objects have nominal underlying X-ray emission. However, a soft stacked effective photon index (Γ~1.8) for this sample disfavors Compton-thick absorption in general. The uniform hard X-ray weakness observed by NuSTAR for this and the pilot samples selected with <10 keV weakness also suggests that the X-ray weakness is intrinsic in at least some of the targets. We conclude that the NuSTAR observations have likely discovered a significant population (>33%) of intrinsically X-ray weak objects among the BAL quasars with significantly weak <10 keV emission. We suggest that intrinsically X-ray weak quasars might be preferentially observed as BAL quasars.

preprint2013arXiv

A statistical relation between the X-ray spectral index and Eddington ratio of active galactic nuclei in deep surveys

We present an investigation into how well the properties of the accretion flow onto a supermassive black hole may be coupled to those of the overlying hot corona. To do so, we specifically measure the characteristic spectral index, Gamma, of a power-law energy distribution, over an energy range of 2 to 10 keV, for X-ray selected, broad-lined radio-quiet AGN up to z~2 in COSMOS and E-CDF-S. We test the previously reported dependence between Gamma and black hole mass, FWHM and Eddington ratio using a sample of AGN covering a broad range in these parameters based on both the Mg ii and H-alpha emission lines with the later afforded by recent near infrared spectroscopic observations using Subaru/FMOS. We calculate the Eddington ratios, lambda_Edd, for sources where a bolometric luminosity (L_Bol) has been presented in the literature, based on SED fitting, or, for sources where these data do not exist, we calculate L_Bol using a bolometric correction to the LX, derived from a relationship between the bolometric correction, and LX/L3000. From a sample of 69 X-ray bright sources (>250 counts), where Gamma can be measured with greatest precision, with an estimate of L_Bol, we find a statistically significant correlation between Gamma and lambda_Edd, which is highly significant with a chance probability of 6.59x10^-8. A statistically significant correlation between Gamma and the FWHM of the optical lines is confirmed, but at lower significance than with lambda_Edd indicating that lambda_Edd is the key parameter driving conditions in the corona. Linear regression analysis reveals that Gamma=(0.32+/-0.05)log10 lambda_Edd+(2.27+/-0.06) and Gamma=(-0.69+/-0.11)log10(FWHM/km/s)+(4.44+/-0.42). Our results on Gamma-lambda_Edd are in very good agreement with previous results. (ABRIDGED)

preprint2013arXiv

A wide search for obscured Active Galactic Nuclei using XMM-Newton and WISE

We use a combination of the XMM-Newton serendipitous X-ray survey with the optical SDSS, and the infrared WISE all-sky survey in order to check the efficiency of the low X-ray to infrared luminosity selection method in finding heavily obscured AGN. We select sources in the 2-8 keV X-ray band which have a redshift determination in the SDSS catalogue. We match this sample with the WISE catalogue, and fit the SEDs of the 2844 sources which have three, or more, photometric data-points in the infrared. We then select the heavily obscured AGN candidates by comparing their 12 micron AGN luminosity to the observed 2-10 keV X-ray luminosity and their expected intrinsic relation. With this approach we find 20 candidates, and we examine their X-ray and optical spectra. Of the 20 initial candidates, we find nine (64%; out of the 14, for which X-ray spectra could be fit) based on the X-ray spectra, and seven (78%; out of the nine detected spectroscopically in the SDSS) based on the [OIII] line fluxes. Combining all criteria, we determine the final number of heavily obscured AGN to be 12-19, and the number of Compton-thick AGN to be 2-5, showing that the method is reliable in finding obscured AGN, but not Compton-thick. However those numbers are smaller than what would be expected from X-ray background population synthesis models, which demonstrates how the optical-infrared selection and the scatter of the L_x-L_MIR relation introduced by observational constraints limit the efficiency of the method. Finally, we test popular obscured AGN selection methods based on mid-infrared colours, and find that the probability of an AGN to be selected by its mid-infrared colours increases with the X-ray luminosity. However, a selection scheme based on a relatively low X-ray luminosity and mid-infrared colours characteristic of QSOs would not select ~25% of the heavily obscured AGN of our sample. (abridged)

preprint2013arXiv

An ALMA survey of submillimeter galaxies in the Extended Chandra Deep Field South: Source catalog and multiplicity

We present an Atacama Large Millimeter/submillimeter Array (ALMA) Cycle 0 survey of 126 submillimeter sources from the LABOCA ECDFS Submillimeter Survey (LESS). Our 870 micron survey with ALMA (ALESS) has produced maps ~3X deeper and with a beam area ~200X smaller than the original LESS observations, doubling the current number of interferometrically-observed submillimeter sources. The high resolution of these maps allows us to resolve sources that were previously blended and accurately identify the origin of the submillimeter emission. We discuss the creation of the ALESS submillimeter galaxy (SMG) catalog, including the main sample of 99 SMGs and a supplementary sample of 32 SMGs. We find that at least 35% (possibly up to 50%) of the detected LABOCA sources have been resolved into multiple SMGs, and that the average number of SMGs per LESS source increases with LESS flux density. Using the (now precisely known) SMG positions, we empirically test the theoretical expectation for the uncertainty in the single-dish source positions. We also compare our catalog to the previously predicted radio/mid-infrared counterparts, finding that 45% of the ALESS SMGs were missed by this method. Our ~1.6" resolution allows us to measure a size of ~9 kpc X 5 kpc for the rest-frame ~300 um emission region in one resolved SMG, implying a star formation rate surface density of 80 M_sol yr^-1 kpc^-2, and we constrain the emission regions in the remaining SMGs to be <10 kpc. As the first statistically reliable survey of SMGs, this will provide the basis for an unbiased multiwavelength study of SMG properties.

preprint2013arXiv

An ALMA Survey of Submillimeter Galaxies in the Extended Chandra Deep Field South: The AGN Fraction and X-ray Properties of Submillimeter Galaxies

The large gas and dust reservoirs of submm galaxies (SMGs) could potentially provide ample fuel to trigger an Active Galactic Nucleus (AGN), but previous studies of the AGN fraction in SMGs have been controversial largely due to the inhomogeneity and limited angular resolution of the available submillimeter surveys. Here we set improved constraints on the AGN fraction and X-ray properties of the SMGs with ALMA and Chandra observations in the Extended Chandra Deep Field-South (E-CDF-S). This study is the first among similar works to have unambiguously identified the X-ray counterparts of SMGs; this is accomplished using the fully submm-identified, statistically reliable SMG catalog with 99 SMGs from the ALMA LABOCA E-CDF-S Submillimeter Survey (ALESS). We found 10 X-ray sources associated with SMGs (median redshift z = 2.3), of which 8 were identified as AGNs using several techniques that enable cross-checking. The other 2 X-ray detected SMGs have levels of X-ray emission that can be plausibly explained by their star-formation activity. 6 of the 8 SMG-AGNs are moderately/highly absorbed, with N_H > 10e23 cm-2. An analysis of the AGN fraction, taking into account the spatial variation of X-ray sensitivity, yields an AGN fraction of 17+16-6% for AGNs with rest-frame 0.5-8 keV absorption-corrected luminosity >7.8x10e42 erg s-1; we provide estimated AGN fractions as a function of X-ray flux and luminosity. ALMA's high angular resolution also enables direct X-ray stacking at the precise positions of SMGs for the first time, and we found 4 potential SMG-AGNs in our stacking sample.

preprint2013arXiv

Concurrent Supermassive Black Hole and Galaxy Growth: Linking Environment and Nuclear Activity in z = 2.23 H-alpha Emitters

We present results from a ~100 ks Chandra observation of the 2QZ Cluster 1004+00 structure at z = 2.23 (hereafter, 2QZ Clus). 2QZ Clus was originally identified as an overdensity of four optically-selected QSOs at z = 2.23 within a 15x15 arcmin^2 region. Narrow-band imaging in the near-IR revealed that the structure contains an additional overdensity of 22 z = 2.23 Halpha-emitting galaxies (HAEs), resulting in 23 unique z = 2.23 HAEs/QSOs. Our Chandra observations reveal that 3 HAEs in addition to the 4 QSOs harbor powerfully accreting supermassive black holes (SMBHs), with 2-10 keV luminosities of ~(8-60) x 10^43 ergs/s and X-ray spectral slopes consistent with unobscured AGN. Using a large comparison sample of 210 HAEs in Chandra-COSMOS (C-COSMOS), we find suggestive evidence that the AGN fraction increases with local HAE galaxy density. The 2QZ Clus HAEs reside in a moderately overdense environment (a factor of ~2 times over the field), and after excluding optically-selected QSOs, we find the AGN fraction is a factor of ~3.5^+3.8_-2.2 times higher than C-COSMOS HAEs in similar environments. Using stacking analyses of the Chandra data and Herschel SPIRE observations at 250 um, we respectively estimate mean SMBH accretion rates (Mdot) and star-formation rates (SFRs) for the 2QZ Clus and C-COSMOS samples. We find that the mean 2QZ Clus HAE stacked 2-10 keV luminosity is QSO-like (~6-10 x 10^43 ergs/s), and the implied Mdot/SFR \sim (1.6-3.2) x10^-3 is broadly consistent with the local MBH/M* relation and z ~ 2 X-ray selected AGN. The C-COSMOS HAEs are on average an order of magnitude less X-ray luminous and have Mdot/SFR \sim (0.2-0.4) x10^-3, comparable to z ~ 1-2 star-forming galaxies with similar mean X-ray luminosities. We estimate that a periodic QSO phase with duty cycle ~2-8% would be sufficient to bring the star-forming galaxies onto the local MBH/M* relation.

preprint2013arXiv

Narrow-line region gas kinematics of 24,264 optically-selected AGN: the radio connection

Using a sample of 24264 optically selected AGNs from the SDSS DR7 database, we characterise how the profile of the [OIII] emission line relates to bolometric luminosity (L_Bol), Eddington ratio, radio loudness, radio luminosity (L_Rad) and optical class (Sy 1/2) to determine what drives the kinematics of this kpc-scale line emitting gas. Spectral stacking is used to characterise how the average [OIII] profile changes as a function of these variables. After accounting for the correlation between L_Bol and L_Rad, we report that L_Rad has the strongest influence on the [OIII] profile, with moderate radio luminosity AGNs (log(L_Rad)=23-25 W/Hz) having the broadest [OIII] profiles. When binned according to Eddington ratio, only AGNs in our highest bin (i.e., R_Edd>0.3) show any signs of having broadened [OIII] profiles, although the small numbers of such extreme AGNs mean we cannot rule out that other processes (e.g., radio jets) are responsible for this broadening. We find no significant difference between the [OIII] profiles of broad and narrow-line Seyfert 1s. The [OIII] profiles of Sy 1 and 2 AGNs show the same trends in terms of line width, but Sy 1 AGNs display a much stronger blue wing, which we interpret as evidence of outflowing ionised gas. Using multi-component line fitting we calculate the proportions of AGNs with broad [OIII] profiles, finding that moderate radio luminosity AGNs are roughly 5 times more likely to have extremely broad [OIII] lines (FWHM>1000km/s) compared to lower L_Rad AGNs. Our results are consistent with the most disturbed gas kinematics being induced by compact radio cores (rather than powerful radio jets), although broadened [OIII] lines are also present, but much rarer, in low L_Rad systems (<10^23 W/Hz). Our catalogue of multi-component fits is freely available as an online resource.

preprint2013arXiv

New insights on the distant AGN population

Current X-ray surveys have proved to be essential tools in order to identify and study AGNs across cosmic time. However, there is evidence that the most heavily obscured AGNs are largely missing even in the deepest surveys. The search for these obscured AGNs is one of the most outstanding issues of extragalactic astronomy, since they are expected to make a major contribution to the high energy peak of the X-ray background (XRB) and might constitute a particularly active and dusty phase of black hole and galaxy evolution. Using a newly developed SED fitting technique to decompose the AGN from the star-formation emission in the infrared band we identify a sample of 17 IR bright quasars at z=1-3. For the majority of these sources the X-ray spectrum is well characterised by an absorbed power law model, revealing that ~25% of the sources are unabsorbed; the remainder of these IR quasars have moderate-to-high absorption (N_H>10^23 cm^-2), and ~25% are not detected in the X-rays and are likely to be Compton thick (N_H>10^24 cm^-2). We therefore find a much higher fraction of obscured than unobscured quasars, indicating that a large fraction of the luminous black hole accretion was very heavily obscured at z~2.

preprint2013arXiv

Reflection-dominated nuclear X-ray emission in the early-type galaxy ESO 565--G019

We present the discovery of a reflection-dominated active galactic nucleus (AGN) in the early-type radio-quiet galaxy ESO 565--G019 with Suzaku and Swift/BAT. The source X-ray spectrum below 10 keV is characteristic of other Compton-thick (CT) AGN, clearly showing an inverted continuum and prodigious fluorescence iron emission above ~3 keV. A Compton shoulder to the neutral Fe Kalpha line also appears to be present. There is evidence for long-term hard X-ray flux variability which we associate with changes in the intrinsic AGN power-law. The increasing sensitivity of ongoing and new hard X-ray surveys means that more such reflection-dominated AGN ought to be uncovered in the near future. ESO 565--G019 is hosted in an early-type galaxy whose morphology has been variously classified as either type E or type S0. Only about 20 bona fide CT-AGN have been identified in the local universe so far, and all exist in host galaxies with late Hubble types (S0 or later). CT columns of nuclear obscuring gas are uncommon in early-type galaxies in the local universe, so confirmation of the exact morphological class of ESO 565--G019 is important. Infrared photometry also shows the presence of large quantities of cool dust in the host, indicative of significant ongoing star-formation. ESO 565--G019 may be the first identified local example of minor-merger driven CT-AGN growth in an early-type host, or may be the result of interaction with its neighboring galaxy ESO 565--G018 in a wide pair.

preprint2013arXiv

The Mean Star Formation Rate of X-ray selected Active Galaxies and its Evolution from z=2.5: Results from PEP-Herschel

(Abridged) We study relationships between the SFR and the nuclear properties of X-ray selected AGNs out to z=2.5, using far-IR data in three extragalactic deep fields as part of the PACS Evolutionary Probe (PEP) program. Guided by studies of intrinsic infra-red AGN SEDs, we show that the majority of the FIR emission in AGNs is produced by cold dust heated by star-formation. We uncover characteristic redshift-dependent trends between the mean FIR luminosity (L_fir) and accretion luminosity (L_agn) of AGNs. At low AGN luminosities, accretion and SFR are uncorrelated at all redshifts, consistent with a scenario where most low-luminosity AGNs are primarily fueled by secular processes in their host galaxies. At high AGN luminosities, a significant correlation is observed between L_fir and L_agn, but only among AGNs at low and moderate redshifts (z<1). We interpret this as a signature of the increasing importance of major-mergers in driving both the growth of super-massive black holes (SMBHs) and global star-formation in their hosts at high AGN luminosities. However, we also find that the enhancement of SFR in luminous AGNs weakens or disappears at high redshifts (z>1). This suggests that the role of mergers in SMBH-galaxy co-evolution is less important at these epochs. At all redshifts, we find essentially no relationship between L_fir and nuclear obscuration across five orders of magnitude in obscuring column density, suggesting that various different mechanisms are likely to be responsible for obscuring X-rays in active galaxies. We explain our results within a scenario in which two different modes of SMBH fueling operate among low- and high-luminosity AGNs. We postulate, guided by emerging knowledge about the properties of high redshift galaxies, that the dominant mode of accretion among high-luminosity AGNs evolves with redshift.

preprint2013arXiv

The NuSTAR Extragalactic Survey: A First Sensitive Look at the High-Energy Cosmic X-ray Background Population

We report on the first ten identifications of sources serendipitously detected by the NuSTAR to provide the first sensitive census of the cosmic X-ray background (CXB) source population at >10 keV. We find that these NuSTAR-detected sources are ~100x fainter than those previously detected at >10 keV and have a broad range in redshift and luminosity (z=0.020-2.923 and L_10-40 keV~4x10^{41}-5x10^{45} erg/s); the median redshift and luminosity are z~0.7 and L_10-40 keV~3x10^{44} erg/s, respectively. We characterize these sources on the basis of broad-band ~0.5-32 keV spectroscopy, optical spectroscopy, and broad-band ultraviolet-to-mid-infrared SED analyzes. We find that the dominant source population is quasars with L_10-40 keV>10^{44} erg/s, of which ~50% are obscured with N_H>10^{22} cm^{-2}. However, none of the ten NuSTAR sources are Compton thick (N_H>10^{24} cm^{-2}) and we place a 90% confidence upper limit on the fraction of Compton-thick quasars (L_10-40 keV>10^{44} erg/s) selected at >10 keV of ~33% over the redshift range z=0.5-1.1. We jointly fitted the rest-frame ~10-40 keV data for all of the non-beamed sources with L_10-40 keV>10^{43} erg/s to constrain the average strength of reflection; we find R<1.4 for Γ=1.8, broadly consistent with that found for local AGNs observed at >10 keV. We also constrain the host galaxy masses and find a median stellar mass of ~10^{11} M_sun, a factor ~5 times higher than the median stellar mass of nearby high-energy selected AGNs, which may be at least partially driven by the order of magnitude higher X-ray luminosities of the NuSTAR sources. Within the low source-statistic limitations of our study, our results suggest that the overall properties of the NuSTAR sources are broadly similar to those of nearby high-energy selected AGNs but scaled up in luminosity and mass.

preprint2013arXiv

Weak Hard X-ray Emission from Two Broad Absorption Line Quasars Observed with NuSTAR: Compton-thick Absorption or Intrinsic X-ray Weakness?

We present NuSTAR hard X-ray observations of two X-ray weak broad absorption line (BAL) quasars, PG 1004+130 (radio loud) and PG 1700+518 (radio quiet). Many BAL quasars appear X-ray weak, probably due to absorption by the shielding gas between the nucleus and the accretion-disk wind. The two targets are among the optically brightest BAL quasars, yet they are known to be significantly X-ray weak at rest-frame 2-10 keV (16-120 times fainter than typical quasars). We would expect to obtain ~400-600 hard X-ray (>10 keV) photons with NuSTAR, provided that these photons are not significantly absorbed (NH<1E24 cm^{-2}). However, both BAL quasars are only detected in the softer NuSTAR bands (e.g., 4-20 keV) but not in its harder bands (e.g., 20-30 keV), suggesting that either the shielding gas is highly Compton-thick or the two targets are intrinsically X-ray weak. We constrain the column densities for both to be NH~7E24 cm^{-2} if the weak hard X-ray emission is caused by obscuration from the shielding gas. We discuss a few possibilities for how PG 1004+130 could have Compton-thick shielding gas without strong Fe K? line emission; dilution from jet-linked X-ray emission is one likely explanation. We also discuss the intrinsic X-ray weakness scenario based on a coronal-quenching model relevant to the shielding gas and disk wind of BAL quasars. Motivated by our NuSTAR results, we perform a Chandra stacking analysis with the Large Bright Quasar Survey BAL quasar sample and place statistical constraints upon the fraction of intrinsically X-ray weak BAL quasars; this fraction is likely 17-40%.

D. M. Alexander

What is connected

Connect this record

See the researcher in context

Building this map preview

140 published item(s)

GOODS-ALMA 2.0: Starbursts in the main sequence reveal compact star formation regulating galaxy evolution prequenching

NuSTAR Observations of Intrinsically X-ray Weak Quasar Candidates: An Obscuration-Only Scenario

Quasar Feedback Survey: Multi-phase outflows, turbulence and evidence for feedback caused by low power radio jets inclined into the galaxy disk

The Black Hole Mass Function Across Cosmic Times II. Heavy Seeds and (Super)Massive Black Holes

Warm molecular and ionized gas kinematics in the Type-2 quasar J0945+1737

GOODS-ALMA 2.0: Source catalog, number counts, and prevailing compact sizes in 1.1 mm galaxies

The post-Herschel view of intrinsic AGN emission: constructing templates for galaxy and AGN emission at IR wavelengths

$NuSTAR$ Observations of Four Nearby X-ray Faint AGN: Low Luminosity or Heavy Obscuration?

A binning-free method reveals a continuous relationship between galaxies' AGN power and offset from main sequence

Fundamental differences in the radio properties of red and blue quasars: enhanced compact AGN emission in red quasars

Fundamental differences in the radio properties of red and blue quasars: Insight from the LOFAR Two-metre Sky Survey (LoTSS)

KASHz: No evidence for ionised outflows instantaneously suppressing star formation in moderate luminosity AGN at $z$$\sim$$1.4$-$2.6$

Local AGN Survey (LASr): I. Galaxy sample, infrared colour selection and predictions for AGN within 100 Mpc

Searching for obscured AGN in z $\sim$ 2 submillimetre galaxies

The Bolometric Quasar Luminosity Function at z = 0-7

The Star-Formation Properties of the Observed and Simulated AGN Universe: BAT vs EAGLE

Do AGN really suppress star formation?

A New Population of Compton-Thick AGN Identified Using the Spectral Curvature Above 10 keV

ALMA observations of a z~3.1 Protocluster: Star Formation from Active Galactic Nuclei and Lyman-Alpha Blobs in an Overdense Environment

ALMA resolves extended star formation in high-z AGN host galaxies

First NuSTAR observations of the BL Lac - type blazar PKS~2155-304: constraints on the jet content and distribution of radiating particles

Herschel protocluster survey: A search for dusty star-forming galaxies in protoclusters at z=2-3

Kiloparsec-scale dust disks in high-redshift luminous submillimeter galaxies

The 2 Ms Chandra Deep Field-North Survey and the 250 ks Extended Chandra Deep Field-South Survey: Improved Point-Source Catalogs

The average submillimetre properties of Lyman-alpha Blobs at z=3

The Chandra COSMOS Legacy survey: overview and point source catalog

The Chandra Deep Field-South Survey: 7 Ms Source Catalogs

The NuSTAR Extragalactic Surveys: The Number Counts of Active Galactic Nuclei and the Resolved Fraction of the Cosmic X-ray Background

The SCUBA-2 Cosmology Legacy Survey: 850um maps, catalogues and number counts

VLA and ALMA Imaging of Intense, Galaxy-Wide Star Formation in z ~ 2 Galaxies

A NuSTAR Survey of Nearby Ultraluminous Infrared Galaxies

A remarkably flat relationship between the average star formation rate and AGN luminosity for distant X-ray AGN

ALMA and Herschel reveal that X-ray selected AGN and main-sequence galaxies have different star formation rate distributions

An enhanced merger fraction within the galaxy population of the SSA22 protocluster at z ~ 3.1

An Extragalactic Spectroscopic Survey of the SSA22 Field

Cold Dust Emission from X-ray AGN in the SCUBA-2 Cosmology Legacy Survey: Dependence on Luminosity, Obscuration & AGN Activity

Cosmic X-ray Surveys of Distant Active Galaxies: The Demographics, Physics, and Ecology of Growing Supermassive Black Holes

Mid-infrared luminous quasars in the GOODS-Herschel fields: a large population of heavily-obscured, Compton-thick quasars at z~2

NuSTAR Observations of the Compton-thick Active Galactic Nucleus and Ultraluminous X-ray Source Candidate in NGC 5643

NuSTAR Reveals Extreme Absorption in z < 0.5 Type 2 Quasars

Rapid Variability of Blazar 3C 279 during Flaring States in 2013-2014 with Joint Fermi-LAT, NuSTAR, Swift, and Ground-Based Multi-wavelength Observations

Storm in a "Teacup": a radio-quiet quasar with ~10kpc radio-emitting bubbles and extreme gas kinematics

The hard X-ray emission of the luminous infrared galaxy NGC 6240 as observed by NuSTAR

The NuSTAR Extragalactic Surveys: First Direct Measurements of the >10 keV X-Ray Luminosity Function for Active Galactic Nuclei at z>0.1

The NuSTAR Extragalactic Surveys: Initial Results and Catalog from the Extended Chandra Deep Field South

The NuSTAR Extragalactic Surveys: Overview and Catalog from the COSMOS Field

{\it NuSTAR} Reveals an Intrinsically X-ray Weak Broad Absorption Line Quasar in the Ultraluminous Infrared Galaxy Markarian 231

A deficit of ultraluminous X-ray sources in luminous infrared galaxies

A submillimeter galaxy illuminating its circumgalactic medium: Ly-alpha scattering in a cold, clumpy outflow

An ALMA Survey of Submillimeter Galaxies in the Extended Chandra Deep Field South: Near-infrared morphologies and stellar sizes

An ALMA survey of submillimetre galaxies in the Extended Chandra Deep Field South: radio properties and the far-infrared/radio correlation

Asymmetries in core-collapse supernovae from maps of radioactive $^{44}$Ti in CassiopeiaA

Kiloparsec-scale outflows are prevalent among luminous AGN: outflows and feedback in the context of the overall AGN population

NuSTAR and XMM-Newton Observations of Luminous, Heavily Obscured, WISE-Selected Quasars at z ~ 2

NuSTAR J033202-2746.8: direct constraints on the Compton reflection in a heavily obscured quasar at z~2

NuSTAR Observations of Heavily Obscured Quasars at z ~ 0.5

NuSTAR unveils a Compton-thick Type 2 quasar in Mrk 34

Photometric Redshifts in the Hawaii-Hubble Deep Field-North (H-HDF-N)

The host galaxies of X-ray selected Active Galactic Nuclei to z=2.5: Structure, star-formation and their relationships from CANDELS and Herschel/PACS

The NuSTAR View of Nearby Compton-thick AGN: The Cases of NGC 424, NGC 1320 and IC 2560

Weak Hard X-ray Emission from Broad Absorption Line Quasars: Evidence for Intrinsic X-ray Weakness

A statistical relation between the X-ray spectral index and Eddington ratio of active galactic nuclei in deep surveys

A wide search for obscured Active Galactic Nuclei using XMM-Newton and WISE

An ALMA survey of submillimeter galaxies in the Extended Chandra Deep Field South: Source catalog and multiplicity

An ALMA Survey of Submillimeter Galaxies in the Extended Chandra Deep Field South: The AGN Fraction and X-ray Properties of Submillimeter Galaxies

Concurrent Supermassive Black Hole and Galaxy Growth: Linking Environment and Nuclear Activity in z = 2.23 H-alpha Emitters

Narrow-line region gas kinematics of 24,264 optically-selected AGN: the radio connection

New insights on the distant AGN population

Reflection-dominated nuclear X-ray emission in the early-type galaxy ESO 565--G019

The Mean Star Formation Rate of X-ray selected Active Galaxies and its Evolution from z=2.5: Results from PEP-Herschel

The NuSTAR Extragalactic Survey: A First Sensitive Look at the High-Energy Cosmic X-ray Background Population

Weak Hard X-ray Emission from Two Broad Absorption Line Quasars Observed with NuSTAR: Compton-thick Absorption or Intrinsic X-ray Weakness?

Deep Silicate absorption features in Compton-thick AGN predominantly arise due to dust in the host galaxy

Energetic galaxy-wide outflows in high-redshift ultra-luminous infrared galaxies hosting AGN activity