Mid-Infrared Colors of Dwarf Galaxies: Young Starbursts Mimicking Active Galactic Nuclei

Searching for active galactic nuclei (AGN) in dwarf galaxies is important for our understanding of the seed black holes that formed in the early Universe. Here, we test infrared selection methods for AGN activity at low galaxy masses. Our parent sample consists of ~18,000 nearby dwarf galaxies (M*< 3 x 10^9 Msun, $z<0.055$) in the Sloan Digital Sky Survey with significant detections in the first three bands of the AllWISE data release from the Wide-field Infrared Survey Explorer (WISE). First, we demonstrate that the majority of optically-selected AGNs in dwarf galaxies are not selected as AGNs using WISE infrared color diagnostics and that the infrared emission is dominated by the host galaxies. We then investigate the infrared properties of optically-selected star-forming dwarf galaxies, finding that the galaxies with the reddest infrared colors are the most compact, with blue optical colors, young stellar ages and large specific star formation rates. These results indicate that great care must be taken when selecting AGNs in dwarf galaxies using infrared colors, as star-forming dwarf galaxies are capable of heating dust in such a way that mimics the infrared colors of more luminous AGNs. In particular, a simple $\mathrm{W1}-\mathrm{W2}$ color cut alone should not be used to select AGNs in dwarf galaxies. With these complications in mind, we present a sample of 41 dwarf galaxies worthy of follow-up observations that fall in WISE infrared color space typically occupied by more luminous AGNs.Comment: 17 pages, 12 figures, accepted by the Astrophysical Journa

The Rest Frame Ultraviolet Spectra of UV-Selected Active Galactic Nuclei at z ~ 2-3

We present new results for a sample of 33 narrow-lined UV-selected active galactic nuclei (AGNs), identified in the course of a spectroscopic survey for star-forming galaxies at z ~ 2-3. The rest-frame UV composite spectrum for our AGN sample shows several emission lines characteristic of AGNs, as well as interstellar absorption features seen in star-forming Lyman Break Galaxies (LBGs). We report a detection of NIV]1486, which has been observed in high-redshift radio galaxies, as well as in rare optically-selected quasars. The UV continuum slope of the composite spectrum is significantly redder than that of a sample of non-AGN UV-selected star forming galaxies. Blueshifted SiIV absorption provides evidence for outflowing highly-ionized gas in these objects at speeds of ~ 10^(3) km/s, quantitatively different from what is seen in the outflows of non-AGN LBGs. Grouping the individual AGNs by parameters such as Ly-alpha equivalent width, redshift, and UV continuum magnitude allows for an analysis of the major spectroscopic trends within the sample. Stronger Ly-alpha emission is coupled with weaker low-ionization absorption, which is similar to what is seen in the non-AGN LBGs, and highlights the role that cool interstellar gas plays in the escape of Ly-alpha photons. However, the AGN composite does not show the same trends between Ly-alpha strength and extinction seen in the non-AGN LBGs. These results represent the first such comparison at high-redshift between star-forming galaxies and similar galaxies that host AGN activity.Comment: 13 pages, 6 figures, accepted by Ap

SALT Long-slit Spectroscopy of Luminous Obscured Quasars: An Upper Limit on the Size of the Narrow-Line Region?

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

Composite Spectral Energy Distributions and Infrared-Optical Colors of Type 1 and Type 2 Quasars

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

The Kinematics of Multiple-Peaked Lyα Emission in Star-Forming Galaxies at z ~ 2-3

We present new results on the Lyα emission-line kinematics of 18 z ~ 2-3 star-forming galaxies with multiple-peaked Lyα profiles. With our large spectroscopic database of UV-selected star-forming galaxies at these redshifts, we have determined that ~30% of such objects with detectable Lyα emission display multiple-peaked emission profiles. These profiles provide additional constraints on the escape of Lyα photons due to the rich velocity structure in the emergent line. Despite recent advances in modeling the escape of Lyα from star-forming galaxies at high redshifts, comparisons between models and data are often missing crucial observational information. Using Keck II NIRSPEC spectra of Hα (z ~ 2) and [O III]λ5007 (z ~ 3), we have measured accurate systemic redshifts, rest-frame optical nebular velocity dispersions, and emission-line fluxes for the objects in the sample. In addition, rest-frame UV luminosities and colors provide estimates of star formation rates and the degree of dust extinction. In concert with the profile sub-structure, these measurements provide critical constraints on the geometry and kinematics of interstellar gas in high-redshift galaxies. Accurate systemic redshifts allow us to translate the multiple-peaked Lyα profiles into velocity space, revealing that the majority (11/18) display double-peaked emission straddling the velocity-field zero point with stronger red-side emission. Interstellar absorption-line kinematics suggest the presence of large-scale outflows for the majority of objects in our sample, with an average measured interstellar absorption velocity offset of (Δv_(abs))=–230 km s^(–1). A comparison of the interstellar absorption kinematics for objects with multiple- and single-peaked Lyα profiles indicate that the multiple-peaked objects are characterized by significantly narrower absorption line widths. We compare our data with the predictions of simple models for outflowing and infalling gas distributions around high-redshift galaxies. While popular "shell" models provide a qualitative match with many of the observations of Lyα emission, we find that in detail there are important discrepancies between the models and data, as well as problems with applying the framework of an expanding thin shell of gas to explain high-redshift galaxy spectra. Our data highlight these inconsistencies, as well as illuminating critical elements for success in future models of outflow and infall in high-redshift galaxies

Gemini Long-Slit Observations of Luminous Obscured Quasars: Further Evidence for an Upper Limit on the Size of the Narrow-Line Region

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

A Spectroscopic Survey of Wise -Selected Obscured Quasars with the Southern African Large Telescope

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

Gemini Long-slit Observations of Luminous Obscured Quasars: Further Evidence for an Upper Limit on the Size of the Narrow-Line Region

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

A Tale of Two Narrow-Line Regions: Ionization, Kinematics, and Spectral Energy Distributions for a Local Pair of Merging Obscured Active Galaxies

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."Comment: 15 pages, 10 figures, accepted by the Astrophysical Journa

The Intrinsic Eddington Ratio Distribution of Active Galactic Nuclei in Star-Forming Galaxies from the Sloan Digital Sky Survey

An important question in extragalactic astronomy concerns the distribution of black hole accretion rates of active galactic nuclei (AGNs). 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 whether or not an intrinsic Eddington ratio distribution that takes the form of a Schechter function is consistent with previous work suggesting 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 AGNs 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 their 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 AGNs 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 X-rays