Source author record

Jenny E. Greene

Jenny E. Greene appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

119works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

119 published item(s)

preprint2026arXiv

The X-Ray Dot: Exotic Dust or a Late-Stage Little Red Dot?

JWST's "Little Red Dots" (LRDs) are increasingly interpreted as active galactic nuclei (AGN) obscured by dense thermalized gas rather than dust as evidenced by their X-ray weakness, blackbody-like continua, and Balmer line profiles. A key question is how LRDs connect to standard UV-luminous AGN and whether transitional phases exist and if they are observable. We present the "X-Ray Dot" (XRD), a compact source at $z=3.28$ observed by the NIRSpec WIDE GTO survey. The XRD exhibits LRD hallmarks: a blackbody-like ($T_{\rm eff} \simeq 6400\,$K) red continuum, a faint but blue rest-UV excess, falling mid-IR emission, and broad Balmer lines ($\rm FWHM \sim 2700-3200\,km\,s^{-1}$). Unlike LRDs, however, it is remarkably X-ray luminous ($L_\textrm{2$-$10$\,$keV} = 10^{44.18}\,$erg$\,$s$^{-1}$) and has a continuum inflection that is bluewards of the Balmer limit. We find that the red rest-optical and blue mid-IR continuum cannot be reproduced by standard dust-attenuated AGN models without invoking extremely steep extinction curves, nor can the weak mid-IR emission be reconciled with well-established X-ray--torus scaling relations. We therefore consider an alternative scenario: the XRD may be an LRD in transition, where the gas envelope dominates the optical continuum but optically thin sightlines allow X-rays to escape. The XRD may thus provide a physical link between LRDs and standard AGN, offering direct evidence that LRDs are powered by supermassive black holes and providing insight into their accretion properties.

preprint2025arXiv

A weak Ly$α$ halo for an extremely bright Little Red Dot: Indications of enshrouded SMBH growth

The abundant population of "Little Red Dots" (LRDs)-compact objects with red UV to optical colors and broad Balmer lines at high redshift-is unveiling new insights into the properties of early active galactic nuclei (AGN). Perhaps the most surprising features of this population are the presence of Balmer absorption and ubiquitous strong Balmer breaks. Recent models link these features to an active supermassive black hole (SMBH) cocooned in very dense gas ($N_{\rm H}\sim10^{24}\,\rm cm^{-2}$). We present a stringent test of such models using VLT/MUSE observations of A2744-45924, the most luminous LRD known to date ($L_{\rm Hα}\approx10^{44}~\rm erg\,s^{-1}$), located behind the Abell-2744 lensing cluster at $z=4.464$ ($μ=1.8$). We detect a moderately extended Ly$α$ nebula ($h\approx5.7$ pkpc), spatially offset from the point-like H$α$ seen by JWST. The Ly$α$ emission is narrow ($\rm FWHM=270\pm 15~km\,s^{-1}$), spatially offset to H$α$, and faint ($\rm Lyα=0.07Hα$) compared to Ly$α$ nebulae typically observed around quasars of similar luminosity. We detect compact N$\,$IV]$λ$1486 emission, spatially aligned with H$α$, and a spatial shift in the far-UV continuum matching the Ly$α$ offset. We discuss that H$α$ and Ly$α$ have distinct physical origins: H$α$ originates from the AGN, while Ly$α$ is powered by star formation. In the environment of A2744-45924, we identify four extended Ly$α$ halos ($Δz<0.02$, $Δr<100$ pkpc). Their Ly$α$ luminosities match expectations based on H$α$ emission, indicating no evidence for radiation from A2744-45924 affecting its surroundings. The lack of strong, compact, and broad Ly$α$ and the absence of a luminous extended halo, suggest that the UV AGN light is obscured by dense gas cloaking the SMBH with covering factor close to unity.

preprint2025arXiv

Picture an Astronomer: Best Practices for Retaining Talent in Astrophysics

Women are consistently underrepresented in astrophysics yet are simultaneously subject to disproportionate attrition at every career stage. This disparity between demonstrated efficacy in job performance and ultimate career outcome was the primary motivation for the Picture an Astronomer series, which included both targeted public outreach to increase representation of women in astrophysics and high-level, solution-oriented discussions among professional astronomers. In March 2025, more than 200 astronomers came together in a hybrid-format symposium focused on the state of the field for female scientists, combining scientific exchange with discussions of policies and practices to strengthen retention of talent in the field. This white paper is the result of those discussions, offering a wide range of recommendations developed in the context of gendered attrition in astrophysics but which ultimately support a healthier climate for all scientists alike.

preprint2024arXiv

Dwarf-Dwarf Interactions Can Both Trigger and Quench Star Formation

It is exceedingly rare to find quiescent low-mass galaxies in the field. UGC5205 is an example of such a quenched field dwarf ($M_\star\sim3\times10^8M_\odot$). Despite a wealth of cold gas ($M_{\rm HI}\sim 3.5 \times 10^8 M_\odot$) and GALEX emission that indicates significant star formation in the past few hundred Myr, there is no detection of H$α$ emission -- star formation in the last $\sim 10$ Myr -- across the face of the galaxy. Meanwhile, the near equal-mass companion of UGC5205, PGC027864, is starbursting ($\rm EW_{\rm Hα}>1000$ Angstrom). In this work, we present new Karl G. Jansky Very Large Array (VLA) 21 cm line observations of UGC5205 that demonstrate that the lack of star formation is caused by an absence of HI in the main body of the galaxy. The HI of UGC5205 is highly disturbed; the bulk of the HI resides in several kpc-long tails, while the HI of PGC027864 is dominated by ordered rotation. We model the stellar populations of UGC5205 to show that, as indicated by the UV-H$α$ emission, the galaxy underwent a coordinated quenching event $\sim\!100-300$ Myr ago. The asymmetry of outcomes for UGC5205 and PGC027864 demonstrate that major mergers can both quench and trigger star formation in dwarfs. However, because the gas remains bound to the system, we suggest that such mergers only temporarily quench star formation. We estimate a total quenched time of $\sim 560$ Myr for UGC5205, consistent with established upper limits on the quenched fraction of a few percent for dwarfs in the field.

preprint2023arXiv

Multiwavelength Scrutiny of X-ray Sources in Dwarf Galaxies: ULXs versus AGN

Owing to their quiet evolutionary histories, nearby dwarf galaxies (stellar masses $M_\star \lesssim 3 \times 10^9 M_\odot$) have the potential to teach us about the mechanism(s) that 'seeded' the growth of supermassive black holes, and also how the first stellar mass black holes formed and interacted with their environments. Here, we present high spatial-resolution observations of three dwarf galaxies in the X-ray (Chandra), the optical/near-infrared (Hubble Space Telescope), and the radio (Karl G. Jansky Very Large Array). These three galaxies were previously identified as hosting candidate active galactic nuclei on the basis of lower resolution X-ray imaging. With our new observations, we find that X-ray sources in two galaxies (SDSS J121326.01+543631.6 and SDSS J122111.29+173819.1) are off nuclear and lack corresponding radio emission, implying they are likely luminous X-ray binaries. The third galaxy (Mrk 1434) contains two X-ray sources (each with $L_{\rm X} \approx 10^{40}$ erg s$^{-1}$) separated by 2".8, has a low metallicity (12 + log (O/H) = 7.8), and emits nebular \ion{He}{II} $λ$4686 line emission. The northern source has spatially coincident point-like radio emission at 9.0 GHz and extended radio emission at 5.5 GHz. We discuss X-ray binary interpretations (where an ultraluminous X-ray source blows a 'radio bubble') and active galactic nucleus interpretations (where a $\approx 4\times10^5 M_\odot$ black hole launches a jet). In either case, we find that the \ion{He}{II} emission cannot be photoionised by the X-ray source, unless the source was $\approx$30-90 times more luminous several hundred years ago.

preprint2023arXiv

UNCOVER: Illuminating the Early Universe -- JWST/NIRSpec Confirmation of $z > 12$ Galaxies

Observations of high-redshift galaxies provide a critical direct test to the theories of early galaxy formation, yet to date, only three have been spectroscopically confirmed at $z>12$. Due to strong gravitational lensing over a wide area, the galaxy cluster field A2744 is ideal for searching for the earliest galaxies. Here we present JWST/NIRSpec observations of two galaxies: a robust detection at $z_{\rm spec} = 12.393^{+0.004}_{-0.001}$, and a plausible candidate at $z_{\rm spec} = 13.079^{+0.013}_{-0.001}$. The galaxies are discovered in JWST/NIRCam imaging and their distances are inferred with JWST/NIRSpec spectroscopy, all from the JWST Cycle 1 UNCOVER Treasury survey. Detailed stellar population modeling using JWST NIRCam and NIRSpec data corroborates the primeval characteristics of these galaxies: low mass ($\sim 10^8~{\rm M_\odot}$), young, rapidly-assembling, metal-poor, and star-forming. Interestingly, both galaxies are spatially resolved, having lensing-corrected rest-UV effective radii on the order of 300-400 pc, which are notably larger than other spectroscopically confirmed systems at similar redshifts. The observed dynamic range of $z \gtrsim 10$ sizes spans over 1 order of magnitude, implying a significant scatter in the size-mass relation at early times. Deep into the epoch of reionization, these discoveries elucidate the emergence of the first galaxies.

preprint2022arXiv

Recovering the star formation histories of recently-quenched galaxies: the impact of model and prior choices

Accurate models of the star formation histories (SFHs) of recently-quenched galaxies can provide constraints on when and how galaxies shut down their star formation. The recent development of "non-parametric" SFH models promises the flexibility required to make these measurements. However, model and prior choices significantly affect derived SFHs, particularly for post-starburst galaxies (PSBs) which have sharp changes in their recent SFH. In this paper, we create mock PSBs, then use the Prospector SED fitting software to test how well four different SFH models recover key properties. We find that a two-component parametric model performs well for our simple mock galaxies, but is sensitive to model mismatches. The fixed- and flexible-bin non-parametric models included in Prospector are able to rapidly quench a major burst of star formation, but systematically underestimate the post-burst age by up to 200 Myr. We develop a custom SFH model that allows for additional flexibility in the recent SFH. Our flexible non-parametric model is able to constrain post-burst ages with no significant offset and just ~90 Myr of scatter. Our results suggest that while standard non-parametric models are able to recover first-order quantities of the SFH (mass, SFR, average age), accurately recovering higher-order quantities (burst fraction, quenching time) requires careful consideration of model flexibility. These mock recovery tests are a critical part of future SFH studies. Finally, we show that our new, public SFH model is able to accurately recover the properties of mock star-forming and quiescent galaxies and is suitable for broader use in the SED fitting community.

preprint2022arXiv

Star Formation Suppresion by Tidal Removal of Cold Molecular Gas from an Intermediate-Redshift Massive Post-Starburst Galaxy

Observations and simulations have demonstrated that star formation in galaxies must be actively suppressed to prevent the formation of over-massive galaxies. Galactic outflows driven by stellar feedback or supermassive black hole accretion are often invoked to regulate the amount of cold molecular gas available for future star formation, but may not be the only relevant quenching processes in all galaxies. We present the discovery of vast molecular tidal features extending up to 64 kpc outside of a massive z=0.646 post-starburst galaxy that recently concluded its primary star-forming episode. The tidal tails contain (1.2 +/- 0.1)x10^10 Msun of molecular gas, 47 +/- 5 % of the total cold gas reservoir of the system. Both the scale and magnitude of the molecular tidal features are unprecedented compared to all known nearby or high-redshift merging systems. We infer that the cold gas was stripped from the host galaxies during the merger, which is most likely responsible for triggering the initial burst phase and the subsequent suppression of star formation. While only a single example, this result shows that galaxy mergers can regulate the cold gas contents in distant galaxies by directly removing a large fraction of the molecular gas fuel, and plausibly suppress star formation directly, a qualitatively different physical mechanism than feedback-driven outflows.

preprint2022arXiv

The Compact Structures of Massive $z\sim0.7$ Post-Starburst Galaxies in the SQuIGG$\vec{L}$E Sample

We present structural measurements of 145 spectroscopically selected intermediate-redshift (z$\sim$0.7), massive ($M_\star \sim 10^{11} \ M_\odot$) post-starburst galaxies from the SQuIGG$\vec{L}$E Sample measured using wide-depth Hyper Suprime-Cam i-band imaging. This deep imaging allows us to probe the sizes and structures of these galaxies, which we compare to a control sample of star forming and quiescent galaxies drawn from the LEGA-C Survey. We find that post-starburst galaxies systematically lie $\sim0.1$ dex below the quiescent mass-size (half-light radius) relation, with a scatter of $\sim0.2$ dex. This finding is bolstered by non-parametric measures, such as the Gini coefficient and the concentration, which also reveal these galaxies to have more compact light profiles than both quiescent and star-forming populations at similar mass and redshift. The sizes of post-starburst galaxies show either negative or no correlation with the time since quenching, such that more recently quenched galaxies are larger or similarly sized. This empirical finding disfavors the formation of post-starburst galaxies via a purely central burst of star formation that simultaneously shrinks the galaxy and shuts off star formation. We show that the central densities of post-starburst and quiescent galaxies at this epoch are very similar, in contrast with their effective radii. The structural properties of z$\sim$0.7 post-starburst galaxies match those of quiescent galaxies that formed in the early universe, suggesting that rapid quenching in the present epoch is driven by a similar mechanism to the one at high redshift.

preprint2022arXiv

The Exploration of Local VolumE Satellites (ELVES) Survey: A Nearly Volume-Limited Sample of Nearby Dwarf Satellite Systems

We present the final results of the Exploration of Local VolumE Satellites (ELVES) Survey, a survey of the dwarf satellites of a nearly volume-limited sample of Milky Way (MW)-like hosts in the Local Volume. Hosts are selected simply via a cut in luminosity ($M_{K_s}<-22.1$ mag) and distance ($D<12$ Mpc). We have cataloged the satellites of 25 of the 31 such hosts, with another five taken from the literature. All hosts are surveyed out to at least 150 projected kpc ($\sim R_\mathrm{vir}/2$) with the majority surveyed to 300 kpc ($\sim R_\mathrm{vir}$). Satellites are detected using a consistent semi-automated algorithm specialized for low surface brightness dwarfs. As shown through extensive tests with injected galaxies, the catalogs are complete to $M_V\sim-9$ mag and $μ_{0,V}\sim26.5$ mag arcsec$^{-2}$. Candidates are confirmed to be real satellites through distance measurements including redshift, tip of the red giant branch, and surface brightness fluctuations. Across all 30 surveyed hosts, there are 338 confirmed satellites with a further 105 candidates awaiting distance measurement. For the vast majority of these, we provide consistent multi-band Sérsic photometry. We show that satellite abundance correlates with host mass, with the MW being quite typical amongst comparable systems, and that satellite quenched fraction rises steeply with decreasing satellite mass, mirroring the quenched fraction for the MW and M31. The ELVES survey represents a massive increase in the statistics of surveyed systems with known completeness, and the provided catalogs are a unique dataset to explore various aspects of small-scale structure and dwarf galaxy evolution.

preprint2022arXiv

The In-situ Origins of Dwarf Stellar Outskirts in FIRE-2

Extended, old, and round stellar halos appear to be ubiquitous around high-mass dwarf galaxies ($10^{8.5}<M_\star/M_\odot<10^{9.6}$) in the observed universe. However, it is unlikely that these dwarfs have undergone a sufficient number of minor mergers to form stellar halos that are composed of predominantly accreted stars. Here, we demonstrate that FIRE-2 (Feedback in Realistic Environments) cosmological zoom-in simulations are capable of producing dwarf galaxies with realistic structure, including both a thick disk and round stellar halo. Crucially, these stellar halos are formed in-situ, largely via the outward migration of disk stars. However, there also exists a large population of "non-disky" dwarfs in FIRE-2 that lack a well-defined disk/halo and do not resemble the observed dwarf population. These non-disky dwarfs tend to be either more gas poor or to have burstier recent star formation histories than the disky dwarfs, suggesting that star formation feedback may be preventing disk formation. Both classes of dwarfs underscore the power of a galaxy's intrinsic shape -- which is a direct quantification of the distribution of the galaxy's stellar content -- to interrogate the feedback implementation in simulated galaxies.

preprint2022arXiv

The MASSIVE Survey -- XVII. A Triaxial Orbit-based Determination of the Black Hole Mass and Intrinsic Shape of Elliptical Galaxy NGC 2693

We present a stellar dynamical mass measurement of a newly detected supermassive black hole (SMBH) at the center of the fast-rotating, massive elliptical galaxy NGC 2693 as part of the MASSIVE survey. We combine high signal-to-noise integral field spectroscopy (IFS) from the Gemini Multi-Object Spectrograph (GMOS) with wide-field data from the Mitchell Spectrograph at McDonald Observatory to extract and model stellar kinematics of NGC 2693 from the central $\sim 150$ pc out to $\sim2.5$ effective radii. Observations from Hubble Space Telescope (HST) WFC3 are used to determine the stellar light distribution. We perform fully triaxial Schwarzschild orbit modeling using the latest TriOS code and a Bayesian search in 6-D galaxy model parameter space to determine NGC 2693's SMBH mass ($M_\text{BH}$), stellar mass-to-light ratio, dark matter content, and intrinsic shape. We find $M_\text{BH} = \left(1.7\pm 0.4\right)\times 10^{9}\ M_\odot$ and a triaxial intrinsic shape with axis ratios $p=b/a=0.902 \pm 0.009$ and $q=c/a=0.721^{+0.011}_{-0.010}$, triaxiality parameter $T = 0.39 \pm 0.04$. In comparison, the best-fit orbit model in the axisymmetric limit and (cylindrical) Jeans anisotropic model of NGC 2693 prefer $M_\text{BH} = \left(2.4\pm 0.6\right)\times 10^{9}\ M_\odot$ and $M_\text{BH} = \left(2.9\pm 0.3\right)\times 10^{9}\ M_\odot$, respectively. Neither model can account for the non-axisymmetric stellar velocity features present in the IFS data.

preprint2022arXiv

Towards a More Complete Optical Census of Active Galactic Nuclei, Via Spatially-Resolved Spectroscopy

While emission-line flux ratio diagnostics are the most common technique for identifying active galactic nuclei (AGNs) in optical spectra, applying this approach to single fiber spectra of galaxies can omit entire subpopulations of AGNs. Here, we use spatially resolved spectroscopy from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey to construct a sample of 10 galaxies where Baldwin-Philips-Terlevich line flux ratio diagnostics classify each galaxy's central $3^{\prime\prime}$ spectrum as LINER or star forming, while $>10\%$ of the spaxels in the galaxy's MaNGA footprint are classified as Seyfert. We obtain Chandra observations of these 10 galaxies with off-nuclear Seyfert regions to determine whether AGNs are actually present in them. Our main result is that 7-10 (depending on strictness of criteria) of the galaxies host one or more X-ray AGNs, even though none of them were classified as AGNs based on their single-fiber optical spectra. We find that these AGNs were not identified in the single-fiber spectra because they are AGNs in the nuclei of companion galaxies, low luminosity AGNs, dust obscured AGNs, and/or flickering AGNs. In summary, we find that off-nuclear AGN signatures may increase the number of known AGNs by a factor of two over what conventional single nuclear fiber spectra identify. Our results show that spatially resolved spectroscopy can be leveraged to reveal a more complete census of AGNs that are traditionally missed by single fiber spectra.

