Hot Dust-Obscured Galaxies (Hot DOGs) are among the most luminous galaxies in the Universe. Powered by highly obscured, possibly Compton-thick, active galactic nuclei (AGNs), Hot DOGs are characterized by SEDs that are very red in the mid-IR yet dominated by the host galaxy stellar emission in the UV and optical. An earlier study identified a sub-sample of Hot DOGs with significantly enhanced UV emission. One target, W0204-0506, was studied in detail and, based on Chandra observations, it was concluded that the enhanced emission was most likely due to either extreme unobscured star-formation (${\rm SFR}>1000~M_{\odot}~\rm yr^{-1}$) or to light from the highly obscured AGN scattered by gas or dust into our line of sight. Here, we present a follow-up study of W0204-0506 as well as two more Hot DOGs with excess UV emission. For the two new objects we obtained Chandra/ACIS-S observations, and for all three targets we obtained HST/WFC3 F555W and F160W imaging. We conclude that the excess UV emission is primarily dominated by light from the central highly obscured, hyper-luminous AGN that has been scattered into our line of sight. We cannot rule out, however, that star-formation may significantly contribute to the UV excess of W0204-0506.

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.

We present a broad-band X-ray study of the effect of superorbital periods on X-ray spectra and pulse profiles in the neutron star X-ray binaries LMC X-4 and SMC X-1. These two sources display periodic or quasi-periodic variations in luminosity on the order of tens of days which are known to be superorbital, and are attributed to warped, precessing accretion disks. Using joint NuSTAR and XMM-Newton observations that span a complete superorbital cycle, we examine the broad-band spectra of these sources and find the shape to be well described by an absorbed power law with a soft blackbody component. Changes in spectral shape and pulse profile shape are periodic with superorbital period, as expected from a precessing disk. We perform X-ray tomography using the changes in pulse profiles to model the geometry and kinematics of the inner accretion disk. Our simple beam and inner disk geometric model indicates that the long term changes in soft pulse shape and phase are consistent with reprocessed emission from a precessing inner disk.

Clustering measurements of obscured and unobscured quasars show that obscured quasars reside in more massive dark matter halos than their unobscured counterparts. These results are inconsistent with simple unified (torus) scenarios, but might be explained by models in which the distribution of obscuring material depends on Eddington ratio or galaxy stellar mass. We test these possibilities by constructing simple physical models to compare to observed AGN populations. We find that previously observed relationships between obscuration and Eddington ratio or stellar mass are not sufficient reproduce the observed quasar clustering results ($\langle \log M_{\text{halo}}/M_{\odot} \rangle = 12.94 ^{+ 0.10}_{- 0.11}$ and $\langle \log M_{\text{halo}}/M_{\odot} \rangle = 12.49 ^{+ 0.08}_{- 0.08}$ for obscured and unobscured populations, respectively) while maintaining the observed fraction of obscured quasars (30-65$\%$). This work suggests that evolutionary models, in which obscuration evolves on the typical timescale for black hole growth, are necessary to understand the observed clustering of mid-IR selected quasars.

We constrain X-ray spectral shapes for the ensemble of AGN based on the shape of the Cosmic X-ray Background (CXB). Specifically, we rule out regions of X-ray spectral parameter space that do not reproduce the CXB in the energy range 1-100 keV. The key X-ray spectral parameters are the photon index, {\Gamma}; the cutoff energy, Ecutoff; and the reflection scaling factor, R. Assuming each parameter follows a Gaussian distribution, we first explore the parameter space using a Bayesian approach and a fixed X-ray luminosity function (XLF). For {\sigma}_E = 36 keV and {\sigma}_R = 0.14, fixed at the observed values from the Swift-BAT 70-month sample, we allow <R>, <Ecutoff > and <{\Gamma}> to vary subject to reproducing the CXB. We report results for {\sigma}_{\Gamma} = 0.1-0.5. In an alternative approach, we define the parameter distributions, then forward model to fit the CXB by perturbing the XLF using a neural network. This approach allows us to rule out parameter combinations that cannot reproduce the CXB for any XLF. The marginalized conditional probabilities for the four free parameters are: <R> = 0.99^{+0.11}_{-0.26}, <Ecutoff> = 118^{+24}_{-23}, {\sigma}_{\Gamma} = 0.101^{+0.097}_{-0.001} and <{\Gamma}> = 1.9^{+0.08}_{-0.09}. We provide an interactive online tool for users to explore any combination of <Ecutoff>, {\sigma}_E, <{\Gamma}>, {\sigma}_{\Gamma}, <R> and {\sigma}_R including different distributions for each absorption bin, subject to the integral CXB constraint. The distributions observed in many AGN samples can be ruled out by our analysis, meaning these samples can not be representative of the full AGN population. The few samples that fall within the acceptable parameter space are hard X-ray-selected, commensurate with their having fewer selection biases.

Ninety per cent of baryons are located outside galaxies, either in the circumgalactic or intergalactic medium. Theory points to galactic winds as the primary source of the enriched and massive circumgalactic medium. Winds from compact starbursts have been observed to flow to distances somewhat greater than ten kiloparsecs, but the circumgalactic medium typically extends beyond a hundred kiloparsecs. Here we report optical integral field observations of the massive but compact galaxy SDSS J211824.06+001729.4. The oxygen [O II] lines at wavelengths of 3726 and 3729 angstroms reveal an ionized outflow spanning 80 by 100 square kiloparsecs, depositing metal-enriched gas at 10,000 kelvin through an hourglass-shaped nebula that resembles an evacuated and limb-brightened bipolar bubble. We also observe neutral gas phases at temperatures of less than 10,000 kelvin reaching distances of 20 kiloparsecs and velocities of around 1,500 kilometres per second. This multi-phase outflow is probably driven by bursts of star formation, consistent with theory.

Active Galactic Nuclei (AGNs) hosting disk water megamasers are well known to be obscured by large amounts of gas, likely due to the presence along the line of sight of an almost edge-on disky structure orbiting the supermassive black hole. Correcting for the high obscuration is crucial to infer parameters intrinsic to the source, like its luminosity. We present a broadband X-ray spectral analysis of a water megamaser AGN in an early merger (NGC 5765B), combining Chandra and NuSTAR data. NGC 5765B is highly Compton-thick and reflection-dominated, following the general trend among disk megamasers. Combining the exquisite black hole mass from masers with our X-ray spectroscopy, the Eddington ratio of the megamaser is estimated to be in the $2-14\%$ range, and its robustness is confirmed through SED fitting.

We present a simple semi-numerical model designed to explore black hole growth and galaxy evolution. This method builds on a previous model for black hole accretion that uses a semi-numerical galaxy formation model and universal Eddington ratio distribution to describe the full AGN population by independently connecting galaxy and AGN growth to the evolution of the host dark matter halos. We fit observed X-ray luminosity functions up to a redshift of z ~ 4, as well as investigate the evolution of the Eddington ratio distributions. We find that the Eddington ratio distribution evolves with redshift such that the slope of the low-Eddington accretion rate distribution increases with cosmic time, consistent with the behavior predicted in hydrodynamical simulations for galaxies with different gas fractions. We also find that the evolution of our average Eddington ratio is correlated with observed star formation histories, supporting a picture in which black holes and galaxies evolve together in a global sense. We further confirm the impact of luminosity limits on observed galaxy and halo properties by applying selection criteria to our fiducial model and comparing to surveys across a wide range of redshifts.

We investigate the infrared contribution from supermassive black hole activity versus host galaxy emission in the mid to far-infrared (IR) spectrum for a large sample of X-ray bright active galactic nuclei (AGN) residing in dusty, star-forming host galaxies. We select 703 AGN with L_X = 10^42-46 ergs/s at 0.1 < z < 5 from the Chandra XBootes X-ray Survey with rich multi-band observations in the optical to far-IR. This is the largest sample to date of X-ray AGN with mid and far-IR detections that uses spectral energy distribution (SED) decomposition to determine intrinsic AGN and host galaxy infrared luminosities. We determine weak or nonexistent relationships when averaging star-formation activity as a function of AGN activity, but see stronger positive trends when averaging L_X in bins of star-forming activity for AGN at low redshifts. We estimate an average dust covering factor of 33% based on infrared SEDs and bolometric AGN luminosity, corresponding to a Type 2 AGN population of roughly a third. We also see a population of AGN that challenge the inclination based unification model with individual dust covering factors that contradict the nuclear obscuration expected from observed X-ray hardness ratios. We see no strong connection between AGN fractions in the IR and corresponding total infrared, 24 um, or X-ray luminosities. The average rest-frame AGN contribution as a function of IR wavelength shows significant (~80%) contributions in the mid-IR that trail off at lambda > 30 um. Additionally, we provide a relation between observed L_X and pure AGN IR output for high-z AGN allowing future studies to estimate AGN infrared contribution using only observed X-ray flux density estimates.

The cosmic X-ray background (CXB), which peaks at an energy of ~30 keV, is produced primarily by emission from accreting supermassive black holes (SMBHs). The CXB therefore serves as a constraint on the integrated SMBH growth in the Universe and the accretion physics and obscuration in active galactic nuclei (AGNs). This paper gives an overview of recent progress in understanding the high-energy (>~10 keV) X-ray emission from AGNs and the synthesis of the CXB, with an emphasis on results from NASA's NuSTAR hard X-ray mission. We then discuss remaining challenges and open questions regarding the nature of AGN obscuration and AGN physics. Finally, we highlight the exciting opportunities for a next-generation, high-resolution hard X-ray mission to achieve the long-standing goal of resolving and characterizing the vast majority of the accreting SMBHs that produce the CXB.

We present an analysis of a sample of 69 local obscured Swift/Burst Alert Telescope active galactic nuclei (AGNs) with X-ray spectra from NuSTAR and infrared (IR) spectral energy distributions from Herschel and WISE. We combine this X-ray and IR phenomenological modeling and find a significant correlation between reflected hard X-ray emission and IR AGN emission, with suggestive indications that this correlation may be stronger than the one between intrinsic hard X-ray and IR emissions. This relation between the IR and reflected X-ray emission suggests that both are the result of the processing of intrinsic emission from the corona and accretion disk by the same structure. We explore the resulting implications on the underlying distribution of covering fraction for all AGNs, by generating mock observables for the reflection parameter and IR luminosity ratio using empirical relations found for the covering fraction with each quantity. We find that the observed distributions of the reflection parameter and IR-to-X-ray ratio are reproduced with broad distributions centered around covering fractions of at least ~40%-50%, whereas narrower distributions match our observations only when centered around covering fractions of ~70%-80%. Our results are consistent with both independent estimates of the covering fractions of individual objects and the typical covering fraction obtained on the basis of obscured fractions for samples of AGNs. These results suggest that the level of reprocessing in AGNs, including X-ray reflection, is related in a relatively straightforward way to the geometry of the obscuring material.

The discoveries made over the past 20 years by Chandra and XMM-Newton surveys in conjunction with multiwavelength imaging and spectroscopic data available in the same fields have significantly changed the view of the supermassive black hole (SMBH) and galaxy connection. These discoveries have opened up several exciting questions that are beyond the capabilities of current X-ray telescopes and will need to be addressed by observatories in the next two decades. As new observatories peer into the early Universe, we will begin to understand the physics and demographics of SMBH infancy (at $z>6$) and investigate the influence of their accretion on the formation of the first galaxies ($\S$ 2.1). We will also be able to understand the accretion and evolution over the cosmic history (at $z\sim$1-6) of the full population of black holes in galaxies, including low accretion rate, heavily obscured AGNs at luminosities beyond the reach of current X-ray surveys ($\S$2.2 and $\S$2.3), enabling us to resolve the connection between SMBH growth and their environment.

Galaxies grow their supermassive black holes in concert with their stars, although the relationship between these major galactic components is poorly understood. Observations of the cosmic growth of stars and black holes in galaxies suffer from disjoint samples and the strong effects of dust attenuation. The thermal infrared holds incredible potential for simultaneously measuring both the star formation and black hole accretion rates in large samples of galaxies covering a wide range of physical conditions. Spitzer demonstrated this potential at low redshift, and by observing some of the most luminous galaxies at z~2. JWST will apply these methods to normal galaxies at these epochs, but will not be able to generate large spectroscopic samples or access the thermal infrared at high-redshift. An order of magnitude gap in our wavelength coverage will persist between JWST and ALMA. A large, cold infrared telescope can fill this gap to determine when (in cosmic time), and where (within the cosmic web), stars and black holes co-evolve, by measuring these processes simultaneously in statistically complete and unbiased samples of galaxies to z>8. A next-generation radio interferometer will have the resolution and sensitivity to measure star-formation and nuclear accretion in even the dustiest galaxies. Together, the thermal infrared and radio can uniquely determine how stars and supermassive blackholes co-evolve in galaxies over cosmic time.

We measure the average deflection of cosmic microwave background photons by quasars at $\langle z \rangle =1.7$. Our sample is selected from the Sloan Digital Sky Survey to cover the redshift range $0.9\leq z\leq2.2$ with absolute i-band magnitudes of $M_i\leq-24$ (K-corrected to z=2). A stack of nearly 200,000 targets reveals an 8$\sigma$ detection of Planck's estimate of the lensing convergence towards the quasars. We fit the signal with a model comprising a Navarro-Frenk-White density profile and a 2-halo term accounting for correlated large scale structure, which dominates the observed signal. The best-fitting model is described by an average halo mass $\log_{10}(M_{\rm h}/h^{-1}M_\odot)=12.6\pm0.2$ and linear bias $b=2.7\pm0.3$ at $z=1.7$, in excellent agreement with clustering studies. We also report of a hint, at a 90% confidence level, of a correlation between the convergence amplitude and luminosity, indicating that quasars brighter than $M_i\lesssim -26$ reside in halos of typical mass ${M_{\rm h}\approx 10^{13}\,h^{-1}M_\odot}$, scaling roughly as ${M_{\rm h}\propto L_{\rm opt}^{3/4}}$ at ${M_i\lesssim-24}$, in good agreement with physically-motivated quasar demography models. Although we acknowledge this luminosity dependence is a marginal result, the observed $M_{\rm h}$-$L_{\rm opt}$ relationship could be interpreted as a reflection of the cutoff in the distribution of black hole accretion rates towards high Eddington ratios: the weak trend of $M_{\rm h}$ with $L_{\rm opt}$ observed at low luminosity becomes stronger for the most powerful quasars, which tend to be accreting close to the Eddington limit.

A primary aim of the Nuclear Spectroscopic Telescope Array (NuSTAR) mission is to find and characterize heavily obscured Active Galactic Nuclei (AGNs). Based on mid-infrared photometry from the Wide-Field Infrared Survey Explorer (WISE) and optical photometry from the Sloan Digital Sky Surveys, we have selected a large population of luminous obscured AGN (i.e., obscured quasars). Here we report NuSTAR observations of four WISE-selected heavily obscured quasars for which we have optical spectroscopy from the Southern African Large Telescope and W. M. Keck Observatory. Optical diagnostics confirm that all four targets are AGNs. With NuSTAR hard X-ray observations, three of the four objects are undetected, while the fourth has a marginal detection. We confirm that these objects have observed hard X-ray (10-40 keV) luminosities at or below ~ 10^43 erg s^-1. We compare X-ray and IR luminosities to obtain estimates of the hydrogen column densities (N_H) based on the suppression of the hard X-ray emission. We estimate N_H of these quasars to be at or larger than 10^25 cm^-2, confirming that WISE and optical selection can identify very heavily obscured quasars that may be missed in X-ray surveys, and do not contribute significantly to the cosmic X-ray background. From the optical Balmer decrements, we found that our three extreme obscured targets lie in highly reddened host environments. This galactic extinction is not adequate to explain the more obscured AGN, but it may imply a different scale of obscuration in the galaxy.

We present Atacama Large Millimeter/submillimeter Array observations of a massive (M_stars~10^11 M_Sun) compact (r_e,UV~100 pc) merger remnant at z=0.66 that is driving a 1000 km/s outflow of cool gas, with no observational trace of an active galactic nucleus (AGN). We resolve molecular gas on scales of approximately 1-2 kpc, and our main finding is the discovery of a wing of blueshifted CO(2-1) emission out to -1000 km/s relative to the stars. We argue that this is the molecular component of a multiphase outflow, expelled from the central starburst within the past 5 Myr through stellar feedback, although we cannot rule out previous AGN activity as a launching mechanism. If the latter is true, then this is an example of a relic multiphase AGN outflow. We estimate a molecular mass outflow rate of approximately 300 M_Sun/yr, or about one third of the 10 Myr-averaged star formation rate. This system epitomizes the multiphase 'blowout' episode following a dissipational major merger - a process that has violently quenched central star formation and supermassive black hole growth.

We report on the results of active galactic nuclei (AGNs) detection by NuSTAR performed in three extragalactic survey fields (COSMOS, UDS, ECDFS) in three hard bands, namely H1 (8-16 keV), H2 (16-24 keV) and VH (35-55 keV). The aggregated area of the surveys is $\sim 2.7$ deg$^2$. While a large number of sources is detected in the H1 band (72 at the $97\%$ level of reliability), the H2 band directly probing close to the peak of the Cosmic X-ray Background (CXB) returns four significant detections, and two tentative, although not significant, detections are found in the VH band. All the sources detected above 16 keV are also detected at lower energies. We compute the integral number counts for sources in such bands, which show broad consistency with population synthesis models of the CXB. We furthermore identify two Compton-thick AGNs, one in the COSMOS field, associated with a hard and faint Chandra source, and one in the UDS field, never detected in the X-ray band before. Both sources are at the same redshift $z \sim 1.25$, which shifts their Compton-hump into the H1 band, and were previously missed in the usually employed NuSTAR bands, confirming the potential of using the H1 band to discover obscured AGNs at $z > 1$ in deep surveys.

We present a sample of accreting supermassive black holes (SMBHs) in dwarf galaxies at $z<1$. We identify dwarf galaxies in the NEWFIRM Medium Band Survey with stellar masses $M_{\star}<3\times 10^{9} M_{\odot}$ that have spectroscopic redshifts from the DEEP2 survey and lie within the region covered by deep (flux limit of $\sim 5\times 10^{-17} - 6\times 10^{-16} \rm{erg \ cm}^{-2} \ \rm{s}^{-1}$) archival Chandra X-ray data. From our sample of $605$ dwarf galaxies, $10$ exhibit X-ray emission consistent with that arising from AGN activity. If black hole mass scales roughly with stellar mass, then we expect that these AGN are powered by SMBHs with masses of $\sim 10^5-10^6 \ M_{\odot}$ and typical Eddington ratios $\sim 5\%$. Furthermore, we find an AGN fraction consistent with extrapolations of other searches of $\sim 0.6-3\%$ for $10^9 \ M_{\odot} \leq M_{\star} \leq 3\times 10^{9} \ M_{\odot}$ and $0.1<z<0.6$. Our AGN fraction is in good agreement with a semi-analytic model, suggesting that as we search larger volumes we may use comparisons between observed AGN fractions and models to understand seeding mechanisms in the early universe.

