The unique combination of superb spatial resolution, wide-area coverage, and deep depth of the optical imaging from the Hyper Suprime-Cam (HSC) Subaru Strategic Program is utilized to search for dual quasar candidates. Using an automated image analysis routine on 34,476 known SDSS quasars, we identify those with two (or more) distinct optical point sources in HSC images covering 796 deg^2. We find 421 candidates out to a redshift of 4.5 of which one hundred or so are more likely after filtering out contaminating stars. Angular separations of 0.6 - 4.0" correspond to projected separations of 3 - 30 kpc, a range relatively unexplored for population studies of luminous dual quasars. Using Keck-I/LRIS and Gemini-N/NIFS, we spectroscopically confirm three dual quasar systems at z < 1, two of which are previously unknown out of eight observed, based on the presence of characteristic broad emission lines in each component, while highlighting that the continuum of one object in one of the pairs is reddened. In all cases, the [OIII]5007 emission lines have mild velocity offsets, thus the joint [OIII] line profile is not double-peaked. We find a dual quasar fraction of 0.26+/-0.18% and no evidence for evolution. A comparison with the Horizon-AGN simulation seems to support the case of no evolution in the dual quasar fraction when broadly matching the quasar selection. These results may indicate a scenario in which the frequency of the simultaneous triggering of luminous quasars is not as sensitive as expected to the cosmic evolution of the merger rate or gas content of galaxies.

We study the incidence of nuclear activity in a large sample of massive post-starburst galaxies at z~0.7 selected from the Sloan Digital Sky Survey, and identify active galactic nuclei based on radio continuum and optical emission lines. Over our mass range of 10^10.6-10^11.5 Msun, the incidence of radio activity is weakly dependent on stellar mass and independent of stellar age, while radio luminosity depends strongly on stellar mass. Optical nuclear activity incidence depends most strongly on the Dn4000 line index, a proxy for stellar age, with an active fraction that is ~ten times higher in the youngest versus oldest post-starburst galaxies. Since a similar trend is seen between age and molecular gas fractions, we argue that, like in local galaxies, the age trend reflects a peak in available fueling rather than feedback from the central black hole on the surrounding galaxy.

Galactic-scale winds driven by active galactic nuclei (AGN) are often invoked to suppress star formation in galaxy evolution models, but the mechanisms driving these outflows are hotly debated. Two key AGN feedback models are (1) radiation pressure accelerating cool gas and (2) a hot outflowing wind entraining the ISM. Highly ionized emission-line diagnostics represent a powerful means of differentiating these scenarios because of their sensitivity to the expected compression of the ISM clouds by the hot wind. Here, we report the first spatially resolved UV emission spectroscopy of a prototypical (radio-quiet) quasar-driven superwind around the obscured quasar SDSSJ1356+1026 at z=0.123. We observe ratios of OVI/CIV, NV/CIV, and CIV/HeII that are remarkably similar for outflowing gas clouds <100 pc and ~10 kpc from the nucleus. Such similarity is expected for clouds with AGN radiation pressure dominated dynamics. Comparing the observed line emission to models of clouds in balance with radiation pressure and/or a hot wind, we rule out the presence of a dynamically important hot wind and constrain the ratio of hot gas pressure to radiation pressure to P_hot/P_rad<0.25 both at <100 pc and ~10 kpc from the nucleus. Moreover, the predictions of the radiation pressure confined cloud models that best fit observed UV line ratios are consistent with the observed diffuse X-ray spectrum. These results indicate that this AGN superwind is driven by radiation pressure or was driven by a hot wind that has since dissipated despite on-going AGN activity.