preprint2022arXiv

Wandering Black Hole Candidates in Dwarf Galaxies at VLBI Resolution

Thirteen dwarf galaxies have recently been found to host radio-selected accreting massive black hole (MBH) candidates, some of which are ``wandering" in the outskirts of their hosts. We present 9 GHz Very Long Baseline Array (VLBA) observations of these sources at milliarcsecond resolution. Our observations have beam solid angles ${\sim}10^4$ times smaller than the previous Very Large Array (VLA) observations at 9 GHz, with comparable point source sensitivities. We detect milliarcsecond-scale radio sources at the positions of the four VLA sources most distant from the photo-centers of their associated dwarf galaxies. These sources have brightness temperatures of ${>}10^6~\mathrm{K}$, consistent with active galactic nuclei (AGNs), but the significance of their preferential location at large distances ($p$-value~$=0.0014$) favors a background AGN interpretation. The VLBA non-detections toward the other 9 galaxies indicate that the VLA sources are resolved out on scales of tens of milliarcseconds, requiring extended radio emission and lower brightness temperatures consistent with either star formation or radio lobes associated with AGN activity. We explore the star formation explanation by calculating the expected radio emission for these nine VLBA non-detections, finding that about 5 have VLA luminosities that are inconsistent with this scenario. Of the remaining four, two are associated with spectroscopically confirmed AGNs that are consistent with being located at their galaxy photo-centers. There are therefore between 5 and 7 wandering MBH candidates out of the 13 galaxies we observed, although we cannot rule out background AGNs for five of them with the data in hand.

preprint2021arXiv

A Chandra and HST View of WISE-Selected AGN Candidates in Dwarf Galaxies

Reliably identifying active galactic nuclei (AGNs) in dwarf galaxies is key to understanding black hole demographics at low masses and constraining models for black hole seed formation. Here we present Chandra X-ray Observatory observations of eleven dwarf galaxies that were chosen as AGN candidates using Wide-field Infrared Survey Explorer (WISE) mid-infrared (mid-IR) color-color selection. Hubble Space Telescope images are also presented for ten of the galaxies. Based on Sloan Digital Sky Survey spectroscopy, six galaxies in our sample have optical evidence for hosting AGNs and five are classified as star-forming. We detect X-ray point sources with luminosities above that expected from X-ray binaries in the nuclei of five of the six galaxies with optical evidence of AGNs. However, the X-ray emission from these AGNs is generally much lower than expected based on AGN scaling relations with infrared and optical tracers. We do not find compelling evidence for AGNs in the five optically-selected star-forming galaxies despite having red mid-IR colors. Only two are detected in X-rays and their properties are consistent with stellar-mass X-ray binaries. Based on this multiwavelength study, we conclude that two-color mid-IR AGN diagnostics at the resolution of WISE cannot be used to reliably select AGNs in optically-star-forming dwarf galaxies. Future observations in the infrared with the James Webb Space Telescope offer a promising path forward.

preprint2021arXiv

A quasar-based supermassive black hole binary population model: implications for the gravitational-wave background

The nanohertz gravitational wave background (GWB) is believed to be dominated by GW emission from supermassive black hole binaries (SMBHBs). Observations of several dual active galactic nuclei (AGN) strongly suggest a link between AGN and SMBHBs, given that these dual AGN systems will eventually form bound binary pairs. Here we develop an exploratory SMBHB population model based on empirically constrained quasar populations, allowing us to decompose the GWB amplitude into an underlying distribution of SMBH masses, SMBHB number density, and volume enclosing the GWB. Our approach also allows us to self-consistently predict the number of local SMBHB systems from the GWB amplitude. Interestingly, we find the local number density of SMBHBs implied by the common-process signal in the NANOGrav 12.5-yr dataset to be roughly five times larger than previously predicted by other models. We also find that at most $\sim 25 \%$ of SMBHBs can be associated with quasars. Furthermore, our quasar-based approach predicts $\gtrsim 95\%$ of the GWB signal comes from $z \lesssim 2.5$, and that SMBHBs contributing to the GWB have masses $\gtrsim 10^8 M_\odot$. We also explore how different empirical galaxy-black hole scaling relations affect the local number density of GW sources, and find that relations predicting more massive black holes decrease the local number density of SMBHBs. Overall, our results point to the important role that a measurement of the GWB will play in directly constraining the cosmic population of SMBHBs, as well as their connections to quasars and galaxy mergers.

preprint2021arXiv

CLIMBER: Galaxy-Halo Connection Constraints from Next-Generation Surveys

In the coming decade, a new generation of massively multiplexed spectroscopic surveys, such as PFS, WAVES, and MOONS, will probe galaxies in the distant universe in vastly greater numbers than was previously possible. In this work, we generate mock catalogs for each of these three planned surveys to help quantify and optimize their scientific output. To assign photometry into the UniverseMachine empirical model, we develop the Calibrating Light: Illuminating Mocks By Empirical Relations (CLIMBER) procedure using UltraVISTA photometry. Using the published empirical selection functions for each aforementioned survey, we quantify the mass completeness of each survey. We compare different targeting strategies by varying the area and targeting completeness, and quantify how these survey parameters affect the uncertainty of the two-point correlation function. We demonstrate that the PFS and MOONS measurements will be primarily dominated by cosmic variance, not shot noise, motivating the need for increasingly large survey areas. On the other hand, the WAVES survey, which covers a much larger area, will strike a good balance between cosmic variance and shot noise. For a fixed number of targets, a 5% increased survey area (and $\sim$5% decreased completeness) would decrease the uncertainty of the correlation function at intermediate scales by 0.15%, 1.2%, and 1.1% for our WAVES, PFS, and MOONS samples, respectively. Meanwhile, for a fixed survey area, 5% increased targeting completeness improves the same constraints by 0.7%, 0.25%, and 0.1%. All of the utilities used to construct our mock catalogs and many of the catalogs themselves are publicly available.

preprint2021arXiv

Deep Realistic Extragalactic Model (DREaM) Galaxy Catalogs: Predictions for a Roman Ultra-Deep Field

In the next decade, deep galaxy surveys from telescopes such as the James Webb Space Telescope and Roman Space Telescope will provide transformational data sets that will greatly enhance the understanding of galaxy formation during the epoch of reionization (EoR). In this work, we present the Deep Realistic Extragalactic Model (DREaM) for creating synthetic galaxy catalogs. Our model combines dark matter simulations, subhalo abundance matching and empirical models, and includes galaxy positions, morphologies, and spectral energy distributions (SEDs). The resulting synthetic catalog extends to redshifts $z \sim 12 $, and galaxy masses $\log_{10}(M/M_{\odot}) = 5 $ covering an area of $1\, {\rm deg}^2$ on the sky. We use DREaM to explore the science returns of a $1\, {\rm deg}^2$ \emph{Roman} UDF, and to provide a resource for optimizing ultra-deep survey designs. We find that a \emph{Roman} UDF to $\sim 30\,m_{\rm AB}$ will potentially detect more than $10^6$ $M_{\rm UV}<-17$ galaxies, with more than $10^4$ at redshifts $z>7$, offering an unparalleled dataset for constraining galaxy properties during the EoR. Our synthetic catalogs and simulated images are made publicly available to provide the community with a tool to prepare for upcoming data.

preprint2021arXiv

ELVES II: GCs and Nuclear Star Clusters of Dwarf Galaxies; The Importance of Environment

We present the properties of the globular clusters (GCs) and nuclear star clusters (NSCs) of low-mass ($10^{5.5}<M_\star<10^{8.5}$ $M_\odot$) early-type satellites of Milky Way-like and small group hosts in the Local Volume (LV) using deep, ground-based data from the ongoing Exploration of Local VolumE Satellites (ELVES) Survey. This sample of 177 dwarfs significantly increases the statistics for studying the star clusters of dwarfs in low-density environments, offering an important comparison to samples from nearby galaxy clusters. The LV dwarfs exhibit significantly lower nucleation fractions at fixed galaxy mass than dwarfs in nearby clusters. The mass of NSCs of LV dwarfs show a similar scaling of $M_{\star,\mathrm{NSC}}\propto M_{\star,\mathrm{gal}}^{0.4}$ as that found in clusters but offset to lower NSC masses. To deal with foreground/background contamination in the GC analysis, we employ both a statistical subtraction and Bayesian approach to infer the average GC system properties from all dwarfs simultaneously. We find that the GC occupation fraction and average abundance are both increasing functions of galaxy stellar mass, and the LV dwarfs show significantly lower average GC abundance at fixed galaxy mass than a comparable sample of Virgo dwarfs analyzed in the same way, demonstrating that GC prevalence also shows an important secondary dependence on the dwarf's environment. This result strengthens the connection between GCs and NSCs in low-mass galaxies. We discuss these observations in the context of modern theories of GC and NSC formation, finding that the environmental dependencies can be well-explained by these models.

preprint2021arXiv

Now you see it, now you don't: Star formation truncation precedes the loss of molecular gas by ~100 Myr in massive post-starburst galaxies at z~0.6

We use ALMA observations of CO(2-1) in 13 massive ($M_{\star}\gtrsim 10^{11} M_{\odot}$) post-starburst galaxies at $z\sim0.6$ to constrain the molecular gas content in galaxies shortly after they quench their major star-forming episode. The post-starburst galaxies in this study are selected from the Sloan Digital Sky Survey spectroscopic samples (DR14) based on their spectral shapes, as part of the SQuIGGLE program. Early results showed that two post-starburst galaxies host large H$_2$ reservoirs despite their low inferred star formation rates. Here we expand this analysis to a larger statistical sample of 13 galaxies. Six of the primary targets (45%) are detected, with $M_{H_2}\gtrsim10^9 M_{\odot}$. Given their high stellar masses, this mass limit corresponds to an average gas fraction of $\langle f_{H_2} \equiv M_{H_2}/M_{\star} \rangle \sim7\%$, or ${\sim}14\%$ using lower stellar masses estimates derived from analytic, exponentially declining star formation histories. The gas fraction correlates with the $D_n4000$ spectral index, suggesting that the cold gas reservoirs decrease with time since burst, as found in local K+A galaxies. Star formation histories derived from flexible stellar population synthesis modeling support this empirical finding: galaxies that quenched $\lesssim 150$ Myr prior to observation host detectable CO(2-1) emission, while older post-starburst galaxies are undetected. The large $\mathrm{H_2}$ reservoirs and low star formation rates in the sample imply that the quenching of star formation precedes the disappearance of the cold gas reservoirs. However, within the following 100-200 Myrs, the SQuIGGLE galaxies require the additional and efficient heating or removal of cold gas to bring their low star formation rates in line with standard $\mathrm{H_2}$ scaling relations.

preprint2021arXiv

Reaching for the Edge I: Probing the Outskirts of Massive Galaxies with HSC, DECaLS, SDSS, and Dragonfly

The outer light (stellar halos) of massive galaxies has recently emerged as a possible low scatter tracer of dark matter halo mass. To test the robustness of outer light measurements across different data sets, we compare the surface brightness profiles of massive galaxies using four independent data sets: the Hyper Suprime-Cam survey (HSC), the Dark Energy Camera Legacy Survey (DECaLS), the Sloan Digital Sky Survey (SDSS), and the Dragonfly Wide Field Survey (Dragonfly). We use customized pipelines for HSC and DECaLS to achieve better sky background subtraction. For galaxies at $z<0.05$, Dragonfly has the best control of systematics, reaching surface brightness levels of $μ_r \sim 30$ mag/arcsec$^{2}$. At $0.19<z<0.50$, HSC can reliably recover surface brightness profiles to $μ_{r} \sim 28.5$ mag/arcsec$^{2}$ reaching $R=100 - 150$ kpc. DECaLS surface brightness profiles show good agreement with HSC but are noisier at large radii. The median profiles of galaxy ensembles in both HSC and DECaLS reach $R > 200$ kpc without significant bias. At $0.19<z<0.50$, DECaLS and HSC measurements of the stellar mass contained within 100 kpc agree within 0.05 dex. Finally, we use weak gravitational lensing to show that measurements of outer light with DECaLS at $0.19<z<0.50$ show a similar promise as HSC as a low scatter proxy of halo mass. The tests and results from this paper represent an important step forward for accurate measurements of the outer light of massive galaxies and demonstrate that outer light measurements from DECam imaging will be a promising method for finding galaxy clusters for DES and DESI.

preprint2021arXiv

SQuIGGLE: Studying Quenching in Intermediate-z Galaxies -- Gas, AnguLar Momentum, and Evolution

We describe the SQuIGGLE survey of intermediate-redshift post-starburst galaxies. We leverage the large sky coverage of the SDSS to select ~1300 recently-quenched galaxies at 0.5<z<~0.9 based on their unique spectral shapes. These bright, intermediate-redshift galaxies are ideal laboratories to study the physics responsible for the rapid quenching of star formation: they are distant enough to be useful analogs for high-redshift quenching galaxies, but low enough redshift that multi-wavelength follow-up observations are feasible with modest telescope investments. We use the Prospector code to infer the stellar population properties and non-parametric star formation histories of all galaxies in the sample. We find that SQuIGGLE galaxies are both very massive (M* ~ 10^11.25 Msun) and quenched, with inferred star formation rates <~1Msun/yr, more than an order of magnitude below the star-forming main sequence. The best-fit star formation histories confirm that these galaxies recently quenched a major burst of star formation: >75% of SQuIGGLE galaxies formed at least a quarter of their total stellar mass in the recent burst, which ended just ~200Myr before observation. We find that SQuIGGLE galaxies are on average younger and more burst-dominated than most other z<~1 post-starburst samples. This large sample of bright post-starburst galaxies at intermediate redshift opens a wide range of studies into the quenching process. In particular, the full SQuIGGLE survey will investigate the molecular gas reservoirs, morphologies, kinematics, resolved stellar populations, AGN incidence, and infrared properties of this unique sample of galaxies in order to place definitive constraints on the quenching process.

preprint2020arXiv

A Second Look at 12 Candidate Dual AGNs using ${\tt BAYMAX}$

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] $λ$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.

preprint2020arXiv

Ionised gas outflow signatures in SDSS-IV MaNGA active galactic nuclei

The prevalence of outflow and feedback signatures in active galactic nuclei (AGN) is a major unresolved question which large integral field unit (IFU) surveys now allow to address. In this paper, we present a kinematic analysis of the ionised gas in 2778 galaxies at z $\sim$ 0.05 observed by SDSS-IV MaNGA. Specifically, we measure the kinematics of the [OIII] λ5007Å emission line in each spatial element and fit multiple Gaussian components to account for possible non-gravitational motions of gas. Comparing the kinematics of the ionised gas between 308 MaNGA-selected AGN that have been previously identified through emission line diagnostics and sources not classified as AGN, we find that while 25% of MaNGA-selected AGN show [OIII] components with emission line widths of $>$ 500 km/s in more than 10% of their spaxels, only 7% of MaNGA non-AGN show a similar signature. Even the AGN that do not show nuclear AGN photoionisation signatures and that were only identified as AGN based on their larger scale photoionisation signatures show similar kinematic characteristics. In addition to obscuration, another possibility is that outflow and mechanical feedback signatures are longer lived than the AGN itself. Our measurements demonstrate that high velocity gas is more prevalent in AGN compared to non-AGN and that outflow and feedback signatures in low-luminosity, low-redshift AGN may so far have been underestimated. We show that higher luminosity MaNGA-selected AGN are able to drive larger scale outflows than lower luminosity AGN. But estimates of the kinetic coupling efficiencies are $<$ 1% and suggest that the feedback signatures probed in this paper are unlikely to have a significant impact on the AGN host galaxies. However, continuous energy injection may still heat a fraction of the cool gas and delay or suppress star formation in individual galaxies even when the AGN is weak.

preprint2020arXiv

Luminosity Functions and Host-to-Host Scatter of Dwarf Satellite Systems in the Local Volume

Low-mass satellites around Milky Way (MW)-like galaxies are important probes of small scale structure and galaxy formation. However, confirmation of satellite candidates with distance measurements remains a key barrier to fast progress in the Local Volume (LV). We measure the surface brightness fluctuation (SBF) distances to recently cataloged candidate dwarf satellites around 10 massive hosts within $D<12$ Mpc to confirm association. The satellite systems of these hosts are complete and mostly cleaned of contaminants down to $M_g{\sim}-9$ to $-10$, within the area of the search footprints. Joining this sample with hosts surveyed to comparable or better completeness in the literature, we explore how well cosmological simulations combined with common stellar to halo mass relations (SHMR) match observed satellite luminosity functions in the classical satellite luminosity regime. Adopting a SHMR that matches hydrodynamic simulations, the predicted overall satellite abundance agrees well with the observations. The MW is remarkably typical in its luminosity function amongst LV hosts. Contrary to recent results, we find that the host-to-host scatter predicted by the model is in close agreement with the scatter between the observed systems, once the different masses of the observed systems are taken into account. However, we find significant evidence that the observed systems have more bright and fewer faint satellites than the SHMR model predicts, necessitating a higher normalization of the SHMR around halo masses of $10^{11}$ \msun\ than present in common SHMRs. These results demonstrate the utility of nearby satellite systems in inferring the galaxy-subhalo connection in the low-mass regime.

preprint2020arXiv

Observational nonstationarity of AGN variability: The only way to go is down!

To gain insights into long-term Active Galactic Nuclei (AGN) variability, we analyze an AGN sample from the Sloan Digital Sky Survey (SDSS) and compare their photometry with observations from the Hyper Suprime-Cam survey (HSC) observed $\langle 14.85 \rangle$ years after SDSS. On average, the AGN are fainter in HSC than SDSS. We demonstrate that the difference is not due to subtle differences in the SDSS versus HSC filters or photometry. The decrease in mean brightness is redshift dependent, consistent with expectations for a change that is a function of the rest-frame time separation between observations. At a given redshift, the mean decrease in brightness is stronger for more luminous AGN and for objects with longer time separation between measurements. We demonstrate that the dependence on redshift and luminosity of measured mean brightness decrease is consistent with simple models of Eddington ratio variability in AGN on long (Myr, Gyr) timescales. We show how our results can be used to constrain the variability and demographic properties of AGN populations.

preprint2020arXiv

Properties of simulated galaxies and supermassive black holes in cosmic voids

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.

preprint2020arXiv

Revealing the intermediate-mass black hole at the heart of the dwarf galaxy NGC 404 with sub-parsec resolution ALMA observations

We estimate the mass of the intermediate-mass black hole at the heart of the dwarf elliptical galaxy NGC 404 using Atacama Large Millimeter/submillimeter Array (ALMA) observations of the molecular interstellar medium at an unprecedented linear resolution of ~0.5 pc, in combination with existing stellar kinematic information. These ALMA observations reveal a central disc/torus of molecular gas clearly rotating around the black hole. This disc is surrounded by a morphologically and kinematically complex flocculent distribution of molecular clouds, that we resolve in detail. Continuum emission is detected from the central parts of NGC 404, likely arising from the Rayleigh-Jeans tail of emission from dust around the nucleus, and potentially from dusty massive star-forming clumps at discrete locations in the disc. Several dynamical measurements of the black hole mass in this system have been made in the past, but they do not agree. We show here that both the observed molecular gas and stellar kinematics independently require a ~5x10$^5$ Msun black hole once we include the contribution of the molecular gas to the potential. Our best estimate comes from the high-resolution molecular gas kinematics, suggesting the black hole mass of this system is 5.5$^{+4.1}_{-3.8}\times$10$^5$ Msun (at the 99% confidence level), in good agreement with our revised stellar kinematic measurement and broadly consistent with extrapolations from the black hole mass - velocity dispersion and black hole mass - bulge mass relations. This highlights the need to accurately determine the mass and distribution of each dynamically important component around intermediate-mass black holes when attempting to estimate their masses.