Two \textit{NuSTAR} observations of the luminous X-ray pulsar LMC X-4 in October and November 2015 captured several bright accretion flares from this source, which has a long history of stable pulse and superorbital behavior. We present a timing analysis of these data in which we detect a rapid pulse turn-on in association with the accretion flares, during which the source reaches super-Eddington luminosities. Pulsations, which are normally seen from this source, are found to only occur for approximately one hour before and during the bright flares. Beyond one hour before and after the flares, we find pulsations to be weak or nonexistent, with fractional RMS amplitudes of less than 0.05. At the onset of the flare, the pulse profiles exhibit a phase shift of 0.25 cycles that could be associated with a change in the emission geometry. This increase in pulse strength occurring well before the flare cannot be explained by the propeller effect, and potentially offers a connection between the magnetic properties of pulsars that accrete close to their Eddington limits and ULX pulsars.

Active Galactic Nuclei (AGN) are powered by the accretion of material onto a supermassive black hole (SMBH), and are among the most luminous objects in the Universe. However, the huge radiative power of most AGN cannot be seen directly, as the accretion is "hidden" behind gas and dust that absorbs many of the characteristic observational signatures. This obscuration presents an important challenge for uncovering the complete AGN population and understanding the cosmic evolution of SMBHs. In this review we describe a broad range of multi-wavelength techniques that are currently employed to identify obscured AGN, and assess the reliability and completeness of each technique. We follow with a discussion of the demographics of obscured AGN activity, explore the nature and physical scales of the obscuring material, and assess the implications of obscured AGN for observational cosmology. We conclude with an outline of the prospects for future progress from both observations and theoretical models, and highlight some of the key outstanding questions.

We present a study of X-ray source populations in M87, the cD galaxy of the Virgo cluster, using 12 archival {\it Chandra} observations with a total exposure of $\sim$680 ks spanning about a decade. A total of 346 point-like sources are detected down to a limiting 0.5--8 keV luminosity of $4\times10^{37}{\rm~erg~s^{-1}}$ and out to a galactocentric radius of $\sim$40 kpc. We cross-correlate the X-ray sources with published catalogs of globular clusters (GCs), derived from the ACS Virgo Cluster Survey and the Next Generation Virgo Cluster Survey. This results in 122 matches, making it one of the largest samples of GC-hosting X-ray sources in an external galaxy. These sources, most likely low-mass X-ray binaries (LMXBs), correspond to $\sim$5\% of all known GCs within the {\it Chandra} field-of-view. Conversely, $\sim$50\% of the detected X-ray sources are found in a GC. Moreover, red (metal-rich) GCs are $\sim$2.2 times more likely to host an X-ray source than blue (metal-poor) GCs. We also examine 76 currently known ultra-compact dwarf galaxies around M87, but find no significant X-ray counterparts. After statistically accounting for the cosmic X-ray background, we identify $\sim$110 field-LMXBs. The GC-LMXBs and field-LMXBs differ in their luminosity function and radial distribution, which indicates that the latter cannot be primarily originated from GCs. Using another set of deep {\it Chandra} observations toward $\sim$100 kpc northwest of the M87 center, we statistically constrain the abundance of field-LMXBs in the stellar halo, which is consistent with that found in the central region. We also identify 40 variable X-ray sources, among which one source is likely a black hole binary residing in a GC.

We explore the kinematics of ionized gas via the [O III] $\lambda$5007 emission lines in active galactic nuclei (AGN) selected on the basis of their mid-infrared (IR) emission, and split into obscured and unobscured populations based on their optical-IR colors. After correcting for differences in redshift distributions, we provide composite spectra of spectroscopically and photometrically defined obscured/Type 2 and unobscured/Type 1 AGN from 3500 to 7000 \AA. The IR-selected obscured sources contain a mixture of narrow-lined Type 2 AGN and intermediate sources that have broad H$\alpha$ emission and significantly narrower H$\beta$. Using both [OIII] luminosities and AGN luminosities derived from optical-IR spectral energy distribution fitting, we find evidence for enhanced large-scale obscuration in the obscured sources. In matched bins of luminosity we find that the obscured population typically has broader, more blueshifted \OIII\ emission than in the unobscured sample, suggestive of more powerful AGN-driven outflows. This trend is not seen in spectroscopically classified samples, and is unlikely to be entirely explained by orientation effects. In addition, outflow velocities increase from small to moderate AGN $E(B-V)$ values, before flattening out (as traced by FWHM) and even decreasing (as traced by blueshift). While difficult to fully interpret in a single physical model, due to both the averaging over populations and the spatially-averaged spectra, these results agree with previous findings that simple geometric unification models are insufficient for the IR-selected AGN population, and may fit into an evolutionary model for obscured and unobscured AGN.

The lack of a strong correlation between AGN X-ray luminosity ($L_X$; a proxy for AGN power) and the star formation rate (SFR) of their host galaxies has recently been attributed to stochastic AGN variability. Studies using population synthesis models have incorporated this by assuming a broad, universal (i.e. does not depend on the host galaxy properties) probability distribution for AGN specific X-ray luminosities (i.e. the ratio of $L_X$ to host stellar mass; a common proxy for Eddington ratio). However, recent studies have demonstrated that this universal Eddington ratio distribution fails to reproduce the observed X-ray luminosity functions beyond z$\sim$1.2. Furthermore, empirical studies have recently shown that the Eddington ratio distribution may instead depend upon host galaxy properties, such as SFR and/or stellar mass. To investigate this further we develop a population synthesis model in which the Eddington ratio distribution is different for star-forming and quiescent host galaxies. We show that, although this model is able to reproduce the observed X-ray luminosity functions out to z$\sim$2, it fails to simultaneously reproduce the observed flat relationship between SFR and X-ray luminosity. We can solve this, however, by incorporating a mass dependency in the AGN Eddington ratio distribution for star-forming host galaxies. Overall, our models indicate that a relative suppression of low Eddington ratios ($\lambda_{\rm Edd}\lesssim$0.1) in lower mass galaxies (M<$10^{10-11}$Msun) is required to reproduce both the observed X-ray luminosity functions and the observed flat SFR/X-ray relationship.

We present the results and the source catalog of the NuSTAR survey in the UKIDSS Ultra Deep Survey (UDS) field, bridging the gap in depth and area between NuSTAR's ECDFS and COSMOS surveys. The survey covers a $\sim 0.6$ deg$^2$ area of the field for a total observing time of $\sim$ 1.75 Ms, to a half-area depth of $\sim$ 155 ks corrected for vignetting at $3-24$ keV, and reaching sensitivity limits at half-area in the full ($3-24$ keV), soft ($3-8$ keV) and hard ($8-24$ keV) bands of $2.2 \times 10^{-14}$ erg cm$^{-2}$ s$^{-1}$, $1.0 \times 10^{-14}$ erg cm$^{-2}$ s$^{-1}$, and $2.7 \times 10^{-14}$ erg cm$^{-2}$ s$^{-1}$, respectively. A total of 67 sources are detected in at least one of the three bands, 56 of which have a robust optical redshift with a median of $\langle z\rangle \sim 1.1$. Through a broadband ($0.5-24$ keV) spectral analysis of the whole sample combined with the NuSTAR hardness ratios, we compute the observed Compton-thick (CT; $N_{\rm H} > 10^{24}$ cm$^{-2}$) fraction. Taking into account the uncertainties on each $N_{\rm H}$ measurement, the final number of CT sources is $6.8\pm1.2$. This corresponds to an observed CT fraction of $(11.5\pm2.0)\%$, providing a robust lower limit to the intrinsic fraction of CT AGN and placing constraints on cosmic X-ray background synthesis models.

We present sensitive 870$\mu$m continuum measurements from our ALMA programmes of 114 X-ray selected AGN in the CDF-S and COSMOS fields. We use these observations in combination with data from Spitzer and Herschel to construct a sample of 86 X-ray selected AGN, 63 with ALMA constraints at $z=1.5-3.2$ with stellar mass $>2\times10^{10}M_{\odot}$. We constructed broad-band spectral energy distributions in the infrared band (8 - 1000$\mu$m) and constrain star-formation rates (SFRs) uncontaminated by the AGN. Using a hierarchical Bayesian method that takes into account the information from upper limits, we fit SFR and specific SFR (sSFR) distributions. We explore these distributions as a function of both X-ray luminosity and stellar mass. We compare our measurements to two versions of the EAGLE hydrodynamical simulations: the reference model with AGN feedback and the model without AGN. We find good agreement between the observations and that predicted by the EAGLE reference model for the modes and widths of the sSFR distributions as a function of both X-ray luminosity and stellar mass; however, we found that the EAGLE model without AGN feedback predicts a significantly narrower width when compared to the data. Overall, from the combination of the observations with the model predictions, we conclude that (1) even with AGN feedback, we expect no strong relationship between the sSFR distribution parameters and instantaneous AGN luminosity and (2) a signature of AGN feedback is a broad distribution of sSFRs for all galaxies (not just those hosting an AGN) with stellar masses above $\approx 10^{10}$M$_{\odot}$.

We present an X-ray spectral and timing analysis of two $NuSTAR$ observations of the transient Be X-ray binary SAX J2103.5+4545 during its April 2016 outburst, which was characterized by the highest flux since $NuSTAR$'s launch. These observations provide detailed hard X-ray spectra of this source during its bright precursor flare and subsequent fainter regular outburst for the first time. In this work, we model the phase-averaged spectra for these observations with a negative and positive power law with an exponential cut-off (NPEX) model and compare the pulse profiles at different flux states. We found that the broad-band pulse profile changes from a three peaked pulse in the first observation to a two peaked pulse in the second observation, and that each of the pulse peaks has some energy dependence. We also perform pulse-phase spectroscopy and fit phase-resolved spectra with NPEX to evaluate how spectral parameters change with pulse phase. We find that while the continuum parameters are mostly constant with pulse phase, a weak absorption feature at ~12 keV that might, with further study, be classified as a cyclotron line, does show strong pulse phase dependence.

Hot Dust-Obscured Galaxies (Hot DOGs) are hyperluminous ($L_{\mathrm{8-1000\,\mu m}}>10^{13}\,\mathrm{L_\odot}$) infrared galaxies with extremely high (up to hundreds of K) dust temperatures. The sources powering both their extremely high luminosities and dust temperatures are thought to be deeply buried and rapidly accreting supermassive black holes (SMBHs). Hot DOGs could therefore represent a key evolutionary phase in which the SMBH growth peaks. X-ray observations can be used to study their obscuration levels and luminosities. In this work, we present the X-ray properties of the 20 most-luminous ($L_{\mathrm{bol}}\gtrsim10^{14}\, L_\odot$) known Hot DOGs at $z=2-4.6$. Five of them are covered by long-exposure ($10-70$ ks) Chandra and XMM-Newton observations, with three being X-ray detected, and we study their individual properties. One of these sources (W0116$-$0505) is a Compton-thick candidate, with column density $N_H=(1.0-1.5)\times10^{24}\,\mathrm{cm^{-2}}$ derived from X-ray spectral fitting. The remaining 15 Hot DOGs have been targeted by a Chandra snapshot (3.1 ks) survey. None of these 15 is individually detected; therefore we applied a stacking analysis to investigate their average emission. From hardness-ratio analysis, we constrained the average obscuring column density and intrinsic luminosity to be log$N_H\,\mathrm{[cm^{-2}]}>23.5$ and $L_X\gtrsim10^{44}\,\mathrm{erg\,cm^{-2}\,s^{-1}}$, which are consistent with results for individually detected sources. We also investigated the $L_X-L_{6\mu\mathrm{m}}$ and $L_X-L_{bol}$ relations, finding hints that Hot DOGs are typically X-ray weaker than expected, although larger samples of luminous obscured QSOs are needed to derive solid conclusions.

We present the results of an optical spectroscopic survey of 46 heavily obscured quasar candidates. Objects are selected using their mid-infrared (mid-IR) colours and magnitudes from the Wide-Field Infrared Survey Explorer (WISE) and their optical magnitudes from the Sloan Digital Sky Survey (SDSS). Candidate Active Galactic Nuclei (AGNs) are selected to have mid-IR colours indicative of quasar activity and lie in a region of mid-IR colour space outside previously published X-ray based selection regions. We obtain optical spectra for our sample using the Robert Stobie Spectrograph on the Southern African Large Telescope. Thirty objects (65%) have identifiable emission lines, allowing for the determination of spectroscopic redshifts. Other than one object at $z\sim2.6$, candidates have moderate redshifts ranging from $z=0.1$ to $0.8$ with a median of 0.3. Twenty-one (70%) of our objects with identified redshift (46% of the whole sample) are identified as AGNs through common optical diagnostics. We model the spectral energy distributions of our sample and found that all require a strong AGN component, with an average intrinsic AGN fraction at 8$\,\mu$m of 0.91. Additionally, the fits require large extinction coefficients with an average $E(B-V)_\textrm{AGN} = 17.8$ (average $A(V)_\textrm{AGN} = 53.4$). By focusing on the area outside traditional mid-IR photometric cuts, we are able to capture and characterise a population of deeply buried quasars that were previously unattainable through X-ray surveys alone.

We present the variability and time lag measurements of PG 0934+013 based on a photometric and spectroscopic monitoring campaign over a two year period. We obtained 46 epochs of data from the spectroscopic campaign, which was carried out using the Southern African Large Telescope with $\sim$1 week cadence over two sets of 4 month-long observing period, while we obtained 80 epochs of \textit{B}-band imaging data using a few 1-m class telescopes. Due to the seven month gap between the two observing periods, we separately measured the time lags of broad emission lines including H$\beta$, by comparing the emission line light curve with the \textit{B}-band continuum light curve using the cross-correlation function techniques. We determined the H$\beta$ lag, $\tau_{\rm cent} = 8.46^{+2.08}_{-2.14}$ days in the observed-frame based on Year 2 data, while the time lag from Year 1 data was not reliably determined. Using the rms spectrum of Year 2 data, we measured the \Hb\ line dispersion \sigmaline = 668 $\pm$ 44 \kms\ after correcting for the spectral resolution. Adopting a virial factor f = 4.47 from Woo et al. 2015, we determined the black hole mass M$_{BH}$ = $3.13 ^{+0.91} _{-0.93} \times 10^{6}$ \msun based on the \Hb\ time lag and velocity.

We present observed mid-infrared and optical colors and composite spectral energy distributions (SEDs) of type 1 (broad-line) and 2 (narrow-line) quasars selected from Sloan Digital Sky Survey (SDSS) spectroscopy. A significant fraction of powerful quasars are obscured by dust, and are difficult to detect in optical photometric or spectroscopic surveys. However these may be more easily identified on the basis of mid-infrared (MIR) colors and SEDs. Using samples of SDSS type 1 type 2 matched in redshift and [OIII] luminosity, we produce composite rest-frame 0.2-15 micron SEDs based on SDSS, UKIDSS, and Wide-Field Infrared Survey Explorer (WISE) photometry and perform model fits using simple galaxy and quasar SED templates. The SEDs of type 1 and 2 quasars are remarkably similar, with the differences explained primarily by the extinction of the quasar component in the type 2 systems. For both types of quasar, the flux of the AGN relative to the host galaxy increases with AGN luminosity (L_[OIII]) and redder observed MIR color, but we find only weak dependencies of the composite SEDs on mechanical jet power as determined through radio luminosity. We conclude that luminous quasars can be effectively selected using simple MIR color criteria similar to those identified previously (W1-W2 > 0.7 [Vega]), although these criteria miss many heavily obscured objects. Obscured quasars can be further identified based on optical-IR colors (for example, (u-W3 [AB]) > 1.4(W1-W2 [Vega])+3.2). These results illustrate the power of large statistical studies of obscured quasars selected on the basis of mid-IR and optical photometry.

Several recent studies have reported different intrinsic correlations between the AGN mid-IR luminosity ($L_{MIR}$) and the rest-frame 2-10 keV luminosity ($L_{X}$) for luminous quasars. To understand the origin of the difference in the observed $L_{X}-L_{MIR}$ relations, we study a sample of 3,247 spectroscopically confirmed type 1 AGNs collected from Bo\"{o}tes, XMM-COSMOS, XMM-XXL-North, and the SDSS quasars in the Swift/XRT footprint spanning over four orders of magnitude in luminosity. We carefully examine how different observational constraints impact the observed $L_{X}-L_{MIR}$ relations, including the inclusion of X-ray non-detected objects, possible X-ray absorption in type 1 AGNs, X-ray flux limits, and star formation contamination. We find that the primary factor driving the different $L_{X}-L_{MIR}$ relations reported in the literature is the X-ray flux limits for different studies. When taking these effects into account, we find that the X-ray luminosity and mid-IR luminosity (measured at rest-frame $6\mu m$, or $L_{6\mu m}$) of our sample of type 1 AGNs follow a bilinear relation in the log-log plane: $\log L_X =(0.84\pm0.03)\times\log L_{6\mu m}/10^{45}{\rm erg\;s^{-1}} + (44.60\pm0.01)$ for $L_{6\mu m} < 10^{44.79}{\rm erg\;s^{-1}} $, and $\log L_X = (0.40\pm0.03)\times\log L_{6\mu m}/10^{45}{\rm erg\;s^{-1}} +(44.51\pm0.01)$ for $L_{6\mu m} \geq 10^{44.79}{\rm erg\;s^{-1}} $. This suggests that the luminous type 1 quasars have a shallower $L_{X}-L_{MIR}$ correlation than the approximately linear relations found in local Seyfert galaxies. This result is consistent with previous studies reporting a luminosity-dependent $L_{X}-L_{MIR}$ relation, and implies that assuming a linear $L_{X}-L_{MIR}$ relation to infer the neutral gas column density for X-ray absorption might overestimate the column densities in luminous quasars.

Active Galactic Nuclei (AGN) are energetic astrophysical sources powered by accretion onto supermassive black holes in galaxies, and present unique observational signatures that cover the full electromagnetic spectrum over more than twenty orders of magnitude in frequency. The rich phenomenology of AGN has resulted in a large number of different "flavours" in the literature that now comprise a complex and confusing AGN "zoo". It is increasingly clear that these classifications are only partially related to intrinsic differences between AGN, and primarily reflect variations in a relatively small number of astrophysical parameters as well the method by which each class of AGN is selected. Taken together, observations in different electromagnetic bands as well as variations over time provide complementary windows on the physics of different sub-structures in the AGN. In this review, we present an overview of AGN multi-wavelength properties with the aim of painting their "big picture" through observations in each electromagnetic band from radio to gamma-rays as well as AGN variability. We address what we can learn from each observational method, the impact of selection effects, the physics behind the emission at each wavelength, and the potential for future studies. To conclude we use these observations to piece together the basic architecture of AGN, discuss our current understanding of unification models, and highlight some open questions that present opportunities for future observational and theoretical progress.