We present an analysis of 12 optically selected dual AGN candidates at $z < 0.34$. Each candidate was originally identified via double-peaked [O III] $\lambda$5007 emission lines, and have received follow-up $Chandra$ and $HST$ observations. Because the X-ray data are low-count ($<100$ counts) with small separations ($<1$"), a robust analysis is necessary for classifying each source. Pairing long-slit [O III] observations with existing $Chandra$ observations, we re-analyze the X-ray observations with ${\tt BAYMAX}$ to determine whether the X-ray emission from each system is more likely a single or dual point source. We find that 4 of the 12 sources are likely dual X-ray point source systems. We examine each point source's spectra via a Monte Carlo method that probabilistically identifies the likely origin of each photon. When doing so, we find that (i) the secondary X-ray point sources in 2 of the systems have $L_{\mathrm{X}}<10^{40}$ erg s$^{-1}$, such that we cannot rule out a non-AGN origin, (ii) one source has a secondary with $L_{\mathrm{X}}>10^{40}$ erg s$^{-1}$ but a spectrum that is too soft to definitively preclude being X-ray emitting diffuse gas that was photoionized by the primary AGN, and (iii) one system (SDSS J1126+2944) is a dual AGN. Additionally, using complementary $HST$ observations, we analyze a sub-sample of systems that are visually identified as merging. Our results suggest that dual AGNs may preferentially reside in mergers with small separations, consistent with both simulations and observations.

Cosmic voids, the under-dense regions of the cosmic web, are widely used to constrain cosmology. Voids contain few, isolated galaxies, presumably expected to be less evolved and preserving memory of the pristine Universe. We use the cosmological hydrodynamical simulation Horizon-AGN coupled to the void finder {\sc \texttt{VIDE}} to investigate properties of galaxies in voids at z=0. We find that, closer to void centers, low-mass galaxies are more common than their massive counterparts. At fixed dark matter halo mass, they have smaller stellar masses than in denser regions. The star formation rate of void galaxies diminishes when approaching void centers, but their sSFR slightly increases, suggesting that void galaxies form stars more efficiently with respect to their stellar mass. We find that this can not only be attributed to the prevalence of low-mass galaxies. The inner region of voids also predominantly host low-mass BHs. However, the BH mass to galaxy mass ratios resemble those of the whole simulation at z=0. Our results suggest that even if the growth channels in cosmic voids are different than in denser environments, voids grow their galaxies and BHs in a similar way. While a large fraction of the BHs have low Eddington ratios, we find that 20\% could be observed as AGN with log10 L=41.5-42.5 erg/s in hard X-ray (2-10 keV). These results pave the way to future work with larger next-generation hydro simulations, aiming to confirm our findings and prepare the application on data from upcoming large surveys such as PFS, Euclid and WFIRST.

Like massive galaxies, dwarf galaxies are expected to undergo major mergers with other dwarfs. However, the end state of these mergers and the role that merging plays in regulating dwarf star formation is uncertain. Using imaging from the Hyper Suprime-Cam Subaru Strategic program, we construct a sample of dwarf-dwarf mergers and examine the star formation and host properties of the merging systems. These galaxies are selected via an automated detection algorithm from a sample of 6875 spectroscopically selected isolated dwarf galaxies at $z<0.12$ and $\log(M_\star/M_\odot)<9.6$ from the Galaxy and Mass Assembly (GAMA) and Sloan Digital Sky Survey (SDSS) spectroscopic campaigns. We find a total tidal feature detection fraction of 3.29% (6.1% when considering only galaxies at $z<0.05$). The tidal feature detection fraction rises strongly as a function of star formation activity; 15%-20% of galaxies with extremely high H$\alpha$ equivalent width (H$\alpha$ EW > 250 Angstrom) show signs of tidal debris. Galaxies that host tidal debris are also systematically bluer than the average galaxy at fixed stellar mass. These findings extend the observed dwarf-dwarf merger sequence with a significant sample of dwarf galaxies, indicating that star formation triggered in mergers between dwarf galaxies continues after coalescence.

The properties of low-mass galaxies hosting central black holes provide clues about the formation and evolution of the progenitors of supermassive black holes. In this letter, we present HSC-XD 52, a spectroscopically confirmed low-mass active galactic nucleus (AGN) at an intermediate redshift of $z\sim0.56$. We detect this object as a very luminous X-ray source coincident with a galaxy observed by the Hyper Suprime-Cam (HSC) as part of a broader search for low-mass AGN. We constrain its stellar mass through spectral energy distribution modeling to be LMC-like at $M_\star \approx 3 \times 10^9 M_\odot$, placing it in the dwarf regime. We estimate a central black hole mass of $M_\mathrm{BH} \sim 10^{6} M_\odot$. With an average X-ray luminosity of $L_X \approx 3.5 \times 10^{43}~\mathrm{erg}~\mathrm{s}^{-1}$, HSC-XD 52 is among the most luminous X-ray selected AGN in dwarf galaxies. The spectroscopic and photometric properties of HSC-XD 52 indicate that it is an intermediate redshift counterpart to local low-mass AGN.

The relationship between galaxies and the state/chemical enrichment of the warm-hot intergalactic medium (WHIM) expected to dominate the baryon budget at low-z provides sensitive constraints on structure formation and galaxy evolution models. We present a deep redshift survey in the field of 1ES1553+113, a blazar with a unique combination of UV+X-ray spectra for surveys of the circum-/intergalactic medium (CGM/IGM). Nicastro et al. 2018 reported the detection of two O VII WHIM absorbers at $z=0.4339$ and $0.3551$ in its spectrum, suggesting that the WHIM is metal-rich and sufficient to close the missing baryons problem. Our survey indicates that the blazar is a member of a $z=0.433$ group and that the higher-$z$ O VII candidate arises from its intragroup medium. The resulting bias precludes its use in baryon censuses. The $z=0.3551$ candidate occurs in an isolated environment 630 kpc from the nearest galaxy (with stellar mass $\log M_*/M_\odot \approx 9.7$) which we show is unexpected for the WHIM. Finally, we characterize the galactic environments of broad H I Ly$\alpha$ absorbers (Doppler widths of $b=40-80$ \kms; $T\lesssim4\times10^5$ K) which provide metallicity independent WHIM probes. On average, broad Ly$\alpha$, absorbers are ${\approx}2\times$ closer to the nearest luminous ($L>0.25 L_*$) galaxy (700 kpc) than narrow ($b<30$ \kms; $T\lesssim4\times10^5$ K) ones (1300 kpc) but ${\approx}2\times$ further than O\,VI absorbers (350 kpc). These observations suggest that gravitational collapse heats portions of the IGM to form the WHIM but with feedback that does not enrich the IGM far beyond galaxy/group halos to levels currently observable in UV/X-ray metal lines.

This paper presents the second data release of the Hyper Suprime-Cam Subaru Strategic Program, a wide-field optical imaging survey on the 8.2 meter Subaru Telescope. The release includes data from 174 nights of observation through January 2018. The Wide layer data cover about 300 deg^2 in all five broadband filters (grizy) to the nominal survey exposure (10min in gr and 20min in izy). Partially observed areas are also included in the release; about 1100 deg^2 is observed in at least one filter and one exposure. The median seeing in the i-band is 0.6 arcsec, demonstrating the superb image quality of the survey. The Deep (26 deg^2) and UltraDeep (4 deg^2) data are jointly processed and the UltraDeep-COSMOS field reaches an unprecedented depth of i~28 at 5 sigma for point sources. In addition to the broad-bands, narrow-band data are also available in the Deep and UltraDeep fields. This release includes a major update to the processing pipeline, including improved sky subtraction, PSF modeling, object detection, and artifact rejection. The overall data quality has been improved, but this release is not without problems; there is a persistent deblender problem as well as new issues with masks around bright stars. The user is encouraged to review the issue list before utilizing the data for scientific explorations. All the image products as well as catalog products are available for download. The catalogs are also loaded to a database, which provides an easy interface for users to retrieve data for objects of interest. In addition to these main data products, detailed galaxy shape measurements withheld from the Public Data Release 1 (PDR1) are now available to the community. The shape catalog is drawn from the S16A internal release, which has a larger area than PDR1 (160 deg^2). All products are available at the data release site, https://hsc-release.mtk.nao.ac.jp/.

Supermassive black hole (SMBH) binaries with masses of ~10^8--10^9 Msun are expected to dominate the contribution to the as-yet undetected gravitational wave background (GWB) signal at the nanohertz frequencies accessible to Pulsar Timing Arrays (PTA). We currently lack firm empirical constraints on the amplitude of the GWB due to the dearth of confirmed SMBH binaries in the required mass range. Using HST/WFC3 images, we have discovered a z~0.2 quasar hosted in a merger remnant with two closely separated (0.13'' or ~430pc) continuum cores at the heart of the galaxy SDSSJ1010+1413. The two cores are spatially coincident with two powerful [OIII]-emitting point sources with quasar-like luminosities (L_AGN ~ 5x10^46 erg/s, suggesting the presence of a bound SMBH system, each with M_BH > 4x10^8 Msun. We place an upper limit on the merging timescale of the SMBH pair of 2.5 billion years, roughly the Universe lookback time at z~0.2. There is likely a population of quasar binaries similar to SDSSJ1010+1413 that contribute to a stochastic GWB that should be detected in the next several years. If the GWB is not detected this could indicate that SMBHs merge only over extremely long timescales, remaining as close separation binaries for many Hubble times, the so-called `final-parsec problem'.

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.

Cosmic voids, the less dense patches of the Universe, are promising laboratories to extract cosmological information. Thanks to their unique low density character, voids are extremely sensitive to diffuse components such as neutrinos and dark energy, and represent ideal environments to study modifications of gravity, where the effects of such modifications are expected to be more prominent. Robust void-related observables, including for example redshift-space distortions (RSD) and weak lensing around voids, are a promising way to chase and test new physics. Cosmological analysis of the large-scale structure of the Universe predominantly relies on the high density regions. Current and upcoming surveys are designed to optimize the extraction of cosmological information from these zones, but leave voids under-exploited. A dense, large area spectroscopic survey with imaging capabilities is ideal to exploit the power of voids fully. Besides helping illuminate the nature of dark energy, modified gravity, and neutrinos, this survey will give access to a detailed map of under-dense regions, providing an unprecedented opportunity to observe and study a so far under-explored galaxy population.

In recent years, machine learning (ML) methods have remarkably improved how cosmologists can interpret data. The next decade will bring new opportunities for data-driven cosmological discovery, but will also present new challenges for adopting ML methodologies and understanding the results. ML could transform our field, but this transformation will require the astronomy community to both foster and promote interdisciplinary research endeavors.

We measure the stellar populations as a function of radius for 90 early-type galaxies (ETGs) in the MASSIVE survey, a volume-limited integral-field spectroscopic (IFS) galaxy survey targeting all northern-sky ETGs with absolute K-band magnitude M_K < -25.3 mag, or stellar mass M* 4x10^11 M_sun, within 108 Mpc. We are able to measure reliable stellar population parameters for individual galaxies out to 10-20 kpc (1-3 R_e) depending on the galaxy. Focusing on ~R_e (~10 kpc), we find significant correlations between the abundance ratios, sigma, and M* at large radius, but we also find that the abundance ratios saturate in the highest-mass bin. We see a strong correlation between the kurtosis of the line of sight velocity distribution (h4) and the stellar population parameters beyond R_e. Galaxies with higher radial anisotropy appear to be older, with metal-poorer stars and enhanced [alpha/Fe]. We suggest that the higher radial anisotropy may derive from more accretion of small satellites. Finally, we see some evidence for correlations between environmental metrics (measured locally and on >5 Mpc scales) and the stellar populations, as expected if satellites are quenched earlier in denser environments.

We present optical long-slit spectroscopy and far-ultraviolet to near-infrared spectral energy distribution fitting of two diffuse dwarf galaxies, LSBG-285 and LSBG-750, which were recently discovered by the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). We measure redshifts using H$\alpha$ line emission, and find that these galaxies are at comoving distances of ${\approx}25$ and ${\approx}41$ Mpc, respectively, after correcting for the local velocity field. They have effective radii of $r_\mathrm{eff}=1.2$ and 1.8 kpc and stellar masses of $M_\star\approx2$-$3\times10^{7}~M_\odot$. There are no massive galaxies ($M_\star>10^{10} M_\odot$) within a comoving separation of at least 1.5 Mpc from LSBG-285 and 2 Mpc from LSBG-750. These sources are similar in size and surface brightness to ultra-diffuse galaxies, except they are isolated, star-forming objects that were optically selected in an environmentally blind survey. Both galaxies likely have low stellar metallicities $[Z_\star/Z_\odot] < -1.0$ and are consistent with the stellar mass-metallicity relation for dwarf galaxies. We set an upper limit on LSBG-750's rotational velocity of ${\sim}50$ km s$^{-1}$, which is comparable to dwarf galaxies of similar stellar mass with estimated halo masses $<10^{11}~M_\odot$. We find tentative evidence that the gas-phase metallicities in both of these diffuse systems are high for their stellar mass, though a statistically complete, optically-selected galaxy sample at very low surface brightness will be necessary to place these results into context with the higher-surface-brightness galaxy population.

Astronomical data often suffer from noise and incompleteness. We extend the common mixtures-of-Gaussians density estimation approach to account for situations with a known sample incompleteness by simultaneous imputation from the current model. The method, called GMMis, generalizes existing Expectation-Maximization techniques for truncated data to arbitrary truncation geometries and probabilistic rejection processes, as long as they can be specified and do not depend on the density itself. The method accounts for independent multivariate normal measurement errors for each of the observed samples and recovers an estimate of the error-free distribution from which both observed and unobserved samples are drawn. It can perform a separation of a mixtures-of-Gaussian signal from a specified background distribution whose amplitude may be unknown. We compare GMMis to the standard Gaussian mixture model for simple test cases with different types of incompleteness, and apply it to observational data from the NASA Chandra X-ray telescope. The python code is released as an open-source package at https://github.com/pmelchior/pyGMMis

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.

We use spatially resolved two-dimensional stellar velocity maps over a $107"\times 107"$ field of view to investigate the kinematic features of 90 early-type galaxies above stellar mass $10^{11.5}M_\odot$ in the MASSIVE survey. We measure the misalignment angle $\Psi$ between the kinematic and photometric axes and identify local features such as velocity twists and kinematically distinct components. We find 46% of the sample to be well aligned ($\Psi < 15^{\circ}$), 33% misaligned, and 21% without detectable rotation (non-rotators). Only 24% of the sample are fast rotators, the majority of which (91%) are aligned, whereas 57% of the slow rotators are misaligned with a nearly flat distribution of $\Psi$ from $15^{\circ}$ to $90^{\circ}$. 11 galaxies have $\Psi \gtrsim 60^{\circ}$ and thus exhibit minor-axis ("prolate") rotation in which the rotation is preferentially around the photometric major axis. Kinematic misalignments occur more frequently for lower galaxy spin or denser galaxy environments. Using the observed misalignment and ellipticity distributions, we infer the intrinsic shape distribution of our sample and find that MASSIVE slow rotators are consistent with being mildly triaxial, with mean axis ratios of $b/a=0.88$ and $c/a=0.65$. In terms of local kinematic features, 51% of the sample exhibit kinematic twists of larger than $20^{\circ}$, and 2 galaxies have kinematically distinct components. The frequency of misalignment and the broad distribution of $\Psi$ reported here suggest that the most massive early-type galaxies are mildly triaxial, and that formation processes resulting in kinematically misaligned slow rotators such as gas-poor mergers occur frequently in this mass range.

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 present 1,201 galaxies at $0.05<z<0.45$ that host tidal features, detected from the first $\sim\! 200$ deg$^2$ of imaging from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). All galaxies in the present sample have spectroscopic observations from the Sloan Digital Sky Survey (SDSS) spectroscopic campaigns, generating a sample of 21208 galaxies. Of these galaxies, we identify 214 shell systems and 987 stream systems. For 575 of these systems, we are additionally able to measure the $(g-i)$ colors of the tidal features. We find evidence for star formation in a subset of the streams, with the exception of streams around massive ellipticals, and find that stream host galaxies span the full range of stellar masses in our sample. Galaxies which host shells are predominantly red and massive: we find that observable shells form more frequently around ellipticals than around disc galaxies of the same stellar mass. Although the majority of the shells in our sample are consistent with being formed by minor mergers, $15\% \pm 4.4\%$ of shell host galaxies have $(g-i)$ colors as red as their host galaxy, consistent with being formed by major mergers. These "red shells" are additionally preferentially aligned with the major axis of the host galaxy, as previously predicted from simulations. We suggest that although the bulk of the observable shell population originates from fairly minor mergers, which preferentially form shells that are not aligned with the major axis of the galaxy, major mergers produce a significant number of observable shells.

We present a catalog of extended low-surface-brightness galaxies (LSBGs) identified in the Wide layer of the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). Using the first ${\sim}$200 deg$^2$ of the survey, we have uncovered 781 LSBGs, spanning red ($g-i\geq0.64$) and blue ($g-i<0.64$) colors and a wide range of morphologies. Since we focus on extended galaxies ($r_\mathrm{eff}=2.5$-$14^{\prime\prime}$), our sample is likely dominated by low-redshift objects. We define LSBGs to have mean surface brightnesses $\bar{\mu}_\mathrm{eff}(g)>24.3$ mag arcsec$^{-2}$, which allows nucleated galaxies into our sample. As a result, the central surface brightness distribution spans a wide range of $\mu_0(g)=18$-$27.4$ mag arcsec$^{-2}$, with 50% and 95% of galaxies fainter than 24.3 and 22 mag arcsec$^{-2}$, respectively. Furthermore, the surface brightness distribution is a strong function of color, with the red distribution being much broader and generally fainter than that of the blue LSBGs, and this trend shows a clear correlation with galaxy morphology. Red LSBGs typically have smooth light profiles that are well-characterized by single-component S\'{e}rsic functions. In contrast, blue LSBGs tend to have irregular morphologies and show evidence for ongoing star formation. We crossmatch our sample with existing optical, HI, and ultraviolet catalogs to gain insight into the physical nature of the LSBGs. We find that our sample is diverse, ranging from dwarf spheroidals and ultra-diffuse galaxies in nearby groups to gas-rich irregulars to giant LSB spirals, demonstrating the potential of the HSC-SSP to provide a truly unprecedented view of the LSBG population.

Narrow-line regions excited by active galactic nuclei (AGN) are important for studying AGN photoionization and feedback. Their strong [O III] lines can be detected with broadband images, allowing morphological studies of these systems with large-area imaging surveys. We develop a new technique to reconstruct the [O III] images using the Subaru Hyper Suprime-Cam (HSC) Survey aided with spectra from the Sloan Digital Sky Survey (SDSS). The technique involves a careful subtraction of the galactic continuum to isolate emission from the [O III]$\lambda$5007 and [O III]$\lambda$4959 lines. Compared to traditional targeted observations, this technique is more efficient at covering larger samples with less dedicated observational resources. We apply this technique to an SDSS spectroscopically selected sample of 300 obscured AGN at redshifts 0.1 - 0.7, uncovering extended emission-line region candidates with sizes up to tens of kpc. With the largest sample of uniformly derived narrow-line region sizes, we revisit the narrow-line region size-luminosity relation. The area and radii of the [O III] emission-line regions are strongly correlated with the AGN luminosity inferred from the mid-infrared (15 $\mu$m rest-frame) with a power-law slope of $0.62^{+0.05}_{-0.06} \pm 0.10$ (statistical and systemic errors), consistent with previous spectroscopic findings. We discuss the implications for the physics of AGN emission-line region and future applications of this technique, which should be useful for current and next-generation imaging surveys to study AGN photoionization and feedback with large statistical samples.

Hyper Suprime-Cam (HSC) is a wide-field imaging camera on the prime focus of the 8.2m Subaru telescope on the summit of Maunakea in Hawaii. A team of scientists from Japan, Taiwan and Princeton University is using HSC to carry out a 300-night multi-band imaging survey of the high-latitude sky. The survey includes three layers: the Wide layer will cover 1400 deg$^2$ in five broad bands ($grizy$), with a $5\,\sigma$ point-source depth of $r \approx 26$. The Deep layer covers a total of 26~deg$^2$ in four fields, going roughly a magnitude fainter, while the UltraDeep layer goes almost a magnitude fainter still in two pointings of HSC (a total of 3.5 deg$^2$). Here we describe the instrument, the science goals of the survey, and the survey strategy and data processing. This paper serves as an introduction to a special issue of the Publications of the Astronomical Society of Japan, which includes a large number of technical and scientific papers describing results from the early phases of this survey.

Quasars may have played a key role in limiting the stellar mass of massive galaxies. Identifying those quasars in the process of removing star formation fuel from their hosts is an exciting ongoing challenge in extragalactic astronomy. In this paper we present X-ray observations of eleven extremely red quasars (ERQs) with $L_{\rm bol}\sim 10^{47}$ erg s$^{-1}$ at $z=1.5-3.2$ with evidence for high-velocity ($v > 1000$ km s$^{-1}$) [OIII]$\lambda$5007\AA\ outflows. X-rays allow us to directly probe circumnuclear obscuration and to measure the instantaneous accretion luminosity. We detect ten out of eleven extremely red quasars available in targeted and archival data. Using a combination of X-ray spectral fitting and hardness ratios, we find that all of the ERQs show signs of absorption in the X-rays with inferred column densities of $N_{\rm H}\approx 10^{23}$ cm$^{-2}$, including four Compton-thick candidates ($N_{\rm H} > 10^{24}$ cm$^{-2}$). We stack the X-ray emission of the seven weakly detected sources, measuring an average column density of $N_{\rm H}\sim 8\times 10^{23}$ cm$^{-2}$. The absorption-corrected (intrinsic) $2-10$ keV X-ray luminosity of the stack is $2.7\times 10^{45}$ erg s$^{-1}$, consistent with X-ray luminosities of type 1 quasars of the same infrared luminosity. Thus, we find that ERQs are a highly obscured, borderline Compton-thick population, and based on optical and infrared data we suggest that these objects are partially hidden by their own equatorial outflows. However, unlike some quasars with known outflows, ERQs do not appear to be intrinsically underluminous in X-rays for their bolometric luminosity. Our observations indicate that low X-rays are not necessary to enable some types of radiatively driven winds.

We discuss the spectral analysis of a sample of 63 Active Galactic Nuclei (AGN) detected above a limiting flux of $S(8-24$ keV$)=7\times10^{-14}$ erg/s/cm$^2$ in the multi-tiered NuSTAR Extragalactic Survey program. The sources span a redshift range z=0-2.1 (median $\langle$z$\rangle=$0.58). The spectral analysis is performed over the broad 0.5-24 keV energy range, combining NuSTAR with Chandra and/or XMM-Newton data and employing empirical and physically motivated models. This constitutes the largest sample of AGN selected at $ >10$ keV to be homogeneously spectrally analyzed at these flux levels. We study the distribution of spectral parameters such as photon index, column density ($N_{\rm H}$), reflection parameter ($R$) and 10-40 keV luminosity ($L_{X}$). Heavily obscured ($\log [N_{\rm H}/cm^{-2}]\ge23$) and Compton Thick (CT; $\log [N_{\rm H}/cm^{-2}]\ge24$) AGN constitute $\sim$25% (15-17 sources) and $\sim$2-3% ( 1-2 sources) of the sample, respectively. The observed $N_{\rm H}$ distribution fairly agrees with predictions of Cosmic X-ray Background population synthesis models (CXBPSM). We estimate the intrinsic fraction of AGN as a function of $N_{H}$, accounting for the bias against obscured AGN in a flux-selected sample. The fraction of CT AGN relative to $\log [N_{\rm H}/cm^{-2}]=20-24$ AGN is poorly constrainted, formally in the range 2-56% (90% upper limit of 66%). We derived a fraction ($f_{abs}$) of obscured AGN ($\log [N_{\rm H}/cm^{-2}]=22-24$) as a function of $L_{X}$ in agreement with CXBPSM and previous z < 1 X-ray determinations. Furthermore $f_{abs}$ at z=0.1-0.5 and $\log (L_{X}/erg$ $s^{-1})\approx43.6-44.3$ agrees with observational measurements/trends obtained over larger redshift intervals. We report a significant anti-correlation of $R$ with $L_{X}$ (confirmed by our companion paper on stacked spectra) with considerable scatter around the median $R$ values.

Collisions and interactions between gas-rich galaxies are thought to be pivotal stages in their formation and evolution, causing the rapid production of new stars, and possibly serving as a mechanism for fueling supermassive black holes (BH). Harnessing the exquisite spatial resolution (~0.5 arcsec) afforded by the first ~170 deg^2 of the Hyper Suprime-Cam (HSC) Survey, we present our new constraints on the importance of galaxy-galaxy major mergers (1:4) in growing BHs throughout the last ~8 Gyrs. Utilizing mid-infrared observations in the WISE All-Sky survey, we robustly select active galactic nuclei (AGN) and mass-matched control galaxy samples, totaling ~140,000 spectroscopically confirmed systems at i<22 mag. We identify galaxy interaction signatures using a novel machine-learning random forest decision tree technique allowing us to select statistically significant samples of major-mergers, minor-mergers/irregular-systems, and non-interacting galaxies. We use these samples to show that galaxies undergoing mergers are a factor ~2-7 more likely to contain luminous obscured AGN than non-interacting galaxies, and this is independent of both stellar mass and redshift to z < 0.9. Furthermore, based on our comparison of AGN fractions in mass-matched samples, we determine that the most luminous AGN population (L_AGN > 10^45 erg/s) systematically reside in merging systems over non-interacting galaxies. Our findings show that galaxy-galaxy interactions do, on average, trigger luminous AGN activity substantially more often than in secularly evolving non-interacting galaxies, and we further suggest that the BH growth rate may be closely tied to the dynamical time of the merger system.

We investigate the relationship between the black hole accretion rate (BHAR) and star-formation rate (SFR) for Milky Way (MW) and Andromeda (M31)-mass progenitors from z = 0.2 - 2.5. We source galaxies from the Ks-band selected ZFOURGE survey, which includes multi-wavelenth data spanning 0.3 - 160um. We use decomposition software to split the observed SEDs of our galaxies into their active galactic nuclei (AGN) and star-forming components, which allows us to estimate BHARs and SFRs from the infrared (IR). We perform tests to check the robustness of these estimates, including a comparison to BHARs and SFRs derived from X-ray stacking and far-IR analysis, respectively. We find as the progenit- ors evolve, their relative black hole-galaxy growth (i.e. their BHAR/SFR ratio) increases from low to high redshift. The MW-mass progenitors exhibit a log-log slope of 0.64 +/- 0.11, while the M31-mass progenitors are 0.39 +/- 0.08. This result contrasts with previous studies that find an almost flat slope when adopting X-ray/AGN-selected or mass-limited samples and is likely due to their use of a broad mixture of galaxies with different evolutionary histories. Our use of progenitor-matched samples highlights the potential importance of carefully selecting progenitors when searching for evolutionary relationships between BHAR/SFRs. Additionally, our finding that BHAR/SFR ratios do not track the rate at which progenitors quench casts doubts over the idea that the suppression of star-formation is predominantly driven by luminous AGN feedback (i.e. high BHARs).

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.

The relation between X-ray luminosity (L_X) and ambient gas temperature (T) among massive galactic systems is an important cornerstone of both observational cosmology and galaxy-evolution modeling. In the most massive galaxy clusters, the relation is determined primarily by cosmological structure formation. In less massive systems, it primarily reflects the feedback response to radiative cooling of circumgalactic gas. Here we present a simple but powerful model for the L_X-T relation as a function of physical aperture R within which those measurements are made. The model is based on the precipitation framework for AGN feedback and assumes that the circumgalactic medium is precipitation-regulated at small radii and limited by cosmological structure formation at large radii. We compare this model with many different data sets and show that it successfully reproduces the slope and upper envelope of the L_X-T-R relation over the temperature range from ~0.2 keV through >10 keV. Our findings strongly suggest that the feedback mechanisms responsible for regulating star formation in individual massive galaxies have much in common with the precipitation-triggered feedback that appears to regulate galaxy-cluster cores.

The Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) is a three-layered imaging survey aimed at addressing some of the most outstanding questions in astronomy today, including the nature of dark matter and dark energy. The survey has been awarded 300 nights of observing time at the Subaru Telescope and it started in March 2014. This paper presents the first public data release of HSC-SSP. This release includes data taken in the first 1.7 years of observations (61.5 nights) and each of the Wide, Deep, and UltraDeep layers covers about 108, 26, and 4 square degrees down to depths of i~26.4, ~26.5, and ~27.0 mag, respectively (5sigma for point sources). All the layers are observed in five broad bands (grizy), and the Deep and UltraDeep layers are observed in narrow bands as well. We achieve an impressive image quality of 0.6 arcsec in the i-band in the Wide layer. We show that we achieve 1-2 per cent PSF photometry (rms) both internally and externally (against Pan-STARRS1), and ~10 mas and 40 mas internal and external astrometric accuracy, respectively. Both the calibrated images and catalogs are made available to the community through dedicated user interfaces and database servers. In addition to the pipeline products, we also provide value-added products such as photometric redshifts and a collection of public spectroscopic redshifts. Detailed descriptions of all the data can be found online. The data release website is https://hsc-release.mtk.nao.ac.jp/.

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 report the discovery of a diffuse stellar cloud with an angular extent $\gtrsim30^{\prime\prime}$, which we term "Sumo Puff", in data from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). While we do not have a redshift for this object, it is in close angular proximity to a post-merger galaxy at redshift $z=0.0431$ and is projected within a few virial radii (assuming similar redshifts) of two other ${\sim}L_\star$ galaxies, which we use to bracket a potential redshift range of $0.0055 < z < 0.0431$. The object's light distribution is flat, as characterized by a low Sersic index ($n\sim0.3$). It has a low central $g$-band surface brightness of ${\sim}26.4$ mag arcsec$^{-2}$, large effective radius of ${\sim}13^{\prime\prime}$ (${\sim}11$ kpc at $z=0.0431$ and ${\sim}1.5$ kpc at $z=0.0055$), and an elongated morphology ($b/a\sim0.4$). Its red color ($g-i\sim1$) is consistent with a passively evolving stellar population and similar to the nearby post-merger galaxy, and we may see tidal material connecting Sumo Puff with this galaxy. We offer two possible interpretations for the nature of this object: (1) it is an extreme, galaxy-size tidal feature associated with a recent merger event, or (2) it is a foreground dwarf galaxy with properties consistent with a quenched, disturbed ultra-diffuse galaxy. We present a qualitative comparison with simulations that demonstrates the feasibility of forming a structure similar to this object in a merger event. Follow-up spectroscopy and/or deeper imaging to confirm the presence of the bridge of tidal material will be necessary to reveal the true nature of this object.

We present the first systematic investigation of the existence, spatial distribution, and kinematics of warm ionized gas as traced by the [O II] 3727AA emission line in 74 of the most massive galaxies in the local Universe. All of our galaxies have deep integral field spectroscopy from the volume- and magnitude-limited MASSIVE survey of early-type galaxies with stellar mass log(M_*/M_sun) > 11.5 (M_K < -25.3 mag) and distance D < 108 Mpc. Of the 74 galaxies in our sample, we detect warm ionized gas in 28, which yields a global detection fraction of 38\pm6% down to a typical [O II] equivalent width limit of 2AA. MASSIVE fast rotators are more likely to have gas than MASSIVE slow rotators with detection fractions of 80\pm10% and 28\pm6%, respectively. The spatial extents span a wide range of radii (0.6 - 18.2 kpc; 0.1 - 4R_e), and the gas morphologies are diverse, with 17/28 = 61\pm9% being centrally concentrated, 8/28 = 29\pm9% exhibiting clear rotation out to several kpc, and 3/28 = 11\pm6% being extended but patchy. Three out of four fast rotators show kinematic alignment between the stars and gas, whereas the two slow rotators with robust kinematic measurements available exhibit kinematic misalignment. Our inferred warm ionized gas masses are roughly ~10^5M_sun. The emission line ratios and radial equivalent width profiles are generally consistent with excitation of the gas by the old underlying stellar population. We explore different gas origin scenarios for MASSIVE galaxies and find that a variety of physical processes are likely at play, including internal gas recycling, cooling out of the hot gaseous halo, and gas acquired via mergers.

To investigate the relationship between thermal and non-thermal components in merger galaxy clusters, we present deep JVLA and Chandra observations of the HST Frontier Fields cluster MACS J0717.5+3745. The Chandra image shows a complex merger event, with at least four components belonging to different merging subclusters. NW of the cluster, $\sim 0.7$ Mpc from the center, there is a ram-pressure-stripped core that appears to have traversed the densest parts of the cluster after entering the ICM from the direction of a galaxy filament to the SE. We detect a density discontinuity NNE of this core which we speculate is associated with a cold front. Our radio images reveal new details for the complex radio relic and radio halo in this cluster. In addition, we discover several new filamentary radio sources with sizes of 100-300 kpc. A few of these seem to be connected to the main radio relic, while others are either embedded within the radio halo or projected onto it. A narrow-angled-tailed (NAT) radio galaxy, a cluster member, is located at the center of the radio relic. The steep spectrum tails of this AGN leads into the large radio relic where the radio spectrum flattens again. This morphological connection between the NAT radio galaxy and relic provides evidence for re-acceleration (revival) of fossil electrons. The presence of hot $\gtrsim 20$ keV ICM gas detected by Chandra near the relic location provides additional support for this re-acceleration scenario.

NGC 1448 is one of the nearest luminous galaxies ($L_{8-1000\mu m} >$ 10$^{9} L_{\odot}$) to ours ($z$ $=$ 0.00390), and yet the active galactic nucleus (AGN) it hosts was only recently discovered, in 2009. In this paper, we present an analysis of the nuclear source across three wavebands: mid-infrared (MIR) continuum, optical, and X-rays. We observed the source with the Nuclear Spectroscopic Telescope Array (NuSTAR), and combined this data with archival Chandra data to perform broadband X-ray spectral fitting ($\approx$0.5-40 keV) of the AGN for the first time. Our X-ray spectral analysis reveals that the AGN is buried under a Compton-thick (CT) column of obscuring gas along our line-of-sight, with a column density of $N_{\rm H}$(los) $\gtrsim$ 2.5 $\times$ 10$^{24}$ cm$^{-2}$. The best-fitting torus models measured an intrinsic 2-10 keV luminosity of $L_{2-10\rm{,int}}$ $=$ (3.5-7.6) $\times$ 10$^{40}$ erg s$^{-1}$, making NGC 1448 one of the lowest luminosity CTAGNs known. In addition to the NuSTAR observation, we also performed optical spectroscopy for the nucleus in this edge-on galaxy using the European Southern Observatory New Technology Telescope. We re-classify the optical nuclear spectrum as a Seyfert on the basis of the Baldwin-Philips-Terlevich diagnostic diagrams, thus identifying the AGN at optical wavelengths for the first time. We also present high spatial resolution MIR observations of NGC 1448 with Gemini/T-ReCS, in which a compact nucleus is clearly detected. The absorption-corrected 2-10 keV luminosity measured from our X-ray spectral analysis agrees with that predicted from the optical [OIII]$\lambda$5007\AA\ emission line and the MIR 12$\mu$m continuum, further supporting the CT nature of the AGN.

We use new precision measurements of black hole masses from water megamaser disks to investigate scaling relations between macroscopic galaxy properties and supermassive black hole (BH) mass. The megamaser-derived BH masses span 10^6-10^8 M_sun, while all the galaxy properties that we examine (including stellar mass, central mass density, central velocity dispersion) lie within a narrow range. Thus, no galaxy property correlates tightly with M_BH in ~L* spiral galaxies. Of them all, stellar velocity dispersion provides the tightest relation, but at fixed sigma* the mean megamaser M_BH are offset by -0.6+/-0.1 dex relative to early-type galaxies. Spiral galaxies with non-maser dynamical BH masses do not show this offset. At low mass, we do not yet know the full distribution of BH mass at fixed galaxy property; the non-maser dynamical measurements may miss the low-mass end of the BH distribution due to inability to resolve the spheres of influence and/or megamasers may preferentially occur in lower-mass BHs.

Studies of the physical properties of local elliptical galaxies (e.g., gas temperatures, halo masses, stellar kinematics) are shedding new light on galaxy formation. Here we present the hot X-ray gas properties of 33 early-type systems within the MASSIVE galaxy survey sample that have archival Chandra X-ray observations. Through careful X-ray spectral modeling, we derive X-ray luminosities (L_X) and plasma temperatures (T_gas) for the diffuse gas components in these galaxies. We combine the MASSIVE sample with 41 galaxies from the ATLAS^3D survey to investigate the X-ray and optical properties of a statistically significant sample of nearby early-type galaxies across a wide-range of environments. We deduce that all early-type galaxies (independent of galaxy mass and rotational support) follow a universal scaling law such that L_X~T_gas^4.5. When X-ray measurements are performed consistently in apertures set by the galaxy stellar content, the wide-scale environment does not contribute to the intrinsic scatter (~0.5dex) within the scaling relation. We further demonstrate that the scatter in L_X around both K-band luminosity (L_K) and the galaxy stellar velocity dispersion is primarily driven by T_gas, with no clear trends with halo mass, radio power, or angular momentum of the stars. It is not trivial to tie the origin of the gas directly to either the stellar mass or the galaxy potential. Indeed, our data require a steeper relation between L_X, L_K, and sigma_e than predicted by standard mass-loss models. Finally, we find a statistically significant correlation between sigma_e and T_gas, suggesting that T_gas is set by the galaxy potential inside the optical effective radius. We conclude that within the inner-most 10-30kpc region, early-types maintain pressure-supported hot gas, with a minimum T_gas set by the virial temperature, but the majority show evidence for some additional heating.

We report on high-resolution JVLA and Chandra observations of the HST Frontier Cluster MACS J0717.5+3745. MACS J0717.5+3745 offers the largest contiguous magnified area of any known cluster, making it a promising target to search for lensed radio and X-ray sources. With the high-resolution 1.0-6.5 GHz JVLA imaging in A and B configuration, we detect a total of 51 compact radio sources within the area covered by the HST imaging. Within this sample we find 7 lensed sources with amplification factors larger than $2$. None of these sources are identified as multiply-lensed. Based on the radio luminosities, the majority of these sources are likely star forming galaxies with star formation rates of 10-50 M$_\odot$ yr$^{-1}$ located at $1 \lesssim z \lesssim 2$. Two of the lensed radio sources are also detected in the Chandra image of the cluster. These two sources are likely AGN, given their $2-10$ keV X-ray luminosities of $\sim 10^{43-44}$ erg s$^{-1}$. From the derived radio luminosity function, we find evidence for an increase in the number density of radio sources at $0.6<z<2.0$, compared to a $z < 0.3$ sample. Our observations indicate that deep radio imaging of lensing clusters can be used to study star forming galaxies, with star formation rates as low as $\sim10$ M$_{\odot}$ yr$^{-1}$, at the peak of cosmic star formation history.

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.

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

The HST Frontier Fields cluster MACS J1149.6+2223 is one of the most complex merging clusters, believed to consist of four dark matter halos. We present results from deep (365 ks) Chandra observations of the cluster, which reveal the most distant cold front (z=0.544) discovered to date. In the cluster outskirts, we also detect hints of a surface brightness edge that could be the bow shock preceding the cold front. The substructure analysis of the cluster identified several components with large relative radial velocities, thus indicating that at least some collisions occur almost along the line of sight. The inclination of the mergers with respect to the plane of the sky poses significant observational challenges at X-ray wavelengths. MACS J1149.6+2223 possibly hosts a steep-spectrum radio halo. If the steepness of the radio halo is confirmed, then the radio spectrum, combined with the relatively regular ICM morphology, could indicate that MACS J1149.6+2223 is an old merging cluster.

We present a catalogue containing the redshifts of 3,660 X-ray selected targets in the XXL southern field. The redshifts were obtained with the AAOmega spectrograph and 2dF fibre positioner on the Anglo-Australian Telescope. The catalogue contains 1,515 broad line AGN, 528 stars, and redshifts for 41 out of the 49 brightest X-ray selected clusters in the XXL southern field.

Merging galaxy clusters leave long-lasting signatures on the baryonic and non-baryonic cluster constituents, including shock fronts, cold fronts, X-ray substructure, radio halos, and offsets between the dark matter and the gas components. Using observations from Chandra, the Jansky Very Large Array, the Giant Metrewave Radio Telescope, and the Hubble Space Telescope, we present a multiwavelength analysis of the merging Frontier Fields cluster MACS J0416.1-2403 (z=0.396), which consists of a NE and a SW subclusters whose cores are separated on the sky by ~250 kpc. We find that the NE subcluster has a compact core and hosts an X-ray cavity, yet it is not a cool core. Approximately 450 kpc south-south west of the SW subcluster, we detect a density discontinuity that corresponds to a compression factor of ~1.5. The discontinuity was most likely caused by the interaction of the SW subcluster with a less massive structure detected in the lensing maps SW of the subcluster's center. For both the NE and the SW subclusters, the dark matter and the gas components are well-aligned, suggesting that MACS J0416.1-2403 is a pre-merging system. The cluster also hosts a radio halo, which is unusual for a pre-merging system. The halo has a 1.4 GHz power of (1.06 +/- 0.09) x 10^{24} W Hz^{-1}, which is somewhat lower than expected based on the X-ray luminosity of the cluster. We suggest that we are either witnessing the birth of a radio halo, or have discovered a rare ultra-steep spectrum halo.

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.

Observational studies of nearby galaxies have demonstrated correlations between the mass of the central supermassive black holes (BHs) and properties of the host galaxies, notably the stellar bulge mass or central stellar velocity dispersion. Motivated by these correlations, the theoretical paradigm has emerged, in which BHs and bulges co-evolve. However, this picture was challenged by observational and theoretical studies, which hinted that the fundamental connection may be between BHs and dark matter halos, and not necessarily with their host galaxies. Based on a study of 3130 elliptical galaxies $-$ selected from the Sloan Digital and ROSAT All Sky Surveys $-$ we demonstrate that the central stellar velocity dispersion exhibits a significantly tighter correlation with the total gravitating mass, traced by the X-ray luminosity of the hot gas, than with the stellar mass. This hints that the central stellar velocity dispersion, and hence the central gravitational potential, may be the fundamental property of elliptical galaxies that is most tightly connected to the larger-scale dark matter halo. Furthermore, using the central stellar velocity dispersion as a surrogate for the BH mass, we find that in elliptical galaxies the inferred BH mass and inferred total gravitating mass within the virial radius (or within five effective radii) can be expressed as $M_{\rm{BH}} \propto M_{\rm tot}^{1.6^{+0.6}_{-0.4}} $ (or $M_{\rm{BH}} \propto M_{\rm{5r_{eff}}}^{1.8^{+0.7}_{-0.6}}$). These results are consistent with a picture in which the BH mass is directly set by the central stellar velocity dispersion, which, in turn, is determined by the total gravitating mass of the system.

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.

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

We present the X-ray and optical properties of the galaxy groups selected in the Chandra X-Bo\"otes survey. Our final sample comprises 32 systems at \textbf{$z<1.75$}, with 14 below $z = 0.35$. For these 14 systems we estimate velocity dispersions ($\sigma_{gr}$) and perform a virial analysis to obtain the radii ($R_{200}$ and $R_{500}$) and total masses ($M_{200}$ and $M_{500}$) for groups with at least five galaxy members. We use the Chandra X-ray observations to derive the X-ray luminosity ($L_X$). We examine the performance of the group properties $\sigma_{gr}$, $L_{opt}$ and $L_X$, as proxies for the group mass. Understanding how well these observables measure the total mass is important to estimate how precisely the cluster/group mass function is determined. Exploring the scaling relations built with the X-Bo\"otes sample and comparing these with samples from the literature, we find a break in the $L_X$-$M_{500}$ relation at approximately $M_{500} = 5\times10^{13}$ M$_\odot$ (for $M_{500} > 5\times10^{13}$ M$_\odot$, $M_{500} \propto L_X^{0.61\pm0.02}$, while for $M_{500} \leq 5\times10^{13}$ M$_\odot$, $M_{500} \propto L_X^{0.44\pm0.05}$). Thus, the mass-luminosity relation for galaxy groups cannot be described by the same power law as galaxy clusters. A possible explanation for this break is the dynamical friction, tidal interactions and projection effects which reduce the velocity dispersion values of the galaxy groups. By extending the cluster luminosity function to the group regime, we predict the number of groups that new X-ray surveys, particularly eROSITA, will detect. Based on our cluster/group luminosity function estimates, eROSITA will identify $\sim$1800 groups ($L_X = 10^{41}-10^{43}$ ergs s$^{-1}$) within a distance of 200 Mpc. Since groups lie in large scale filaments, this group sample will map the large scale structure of the local universe.

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

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

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.

In the past few years more and more pieces of evidence have been presented for a revision of the widely accepted Unified Model of Active Galactic Nuclei. A model based solely on orientation cannot explain all the observed phenomenology. In the following, we will present evidence that accretion rate is also a key parameter for the presence of Hidden Broad Line Regions in Seyfert 2 galaxies. Our sample consists of 21 sources with polarized Hidden Broad Lines and 18 sources without Hidden Broad Lines. We use stellar velocity dispersions from several studies on the CaII and Mg b triplets in Seyfert 2 galaxies, to estimate the mass of the central black holes via the Mbh-{\sigma}\ast relation. The ratio between the bolometric luminosity, derived from the intrinsic (i.e. unabsorbed) X-ray luminosity, and the Eddington luminosity is a measure of the rate at which matter accretes onto the central supermassive black hole. A separation between Compton-thin HBLR and non-HBLR sources is clear, both in accretion rate (log Lbol/LEdd = -1.9) and in luminosity (log Lbol = 43.90). When, properly luminosity-corrected, Compton-thick sources are included, the separation between HBLR and non-HBLR is less sharp but no HBLR source falls below the Eddington ratio threshold. We speculate that non-HBLR Compton-thick sources with accretion rate higher than the threshold, do possess a BLR, but something, probably related to their heavy absorption, is preventing us from observing it even in polarized light. Our results for Compton-thin sources support theoretical expectations. In a model presented by Nicastro (2000), the presence of broad emission lines is intrinsically connected with disk instabilities occuring in proximity of a transition radius, which is a function of the accretion rate, becoming smaller than the innermost stable orbit for very low accretion rates and therefore luminosities.

Surveys have revealed a class of object displaying both high X-ray luminosities (Lx > 10^42 erg/s), and a lack of a discernible active galactic nucleus (AGN) in the optical band. If these sources are powered by star formation activity alone, they would be the most extreme X-ray luminosity star forming galaxies known. We have investigated the mechanism driving the X-ray luminosities of such galaxies by studying the X-ray emission of three moderate redshift (z ~ 0.1) examples of this class, selected from a cross-correlation of the SDSS-DR5 and 2XMMp-DR0 catalogues. X-ray spatial and long-term variability diagnostics of these sources suggest that they are compact X-ray emitters. This result is supported by the detection of rapid short term variability in an observation of one of the sources. The X-ray spectra of all three sources are best fitted with a simple absorbed power-law model, thus betraying no significant signs of star formation. These results indicate that the X-ray emission is powered by AGN activity. But why do these sources not display optical AGN signatures? We show that the most likely explanation is that the optical AGN emission lines are being diluted by star formation signatures from within their host 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 preliminary results of two Gemini campaigns to constrain the mass of the black hole in an ultraluminous X-ray source (ULX) via optical spectroscopy. Pilot studies of the optical counterparts of a number of ULXs revealed two candidates for further detailed study, based on the presence of a broad He II 4686 Angstrom emission line. A sequence of 10 long-slit spectra were obtained for each object, and the velocity shift of the ULX counterpart measured. Although radial velocity variations are observed, they are not sinusoidal, and no mass function is obtained. However, the broad He II line is highly variable on timescales shorter than a day. If associated with the reprocessing of X-rays in the accretion disc, its breadth implies that the disc must be close to face-on.

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 new Chandra observations that complete a sample of seventeen (17) luminous infrared galaxies (LIRGs) with D < 60 Mpc and low Galactic column densities of N_H < 5 X 10^20 cm^-2. The LIRGs in our sample have total infrared (8-1000um) luminosities in the range of L_IR ~ (1-8) X 10^11 L_sol. The high-resolution imaging and X-ray spectral information from our Chandra observations allow us to measure separately X-ray contributions from active galactic nuclei (AGNs) and normal galaxy processes (e.g., X-ray binaries and hot gas). We utilized total infrared plus UV luminosities to estimate star-formation rates (SFRs) and K-band luminosities and optical colors to estimate stellar masses (M*) for the sample. Under the assumption that the galaxy-wide 2-10 keV luminosity (LX) traces the combined emission from high mass X-ray binaries (HMXBs) and low mass X-ray binaries (LMXBs), and that the power output from these components are linearly correlated with SFR and M*, respectively, we constrain the relation LX = alpha M* + beta SFR. To achieve this, we construct a Chandra-based data set composed of our new LIRG sample combined with additional samples of less actively star-forming normal galaxies and more powerful LIRGs and ultraluminous infrared galaxies (ULIRGs) from the literature. Using these data, we measure best-fit values of alpha = (9.05 +/- 0.37) X 10^28 ergs s^-1 Msol^-1 and beta = (1.62 +/- 0.22) X 10^39 ergs s^-1 (Msol yr^-1)^-1. This scaling provides a more physically meaningful estimate of LX, with ~0.1-0.2 dex less scatter, than a direct linear scaling with SFR (abridged).

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

(Abridged) We investigate the far-infrared properties of X-ray sources detected in the Chandra Deep Field-South (CDF-S) survey using the ultra-deep 70um and 24um Spitzer observations taken in this field. We rely on stacking analyses of the 70um data to characterise the average 70um properties of the X-ray sources. Using Spitzer-IRS data of the Swift-BAT sample of z~0 active galactic nuclei (hereafter, AGNs), we show that the 70um/24um flux ratio can distinguish between AGN-dominated and starburst-dominated systems out to z~1.5. From stacking analysis we find that both high redshift and z~0 AGNs follow the same tendency toward warmer 70um/24um colours with increasing X-ray luminosity (L_X). We also show that the 70um flux can be used to determine the infrared (8-1000um) luminosities of high redshift AGNs. We use this information to show that L_X=10^{42-43} erg/s AGNs at high redshifts (z=1-2) have infrared to X-ray luminosity ratios (hereafter, L_IR/L_X) that are, on average, 4.7_{-2.0}^{+10.2} and 12.7^{+7.1}_{-2.6} times higher than AGNs with similar X-ray luminosities at z=0.5-1 and z~0, respectively. By comparison, we find that the L_IR/L_X ratios of L_X=10^{43-44} erg/s AGNs remain largely unchanged across this same redshift interval. We explore the consequences that these results may have on the identification of distant, potentially Compton thick AGNs using L_IR/L_X ratios. We discuss possible scenarios for the observed increase in the L_IR/L_X ratio with redshift, including changes in the dust covering factor of AGNs and/or the star formation rates of their host galaxies. Finally, we show how deep observations to be undertaken by the Herschel Space Observatory will enable us to discriminate between these proposed scenarios and also identify Compton-thick AGNs at high redshifts.

Using Spitzer-IRS spectroscopy, we investigate the ubiquity of Active Galactic Nuclei (AGN) in a complete (~94%), volume-limited sample of the most bolometrically-luminous galaxies (L_IR > (0.3-20) x 10^10 L_sun) to D < 15 Mpc. Our analyses are based on the detection of the high-excitation emission line [NeV](\lambda 14.32 um; 97.1 eV) to unambiguously identify AGN activity. We find that 17 of the 64 IR-bright galaxies in our sample host AGN activity (~27^{+8}_{-6}%), >50% of which are not identified as AGNs using optical spectroscopy. The large AGN fraction indicates a tighter connection between AGN activity and IR luminosity for galaxies in the local Universe than previously found, potentially indicating a close association between AGN activity and star formation. The optically unidentified AGNs span a wide range of galaxy type (S0-Ir) and are typically starburst-dominated systems hosting modest-luminosity AGN activity (L_[NeV] ~ 10^37 - 10^39 erg s^-1). The non-identification of optical AGN signatures in the majority of these galaxies appears to be due to extinction towards the AGN, rather than intrinsically low-luminosity AGN activity. Examination of optical images shows that the optically unidentified AGNs with evidence for extinction are hosted in either highly inclined galaxies or galaxies with dust lanes, indicating that obscuration of the AGN is not necessarily due to an obscuring torus. We therefore conclude that optical spectroscopic surveys miss approximately half of the AGN population simply due to extinction through the host galaxy.