preprint2020arXiv

Spatially-resolved UV diagnostics of AGN feedback: radiation pressure dominates in a prototypical quasar-driven superwind

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.

preprint2020arXiv

The MASSIVE Survey -- XV. A Stellar Dynamical Mass Measurement of the Supermassive Black Hole in Massive Elliptical Galaxy NGC 1453

We present a new stellar dynamical mass measurement ($M_\mathrm{BH}$) of the supermassive black hole (SMBH) in NGC 1453, a fast-rotating massive elliptical galaxy in the MASSIVE survey. We measure stellar kinematics in 135 spatial bins in the central 1.5 kpc by 2 kpc region of the galaxy using high signal-to-noise ($S/N \sim 130$) spectra from the Gemini-North GMOS integral-field spectrograph (IFS). Combining with wide-field IFS kinematics out to $\sim 3$ effective radii and stellar light distributions from Hubble Space Telescope (HST) WFC3 images, we perform Schwarzschild orbit-based mass modeling in the axisymmetric limit to constrain the mass components in NGC 1453. The best-fit black hole mass is $M_\mathrm{BH} =(2.9 \pm 0.4) 10^9 M_\odot$; the mass models without a central black hole are excluded at the $8.7σ$ level. The NGC 1453 black hole lies within the intrinsic scatter of the SMBH and galaxy scaling relations, unlike three other galaxies hosting $\gtrsim 10^{10} M_\odot$ SMBHs in the MASSIVE sample. The high-$S/N$ GMOS spectra enable us to determine 8 moments of the Gauss-Hermite expansion of the line-of-sight velocity distributions (LOSVDs), which are used as constraints in the orbit modeling. The stellar orbits in the mass models are further constrained to produce negligible $h_9$ through $h_{12}$ to minimize spurious behavior in the LOSVDs. We show that truncating the series at $h_4$, as was often done in prior work, leads to a much weaker constraint on the inferred $M_\mathrm{BH}$ for NGC 1453. Furthermore, we discuss precautions and modifications that are needed to achieve axisymmetry in triaxial orbit codes that use the Schwarzschild method to sample the start space of stellar orbits in triaxial gravitational potentials.

preprint2020arXiv

The MASSIVE Survey XIV -- Stellar Velocity Profiles and Kinematic Misalignments from 200 pc to 20 kpc in Massive Early-type Galaxies

We use high spatial resolution stellar velocity maps from the Gemini GMOS integral-field spectrograph (IFS) and wide-field velocity maps from the McDonald Mitchell IFS to study the stellar velocity profiles and kinematic misalignments from $\sim 200$ pc to $\sim 20$ kpc in 20 early-type galaxies with stellar mass $M_* > 10^{11.7} M_\odot$ in the MASSIVE survey. While 80% of the galaxies have low spins ($λ< 0.1$) and low rotational velocities ($< 50$ km/s) in both the central region and the main body, we find a diverse range of velocity features and misalignment angles. For the 18 galaxies with measurable central kinematic axes, 10 have well aligned kinematic axis and photometric major axis, and the other 8 galaxies have misalignment angles that are distributed quite evenly from $15^\circ$ to the maximal value of $90^\circ$. There is a strong correlation between central kinematic misalignment and galaxy spin, where all 4 galaxies with significant spins have well aligned kinematic and photometric axes, but only 43% of the low-spin galaxies are well aligned. The central and main-body kinematic axes within a galaxy are not always aligned. When the two kinematic axes are aligned ($\sim 60$% of the cases), they are either also aligned with the photometric major axis or orthogonal to it. We find 13 galaxies to also exhibit noticeable local kinematic twists, and one galaxy to have a counter-rotating core. A diverse assembly history consisting of multiple gas-poor mergers of a variety of progenitor mass ratios and orbits is likely to be needed to account for the predominance of low spins and the wide range of central and main-body velocity features reported here for local massive ETGs.

preprint2020arXiv

The Role of Active Galactic Nuclei in the Quenching of Massive Galaxies in the SQuiGGLE Survey

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.

preprint2020arXiv

Tracing the Intrinsic Shapes of Dwarf Galaxies out to Four Effective Radii: Clues to Low-Mass Stellar Halo Formation

Though smooth, extended spheroidal stellar outskirts have long been observed around nearby dwarf galaxies, it is unclear whether dwarfs generically host an extended stellar halo. We use imaging from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) to measure the shapes of dwarf galaxies out to four effective radii for a sample of dwarfs at 0.005<z<0.2 and 10^7<M_star/M_sun<10^9.6. We find that dwarfs are slightly triaxial, with a <B/A> >~ 0.75 (where the ellipsoid is characterized by three principle semi-axes constrained by C<=B<=A). At M_star>10^8.5 M_sun, the galaxies grow from thick disk-like near their centers towards the spheroidal extreme at four effective radii. We also see that although blue dwarfs are, on average, characterized by thinner discs than red dwarfs, both blue and red dwarfs grow more spheroidal as a function of radius. This relation also holds true for a comparison between field and satellite dwarfs. This uniform trend towards relatively spheroidal shapes as a function of radius is consistent with an in-situ formation mechanism for stellar outskirts around low-mass galaxies, in agreement with proposed models where star formation feedback produces round stellar outskirts around dwarfs.

preprint2019arXiv

Star Formation in Isolated Dwarf Galaxies Hosting Tidal Debris: Extending the Dwarf-Dwarf Merger Sequence

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$α$ equivalent width (H$α$ 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.

preprint2019arXiv

Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) VIII. A less biased view of the early co-evolution of black holes and host galaxies

We present ALMA [CII] line and far-infrared (FIR) continuum observations of three $z > 6$ low-luminosity quasars ($M_{\rm 1450} > -25$) discovered by our Subaru Hyper Suprime-Cam (HSC) survey. The [CII] line was detected in all three targets with luminosities of $(2.4 - 9.5) \times 10^8~L_\odot$, about one order of magnitude smaller than optically luminous ($M_{\rm 1450} \lesssim -25$) quasars. The FIR continuum luminosities range from $< 9 \times 10^{10}~L_\odot$ (3$σ$ limit) to $\sim 2 \times 10^{12}~L_\odot$, indicating a wide range in star formation rates in these galaxies. Most of the HSC quasars studied thus far show [CII]/FIR luminosity ratios similar to local star-forming galaxies. Using the [CII]-based dynamical mass ($M_{\rm dyn}$) as a surrogate for bulge stellar mass ($M_{\rm bulge}$), we find that a significant fraction of low-luminosity quasars are located on or even below the local $M_{\rm BH} - M_{\rm bulge}$ relation, particularly at the massive end of the galaxy mass distribution. In contrast, previous studies of optically luminous quasars have found that black holes are overmassive relative to the local relation. Given the low luminosities of our targets, we are exploring the nature of the early co-evolution of supermassive black holes and their hosts in a less biased way. Almost all of the quasars presented in this work are growing their black hole mass at much higher pace at $z \sim 6$ than the parallel growth model, in which supermassive black holes and their hosts grow simultaneously to match the local $M_{\rm BH} - M_{\rm bulge}$ relation at all redshifts. As the low-luminosity quasars appear to realize the local co-evolutionary relation even at $z \sim 6$, they should have experienced vigorous starbursts prior to the currently observed quasar phase to catch up with the relation.

preprint2019arXiv

Surprisingly Strong K-band Emission Found in Low Luminosity Active Galactic Nuclei

We examine the near-infrared (NIR) emission from low-luminosity AGNs (LLAGNs). Our galaxy sample includes 15 objects with detected 2-10 keV X-ray emission, dynamical black hole mass estimates from the literature, and available Gemini/NIFS integral field spectroscopy (IFU) data. We find evidence for red continuum components at the center of most galaxies, consistent with the hot dust emission seen in higher luminosity AGN. We decompose the spectral data cubes into a stellar and continuum component, assuming the continuum component comes from thermal emission from hot dust. We detect nuclear thermal emission in 14 out of 15 objects. This emission causes weaker CO absorption lines and redder continuum ($2.05-2.28\:μ$m) in our $K$-band data, as expected from hot dust around an AGN. The NIR emission is clearly correlated with the 2-10 keV X-ray flux, with a Spearman coefficient of $r_{spearman}=0.69$ suggesting a $>99\%$ significance of correlation, providing further evidence of an AGN origin. Our sample has typical X-ray and NIR fluxes $3-4$ orders of magnitude less luminous than previous work studying the NIR emission from AGN. We find that the ratio of NIR to X-ray emission increases towards lower Eddington ratios. The NIR emission in our sample is often brighter than the X-ray emission, with our $K$-band AGN luminosities comparable to or greater than the 2-10 keV X-ray luminosities in all objects with Eddington ratios below $0.01\%$. The nature of this LLAGN NIR emission remains unclear, with one possibility being an increased contribution from jet emission at these low luminosities. These observations suggest JWST will be a useful tool for detecting the lowest luminosity AGN.

preprint2019arXiv

The MBHBM* Project -- I: Measurement of the Central Black Hole Mass in Dwarf Galaxy NGC 3504 Using Molecular Gas Kinematics

We present the first measurement of the mass of a supermassive black hole (SMBH) in the nearby double-barred spiral galaxy NGC 3504 as part of the Measuring Black Holes Below the Milky Way ($M_{\star}$) mass galaxies (MBHBM$_{\star}$) Project. Our analysis is based on Atacama Large Millimeter/submillimeter Array (ALMA) Cycle-5 observations of the ${\rm ^{12}CO(2-1)}$ emission line. NGC 3504 has a circumnuclear gas disk (CND), which has a relatively high-velocity dispersion of 30 \kms. Our dynamical models of the CND yield a \Mbh~of $M_{\rm BH}=1.02^{+0.18}_{-0.15}\times10^7$\Msun and a mass-to-light ratio in $H$-band of \ml$_{\rm H}=0.66^{+1.44}_{-0.65}$ (\Msun/\Lsun). This black hole (BH) mass is consistent with BH--galaxy scaling relations. We also detect a central deficiency in the ${\rm ^{12}CO(2-1)}$ integrated intensity map with a diameter of 2.7 pc at the putative position of the SMBH. However, this hole is filled by a dense gas tracer ${\rm CS(5-4)}$ that peaks at the galaxy center found in one of the three low-velocity-resolution continuum spectral correlators. The ${\rm CS(5-4)}$ line has the same kinematics with the ${\rm ^{12}CO(2-1)}$ line within the CND, suggesting that it is also an alternative transition for measuring the central \Mbh~in NGC 3504 probably more accurately than the current commonly used of ${\rm ^{12}CO(2-1)}$ due to its centralization.

preprint2019arXiv

Wide-Field Survey of Dwarf Satellite Systems Around 10 Hosts in the Local Volume

We present the results of an extensive search for dwarf satellite galaxies around 10 primary host galaxies in the Local Volume (D$<$12 Mpc) using archival CFHT/MegaCam imaging data. The hosts span a wide range in properties, with stellar masses ranging from that of the LMC to ${\sim}3$ times that of the Milky Way (MW). The surveyed hosts are: NGC 1023, NGC 1156, NGC 2903, NGC 4258, NGC 4565, NGC 4631, NGC 5023, M51, M64, and M104. We detect satellite candidates using a consistent semi-automated detection algorithm that is optimized for the detection of low surface brightness objects. Depending on the host, our completeness limit is $M_g{\sim}-8$ to $-10$ (assuming the distance of the host). We detect objects with surface brightness down to $μ_{0,g}{\sim}26$ mag arcsec$^{-2}$ at $\gtrsim90\%$ completeness. The survey areas of the six best-surveyed hosts cover most of the inner projected $R<150$ kpc area, which roughly doubles the number of MW-mass hosts surveyed at this level of area and luminosity completeness. The number of detected candidates range from 1 around M64 to 33 around NGC 4258. In total, 153 candidates are found, of which 93 are new. While we defer an analysis of the satellite luminosity functions of the hosts until distance information is available for the candidates, we do show that the candidates are primarily red, spheroid systems with properties roughly consistent with known satellites in the Local Group.

preprint2016arXiv

A 17-billion-solar-mass black hole in a group galaxy with a diffuse core

Quasars are associated with and powered by the accretion of material onto massive black holes; the detection of highly luminous quasars with redshifts greater than z = 6 suggests that black holes of up to ten billion solar masses already existed 13 billion years ago. Two possible present-day dormant descendants of this population of active black holes have been found in the galaxies NGC 3842 and NGC 4889 at the centres of the Leo and Coma galaxy clusters, which together form the central region of the Great Wall - the largest local structure of galaxies. The most luminous quasars, however, are not confined to such high-density regions of the early Universe; yet dormant black holes of this high mass have not yet been found outside of modern-day rich clusters. Here we report observations of the stellar velocity distribution in the galaxy NGC 1600 - a relatively isolated elliptical galaxy near the centre of a galaxy group at a distance of 64 Mpc from Earth. We use orbit superposition models to determine that the black hole at the centre of NGC 1600 has a mass of 17 billion solar masses. The spatial distribution of stars near the centre of NGC 1600 is rather diffuse. We find that the region of depleted stellar density in the cores of massive elliptical galaxies extends over the same radius as the gravitational sphere of influence of the central black holes, and interpret this as the dynamical imprint of the black holes.

preprint2016arXiv

A Comprehensive Archival Chandra Search for X-ray Emission from Ultracompact Dwarf Galaxies

We present the first comprehensive archival study of the X-ray properties of ultracompact dwarf (UCD) galaxies, with the goal of identifying weakly-accreting central black holes in UCDs. Our study spans 578 UCDs distributed across thirteen different host systems, including clusters, groups, fossil groups, and isolated galaxies. Of the 336 spectroscopically-confirmed UCDs with usable archival Chandra imaging observations, 21 are X-ray-detected. Imposing a completeness limit of $L_X>2\times10^{38}$ erg s$^{-1}$, the global X-ray detection fraction for the UCD population is $\sim3\%$. Of the 21 X-ray-detected UCDs, seven show evidence of long-term X-ray time variability on the order of months to years. X-ray-detected UCDs tend to be more compact than non-X-ray-detected UCDs, and we find tentative evidence that the X-ray detection fraction increases with surface luminosity density and global stellar velocity dispersion. The X-ray emission of UCDs is fully consistent with arising from a population of low-mass X-ray binaries (LMXBs). In fact, there are fewer X-ray sources than expected using a naive extrapolation from globular clusters. Invoking the fundamental plane of black hole activity for SUCD1 near the Sombrero galaxy, for which archival Jansky Very Large Array imaging at 5 GHz is publicly available, we set an upper limit on the mass of a hypothetical central black hole in that UCD to be $\lesssim10^5M_{\odot}$. While the majority of our sources are likely LMXBs, we cannot rule out central black holes in some UCDs based on X-rays alone, and so we address the utility of follow-up radio observations to find weakly-accreting central black holes.

preprint2016arXiv

Chandra X-ray and Hubble Space Telescope Imaging of Optically Selected Kiloparsec-Scale Binary Active Galactic Nuclei II: Host Galaxy Morphology and AGN Activity

Binary active galactic nuclei (AGNs) provide clues to how gas-rich mergers trigger and fuel AGNs and how supermassive black hole (SMBH) pairs evolve in a gas-rich environment. While significant effort has been invested in their identification, the detailed properties of binary AGNs and their host galaxies are still poorly constrained. In a companion paper, we examined the nature of ionizing sources in the double nuclei of four kpc-scale binary AGNs with redshifts between 0.1~0.2. Here, we present their host galaxy morphology based on F336W (U-band) and F105W (Y-band) images taken by the Wide Field Camera 3 (WFC3) onboard the Hubble Space Telescope. Our targets have double-peaked narrow emission lines and were confirmed to host binary AGNs with follow up observations. We find that kpc-scale binary AGNs occur in galaxy mergers with diverse morphological types. There are three major mergers with intermediate morphologies and a minor merger with a dominant disk component. We estimate the masses of the SMBHs from their host bulge stellar masses and obtain Eddington ratios for each AGN. Compared with a representative control sample drawn at the same redshift and stellar mass, the AGN luminosities and Eddington ratios of our binary AGNs are similar to those of single AGNs. The U-Y color maps indicate that clumpy star forming regions could significantly affect the X-ray detection of binary AGNs, e.g., the hardness ratio. Considering the weak X-ray emission in AGNs triggered in merger systems, we suggest that samples of X-ray selected AGNs may be biased against gas-rich mergers.

preprint2016arXiv

Discovery of extreme [OIII]5007A outflows in high-redshift red quasars

Black hole feedback is now a standard component of galaxy formation models. These models predict that the impact of black hole activity on its host galaxy likely peaked at z=2-3, the epoch of strongest star formation activity and black hole accretion activity in the Universe. We used XShooter on the Very Large Telescope to measure rest-frame optical spectra of four z~2.5 extremely red quasars with infrared luminosities ~10^47 erg/sec. We present the discovery of very broad (full width at half max= 2600-5000 km/sec), strongly blue-shifted (by up to 1500 km/sec) [OIII]5007A emission lines in these objects. In a large sample of obscured and red quasars, [OIII] kinematics are positively correlated with infrared luminosity, and the four objects in our sample are on the extreme end both in [OIII] kinematics and infrared luminosity. We estimate that ~3% of the bolometric luminosity in these objects is being converted into the kinetic power of the observed wind. Photo-ionization estimates suggest that the [OIII] emission might be extended on a few kpc scales, which would suggest that the extreme outflow is affecting the entire host galaxy of the quasar. These sources may be the signposts of the most extreme form of quasar feedback at the peak epoch of galaxy formation, and may represent an active "blow-out" phase of quasar evolution.

preprint2016arXiv

Extremely Red Quasars from SDSS, BOSS and WISE: Classification of Optical Spectra

Quasars with extremely red infrared-to-optical colours are an interesting population that can test ideas about quasar evolution as well as orientation, obscuration and geometric effects in the so-called AGN unified model. To identify such a population we match the quasar catalogues of the Sloan Digital Sky Survey (SDSS), the Baryon Oscillation Spectroscopic Survey (BOSS) to the Wide-Field Infrared Survey Explorer (WISE) to identify quasars with extremely high infrared-to-optical ratios. We identify 65 objects with r(AB)-W4(Vega)>14 mag (i.e., F_nu(22um)/F_nu(r) > ~1000). This sample spans a redshift range of 0.28<z<4.36 and has a bimodal distribution, with peaks at z~0.8 and z~2.5. It includes three z>2.6 objects that are detected in the W4-band but not W1 or W2 (i.e., W1W2-dropouts). The SDSS/BOSS spectra show that the majority of the objects are reddened Type 1 quasars, Type 2 quasars (both at low and high redshift) or objects with deep low-ionization broad absorption lines (BALs) that suppress the observed r-band flux. In addition, we identify a class of Type 1 permitted broad-emission line objects at z~2-3 which are characterized by emission line rest-frame equivalent widths (REWs) of >~150Ang , much larger than those of typical quasars. In particular, 55% (45%) of the non-BAL Type 1s with measurable CIV in our sample have REW(CIV) > 100 (150)Ang, compared to only 5.8% (1.3%) for non-BAL quasars in BOSS. These objects often also have unusual line ratios, such as very high NV/Ly-alpha ratios. These large REWs might be caused by suppressed continuum emission analogous to Type 2 quasars; however, there is no obvious mechanism in standard Unified Models to suppress the continuum without also obscuring the broad emission lines.