In studies of the connection between active galactic nuclei (AGN) and their host galaxies there is widespread disagreement on some key aspects stemming largely from a lack of understanding of the nature of the full underlying AGN population. Recent attempts to probe this connection utilize both observations and simulations to correct for a missed population, but presently are limited by intrinsic biases and complicated models. We take a simple simulation for galaxy evolution and add a new prescription for AGN activity to connect galaxy growth to dark matter halo properties and AGN activity to star formation. We explicitly model selection effects to produce an "observed" AGN population for comparison with observations and empirically motivated models of the local universe. This allows us to bypass the difficulties inherent in many models which attempt to infer the AGN population by inverting selection effects. We investigate the impact of selecting AGN based on thresholds in luminosity or Eddington ratio on the "observed" AGN population. By limiting our model AGN sample in luminosity, we are able to recreate the observed local AGN luminosity function and specific star formation-stellar mass distribution, and show that using an Eddington ratio threshold introduces less bias into the sample by selecting the full range of growing black holes, despite the challenge of selecting low mass black holes. We find that selecting AGN using these various thresholds yield samples with different AGN host galaxy properties.

Recent work has found evidence for a difference in the bias and dark matter halo masses of WISE-selected obscured and unobscured quasars, implying a distinction between these populations beyond random line-of-sight effects. However, the significance of this difference in the most up-to-date measurements is relatively weak, at $\sim$2$\sigma$ for individual measurements but bolstered by agreement from different techniques, including angular clustering and cross-correlations with cosmic microwave background (CMB) lensing maps. Here, we expand the footprint of previous work, aiming to improve the precision of both methods. In this larger area we correct for position dependent selection effects, in particular fluctuations of the WISE-selected quasar density as a function of Galactic latitude. We also measure the cross-correlation of the obscured and unobscured samples and confirm that they are well-matched in redshift, both centred at $z=1$. Combined with very similar detection fractions and magnitude distributions in the long-wavelength WISE bands, this redshift match strongly supports the fact that IR selection identifies obscured and unobscured quasars of similar bolometric luminosity. Finally, we perform cross-correlations with confirmed spectroscopic quasars, again confirming the results from other methods --- obscured quasars reside in haloes a factor of 3 times more massive than unobscured quasars. This difference is significant at the $\sim$5$\sigma$ level when the measurements are combined, strong support for the idea that obscuration in at least some quasars is tied to the larger environment, and may have an evolutionary component.

Galaxy-scale bars are expected to provide an effective means for driving material towards the central region in spiral galaxies, and possibly feeding supermassive black holes (BHs). Here we present a statistically-complete study of the effect of bars on average BH accretion. From a well-selected sample of 50,794 spiral galaxies (with M* ~ 0.2-30 x 10^10 Msun) extracted from the Sloan Digital Sky Survey Galaxy Zoo 2 project, we separate those sources considered to contain galaxy-scale bars from those that do not. Using archival data taken by the Chandra X-ray Observatory, we identify X-ray luminous (L_X >~ 10^41 erg/s) active galactic nuclei (AGN) and perform an X-ray stacking analysis on the remaining X-ray undetected sources. Through X-ray stacking, we derive a time-averaged look at accretion for galaxies at fixed stellar mass and star formation rate, finding that the average nuclear accretion rates of galaxies with bar structures are fully consistent with those lacking bars (Mdot_acc ~ 3 x 10^-5 Msun/yr). Hence, we robustly conclude that large-scale bars have little or no effect on the average growth of BHs in nearby (z < 0.15) galaxies over gigayear timescales.

We present a sample of 10 low-mass active galactic nuclei (AGNs) selected from the 40-month NuSTAR serendipitous survey. The sample is selected to have robust NuSTAR detections at $3 - 24$~keV, to be at $z < 0.3$, and to have optical r-band magnitudes at least 0.5~mag fainter than an $L_\star$ galaxy at its redshift. The median values of absolute magnitude, stellar mass and 2--10 X-ray luminosity of our sample are $\langle M_r\rangle = -20.03$, $\langle M_\star\rangle = 4.6\times10^{9}M_\odot$, and $\langle L_{2-10\mathrm{keV}}\rangle = 3.1\times10^{42}$ erg s$^{-1}$, respectively. Five objects have detectable broad H$\alpha$ emission in their optical spectra, indicating black-hole masses of $(1.1-10.4)\times 10^6 M_\odot$. We find that $30^{+17}_{-10}\%$ of the galaxies in our sample do not show AGN-like optical narrow emission lines, and one of the ten galaxies in our sample, J115851+4243.2, shows evidence for heavy X-ray absorption. This result implies that a non-negligible fraction of low-mass galaxies might harbor accreting massive black holes that are missed by optical spectroscopic surveys and $<10$ keV X-ray surveys. The mid-IR colors of our sample also indicate these optically normal low-mass AGNs cannot be efficiently identified with typical AGN selection criteria based on WISE colors. While the hard ($>10$ keV) X-ray selected low-mass AGN sample size is still limited, our results show that sensitive NuSTAR observations are capable of probing faint hard X-ray emission originating from the nuclei of low-mass galaxies out to moderate redshift ($z<0.3$), thus providing a critical step in understanding AGN demographics in low-mass galaxies.

We present the first full catalog and science results for the NuSTAR serendipitous survey. The catalog incorporates data taken during the first 40 months of NuSTAR operation, which provide ~20Ms of effective exposure time over 331 fields, with an areal coverage of 13 sq deg, and 497 sources detected in total over the 3-24 keV energy range. There are 276 sources with spectroscopic redshifts and classifications, largely resulting from our extensive campaign of ground-based spectroscopic followup. We characterize the overall sample in terms of the X-ray, optical, and infrared source properties. The sample is primarily comprised of active galactic nuclei (AGNs), detected over a large range in redshift from z = 0.002 - 3.4 (median of <z> = 0.56), but also includes 16 spectroscopically confirmed Galactic sources. There is a large range in X-ray flux, from log( f_3-24keV / erg s^-1 cm^-2 ) ~ -14 to -11, and in rest-frame 10-40 keV luminosity, from log( L_10-40keV / erg s^-1 ) ~ 39 to 46, with a median of 44.1. Approximately 79% of the NuSTAR sources have lower energy (<10 keV) X-ray counterparts from XMM-Newton, Chandra, and Swift/XRT. The mid-infrared (MIR) analysis, using WISE all-sky survey data, shows that MIR AGN color selections miss a large fraction of the NuSTAR-selected AGN population, from ~15% at the highest luminosities (Lx > 10^44 erg s^-1) to ~80% at the lowest luminosities (Lx < 10^43 erg s^-1). Our optical spectroscopic analysis finds that the observed fraction of optically obscured AGNs (i.e., the Type 2 fraction) is F_Type2 = 53(+14-15)%, for a well-defined subset of the 8-24 keV selected sample. This is higher, albeit at a low significance level, than the Type 2 fraction measured for redshift- and luminosity-matched AGNs selected by <10 keV X-ray missions.

We analyze the relationship between star formation (SF), substructure, and supercluster environment in a sample of 107 nearby galaxy clusters using data from the Sloan Digital Sky Survey. Previous works have investigated the relationships between SF and cluster substructure, and cluster substructure and supercluster environment, but definitive conclusions relating all three of these variables has remained elusive. We find an inverse relationship between cluster SF fraction (f_SF) and supercluster environment density, calculated using the galaxy luminosity density field at a smoothing length of 8 h^-1 Mpc (D8). The slope of f_SF vs. D8 is -0.008 +/- 0.002. The f_SF of clusters located in low-density large-scale environments, 0.244 +/- 0.011, is higher than for clusters located in high-density supercluster cores, 0.202 +/- 0.014. We also divide superclusters, according to their morphology, into filament- and spider-type systems. The inverse relationship between cluster f_SF and large-scale density is dominated by filament- rather than spider-type superclusters. In high-density cores of superclusters, we find a higher f_SF in spider-type superclusters, 0.229 +/- 0.016, than in filament-type superclusters, 0.166 +/- 0.019. Using principal component analysis, we confirm these results and the direct correlation between cluster substructure and SF. These results indicate that cluster SF is affected by both the dynamical age of the cluster (younger systems exhibit higher amounts of SF); the large-scale density of the supercluster environment (high-density core regions exhibit lower amounts of SF); and supercluster morphology (spider-type superclusters exhibit higher amounts of SF at high densities).

Our understanding of the cosmic evolution of supermassive black holes (SMBHs) has been revolutionized by the advent of large multiwavelength extragalactic surveys, which have enabled detailed statistical studies of the host galaxies and large-scale structures of active galactic nuclei (AGN). We give an overview of some recent results on SMBH evolution, including the connection between AGN activity and star formation in galaxies, the role of galaxy mergers in fueling AGN activity, the nature of luminous obscured AGN, and the connection between AGN and their host dark matter halos. We conclude by looking to the future of large-scale extragalactic X-ray and spectroscopic surveys.

Recent measurements of the dark matter halo masses of infrared-selected obscured quasars are in tension --- some indicate that obscured quasars have higher halo mass compared to their unobscured counterparts, while others find no difference. The former result is inconsistent with the simplest models of quasar unification that rely solely on viewing angle, while the latter may support such models. Here, using empirical relationships between dark matter halo and supermassive black hole masses, we provide a simple evolutionary picture that naturally explains these findings and is motivated by more sophisticated merger-driven quasar fueling models. The model tracks the growth rate of haloes, with the black hole growing in spurts of quasar activity in order to "catch-up" with the $M_{\textrm{BH}}$ - $M_{\textrm{star}}$ - $M_{\textrm{halo}}$ relationship. The first part of the quasar phase is obscured and is followed by an unobscured phase. Depending on the luminosity limit of the sample, driven by observational selection effects, a difference in halo masses may or may not be significant. For high luminosity samples, the difference can be large (a few to 10 times higher masses in obscured quasars), while for lower luminosity samples the halo mass difference is very small, much smaller than current observational constraints. Such a simple model provides a qualitative explanation for the higher mass haloes of obscured quasars, as well as rough quantitative agreement with seemingly disparate results.

We present follow-up X-ray observations of the candidate massive black hole (BH) in the nucleus of the low-mass, compact starburst galaxy Henize 2-10. Using new high-resolution observations from the Chandra X-ray Observatory totaling 200 ks in duration, as well as archival Chandra observations from 2001, we demonstrate the presence of a previously unidentified X-ray point source that is spatially coincident with the known nuclear radio source in Henize 2-10 (i.e., the massive BH). We show that the hard X-ray emission previously identified in the 2001 observation is dominated by a source that is distinct from the nucleus, with the properties expected for a high-mass X-ray binary. The X-ray luminosity of the nuclear source suggests the massive BH is radiating significantly below its Eddington limit (~10^-6 L_Edd), and the soft spectrum resembles other weakly accreting massive BHs including Sagittarius A*. Analysis of the X-ray light curve of the nucleus reveals the tentative detection of a ~9-hour periodicity, although additional observations are required to confirm this result. Our study highlights the need for sensitive high-resolution X-ray observations to probe low-level accretion, which is the dominant mode of BH activity throughout the Universe.

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.

We use time-domain optical spectroscopy to distinguish between broad emission lines powered by accreting black holes (BHs) or stellar processes (i.e., supernovae) for 16 galaxies identified as AGN candidates by Reines \etal (2013). Our study is primarily focused on those objects with narrow emission-line ratios dominated by star formation. Based on follow-up spectra taken with the Magellan Echellette Spectrograph (MagE), the Dual Imaging Spectrograph, and the Ohio State Multi-Object Spectrograph, we find that the broad H$\alpha$ emission has faded or was ambiguous for all of the star-forming objects (14/16) over baselines ranging from 5 to 14 years. For the two objects in our follow-up sample with narrow-line AGN signatures (RGG 9 and RGG 119), we find persistent broad H$\alpha$ emission consistent with an AGN origin. Additionally, we use our MagE observations to measure stellar velocity dispersions for 15 objects in the Reines et al. (2013) sample, all with narrow-line ratios indicating the presence of an AGN. Stellar masses range from $\sim5\times10^{8}$ to $3\times10^{9}$~\msun, and we measure $\sigma_{\ast}$ ranging from $28-71~{\rm km~s^{-1}}$. These $\sigma_{\ast}$ correspond to some of the lowest-mass galaxies with optical signatures of AGN activity. We show that RGG 119, the one object which has both a measured $\sigma_{\ast}$ and persistent broad H$\alpha$ emission, falls near the extrapolation of the $\rm M_{BH}-\sigma_{\star}$ relation to the low-mass end.

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.

There is mounting evidence that supermassive black holes form and grow in bulgeless galaxies. However, a robust determination of the fraction of AGNs in bulgeless galaxies, an important constraint to models of supermassive black hole seed formation and merger-free models of AGN fueling, is unknown, since optical studies have been shown to be incomplete for low mass AGNs. In a recent study using the Wide-field Infrared Survey Explorer, we discovered hundreds of bulgeless galaxies that display mid-infrared signatures of extremely hot dust suggestive of powerful accreting massive black holes, despite having no signatures of black hole activity at optical wavelengths. Here we report X-ray follow-up observations of J122434.66+555522.3, a nearby (z=0.052) isolated bulgeless galaxy that contains an unresolved X-ray source detected at the 3 sigma level by XMM-Newton with an observed luminosity uncorrected for intrinsic absorption of L{2-10~keV}=1.1+/-0.4 10^40 ergs/s. Ground-based near-infrared spectroscopy with the Large Binocular Telescope together with multiwavelength observations from ultraviolet to millimeter wavelengths together suggest that J1224+5555 harbors a highly absorbed AGN with an intrinsic absorption of ~NH >10^24 cm^-2. The hard X-ray luminosity of the putative AGN corrected for absorption is L{2-10~keV}~3x10^42 ergs/s, which, depending on the bolometric correction factor, corresponds to a bolometric luminosity of the AGN of 6x10^43ergs/s - 3x10^44 erg/s, and a lower mass limit for the black hole of MBH~2x10^6 Msun, based on the Eddington limit.

An important question in extragalactic astronomy concerns the distribution of black hole accretion rates of active galactic nuclei (AGN). Based on observations at X-ray wavelengths, the observed Eddington ratio distribution appears as a power law, while optical studies have often yielded a lognormal distribution. There is increasing evidence that these observed discrepancies may be due to contamination by star formation and other selection effects. Using a sample of galaxies from the Sloan Digital Sky Survey Data Release 7, we test if an intrinsic Eddington ratio distribution that takes the form of a Schechter function is consistent with previous work that suggests that young galaxies in optical surveys have an observed lognormal Eddington ratio distribution. We simulate the optical emission line properties of a population of galaxies and AGN using a broad instantaneous luminosity distribution described by a Schechter function near the Eddington limit. This simulated AGN population is then compared to observed galaxies via the positions on an emission line excitation diagram and Eddington ratio distributions. We present an improved method for extracting the AGN distribution using BPT diagnostics that allows us to probe over one order of magnitude lower in Eddington ratio counteracting the effects of dilution by star formation. We conclude that for optically selected AGN in young galaxies, the intrinsic Eddington ratio distribution is consistent with a possibly universal, broad power law with an exponential cutoff, as this distribution is observed in old optically selected galaxies and in X-rays.

Recent studies have found that obscured quasars cluster more strongly and are thus hosted by dark matter haloes of larger mass than their unobscured counterparts. These results pose a challenge for the simplest unification models, in which obscured objects are intrinsically the same as unobscured sources but seen through a dusty line of sight. There is general consensus that a structure like a "dusty torus" exists, meaning that this intrinsic similarity is likely the case for at least some subset of obscured quasars. However, the larger host halo masses of obscured quasars implies that there is a second obscured population that has an even higher clustering amplitude and typical halo mass. Here, we use simple assumptions about the host halo mass distributions of quasars, along with analytical methods and cosmological $N$-body simulations to isolate the signal from this population. We provide values for the bias and halo mass as a function of the fraction of the "non-torus obscured" population. Adopting a reasonable value for this fraction of $\sim$25% implies a non-torus obscured quasar bias that is much higher than the observed obscured quasar bias, because a large fraction of the obscured population shares the same clustering strength as the unobscured objects. For this non-torus obscured population, we derive a bias of $\sim$3, and typical halo masses of $\sim3\times10^{13}$ M$_{\odot}/h$ at $z=1$. These massive haloes are likely the descendants of high-mass unobscured quasars at high redshift, and will evolve into members of galaxy groups at $z=0$.

We explore the gas ionization and kinematics, as well as the optical--IR spectral energy distributions for UGC 11185, a nearby pair of merging galaxies hosting obscured active galactic nuclei (AGNs), also known as SDSS J181611.72+423941.6 and J181609.37+423923.0 (J1816NE and J1816SW, $z \approx 0.04$). Due to the wide separation between these interacting galaxies ($\sim 23$ kpc), observations of these objects provide a rare glimpse of the concurrent growth of supermassive black holes at an early merger stage. We use BPT line diagnostics to show that the full extent of the narrow line emission in both galaxies is photoionized by an AGN and confirm the existence of a 10-kpc-scale ionization cone in J1816NE, while in J1816SW the AGN narrow-line region is much more compact (1--2 kpc) and relatively undisturbed. Our observations also reveal the presence of ionized gas that nearly spans the entire distance between the galaxies which is likely in a merger-induced tidal stream. In addition, we carry out a spectral analysis of the X-ray emission using data from {\em XMM-Newton}. These galaxies represent a useful pair to explore how the [\ion{O}{3}] luminosity of an AGN is dependent on the size of the region used to explore the extended emission. Given the growing evidence for AGN "flickering" over short timescales, we speculate that the appearances and impact of these AGNs may change multiple times over the course of the galaxy merger, which is especially important given that these objects are likely the progenitors of the types of systems commonly classified as "dual AGNs."

Some reddened quasars appear to be transitional objects in the merger-induced black hole growth/galaxy evolution paradigm, where a heavily obscured nucleus starts to be unveiled by powerful quasar winds evacuating the surrounding cocoon of dust and gas. Hard X-ray observations are able to peer through this gas and dust, revealing the properties of circumnuclear obscuration. Here, we present NuSTAR and XMM-Newton/Chandra observations of FIRST-2MASS selected red quasars F2M 0830+3759 and F2M 1227+3214. We find that though F2M 0830+3759 is moderately obscured ($N_{\rm H,Z} = 2.1\pm0.2 \times10^{22}$ cm$^{-2}$) and F2M 1227+3214 is mildly absorbed ($N_{\rm H,Z} = 3.4^{+0.8}_{-0.7}\times10^{21}$ cm$^{-2}$) along the line-of-sight, heavier global obscuration may be present in both sources, with $N_{\rm H,S} = 3.7^{+4.1}_{-2.6} \times 10^{23}$ cm$^{-2}$ and $< 5.5\times10^{23}$ cm$^{-2}$, for F2M 0830+3759 and F2M 1227+3214, respectively. F2M 0830+3759 also has an excess of soft X-ray emission below 1 keV which is well accommodated by a model where 7% of the intrinsic AGN X-ray emission is scattered into the line-of-sight. While F2M 1227+3214 has a dust-to-gas ratio ($E(B-V)$/$N_{\rm H}$) consistent with the Galactic value, the $E(B-V)$/$N_{\rm H}$ value for F2M 0830+3759 is lower than the Galactic standard, consistent with the paradigm that the dust resides on galactic scales while the X-ray reprocessing gas originates within the dust-sublimation zone of the broad-line-region. The X-ray and 6.1$\mu$m luminosities of these red quasars are consistent with the empirical relations derived for high-luminosity, unobscured quasars, extending the parameter space of obscured AGN previously observed by NuSTAR to higher luminosities.