preprint2016arXiv

Megamaser Disks Reveal a Broad Distribution of Black Hole Mass in Spiral Galaxies

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.

preprint2016arXiv

Prime Focus Spectrograph (PFS) for the Subaru Telescope: Overview, recent progress, and future perspectives

PFS (Prime Focus Spectrograph), a next generation facility instrument on the 8.2-meter Subaru Telescope, is a very wide-field, massively multiplexed, optical and near-infrared spectrograph. Exploiting the Subaru prime focus, 2394 reconfigurable fibers will be distributed over the 1.3 deg field of view. The spectrograph has been designed with 3 arms of blue, red, and near-infrared cameras to simultaneously observe spectra from 380nm to 1260nm in one exposure at a resolution of ~1.6-2.7A. An international collaboration is developing this instrument under the initiative of Kavli IPMU. The project is now going into the construction phase aiming at undertaking system integration in 2017-2018 and subsequently carrying out engineering operations in 2018-2019. This article gives an overview of the instrument, current project status and future paths forward.

preprint2016arXiv

Sensitive Radio Survey of Obscured Quasar Candidates

We study the radio properties of moderately obscured quasars over a range of redshifts to understand the role of radio activity in accretion using the Jansky Very Large Array (JVLA) at 6.0GHz and 1.4GHz. Our z~2.5 sample consists of optically-selected obscured quasar candidates, all of which are radio-quiet, with typical radio luminosities of $νL_ν$[1.4 GHz] < $10^{40}$ erg s$^{-1}$. Only a single source is individually detected in our deep (rms~10 $μ$Jy) exposures. This population would not be identified by radio-based selection methods used for distinguishing dusty star-forming galaxies and obscured active nuclei. In our pilot A-array study of z~0.5 radio-quiet quasars, we spatially resolve four of five objects on scales ~ 5 kpc and find they have steep spectral indices. Therefore, radio emission in these sources could be due to jet-driven or radiatively driven bubbles interacting with interstellar material on the scale of the host galaxy. Finally, we also study the population of ~ 200 faint (~40 $μ$Jy - 40 mJy) radio sources observed over ~ 120 arcmin$^2$ of our data. 60% of these detections are matched in the SDSS and/or WISE and are, in roughly equal shares, active nuclei at a broad range of redshifts, passive galaxies with no other signs of nuclear activity and IR-bright but optically faint sources. Spectroscopically or photometrically confirmed star-forming galaxies constitute only a small minority of the matches. Such sensitive radio surveys allow us to address important questions of AGN evolution and evaluate the AGN contribution to the radio-quiet sky.

preprint2016arXiv

Spatially Offset Active Galactic Nuclei I: Selection and Spectroscopic Properties

We present a sample of 18 optically-selected and X-ray detected spatially offset active galactic nuclei (AGN) from the Sloan Digital Sky Survey (SDSS). In 9 systems, the X-ray AGN is spatially offset from the galactic stellar core that is located within the 3'' diameter SDSS spectroscopic fiber. In 11 systems, the X-ray AGN is spatially offset from a stellar core that is located outside the fiber, with an overlap of 2. To build the sample, we cross-matched Type II AGN selected from the SDSS galaxy catalogue with archival Chandra imaging and employed our custom astrometric and registration procedure. The projected angular (physical) offsets span a range of 0."6 (0.8 kpc) to 17."4 (19.4 kpc), with a median value of 2."7 (4.6 kpc). The offset nature of an AGN is an unambiguous signature of a galaxy merger, and these systems can be used to study the properties of AGN in galaxy mergers without the biases introduced by morphological merger selection techniques. In this paper (Paper I), we use our sample to assess the kinematics of AGN photoionized gas in galaxy mergers. We find that spectroscopic offset AGN selection may be up to 89% incomplete due to small projected velocity offsets. We also find that the magnitude of the velocity offsets are generally larger than expected if our spatial selection introduces a bias toward face-on orbits, suggesting the presence of complex kinematics in the emission line gas of AGN in galaxy mergers.

preprint2016arXiv

Suppressing star formation in quiescent galaxies with supermassive black hole winds

Quiescent galaxies with little or no ongoing star formation dominate the galaxy population above $M_{*}\sim 2 \times 10^{10}~M_{\odot}$, where their numbers have increased by a factor of $\sim25$ since $z\sim2$. Once star formation is initially shut down, perhaps during the quasar phase of rapid accretion onto a supermassive black hole, an unknown mechanism must remove or heat subsequently accreted gas from stellar mass loss or mergers that would otherwise cool to form stars. Energy output from a black hole accreting at a low rate has been proposed, but observational evidence for this in the form of expanding hot gas shells is indirect and limited to radio galaxies at the centers of clusters, which are too rare to explain the vast majority of the quiescent population. Here we report bisymmetric emission features co-aligned with strong ionized gas velocity gradients from which we infer the presence of centrally-driven winds in typical quiescent galaxies that host low-luminosity active nuclei. These galaxies are surprisingly common, accounting for as much as $10\%$ of the population at $M_* \sim 2 \times 10^{10}~ M_{\odot}$. In a prototypical example, we calculate that the energy input from the galaxy's low-level active nucleus is capable of driving the observed wind, which contains sufficient mechanical energy to heat ambient, cooler gas (also detected) and thereby suppress star formation.

preprint2016arXiv

The Black Hole - Bulge Mass Relation in Megamaser Host Galaxies

We present HST images for nine megamaser disk galaxies with the primary goal of studying photometric BH-galaxy scaling relations. The megamaser disks provide the highest-precision extragalactic BH mass measurements, while our high-resolution HST imaging affords us the opportunity to decompose the complex nuclei of their late-type hosts in detail. Based on the morphologies and shapes of the galaxy nuclei, we argue that most of these galaxies' central regions contain secularly evolving components (pseudo-bulges), and in many cases we photometrically identify co-existing "classical" bulge components as well. Using these decompositions, we draw the following conclusions: (1) The megamaser BH masses span two orders of magnitude ($10^6$ -- $10^8 M_\odot$) while the stellar mass of their spiral host galaxies are all $\sim 10^{11} M_\odot$ within a factor of three; (2) the BH masses at a given bulge mass or total stellar mass in the megamaser host spiral galaxies tend to be lower than expected, when compared to an extrapolation of the BH-bulge relation based on early-type galaxies; (3) the observed large intrinsic scatter of BH masses in the megamaser host galaxies raises the question of whether scaling relations exist in spiral galaxies.

preprint2016arXiv

The MASSIVE Survey - V. Spatially-Resolved Stellar Angular Momentum, Velocity Dispersion, and Higher Moments of the 41 Most Massive Local Early-Type Galaxies

We present spatially-resolved two-dimensional stellar kinematics for the 41 most massive early-type galaxies (MK <~ -25.7 mag, stellar mass M* >~ 10^11.8 Msun) of the volume-limited (D < 108 Mpc) MASSIVE survey. For each galaxy, we obtain high-quality spectra in the wavelength range of 3650 to 5850 A from the 246-fiber Mitchell integral-field spectrograph (IFS) at McDonald Observatory, covering a 107" x 107" field of view (often reaching 2 to 3 effective radii). We measure the 2-D spatial distribution of each galaxy's angular momentum (lambda and fast or slow rotator status), velocity dispersion (sigma), and higher-order non-Gaussian velocity features (Gauss-Hermite moments h3 to h6). Our sample contains a high fraction (~80% ) of slow and non-rotators with lambda <~ 0.2. When combined with the lower-mass ETGs in the ATLAS3D survey, we find the fraction of slow-rotators to increase dramatically with galaxy mass, reaching ~50% at MK ~ -25.5 mag and ~90% at MK <~ -26 mag. All of our fast rotators show a clear anti-correlation between h3 and V/sigma, and the slope of the anti-correlation is steeper in more round galaxies. The radial profiles of sigma show a clear luminosity and environmental dependence: the 12 most luminous galaxies in our sample (MK <~ -26 mag) are all brightest cluster/group galaxies (except NGC 4874) and all have rising or nearly flat sigma profiles, whereas five of the seven "isolated" galaxies are all fainter than MK = -25.8 mag and have falling sigma. All of our galaxies have positive average h4; the most luminous galaxies have average h4 ~ 0.05 while less luminous galaxies have a range of values between 0 and 0.05. Most of our galaxies show positive radial gradients in h4, and those galaxies also tend to have rising sigma profiles. We discuss the implications for the relationship among dynamical mass, sigma, h4, and velocity anisotropy for these massive galaxies.

preprint2016arXiv

The MASSIVE Survey IV.: The X-ray halos of the most massive early-type galaxies in the nearby Universe

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.

preprint2016arXiv

The Sloan Digital Sky Survey Reverberation Mapping Project: Velocity Shifts of Quasar Emission Lines

Quasar emission lines are often shifted from the systemic velocity due to various dynamical and radiative processes in the line-emitting region. The level of these velocity shifts depends both on the line species and on quasar properties. We study velocity shifts for the line peaks of various narrow and broad quasar emission lines relative to systemic using a sample of 849 quasars from the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project. The coadded (from 32 epochs) spectra of individual quasars have sufficient signal-to-noise ratio (SNR) to measure stellar absorption lines to provide reliable systemic velocity estimates, as well as weak narrow emission lines. The sample also covers a large dynamic range in quasar luminosity (~2 dex), allowing us to explore potential luminosity dependence of the velocity shifts. We derive average line peak velocity shifts as a function of quasar luminosity for different lines, and quantify their intrinsic scatter. We further quantify how well the peak velocity can be measured for various lines as a function of continuum SNR, and demonstrate there is no systematic bias in the line peak measurements when the spectral quality is degraded to as low as SNR~3 per SDSS pixel. Based on the observed line shifts, we provide empirical guidelines on redshift estimation from [OII]3728, [OIII]5008, [NeV]3426, MgII, CIII], HeII1640, broad Hbeta, CIV, and SiIV, which are calibrated to provide unbiased systemic redshifts in the mean, but with increasing intrinsic uncertainties of 46, 56, 119, 205, 233, 242, 400, 415, and 477 km/s, in addition to the measurement uncertainties. These more realistic redshift uncertainties are generally much larger than the formal uncertainties reported by the redshift pipelines for spectroscopic quasar surveys, and demonstrate the infeasibility of measuring quasar redshifts to better than ~200 km/s with only broad lines.

preprint2016arXiv

Towards an Understanding of Changing-Look Quasars: An Archival Spectroscopic Search in SDSS

The uncertain origin of the recently-discovered `changing-looking' quasar phenomenon -- in which a luminous quasar dims significantly to a quiescent state in repeat spectroscopy over ~10 year timescales -- may present unexpected challenges to our understanding of quasar accretion. To better understand this phenomenon, we take a first step to building a sample of changing-look quasars with a systematic but simple archival search for these objects in the Sloan Digital Sky Survey Data Release 12. By leveraging the >10 year baselines for objects with repeat spectroscopy, we uncover two new changing-look quasars, and a third discovered previously. Decomposition of the multi-epoch spectra and analysis of the broad emission lines suggest that the quasar accretion disk emission dims due to rapidly decreasing accretion rates (by factors of >2.5), while disfavoring changes in intrinsic dust extinction for the two objects where these analyses are possible. Broad emission line energetics also support intrinsic dimming of quasar emission as the origin for this phenomenon rather than transient tidal disruption events or supernovae. Although our search criteria included quasars at all redshifts and transitions from either quasar-like to galaxy-like states or the reverse, all of the clear cases of changing-look quasars discovered were at relatively low-redshift (z ~ 0.2 - 0.3) and only exhibit quasar-like to galaxy-like transitions.

preprint2016arXiv

Towards Precision Supermassive Black Hole Masses using Megamaser Disks

Megamaser disks provide the most precise and accurate extragalactic supermassive black hole masses. Here we describe a search for megamasers in nearby galaxies using the Green Bank Telescope (GBT). We focus on galaxies where we believe that we can resolve the gravitational sphere of influence of the black hole and derive a stellar or gas dynamical measurement with optical or NIR observations. Since there are only a handful of super massive black holes (SMBH) that have direct black hole mass measurements from more than one method, even a single galaxy with a megamaser disk and a stellar dynamical black hole mass would provide necessary checks on the stellar dynamical methods. We targeted 87 objects from the Hobby-Eberly Telescope Massive Galaxy Survey, and detected no new maser disks. Most of the targeted objects are elliptical galaxies with typical stellar velocity dispersions of 250 km/s and distances within 130 Mpc. We discuss the implications of our non-detections, whether they imply a threshold X-ray luminosity required for masing, or possibly reflect the difficulty of maintaining a masing disk around much more massive (>10^8 Msun) black holes at low Eddington ratio. Given the power of maser disks at probing black hole accretion and demographics, we suggest that future maser searches should endeavour to remove remaining sample biases, in order to sort out the importance of these covariant effects.

preprint2015arXiv

Merger-driven Fueling of Active Galactic Nuclei: Six Dual and Offset Active Galactic Nuclei Discovered with Chandra and Hubble Space Telescope Observations

Dual active galactic nuclei (AGNs) and offset AGNs are kpc-scale separation supermassive black holes pairs created during galaxy mergers, where both or one of the black holes are AGNs, respectively. These dual and offset AGNs are valuable probes of the link between mergers and AGNs but are challenging to identify. Here we present Chandra/ACIS observations of 12 optically-selected dual AGN candidates at z < 0.34, where we use the X-rays to identify AGNs. We also present HST/WFC3 observations of 10 of these candidates, which reveal any stellar bulges accompanying the AGNs. We discover a dual AGN system with separation of 2.2 kpc, where the two stellar bulges have coincident [O III] and X-ray sources. This system is an extremely minor merger (460:1) that may include a dwarf galaxy hosting an intermediate mass black hole. We also find six single AGNs, and five systems that are either dual or offset AGNs with separations < 10 kpc. Four of the six dual AGNs and dual/offset AGNs are in ongoing major mergers, and these AGNs are 10 times more luminous, on average, than the single AGNs in our sample. This hints that major mergers may preferentially trigger higher luminosity AGNs. Further, we find that confirmed dual AGNs have hard X-ray luminosities that are half of those of single AGNs at fixed [O III] luminosity, on average. This could be explained by high densities of gas funneled to galaxy centers during mergers, and emphasizes the need for deeper X-ray observations of dual AGN candidates.

preprint2015arXiv

Now You See It, Now You Don't: The Disappearing Central Engine of the Quasar J1011+5442

We report the discovery of a new "changing-look" quasar, SDSS J101152.98+544206.4, through repeat spectroscopy from the Time Domain Spectroscopic Survey. This is an addition to a small but growing set of quasars whose blue continua and broad optical emission lines have been observed to decline by a large factor on a time scale of approximately a decade. The 5100 Angstrom monochromatic continuum luminosity of this quasar drops by a factor of > 9.8 in a rest-frame time interval of < 9.7 years, while the broad H-alpha luminosity drops by a factor of 55 in the same amount of time. The width of the broad H-alpha line increases in the dim state such that the black hole mass derived from the appropriate single-epoch scaling relation agrees between the two epochs within a factor of 3. The fluxes of the narrow emission lines do not appear to change between epochs. The light curve obtained by the Catalina Sky Survey suggests that the transition occurs within a rest-frame time interval of approximately 500 days. We examine three possible mechanisms for this transition suggested in the recent literature. An abrupt change in the reddening towards the central engine is disfavored by the substantial difference between the timescale to obscure the central engine and the observed timescale of the transition. A decaying tidal disruption flare is consistent with the decay rate of the light curve but not with the prolonged bright state preceding the decay, nor can this scenario provide the power required by the luminosities of the emission lines. An abrupt drop in the accretion rate onto the supermassive black hole appears to be the most plausible explanation for the rapid dimming.

preprint2015arXiv

On R-W1 as A Diagnostic to Discover Obscured Active Galactic Nuclei in Wide-Area X-ray Surveys

Capitalizing on the all-sky coverage of {\it WISE}, and the 35\% and 50\% sky coverage from SDSS and Pan-STARRS, respectively, we explore the efficacy of $m_{R}$ (optical) - $m_{3.4 μm}$ (mid-infrared), hereafter $R-W1$, as a color diagnostic to identify obscured supermassive black hole accretion in wide-area X-ray surveys. We use the $\sim$16.5 deg$^2$ Stripe 82 X-ray survey data as a test-bed to compare $R-W1$ with $R-K$, an oft-used obscured AGN selection criterion, and examine where different classes of objects lie in this parameter space. Most stars follow a well-defined path in $R-K$ vs. $R-W1$ space. We demonstrate that optically normal galaxies hosting X-ray AGN at redshifts $0.5<z<1$ can be recovered with an $R-W1>4$ color-cut, while they typically are not selected as AGN based on their $W1-W2$ colors. Additionally, different observed X-ray luminosity bins favor different regions in $R-W1$ parameter space: moderate luminosity AGN ($10^{43}$ erg s$^{-1} < L_{\rm 0.5-10 keV} < 10^{44}$ erg s$^{-1}$) tend to have red colors while the highest luminosity AGN ($L_{\rm 0.5-10 keV} > 10^{45}$ erg s$^{-1}$) have bluer colors; higher spectroscopic completeness of the Stripe 82X sample is needed to determine whether this is a selection effect or an intrinsic property. Finally, we parameterize X-ray obscuration of Stripe 82X AGN by calculating their hardness ratios (HRs) and find no clear trends between HR and optical reddening. Our results will help inform best-effort practices in following-up obscured AGN candidates in current and future wide-area, shallow X-ray surveys, including the all-sky {\it eROSITA} mission.

preprint2015arXiv

Similarity of ionized gas nebulae around unobscured and obscured quasars

Quasar feedback is suspected to play a key role in the evolution of massive galaxies, by removing or reheating gas in quasar host galaxies and thus limiting the amount of star formation. In this paper we continue our investigation of quasar-driven winds on galaxy-wide scales. We conduct Gemini Integral Field Unit spectroscopy of a sample of luminous unobscured (type 1) quasars, to determine the morphology and kinematics of ionized gas around these objects, predominantly via observations of the [O III]5007 emission line. We find that ionized gas nebulae extend out to ~13 kpc from the quasar, that they are smooth and round, and that their kinematics are inconsistent with gas in dynamical equilibrium with the host galaxy. The observed morphological and kinematic properties are strikingly similar to those of ionized gas around obscured (type 2) quasars with matched [O III] luminosity, with marginal evidence that nebulae around unobscured quasars are slightly more compact. Therefore in samples of obscured and unobscured quasars carefully matched in [O III] luminosity we find support for the standard geometry-based unification model of active galactic nuclei, in that the intrinsic properties of quasars, of their hosts and of their ionized gas appear to be very similar. Given the apparent ubiquity of extended ionized regions, we are forced to conclude that either the quasar is at least partially illuminating pre-existing gas or that both samples of quasars are seen during advanced stages of quasar feedback. In the latter case, we may be biased by our [O III]-based selection against quasars in the early "blow-out" phase, for example due to dust obscuration.

preprint2015arXiv

Star formation in quasar hosts and the origin of radio emission in radio-quiet quasars

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

preprint2015arXiv

The MASSIVE Survey II: Stellar Population Trends Out to Large Radius in Massive Early Type Galaxies

We examine stellar population gradients in ~100 massive early type galaxies spanning 180 < sigma* < 370 km/s and M_K of -22.5 to -26.5 mag, observed as part of the MASSIVE survey (Ma et al. 2014). Using integral-field spectroscopy from the Mitchell Spectrograph on the 2.7m telescope at McDonald Observatory, we create stacked spectra as a function of radius for galaxies binned by their stellar velocity dispersion, stellar mass, and group richness. With excellent sampling at the highest stellar mass, we examine radial trends in stellar population properties extending to beyond twice the effective radius (~2.5 R_e). Specifically, we examine trends in age, metallicity, and abundance ratios of Mg, C, N, and Ca, and discuss the implications for star formation histories and elemental yields. At a fixed physical radius of 3-6 kpc (the likely size of the galaxy cores formed at high redshift) stellar age and [alpha/Fe] increase with increasing sigma* and depend only weakly on stellar mass, as we might expect if denser galaxies form their central cores earlier and faster. If we instead focus on 1-1.5 R_e, the trends in abundance and abundance ratio are washed out, as might be expected if the stars at large radius were accreted by smaller galaxies. Finally, we show that when controlling for \sigmastar, there are only very subtle differences in stellar population properties or gradients as a function of group richness; even at large radius internal properties matter more than environment in determining star formation history.

preprint2015arXiv

The Origin of Double-Peaked Narrow Lines in Active Galactic Nuclei I: Very Large Array Detections of Dual AGNs and AGN Outflows

We have examined a subset of 18 active galactic nuclei (AGNs) drawn from a sample of 81 galaxies that possess double-peaked narrow optical emission line spectra in the Sloan Digital Sky Survey, have two optical AGN emission components separated by >0.2", and are detected in the Faint Images of the Radio Sky at Twenty-centimeters survey. Without follow-up observations, the sources of the double-peaked narrow emission lines are uncertain, and may be produced by kpc-scale separation dual active supermassive black holes, AGN outflows, or disk rotation. In this work, we propose a new methodology to characterize double-peaked narrow emission-line galaxies based on optical long-slit spectroscopy and high resolution multi-band Very Large Array observations. The nature of the radio emission in the sample galaxies is varied. Of the 18 galaxies, we detect two compact flat-spectrum radio cores with projected spatial separations on the sky between 0.6-1.6 kpc in three galaxies: J1023+3243, J1158+3231, and J1623+0808. The two radio sources are spatially coincident with the two optical components of ionized gas with AGN-like line ratios, which confirms the presence of dual AGNs in these three galaxies. Dual AGNs account for only ~15% (3/18) of the double-peaked AGNs in our sample. Gas kinematics produce ~75% (13/18) of the double-peaked narrow emission lines, distributed in the following way: 7 AGN wind-driven outflows, 5 radio-jet driven outflows, and one rotating narrow-line region. The remaining ~10% (2/18) are ambiguous cases. Our method demonstrates the power of spatially resolved spectroscopy and high resolution radio observations for the identification of AGN outflows and AGN pairs with angular separations as small as 0.18".

preprint2015arXiv

The Sloan Digital Sky Survey Reverberation Mapping Project: First Broad-line Hbeta and MgII Lags at z>~0.3 from six-Month Spectroscopy

Reverberation mapping (RM) measurements of broad-line region (BLR) lags in z>0.3 quasars are important for directly measuring black hole masses in these distant objects, but so far there have been limited attempts and success given the practical difficulties of RM in this regime. Here we report preliminary results of 15 BLR lag measurements from the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project, a dedicated RM program with multi-object spectroscopy designed for RM over a wide redshift range. The lags are based on the 2014 spectroscopic light curves alone (32 epochs over 6 months) and focus on the Hbeta and MgII broad lines in the 100 lowest-redshift (z<0.8) quasars included in SDSS-RM; they represent a small subset of the lags that SDSS-RM (including 849 quasars to z~4.5) is expected to deliver. The reported preliminary lag measurements are for intermediate-luminosity quasars at 0.3<~z<0.8, including 9 Hbeta lags and 6 MgII lags, for the first time extending RM results to this redshift-luminosity regime and providing direct quasar black hole mass estimates over ~ half of cosmic time. The MgII lags also increase the number of known MgII lags by several-fold, and start to explore the utility of MgII for RM at high redshift. The location of these new lags at higher redshifts on the observed BLR size-luminosity relationship is statistically consistent with previous Hbeta results at z<0.3. However, an independent constraint on the relationship slope at z>0.3 is not yet possible due to the limitations in our current sample. Our results demonstrate the general feasibility and potential of multi-object RM for z>0.3 quasars.

preprint2015arXiv

The Sloan Digital Sky Survey Reverberation Mapping Project: No Evidence for Evolution in the M-sigma Relation to z~1

We present host stellar velocity dispersion measurements for a sample of 88 broad-line quasars at 0.1<z<1 (46 at z>0.6) from the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project. High signal-to-noise ratio coadded spectra (average S/N~30 per 69 km/s pixel) from SDSS-RM allowed decomposition of the host and quasar spectra, and measurement of the host stellar velocity dispersions and black hole (BH) masses using the single-epoch (SE) virial method. The large sample size and dynamic range in luminosity (L5100=10^(43.2-44.7) erg/s) lead to the first clear detection of a correlation between SE virial BH mass and host stellar velocity dispersion far beyond the local universe. However, the observed correlation is significantly flatter than the local relation, suggesting that there are selection biases in high-z luminosity-threshold quasar samples for such studies. Our uniform sample and analysis enable an investigation of the redshift evolution of the M-sigma relation free of caveats by comparing different samples/analyses at disjoint redshifts. We do not observe evolution of the M-sigma relation in our sample, up to z~1, but there is an indication that the relation flattens towards higher redshifts. Coupled with the increasing threshold luminosity with redshift in our sample, this again suggests certain selection biases are at work, and simple simulations demonstrate that a constant M-sigma relation is favored to z~1. Our results highlight the scientific potential of deep coadded spectroscopy from quasar monitoring programs, and offer a new path to probe the co-evolution of BHs and galaxies at earlier times.

preprint2015arXiv

The Sloan Digital Sky Survey Reverberation Mapping Project: Post-Starburst Signatures in Quasar Host Galaxies at z < 1

Quasar host galaxies are key for understanding the relation between galaxies and the supermassive black holes (SMBHs) at their centers. We present a study of 191 broad-line quasars and their host galaxies at z < 1, using high signal-to-noise ratio (SNR) spectra produced by the Sloan Digital Sky Survey Reverberation Mapping project. Clear detection of stellar absorption lines allows a reliable decomposition of the observed spectra into nuclear and host components, using spectral models of quasar and stellar radiations as well as emission lines from the interstellar medium. We estimate age, mass (M*), and velocity dispersion (sigma*) of the host stars, the star formation rate (SFR), quasar luminosity, and SMBH mass (Mbh), for each object. The quasars are preferentially hosted by massive galaxies with M* ~ 10^{11} Msun characterized by stellar ages around a billion years, which coincides with the transition phase of normal galaxies from the blue cloud to the red sequence. The host galaxies have relatively low SFRs and fall below the main sequence of star-forming galaxies at similar redshifts. These facts suggest that the hosts have experienced an episode of major star formation sometime in the past billion years, which was subsequently quenched or suppressed. The derived Mbh - sigma* and Mbh - M* relations agree with our past measurements and are consistent with no evolution from the local Universe. The present analysis demonstrates that reliable measurements of stellar properties of quasar host galaxies are possible with high-SNR fiber spectra, which will be acquired in large numbers with future powerful instruments such as the Subaru Prime Focus Spectrograph.

preprint2015arXiv

The Structure of Nuclear Star Clusters in Nearby Late-type Spiral Galaxies from Hubble Space Telescope Wide Field Camera 3 Imaging

We obtained Hubble Space Telescope/Wide Field Camera 3 imaging of a sample of ten of the nearest and brightest nuclear clusters residing in late-type spiral galaxies, in seven bands that span the near-ultraviolet to the near-infrared. Structural properties of the clusters were measured by fitting two-dimensional surface brightness profiles to the images using GALFIT. The clusters exhibit a wide range of structural properties. For six of the ten clusters in our sample, we find changes in the effective radius with wavelength, suggesting radially varying stellar populations. In four of the objects, the effective radius increases with wavelength, indicating the presence of a younger population which is more concentrated than the bulk of the stars in the cluster. However, we find a general decrease in effective radius with wavelength in two of the objects in our sample, which may indicate extended, circumnuclear star formation. We also find a general trend of increasing roundness of the clusters at longer wavelengths, as well as a correlation between the axis ratios of the NCs and their host galaxies. These observations indicate that blue disks aligned with the host galaxy plane are a common feature of nuclear clusters in late-type galaxies, but are difficult to detect in galaxies that are close to face-on. In color-color diagrams spanning the near-UV through the near-IR, most of the clusters lie far from single-burst evolutionary tracks, showing evidence for multi-age populations. Most of the clusters have integrated colors consistent with a mix of an old population (> 1 Gyr) and a young population (~100-300 Myr). The wide wavelength coverage of our data provides a sensitivity to populations with a mix of ages that would not be possible to achieve with imaging in optical bands only.

preprint2015arXiv

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

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

preprint2014arXiv

ALMA Observations of a Candidate Molecular Outflow in an Obscured Quasar

We present Atacama Large Millimeter/Submillimeter Array (ALMA) CO (1-0) and CO (3-2) observations of SDSS J135646.10+102609.0, an obscured quasar and ultra-luminous infrared galaxy (ULIRG) with two merging nuclei and a known 20-kpc-scale ionized outflow. The total molecular gas mass is M_{mol} ~ 9^{+19}_{-6} x 10^8 Msun, mostly distributed in a compact rotating disk at the primary nucleus (M_{mol} ~ 3 x 10^8 Msun) and an extended tidal arm (M_{mol} ~ 5 x 10^8 Msun). The tidal arm is one of the most massive molecular tidal features known; we suggest that it is due to the lower chance of shock dissociation in this elliptical/disk galaxy merger. In the spatially resolved CO (3-2) data, we find a compact (r ~ 0.3 kpc) high velocity (v ~ 500 km/s) red-shifted feature in addition to the rotation at the N nucleus. We propose a molecular outflow as the most likely explanation for the high velocity gas. The outflowing mass of M_{mol} ~ 7 x 10^7 Msun and the short dynamical time of t_{dyn} ~ 0.6 Myr yield a very high outflow rate of \dot{M}_{mol} ~ 350 Msun/yr and can deplete the gas in a million years. We find a low star formation rate (< 16 Msun/yr from the molecular content and < 21 Msun/yr from the far-infrared spectral energy distribution decomposition) that is inadequate to supply the kinetic luminosity of the outflow (\dot{E} ~ 3 x 10^43 erg/s). Therefore, the active galactic nucleus, with a bolometric luminosity of 10^46 erg/s, likely powers the outflow. The momentum boost rate of the outflow (\dot{p}/(Lbol/c) ~ 3) is lower than typical molecular outflows associated with AGN, which may be related to its compactness. The molecular and ionized outflows are likely two distinct bursts induced by episodic AGN activity that varies on a time scale of 10^7 yr.

preprint2014arXiv

Circumnuclear Molecular Gas in Megamaser Disk Galaxies NGC 4388 and NGC 1194

We explore the warm molecular and ionized gas in the centers of two megamaser disk galaxies using K-band spectroscopy. Our ultimate goal is to determine how gas is funneled onto the accretion disk, here traced by megamaser spots on sub-pc scales. We present NIR IFU data with a resolution of ~50 pc for two galaxies: NGC 4388 with VLT/SINFONI and NGC 1194 with Keck/OSIRIS+AO. The high spatial resolution and rich spectral diagnostics allow us to study both the stellar and gas kinematics as well as gas excitation on scales only an order of magnitude larger than the maser disk. We find a drop in the stellar velocity dispersion in the inner ~100 pc of NGC 4388, a common signature of a dynamically cold central component seen in many active nuclei. We also see evidence for non-circular gas motions in the molecular hydrogen on similar scales, with the gas kinematics on 100-pc scales aligned with the megamaser disk. In contrast, the high ionization lines and Br-gamma trace outflow along the 100 pc-scale jet. In NGC 1194, the continuum from the accreting black hole is very strong, making it difficult to measure robust two-dimensional kinematics, but the spatial distribution and line ratios of the molecular hydrogen and Br-gamma have consistent properties between the two galaxies.

preprint2014arXiv

Extended X-ray Emission From a Quasar-Driven Superbubble

We present observations of extended, 20-kpc scale soft X-ray gas around a luminous obscured quasar hosted by an ultra-luminous infrared galaxy caught in the midst of a major merger. The extended X-ray emission is well fit as a thermal gas with a temperature of kT ~ 280 eV and a luminosity of L_X ~ 10^42 erg/s and is spatially coincident with a known ionized gas outflow. Based on the X-ray luminosity, a factor of ~10 fainter than the [OIII] emission, we conclude that the X-ray emission is either dominated by photoionization, or by shocked emission from cloud surfaces in a hot quasar-driven wind.

preprint2014arXiv

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érsic or Voigt profiles convolved with the seeing. These data, which span a full order of magnitude in IR luminosity ($\log{(L_{8 μ\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.

preprint2014arXiv

Offset Active Galactic Nuclei as Tracers of Galaxy Mergers and Supermassive Black Hole Growth

Offset active galactic nuclei (AGNs) are AGNs that are in ongoing galaxy mergers, which produce kinematic offsets in the AGNs relative to their host galaxies. Offset AGNs are also close relatives of dual AGNs. We conduct a systematic search for offset AGNs in the Sloan Digital Sky Survey, by selecting AGN emission lines that exhibit statistically significant line-of-sight velocity offsets relative to systemic. From a parent sample of 18314 Type 2 AGNs at z<0.21, we identify 351 offset AGN candidates with velocity offsets of 50 km/s < |v| < 410 km/s. When we account for projection effects in the observed velocities, we estimate that 4% - 8% of AGNs are offset AGNs. We designed our selection criteria to bypass velocity offsets produced by rotating gas disks, AGN outflows, and gravitational recoil of supermassive black holes, but follow-up observations are still required to confirm our candidates as offset AGNs. We find that the fraction of AGNs that are offset candidates increases with AGN bolometric luminosity, from 0.7% to 6% over the luminosity range 43 < log(L_bol) [erg/s] < 46. If these candidates are shown to be bona fide offset AGNs, then this would be direct observational evidence that galaxy mergers preferentially trigger high-luminosity AGNs. Finally, we find that the fraction of AGNs that are offset AGN candidates increases from 1.9% at z=0.1 to 32% at z=0.7, in step with the growth in the galaxy merger fraction over the same redshift range.

preprint2014arXiv

Quasar feedback and the origin of radio emission in radio-quiet quasars

We conduct kinematic analysis of the SDSS spectra of 568 obscured luminous quasars, with the emphasis on the kinematic structure of the [OIII]5007 emission line. [OIII] emission tends to show blueshifts and blue excess, which indicates that at least part of the narrow-line gas is undergoing an organized outflow. The velocity width containing 90% of line power ranges from 370 to 4780 km/sec, suggesting outflow velocities up to 2000 km/sec. The velocity width of the [OIII] emission is positively correlated with the radio luminosity among the radio-quiet quasars. We propose that radio emission in radio-quiet quasars is due to relativistic particles accelerated in the shocks within the quasar-driven outflows; star formation in quasar hosts is insufficient to explain the observed radio emission. The median radio luminosity of the sample of nu L_nu[1.4GHz] = 10^40 erg/sec suggests a median kinetic luminosity of the quasar-driven wind of L_wind=3x10^44 erg/sec, or about 4% of the estimated median bolometric luminosity L_bol=8x10^45 erg/sec. Furthermore, the velocity width of [OIII] is positively correlated with mid-infrared luminosity, which suggests that outflows are ultimately driven by the radiative output of the quasar. As the outflow velocity increases, some emission lines characteristic of shocks in quasi-neutral medium increase as well, which we take as further evidence of quasar-driven winds propagating into the interstellar medium of the host galaxy. None of the kinematic components show correlations with the stellar velocity dispersions of the host galaxies, so there is no evidence that any of the gas in the narrow-line region of quasars is in dynamical equilibrium with the host galaxy. Quasar feedback appears to operate above the threshold luminosity of L_bol=3x10^45 erg/sec.

preprint2014arXiv

The MASSIVE Survey - I. A Volume-Limited Integral-Field Spectroscopic Study of the Most Massive Early-Type Galaxies within 108 Mpc

Massive early-type galaxies represent the modern-day remnants of the earliest major star formation episodes in the history of the universe. These galaxies are central to our understanding of the evolution of cosmic structure, stellar populations, and supermassive black holes, but the details of their complex formation histories remain uncertain. To address this situation, we have initiated the MASSIVE Survey, a volume-limited, multi-wavelength, integral-field spectroscopic (IFS) and photometric survey of the structure and dynamics of the ~100 most massive early-type galaxies within a distance of 108 Mpc. This survey probes a stellar mass range M* > 10^{11.5} Msun and diverse galaxy environments that have not been systematically studied to date. Our wide-field IFS data cover about two effective radii of individual galaxies, and for a subset of them, we are acquiring additional IFS observations on sub-arcsecond scales with adaptive optics. We are also acquiring deep K-band imaging to trace the extended halos of the galaxies and measure accurate total magnitudes. Dynamical orbit modeling of the combined data will allow us to simultaneously determine the stellar, black hole, and dark matter halo masses. The primary goals of the project are to constrain the black hole scaling relations at high masses, investigate systematically the stellar initial mass function and dark matter distribution in massive galaxies, and probe the late-time assembly of ellipticals through stellar population and kinematical gradients. In this paper, we describe the MASSIVE sample selection, discuss the distinct demographics and structural and environmental properties of the selected galaxies, and provide an overview of our basic observational program, science goals and early survey results.

preprint2014arXiv

The Sloan Digital Sky Survey Reverberation Mapping Project: Technical Overview

The Sloan Digital Sky Survey Reverberation Mapping project (SDSS-RM) is a dedicated multi-object RM experiment that has spectroscopically monitored a sample of 849 broad-line quasars in a single 7 deg$^2$ field with the SDSS-III BOSS spectrograph. The RM quasar sample is flux-limited to i_psf=21.7 mag, and covers a redshift range of 0.1<z<4.5. Optical spectroscopy was performed during 2014 Jan-Jul dark/grey time, with an average cadence of ~4 days, totaling more than 30 epochs. Supporting photometric monitoring in the g and i bands was conducted at multiple facilities including the CFHT and the Steward Observatory Bok telescopes in 2014, with a cadence of ~2 days and covering all lunar phases. The RM field (RA, DEC=14:14:49.00, +53:05:00.0) lies within the CFHT-LS W3 field, and coincides with the Pan-STARRS 1 (PS1) Medium Deep Field MD07, with three prior years of multi-band PS1 light curves. The SDSS-RM 6-month baseline program aims to detect time lags between the quasar continuum and broad line region (BLR) variability on timescales of up to several months (in the observed frame) for ~10% of the sample, and to anchor the time baseline for continued monitoring in the future to detect lags on longer timescales and at higher redshift. SDSS-RM is the first major program to systematically explore the potential of RM for broad-line quasars at z>0.3, and will investigate the prospects of RM with all major broad lines covered in optical spectroscopy. SDSS-RM will provide guidance on future multi-object RM campaigns on larger scales, and is aiming to deliver more than tens of BLR lag detections for a homogeneous sample of quasars. We describe the motivation, design and implementation of this program, and outline the science impact expected from the resulting data for RM and general quasar science.

preprint2014arXiv

X-ray constraints on the local supermassive black hole occupation fraction

Distinct seed formation mechanisms are imprinted upon the fraction of dwarf galaxies currently containing a central supermassive black hole. Seeding by Pop III remnants is expected to produce a higher occupation fraction than is generated with direct gas collapse precursors. Chandra observations of nearby early-type galaxies can directly detect even low-level supermassive black hole activity, and the active fraction immediately provides a firm lower limit to the occupation fraction. Here, we use the volume-limited AMUSE surveys of ~200 optically-selected early-type galaxies to characterize simultaneously, for the first time, the occupation fraction and the scaling of nuclear X-ray luminosity with stellar mass, accounting for intrinsic scatter, measurement uncertainties, and X-ray limits. For early-type galaxies with log(M_star/M_sun)<10, we obtain a lower limit to the occupation fraction of >20% (at 95% confidence), but full occupation cannot be excluded. The preferred dependence of log(L_X) upon log(M_star) has a slope of about 0.7-0.8, consistent with the "downsizing" trend previously identified from the AMUSE dataset, and a uniform Eddington efficiency is disfavored at ~2 sigma. We provide guidelines for the future precision with which these parameters may be refined with larger or more sensitive samples.

preprint2013arXiv

Candidate Type II Quasars at 2 < z < 4.3 in the Sloan Digital Sky Survey III

At low redshifts, dust-obscured quasars often have strong yet narrow permitted lines in the rest-frame optical and ultraviolet, excited by the central active nucleus, earning the designation Type II quasars. We present a sample of 145 candidate Type II quasars at redshifts between 2 and 4.3, encompassing the epoch at which quasar activity peaked in the universe. These objects, selected from the quasar sample of the Baryon Oscillation Spectroscopic Survey of the Sloan Digital Sky Survey III, are characterized by weak continuum in the rest-frame ultraviolet (typical continuum magnitude of i \approx 22) and strong lines of CIV and Ly α, with Full Width at Half Maximum less than 2000 kms-1. The continuum magnitudes correspond to an absolute magnitude of -23 or brighter at redshift 3, too bright to be due exclusively to the host galaxies of these objects. Roughly one third of the objects are detected in the shorter-wavelength bands of the WISE survey; the spectral energy distributions (SEDs) of these objects appear to be intermediate between classic Type I and Type II quasars seen at lower redshift. Five objects are detected at rest frame 6μm by Spitzer, implying bolometric luminosities of several times 10^46 erg s-1. We have obtained polarization measurements for two objects; they are roughly 3% polarized. We suggest that these objects are luminous quasars, with modest dust extinction (A_V ~ 0.5 mag), whose ultraviolet continuum also includes a substantial scattering contribution. Alternatively, the line of sight to the central engines of these objects may be partially obscured by optically thick material.

preprint2013arXiv

Dual Supermassive Black Hole Candidates in the AGN and Galaxy Evolution Survey

Dual supermassive black holes (SMBHs) with kiloparsec scale separations in merger-remnant galaxies are informative tracers of galaxy evolution, but the avenue for identifying them in large numbers for such studies is not yet clear. One promising approach is to target spectroscopic signatures of systems where both SMBHs are fueled as dual active galactic nuclei (AGNs), or where one SMBH is fueled as an offset AGN. Dual AGNs may produce double-peaked narrow AGN emission lines, while offset AGNs may produce single-peaked narrow AGN emission lines with line-of-sight velocity offsets relative to the host galaxy. We search for such dual and offset systems among 173 Type 2 AGNs at z<0.37 in the AGN and Galaxy Evolution Survey (AGES), and we find two double-peaked AGNs and five offset AGN candidates. When we compare these results to a similar search of the DEEP2 Galaxy Redshift Survey and match the two samples in color, absolute magnitude, and minimum velocity offset, we find that the fraction of AGNs that are dual SMBH candidates increases from z=0.25 to z=0.7 by a factor of ~6 (from 2/70 to 16/91, or 2.9% to 18%). This may be associated with the rise in the galaxy merger fraction over the same cosmic time. As further evidence for a link with galaxy mergers, the AGES offset and dual AGN candidates are tentatively ~3 times more likely than the overall AGN population to reside in a host galaxy that has a companion galaxy (from 16/173 to 2/7, or 9% to 29%). Follow-up observations of the seven offset and dual AGN candidates in AGES will definitively distinguish velocity offsets produced by dual SMBHs from those produced by narrow-line region kinematics, and will help sharpen our observational approach to detecting dual SMBHs.

preprint2013arXiv

Dwarf Galaxies with Optical Signatures of Active Massive Black Holes

We present a sample of 151 dwarf galaxies (10^8.5 < M_stellar < 10^9.5 Msun) that exhibit optical spectroscopic signatures of accreting massive black holes (BHs), increasing the number of known active galaxies in this stellar mass range by more than an order of magnitude. Utilizing data from the Sloan Digital Sky Survey Data Release 8 and stellar masses from the NASA-Sloan Atlas, we have systematically searched for active BHs in ~25,000 emission-line galaxies with stellar masses comparable to the Magellanic Clouds and redshifts z<0.055. Using the narrow-line [OIII]/H-beta versus [NII]/H-alpha diagnostic diagram, we find photoionization signatures of BH accretion in 136 galaxies, a small fraction of which also exhibit broad H-alpha emission. For these broad-line AGN candidates, we estimate BH masses using standard virial techniques and find a range of 10^5 < M_BH < 10^6 Msun and a median of M_BH ~ 2 x 10^5 Msun. We also detect broad H-alpha in 15 galaxies that have narrow-line ratios consistent with star-forming galaxies. Follow-up observations are required to determine if these are true type 1 AGN or if the broad H-alpha is from stellar processes. The median absolute magnitude of the host galaxies in our active sample is Mg = -18.1 mag, which is ~1-2 magnitudes fainter than previous samples of AGN hosts with low-mass BHs. This work constrains the smallest galaxies that can form a massive BH, with implications for BH feedback in low-mass galaxies and the origin of the first supermassive BH seeds.

preprint2013arXiv

Extragalactic Science, Cosmology and Galactic Archaeology with the Subaru Prime Focus Spectrograph (PFS)

The Subaru Prime Focus Spectrograph (PFS) is a massively-multiplexed fiber-fed optical and near-infrared 3-arm spectrograph (N_fiber=2400, 380<lambda<1260nm, 1.3 degree diameter FoV), offering unique opportunities in survey astronomy. Here we summarize the science case feasible for a survey of Subaru 300 nights. We describe plans to constrain the nature of dark energy via a survey of emission line galaxies spanning a comoving volume of 9.3 (Gpc/h)^3 in the redshift range 0.8<z<2.4. In each of 6 redshift bins, the cosmological distances will be measured to 3% precision via BAO, and redshift-space distortions will be used to constrain structure growth to 6% precision. In the GA program, radial velocities and chemical abundances of stars in the Milky Way and M31 will be used to infer the past assembly histories of spiral galaxies and the structure of their dark matter halos. Data will be secured for 10^6 stars in the Galactic thick-disk, halo and tidal streams as faint as V~22, including stars with V < 20 to complement the goals of the Gaia mission. A medium-resolution mode with R = 5000 to be implemented in the red arm will allow the measurement of multiple alpha-element abundances and more precise velocities for Galactic stars, elucidating the detailed chemo-dynamical structure and evolution of each of the main stellar components of the Milky Way Galaxy and of its dwarf spheroidal galaxies. For the extragalactic program, our simulations suggest the wide avelength range will be powerful in probing the galaxy population and its clustering over a wide redshift range. We propose to conduct a color-selected survey of 1<z<2 galaxies and AGN over 16 deg^2 to J~23.4, yielding a fair sample of galaxies with stellar masses above ~10^{10}Ms at z~2. A two-tiered survey of higher redshift LBGs and LAEs will quantify the properties of early systems close to the reionization epoch.

preprint2013arXiv

Observations of Feedback from Radio-Quiet Quasars - II. Kinematics of Ionized Gas Nebulae

The prevalence and energetics of quasar feedback is a major unresolved problem in galaxy formation theory. In this paper, we present Gemini Integral Field Unit observations of ionized gas around eleven luminous, obscured, radio-quiet quasars at z~0.5 out to ~15 kpc from the quasar; specifically, we measure the kinematics and morphology of [O III]5007 emission. The round morphologies of the nebulae and the large line-of-sight velocity widths (with velocities containing 80% of the emission as high as 1000 km/s combined with relatively small velocity difference across them (from 90 to 520 km/s) point toward wide-angle quasi-spherical outflows. We use the observed velocity widths to estimate a median outflow velocity of 760 km/s, similar to or above the escape velocities from the host galaxies. The line-of-sight velocity dispersion declines slightly toward outer parts of the nebulae (by 3% per kpc on average). The majority of nebulae show blueshifted excesses in their line profiles across most of their extents, signifying gas outflows. For the median outflow velocity, we find a kinetic energy flow between 4x10^{44} and 3x10^{45} erg/s and mass outflow rate between 2000 and 20000 Msun/yr. These values are large enough for the observed quasar winds to have a significant impact on their host galaxies. The median rate of converting bolometric luminosity to kinetic energy of ionized gas clouds is ~2%. We report four new candidates for "super-bubbles" -- outflows that may have broken out of the denser regions of the host galaxy.

preprint2013arXiv

Observations of Feedback from Radio-Quiet Quasars: I. Extents and Morphologies of Ionized Gas Nebulae

Black hole feedback -- the strong interaction between the energy output of supermassive black holes and their surrounding environments -- is routinely invoked to explain the absence of overly luminous galaxies, the black hole vs. bulge correlations and the similarity of black hole accretion and star formation histories. Yet direct probes of this process in action are scarce and limited to small samples of active nuclei. We present Gemini IFU observations of the distribution of ionized gas around luminous, obscured, radio-quiet (RQ) quasars at z~0.5. We detect extended ionized gas nebulae via [O III]5007 emission in every case, with a mean diameter of 28 kpc. These nebulae are nearly perfectly round. The regular morphologies of nebulae around RQ quasars are in striking contrast with lumpy or elongated nebulae seen around radio galaxies at low and high redshifts. We present the uniformly measured size-luminosity relationship of [O III] nebulae around Seyfert 2 galaxies and type 2 quasars spanning 6 orders of magnitude in luminosity and confirm the flat slope of the correlation (R ~ L^{0.25+/-0.02}). We find a universal behavior of the [O III]/H-beta ratio in our entire RQ quasar sample: it persists at a constant value (~10) in the central regions, until reaching a "break" isophotal radius ranging from 4 to 11 kpc where it starts to decrease. We propose a model of clumpy nebulae in which clouds that produce line emission transition from being ionization-bounded at small distances from the quasar to being matter-bounded in the outer parts of the nebula, which qualitatively explains the observed line ratio and surface brightness profiles. It is striking that we see such smooth and round large-scale gas nebulosities in this sample, which are inconsistent with illuminated merger debris and which we suggest may be the signature of accretion energy from the nucleus reaching gas at large scales.

preprint2013arXiv

Refining the M_BH-V_c scaling relation with HI rotation curves of water megamaser galaxies

Black hole - galaxy scaling relations provide information about the coevolution of supermassive black holes and their host galaxies. We compare the black hole mass - circular velocity (MBH - Vc) relation with the black hole mass - bulge stellar velocity dispersion (MBH - sigma) relation, to see whether the scaling relations can passively emerge from a large number of mergers, or require a physical mechanism, such as feedback from an active nucleus. We present VLA H I observations of five galaxies, including three water megamaser galaxies, to measure the circular velocity. Using twenty-two galaxies with dynamical MBH measurements and Vc measurements extending to large radius, our best-fit MBH - Vc relation, log MBH = alpha + beta log(Vc /200 km s^-1), yields alpha = 7.43+/-0.13, beta = 3.68+1.23/-1.20, and intrinsic scatter epsilon_int = 0.51+0.11/-0.09. The intrinsic scatter may well be higher than 0.51, as we take great care to ascribe conservatively large observational errors. We find comparable scatter in the MBH - sigma relations, epsilon_int = 0.48+0.10/-0.08, while pure merging scenarios would likely result in a tighter scaling with the dark halo (as traced by Vc) than baryonic (sigma) properties. Instead, feedback from the active nucleus may act on bulge scales to tighten the MBH - sigma relation with respect to the MBH - Vc relation, as observed.

preprint2013arXiv

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.

preprint2013arXiv

Search for Supermassive Black Hole Binaries in the Sloan Digital Sky Survey Spectroscopic Sample

Supermassive black hole (SMBH) binaries are expected in a Lambda CDM cosmology given that most (if not all) massive galaxies contain a massive black hole at their center. So far, however, direct evidence for such binaries has been elusive. We use cross-correlation to search for temporal velocity shifts in the MgII broad emission lines of 0.36 < z < 2 quasars with multiple observations in the Sloan Digital Sky Survey. For ~ 10^9 Msun BHs in SMBH binaries, we are sensitive to velocity drifts for binary separations of ~ 0.1 pc with orbital periods of ~100 years. We find seven candidate sub-pc--scale binaries with velocity shifts > 3.4 sigma ~ 280 km/s, where sigma is our systematic error. Comparing the detectability of SMBH binaries with the number of candidates (N < 7), we can rule out that most 10^9 Msun BHs exist in ~ 0.03-0.2 pc scale binaries, in a scenario where binaries stall at sub-pc scales for a Hubble time. We further constrain that < one-third of quasars host SMBH binaries after considering gas-assisted sub-pc evolution of SMBH binaries, although this result is very sensitive to the assumed size of the broad line region. We estimate the detectability of SMBH binaries with ongoing or next-generation surveys (e.g., BOSS, Subaru Prime Focus Spectrograph), taking into account the evolution of the sub-parsec binary in circumbinary gas disks. These future observations will provide longer time baselines for searches similar to ours and may in turn constrain the evolutionary scenarios of SMBH binaries.

preprint2013arXiv

The Effect of Spatial Gradients in Stellar Mass-to-Light Ratio on Black Hole Mass Measurements

We have tested the effect of spatial gradients in stellar mass-to-light ratio (Y) on measurements of black hole masses (MBH) derived from stellar orbit superposition models. Such models construct a static gravitational potential for a galaxy and its central black hole, but typically assume spatially uniform Y. We have modeled three giant elliptical galaxies with gradients alpha = d(log Y)/d(log r) from -0.2 to +0.1. Color and line strength gradients suggest mildly negative alpha in these galaxies. Introducing a negative (positive) gradient in Y increases (decreases) the enclosed stellar mass near the center of the galaxy and leads to systematically smaller (larger) MBH measurements. For models with alpha = -0.2, the best-fit values of MBH are 28%, 27%, and 17% lower than the constant-Y case, in NGC 3842, NGC 6086, and NGC 7768, respectively. For alpha = +0.1, MBH are 14%, 22%, and 17% higher than the constant-Y case for the three respective galaxies. For NGC 3842 and NGC 6086, this bias is comparable to the statistical errors from individual modeling trials. At larger radii, negative (positive) gradients in Y cause the total stellar mass to decrease (increase) and the dark matter fraction within one effective radius to increase (decrease).

preprint2013arXiv

The Lick AGN Monitoring Project 2011: Fe II Reverberation from the Outer Broad-Line Region

The prominent broad Fe II emission blends in the spectra of active galactic nuclei have been shown to vary in response to continuum variations, but past attempts to measure the reverberation lag time of the optical Fe II lines have met with only limited success. Here we report the detection of Fe II reverberation in two Seyfert 1 galaxies, NGC 4593 and Mrk 1511, based on data from a program carried out at Lick Observatory in Spring 2011. Light curves for emission lines including H-beta and Fe II were measured by applying a fitting routine to decompose the spectra into several continuum and emission-line components, and we use cross-correlation techniques to determine the reverberation lags of the emission lines relative to V-band light curves. In both cases the measured lag (t_cen) of Fe II is longer than that of H-beta, although the inferred lags are somewhat sensitive to the choice of Fe II template used in the fit. For spectral decompositions done using the Fe II template of Veron-Cetty et al. (2004), we find t_cen(Fe II)/t_cen(H-beta) = 1.9+-0.6 in NGC 4593 and 1.5+-0.3 in Mrk 1511. The detection of highly correlated variations between Fe II and continuum emission demonstrates that the Fe II emission in these galaxies originates in photoionized gas, located predominantly in the outer portion of the broad-line region.

preprint2013arXiv

The Low-Luminosity End of the Radius-Luminosity Relationship for Active Galactic Nuclei

We present an updated and revised analysis of the relationship between the Hbeta broad-line region (BLR) radius and the luminosity of the active galactic nucleus (AGN). Specifically, we have carried out two-dimensional surface brightness decompositions of the host galaxies of 9 new AGNs imaged with the Hubble Space Telescope Wide Field Camera 3. The surface brightness decompositions allow us to create "AGN-free" images of the galaxies, from which we measure the starlight contribution to the optical luminosity measured through the ground-based spectroscopic aperture. We also incorporate 20 new reverberation-mapping measurements of the Hbeta time lag, which is assumed to yield the average Hbeta BLR radius. The final sample includes 41 AGNs covering four orders of magnitude in luminosity. The additions and updates incorporated here primarily affect the low-luminosity end of the R-L relationship. The best fit to the relationship using a Bayesian analysis finds a slope of alpha = 0.533 (+0.035/-0.033), consistent with previous work and with simple photoionization arguments. Only two AGNs appear to be outliers from the relationship, but both of them have monitoring light curves that raise doubt regarding the accuracy of their reported time lags. The scatter around the relationship is found to be 0.19(+/-0.02) dex, but would be decreased to 0.13 dex by the removal of these two suspect measurements. A large fraction of the remaining scatter in the relationship is likely due to the inaccurate distances to the AGN host galaxies. Our results help support the possibility that the R-L relationship could potentially be used to turn the BLRs of AGNs into standardizable candles. This would allow the cosmological expansion of the Universe to be probed by a separate population of objects, and over a larger range of redshifts.

preprint2013arXiv

Using Megamaser Disks to Probe Black Hole Accretion

We examine the alignment between H_2O megamaser disks on sub-pc scales with circumnuclear disks and bars on <500 pc scales observed with HST/WFC3. The HST imaging reveals young stars, indicating the presence of gas. The megamaser disks are not well aligned with the circumnuclear bars or disks as traced by stars in the HST images. We speculate on the implications of the observed misalignments for fueling supermassive black holes in gas-rich spiral galaxies. In contrast, we find a strong preference for the rotation axes of the megamaser disks to align with radio continuum jets observed on >50 pc scales, in those galaxies for which radio continuum detections are available. Sub-arcsecond observations of molecular gas with ALMA will enable a more complete understanding of the interplay between circumnuclear structures.

preprint2013arXiv

X-ray properties of K-selected galaxies at 0.5<z<2.0: Investigating trends with stellar mass, redshift and spectral type

We examine how the total X-ray luminosity correlates with stellar mass, stellar population, and redshift for a K-band limited sample of ~3500 galaxies at 0.5<z<2.0 from the NEWFIRM Medium Band Survey in the COSMOS field. The galaxy sample is divided into 32 different galaxy types, based on similarities between the spectral energy distributions. For each galaxy type, we further divide the sample into bins of redshift and stellar mass, and perform an X-ray stacking analysis using the Chandra COSMOS (C-COSMOS) data. We find that full band X-ray luminosity is primarily increasing with stellar mass, and at similar mass and spectral type is higher at larger redshifts. When comparing at the same stellar mass, we find that the X-ray luminosity is slightly higher for younger galaxies (i.e., weaker 4000Åbreaks), but the scatter in this relation is large. We compare the observed X-ray luminosities to those expected from low and high mass X-ray binaries (XRBs). For blue galaxies, XRBs can almost fully account for the observed emission, while for older galaxies with larger 4000Åbreaks, active galactic nuclei (AGN) or hot gas dominate the measured X-ray flux. After correcting for XRBs, the X-ray luminosity is still slightly higher in younger galaxies, although this correlation is not significant. AGN appear to be a larger component of galaxy X-ray luminosity at earlier times, as the hardness ratio increases with redshift. Together with the slight increase in X-ray luminosity this may indicate more obscured AGNs or higher accretion rates at earlier times.

preprint2012arXiv

Chandra X-ray and Hubble Space Telescope Imaging of Optically Selected kiloparsec-Scale Binary Active Galactic Nuclei I. Nature of the Nuclear Ionizing Sources

Kiloparsec-scale binary active galactic nuclei (AGNs) signal active supermassive black hole (SMBH) pairs in merging galaxies. Despite their significance, unambiguously confirmed cases remain scarce and most have been discovered serendipitously. In a previous systematic search, we optically identified four kpc-scale binary AGNs from candidates selected with double-peaked narrow emission lines at redshifts of 0.1--0.2. Here we present Chandra and Hubble Space Telescope Wide Field Camera 3 (WFC3) imaging of these four systems. We critically examine and confirm the binary-AGN scenario for two of the four targets, by combining high angular resolution X-ray imaging spectroscopy with Chandra ACIS-S, better nuclear position constraints from WFC3 F105W imaging, and direct starburst estimates from WFC3 F336W imaging; for the other two targets, the existing data are still consistent with the binary-AGN scenario, but we cannot rule out the possibility of only one AGN ionizing gas in both merging galaxies. We find tentative evidence for a systematically smaller X-ray-to-[O III] luminosity ratio and/or higher Compton-thick fraction in optically selected kpc-scale binary AGNs than in single AGNs, possibly caused by a higher nuclear gas column due to mergers and/or a viewing angle bias related to the double-peak narrow line selection. While our result lends some further support to the general approach of optically identifying kpc-scale binary AGNs, it also highlights the challenge and ambiguity of X-ray confirmation.

preprint2012arXiv

Detectors and cryostat design for the SuMIRe Prime Focus Spectrograph (PFS)

We describe the conceptual design of the camera cryostats, detectors, and detector readout electronics for the SuMIRe Prime Focus Spectrograph (PFS) being developed for the Subaru telescope. The SuMIRe PFS will consist of four identical spectrographs, each receiving 600 fibers from a 2400 fiber robotic positioner at the prime focus. Each spectrograph will have three channels covering wavelength ranges 3800 Å - 6700 Å, 6500 Å - 10000 Å, and 9700 Å - 13000 Å, with the dispersed light being imaged in each channel by a f/1.10 vacuum Schmidt camera. In the blue and red channels a pair of Hamamatsu 2K x 4K edge-buttable CCDs with 15 um pixels are used to form a 4K x 4K array. For the IR channel, the new Teledyne 4K x 4K, 15 um pixel, mercury-cadmium-telluride sensor with substrate removed for short-wavelength response and a 1.7 um cutoff will be used. Identical detector geometry and a nearly identical optical design allow for a common cryostat design with the only notable difference being the need for a cold radiation shield in the IR camera to mitigate thermal background. This paper describes the details of the cryostat design and cooling scheme, relevant thermal considerations and analysis, and discusses the detectors and detector readout electronics.

preprint2012arXiv

Physical Properties of the Narrow-Line Region of Low-Mass Active Galaxies

We present spectroscopic observations of 27 active galactic nuclei (AGN) with some of the lowest black hole (BH) masses known. We use the high spectral resolution and small aperture of our Keck data, taken with the Echellette Spectrograph and Imager, to isolate the narrow-line regions (NLRs) of these low-mass BHs. We investigate their emission-line properties and compare them with those of AGN with higher-mass black holes. While we are unable to determine absolute metallicities, some of our objects plausibly represent examples of the low-metallicity AGN described by Groves et al. (2006), based on their [N II]/H_alpha ratios and their consistency with the Kewley & Ellison (2008) mass-metallicity relation. We find tentative evidence for steeper far-UV spectral slopes in lower-mass systems. Overall, NLR emission lines in these low-mass AGN exhibit trends similar to those seen in AGN with higher-mass BHs, such as increasing blueshifts and broadening with increasing ionization potential. Additionally, we see evidence of an intermediate line region whose intensity correlates with L/L_Edd, as seen in higher-mass AGN. We highlight the interesting trend that, at least in these low-mass AGN, the [O III] equivalent width (EW) is highest in symmetric NLR lines with no blue wing. This trend of increasing [O III] EW with line symmetry could be explained by a high covering factor of lower ionization gas in the NLR. In general, low-mass AGN preserve many well-known trends in the structure of the NLR, while exhibiting steeper ionizing continuum slopes and somewhat lower gas-phase metallicities.

preprint2012arXiv

Simultaneous Ultraviolet and Optical Emission-line Profiles of Quasars: Implications for Black Hole Mass Determination

The X-shooter instrument on the VLT was used to obtain spectra of seven moderate-redshift quasars simultaneously covering the spectral range 3000 Ang to 2.5 microns. At z ~ 1.5, most of the prominent broad emission lines in the ultraviolet to optical region are captured in their rest frame. We use this unique dataset, which mitigates complications from source variability, to intercompare the line profiles of C IV 1549, C III] 1909, Mg II 2800, and Halpha and evaluate their implications for black hole mass estimation. We confirm that Mg II and the Balmer lines share similar kinematics and that they deliver mutually consistent black hole mass estimates with minimal internal scatter (< 0.1 dex) using the latest virial mass estimators. Although no virial mass formalism has yet been calibrated for C III], this line does not appear promising for such an application because of the large spread of its velocity width compared to lines of both higher and lower ionization; part of the discrepancy may be due to the difficulty of deblending C III] from its neighboring lines. The situation for C IV is complex and, because of the limited statistics of our small sample, inconclusive. On the one hand, slightly more than half of our sample (4/7) have C IV line widths that correlate reasonably well with Halpha line widths, and their respective black hole mass estimates agree to within ~0.15 dex. The rest, on the other hand, exhibit exceptionally broad C IV profiles that overestimate virial masses by factors of 2-5 compared to Halpha. As C IV is widely used to study black hole demographics at high redshifts, we urgently need to revisit our analysis with a larger sample.

preprint2012arXiv

Stellar Populations of UV-Selected Active Galactic Nuclei Host Galaxies at z ~ 2 - 3

We use stellar population synthesis modeling to analyze the host galaxy properties of a sample of 33 UV-selected, narrow-lined active galactic nuclei (AGNs) at z ~ 2 - 3. In order to quantify the contribution of AGN emission to host galaxy broadband spectral energy distributions (SEDs), we use the subsample of 11 AGNs with photometric coverage spanning from rest-frame UV through near-IR wavelengths. Modeling the SEDs of these objects with a linear combination of stellar population and AGN templates, we infer the effect of the AGN on derived stellar population parameters. We also estimate the typical bias in derived stellar populations for AGNs lacking rest-frame near-IR wavelength coverage, and develop a method for inferring the true host galaxy properties. We compare AGN host galaxy properties to those of a sample of UV-selected, star-forming non-AGNs in the same redshift range, including a subsample carefully matched in stellar mass. Although the AGNs have higher masses and SFRs than the full non-active sample, their stellar population properties are consistent with those of the mass-selected sample, suggesting that the presence of an AGN is not connected with the cessation of star-formation activity in star-forming galaxies at z ~ 2 - 3. We suggest that a correlation between M_BH and galaxy stellar mass is already in place at this epoch. Assuming a roughly constant Eddington ratio for AGNs at all stellar masses, we are unable to detect the AGNs in low-mass galaxies because they are simply too faint.

preprint2012arXiv

The Lick AGN Monitoring Project 2011: Dynamical Modeling of the Broad Line Region in Mrk 50

We present dynamical modeling of the broad line region (BLR) in the Seyfert 1 galaxy Mrk 50 using reverberation mapping data taken as part of the Lick AGN Monitoring Project (LAMP) 2011. We model the reverberation mapping data directly, constraining the geometry and kinematics of the BLR, as well as deriving a black hole mass estimate that does not depend on a normalizing factor or virial coefficient. We find that the geometry of the BLR in Mrk 50 is a nearly face-on thick disk, with a mean radius of 9.6(+1.2,-0.9) light days, a width of the BLR of 6.9(+1.2,-1.1) light days, and a disk opening angle of 25\pm10 degrees above the plane. We also constrain the inclination angle to be 9(+7,-5) degrees, close to face-on. Finally, the black hole mass of Mrk 50 is inferred to be log10(M(BH)/Msun) = 7.57(+0.44,-0.27). By comparison to the virial black hole mass estimate from traditional reverberation mapping analysis, we find the normalizing constant (virial coefficient) to be log10(f) = 0.78(+0.44,-0.27), consistent with the commonly adopted mean value of 0.74 based on aligning the M(BH)-σ* relation for AGN and quiescent galaxies. While our dynamical model includes the possibility of a net inflow or outflow in the BLR, we cannot distinguish between these two scenarios.

preprint2012arXiv

The Stellar Halos of Massive Elliptical Galaxies

We use the Mitchell Spectrograph (formerly VIRUS-P) on the McDonald Observatory 2.7m Harlan J. Smith Telescope to search for the chemical signatures of massive elliptical galaxy assembly. The Mitchell Spectrograph is an integral-field spectrograph with a uniquely wide field of view (107x107 sq arcsec), allowing us to achieve remarkably high signal-to-noise ratios of ~20-70 per pixel in radial bins of 2-2.5 times the effective radii of the eight galaxies in our sample. Focusing on a sample of massive elliptical galaxies with stellar velocity dispersions sigma* > 150 km/s, we study the radial dependence in the equivalent widths (EWs) of key metal absorption lines. By twice the effective radius, the Mgb EWs have dropped by ~50%, and only a weak correlation between sigma* and Mgb EW remains. The Mgb EWs at large radii are comparable to those seen in the centers of elliptical galaxies that are approximately an order of magnitude less massive. We find that the well-known metallicity gradients often observed within an effective radius continue smoothly to 2.5R_e, while the abundance ratio gradients remain flat. Much like the halo of the Milky Way, the stellar halos of our galaxies have low metallicities and high alpha-abundance ratios, as expected for very old stars formed in small stellar systems. Our observations support a picture in which the outer parts of massive elliptical galaxies are built by the accretion of much smaller systems whose star formation history was truncated at early times.

preprint2012arXiv

X-ray Properties of Intermediate-mass Black Holes in Active Galaxies. III. Spectral Energy Distribution and Possible Evidence for Intrinsically X-ray-weak AGNs

We present a systematic X-ray study, the third in a series, of 49 active galactic nuclei with intermediate-mass black holes (IMBH; ~10^5-10^6 M_sun) using Chandra observations. We detect 42 out of 49 targets with a 0.5-2 keV X-ray luminosity 10^41-10^43 erg/s. We perform spectral fitting for the 10 objects with enough counts (>200), and they are all well fit by a simple power-law model modified by Galactic absorption, with no sign of significant intrinsic absorption. While we cannot fit the X-ray spectral slope directly for the rest of the sample, we estimate it from the hardness ratio and find a range of photon indices consistent with those seen in more luminous and massive objects. The X-ray-to-optical spectral slope (alphaox) of our IMBH sample is systematically flatter than in active galaxies with more massive black holes, consistent with the well-known correlation between alphaox and UV luminosity. Thanks to the wide dynamic range of our sample, we find evidence that alphaox increases with decreasing M_BH as expected from accretion disk models, where the UV emission systematically decreases as M_BH decreases and the disk temperature increases. We also find a long tail toward low alphaox values. While some of these sources may be obscured, given the high L_bol/L_Eddington values in the sample, we argue that some may be intrinsically X-ray-weak, perhaps owing to a rare state that radiates very little coronal emission.

preprint2011arXiv

A spectacular outflow in an obscured quasar

SDSS J1356+1026 is a pair of interacting galaxies at redshift z=0.123 that hosts a luminous obscured quasar in its northern nucleus. Here we present two long-slit Magellan LDSS-3 spectra that reveal a pair of symmetric ~10 kpc-size outflows emerging from this nucleus, with observed expansion velocities of ~250 km/s in projection. We present a kinematic model of these outflows and argue that the deprojected physical velocities of expansion are likely ~1000 km/s and that the kinetic energy of the expanding shells is likely 10^44-10^45 erg/s, with an absolute minimum of >10^42 erg/s. Although a radio counterpart is detected at 1.4GHz, it is faint enough that the quasar is considered to be radio-quiet by all standard criteria, and there is no evidence of extended emission due to radio lobes, whether aged or continuously powered by an ongoing jet. We argue that the likely level of star formation is probably insufficient to power the observed energetic outflow and that SDSS J1356+1026 makes a strong case for radio-quiet quasar feedback. In further support of this hypothesis, polarimetric observations show that the direction of quasar illumination is coincident with the direction of the outflow.

preprint2011arXiv

Black Hole Mass and Bulge Luminosity for Low-mass Black Holes

We study the scaling between bulge magnitude and central black hole (BH) mass in galaxies with virial BH masses < 10^6 solar mass. Based on careful image decomposition of a snapshot Hubble Space Telescope I-band survey, we found that these BHs are found predominantly in galaxies with pseudobulges. Here we show that the \mbulge\ relation for the pseudobulges at low mass is significantly different from classical bulges with BH masses >10^7 solar mass. Specfically, bulges span a much wider range of bulge luminosity, and on average the luminosity is larger, at fixed black hole mass. The trend holds both for the active galaxies from Bentz et al. and the inactive sample of Gultekin et al. and cannot be explained by differences in stellar populations, as it persists when we use dynamical bulge masses. Put another way, the ratio between bulge and BH mass is much larger than $\sim 1000$ for our sample. This is consistent with recent suggestions that black hole mass does not scale with the pseudobulge luminosity. The low-mass scaling relations appear to flatten, consistent with predictions from Volonteri & Natarajan for massive seed BHs.

preprint2011arXiv

Exploring the Low-Mass End of The M-Sigma Relation with Active Galaxies

We present new measurements of stellar velocity dispersions, using spectra obtained with the Keck Echellette Spectrograph and Imager (ESI) and the Magellan Echellette (MagE), for 76 Seyfert 1 galaxies from the recent catalogue of Greene & Ho. These objects were selected from the Sloan Digital Sky Survey (SDSS) to have estimated black hole (BH) masses below 2\times10^6 M\odot. Combining our results with previous ESI observations of similar objects, we obtain an expanded sample of 93 galaxies and examine the relation between BH mass and velocity dispersion (the M-Sigma relation) for active galaxies with low BH masses. The low-mass active galaxies tend to follow the extrapolation of the M-Sigma relation of inactive galaxies. Including results for active galaxies of higher BH mass from the literature, we find a zero pointα= 7.68\pm0.08 and slope of β= 3.32\pm0.22 for the M-Sigma relation [in the form log Mbh=α+βlog(σ*/200 km s-1)], with intrinsic scatter of 0.46\pm0.03 dex. This result is consistent, within the uncertainties, with the slope of the M-Sigma relation for reverberation-mapped active galaxies with BH masses from 10^6 to 10^9 M\odot. For the subset of our sample having morphological information from Hubble Space Telescope images, we examine the slope of the M-Sigma relation separately for subsamples of barred and unbarred host galaxies, and find no significant evidence for a difference in slope. We do find a mild offset between low-inclination and high-inclination disk galaxies, such that more highly inclined galaxies tend to have larger σ* at a given value of BH mass, presumably due to the contribution of disk rotation within the spectroscopic aperture. We also find that the velocity dispersion of the ionized gas trace the stellar velocity dispersion well for this large sample of low-mass Seyfert 1 galaxies.

preprint2011arXiv

Feedback in Luminous Obscured Quasars

We use spatially resolved long-slit spectroscopy from Magellan to investigate the extent, kinematics, and ionization structure in the narrow-line regions of 15 luminous, obscured quasars with z<0.5. Increasing the dynamic range in luminosity by an order of magnitude, as well as improving the depth of existing observations by a similar factor, we revisit relations between narrow-line region size and the luminosity and linewidth of the narrow emission lines. We find a slope of 0.22 +/- 0.04 for the power-law relationship between size and luminosity, suggesting that the nebulae are limited by availability of gas to ionize at these luminosities. In fact, we find that the active galactic nucleus is effectively ionizing the interstellar medium over the full extent of the host galaxy. Broad (~300-1000 km/s) linewidths across the galaxies reveal that the gas is kinematically disturbed. Furthermore, the rotation curves and velocity dispersions of the ionized gas remain constant out to large distances, in striking contrast to normal and starburst galaxies. We argue that the gas in the entire host galaxy is significantly disturbed by the central active galactic nucleus. While only ~10^7-10^8 M_sun worth of gas are directly observed to be leaving the host galaxies at or above their escape velocities, these estimates are likely lower limits because of the biases in both mass and outflow velocity measurements and may in fact be in accord with expectations from recent feedback models. Additionally, we report the discovery of two dual obscured quasars, one of which is blowing a large-scale (~10 kpc) bubble of ionized gas into the intergalactic medium.

preprint2011arXiv

The Host Galaxies of Low-mass Black Holes

Using HST observations of 147 host galaxies of low-mass black holes (BHs), we systematically study the structures and scaling relations of these active galaxies. Our sample is selected to have central BHs with virial masses ~10^5-10^6 solar mass. The host galaxies have total I-band magnitudes of -23.2<M_I<-18.8 mag and bulge magnitudes of -22.9<M_I<-16.1 mag. Detailed bulge-disk-bar decompositions with GALFIT show that 93% of the galaxies have extended disks, 39% have bars and 5% have no bulges at all at the limits of our observations. Based on the Sersic index and bulge-to-total ratio, we conclude that the majority of the galaxies with disks are likely to contain pseudobulges and very few of these low-mass BHs live in classical bulges. The fundamental plane of our sample is offset from classical bulges and ellipticals in a way that is consistent with the scaling relations of pseudobulges. The sample has smaller velocity dispersion at fixed luminosity in the Faber-Jackson plane, compared with classical bulges and elliptical galaxies. The galaxies without disks are structurally more similar to spheroidals than to classical bulges according to their positions in the fundamental plane, especially the Faber-Jackson projection. Overall, we suggest that BHs with mass < 10^6 solar mass live in galaxies that have evolved secularly over the majority of their history. A classical bulge is not a prerequisite to host a black hole.

preprint2011arXiv

The Lick AGN Monitoring Project: Recalibrating Single-Epoch Virial Black Hole Mass Estimates

We investigate the calibration and uncertainties of black hole mass estimates based on the single-epoch (SE) method, using homogeneous and high-quality multi-epoch spectra obtained by the Lick Active Galactic Nucleus (AGN) Monitoring Project for 9 local Seyfert 1 galaxies with black hole masses < 10^8 M_sun. By decomposing the spectra into their AGN and stellar components, we study the variability of the single-epoch Hbeta line width (full width at half-maximum intensity, FWHM_Hbeta; or dispersion, sigma_Hbeta) and of the AGN continuum luminosity at 5100A (L_5100). From the distribution of the "virial products" (~ FWHM_Hbeta^2 L_5100^0.5 or sigma_Hbeta^2 L_5100^0.5) measured from SE spectra, we estimate the uncertainty due to the combined variability as ~ 0.05 dex (12%). This is subdominant with respect to the total uncertainty in SE mass estimates, which is dominated by uncertainties in the size-luminosity relation and virial coefficient, and is estimated to be ~ 0.46 dex (factor of ~ 3). By comparing the Hbeta line profile of the SE, mean, and root-mean-square (rms) spectra, we find that the Hbeta line is broader in the mean (and SE) spectra than in the rms spectra by ~ 0.1 dex (25%) for our sample with FWHM_Hbeta < 3000 km/s. This result is at variance with larger mass black holes where the difference is typically found to be much less than 0.1 dex. To correct for this systematic difference of the Hbeta line profile, we introduce a line-width dependent virial factor, resulting in a recalibration of SE black hole mass estimators for low-mass AGNs.

preprint2011arXiv

The Mass of the Black Hole in Arp 151 from Bayesian Modeling of Reverberation Mapping Data

Supermassive black holes are believed to be ubiquitous at the centers of galaxies. Measuring their masses is extremely challenging yet essential for understanding their role in the formation and evolution of cosmic structure. We present a direct measurement of the mass of a black hole in an active galactic nucleus (Arp 151) based on the motion of the gas responsible for the broad emission lines. By analyzing and modeling spectroscopic and photometric time series, we find that the gas is well described by a disk or torus with an average radius of 3.99 +- 1.25 light days and an opening angle of 68.9 (+21.4, -17.2) degrees, viewed at an inclination angle of 67.8 +- 7.8 degrees (that is, closer to face-on than edge-on). The black hole mass is inferred to be 10^(6.51 +- 0.28) solar masses. The method is fully general and can be used to determine the masses of black holes at arbitrary distances, enabling studies of their evolution over cosmic time.

preprint2011arXiv

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.

preprint2011arXiv

Type 2 Active Galactic Nuclei with Double-Peaked [OIII] Lines. II. Single AGNs with Complex Narrow-Line Region Kinematics are More Common than Binary AGNs

(Abridged) Approximately 1% of low redshift (z<0.3) optically-selected type 2 AGNs show a double-peaked [OIII] narrow emission line profile in their spatially-integrated spectra. Such features are usually interpreted as due either to kinematics, such as biconical outflows and/or disk rotation of the narrow line region (NLR) around single black holes, or to the relative motion of two distinct NLRs in a merging pair of AGNs. Here we report follow-up near infrared (NIR) imaging and optical slit spectroscopy of 31 double-peaked [OIII] type 2 AGNs drawn from the SDSS parent sample presented in Liu et al (2010). These data reveal a mixture of origins for the double-peaked feature. Roughly 10% of our objects are best explained by binary AGNs at (projected) kpc-scale separations, where two stellar components with spatially coincident NLRs are seen. ~ 50% of our objects have [OIII] emission offset by a few kpc, corresponding to the two velocity components seen in the SDSS spectra, but there are no corresponding double stellar components seen in the NIR imaging. For those objects with sufficiently high quality slit spectra, we see velocity and/or velocity dispersion gradients in [OIII] emission, suggestive of the kinematic signatures of a single NLR. The remaining ~40% of our objects are ambiguous, and will need higher spatial resolution observations to distinguish between the two scenarios. Our observations therefore favor the kinematics scenario with a single AGN for the majority of these double-peaked [OIII] type 2 AGNs. We emphasize the importance of combining imaging and slit spectroscopy in identifying kpc binary AGNs, i.e., in no cases does one of these alone allow an unambiguous identification. We estimate that ~ 0.5-2.5% of the z<0.3 type 2 AGNs are kpc-scale binary AGNs of comparable luminosities, with a relative orbital velocity >~150 km/s.

preprint2010arXiv

Discovery of Four kpc-Scale Binary AGNs

We report the discovery of four kpc-scale binary AGNs. These objects were originally selected from the Sloan Digital Sky Survey based on double-peaked [O III] 4959,5007 emission lines in their fiber spectra. The double peaks could result from pairing active supermassive black holes (SMBHs) in a galaxy merger, or could be due to bulk motions of narrow-line region gas around a single SMBH. Deep near-infrared (NIR) images and optical slit spectra obtained from the Magellan 6.5 m and the APO 3.5 m telescopes strongly support the binary SMBH scenario for the four objects. In each system, the NIR images reveal tidal features and double stellar components with a projected separation of several kpc, while optical slit spectra show two Seyfert 2 nuclei spatially coincident with the stellar components, with line-of-sight velocity offsets of a few hundred km/s. These objects were drawn from a sample of only 43 objects, demonstrating the efficiency of this technique to find kpc-scale binary AGNs.

preprint2010arXiv

The Lick AGN Monitoring Project: Alternate Routes to a Broad-line Region Radius

It is now possible to estimate black hole masses across cosmic time, using broad emission lines in active galaxies. This technique informs our views of how galaxies and their central black holes coevolve. Unfortunately, there are many outstanding uncertainties associated with these "virial" mass estimates. One of these comes from using the accretion luminosity to infer a size for the broad-line region. Incorporating the new sample of low-luminosity active galaxies from our recent monitoring campaign at Lick Observatory, we recalibrate the radius-luminosity relation with tracers of the accretion luminosity other than the optical continuum. We find that the radius of the broad-line region scales as the square root of the X-ray and Hbeta luminosities, in agreement with recent optical studies. On the other hand, the scaling appears to be marginally steeper with narrow-line luminosities. This is consistent with a previously observed decrease in the ratio of narrow-line to X-ray luminosity with increasing total luminosity. The radius of the broad-line region correlates most tightly with Hbeta luminosity, while the X-ray and narrow-line relations both have comparable scatter of a factor of two. These correlations provide useful alternative virial BH masses in objects with no detectable optical/UV continuum emission, such as high-redshift galaxies with broad emission lines, radio-loud objects, or local active galaxies with galaxy-dominated continua.

preprint2010arXiv

The Lick AGN Monitoring Project: Reverberation Mapping of Optical Hydrogen and Helium Recombination Lines

We have recently completed a 64-night spectroscopic monitoring campaign at the Lick Observatory 3-m Shane telescope with the aim of measuring the masses of the black holes in 12 nearby (z < 0.05) Seyfert 1 galaxies with expected masses in the range ~10^6-10^7M_sun and also the well-studied nearby active galactic nucleus (AGN) NGC 5548. Nine of the objects in the sample (including NGC 5548) showed optical variability of sufficient strength during the monitoring campaign to allow for a time lag to be measured between the continuum fluctuations and the response to these fluctuations in the broad Hbeta emission, which we have previously reported. We present here the light curves for the Halpha, Hgamma, HeII 4686, and HeI 5876 emission lines and the time lags for the emission-line responses relative to changes in the continuum flux. Combining each emission-line time lag with the measured width of the line in the variable part of the spectrum, we determine a virial mass of the central supermassive black hole from several independent emission lines. We find that the masses are generally consistent within the uncertainties. The time-lag response as a function of velocity across the Balmer line profiles is examined for six of the AGNs. Finally we compare several trends seen in the dataset against the predictions from photoionization calculations as presented by Korista & Goad. We confirm several of their predictions, including an increase in responsivity and a decrease in the mean time lag as the excitation and ionization level for the species increases. Further confirmation of photoionization predictions for broad-line gas behavior will require additional monitoring programs for these AGNs while they are in different luminosity states. [abridged]

preprint2009arXiv

Host Galaxies of Luminous Type 2 Quasars at z ~ 0.5

We present deep Gemini GMOS optical spectroscopy of nine luminous quasars at redshifts z ~ 0.5, drawn from the SDSS type 2 quasar sample. Our targets were selected to have high intrinsic luminosities (M_V < -26 mag) as indicated by the [O III] 5007 A emission-line luminosity (L_[O III]). Our sample has a median black hole mass of ~ 10^8.8 M_sun inferred assuming the local M_BH-sigma_* relation and a median Eddington ratio of ~ 0.7, using stellar velocity dispersions sigma_* measured from the G band. We estimate the contamination of the stellar continuum from scattered quasar light based on the strength of broad H-beta, and provide an empirical calibration of the contamination as a function of L_[O III]; the scattered light fraction is ~ 30% of L_5100 for objects with L_[O III] = 10^9.5 L_sun. Population synthesis indicates that young post-starburst populations (< 0.1 Gyr) are prevalent in luminous type 2 quasars, in addition to a relatively old population (> 1 Gyr) which dominates the stellar mass. Broad emission complexes around He II 4686 A with luminosities up to 10^8.3 L_sun are unambiguously detected in three out of the nine targets, indicative of Wolf-Rayet populations. Population synthesis shows that ~ 5-Myr post-starburst populations contribute substantially to the luminosities (> 50% of L_5100) of all three objects with Wolf-Rayet detections. We find two objects with double cores and four with close companions. Our results may suggest that luminous type 2 quasars trace an early stage of galaxy interaction, perhaps responsible for both the quasar and the starburst activity.

preprint2009arXiv

Redshift Evolution in Black Hole-Bulge Relations: Testing CIV-based Black Hole Masses

We re-examine claims of redshift evolution in black hole-bulge scaling relations based on lensed quasars. In particular, we refine the black hole mass estimates using measurements of Balmer lines from near-infrared spectroscopy obtained with Triplespec at Apache Point Observatory. In support of previous work, we find a large scatter between Balmer and UV line widths, both MgII 2796, 2803 and CIV 1548, 1550. There is tentative evidence that CIII] 1909, despite being a blend of multiple transitions, may correlate well with MgII, although a larger sample is needed for a real calibration. Most importantly, we find no systematic changes in the estimated BH masses for the lensed sample based on Balmer lines, providing additional support to the interpretation that black holes were overly massive compared to their host galaxies at high redshift.

preprint2009arXiv

Strong [O III] Objects Among SDSS Broad-Line Active Galaxies

We present the results of a spectral principal component analysis on 9046 broad-line AGN from the Sloan Digital Sky Survey. We examine correlations between spectral regions within various eigenspectra (e.g., between Fe II strength and H$β$ width) and confirm that the same trends are apparent in spectral measurements, as validation of our technique. Because we found that our sample had a large range in the equivalent width of [O III] $λ$5007, we divided the data into three subsets based on [O III] strength. Of these, only in the sample with the weakest equivalent width of [O III] were we able to recover the known correlation between [O III] strength and full width at half maximum of H$β$ and their anticorrelation with Fe II strength. At the low luminosities considered here ($L_{5100 Å}$ of $10^{42}-10^{46}$ erg s$^{-1}$), interpretation of the principal components is considerably complicated particularly because of the wide range in [O III] equivalent width. We speculate that variations in covering factor are responsible for this wide range in [O III] strength.

preprint2009arXiv

The Lick AGN Monitoring Project: Broad-Line Region Radii and Black Hole Masses from Reverberation Mapping of Hbeta

We have recently completed a 64-night spectroscopic monitoring campaign at the Lick Observatory 3-m Shane telescope with the aim of measuring the masses of the black holes in 12 nearby (z < 0.05) Seyfert 1 galaxies with expected masses in the range ~10^6-10^7 M_sun and also the well-studied nearby active galactic nucleus (AGN) NGC 5548. Nine of the objects in the sample (including NGC 5548) showed optical variability of sufficient strength during the monitoring campaign to allow for a time lag to be measured between the continuum fluctuations and the response to these fluctuations in the broad Hbeta emission. We present here the light curves for the objects in this sample and the subsequent Hbeta time lags for the nine objects where these measurements were possible. The Hbeta lag time is directly related to the size of the broad-line region, and by combining the lag time with the measured width of the Hbeta emission line in the variable part of the spectrum, we determine the virial mass of the central supermassive black hole in these nine AGNs. The absolute calibration of the black hole masses is based on the normalization derived by Onken et al. We also examine the time lag response as a function of velocity across the Hbeta line profile for six of the AGNs. The analysis of four leads to ambiguous results with relatively flat time lags as a function of velocity. However, SBS 1116+583A exhibits a symmetric time lag response around the line center reminiscent of simple models for circularly orbiting broad-line region (BLR) clouds, and Arp 151 shows an asymmetric profile that is most easily explained by a simple gravitational infall model. Further investigation will be necessary to fully understand the constraints placed on physical models of the BLR by the velocity-resolved response in these objects.

preprint2009arXiv

The Lick AGN Monitoring Project: Photometric Light Curves and Optical Variability Characteristics

The Lick AGN Monitoring Project targeted 13 nearby Seyfert 1 galaxies with the intent of measuring the masses of their central black holes using reverberation mapping. The sample includes 12 galaxies selected to have black holes with masses roughly in the range 10^6-10^7 solar masses, as well as the well-studied AGN NGC 5548. In conjunction with a spectroscopic monitoring campaign, we obtained broad-band B and V images on most nights from 2008 February through 2008 May. The imaging observations were carried out by four telescopes: the 0.76-m Katzman Automatic Imaging Telescope (KAIT), the 2-m Multicolor Active Galactic Nuclei Monitoring (MAGNUM) telescope, the Palomar 60-in (1.5-m) telescope, and the 0.80-m Tenagra II telescope. Having well-sampled light curves over the course of a few months is useful for obtaining the broad-line reverberation lag and black hole mass, and also allows us to examine the characteristics of the continuum variability. In this paper, we discuss the observational methods and the photometric measurements, and present the AGN continuum light curves. We measure various variability characteristics of each of the light curves. We do not detect any evidence for a time lag between the B- and V-band variations, and we do not find significant color variations for the AGNs in our sample.

preprint2009arXiv

Type 2 AGNs with Double-Peaked [O III] Lines: Narrow Line Region Kinematics or Merging Supermassive Black Hole Pairs?

We present a sample of 167 type 2 AGNs with double-peaked [O III] 4959,5007 narrow emission lines, selected from the Seventh Data Release of the Sloan Digital Sky Survey. The double-peaked profiles can be well modeled by two velocity components, blueshifted and redshifted from the systemic velocity. Half of these objects have a more prominent redshifted component. In cases where the H-beta emission line is strong, it also shows two velocity components whose line-of-sight (LOS) velocity offsets are consistent with those of [O III]. The relative LOS velocity offset between the two components is typically a few hundred km/s, larger by a factor of ~ 1.5 than the line full width at half maximum of each component. The offset correlates with the host stellar velocity dispersion sigma*. The host galaxies of this sample show systematically larger sigma*, stellar masses, and concentrations, and older luminosity-weighted mean stellar ages than a regular type 2 AGN sample matched in redshift, [O III] 5007 equivalent width and luminosity; they show no significant difference in radio properties. These double-peaked features could be due to narrow-line region kinematics, or binary black holes. The statistical properties do not show strong preference for or against either scenario, and spatially resolved optical imaging, spectroscopy, radio or X-ray followup are needed to draw firm conclusions.

preprint2009arXiv

X-ray Properties of Intermediate-Mass Black Holes in Active Galaxies. II. X-ray-Bright Accretion and Possible Evidence for Slim Disks

We present X-ray properties of optically-selected intermediate-mass (~10^5--10^6 M_Sun) black holes (BHs) in active galaxies (AGNs), using data from the Chandra X-Ray Observatory. Our observations are a continuation of a pilot study by Greene & Ho (2007). Of the 8 objects observed, 5 are detected with X-ray luminosities in the range L_0.5-2 keV = 10^41--10^43 erg s^-1, consistent with the previously observed sample. Objects with enough counts to extract a spectrum are well fit by an absorbed power law. We continue to find a range of soft photon indices 1 < Γ_s < 2.7, where N(E) \propto E^-Γ_s, consistent with previous AGN studies, but generally flatter than other narrow-line Seyfert 1 active nuclei (NLS1s). The soft photon index correlates strongly with X-ray luminosity and Eddington ratio, but does not depend on BH mass. There is no justification for the inclusion of any additional components, such as a soft excess, although this may be a function of the relative inefficiency of detecting counts above 2 keV in these relatively shallow observations. As a whole, the X-ray-to-optical spectral slope α_ox is flatter than in more massive systems, even other NLS1s. Only X-ray-selected NLS1s with very high Eddington ratios share a similar α_ox. This is suggestive of a physical change in the accretion structure at low masses and at very high accretion rates, possibly due to the onset of slim disks. Although the detailed physical explanation for the X-ray loudness of these intermediate-mass BHs is not certain, it is very striking that targets selected on the basis of optical properties should be so distinctly offset in their broader spectral energy distributions.

preprint2008arXiv

Galaxies Correlating with Ultra-high Energy Cosmic Rays

The Pierre Auger Observatory reports that 20 of the 27 highest energy cosmic rays have arrival directions within 3.2 deg of a nearby galaxy in the Veron-Cetty & Veron Catalog of Quasars and Active Galactic Nuclei (12th Ed.), with ~5 of the correlations expected by chance. In this paper we examine the correlated galaxies to gain insight into the possible UHECR sources. We find that 14 of the 21 correlated VCV galaxies are AGNs and we determine their bolometric luminosities. The remaining 7 are primarily star-forming galaxies. The bolometric luminosities of the correlated AGNs are all greater than 5 x 10^{42} erg/s, which may explain the absence of UHECRs from the Virgo region in spite of the large number of VCV galaxies in Virgo, since most of the VCV galaxies in the Virgo region are low luminosity AGNs. Interestingly, the bolometric luminosities of most of the AGNs are significantly lower than required to satisfy the minimum condition for UHECR acceleration in a continuous jet. If a UHECR-AGN correlation is substantiated with further statistics, our results lend support to the recently proposed ``giant AGN flare" mechanism for UHECR acceleration.

preprint2007arXiv

The Mass Function of Active Black Holes in the Local Universe

We present the first measurement of the black hole (BH) mass function for broad-line active galaxies in the local Universe. Using the ~9000 broad-line active galaxies from the Fourth Data Release of the Sloan Digital Sky Survey, we construct a broad-line luminosity function that agrees very well with the local soft X-ray luminosity function. Using standard virial relations, we then convert observed broad-line luminosities and widths into BH masses. A mass function constructed in this way has the unique capability to probe the mass region <10^6 M_sun, which, while insignificant in terms of total BH mass density, nevertheless may place important constraints on the mass distribution of seed BHs in the early Universe. The characteristic local active BH has a mass of ~10^7 M_sun radiating at 10% of the Eddington rate. The active fraction is a strong function of BH mass; at both higher and lower masses the active mass function falls more steeply than one would infer from the distribution of bulge luminosity. The deficit of local massive radiating BHs is a well-known phenomenon, while we present the first robust measurement of a decline in the space density of active BHs at low mass.