Hot Dust-Obscured Galaxies (Hot DOGs) are a population of hyper-luminous infrared galaxies identified by the WISE mission from their very red mid-IR colors, and characterized by hot dust temperatures ($T>60~\rm K$). Several studies have shown clear evidence that the IR emission in these objects is powered by a highly dust-obscured AGN that shows close to Compton-thick absorption at X-ray wavelengths. Thanks to the high AGN obscuration, the host galaxy is easily observable, and has UV/optical colors usually consistent with those of a normal galaxy. Here we discuss a sub-population of 8 Hot DOGs that show enhanced rest-frame UV/optical emission. We discuss three scenarios that might explain the excess UV emission: (i) unobscured light leaked from the AGN by reflection over the dust or by partial coverage of the accretion disk; (ii) a second unobscured AGN in the system; or (iii) a luminous young starburst. X-ray observations can help discriminate between these scenarios. We study in detail the blue excess Hot DOG WISE J020446.13-050640.8, which was serendipitously observed by Chandra/ACIS-I for 174.5 ks. The X-ray spectrum is consistent with a single, hyper-luminous, highly absorbed AGN, and is strongly inconsistent with the presence of a secondary unobscured AGN. Based on this, we argue that the excess blue emission in this object is most likely either due to reflection or a co-eval starburst. We favor the reflection scenario as the unobscured star-formation rate needed to power the UV/optical emission would be $\gtrsim 1000~\rm M_{\odot}~\rm yr^{-1}$. Deep polarimetry observations could confirm the reflection hypothesis.

It is thought that a class of broad absorption line (BAL) QSOs, characterised by Fe absorption features in their UV spectra (called `FeLoBALs'), could mark a transition stage between the end of an obscured starburst event and a youthful QSO beginning to shed its dust cocoon, where Fe has been injected into the interstellar medium by the starburst. To test this hypothesis we have undertaken deep SCUBA-2 850 $\mu$m observations of a sample of 17 FeLoBAL QSOs with 0.89 $\leq$ z $\leq$ 2.78 and -23.31 $\leq$ M$_{B}$ $\leq$-28.50 to directly detect an excess in the thermal emission of the dust which would probe enhanced star-formation activity. We find that FeLoBALs are not luminous sources in the submillimetre, none of them are individually detected at 850 $\mu$m, nor as a population through stacking ($F_{s}=1.14\pm0.58$ mJy). Statistical and survival analyses reveal that FeLoBALs have sub-mm properties consistent with BAL and non-BAL QSOs with matched redshifts and magnitudes. An SED fitting analysis shows that the FIR emission is dominated by AGN activity, and a starburst component is required only in 6/17 sources of our sample; moreover the integrated total luminosity of 16/17 sources is L$\geq$10$^{12}$L$_{\odot}$, high enough to classify FeLoBALs as infrared luminous. In conclusion, we do not find any evidence in support of FeLoBAL QSOs being a transition population between a ULIRG and an unobscured QSO; in particular, FeLoBALs are not characterized by a cold starburst which would support this hypothesis.

Radio emission from radio-quiet quasars may be due to star formation in the quasar host galaxy, to a jet launched by the supermassive black hole, or to relativistic particles accelerated in a wide-angle radiatively-driven outflow. In this paper we examine whether radio emission from radio-quiet quasars is a byproduct of star formation in their hosts. To this end we use infrared spectroscopy and photometry from Spitzer and Herschel to estimate or place upper limits on star formation rates in hosts of ~300 obscured and unobscured quasars at z<1. We find that low-ionization forbidden emission lines such as [NeII] and [NeIII] are likely dominated by quasar ionization and do not provide reliable star formation diagnostics in quasar hosts, while PAH emission features may be suppressed due to the destruction of PAH molecules by the quasar radiation field. While the bolometric luminosities of our sources are dominated by the quasars, the 160 micron fluxes are likely dominated by star formation, but they too should be used with caution. We estimate median star formation rates to be 6-29 Msun/year, with obscured quasars at the high end of this range. This star formation rate is insufficient to explain the observed radio emission from quasars by an order of magnitude, with log(L_radio, observed/L_radio, SF)=0.6-1.3 depending on quasar type and star formation estimator. Although radio-quiet quasars in our sample lie close to the 8-1000 micron infrared / radio correlation characteristic of the star-forming galaxies, both their infrared emission and their radio emission are dominated by the quasar activity, not by the host galaxy.

Using the most recent releases of WISE and Planck data, we perform updated measurements of the bias and typical dark matter halo mass of infrared-selected obscured and unobscured quasars, using the angular autocorrelation function and cosmic microwave background (CMB) lensing cross-correlations. Since our recent work of this kind, the WISE Allwise catalogue was released with improved photometry, and the Planck mission was completed and released improved products. These new data provide a more reliable measurement of the quasar bias and provide an opportunity to explore the role of changing survey pipelines in results downstream. We present a comparison of IR color-selected quasars, split into obscured and unobscured populations based on optical-IR colors, selected from two versions of the WISE data. Which combination of data is used impacts the final results, particularly for obscured quasars, both because of mitigation of some systematics and because the newer catalogue provides a slightly different sample. We show that Allwise data is superior in several ways, though there may be some systematic trends with Moon contamination that were not present in the previous catalogue. We opt currently for the most conservative sample that meet our selection criteria in both the previous and new WISE catalogues. We measure a higher bias and halo mass for obscured quasars ($b_{\textrm{obsc}} \sim 2.1$, $b_{\textrm{unob}} \sim 1.8$) --- at odds with simple orientation models --- but at a reduced significance ($\sim$1.5$\sigma$) as compared to our work with previous survey data.

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.

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.

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.

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.

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.

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.

We present results on gas flows in the halo of a Milky Way-like galaxy at z=0.413 based on high-resolution spectroscopy of a background galaxy. This is the first study of circumgalactic gas at high spectral resolution towards an extended background source (i.e., a galaxy rather than a quasar). Using longslit spectroscopy of the foreground galaxy, we observe spatially extended H alpha emission with circular rotation velocity v=270 km/s. Using echelle spectroscopy of the background galaxy, we detect Mg II and Fe II absorption lines at impact parameter rho=27 kpc that are blueshifted from systemic in the sense of the foreground galaxy's rotation. The strongest absorber EW(2796) = 0.90 A has an estimated column density (N_H>10^19 cm-2) and line-of-sight velocity dispersion (sigma=17 km/s) that are consistent with the observed properties of extended H I disks in the local universe. Our analysis of the rotation curve also suggests that this r=30 kpc gaseous disk is warped with respect to the stellar disk. In addition, we detect two weak Mg II absorbers in the halo with small velocity dispersions (sigma<10 km/s). While the exact geometry is unclear, one component is consistent with an extraplanar gas cloud near the disk-halo interface that is co-rotating with the disk, and the other is consistent with a tidal feature similar to the Magellanic Stream. We can place lower limits on the cloud sizes (l>0.4 kpc) for these absorbers given the extended nature of the background source. We discuss the implications of these results for models of the geometry and kinematics of gas in the circumgalactic medium.

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.

We investigate whether the fuelling of low excitation radio galaxies (LERGs) is linked to major galaxy interactions. Our study utilizes a sample of 10,800 spectroscopic galaxy pairs and 97 post-mergers selected from the Sloan Digital Sky Survey with matches to multi-wavelength datasets. The LERG fraction amongst interacting galaxies is a factor of 3.5 higher than that of a control sample matched in local galaxy density, redshift and stellar mass. However, the LERG excess in pairs does not depend on projected separation and remains elevated out to at least 500 kpc, suggesting that major mergers are not their main fuelling channel. In order to identify the primary fuelling mechanism of LERGs, we compile samples of control galaxies that are matched in various host galaxy and environmental properties. The LERG excess is reduced, but not completely removed, when halo mass or D4000 are included in the matching parameters. However, when BOTH M_halo and D4000 are matched, there is no LERG excess and the 1.4 GHz luminosities (which trace jet mechanical power) are consistent between the pairs and control. In contrast, the excess of optical and mid-IR selected AGN in galaxy pairs is unchanged when the additional matching parameters are implemented. Our results suggest that whilst major interactions may trigger optically and mid-IR selected AGN, the gas which fuels the LERGs has two secular origins: one associated with the large scale environment, such as accretion from the surrounding medium or minor mergers, plus an internal stellar mechanism, such as winds from evolved stars.

We investigate the relationship between star formation (SF) and level of relaxation in a sample of 379 galaxy clusters at z < 0.2. We use data from the Sloan Digital Sky Survey to measure cluster membership and level of relaxation, and to select star-forming galaxies based on mid-infrared emission detected with the Wide-Field Infrared Survey Explorer. For galaxies with absolute magnitudes M_r < -19.5, we find an inverse correlation between SF fraction and cluster relaxation: as a cluster becomes less relaxed, its SF fraction increases. Furthermore, in general, the subtracted SF fraction in all unrelaxed clusters (0.117 +/- 0.003) is higher than that in all relaxed clusters (0.097 +/- 0.005). We verify the validity of our SF calculation methods and membership criteria through analysis of previous work. Our results agree with previous findings that a weak correlation exists between cluster SF and dynamical state, possibly because unrelaxed clusters are less evolved relative to relaxed clusters.

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

Between July 2012 and February 2013, NuSTAR and XMM-Newton performed four long-look joint obser- vations of the type 1.8 Seyfert, NGC 1365. We have analyzed the variable absorption seen in these observations in order to characterize the geometry of the absorbing material. Two of the observations caught NGC 1365 in an unusually low absorption state, revealing complexity in the multi-layer absorber which had previously been hidden. We find the need for three distinct zones of neutral absorption in addition to the two zones of ionized absorption and the Compton-thick torus previously seen in this source. The most prominent absorber is likely associated with broad line region clouds with column densities of around $\sim\,$10$^{23}$ cm$^{-2}$ and a highly clumpy nature as evidenced by an occultation event in February 2013. We also find evidence of a patchy absorber with a variable column around $\sim\,10^{22}$ cm$^{-2}$ and a line of sight covering fraction of 0.3$-$0.9 which responds directly to the intrinsic source flux, possibly due to a wind geometry. A full-covering, constant absorber with a low column density of $\sim\,1 \times$ 10$^{22}$ cm$^{-2}$ is also present, though the location of this low density haze is unknown.

The covering factor of Compton-thick obscuring material associated with the torus in active galactic nuclei (AGN) is at present best understood through the fraction of sources exhibiting Compton-thick absorption along the line of sight ($N_{H}>1.5\times10^{24}$ cm$^{-2}$) in the X-ray band, which reveals the average covering factor. Determining this Compton-thick fraction is difficult however, due to the extreme obscuration. With its spectral coverage at hard X-rays ($>$10 keV), NuSTAR is sensitive to the AGN covering factor since Compton scattering of X-rays off optically thick material dominates at these energies. We present a spectral analysis of 10 AGN observed with NuSTAR where the obscuring medium is optically thick to Compton scattering, so called Compton-thick (CT) AGN. We use the torus models of Brightman & Nandra which predict the X-ray spectrum from reprocessing in a torus and include the torus opening angle as a free parameter and aim to determine the covering factor of the Compton-thick gas in these sources individually. Across the sample we find mild to heavy Compton-thick columns, with $N_{H}$ measured from $10^{24}-10^{26}$ cm$^{-2}$, and a wide range of covering factors, where individual measurements range from 0.2-0.9. We find that the covering factor, $f_{c}$, is a strongly decreasing function of the intrinsic 2-10 keV luminosity, $L_{X}$, where $f_{c}=(-0.41\pm0.13)$log$_{10}$($L_{X}$/erg s$^{-1}$)$+18.31\pm5.33$, across more than two orders of magnitude in $L_{X}$ (10$^{41.5}-10^{44}$ erg s$^{-1}$). The covering factors measured here agree well with the obscured fraction as a function of $L_{X}$ as determined by studies of local AGN with $L_{X}>10^{42.5}$ erg s$^{-1}$.

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

Recent observations have revealed that starburst galaxies can drive molecular gas outflows through stellar radiation pressure. Molecular gas is the phase of the interstellar medium from which stars form, so these outflows curtail stellar mass growth in galaxies. Previously known outflows, however, involve small fractions of the total molecular gas content and are restricted to sub-kiloparsec scales. It is also apparent that input from active galactic nuclei is in at least some cases dynamically important, so pure stellar feedback has been considered incapable of aggressively terminating star formation on galactic scales. Extraplanar molecular gas has been detected in the archetype starburst galaxy M82, but so far there has been no evidence that starbursts can propel significant quantities of cold molecular gas to the same galactocentric radius (~10 kpc) as the warmer gas traced by metal absorbers. Here we report observations of molecular gas in a compact (effective radius 100 pc) massive starburst galaxy at z=0.7, which is known to drive a fast outflow of ionized gas. We find that 35 per cent of the total molecular gas is spatially extended on a scale of approximately 10 kpc, and one third of this has a velocity of up to 1000 km/s. The kinetic energy associated with this high-velocity component is consistent with the momentum flux available from stellar radiation pressure. This result demonstrates that nuclear bursts of star formation are capable of ejecting large amounts of cold gas from the central regions of galaxies, thereby strongly affecting their evolution.

We cross-correlate a cosmic microwave background (CMB) lensing map with the projected space densities of quasars to measure the bias and halo masses of a quasar sample split into obscured and unobscured populations, the first application of this method to distinct quasar subclasses. Several recent studies of the angular clustering of obscured quasars have shown that these objects likely reside in higher-mass halos compared to their unobscured counterparts. This has important implications for models of the structure and geometry of quasars, their role in growing supermassive black holes, and mutual quasar/host galaxy evolution. However, the magnitude and significance of this difference has varied from study to study. Using data from \planck, \wise, and SDSS, we follow up on these results using the independent method of CMB lensing cross-correlations. The region and sample are identical to that used for recent angular clustering measurements, allowing for a direct comparison of the CMB-lensing and angular clustering methods. At $z \sim 1$, we find that the bias of obscured quasars is $b_q = 2.57 \pm 0.24$, while that of unobscured quasars is $b_q = 1.89 \pm 0.19$. This corresponds to halo masses of $\log (M_h / M_{\odot} h^{-1}) = 13.24_{-0.15}^{+0.14}$ (obscured) and $\log (M_h / M_{\odot} h^{-1}) = 12.71_{-0.13}^{+0.15}$ (unobscured). These results agree well with with those from angular clustering (well within $1\sigma$), and confirm that obscured quasars reside in host halos $\sim$3 times as massive as halos hosting unobscured quasars. This implies that quasars spend a significant portion of their lifetime in an obscured state, possibly more than one half of the entire active phase.

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

We present the results of an optical spectroscopic survey of a sample of 40 candidate obscured quasars identified on the basis of their mid-infrared emission detected by the Wide-Field Infrared Survey Explorer (WISE). Optical spectra for this survey were obtained using the Robert Stobie Spectrograph (RSS) on the Southern African Large Telescope (SALT). Our sample was selected with WISE colors characteristic of AGNs, as well as red optical to mid-IR colors indicating that the optical/UV AGN continuum is obscured by dust. We obtain secure redshifts for the majority of the objects that comprise our sample (35/40), and find that sources that are bright in the WISE W4 (22$\mu$m) band are typically at moderate redshift (<z> = 0.35$) while sources fainter in W4 are at higher redshifts (<z> = 0.73$). The majority of the sources have narrow emission lines, with optical colors and emission line ratios of our WISE-selected sources that are consistent with the locus of AGN on the rest-frame $g-z$ color vs. [NeIII]$\lambda$3869 / [OII]$\lambda\lambda$3726+3729 line ratio diagnostic diagram. We also use empirical AGN and galaxy templates to model the spectral energy distributions (SEDs) for the objects in our sample, and find that while there is significant variation in the observed SEDs for these objects, the majority require a strong AGN component. Finally, we use the results from our analysis of the optical spectra and the SEDs to compare our selection criteria to alternate criteria presented in the literature. These results verify the efficacy of selecting luminous obscured AGNs based on their WISE colors.

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.

We explore the connection between absorption by neutral gas and extinction by dust in mid-infrared (IR) selected luminous quasars. We use a sample of 33 quasars at redshifts 0.7 < z < 3 in the 9 deg^2 Bo\"otes multiwavelength survey field that are selected using Spitzer Space Telescope Infrared Array Camera colors and are well-detected as luminous X-ray sources (with >150 counts) in Chandra observations. We divide the quasars into dust-obscured and unobscured samples based on their optical to mid-IR color, and measure the neutral hydrogen column density N_H through fitting of the X-ray spectra. We find that all subsets of quasars have consistent power law photon indices equal to 1.9 that are uncorrelated with N_H. We classify the quasars as gas-absorbed or gas-unabsorbed if N_H > 10^22 cm^-2 or N_H < 10^22 cm^-2, respectively. Of 24 dust-unobscured quasars in the sample, only one shows clear evidence for significant intrinsic N_H, while 22 have column densities consistent with N_H < 10^22 cm^-2. In contrast, of the nine dust-obscured quasars, six show evidence for intrinsic gas absorption, and three are consistent with N_H < 10^22 cm^-2. We conclude that dust extinction in IR-selected quasars is strongly correlated with significant gas absorption as determined through X-ray spectral fitting. These results suggest that obscuring gas and dust in quasars are generally co-spatial, and confirm the reliability of simple mid-IR and optical photometric techniques for separating quasars based on obscuration.

Recent studies of luminous infrared-selected active galactic nuclei (AGN) suggest that the reddest, most obscured objects display a higher angular clustering amplitude, and thus reside in higher-mass dark matter halos. This is a direct contradiction to the prediction of the simplest unification-by-orientation models of AGN and quasars. However, clustering measurements depend strongly on the "mask" that removes low-quality data and describes the sky and selection function. We find that applying a robust, conservative mask to WISE-selected quasars yields a weaker but still significant difference in the bias between obscured and unobscured quasars. These findings are consistent with results from previous Spitzer surveys, and removes any scale dependence of the bias. For obscured quasars with $\langle z \rangle = 0.99$ we measure a bias of $b_q = 2.67 \pm 0.16$, corresponding to a halo mass of $\log (M_h / M_{\odot} h^{-1}) = 13.3 \pm 0.1$, while for unobscured sources with $\langle z \rangle = 1.04$ we find $b_q = 2.04 \pm 0.17$ with a halo mass $\log (M_h / M_{\odot} h^{-1} )= 12.8 \pm 0.1$. This improved measurement indicates that WISE-selected obscured quasars reside in halos only a few times more massive than the halos of their unobscured counterparts, a reduction in the factor of $\sim$10 larger halo mass as has been previously reported using WISE-selected samples. Additionally, an abundance matching analysis yields lifetimes for both obscured and unobscured quasar phases on the order of a few 100 Myr ($\sim$ 1\% of the Hubble time) --- however, the obscured phase lasts roughly twice as long, in tension with many model predictions.

We examine the spatial extent of the narrow-line regions (NLRs) of a sample of 30 luminous obscured quasars at $0.4 < z < 0.7$ observed with spatially resolved Gemini-N GMOS long-slit spectroscopy. Using the [OIII]$\lambda5007$ emission feature, we estimate the size of the NLR using a cosmology-independent measurement: the radius where the surface brightness falls to 10$^{-15}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$. We then explore the effects of atmospheric seeing on NLR size measurements and conclude that direct measurements of the NLR size from observed profiles are too large by 0.1 - 0.2 dex on average, as compared to measurements made to best-fit S\'{e}rsic or Voigt profiles convolved with the seeing. These data, which span a full order of magnitude in IR luminosity ($\log{(L_{8 \mu \mathrm{m}} / \mathrm{erg\, s}^{-1})} = 44.4 - 45.4$) also provide strong evidence that there is a flattening of the relationship between NLR size and AGN luminosity at a seeing-corrected size of $\sim 7$ kpc. The objects in this sample have high luminosities which place them in a previously under-explored portion of the size-luminosity relationship. These results support the existence of a maximal size of the narrow-line region around luminous quasars; beyond this size either there is not enough gas, or the gas is over-ionized and does not produce enough [OIII]$\lambda5007$ emission.

We investigate the process of rapid star formation quenching in a sample of 12 massive galaxies at intermediate redshift (z~0.6) that host high-velocity ionized gas outflows (v>1000 km/s). We conclude that these fast outflows are most likely driven by feedback from star formation rather than active galactic nuclei (AGN). We use multiwavelength survey and targeted observations of the galaxies to assess their star formation, AGN activity, and morphology. Common attributes include diffuse tidal features indicative of recent mergers accompanied by bright, unresolved cores with effective radii less than a few hundred parsecs. The galaxies are extraordinarily compact for their stellar mass, even when compared with galaxies at z~2-3. For 9/12 galaxies, we rule out an AGN contribution to the nuclear light and hypothesize that the unresolved core comes from a compact central starburst triggered by the dissipative collapse of very gas-rich progenitor merging disks. We find evidence of AGN activity in half the sample but we argue that it accounts for only a small fraction (<10%) of the total bolometric luminosity. We find no correlation between AGN activity and outflow velocity and we conclude that the fast outflows in our galaxies are not powered by on-going AGN activity, but rather by recent, extremely compact starbursts.

Interactions between galaxies are predicted to cause gas inflows that can potentially trigger nuclear activity. Since the inflowing material can obscure the central regions of interacting galaxies, a potential limitation of previous optical studies is that obscured Active Galactic Nuclei (AGNs) can be missed at various stages along the merger sequence. We present the first large mid-infrared study of AGNs in mergers and galaxy pairs, in order to quantify the incidence of obscured AGNs triggered by interactions. The sample consists of galaxy pairs and post-mergers drawn from the Sloan Digital Sky Survey that are matched to detections by the Wide Field Infrared Sky Explorer (WISE). We find that the fraction of AGN in the pairs, relative to a mass-, redshift- and environment-matched control sample, increases as a function of decreasing projected separation. This enhancement is most dramatic in the post-merger sample, where we find a factor of 10-20 excess in the AGN fraction compared with the control. Although this trend is in qualitative agreement with results based on optical AGN selection, the mid-infrared selected AGN excess increases much more dramatically in the post-mergers than is seen for optical AGN. Our results suggest that energetically dominant optically obscured AGNs become more prevalent in the most advanced mergers, consistent with theoretical predictions.

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.

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.

We classify the spectral energy distributions (SEDs) of 431,038 sources in the 9 sq. deg Bootes field of the NOAO Deep Wide-Field Survey (NDWFS). There are up to 17 bands of data available per source, including ultraviolet (GALEX), optical (NDWFS), near-IR (NEWFIRM), and mid-infrared (IRAC/MIPS) data, as well as spectroscopic redshifts for ~20,000 objects, primarily from the AGN and Galaxy Evolution Survey (AGES). We fit galaxy, AGN, stellar, and brown dwarf templates to the observed SEDs, which yield spectral classes for the Galactic sources and photometric redshifts and galaxy/AGN luminosities for the extragalactic sources. The photometric redshift precision of the galaxy and AGN samples are sigma/(1+z)=0.040 and sigma/(1+z)=0.169, respectively, with the worst 5% outliers excluded. Based on the reduced chi-squared of the SED fit for each SED model, we are able to distinguish between Galactic and extragalactic sources for sources brighter than I=23.5. We compare the SED fits for a galaxy-only model and a galaxy+AGN model. Using known X-ray and spectroscopic AGN samples, we confirm that SED fitting can be successfully used as a method to identify large populations of AGN, including spatially resolved AGN with significant contributions from the host galaxy and objects with the emission line ratios of "composite" spectra. We also use our results to compare to the X-ray, mid-IR, optical color and emission line ratio selection techniques. For an F-ratio threshold of F>10 we find 16,266 AGN candidates brighter than I=23.5 and a surface density of ~1900 AGN per deg^2.

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 ($\alpha_{\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.

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.

We investigate the relationship between star formation (SF) and substructure in a sample of 107 nearby galaxy clusters using data from the Sloan Digital Sky Survey (SDSS). Several past studies of individual galaxy clusters have suggested that cluster mergers enhance cluster SF, while others find no such relationship. The SF fraction in multi-component clusters (0.228 +/- 0.007) is higher than that in single-component clusters (0.175 +/- 0.016) for galaxies with M^0.1_r < -20.5. In both single- and multi-component clusters, the fraction of star-forming galaxies increases with clustercentric distance and decreases with local galaxy number density, and multi-component clusters show a higher SF fraction than single-component clusters at almost all clustercentric distances and local densities. Comparing the SF fraction in individual clusters to several statistical measures of substructure, we find weak, but in most cases significant at greater than 2 sigma, correlations between substructure and SF fraction. These results could indicate that cluster mergers may cause weak but significant SF enhancement in clusters, or unrelaxed clusters exhibit slightly stronger SF due to their less evolved states relative to relaxed clusters.

We present the results of recent studies on the co-evolution of galaxies and the supermassive black holes (SMBHs) using Herschel far-infrared and Chandra X-ray observations in the Bo\"otes survey region. For a sample of star-forming (SF) galaxies, we find a strong correlation between galactic star formation rate and the average SMBH accretion rate in SF galaxies. Recent studies have shown that star formation and AGN accretion are only weakly correlated for individual AGN, but this may be due to the short variability timescale of AGN relative to star formation. Averaging over the full AGN population yields a strong linear correlation between accretion and star formation, consistent with a simple picture in which the growth of SMBHs and their host galaxies are closely linked over galaxy evolution time scales.

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

We present the results of a combined galaxy population analysis for the host galaxies of active galactic nuclei (AGN) identified at 0 < z < 1.4 within the SDSS, Bootes and DEEP2 surveys. We identified AGN in a uniform and unbiased manner at X-ray, infrared and radio wavelengths. Supermassive black holes undergoing radiatively-efficient accretion (detected as X-ray and/or infrared AGN) appear to be hosted in a separate and distinct galaxy population than AGN undergoing powerful mechanically dominated accretion (radio AGN). Consistent with some previous studies, radiatively efficient AGN appear to be preferentially hosted in modest star-forming galaxies, with little dependence on AGN or galaxy luminosity. AGN exhibiting radio-emitting jets due to mechanically-dominated accretion are almost exclusively observed in massive, passive galaxies. Crucially, we now provide strong evidence that the observed host-galaxy trends are independent of redshift. In particular, these different accretion-mode AGN have remained as separate galaxy populations throughout the last 9 Gyr. Furthermore, it appears that galaxies hosting AGN have evolved along the same path as galaxies that are not hosting AGN with little evidence for distinctly separate evolution.

We measure the cross-power spectrum of the projected mass density as traced by the convergence of the cosmic microwave background lensing field from the South Pole Telescope (SPT) and a sample of Type 1 and 2 (unobscured and obscured) quasars at z~1 selected with the Wide-field Infrared Survey Explorer, over 2500 deg^2. The cross-power spectrum is detected at ~7-sigma, and we measure a linear bias b=1.67+/-0.24, consistent with clustering analyses. Using an independent lensing map, derived from Planck observations, to measure the cross-spectrum, we find excellent agreement with the SPT analysis. The bias of the combined sample of Type 1 and 2 quasars determined in this work is similar to that previously determined for Type 1 quasars alone; we conclude that that obscured and unobscured quasars must trace the matter field in a similar way. This result has implications for our understanding of quasar unification and evolution schemes.

We present a detailed spectral analysis of new XMM-Newton data of the source CXOC J100043.1+020637, also known as CID-42, detected in the COSMOS survey at z = 0.359. Previous works suggested that CID-42 is a candidate recoiling supermassive black holes showing also an inverted P-Cygni profile in the X- ray spectra at ~6 keV (rest) with an iron emission line plus a redshifted absorption line (detected at 3sigma in previous XMM-Newton and Chandra observations). Detailed analysis of the absorption line suggested the presence of ionized material inflowing into the black hole at high velocity. In the new long XMM-Newton observation, while the overall spectral shape remains constant, the continuum 2-10 keV flux decreased of ~20% with respect to previous observation and the absorption line is undetected. The upper limit on the intensity of the absorption line is EW<162 keV. Extensive Monte Carlo simulations show that the non detection of the line is solely due to variation in the properties of the inflowing material, in agreement with the transient nature of these features, and that the intensity of the line is lower than the previously measured with a probability of 98.8%. In the scenario of CID-42 as recoiling SMBH, the absorption line can be interpreted as due to inflow of gas with variable density and located in the proximity of the SMBH and recoiling with it. New monitoring observations will be requested to further characterize this line.

We perform halo occupation distribution (HOD) modeling of the projected two-point correlation function (2PCF) of high-redshift (z~1.2) X-ray-bright active galactic nuclei (AGN) in the XMM-COSMOS field measured by Allevato et al. The HOD parameterization is based on low-luminosity AGN in cosmological simulations. At the median redshift of z~1.2, we derive a median mass of (1.02+0.21/-0.23)x10^{13} Msun/h for halos hosting central AGN and an upper limit of ~10% on the AGN satellite fraction. Our modeling results indicate (at the 2.5-sigma level) that X-ray AGN reside in more massive halos compared to more bolometrically luminous, optically-selected quasars at similar redshift. The modeling also yields constraints on the duty cycle of the X-ray AGN, and we find that at z~1.2 the average duration of the X-ray AGN phase is two orders of magnitude longer than that of the quasar phase. Our inferred mean occupation function of X-ray AGN is similar to recent empirical measurements with a group catalog and suggests that AGN halo occupancy increases with increasing halo mass. We project the XMM-COSMOS 2PCF measurements to forecast the required survey parameters needed in future AGN clustering studies to enable higher precision HOD constraints and determinations of key physical parameters like the satellite fraction and duty cycle. We find that N^{2}/A~5x10^{6} deg^{-2} (with N the number of AGN in a survey area of A deg^{2}) is sufficient to constrain the HOD parameters at the 10% level, which is easily achievable by upcoming and proposed X-ray surveys.

We present spatially resolved long-slit spectroscopy from the Southern African Large Telescope (SALT) to examine the spatial extent of the narrow-line regions (NLRs) of a sample of 8 luminous obscured quasars at 0.10 < z < 0.43. Our results are consistent with an observed shallow slope in the relationship between NLR size and L_[OIII], which has been interpreted to indicate that NLR size is limited by the density and ionization state of the NLR gas rather than the availability of ionizing photons. We also explore how the NLR size scales with a more direct measure of instantaneous AGN power using mid-IR photometry from WISE, which probes warm to hot dust near the central black hole and so, unlike [OIII], does not depend on the properties of the NLR. Using our results as well as samples from the literature, we obtain a power-law relationship between NLR size and L_8micron that is significantly steeper than that observed for NLR size and L_[OIII]. We find that the size of the NLR goes approximately as L^(1/2)_8micron, as expected from the simple scenario of constant-density clouds illuminated by a central ionizing source. We further see tentative evidence for a flattening of the relationship between NLR size and L_8micron at the high luminosity end, and propose that we are seeing a limiting NLR size of 10 - 20 kpc, beyond which the availability of gas to ionize becomes too low. We find that L_[OIII] ~ L_8micron^(1.4), consistent with a picture in which the L_[OIII] is dependent on the volume of the NLR. These results indicate that high-luminosity quasars have a strong effect in ionizing the available gas in a galaxy.

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 \Gamma=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.

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

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

We present optical and near-infrared observations of SN 2012au, a slow-evolving supernova (SN) with properties that suggest a link between subsets of energetic and H-poor SNe and superluminous SNe. SN 2012au exhibited conspicuous SN Ib-like He I lines and other absorption features at velocities reaching 2 x 10^4 km/s in its early spectra, and a broad light curve that peaked at M_B = -18.1 mag. Models of these data indicate a large explosion kinetic energy of 10^{52} erg and 56Ni mass ejection of 0.3 Msolar on par with SN 1998bw. SN 2012au's spectra almost one year after explosion show a blend of persistent Fe II P-Cyg absorptions and nebular emissions originating from two distinct velocity regions. These late-time emissions include strong [Fe II], [Ca II], [O I], Mg I], and Na I lines at velocities > 4500 km/s, as well as O I and Mg I lines at noticeably smaller velocities of 2000 km/s. Many of the late-time properties of SN 2012au are similar to the slow-evolving hypernovae SN 1997dq and SN 1997ef, and the superluminous SN 2007bi. Our observations suggest that a single explosion mechanism may unify all of these events that span -21 < M_B < -17 mag. The aspherical and possibly jetted explosion was most likely initiated by the core collapse of a massive progenitor star and created substantial high-density, low-velocity Ni-rich material.

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.

We report observations of the CO(2-1) emission of SDSSJ1506+54, a compact (r_e~135pc) starburst galaxy at z=0.6. SDSSJ1506+54 appears to be forming stars close to the limit allowed by stellar radiation pressure feedback models: the measured L_IR/L'_CO 1500 is one of the highest measured for any galaxy. With its compact optical morphology but extended low surface brightness envelope, post-starburst spectral features, high infrared luminosity (L_IR>10^12.5 L_Sun), low gas fraction (M_H2/M_stars~15%), and short gas depletion time (tens of Myr), we speculate that this is a feedback- limited central starburst episode at the conclusion of a major merger. Taken as such, SDSSJ1504+54 epitomizes the brief closing stage of a classic model of galaxy growth: we are witnessing a key component of spheroid formation during what we term a 'redline' starburst.

The fraction of cluster galaxies that host luminous AGN is an important probe of AGN fueling processes, the cold ISM at the centers of galaxies, and how tightly black holes and galaxies co-evolve. We present a new measurement of the AGN fraction in a sample of 13 clusters of galaxies (M >= 10^{14} Msun) at 1<z<1.5 selected from the Spitzer/IRAC Shallow Cluster Survey, as well as the field fraction in the immediate vicinity of these clusters, and combine these data with measurements from the literature to quantify the relative evolution of cluster and field AGN from the present to z~3. We estimate that the cluster AGN fraction at 1<z<1.5 is f_A = 3.0^{+2.4}_{-1.4}% for AGN with a rest-frame, hard X-ray luminosity greater than L_{X,H} >= 10^{44} erg/s. This fraction is measured relative to all cluster galaxies more luminous than M*_{3.6}(z)+1, where M*_{3.6}(z) is the absolute magnitude of the break in the galaxy luminosity function at the cluster redshift in the IRAC 3.6um bandpass. The cluster AGN fraction is 30 times greater than the 3sigma upper limit on the value for AGN of similar luminosity at z~0.25, as well as more than an order of magnitude greater than the AGN fraction at z~0.75. AGN with L_{X,H} >= 10^{43} erg/s exhibit similarly pronounced evolution with redshift. In contrast with the local universe, where the luminous AGN fraction is higher in the field than in clusters, the X-ray and MIR-selected AGN fractions in the field and clusters are consistent at 1<z<1.5. This is evidence that the cluster AGN population has evolved more rapidly than the field population from z~1.5 to the present. This environment-dependent AGN evolution mimics the more rapid evolution of star-forming galaxies in clusters relative to the field.

Many theoretical models require powerful active galactic nuclei (AGNs) to suppress star formation in distant galaxies and reproduce the observed properties of today's massive galaxies. A recent study based on Herschel-SPIRE submillimetre observations claimed to provide direct support for this picture, reporting a significant decrease in the mean star-formation rates (SFRs) of the most luminous AGNs (Lx>10^44 erg/s) at z=1-3 in the Chandra Deep Field-North (CDF-N). In this letter we extend these results using Herschel-SPIRE 250um data in the COSMOS and CDF-S fields to achieve an order of magnitude improvement in the number of sources at Lx>10^44 erg/s. On the basis of our analysis, we find no strong evidence for suppressed star formation in Lx>10^44 erg/s AGNs at z=1-3. The mean SFRs of the AGNs are constant over the broad X-ray luminosity range of Lx~10^43-10^45 erg/s (with mean SFRs consistent with typical star-forming galaxies at z~2; <SFRs>~100-200 Msol/yr). We suggest that the previous CDF-N results were likely due to low number statistics. We discuss our results in the context of current theoretical models.

We present here a new spectral energy distribution (SED) fitting approach that we adopt to select radio-excess sources amongst distant star-forming galaxies in the GOODS-Herschel (North) field and to reveal the presence of hidden, highly obscured AGN. Through extensive SED analysis of 458 galaxies with radio 1.4 GHz and mid-IR 24 um detections using some of the deepest Chandra X-ray, Spitzer and Herschel infrared, and VLA radio data available to date, we have robustly identified a sample of 51 radio-excess AGN (~1300 deg^-2) out to redshift z~3. These radio-excess AGN have a significantly lower far-IR/radio ratio (q<1.68) than the typical relation observed for star-forming galaxies (q~2.2). We find that ~45% of these radio-excess sources have a dominant AGN component in the mid-IR band, while for the remainders the excess radio emission is the only indicator of AGN activity. The fraction of radio-excess AGN increases with X-ray luminosity reaching ~60% at Lx~10^44-10^45 erg/s, making these sources an important part of the total AGN population. However, almost half (24/51) of these radio-excess AGN are not detected in the deep Chandra X-ray data, suggesting that some of these sources might be heavily obscured. We also find that the specific star formation rates (sSFRs) of the radio-excess AGN are on average lower that those observed for X-ray selected AGN hosts, indicating that our sources are forming stars more slowly than typical AGN hosts, and possibly their star formation is progressively quenching.

We present integral field spectroscopy observations, covering the [O III]4959,5007 emission-line doublet of eight high-redshift (z=1.4-3.4) ultra-luminous infrared galaxies (ULIRGs) that host Active Galactic Nuclei (AGN) activity, including known sub-millimetre luminous galaxies (SMGs). The targets have moderate radio luminosities that are typical of high-redshift ULIRGs (L(1.4GHz)=10^(24)-10^(25)W/Hz) and therefore are not radio-loud AGN. We de-couple kinematic components due to the galaxy dynamics and mergers from those due to outflows. We find evidence in the four most luminous systems (L([O III])>~10^(43)erg/s) for the signatures of large-scale energetic outflows: extremely broad [O III] emission (FWHM ~ 700-1400km/s) across ~4-15kpc, with high velocity offsets from the systemic redshifts (up to ~850km/s). The four less luminous systems have lower quality data displaying weaker evidence for spatially extended outflows. We estimate that these outflows are potentially depositing energy into their host galaxies at considerable rates (~10^(43)-10^(45)erg/s); however, due to the lack of constraints on the density of the outflowing material and the structure of the outflow, these estimates should be taken as illustrative only. Based on the measured maximum velocities (v(max)~400-1400km/s) the outflows observed are likely to unbind some fraction of the gas from their host galaxies, but are unlikely to completely remove gas from the galaxy haloes. By using a combination of energetic arguments and a comparison to ULIRGs without clear evidence for AGN activity, we show that the AGN activity could be the dominant power source for driving all of the observed outflows, although star formation may also play a significant role in some of the sources.

We estimate the relative contributions of the supermassive black hole (SMBH) accretion disk, corona, and obscuring torus to the bolometric luminosity of Seyfert galaxies, using Spizter mid-infrared (MIR) observations of a complete sample of 68 nearby active galactic nuclei from the INTEGRAL all-sky hard X-ray (HX) survey. This is the first HX-selected (above 15 keV) sample of AGNs with complementary high angular resolution, high signal to noise, MIR data. Correcting for the host galaxy contribution, we find a correlation between HX and MIR luminosities: L_MIR L_HX^(0.74+/-0.06). Assuming that the observed MIR emission is radiation from an accretion disk reprocessed in a surrounding dusty torus that subtends a solid angle decreasing with increasing luminosity (as inferred from the declining fraction of obscured AGNs), the intrinsic disk luminosity, L_D, is approximately proportional to the luminosity of the corona in the 2-300 keV energy band, L_C, with the L_D/L_C ratio varying by a factor of 2.1 around a mean value of 1.6. This ratio is a factor of ~2 smaller than for typical quasars producing the cosmic X-ray background (CXB). Therefore, over three orders of magnitude in luminosity, HX radiation carries a large, and roughly comparable, fraction of the bolometric output of AGNs. We estimate the cumulative bolometric luminosity density of local AGNs at ~(1-3)x10^40 erg/s/Mpc^3. Finally, the Compton temperature ranges between kT_c~2 and ~6 keV for nearby AGNs, compared to kT_c~2 keV for typical quasars, confirming that radiative heating of interstellar gas can play an important role in regulating SMBH growth.

We present the discovery of compact, obscured star formation in galaxies at z ~ 0.6 that exhibit >1000 km/s outflows. Using optical morphologies from the Hubble Space Telescope and infrared photometry from the Wide-field Infrared Survey Explorer, we estimate star formation rate (SFR) surface densities that approach Sigma_SFR ~ 3000 Msun/yr/kpc^2, comparable to the Eddington limit from radiation pressure on dust grains. We argue that feedback associated with a compact starburst in the form of radiation pressure from massive stars and ram pressure from supernovae and stellar winds is sufficient to produce the high-velocity outflows we observe, without the need to invoke feedback from an active galactic nucleus.

We present the X-ray point-source catalog produced from the Chandra Advanced CCD Imaging Spectrometer (ACIS-I) observations of the combined \sim3.2 deg2 DEEP2 (XDEEP2) survey fields, which consist of four ~0.7-1.1 deg2 fields. The combined total exposures across all four XDEEP2 fields range from ~10ks-1.1Ms. We detect X-ray point-sources in both the individual ACIS-I observations and the overlapping regions in the merged (stacked) images. We find a total of 2976 unique X-ray sources within the survey area with an expected false-source contamination of ~30 sources (~1%). We present the combined logN-logS distribution of sources detected across the XDEEP2 survey fields and find good agreement with the Extended Chandra Deep Field and Chandra-COSMOS fields to f_{X,0.5-2keV}\sim2x10^{-16} erg/cm^2/s. Given the large survey area of XDEEP2, we additionally place relatively strong constraints on the logN-logS distribution at high fluxes (f_{X,0.5-2keV}\sim3x10^{-14} erg/cm^2/s), and find a small systematic offset (a factor ~1.5) towards lower source numbers in this regime, when compared to smaller area surveys. The number counts observed in XDEEP2 are in close agreement with those predicted by X-ray background synthesis models. Additionally, we present a Bayesian-style method for associating the X-ray sources with optical photometric counterparts in the DEEP2 catalog (complete to R_AB < 25.2) and find that 2126 (~71.4\pm2.8%) of the 2976 X-ray sources presented here have a secure optical counterpart with a <6% contamination fraction. We provide the DEEP2 optical source properties (e.g., magnitude, redshift) as part of the X-ray-optical counterpart catalog.

(abridged) We present a clustering analysis of 370 high-confidence H-alpha emitters (HAEs) at z=2.23. The HAEs are detected in the Hi-Z Emission Line Survey (HiZELS), a large-area blank field 2.121um narrowband survey using the United Kingdom Infrared Telescope (UKIRT) Wide Field Camera (WFCAM). Averaging the two-point correlation function of HAEs in two ~1 degree scale fields (United Kingdom Infrared Deep Sky Survey/Ultra Deep Survey [UDS] and Cosmological Evolution Survey [COSMOS] fields) we find a clustering amplitude equivalent to a correlation length of r_0=3.7+/-0.3 Mpc/h for galaxies with star formation rates of >7 M_sun/yr. The data are well-fitted by the expected correlation function of Cold Dark Matter, scaled by a bias factor: omega_HAE=b^2 omega_DM where b=2.4^{+0.1}_{-0.2}. The corresponding 'characteristic' mass for the halos hosting HAEs is log(M_h/[M_sun/h])=11.7+/-0.1. Comparing to the latest semi-analytic GALFORM predictions for the evolution of HAEs in a LCDM cosmology, we find broad agreement with the observations, with GALFORM predicting a HAE correlation length of ~4 Mpc/h. Motivated by this agreement, we exploit the simulations to construct a parametric model of the halo occupation distribution (HOD) of HAEs, and use this to fit the observed clustering. Our best-fitting HOD can adequately reproduce the observed angular clustering of HAEs, yielding an effective halo mass and bias in agreement with that derived from the scaled omega_DM fit, but with the relatively small sample size the current data provide a poor constraint on the HOD. Our results support the broad picture that 'typical' (~L*) star-forming galaxies have been hosted by dark matter haloes with M_h<10^12 M_sun/h since z~2, but with a broad occupation distribution and clustering that is likely to be a strong function of luminosity.

Using X-ray stacking analyses we estimate the average amounts of supermassive black hole (SMBH) growth taking place in star-forming galaxies (SFGs) at z~1 and z~2 as a function of galaxy stellar mass (M*). We find the average SMBH growth rate follows remarkably similar trends with M* and redshift as the average star-formation rates (SFRs) of their host galaxies (i.e., (i.e., dM_BH/dt ~ M*^(0.86+/-0.39) for the z~1 sample and dM_BH/dt ~ M*^(1.05+/-0.36) for the z~2 sample). It follows that the ratio of SMBH growth rate to SFR is (a) flat with respect to M* (b) not evolving with redshift and (c) close to the ratio required to maintain/establish a SMBH to M* ratio of ~10^{-3} as also inferred from today's M_BH-M_bulge relationship. We interpret this as evidence that SMBHs have, on average, grown in-step with their host galaxies since at least z~2, irrespective of host galaxy mass and AGN triggering mechanism. As such, we suggest that the same secular processes that drive the bulk of star formation are also responsible for the majority of SMBH growth. From this, we speculate that it is the availability of gas reservoirs that regulate both cosmological SMBH growth and star formation.

We explore the origin of mid-infrared (mid-IR) dust extinction in all 20 nearby (z < 0.05) bona-fide Compton-thick (N_H > 1.5 x 10^24 cm^-2) AGN with hard energy (E > 10 keV) X-ray spectral measurements. We accurately measure the silicate absorption features at lambda~9.7um in archival low-resolution (R~57-127) Spitzer Infrared Spectrograph (IRS) spectroscopy, and show that only a minority (~45%) of nearby Compton-thick AGN have strong Si-absorption features (S_9.7 = ln(f_{int}/f_{obs}) > 0.5) which would indicate significant dust attenuation. The majority (~60%) are star-formation dominated (AGN:SB<0.5) at mid-IR wavelengths and lack the spectral signatures of AGN activity at optical wavelengths, most likely because the AGN emission-lines are optically-extinguished. Those Compton-thick AGN hosted in low-inclination angle galaxies exhibit a narrow-range in Si-absorption (S_9.7 ~ 0-0.3), which is consistent with that predicted by clumpy-torus models. However, on the basis of the IR spectra and additional lines of evidence, we conclude that the dominant contribution to the observed mid-IR dust extinction is dust located in the host galaxy (i.e., due to disturbed morphologies; dust-lanes; galaxy inclination angles) and not necessarily a compact obscuring torus surrounding the central engine.

Using a sample of 123 X-ray clusters and groups drawn from the XMM-Cluster Survey first data release, we investigate the interplay between the brightest cluster galaxy (BCG), its black hole, and the intra-cluster/group medium (ICM). It appears that for groups and clusters with a BCG likely to host significant AGN feedback, gas cooling dominates in those with Tx > 2 keV while AGN feedback dominates below. This may be understood through the sub-unity exponent found in the scaling relation we derive between the BCG mass and cluster mass over the halo mass range 10^13 < M500 < 10^15Msol and the lack of correlation between radio luminosity and cluster mass, such that BCG AGN in groups can have relatively more energetic influence on the ICM. The Lx - Tx relation for systems with the most massive BCGs, or those with BCGs co-located with the peak of the ICM emission, is steeper than that for those with the least massive and most offset, which instead follows self-similarity. This is evidence that a combination of central gas cooling and powerful, well fuelled AGN causes the departure of the ICM from pure gravitational heating, with the steepened relation crossing self-similarity at Tx = 2 keV. Importantly, regardless of their black hole mass, BCGs are more likely to host radio-loud AGN if they are in a massive cluster (Tx > 2 keV) and again co-located with an effective fuel supply of dense, cooling gas. This demonstrates that the most massive black holes appear to know more about their host cluster than they do about their host galaxy. The results lead us to propose a physically motivated, empirical definition of 'cluster' and 'group', delineated at 2 keV.

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

X-ray reflection signatures are observed around multiple classes of accreting compact objects. Modelling these features yield important constraints on the physics of accretion disks, motivating the development of X-ray reflection models appropriate for a variety of systems and illumination conditions. Here, constant density ionized X-ray reflection models are presented for a disk irradiated with a very hard power-law X-ray spectrum (\Gamma < 1) and a variable high-energy cutoff. These models are then applied to the Suzaku data of the accreting X-ray pulsar LMC X-4, where very good fits are obtained with a highly ionized reflector responsible for both the broad Fe K line and the soft excess. The ionized reflector shows strong evidence for significant Doppler broadening and is redshifted by ~10^4 km/s. These features indicate that the reflecting material is associated with the complex dynamics occurring at the inner region of the magnetically-truncated accretion disk. Thus, reflection studies of X-ray pulsar spectra may give important insights into the accretion physics at the magnetospheric radius.

We present the galaxy optical luminosity function for the redshift range 0.05<z<0.75 from the AGN and Galaxy Evolution Survey (AGES), a spectroscopic survey of 7.6 sq. deg. in the Bootes field of the NOAO Deep Wide-Field Survey. Our statistical sample is comprised of 12,473 galaxies with known redshifts down to I=20.4 (AB). Our results at low redshift are consistent with those from SDSS; at higher redshift, we find strong evidence for evolution in the luminosity function, including differential evolution between blue and red galaxies. We find that the luminosity density evolves as (1+z)^(0.54+/-0.64) for red galaxies and (1+z)^(1.64+/-0.39) for blue galaxies.

We present a measurement of the spatial clustering of submillimetre galaxies (SMGs) at z = 1-3. Using data from the 870 micron LESS survey, we employ a novel technique to measure the cross-correlation between SMGs and galaxies, accounting for the full probability distributions for photometric redshifts of the galaxies. From the observed projected two-point cross-correlation function we derive the linear bias and characteristic dark matter (DM) halo masses for the SMGs. We detect clustering in the cross-correlation between SMGs and galaxies at the > 4 sigma level. For the SMG autocorrelation we obtain r_0 = 7.7 (+1.8,-2.3) h^-1 Mpc, and derive a corresponding DM halo mass of log(M_halo [h^-1 M_sun]) = 12.8 (+0.3,-0.5). Based on the evolution of DM haloes derived from simulations, we show that that the z = 0 descendants of SMGs are typically massive (~2-3 L*) elliptical galaxies residing in moderate- to high-mass groups (log(M_halo [h^-1 M_sun]) = 13.3 (+0.3,-0.5). From the observed clustering we estimate an SMG lifetime of ~100 Myr, consistent with lifetimes derived from gas consumption times and star-formation timescales, although with considerable uncertainties. The clustering of SMGs at z ~ 2 is consistent with measurements for optically-selected quasi-stellar objects (QSOs), supporting evolutionary scenarios linking starbursts and QSOs. Given that SMGs reside in haloes of characteristic mass ~ 6 x 10^12 h^-1 M_sun, we demonstrate that the redshift distribution of SMGs can be described remarkably well by the combination of two effects: the cosmological growth of structure and the evolution of the molecular gas fraction in galaxies. We conclude that the powerful starbursts in SMGs likely represent a short-lived but universal phase in massive galaxy evolution, associated with the transition between cold gas-rich, star-forming galaxies and passively evolving systems. [Abridged]

The X-ray source CXOXBJ142607.6+353351 (CXOJ1426+35), which was identified in a 172 ks Chandra image in the Bootes field, shows double-peaked rest-frame optical/UV emission lines, separated by 0.69" (5.5 kpc) in the spatial dimension and by 690 km s^-1 in the velocity dimension. The high excitation lines and emission line ratios indicate both systems are ionized by an AGN continuum, and the double-peaked profile resembles that of candidate dual AGN. At a redshift of z=1.175, this source is the highest redshift candidate dual AGN yet identified. However, many sources have similar emission line profiles for which other interpretations are favored. We have analyzed the substantial archival data available in this field, as well as acquired near-infrared (NIR) adaptive optics (AO) imaging and NIR slit spectroscopy. The X-ray spectrum is hard, implying a column density of several 10^23 cm^-2. Though heavily obscured, the source is also one of the brightest in the field, with an absorption-corrected 2-10 keV luminosity of ~10^45 erg s^-1. Outflows driven by an accretion disk may produce the double-peaked lines if the central engine accretes near the Eddington limit. However, we may be seeing the narrow line regions of two AGN following a galactic merger. While the AO image reveals only a single source, a second AGN would easily be obscured by the significant extinction inferred from the X-ray data. Understanding the physical processes producing the complex emission line profiles seen in CXOJ1426+35 and related sources is important for interpreting the growing population of dual AGN candidates.

We present a study of the infrared properties of X-ray selected, moderate luminosity (Lx=10^{42}-10^{44}ergs/s) active galactic nuclei (AGNs) up to z~3, to explore the links between star formation in galaxies and accretion onto their central black holes. We use 100um and 160um fluxes from GOODS-Herschel -the deepest survey yet undertaken by the Herschel telescope- and show that in >94 per cent of cases these fluxes are dominated by the host. We find no evidence of any correlation between the X-ray and infrared luminosities of moderate AGNs at any redshift, suggesting that star-formation is decoupled from nuclear (AGN) activity. The star formation rates of AGN hosts increase strongly with redshift; by a factor of 43 from z<0.1 to z=2-3 for AGNs with the same X-ray luminosities. This increase is consistent with the factor of 25-50 increase in the specific star formation rates (SSFRs) of normal, star-forming (main-sequence) galaxies. Indeed, the average SSFRs of AGN hosts are only marginally (20 per cent) lower than those of main-sequence galaxies, with this small deficit being due to a fraction of AGNs residing in quiescent (low-SSFR) galaxies. We estimate 79+/-10 per cent of moderate AGNs are hosted in main-sequence galaxies, 15+/-7 per cent in quiescent galaxies and <10 per cent in strongly starbursting galaxies. The fractions of all main sequence galaxies at z<2 experiencing a period of moderate nuclear activity is strongly dependent on galaxy stellar mass (Mstars); rising from a few per cent at Mstars~10^{10}Msun to >20 per cent at Mstars>10^{11}Msun. Our results indicate that it is galaxy stellar mass that is most important in dictating whether a galaxy hosts a moderate luminosity AGN. We argue that the majority of moderate nuclear activity is fuelled by internal mechanisms rather than violent mergers, suggesting that disk instabilities could be an important AGN feeding mechanism.

We present one of the most precise measurement to date of the spatial clustering of X-ray selected AGNs using a sample derived from the Chandra X-ray Observatory survey in the Bootes field. The real-space two-point correlation function over a redshift interval from z=0.17 to z~3 is well described by the power law, xi(r)=(r/r0)^-gamma, for comoving separations r<~20h^-1 Mpc. We find gamma=1.84+-0.12 and r0 consistent with no redshift trend within the sample (varying between r0=5.5+-0.6 h^-1 Mpc for <z>=0.37 and r0=6.9+-1.0 h^-1 Mpc for <z>=1.28). Further, we are able to measure the projections of the two-point correlation function both on the sky plane and in the line of sight. We use these measurements to show that the Chandra/Bootes AGNs are predominantly located at the centers of dark matter halos with the circular velocity Vmax>320 km/s or M_200 > 4.1e12 h^-1 Msun, and tend to avoid satellite galaxies in halos of this or higher mass. The halo occupation properties inferred from the clustering properties of Chandra/Bootes AGNs --- the mass scale of the parent dark matter halos, the lack of significant redshift evolution of the clustering length, and the low satellite fraction --- are broadly consistent with the Hopkins et al. scenario of quasar activity triggered by mergers of similarly-sized galaxies.

We use the 4Ms CDF-S survey to place direct X-ray constraints on the ubiquity of z~2 heavily obscured AGNs in K<22 BzK galaxies. Forty seven of the 222 BzK galaxies in the central region of the CDF-S are detected at X-ray energies, 11 of which have hard X-ray spectral slopes (Gamma<1) indicating the presence of heavily obscured AGN activity. The other 36 X-ray detected BzK galaxies appear to be relatively unobscured AGNs and starburst galaxies; we use X-ray variability analyses over a rest-frame baseline of ~3 years to further confirm the presence of AGN activity in many of these systems. The majority (7 out of 11) of the heavily obscured AGNs have excess IR emission over that expected from star formation (termed "IR-excess galaxies"). However, we find that X-ray detected heavily obscured AGNs only comprise ~25% of the IR-excess galaxy population, which is otherwise composed of relatively unobscured AGNs and starburst galaxies. We find that the typical X-ray spectrum of the heavily obscured AGNs is better characterized by a pure reflection model than an absorbed power-law model, suggesting extreme Compton-thick absorption in some systems. We verify this result by producing a composite rest-frame 2-20 keV spectrum, which has a similar shape as a reflection-dominated X-ray spectrum and reveals an emission feature at rest-frame energy ~6.4 keV, likely to be due to Fe K. These heavily obscured AGNs are likely to be the distant analogs of the reflection-dominated AGNs recently identified at z~0 with >10 keV observatories. On the basis of these analyses we estimate the space density for typical (intrinsic X-ray luminosities of L_X>1E43 erg/s) heavily obscured and Compton-thick AGNs at z~2. Our space-density constraints are conservative lower limits but they are already consistent with the range of predictions from X-ray background models.

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

We use infrared spectroscopy and photometry to empirically define the intrinsic, thermal infrared spectral energy distribution (i.e., 6-100 um SED) of typical active galactic nuclei (i.e., 2-10 keV luminosity, Lx=10^{42}-10^{44} ergs/s AGNs). On average, the infrared SED of typical AGNs is best described as a broken power-law at <40 um that falls steeply at >40um (i.e., at far-infrared wavelengths). Despite this fall-off at long wavelengths, at least 3 of the 11 AGNs in our sample have observed SEDs that are AGN-dominated even at 60 um, demonstrating the importance of accounting for possible AGN contribution even at far-infrared wavelengths. Our results also suggest that the average intrinsic AGN 6-100 um SED gets bluer with increasing X-ray luminosity, a trend seen both within our sample and also when we compare against the intrinsic SEDs of more luminous quasars (i.e., Lx>10^{44} ergs/s). We compare our intrinsic AGN SEDs with predictions from dusty torus models and find they are more closely matched by clumpy, rather than continuous, torus models. Next, we use our intrinsic AGN SEDs to define a set of correction factors to convert either monochromatic infrared or X-ray luminosities into total intrinsic AGN infrared (i.e., 8-1000 um) luminosities. Finally, we outline a procedure that uses our newly defined intrinsic AGN infrared SEDs, in conjunction with a selection of host-galaxy templates, to fit the infrared photometry of composite galaxies and measure the AGN contribution to their total infrared output. We verify the accuracy of our SED fitting procedure by comparing our results to two independent measures of AGN contribution. Our SED fitting procedure opens up the possibility of measuring the intrinsic AGN luminosities of large numbers of galaxies with well-sampled infrared data (e.g., IRAS, ISO, Spitzer and Herschel).

We present the J-band luminosity function of 1838 mid-infrared and X-ray selected AGNs in the redshift range 0<z<5.85. These luminosity functions are constructed by combining the deep multi-wavelength broad-band observations from the UV to the mid-IR of the NDWFS Bootes field with the X-ray observations of the XBootes survey and the spectroscopic observations of the same field by AGES. Our sample is primarily composed of IRAC-selected AGNs, targeted using modifications of the Stern et al.(2005) criteria, complemented by MIPS 24 microns and X-ray selected AGNs to alleviate the biases of IRAC mid-IR selection against z~4.5 quasars and AGNs faint with respect to their hosts. This sample provides an accurate link between low and high redshift AGN luminosity functions and does not suffer from the usual incompleteness of optical samples at z~3. We find that the space density of the brightest quasars strongly decreases from z=3 to z=0, while the space density of faint quasars is at least flat, and possibly increasing, over the same redshift range. At z>3 we observe a decrease in the space density of quasars of all brightnesses. We model the luminosity function by a double power-law and find that its evolution cannot be described by either pure luminosity or pure density evolution, but must be a combination of both. Our best-fit model has bright and faint power-law indices consistent with the low redshift measurements based on the 2QZ and 2SLAQ surveys and it generally agrees with the number of bright quasars predicted by other LFs at all redshifts. If we construct the QSO luminosity function using only the IRAC-selected AGNs, we find that the biases inherent to this selection method significantly modify the behavior of phi*(z) only for z<1 and have no significant impact upon the characteristic magnitude M*_J(z).

We report the X-ray detection of two z>1.4 infrared-selected galaxy clusters from the IRAC Shallow Cluster Survey (ISCS). We present new data from the Hubble Space Telescope and the W. M. Keck Observatory that spectroscopically confirm cluster ISCS J1432.4+3250 at z=1.49, the most distant of 18 confirmed z>1 clusters in the ISCS to date. We also present new spectroscopy for ISCS J1438.1+3414, previously reported at z = 1.41, and measure its dynamical mass. Clusters ISCS J1432.4+3250 and ISCS J1438.1+3414 are detected in 36ks and 143ks Chandra exposures at significances of 5.2 sigma and 9.7 sigma, from which we measure total masses of log(M_{200,Lx}/Msun) = 14.4 +/- 0.2 and 14.35^{+0.14}_{-0.11}, respectively. The consistency of the X-ray and dynamical properties of these high redshift clusters further demonstrates that the ISCS is robustly detecting massive clusters to at least z = 1.5.

We use a coupled model of the formation and evolution of galaxies and black holes (BH) to study the evolution of active galactic nuclei (AGN) in a cold dark matter universe. The model predicts the BH mass, spin and mass accretion history. BH mass grows via accretion triggered by discs becoming dynamically unstable or galaxy mergers (called the starburst mode) and accretion from quasi-hydrostatic hot gas haloes (called the hot-halo mode). By taking into account AGN obscuration, we obtain a very good fit to the observed luminosity functions (LF) of AGN (optical, soft and hard X-ray, and bolometric) for a wide range of redshifts (0<z<6). The model predicts a hierarchical build up of BH mass, with the typical mass of actively growing BHs increasing with decreasing redshift. Remarkably, despite this, we find downsizing in the AGN population, in terms of the differential growth with redshift of the space density of faint and bright AGN. This arises naturally from the interplay between the starburst and hot-halo accretion modes. The faint end of the LF is dominated by massive BHs experiencing quiescent accretion via a thick disc, primarily during the hot-halo mode. The bright end of the LF, on the other hand, is dominated by AGN which host BHs accreting close to or in excess of the Eddington limit during the starburst mode. The model predicts that the comoving space density of AGN peaks at z~3, similar to the star formation history. However, when taking into account obscuration, the space density of faint AGN peaks at lower redshift (z<2) than that of bright AGN (z~2-3). This implies that the cosmic evolution of AGN is shaped in part by obscuration.

Using a suite of X-ray, mid-IR and optical active galactic nuclei (AGN) luminosity indicators, we search for Compton-thick (CT) AGNs with intrinsic L_X>10^42erg/s at z~0.03-0.2, a region of parameter space which is currently poorly constrained by deep narrow-field and high-energy (E>10keV) all-sky X-ray surveys. We have used the widest XMM-Newton survey (the serendipitous source catalogue) to select a representative sub-sample (14; ~10%) of the 147 X-ray undetected candidate CT AGNs in the Sloan Digital Sky Survey (SDSS) with f_X/f_[OIII]<1; the 147 sources account for ~50% of the overall Type-2 AGN population in the SDSS-XMM overlap region. We use mid-IR spectral decomposition analyses and emission-line diagnostics, determined from pointed Spitzer-IRS spectroscopic observations of these candidate CT AGNs, to estimate the intrinsic AGN emission (predicted L_X,2-10keV (0.2-30)x10^42erg/s). On the basis of the optical [OIII], mid-IR [OIV] and 6um AGN continuum luminosities we conservatively find that the X-ray emission in at least 6/14 (>43%) of our sample appear to be obscured by CT material with N_H>1.5x10^24cm^-2. Under the reasonable assumption that our 14 AGNs are representative of the overall X-ray undetected AGN population in the SDSS-XMM parent sample, we find that >20% of the optical Type-2 AGN population are likely to be obscured by CT material. This implies a space-density of log(Phi) >-4.9Mpc^-3 for CT AGNs with L_X>10^42erg/s at z~0.1, which we suggest may be consistent with that predicted by X-ray background synthesis models. Furthermore, using the 6um continuum luminosity to infer the intrinsic AGN luminosity and the stellar velocity dispersion to estimate M_BH, we find that the most conservatively identified CT AGNs in this sample may harbour some of the most rapidly growing black holes (median M_BH~3x10^7M_o) in the nearby Universe, with a median Eddington ratio of ~0.2.

We present results from spectral and temporal analyses of Suzaku and RXTE observations of the high mass X-ray binary LMC X-4. Using the full 13 years of available RXTE/ASM data, we apply the ANOVA and Lomb normalized Periodogram methods to obtain an improved superorbital period measurement of 30.32 +/- 0.04 days. The phase-averaged X-ray spectra from Suzaku observations during the high state of the superorbital period can be modeled in the 0.6--50 keV band as the combination of a power-law with Gamma ~ 0.6 and a high-energy cutoff at ~ 25 keV, a blackbody with kT_BB ~ 0.18 keV, and emission lines from Fe K_alpha, O VIII, and Ne IX (X Lyalpha). Assuming a distance of 50 kpc, The source has luminosity L_X ~ 3 x 10^38 ergs s^-1 in the 2--50 keV band, and the luminosity of the soft (blackbody) component is L_BB ~ 1.5 x 10^37 ergs s^-1. The energy resolved pulse profiles show single-peaked soft (0.5-1 keV) and hard (6-10 keV) pulses but a more complex pattern of medium (2-10 keV) pulses; cross-correlation of the hard with the soft pulses shows a phase shift that varies between observations. We interpret these results in terms of a picture in which a precessing disk reprocesses the hard X-rays and produces the observed soft spectral component, as has been suggested for the similar sources Her X-1 and SMC X-1.

Combined Spitzer, Chandra, XMM-Newton, and VLA observations of the giant radio galaxy NGC 1316 (Fornax A) show a radio jet and X-ray cavities from AGN outbursts likely triggered by a merger with a late-type galaxy at least 0.4 Gyr ago. We detect a weak nucleus with an SED typical of a low-luminosity AGN with a bolometric luminosity of 2.4x10^42 erg/s. The Spitzer IRAC and MIPS images show dust emission strongest in regions with little or no radio emission. The large infrared luminosity relative to the galaxy's K-band luminosity implies an external origin for the dust. The dust mass implies that the merger spiral galaxy had a stellar mass of 1-6x10^10 M_sun and a gas mass of 2-4x10^9 M_sun. Chandra images show a small ~15"(1.6 kpc) cavity coincident with the radio jet, while the XMM-Newton image shows two large X-ray cavities lying 320"(34.8 kpc) east and west of the nucleus, each approximately 230"(25 kpc) in radius. The radio lobes lie at radii of 14.3'(93 kpc) and 15.6'(101 kcp), more distant from the nucleus than the detected X-ray cavities. The relative morphology of the large scale 1.4 GHz and X-ray emission suggests they were products of two distinct outbursts, an earlier one creating the radio lobes and a later one producing the X-ray cavities. Alternatively, if a single outburst created both the X-ray cavities and the radio lobes, this would require that the morphology is not fully defined by the 1.4 GHz emission. For the more likely two outburst scenario, we use the buoyancy rise times to estimate an age for the more recent outburst of 0.1 Gyr and use the associated PV work done by the expanding plasma to create the X-ray cavities to estimate the outburst's energy of 10^58 ergs. The present size and location of the large radio lobes implies that the AGN outburst that created them happened ~0.4 Gyr ago and released ~5x10^58 ergs. (abridged)

We present a set of low resolution empirical SED templates for AGNs and galaxies in the wavelength range from 0.03 to 30 microns based on the multi-wavelength photometric observations of the NOAO Deep-Wide Field Survey Bootes field and the spectroscopic observations of the AGN and Galaxy Evolution Survey. Our training sample is comprised of 14448 galaxies in the redshift range 0<~z<~1 and 5347 likely AGNs in the range 0<~z<~5.58. We use our templates to determine photometric redshifts for galaxies and AGNs. While they are relatively accurate for galaxies, their accuracies for AGNs are a strong function of the luminosity ratio between the AGN and galaxy components. Somewhat surprisingly, the relative luminosities of the AGN and its host are well determined even when the photometric redshift is significantly in error. We also use our templates to study the mid-IR AGN selection criteria developed by Stern et al.(2005) and Lacy et al.(2004). We find that the Stern et al.(2005) criteria suffers from significant incompleteness when there is a strong host galaxy component and at z =~ 4.5, when the broad Halpha emission line is redshifted into the [3.6] band, but that it is little contaminated by low and intermediate redshift galaxies. The Lacy et al.(2004) criterion is not affected by incompleteness at z =~ 4.5 and is somewhat less affected by strong galaxy host components, but is heavily contaminated by low redshift star forming galaxies. Finally, we use our templates to predict the color-color distribution of sources in the upcoming WISE mission and define a color criterion to select AGNs analogous to those developed for IRAC photometry. We estimate that in between 640,000 and 1,700,000 AGNs will be identified by these criteria, but will have serious completeness problems for z >~ 3.4.

Galaxy groups are key tracers of galaxy evolution, cluster evolution, and structure formation, yet they are difficult to study at even moderate redshift. We have undertaken a project to observe a flux-limited sample of intermediate-redshift (0.1 < z < 0.5) group candidates identified by the XBootes Chandra survey. When complete, this project will nearly triple the current number of groups with measured temperatures in this redshift range. Here we present deep Suzaku/XIS and Chandra/ACIS follow-up observations of the first 10 targets in this project; all are confirmed sources of diffuse, thermal emission with derived temperatures and luminosities indicative of rich groups/poor clusters. By exploiting the multi-wavelength coverage of the XBootes/NOAO Deep Wide Field Survey (NDWFS) field, we aim to (1) constrain non-gravitational effects that alter the energetics of the intragroup medium, and (2) understand the physical connection between the X-ray and optical properties of groups. We discuss the properties of the current group sample in the context of observed cluster scaling relations and group and cluster evolution and outline the future plans for this project.

Large extragalactic surveys allow us to trace, in a statistical sense, how supermassive black holes, their host galaxies, and their dark matter halos evolve together over cosmic time, and so explore the consequences of AGN feedback on galaxy evolution. Recent studies have found significant links between the accretion states of black holes and galaxy stellar populations, local environments, and obscuration by gas and dust. This article describes some recent results and shows how such studies may provide new constraints on models of the co-evolution of galaxies and their central SMBHs. Finally, I discuss observational prospects for the proposed Wide-Field X-ray Telescope mission.

At low Eddington ratio (mdot), two effects make it harder to detect AGN given some selection criteria. First, even with fixed accretion physics, AGN are diluted/less luminous relative to their hosts; the magnitude of this depends on host properties and so on luminosity and redshift. Second, they may transition to a radiatively inefficient state, changing SED shape and dramatically decreasing in optical/IR luminosity. These effects lead to differences in observed AGN samples, even at fixed bolometric luminosity and after correction for obscuration. The true Eddington ratio distribution may depend strongly on luminosity, but this will be seen only in surveys robust to dilution and radiative inefficiency (X-ray or narrow-line samples); selection effects imply that AGN in optical samples will have uniformly high mdot. This also implies that different selection methods yield systems with different hosts: the clustering of faint optical/IR sources will be weaker than that of X-ray sources, and optical/IR Seyferts will reside in more disk-dominated galaxies while X-ray selected Seyferts will preferentially occupy early-type systems. If observed mdot distributions are correct, a large fraction of low-luminosity AGN currently classified as 'obscured' are in fact diluted and/or radiatively inefficient, not obscured by gas or dust. This is equally true if X-ray hardness is used as a proxy for obscuration, since radiatively inefficient SEDs near mdot~0.01 are X-ray hard. These effects can explain most of the claimed luminosity/redshift dependence in the 'obscured' AGN population, with the true obscured fraction as low as 20%.

We have searched for active galactic nuclei (AGNs) in K+A galaxies, using multiwavelength imaging and spectroscopy in the Bootes field of the NOAO Deep Wide-Field Survey. The K+A galaxies, which have had their star formation rapidly truncated, are selected via their strong Balmer absorption lines and weak H-alpha emission. Our sample consists of 24 K+A galaxies selected from 6594 0.10<z<0.35 galaxies brighter than I=20 with optical spectroscopy from the AGN and Galaxy Evolution Survey. Two thirds of the K+A galaxies are likely ongoing galaxy mergers, with nearby companion galaxies or tidal tails. Galaxy mergers may be responsible for the truncation of star formation, or we are observing the aftermath of merger triggered starbursts. As expected, the optical colors of K+A galaxies largely fall between blue galaxies with ongoing star formation and red passive galaxies. However, only 1% of the galaxies with colors between the red and blue populations are K+A galaxies, and we conclude that the truncation of star formation in K+A galaxies must have been unusually abrupt (<100 Myr). We examined the AGN content of K+A galaxies with both optical emission-line ratios (BPT diagrams) and Chandra X-ray imaging. At least half of all K+A galaxies display the optical emission-line ratios of AGNs, and a third of M_R<-22 K+A galaxies host AGNs with X-ray luminosities of 10^{42} erg/s. The faintest K+A galaxies do not show clear evidence for hosting AGNs, having emission-line ratios consistent with photoionization by massive stars and few X-ray detections. We speculate that two mechanisms may be responsible for the truncation of star formation in K+A galaxies, with AGN feedback only playing a role in M_R<-20.5 galaxies.

Extragalactic X-ray surveys are exceptionally powerful tools for studying the evolution of supermassive black holes and their host galaxies, by detecting large numbers of active galactic nuclei (AGN) and star-forming galaxies over a wide range of redshifts and cosmic environments. With its sensitivity and superb angular resolution, Chandra has been at the forefront of recent extragalactic surveys. This article provides a brief overview of Chandra surveys, and highlights a few recent results on the composition of the cosmic X-ray background, insights on X-ray source populations, and links between AGN and galaxy evolution, as well as discussing prospects for surveys with future X-ray missions.

We explore the connection between different classes of active galactic nuclei (AGNs) and the evolution of their host galaxies, by deriving host galaxy properties, clustering, and Eddington ratios of AGNs selected in the radio, X-ray, and infrared. We study a sample of 585 AGNs at 0.25 < z < 0.8 using redshifts from the AGN and Galaxy Evolution Survey (AGES) and data in the radio (WSRT 1.4 GHz), X-rays (Chandra XBootes), and mid-IR (IRAC Shallow Survey). The radio, X-ray, and IR AGN samples show modest overlap, indicating that to the flux limits of the survey, they represent largely distinct classes of AGNs. We derive host galaxy colors and luminosities, as well as Eddington ratios (lambda), for obscured or optically faint AGNs. We also measure the two-point cross-correlation between AGNs and galaxies on scales of 0.3-10 h^-1 Mpc, and derive typical dark matter halo masses. We find that: (1) radio AGNs are mainly found in luminous red galaxies, are strongly clustered (with M_halo ~ 3x10^13 h^-1 M_sun), and have very low lambda <~ 10^-3; (2) X-ray-selected AGNs are preferentially found in galaxies in the "green valley" of color-magnitude space and are clustered similarly to typical AGES galaxies (M_halo ~ 10^13 h^-1 M_sun), with 10^-3 <~ lambda <~ 1; (3) IR AGNs reside in slightly bluer, less luminous galaxies than X-ray AGNs, are weakly clustered (M_halo <~ 10^12 h^-1 M_sun), and have lambda > 10^-2. We interpret these results in terms of a simple model of AGN and galaxy evolution, whereby a "quasar" phase and the growth of the stellar bulge occurs when a galaxy's dark matter halo reaches a critical mass between ~10^12 and 10^13 M_sun. Subsequently, star formation ceases and AGN accretion shifts from radiatively efficient (optical- and IR- bright) to radiatively inefficient (optically-faint, radio-bright) modes.

We discuss the central role played by X-ray studies to reconstruct the past history of formation and evolution of supermassive Black Holes (BHs), and the role they played in shaping the properties of their host galaxies. We shortly review the progress in this field contributed by the current X-ray and multiwavelength surveys. Then, we focus on the outstanding scientific questions that have been opened by observations carried out in the last years and that represent the legacy of Chandra and XMM, as for X-ray observations, and the legacy of the SDSS, as for wide area surveys: 1) When and how did the first supermassive black holes form? 2) How does cosmic environment regulate nuclear activity (and star formation) across cosmic time? 3) What is the history of nuclear activity in a galaxy lifetime? We show that the most efficient observational strategy to address these questions is to carry out a large-area X-ray survey, reaching a sensitivity comparable to that of deep Chandra and XMM pointings, but extending over several thousands of square degrees. Such a survey can only be carried out with a Wide-Field X-ray Telescope (WFXT) with a high survey speed, due to the combination of large field of view and large effective area, i.e., grasp, and sharp PSF. We emphasize the important synergies that WFXT will have with a number of future groundbased and space telescopes, covering from the radio to the X-ray bands and discuss the immense legacy value that such a mission will have for extragalactic astronomy at large.

We combine IR, optical and X-ray data from the overlapping, 9.3 square degree NOAO Deep Wide-Field Survey (NDWFS), AGN and Galaxy Evolution Survey (AGES), and XBootes Survey to measure the X-ray evolution of 6146 normal galaxies as a function of absolute optical luminosity, redshift, and spectral type over the largely unexplored redshift range 0.1 < z < 0.5. Because only the closest or brightest of the galaxies are individually detected in X-rays, we use a stacking analysis to determine the mean properties of the sample. Our results suggest that X-ray emission from spectroscopically late-type galaxies is dominated by star formation, while that from early-type galaxies is dominated by a combination of hot gas and AGN emission. We find that the mean star formation and supermassive black hole accretion rate densities evolve like (1+z)^3, in agreement with the trends found for samples of bright, individually detectable starburst galaxies and AGN. Our work also corroborates the results of many previous stacking analyses of faint source populations, with improved statistics.

We explore the connection between different classes of active galactic nuclei (AGNs) and the evolution of their host galaxies, by deriving host galaxy properties, clustering, and Eddington ratios of AGNs selected in the radio, X-ray, and infrared (IR) wavebands from the wide-field (9 deg^2) Bootes survey. We study a sample of 585 AGNs at 0.25 < z < 0.8 using redshifts from the AGN and Galaxy Evolution Survey (AGES). We find that radio and X-ray AGNs reside in relatively large dark matter halos (M_halo ~ 3x10^13 and 10^13 h^-1 M_sun, respectively) and are found in galaxies with red and "green" colors. In contrast, IR AGNs are in less luminous galaxies, have higher Eddington ratios, and reside in halos with M_halo < 10^12 M_sun. We interpret these results in terms of a general picture for AGNs and galaxy evolution, in which quasar activity is triggered when M_halo ~ 10^12 - 10^13 M_sun, after which star formation ceases and AGN accretion shifts to optically-faint, X-ray and radio-bright modes.

We have undertaken a long-term project, Planets in Stellar Clusters Extensive Search (PISCES), to search for transiting planets in open clusters. In this paper we present the results for NGC 188, an old, rather populous cluster. We have monitored the cluster for more than 87 hours, spread over 45 nights. We have not detected any good transiting planet candidates. We have discovered 18 new variable stars in the cluster, bringing the total number of identified variables to 46, and present for them high precision light curves, spanning 15 months.

We identify a population of 640 obscured and 839 unobscured AGNs at redshifts 0.7<z<~3 using multiwavelength observations of the 9 deg^2 NOAO Deep Wide-Field Survey (NDWFS) region in Bootes. We select AGNs on the basis of Spitzer IRAC colors obtained by the IRAC Shallow Survey. Redshifts are obtained from optical spectroscopy or photometric redshift estimators. We classify the IR-selected AGNs as IRAGN 1 (unobscured) and IRAGN 2 (obscured) using a simple criterion based on the observed optical to mid-IR color, with a selection boundary of R-[4.5]=6.1, where R and [4.5] are the Vega magnitudes in the R and IRAC 4.5 micron bands, respectively. We verify this selection using X-ray stacking analyses with data from the Chandra XBootes survey, as well as optical photometry from NDWFS and spectroscopy from MMT/AGES. We show that (1) these sources are indeed AGNs, and (2) the optical/IR color selection separates obscured sources (with average N_H~3x10^22 cm^-2 obtained from X-ray hardness ratios, and optical colors and morphologies typical of galaxies) and unobscured sources (with no X-ray absorption, and quasar colors and morphologies), with a reliability of >~80%. The observed numbers of IRAGNs are comparable to predictions from previous X-ray, optical, and IR luminosity functions, for the given redshifts and IRAC flux limits. We observe a bimodal distribution in R-[4.5] color, suggesting that luminous IR-selected AGNs have either low or significant dust extinction, which may have implications for models of AGN obscuration.

The deepest extragalactic X-ray observation, the 2 Ms Chandra Deep Field North (CDF-N), resolves ~80% of the total extragalactic cosmic X-ray background (CXB) in the 1-2 keV band. Recent work has shown that 70% of the remaining CXB flux is associated with sources detected by the Hubble Space Telescope (HST). This paper uses the existing CDF-N data to constrain the X-ray flux distribution of these X-ray undetected HST sources, by comparing the number of 0.5-2 keV X-ray counts at the HST positions to those expected for model flux distributions. In the simple case where all the undetected HST X-ray sources have the same 0.5-2 keV flux, the data are best fit by 1.5-3 counts per source in 2 Ms, compared to a detection limit (at 10% completeness) of 9 counts. Assuming a more realistic power-law logN-logS distribution [N(>S) S^-alpha], the data favor a relatively steep flux distribution, with alpha=1.1^+0.5_-0.3 (limits are 99% confidence). This slope is very similar to that previously found for faint normal and starburst galaxies in the CDF-N. These results suggest deeper Chandra observations will detect a new population of faint X-ray sources, but extremely deep exposures are needed to resolve the remainder of the soft CXB. In the most optimistic scenario, when the HST sources have the flattest allowed flux distribution and all the sources without HST counterparts are detected, observations 5 times more sensitive than the existing ones would resolve at most ~60% of the remaining soft CXB.

We present a measurement of the surface brightness of the cosmic X-ray background (CXB) in the Chandra Deep Fields, after excluding all detected X-ray, optical and infrared sources. The work is motivated by a recent X-ray stacking analysis by Worsley and collaborators, which showed that galaxies detected by HST but not by Chandra may account for most of the unresolved CXB at E>1 keV. We find that after excluding HST and Spitzer IRAC sources, some CXB still remains, but it is marginally significant: (3.4+/-1.4)x10^-13 ergs cm^-2 s^-1 deg^2 in the 1-2 keV band and (4+/-9)x10^-13 ergs cm^-2 s^-1 deg^2 in the 2-5 keV band, or 7%+/-3% and 4%+/-9% of the total CXB, respectively. Of the 1-2 keV signal resolved by the HST sources, 34%+/-2% comes from objects with optical colors typical of ``normal'' galaxies (which make up 25% of the HST sources), while the remaining flux comes from objects with colors of starburst and irregular galaxies. In the 0.65-1 keV band (just above the bright Galactic O VII line) the remaining diffuse intensity is (1.0+/-0.2)x10^-12 ergs cm^-2 s^-1 deg^2. This flux includes emission from the Galaxy as well as from the hypothetical warm-hot intergalactic medium (WHIM), and provides a conservative upper limit on the WHIM signal that comes interestingly close to theoretical predictions.

We present upper limits on line emission in the Cosmic X-ray background (CXB) that would be produced by decay of sterile neutrino dark matter. We employ the spectra of the unresolved component of the CXB in the Chandra Deep Fields North and South obtained with the Chandra CCD detector in the E=0.8-9 keV band. The expected decay flux comes from the dark matter on the lines of sight through the Milky Way galactic halo. Our constraints on the sterile neutrino decay rate are sensitive to the modeling of the Milky Way halo. The highest halo mass estimates provide a limit on the sterile neutrino mass of m_s<2.9 keV in the Dodelson-Widrow production model, while the lowest halo mass estimates provide the conservative limit of m_s<5.7 keV (2-sigma). We also discuss constraints from a short observation of the softer (E<1 keV) X-ray background with a rocket-borne calorimeter by McCammon and collaborators.

We use observations from the 9 square degree multiwavelength survey in Bootes to identify hundreds of obscured active galactic nuclei (AGN) with high redshifts (z > 0.7), luminosities (L_bol > 10^45 ergs/s), and moderate obscuring columns (N_H > 10^22 cm^-2), and to measure the clustering properties of X-ray AGN at z > 1. In the Bootes region, shallow (5 ks) Chandra X-ray observations have detected ~4,000 X-ray sources, and the same region has been mapped with deep optical imaging and by Spitzer IRAC, which detects ~300,000 point sources, of which ~30,000 have detections in all four IRAC bands, for which we can select AGN on the basis of their mid-IR colors. With the MMT/Hectospec we have obtained modest resolution optical spectra for about half the X-ray sources (out to z > 3) and ~20,000 galaxies (out to z = 0.7). With this multiwavelength data we select >400 AGN per square degree (compared to 12 per square degree from SDSS). Among a sample of IRAC-selected AGN we identify 641 candidate obscured objects based on their R band and IRAC luminosities. We use X-ray stacking techniques to verify that they are obscured AGN and measure their absorbing column densities. We also measure the three-dimensional two-point correlation function for X-ray selected AGN.

We present the absolute measurement of the unresolved 0.5-8 keV cosmic X-ray background (CXB) in the Chandra Deep Fields (CDFs) North and South, the longest observations with Chandra (2 Ms and 1 Ms, respectively). We measure the unresolved CXB intensity by extracting spectra of the sky, removing all point and extended sources detected in the CDF. To model and subtract the instrumental background, we use observations obtained with ACIS in stowed position, not exposed to the sky. The unresolved signal in the 0.5-1 keV band is dominated by diffuse Galactic and local thermal-like emission. In the 1-8 keV band, the unresolved spectrum is adequately described by a power law with a photon index 1.5. We find unresolved CXB intensities of (1.04+/-0.14)x10^-12 ergs cm^-2 s^-1 deg^-2 for the 1-2 keV band and (3.4+/-1.7)x10^-12 ergs cm^-2 s^-1 deg^-2 for the 2-8 keV band. Our detected unresolved intensities in these bands significantly exceed the expected flux from sources below the CDF detection limits, if one extrapolates the logN/logS curve to zero flux. Thus these background intensities imply either a genuine diffuse component, or a steepening of the logN/logS curve at low fluxes, most significantly for energies <2 keV. Adding the unresolved intensity to the total contribution from sources detected in these fields and wider-field surveys, we obtain a total intensity of the extragalactic CXB of (4.6+/-0.3)x10^-12 ergs cm^-2 s^-1 deg^-2 for 1-2 keV and (1.7+/-0.2)x10^-11 ergs cm^-2 s^-1 deg^-2 for 2-8 keV. These totals correspond to a CXB power law normalization (for photon index 1.4) of 10.9 photons cm^-2 s^-1 keV^-1 sr^-1 at 1 keV. This corresponds to resolved fracations of 77+/-3% and 80+/-8% for 1-2 and 2-8 keV, respectively.

We examine the X-ray and infrared properties of galaxies and AGNs in the 9 square degree Bootes field, using data from the Chandra XBootes and Spitzer IRAC Shallow Surveys, as well as optical spectroscopy from the AGES survey. A sample of ~30,000 objects are detected in all four IRAC bands, of which ~2,000 are associated with X-ray sources. We also study X-ray fainter sources using stacking techniques, and find that X-ray fluxes are highest for objects with IRAC colors that are known to be characteristic of AGNs. Because these are shallow, wide-field surveys, they probe the bright end of the AGNluminosity function out to spectroscopic redshifts as high as z=3-4. We can use this multiwavelength dataset to explore the properties and redshift evolution of a large sample of luminous active galaxies.

We have undertaken a long-term project, Planets in Stellar Clusters Extensive Search (PISCES), to search for transiting planets in open clusters. In this paper we present the results for NGC 2158, an intermediate age, populous cluster. We have monitored the cluster for over 260 hours, spread over 59 nights. We have detected one candidate transiting low luminosity object, with eclipse depth of 3.7% in the R-band. If the host star is a member of the cluster, the eclipse depth is consistent with a 1.7 R_J object. Cluster membership of the host is supported by its location on the cluster main sequence (MS) and its close proximity to the cluster center (2'). We have discovered two other stars exhibiting low-amplitude (4-5%) transits, V64 and V70, but they are most likely blends or field stars. Given the photometric precision and temporal coverage of our observations, and current best estimates for the frequency and radii of short-period planets, the expected number of detectable transiting planets in our sample is 0.13. We have observed four outbursts for the candidate cataclysmic variable V57. We have discovered 40 new variable stars in the cluster, bringing the total number of identified variables to 97, and present for them high precision light curves, spanning 13 months.

We present pulse-phase X-ray spectroscopy of the high-mass X-ray pulsar SMC X-1 from five different epochs, using Chandra X-ray Observatory ACIS-S and XMM-Newton EPIC-pn data. The X-ray spectrum consistently shows two distinct components, a hard power law and a soft blackbody with kT_BB ~ 0.18 keV. For all five epochs the hard component shows a simple double-peaked pulse shape, and also a variation in the power law slope, which becomes harder at maximum flux and softer at minimum flux. For the soft component, the pulse profile changes between epochs, in both shape and phase relative to the power law pulses. The soft component is likely produced by reprocessing of the hard X-ray pulsar beam by the inner accretion disk. We use a model of a twisted inner disk, illuminated by the rotating X-ray pulsar beam, to simulate pulsations in the soft component due to this reprocessing. We find that for some disk and beam geometries, precession of an illuminated accretion disk can roughly reproduce the observed long-term changes in the soft pulse profiles.

We present high-precision photometry of the hypervelocity star SDSS J090745.0+024507 (HVS), which has a Galactic rest-frame radial velocity of v=709 km/s, and so has likely been ejected from the supermassive black hole in the Galactic center. Our data were obtained on two nights using the MMT 6.5m telescope, and is supplemented by lower precision photometry obtained on four nights using the FLWO 1.2m telescope. The high-precision photometry indicates that the HVS is a short-period, low-amplitude variable, with period P=0.2-2 days and amplitude A = 2-10%. Together with the known effective temperature of T_eff ~ 10,500 K (spectral type B9), this variability implies that the HVS is a member of the class of slowly pulsating B-type main sequence stars, thus resolving the previously-reported two-fold degeneracy in the luminosity and distance of the star. The HVS has a heliocentric distance of 71 kpc, and an age of ~0.35 Gyr. The time of ejection from the center of the Galaxy is < 100 Myr, and thus the existence of the OS constitutes observational evidence of a population of young stars in the proximity of the central supermassive black hole ~0.1 Gyr ago. It is possible that the HVS was a member of a binary that was tidally disrupted by the central black hole; we discuss constraints on the properties of the companion's orbit.

We present the results of timing and spectral analysis of X-ray high state observations of the high-mass X-ray pulsar SMC X-1 with Chandra, XMM-Newton, and ROSAT, taken between 1991 and 2001. The source has L_X ~ 3-5 x 10^38 ergs/s, and the spectra can be modeled as a power law plus blackbody with kT_BB \~ 0.18 keV and reprocessed emission radius R_BB ~ 2 x 10^8 cm, assuming a distance of 60 kpc to the source. Energy-resolved pulse profiles show several distinct forms, more than half of which include a second pulse in the soft profile, previously documented only in hard energies. We also detect significant variation in the phase shift between hard and soft pulses, as has recently been reported in Her X-1. We suggest an explanation for the observed characteristics of the soft pulses in terms of precession of the accretion disk.

The spectra of many X-ray pulsars show, in addition to a power law, a low-energy component that has often been modeled as a blackbody with kT ~ 0.1 keV. However the physical origin of this soft excess has remained a mystery. We examine a sample of well-studied, bright X-ray pulsars, which have been observed using ROSAT, ASCA, Ginga, RXTE, BeppoSAX, Chandra, and XMM-Newton. In particular we consider the Magellanic Cloud pulsars SMC X-1, LMC X-4, XTE J0111.2-7317, and RX J0059.2-7138 and the Galactic sources Her X-1, 4U 1626-67, Cen X-3, and Vela X-1. We show that the soft excess is a very common if not ubiquitous feature intrinsic to X-ray pulsars. We evaluate several possible mechanisms for the soft emission, using theoretical arguments as well as observational clues such as spectral shapes, eclipses, pulsations of the soft component, and superorbital modulation of the source flux. We find that reprocessing of hard X-rays from the neutron star by the inner region of the accretion disk is the only process that can explain the soft excess in all the pulsars with Lx > 10^38 ergs/s. Other mechanisms, such as emission from diffuse gas in the system, are important in less luminous objects.