Black Hole - Galaxy Scaling Relationships for Active Galactic Nuclei with Reverberation Masses
We have utilized high-resolution optical Hubble Space Telescope images and deep, ground-based near-infrared images to examine the host-galaxies of 37 active galactic nuclei with reverberation-based black hole masses. Using two-dimensional image decompositions, we have separated the host galaxy from the bright central AGN, allowing a re-examination of the and relationships and the and relationships using V-H color to constrain the stellar mass-to-light ratio. We find clear correlations for all of these scaling relationships, and the best-fit correlations are generally in good agreement with the sample of early-type galaxies with from dynamical modeling and the sample of megamasers. We also find good agreement with the expectations from the Illustris simulations, although the agreement with other simulations is less clear because of the different black hole mass ranges that are probed. is found to have the least scatter, and is therefore the best predictor of among the relationships examined here. Large photometric surveys that rely on automated analysis and forego bulge-to-disk decompositions will achieve more accurate predictions if they rely on rather than . Finally, we have examined and find a clear trend with black hole mass but not galaxy mass. This trend is also exhibited by galaxies with from dynamical modeling and megamaser galaxies, as well as simulated galaxies from Illustris, rising from % at M to % at M. This scaling should be taken into account when comparing galaxy samples that are not matched in .
Subject headings:galaxies: active — galaxies: photometry — galaxies: Seyfert — galaxies: supermassive black holes
The discovery that nearly every massive galaxy hosts a supermassive black hole in its nucleus is one of the lasting legacies of the Hubble Space Telescope (HST). Direct measurements of the masses of these black holes (), using luminous tracers inside the gravitational potential of the invisible central massive object, have led to the discovery of scaling relationships between the black holes and other characteristics of their host galaxies. This is true both for the sample of mostly-quiescent galaxies with measurements of from dynamical modeling of stars or gas (e.g., Kormendy & Ho 2013) and for the sample of active galaxies that have measurements of from reverberation mapping (e.g., Bentz & Katz 2015).
Direct black hole mass measurements are time and resource intensive, and they are generally only applicable to galaxies that meet a specific set of criteria. For instance, reverberation mapping is only applicable to broad-lined active galactic nuclei (AGNs), which are rare in the local universe, whereas dynamical modeling is only applicable when the black hole sphere of influence is resolved or nearly so, which is generally only possible out to Mpc. The resource-intensive nature of these measurements as well as the limitations on the applicability of each technique mean that, in practical terms, the number of direct measurements that may be accumulated over time is necessarily limited. The scaling relationships derived from these direct measurements, however, provide valuable shortcuts for estimating based on less resource-intensive measurements, such as the bulge stellar velocity dispersion (the relationship; Ferrarese & Merritt 2000; Gebhardt et al. 2000). As such, direct measurements and the scaling relationships that are based on them provide the foundation for all other determinations, thereby providing avenues to amass large samples for studying black hole and galaxy co-evolution across galaxy types and at different look-back times (e.g., Lapi et al. 2014; Heckman & Best 2014; Kelly & Merloni 2012 and references therein).
Scaling relationships between the central black hole and the host galaxy have also become important tools for critical testing of cosmological simulations of dark matter halo mergers (e.g., Oogi et al. 2016; Degraf et al. 2011; Hopkins et al. 2010), numerical investigations of candidate seed black holes (e.g., Shirakata et al. 2016; Volonteri & Natarajan 2009; Lippai et al. 2009), cosmological modeling of galaxy and black hole growth (e.g., DeGraf et al. 2015; Kim et al. 2011; Bonoli et al. 2009; Miller et al. 2006), and investigations into black hole feedback mechanisms (e.g., Steinborn et al. 2015; Kaviraj et al. 2011; Shabala et al. 2011; Ostriker et al. 2010). Accurate measurements of the host-galaxy characteristics of black holes with direct measurements are therefore necessary and valuable. Uncorrected biases or unmitigated scatter in the galaxy measurements can adversely affect the accurate and precise calibration of widely-used scaling relationships.
In this work, we focus on characterization of the host galaxies of AGNs with reverberation-based measurements. Using high-resolution HST optical images and deep, ground-based near-infrared images, we characterize the photometric properties of the galaxies through two-dimensional image decompositions. We estimate stellar masses based on photometric colors and widely-used prescriptions. These results then allow us to recalibrate several black hole-galaxy scaling relationships, and to investigate the black hole mass to stellar mass fraction across the sample.
Throughout this work, we adopt a standard CDM cosmology of km s Mpc with and .
As part of our ongoing work with the reverberation sample of AGNs, high-resolution medium-band observations were obtained with the Hubble Space Telescope (HST). We have also recently collected deep, ground-based near-infrared imaging for a number of these galaxies at the WIYN observatory. We restrict our analysis here to the sample of 37 galaxies for which we have imaging in both the optical and near infrared. Table 1 lists the sample and details of the observations, which we describe below.
2.1. Optical Imaging
HST imaging of the galaxies in our sample was acquired with the following instrument configurations: the Advanced Camera for Surveys (ACS) High Resolution Channel (HRC) through the F550M filter, the Wide Field Planetary Camera 2 (WFPC2) with the F547M filter, and the Wide Field Camera 3 (WFC3) through the F547M filter. The medium-band filters were specifically chosen to avoid strong emission lines from the AGN and to sample a flat portion of the underlying host-galaxy spectrum. The details of these observations and the post-processing are described by Bentz et al. (2009a, 2013).
We also present here new WFC3 F547M images of eight galaxies in the sample (HST GO-11661 and GO-13816, PI Bentz). Three had not been previously observed, while prior imaging of five galaxies with ACS HRC provided a field of view () that was too narrow to capture their extended morphologies. WFC3 provides a field of view that is well matched to the galaxies in our sample, and a high spatial resolution with a pixel scale of . Each galaxy was observed for a single orbit, with a 2-point dither pattern to fill in the gap between the detectors. At each point in the dither, a short and long exposure were obtained. The short exposures ensure an unsaturated measurement of the bright central AGN at each position, while the long exposures provide more depth for resolving the fainter, extended host galaxy.
The pipeline-reduced images were corrected for cosmic rays with LACosmic (van Dokkum, 2001). Taking advantage of the linear nature of CCDs, we corrected for saturation of the AGN in the long exposures by clipping out the saturated pixels in the long exposures and replacing them with the same pixels from the short exposures taken at the same dither position, scaled up by the exposure time ratio. The individual exposures were then drizzled to a common reference and combined with AstroDrizzle.
2.2. Near-Infrared Imaging
Near-infrared imaging of 29 reverberation-mapped AGN host galaxies was obtained between fall 2011 and spring 2013 with the WIYN High-Resolution Infrared Camera (WHIRC) at the WIYN 3.5-m telescope111The WIYN Observatory is a joint facility of the University of Wisconsin-Madison, Indiana University, the National Optical Astronomy Observatory and the University of Missouri.. The camera is a Raytheon Virgo HgCdTe with a pixel scale of and a field of view of . While broad-band , , and images were obtained for a subset of the sample, the majority of the observations were devoted to -band images and we report those here.
The typical observing sequence involved many short observations of each target with a generous dither pattern between observations. This allowed for the removal of strong fringing in the band, as well as bad pixels and cosmic rays.
Images were reduced in IRAF222IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Association of Universities for Research in Astronomy (AURA) under a cooperative agreement with the National Science Foundation. following standard procedures. Strong fringing is a known problem for band images taken with WHIRC. We were able to correct for this effect by first median-combining a large number of dithered observations of a target, with each image scaled by the median sky level. Then we created a fringe mask from this combined image with the IRAF task objmasks. Finally, we used the mask with the rmfringe task to correct each image. After correcting for fringing, we computed the pixel offsets between dithered images, subtracted the mean sky background, and shifted and combined all of the images. For the final image of each object, we added back the average sky background that had been subtracted in the previous step, to ensure that the image statistics would be properly handled in the fitting process. In Figure 1 we show the final -band images for three of our targets in comparison to the Two Micron All Sky Survey (2MASS;Skrutskie et al. 2006) -band images for the same galaxies. The improvement in depth and spatial resolution provided by the WHIRC images is immediately apparent, allowing for detection and characterization of faint surface brightness features, as well as better separation of distinct photometric components.
We supplemented this sample with HST Near-Infrared Camera and Multi-Object Spectrometer (NICMOS) observations of eight additional PG quasars with the NIC2 camera through the F160W filter. The details of these observations are described by Veilleux et al. (2009). Drizzled and combined images were downloaded from MAST. For each image, we added back the subtracted sky background as recorded in the header, and then multiplied each image by the exposure time to return the image units to counts.
3. Surface Brightness Fits
Two-dimensional surface brightness fits to the AGN host galaxy images were carried out using the software Galfit (Peng et al., 2002, 2010). Galfit allows the user to model surface brightness features with a variety of analytical models. We utilized the general Sérsic (1968) profile to fit the various photometric components of each galaxy. This particular function has the form
where is the pixel surface brightness at the effective radius . An exponential disk profile is simply a Sérsic profile with an index of . Bulges typically have , with the de Vaucouleurs (1948) profile being a special case with . Bars, on the other hand, typically have . In the few cases where a galaxy displayed a ring or a strong dust lane, we utilized the truncation function to truncate the inner and/or outer regions of a Sérsic profile with to represent the ring. For these profiles, there are two quoted radii for each truncation function, which are the break radius and the softening length.
Fits to many of the optical HST images have already been published by Bentz et al. (2009a, 2013). Fits to the new WFC3 images were carried out following similar procedures. The point spread function (PSF) was modeled by StarFit (Hamilton, 2014) in an attempt to better account for slight changes in the PSF width due to telescope breathing. StarFit begins with a TinyTim PSF model (Krist, 1993) and attempts to match the telescope focus by fitting the PSF to a source in the field. Most of the galaxies did not have a suitable field star in the frame to be used as a PSF model, so we used the StarFit model derived from a star in the field of NGC 3516 as the PSF model for all eight galaxies. While this provided a slight improvement over using basic TinyTim PSF models, we still found that in several cases we needed to supplement the PSF model with a narrow Sérsic profile to properly model the AGN in each galaxy nucleus. Without the addition of this component, the Sérsic profiles for the bulge would run up to an unrealistic index of , and would often reach the default maximum value of . Such profiles are extremely peaky with very broad wings, effectively mimicking an unresolved point source and the background sky. Whenever this occurred, we added a narrow (FWHM pixel) Sérsic on top of the PSF model at the location of the AGN. The addition of this profile always resulted in realistic values for all other model components in the image, allowing all the model components to remain unfixed during the fitting. The surface brightness profile of each galaxy was then fit with a bulge and a disk model, with additional model parameters (such as a bar, barlens, or ring) being added when necessary based on inspection of the image and the residuals of the model.
For the near-infrared WHIRC images, our fitting process began by constructing a point spread function (PSF) image from an isolated field star. This first step involved analyzing a small portion of each image centered on the star. The background sky was modeled as a tilted plane, and we fit multiple Gaussians to the star (typically 4-5) with unrestricted shape parameters and initial conditions of widths graduated in size. We also allowed for a single Fourier term to provide an asymmetry in the light distribution of each Gaussian, although we first arrived at a good set of model parameters before turning this option on in the final fitting step. The end result of each PSF model image is a residual pattern (image minus model) that does not retain any “bulls-eye” or other regular pattern, and is simply consistent with noise. These models, excluding the sky component, were then used as the PSF images for fitting the galaxy.
The fits of the WHIRC -band images were guided by the solutions determined from the optical HST images because of their superior sensitivity and spatial resolution and lower sky background. The host-galaxy components in the WHIRC images were constrained to have the same characteristic radii as had been found in the optical images, scaled by the difference in the pixel scales. Furthermore, the indices of the Sérsic profiles were constrained to the values determined from the optical images. In a few cases, unsatisfactory fits of the WHIRC images led us to revisit and refine our previously-published fits to an HST image, and these new fits to the optical images were then used to guide the fits to the WHIRC images. In all cases, the final adopted fits in both the optical and near-infrared bands agree, both in the number of photometric models and their shape parameters. The sky background was again fit as a tilted plane, and the AGN and multiple field stars were fit with the PSF image. Field stars that were not fit were masked out. Because of the non-photometric conditions throughout most of our WHIRC observations, we adopted -band magnitudes for as many field stars as possible in each image from 2MASS. The final zeropoint of each WHIRC image was set by minimizing the differences between the reported 2MASS magnitudes and the Galfit magnitudes of the field stars. Figure 2 shows the HST and WHIRC images for a typical galaxy in our sample, Mrk 6, as well as the surface brightness models fit to each image and the residuals of the fits.
While the Sérsic indices for the surface brightness components fit to the optical HST images were generally allowed to remain free parameters, we followed a slightly different procedure for the analysis of the eight PG quasars included in the sample of Veilleux et al. (2009). At first, we intended to match the procedure described by Veilleux et al. (2009) so that we could adopt the galaxy magnitudes they report, but in the end we found that we preferred a modified version of the procedure, and we thus re-fit all the NICMOS images ourselves.
As with the other galaxies, we began with the optical images. These were all WFPC2 or ACS HRC images, and the small field of view did not allow for a StarFit PSF model to be built from a suitably bright field star. The PSF was instead modeled with TinyTim (Krist, 1993) and we again found that we needed to add an additional narrow Sérsic component on top of the PSF to help avoid mismatch from spacecraft breathing. Each galaxy was fit with either a single Sérsic component with , or with an exponential disk and a Sérsic component with , depending on whether our previous fits and those reported by Veilleux et al. (2009) found evidence for a disk or not. Faint field stars and galaxies in the images were fit simultaneously with the AGN and its host galaxy, rather than masked out. Once a good fit was obtained in the optical image (all of which have finer pixel scales and marginally larger fields of view than the NIC2 camera), we turned to the fitting of the NICMOS images. The NICMOS PSF was modeled by TinyTim and was subsampled by a factor of 5, and the galaxy model components were adopted from the optical images, scaled to the proper size and held fixed during the fitting of the NICMOS images, as we did with the WHIRC images. The final parameters for the adopted fits are listed in Table LABEL:tab:fits.
4. Galaxy Characteristics
The Galfit models described in the previous section constrain the observed magnitudes of the individual photometric components of the host galaxies, which are generally, but not always, related to kinematic components of the galaxy. The fitting process that we adopted allows for a direct comparison of the colors of individual components (e.g., disks, bulges, bars), or the components can be combined to investigate the total magnitude of a galaxy in each passband as well as the overall galaxy color.
4.1. Final Photometry
The optical HST magnitudes represent a few different medium-band filters rather than a true broad-band Johnson filter. For each object, we used synphot and a reddened, redshifted galaxy spectrum to determine the color difference between the filter used for the optical HST observations and a broad-band filter. We adopted the elliptical galaxy template spectrum of Kinney et al. (1996) for these calculations, but the use of an Sa or Sc galaxy spectrum does not significantly change our results. The color differences are small, mag. Similarly, we determined the difference between the magnitudes from the NICMOS images and a “true” -band magnitude using synphot and the tabulated passband for the filter provided on the WHIRC filters webpage333https://www.noao.edu/kpno/manuals/whirc/filters.html. These corrections were slightly larger in magnitude and were always in the same direction, showing a slight bias between the two filters, mag.
The and magnitudes were then corrected for Galactic extinction along the line of sight based on the values determined from the Schlafly & Finkbeiner (2011) recalibration of the Schlegel et al. (1998) dust map of the Milky Way. Table LABEL:tab:mags gives the extinction-corrected and equivalent apparent magnitudes for the integrated galaxies and for their bulges.
Based on our previous experience with Galfit as well as comparison of our fitting results with those of Veilleux et al. (2009) for several of the PG objects, we assume a typical uncertainty of mag for the integrated magnitudes of the galaxies. We also assume a typical uncertainty of mag for the integrated magnitudes of the bulge components on their own. This is not to say that all the uncertainty is in the bulges, but rather that Galfit does a good job of recovering the total galaxy flux, even if there is some abiguity in how the light is divided between the various photometric components.
The galaxies in our sample span a range of distances covering . We used synphot to determine -corrections for each galaxy in and so that we could compare -equivalent photometry and galaxy colors. In , the relatively flat spectral energy distribution (SED) of a galaxy gives rise to small corrections ranging from mag with a median of mag. In , however, there is significantly more structure to a galaxy SED, so the galaxies with the largest redshifts in our sample () have mag. The median is significantly smaller, however, at mag.
After applying the -corrections, we derive the colors, both for the integrated light of the galaxy, which we report in Table LABEL:tab:mass, as well as for each individual component (e.g., bulge, bar, disk). The range of integrated galaxy colors is , with a median color of . Most of the sample is comprised of disk galaxies, and as expected, the colors of the galaxy bulges are generally larger (redder) than the integrated colors of the entire galaxies by a median value of 0.75 mag.
To convert the observed magnitudes to luminosities, we estimated based on the apparent redshifts of the galaxies. We conservatively adopt an uncertainty of 500 km s in peculiar velocities for each distance estimate, based on the distribution of peculiar velocities derived by Tully et al. (2008). This works out to a 17% uncertainty at , and decreases as increases. We caution that this uncertainty may still be significantly underestimated in the case of individual galaxies.
Four of the nearest galaxies in the sample have distances in the literature that were derived through other techniques. We summarize these distance measurements and their potential uncertainties in Bentz et al. (2013), but in brief, they were generally retrieved from the Extragalactic Distance Database (Tully et al., 2009) and derived from an average of the distance moduli for galaxies within the same group. The exception is NGC 3227, where the distance measurement comes from an analysis of the surface brightness fluctuations of NGC 3226 (Tonry et al., 2001), with which it is interacting.
Adopted distances and their uncertainties are listed in Table LABEL:tab:mags.
4.4. Stellar Masses
We estimated the stellar mass-to-light ratio () of each galaxy using the color and the relationships tabulated by Bell & de Jong (2001) in their Table 1. Following their work, we assume solar absolute magnitudes of (Cox, 2000) and (Worthey, 1994) and we derive the expected stellar mass in the and passbands, which are identical and are listed as “ (Bd01)” in Table LABEL:tab:mass. The uncertainties on the stellar masses are based on the propagated uncertainties in the photometry and the distances. The steeper dependency of on bluer colors leads to uncertainties in the stellar masses that are roughly twice as large when based on the luminosity in as for . To be conservative, we adopt the larger uncertainties based on .
We also estimated the stellar mass-to-light ratio of each galaxy using the relationships tabulated by Into & Portinari (2013). Into & Portinari (2013) used updated population synthesis models and applied prescriptions to more accurately account for thermally pulsing asymptotic giant branch stars, which can strongly affect near-infrared photometry of galaxies. We apply their “dusty” models (their Table 6) because we have not attempted to correct for dust internal to the galaxies, and the predominantly late types of the galaxies in our sample mean that they cannot be considered “dust free”. Adopting the solar absolute magnitudes of and derived by Into & Portinari (2013) for consistency, we determined the expected stellar masses in the and passbands. Unlike the masses estimated from the relationships of Bell & de Jong (2001), the two passbands do not predict exactly the same stellar masses, but the differences are only at the 1% level. We adopt the stellar masses based on the luminosity in , again because their larger propagated uncertainties make them a more conservative choice, and we list these in Table LABEL:tab:mass as “ (IP13)”.
The stellar masses predicted by the relationships of Into & Portinari (2013) are typically a factor of 2.4 times smaller than those derived from the relationships tabulated by Bell & de Jong (2001), although for the smallest colors in the sample, they can disagree by a factor of . Some of this difference can be attributed to the choice of a “diet” Salpeter initial-mass-function (IMF) used by Bell & de Jong (2001) and a Kroupa IMF adopted by Into & Portinari (2013). If we adjust the prescription of Bell & de Jong (2001) by dex to better match a Kroupa IMF (Bell et al., 2003), then the agreement is better, although the stellar masses predicted by the Into & Portinari (2013) relations are still typically 1.7 times smaller than those predicted by Bell & de Jong (2001). This difference agrees with the factor-of-two lighter masses found by Into & Portinari (2013) when comparing their new stellar population models with previous models, which they attribute to updates like their treatment of thermally pulsing asymptotic giant branch stars. Kormendy & Ho (2013) also find that the ratios predicted by Into & Portinari (2013) are, on average, a factor of 1.34 smaller than those predicted by the dynamics of the galaxy, although it is unclear how much of the discrepancy may be attributed to dark matter. Given the uncertainties in the methods, we report the results using both the Bell & de Jong (2001) and the Into & Portinari (2013) prescriptions throughout this work.
4.5. Black Hole Masses
Black hole masses for all galaxies were drawn from the compilation of reverberation-based masses in the AGN Black Hole Mass Database (Bentz & Katz, 2015). The basic technique of reverberation mapping (Blandford & McKee, 1982; Peterson, 1993) involves time-resolved spectrophotometry collected over a long time baseline and with dense time sampling (for nearby Seyferts, this typically amounts to daily sampling over a baseline of a few months). Variations in the continuum flux are “echoed” in the broad emission lines, and the time delay between the two is based on the light-travel time between the two regions where the signals arise, namely the accretion disk and the broad line region.
The black hole masses are determined as
where is the measured time delay for a broad emission line, is the velocity width of that same emission line, and is the gravitational constant. The factor is an order-unity scaling factor that is necessary to account for the generally unknown geometry and detailed kinematics of the broad line region in the AGNs. The value of ranges from in the literature, with most current studies finding . We adopted the scaling factor of determined by Grier et al. (2013).
With the measurements of luminosities and masses derived in the previous sections, we examine several black hole scaling relationships here. Linear regressions were carried out with a Bayesian approach using the linmix_err algorithm (Kelly, 2007), which includes measurement errors in both coordinates and a component of intrinsic, random scatter. The values and uncertainties that we report for the slope, intercept, and scatter of each relationship are the median values and widths of a large number of draws from the posterior probability distribution for each term.
5.1. Black Hole Mass – Bulge Luminosity Relationship
The relationship between black hole mass and bulge luminosity, was one of the first black hole scaling relationships to be discovered (Kormendy & Richstone, 1995). However, it was soon eclipsed by the relationship (Ferrarese & Merritt, 2000; Gebhardt et al., 2000), which was initially reported to have a smaller intrinsic scatter and was therefore viewed as being the more fundamental scaling relationship. However, improvements in the black hole mass measurements, in particular, have led to much tighter relationships in recent years where the reported scatter is similar to that of the relationship (Marconi & Hunt, 2003; Gültekin et al., 2009). These studies have tended to focus on bulge-dominated galaxies while neglecting the late-type galaxies common among local Seyfert hosts.
A notable exception, however, is Wandel (2002), who drew photometry from the literature to investigate the relationship for AGN host galaxies with black hole masses from reverberation mapping. A homogeneous reanalysis of the AGN black hole masses by Peterson et al. (2004) combined with consistent bulge photometry derived from high quality HST imaging and galaxy photometric decompositions allowed Bentz et al. (2009b) to update the results of Wandel (2002), finding that for disk-dominated active galaxies is similar in form and scatter to that of bulge-dominated galaxies with predominantly quiescent black holes and masses derived from dynamical modeling.
Here, we are able to improve upon the results of Bentz et al. (2009b) by extending the sample to lower black hole masses, increasing the number of galaxies included in the fit by 40%, and by examining the relationship in both the optical and the near-infrared. This last point is an important addition because it allows for the effects of dust and recent star formation on the photometry to be mitigated.
For each galaxy, we identified the photometric component most consistent with the expected properties of a bulge. In particular, we looked for a round () photometric component with Sérsic index and . In one instance (Mrk 509), there was no such model component and so we do not include it here in the analysis of galaxy bulges. Mrk 509 is thus consistent with either a bulgeless disk galaxy or a disk galaxy with a compact bulge that we could not separate from the central AGN. Some of the PG quasars, on the other hand, were modeled by a single spheroidal component which we include as a “bulge” here. We do not attempt to discriminate between pseudobulges and classical bulges because we have limited kinematic information regarding the bulges of these galaxies. Numerous studies have shown that pseudobulge identification can be extremely uncertain when it is based solely on photometric information (e.g., Läsker et al. 2014a; Kormendy & Kennicutt 2004). In the band, we find the best-fit relationship between the black hole mass and bulge luminosity to be:
with a typical scatter of dex. This is similar to the slope found by Bentz et al. (2009b) using a smaller number of galaxies in the reverberation sample and covering a smaller range of black hole masses. The scatter is much decreased, however, from dex to dex.
In the band, we find a best-fit relationship of:
with a typical scatter of dex. Surprisingly, the scatter in the near-infrared relationship is statistically equivalent to that of the optical relationship, suggesting that dust and/or recent star formation are not strong contributors to the intrinsic scatter in the relationship. As previously mentioned, however, there is still room for improvement in the distances, so it is likely that the scatter in both the optical and near-infrared relationships could be further decreased in the future through efforts to determine distances that do not rely on the galaxy redshift.
We display these relationships in Figure 3. The solid line shows the best fit, while the gray shaded regions show the uncertainties in the fit. We denote broad-line Seyferts 1s (BLS1s) with filled circles and narrow-line Seyferts 1s (NLS1s) with open circles. We follow the original definition of Osterbrock & Pogge (1985) and select NLS1s in cases where the broad H emission line has FWHM km s. While the NLS1s tend to be associated with lower-mass black holes in lower-luminosity bulges, they exhibit the same scatter and general scaling relationship as the BLS1s. Some studies of NLS1s with black hole estimates have shown them to be significantly undermassive relative to BLS1s (e.g., Mathur et al. 2012), but we see no strong tendency for NLS1s to be undermassive relative to the other reverberation-mapped AGNs included here.
Kormendy & Ho (2013) report a near-infrared relationship in the 2MASS band for quiescent galaxies that are ellipticals or contain classical bulges, and for which black hole masses have been determined through dynamical modeling. They find a slightly steeper slope of and a scatter of dex, both of which are consistent within the errors with our finding for the active galaxy sample in . The slightly higher intercept for their sample compared to ours is increased by the color difference between the and bands, given that galaxies are typically somewhat brighter in than .
While Kormendy & Ho (2013) do not report a fit to the relationship in , they do tabulate bulge absolute magnitudes in . We fit the -band relationship matching their accepted sample and adopted uncertainties and find a slope that agrees with their value reported for the band, which is steeper than the slope that we find in for the active galaxy sample. The intercept is also somewhat higher, although the fit to their sample agrees with our findings for the active galaxy sample at the low-mass end. The fits to the Kormendy & Ho (2013) sample are shown as black dashed lines in Figure 3. It is important to note that the active and quiescent samples primarily probe different regions of parameter space in this plot: the active galaxy sample is heavily dominated by galaxies with , while the vast majority of galaxies in the quiescent sample have .
Läsker et al. (2016) report deep -band imaging and surface brightness decompositions for a sample of 9 megamaser galaxies with accurate black hole masses. We find that the megamasers are contained wholly within the scatter of the active galaxy sample presented here in the band. With the good agreement between the active galaxies, the megamasers, and the quiescent galaxy sample, we therefore refit the relationship in with all three samples combined. Based on the typical galaxy properties in 2MASS reported by Jarrett (2000), we adopt mag for the quiescent sample, which should account for any average color offset between the two filters (although we note that the scatter in values is typically mag, even for galaxies with a specific morphological type). The best fit is:
with a typical scatter of dex. While Läsker et al. (2016) do not report -band measurements for the megamaser sample, we can investigate the relationship in for the active and quiescent samples combined. When we do, we find a best fit of:
with a typical scatter of dex. These fits are displayed as the red long-dashed lines in Figure 3. In both the and bands, the best fit for the combined sample has an almost identical scatter to that found for the active sample alone, even though the combination of the samples more than doubles the number of points being fit and extends the range of by an order of magnitude. This may indicate that the galaxies in all three samples are drawn from the same parent population.
5.2. Black Hole Mass – Galaxy Luminosity Relationship
We also examined the relationship between black hole mass and total luminosity of the host galaxy. We find a clear correlation between these two measurements, in both the optical and the near-infrared. The best-fit relationships are found to be:
in the band, with a typical scatter of dex, and
in the band, with a typical scatter of dex.
While the scatter is somewhat higher than that of the relationship, the fact that there is still a relatively tight relationship found when the total galaxy luminosity is used (see Figure 4) suggests that bulge/disk decompositions can be avoided when estimating black hole masses from broad-band photometry of disk galaxies, but with a loss of some accuracy. This may be of particular interest for large photometric surveys that are operational or coming online soon (e.g., LSST), where automated measurements will be key to making sense of the large datasets that will be produced.
Our best-fit relationships for active galaxies may again be compared to the Kormendy & Ho (2013) sample of quiescent galaxies. The best-fit relationships based on their tabulated measurements in and have similar slopes and scatter to our findings, but their intercepts are significantly higher. This appears to stem from the differences in morphology among their sample and ours, as well as the different ranges of between the two samples. While the intercepts for the relationships traced by the active galaxy sample show good agreement with the quiescent galaxies, 2/3 of the galaxies in the Kormendy & Ho (2013) sample are ellipticals. Thus, the relationships for their sample are very similar to the relationships, because 2/3 of the points between them are exactly the same. On the other hand, the active galaxy sample is dominated by later-type galaxies where the bulge contributes a smaller fraction of the integrated galaxy light, and so the best-fit and relationships that we find for the active galaxies are quite different from each other.
We looked at the bulge-to-total ratios for the active galaxy sample and investigated whether splitting the sample into “early” (B/T ¿ 0.5) and “late” (B/T ¡ 0.5) types uncovered any offsets or separations among the sample that may lead to better agreement with the quiescent galaxy sample. The only obvious difference between these two subsamples is that the “early” types have more massive black holes than the “late” types, and so a cut in B/T is similar to a cut in and does not improve the agreement. As before, we also investigated the location of the Läsker et al. (2016) megamasers and find that they are wholly contained within the -band scatter of the active galaxy sample. If we again combine the active, quiescent, and megamaser samples as before, we find best fits of:
in the band, with a typical scatter of dex, and
in the band, with a typical scatter of dex. Thus, while the scatter is significantly increased when galaxies of all morphological types are treated equally, it is likely more representative of the true uncertainty on black hole mass estimates from the total galaxy luminosity.
5.3. Black Hole Mass – Bulge Stellar Mass Relationship
The relationship between black hole mass and bulge stellar mass is expected to be the physical basis for the relationship, where bulge light traces mass. A variety of methods have been used to investigate this relationship in the past, often with the aim of decoupling the relationship from any dependence on the relationship so they can be studied independently.
For example, Magorrian et al. (1998) carried out axisymmetric dynamical models to constrain the bulge mass and the black hole mass simultaneously. Marconi & Hunt (2003) measured effective bulge radii from 2MASS imaging for quiescent galaxies with dynamical black hole masses. The bulge radii were combined with to predict under the assumption that bulges behave similarly to isothermal spheres. Häring & Rix (2004), on the other hand, numerically solved the spherical Jeans equation while matching published luminosity and profiles for quiescent galaxies with dynamical black hole masses.
We can examine this relationship for active galaxies by estimating the bulge stellar mass from its opticalnear-infrared color and the prescriptions described above. The best-fit relationship between the black hole mass and the stellar mass of the bulge, based on the Bell & de Jong (2001) predictions, is found to be:
with a typical scatter of dex.
If we estimate using the prescriptions of Into & Portinari (2013), we find the best fit to be:
with a typical scatter of dex. These relationships are displayed in Figure 5.
For a direct comparison with the quiescent galaxy sample, we recalculated the bulge masses based on the absolute magnitudes of the bulges and the colors tabulated by Kormendy & Ho (2013) with the prescriptions of both Bell & de Jong (2001) and Into & Portinari (2013). Because Kormendy & Ho (2013) only provide an integrated color for each galaxy, we note that we would expect there to be a bias in the bulge masses derived for the disk galaxies in their sample because of the different colors of bulges and disks. The best-fit relationships for the quiescent galaxies are shown as the black dashed lines in Figure 5.
While the active galaxy sample displays a linear relationship between and bulge stellar mass, the quiescent galaxy relationships are quite a bit steeper. The two samples agree better using the prescriptions of Bell & de Jong (2001), although both prescriptions show agreement between the samples at the low mass end.
The megamaser sample of Läsker et al. (2016) reports based on near-infrared HST and ground-based imaging and the prescriptions of Bell et al. (2003), which allows for a simple comparison with our results. We again find that all 9 megamasers are contained wholly within the scatter of the active galaxy sample, with no apparent offsets in bulge mass or black hole mass.
Noting that there is good agreement between the active, quiescent, and megamaser samples, we also fit the relationship with all three samples combined. Assuming the Bell & de Jong (2001) prescriptions, the best-fit relationship is:
with a scatter of () dex. The red long-short dashed line in the left panel of Figure 5 displays this fit. While the slope is quite a bit steeper than that found for the active galaxy sample, they only disagree at the level. Furthermore, there is good agreement with the gray shaded region (which denotes the uncertainty on the fit to the active sample alone) over the range sampled by the active galaxies. This appears to indicate that all three samples may be drawn from the same parent population of galaxies.
We also compared our results to those of simulated galaxies. Some caution must be taken when interpreting such comparisons, because cosmological galaxy simulations are generally tuned to match a set of observables. For example, the slope of the relationship is not expected to be affected by such tuning, but the intercept is. Furthermore, there is no agreement on the best way to separate the bulges of late-type galaxies from their disks in simulations, where the resolution is often a limiting factor, so the simulated galaxies are either compared to samples of massive early-type galaxies where (e.g., Steinborn et al. 2015; Schaye et al. 2015) or a prescription is applied to estimate the bulge contribution.
Sijacki et al. (2015) used the high-resolution hydrodynamical Illustris simulations to explore the predicted relationship for galaxies. The total stellar mass within the stellar half-mass radius was used as a proxy for the bulge mass. This simplification does not take into account different bulge mass fractions of galaxies, nor the fact that some galaxies may not have a bulge at all. Additionally, the Illustris simulations assumed a Chabrier IMF, which can be compared to a “diet” Salpeter like that employed by Bell & de Jong (2001) by adding 0.093 dex (Gallazzi et al., 2008). To compare with a Kroupa IMF like that employed by Into & Portinari (2013), on the other hand, we subtracted 0.057 dex (Bell et al., 2003; Herrmann et al., 2016). The best-fit relationship of Sijacki et al. (2015), with the IMF scaled appropriately, is displayed as the blue long-dashed lines in Figure 5. With a reported slope of , it is in good agreement with our findings, especially when we adopt the Bell & de Jong (2001) prescriptions.
The large-volume Horizon-AGN simulations, which adopt a Salpeter IMF, were analyzed by Volonteri et al. (2016). To separate the bulge contribution, they tried various prescriptions, including examining the kinematics and also adopting a double Sérsic model for each galaxy, where the indices for the two Sérsic profiles were chosen to be [1.0, 1.0], [1.0, 4.0], or [1.0, 1.0 or 4.0]. The slope of the relationship based on these various prescriptions ranges from , which is in good agreement with our findings for the active galaxies using either the Bell & de Jong (2001) or Into & Portinari (2013) prescriptions, although it is somewhat in tension with our results for the combined active, quiescent, and megamaser samples. This tension may result from incompleteness in the Horizon-AGN simulation for black holes with M, which is the region probed by many of the active galaxies.
5.4. Black Hole Mass – Galaxy Stellar Mass Relationship
In the same way, we can examine the best-fit relationship between the black hole mass and the total stellar mass of the galaxy. When we adopt the Bell & de Jong (2001) predictions, we find a best fit of:
with a typical scatter of dex.
If we instead estimate using the prescriptions of Into & Portinari (2013), we find the best fit to be:
with a typical scatter of dex. These relationships are displayed in Figure 6.
Interestingly, the relationship based on the Into & Portinari (2013) values is similar to that found by Reines & Volonteri (2015) for inactive black holes residing in ellipticals and classical bulges. This would seem to contradict their finding that active galaxies lie below the relationship defined by local quiescent galaxies, although a direct comparison is somewhat difficult given that they used prescriptions of Zibetti et al. (2009), who employ a different initial mass function than Into & Portinari (2013).
We therefore recalculated the relationship for local quiescent galaxies based on the absolute magnitudes and the colors tabulated by Kormendy & Ho (2013), using both the Bell & de Jong (2001) and Into & Portinari (2013) prescriptions for direct comparison with the active galaxies in our sample. The best fits are shown as the dashed lines in Figure 6. Using the Bell & de Jong (2001) prescription, we find a nearly identical slope for the quiescent galaxies compared to the active galaxies, but an intercept that is 0.75 dex higher, supporting the findings of Reines & Volonteri (2015) that active galaxies fall below quiescent galaxies in this parameter space. However, using the Into & Portinari (2013) prescription instead, we find a slightly steeper slope for the quiescent galaxies which, when coupled with the intercept, show the two samples to be in general agreement at the low-mass end while diverging at the high mass end.
If we again combine the active sample with the quiescent galaxies and the megamasers, we find a best-fit relationship of
with a scatter of () dex. This fit is denoted with the red long-short dashed line in the left panel of Figure 6. Once again, the consistency with the results derived solely from the active galaxies seems to indicate that all of these subsamples may be drawn from the same parent population.
Recently, Shankar et al. (2016) investigated the potential for selection bias among the quiescent galaxy sample using Monte Carlo simulations and a large sample of galaxies drawn from the Sloan Digital Sky Survey. They concluded that the quiescent galaxy sample is selected from an upper “ridgeline” in the distribution of normal galaxy properties, leading to a bias of a factor of in the normalization of the relationship. If such a bias exists, that would argue against our choice of scaling factor for reverberation-based masses in the active galaxy sample, and would instead argue for . Interestingly, when we adopt for the scaling of reverberation-based , and examine the relationship based on the prescriptions of Bell & de Jong (2001), we find that the active galaxy sample closely follows the predicted unbiased relationship in Equation 6 of Shankar et al. (2016). The stellar masses predicted by Into & Portinari (2013), however, are undermassive compared to the predicted relationship, even when accounting for the slight differences in assumed IMF. However, it appears that bars will affect the measurements of effective radii (Meert et al., 2015) and possibly velocity dispersion (Batiste et al., 2017) that are adopted for the “unbiased” SDSS sample considered by Shankar et al. (2016). These effects will be strongest at the low-mass end, where most of the active galaxies in our sample are found, thus complicating the interpretation for reverberation-based masses. Furthermore, forward modeling of velocity-resolved reverberation signals by Pancoast et al. (2014) and Grier et al. (2017) has constrained the geometry and kinematics of the broad line region and the black hole mass for 9 AGNs, independent of any factor. Both studies recover modeling-based black hole masses that agree well with values derived from traditional reverberation analysis and the use of (as described in Section 4.5). These findings argue against the use of for the proper scaling of reverberation masses, but do not rule out that there may be biases present in the quiescent galaxy sample.
Unlike for the relationship, comparisons with simulated galaxies are much simpler when the entire galaxy stellar mass is used because the issues with bulge-disk decompositions are avoided, although the caveats related to the tuning of parameters in the simulations remain. Mutlu-Pakdil et al. (2018) recently examined the relationships between black holes and large-scale galaxy properties for spiral galaxies in the Illustris simulations. Using the same IMF corrections described in the previous section, we compared our best-fit relationships to theirs (blue long-dashed lines in Figure 6) and we find excellent agreement, especially when we adopt the Bell & de Jong (2001) prescriptions, which may argue against any potential bias in the reverberation-based scaling. Volonteri et al. (2016) examined the relationship for galaxies from the Horizon-AGN simulation, and find a slope that is somewhat shallower than we have found, although the low-mass end of their relationship may be biased by incompleteness. Steinborn et al. (2015) used the Magneticum Pathfinder Simulations to examine the relationship, excluding simulated galaxies for which M. Perhaps unsurprisingly, their reported best-fit relationship (with a slope of 1.09) agrees with the most massive black holes in the active galaxy sample ( M), but diverges at lower black hole masses, predicting a larger at fixed , similar to the findings of Volonteri et al. (2016).
Many large photometric surveys that are currently in operation or are upcoming will collect photometry in multiple filters. When considering that these surveys that may need to be treated in an automated fashion, the stellar mass of the galaxy based on its color appears to be a better predictor of black hole mass than the total galaxy luminosity in a single filter. This can be seen from the decreased scatter in the relationship for the combined active, quiescent, and megamaser samples ( dex) relative to the relationship ( dex).
5.5. Black Hole Mass Fraction
Finally, we also investigated the typical fraction of black hole mass to stellar mass for each galaxy. We find a median value of , however we also find a very clear relationship between and , while there appears to be no obvious relationship between and (see Figure 7).
For comparison, we derived the black hole mass fractions for the quiescent galaxy sample of Kormendy & Ho (2013) and find a median value of . At first glance, this would appear to demonstrate that active galaxies host undermassive black holes compared to quiescent galaxies. However, the samples cover different ranges, with the active galaxy sample skewed toward lower , while the quiescent galaxy sample is skewed to higher (see Figure 8), and seems to depend strongly on . For a better comparison between the two samples, we binned the galaxies in each sample by with bins of width 0.5 dex. For each bin with three or more objects, we computed the median black hole mass to stellar mass fraction. Figure 8 shows the median as a function of for the two samples, with the active sample in red and the quiescent sample in black. The majority of the overlap between the samples exists within the range , with the range extending to lower black hole masses in the active galaxy sample, and extending to higher black hole masses in the quiescent galaxy sample. We have adopted based on the predictions of Bell & de Jong (2001) in Figure 8, but while the values slightly change, the overall trend is the same if we adopt based on the Into & Portinari (2013) values. The two samples show broad agreement, both in the overall trend – with more massive black holes comprising larger mass fractions of their galaxies – and with the typical values for the black hole mass fraction at a fixed value of . While there seems to be a tendency for the active galaxies to lie slightly below the quiescent galaxies in the expected black hole mass fraction at a fixed black hole mass, the values agree within the standard deviation for each bin, and the small and uneven number of objects in each bin make it difficult to draw firm conclusions about any apparent offset between the two samples. Notably, the trend appears to continue across the full range of black hole masses probed by either sample.
We also examined the megamaser sample of Läsker et al. (2016) for comparison. Adopting the same bins for the megamaser sample as for the above two samples, we show the median in blue in Figure 8. There is no apparent offset between the megamaser sample and the reverberation sample, nor with the extension of the quiescent sample to lower black hole masses. Läsker et al. (2016) noted in their study that the megamaser galaxies appeared to probe a lower at fixed galaxy mass than the reverberation sample (as reported by Bentz et al. 2009a), but this discrepancy has been completely erased with the larger sample and extended range of and galaxy properties presented here.
The scaling of as a function of was previously noticed by Trakhtenbrot & Netzer (2010). Using large samples of local non-AGN galaxies and AGN galaxies at and scaling relationships to predict and , they found that . A formal fit to the active, quiescent, and maser galaxies examined here finds:
with a typical scatter of dex, which agrees well with a formal fit to the active galaxies alone, and to the estimated slope reported by Trakhtenbrot & Netzer (2010).
Interestingly, we find the same scaling between and among simulated galaxies from Illustris. Vogelsberger et al. (2014) provide black hole masses and galaxy stellar masses for two subsamples of representative “red” and “blue” galaxies from the Illustris simulation. The “blue” galaxies preferentially occupy the lower range that is probed here by the active and megamaser samples, and the “red” galaxies preferentially occupy the upper range probed by the quiescent galaxies. The scaling in as a function of in the simulated galaxies matches the observed galaxies extremely well in both slope and offset.
It is clear from these studies that the commonly-used assumption of a constant is incorrect in the local universe and possibly up to . Attempts to search for cosmic evolution of black holes and host galaxies, or to search for differences in the evolutionary paths of distinct galaxy samples, should be careful to account for this scaling when the samples are not matched in .
We suggest that the physical meaning of this scaling may be related to differences in feedback efficiency as a function of galaxy mass. High-resolution and zoom-in simulations of individual galaxies show that supernova feedback is extremely effective at prohibiting black hole growth at early times (e.g., Dubois et al. 2015; Trebitsch et al. 2018; Anglés-Alcázar et al. 2017). Once the host galaxy reaches a critical mass (; Dubois et al. 2015; Anglés-Alcázar et al. 2017), supernova feedback can no longer restrict the gas flow to the nucleus and the black hole will undergo a period of rapid growth, effectively “catching up” with the galaxy. This period of rapid growth is short-lived, however, because AGN feedback soon becomes important and the black hole then regulates its own growth and the continued growth of the galaxy (e.g., Dubois et al. 2015; McAlpine et al. 2017). In this scenario, we may currently be witnessing the rapid growth phase for low-mass black holes in the local universe.
Using high-resolution optical HST images and deep, ground-based near-infrared images, we have constrained the photometric properties of 37 active galaxies hosting black holes with reverberation-based measurements. We have compared our results with those of megamaser galaxies and of quiescent galaxies with black hole masses from dynamical modeling, and we have re-examined several black hole-galaxy scaling relationships. In general, we find that megamasers behave as a subset of the active galaxy sample, and there is evidence that the active and megamaser samples may be drawn from the same parent population as the quiescent galaxies. We also find the following:
The relationship for active galaxies is slightly steeper in the near-infrared than the optical, and both bandpasses exhibit similar scatter. There is general agreement with our results and those found for quiescent galaxies by Kormendy & Ho (2013) and the megamaser sample of Läsker et al. (2016). is found to have the tightest correlation with of the relationships examined here, and will provide the least biased estimates from photometry.
The relationship for active galaxies is only slightly less well defined than the relationship, but when combined with the megamaser and quiescent galaxy samples, the scatter increases significantly. Large photometric surveys may forego bulge-disk decompositions and estimate unbiased black hole masses more quickly with total galaxy luminosity, ignoring galaxy morphology, but with a loss of accuracy.
The relationship for active galaxies is linear, while the quiescent galaxy sample displays a steeper slope. Both samples agree at the low mass end, and the agreement is better when the prescriptions of Bell & de Jong (2001) are used rather than those of Into & Portinari (2013). The best-fit relationship for the combined active, megamaser, and quiescent samples agrees well with the relationship for the active galaxies alone, which also agrees well with the expectations from the high-resolution Illustris hydrodynamical simulations. Agreement with other simulations is less clear because of incompleteness at M.
The active galaxy relationship tends to lie slightly below that of the quiescent galaxy sample, but there is excellent agreement with the best fit for the combined active, quiescent, and megamaser samples. There is also excellent agreement between the best-fit relationship for the active galaxies and the expectations from the high-resolution Illustris hydrodynamical simulations, but incompleteness affects comparsions with other simulations. Large photometric surveys with multiple filters will achieve better accuracy in predicted black hole masses using the stellar mass of the galaxy (based on a two filter color) than the galaxy luminosity in a single filter.
The fraction of the black hole mass to the galaxy stellar mass is a strong function of black hole mass (but not stellar mass), with more massive black holes occupying larger fractions of . The same trend is seen in the quiescent galaxy and megamaser samples, and the median black hole mass fractions at fixed black hole mass are similar between all three samples. The median value of the black hole mass fraction ranges from % at M to % at M and follows the form . Studies that seek to compare different galaxy samples should be careful to account for this effect if the samples are not matched in black hole mass.
- Anglés-Alcázar et al. (2017) Anglés-Alcázar, D., Faucher-Giguère, C.-A., Quataert, E., Hopkins, P. F., Feldmann, R., Torrey, P., Wetzel, A., & Kereš, D. 2017, MNRAS, 472, L109
- Batiste et al. (2017) Batiste, M., Bentz, M. C., Manne-Nicholas, E. R., Onken, C. A., & Bershady, M. A. 2017, ApJ, 835, 271
- Bell & de Jong (2001) Bell, E. F., & de Jong, R. S. 2001, ApJ, 550, 212
- Bell et al. (2003) Bell, E. F., McIntosh, D. H., Katz, N., & Weinberg, M. D. 2003, ApJS, 149, 289
- Bentz et al. (2013) Bentz, M. C., Denney, K. D., Grier, C. J., et al. 2013, ApJ, 767, 149
- Bentz & Katz (2015) Bentz, M. C., & Katz, S. 2015, PASP, 127, 67
- Bentz et al. (2009a) Bentz, M. C., Peterson, B. M., Netzer, H., Pogge, R. W., & Vestergaard, M. 2009a, ApJ, 697, 160
- Bentz et al. (2009b) Bentz, M. C., Peterson, B. M., Pogge, R. W., & Vestergaard, M. 2009b, ApJ, 694, L166
- Blandford & McKee (1982) Blandford, R. D., & McKee, C. F. 1982, ApJ, 255, 419
- Bonoli et al. (2009) Bonoli, S., Marulli, F., Springel, V., et al. 2009, MNRAS, 396, 423
- Cox (2000) Cox, A. N. 2000, Allen’s astrophysical quantities (New York: Springer, 4th ed.)
- de Vaucouleurs (1948) de Vaucouleurs, G. 1948, Annales d’Astrophysique, 11, 247
- Degraf et al. (2011) Degraf, C., Di Matteo, T., & Springel, V. 2011, MNRAS, 413, 1383
- DeGraf et al. (2015) DeGraf, C., Di Matteo, T., Treu, T., et al. 2015, MNRAS, 454, 913
- Dubois et al. (2015) Dubois, Y., Volonteri, M., Silk, J., Devriendt, J., Slyz, A., & Teyssier, R. 2015, MNRAS, 452, 1502
- Ferrarese & Merritt (2000) Ferrarese, L., & Merritt, D. 2000, ApJ, 539, L9
- Gallazzi et al. (2008) Gallazzi, A., Brinchmann, J., Charlot, S., & White, S. D. M. 2008, MNRAS, 383, 1439
- Gebhardt et al. (2000) Gebhardt, K., Bender, R., Bower, G., et al. 2000, ApJ, 539, L13
- Grier et al. (2013) Grier, C. J., Martini, P., Watson, L. C., et al. 2013, ApJ, 773, 90
- Grier et al. (2017) Grier, C. J., Pancoast, A., Barth, A. J., Fausnaugh, M. M., Brewer, B. J., Treu, T., & Peterson, B. M. 2017, ApJ, 849, 146
- Gültekin et al. (2009) Gültekin, K., Richstone, D. O., Gebhardt, K., et al. 2009, ApJ, 698, 198
- Hamilton (2014) Hamilton, T. S. 2014, in American Astronomical Society Meeting Abstracts, Vol. 223, American Astronomical Society Meeting Abstracts #223, 145.02
- Häring & Rix (2004) Häring, N., & Rix, H.-W. 2004, ApJ, 604, L89
- Heckman & Best (2014) Heckman, T. M., & Best, P. N. 2014, ARA&A, 52, 589
- Herrmann et al. (2016) Herrmann, K. A., Hunter, D. A., Zhang, H.-X., & Elmegreen, B. G. 2016, AJ, 152, 177
- Hopkins et al. (2010) Hopkins, P. F., Bundy, K., Croton, D., et al. 2010, ApJ, 715, 202
- Into & Portinari (2013) Into, T., & Portinari, L. 2013, MNRAS, 430, 2715
- Jarrett et al. (2000) Jarrett, T. H., Chester, T., Cutri, R., et al. 2000, AJ, 119, 2498
- Jarrett (2000) Jarrett, T. H. 2000, PASP, 112, 1008
- Kaviraj et al. (2011) Kaviraj, S., Schawinski, K., Silk, J., & Shabala, S. S. 2011, MNRAS, 415, 3798
- Kelly (2007) Kelly, B. C. 2007, ApJ, 665, 1489
- Kelly & Merloni (2012) Kelly, B. C., & Merloni, A. 2012, Advances in Astronomy, 2012, 970858
- Kim et al. (2011) Kim, J.-h., Wise, J. H., Alvarez, M. A., & Abel, T. 2011, ApJ, 738, 54
- Kinney et al. (1996) Kinney, A. L., Calzetti, D., Bohlin, R. C., et al. 1996, ApJ, 467, 38
- Kormendy & Ho (2013) Kormendy, J., & Ho, L. C. 2013, ARA&A, 51, 511
- Kormendy & Kennicutt (2004) Kormendy, J., & Kennicutt, Jr., R. C. 2004, ARA&A, 42, 603
- Kormendy & Richstone (1995) Kormendy, J., & Richstone, D. 1995, ARA&A, 33, 581
- Krist (1993) Krist, J. 1993, in ASP Conf. Ser. 52: Astronomical Data Analysis Software and Systems II, 536
- Lapi et al. (2014) Lapi, A., Raimundo, S., Aversa, R., et al. 2014, ApJ, 782, 69
- Läsker et al. (2014a) Läsker, R., Ferrarese, L., & van de Ven, G. 2014a, ApJ, 780, 69
- Läsker et al. (2016) Läsker, R., Greene, J.E., Seth, A., et al. 2016, ApJ, 825, 3
- Lippai et al. (2009) Lippai, Z., Frei, Z., & Haiman, Z. 2009, ApJ, 701, 360
- Magorrian et al. (1998) Magorrian, J., Tremaine, S., Richstone, D., et al. 1998, AJ, 115, 2285
- Marconi & Hunt (2003) Marconi, A., & Hunt, L. K. 2003, ApJ, 589, L21
- Mathur et al. (2012) Mathur, S., Fields, D., Peterson, B. M., & Grupe, D. 2012, ApJ, 754, 146
- McAlpine et al. (2017) McAlpine, S., Bower, R. G., Harrison, C. M., Crain, R. A., Schaller, M., Schaye, J., & Theuns, T. 2017, MNRAS, 468, 3395
- Meert et al. (2015) Meert, A., Vikram, V., & Bernardi, M. 2015, MNRAS, 446, 3943
- Miller et al. (2006) Miller, L., Percival, W. J., Croom, S. M., & Babić, A. 2006, A&A, 459, 43
- Mutlu-Pakdil et al. (2018) Mutlu-Pakdil, B., Seigar, M. S., Hewitt, I. B., et al. 2018, MNRAS, 474, 2594
- Oogi et al. (2016) Oogi, T., Enoki, M., Ishiyama, T., et al. 2016, MNRAS, 456, L30
- Osterbrock & Pogge (1985) Osterbrock, D. E., & Pogge, R. W. 1985, ApJ, 297, 166
- Ostriker et al. (2010) Ostriker, J. P., Choi, E., Ciotti, L., Novak, G. S., & Proga, D. 2010, ApJ, 722, 642
- Pancoast et al. (2014) Pancoast, A., Brewer, B. J., Treu, T., et al. 2014, ApJ, 445, 3073
- Peng et al. (2002) Peng, C. Y., Ho, L. C., Impey, C. D., & Rix, H. 2002, AJ, 124, 266
- Peng et al. (2010) Peng, C. Y., Ho, L. C., Impey, C. D., & Rix, H.-W. 2010, AJ, 139, 2097
- Peterson (1993) Peterson, B. M. 1993, PASP, 105, 247
- Peterson et al. (2004) Peterson, B. M., Ferrarese, L., Gilbert, K. M., et al. 2004, ApJ, 613, 682
- Reines & Volonteri (2015) Reines, A. E., & Volonteri, M. 2015, ApJ, 813, 82
- Schaye et al. (2015) Schaye, J., Crain, R. A., Bower, R.G., et al. 2015, MNRAS, 446, 521
- Schlafly & Finkbeiner (2011) Schlafly, E. F., & Finkbeiner, D. P. 2011, ApJ, 737, 103
- Schlegel et al. (1998) Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 525
- Sérsic (1968) Sérsic, J. L. 1968, Atlas de Galaxias Australes (Cordoba, Argentina: Observatorio Astronomico, 1968)
- Shabala et al. (2011) Shabala, S. S., Kaviraj, S., & Silk, J. 2011, MNRAS, 413, 2815
- Shankar et al. (2016) Shankar, F., Bernardi, M., Sheth, R. K., et al. 2016, MNRAS, 460, 3119
- Shirakata et al. (2016) Shirakata, H., Kawaguchi, T., Okamoto, T., et al. 2016, MNRAS, 461, 4389
- Sijacki et al. (2015) Sijacki, D., Vogelsberger, M., Genel, S., et al. 2015, MNRAS, 452, 575
- Skrutskie et al. (2006) Skrutskie, M. F., Cutri, R. M., Stiening, R., et al. 2006, AJ, 131, 1163
- Steinborn et al. (2015) Steinborn, L. K., Dolag, K., Hirschmann, M., Prieto, M. A., & Remus, R.-S. 2015, MNRAS, 448, 1504
- Tonry et al. (2001) Tonry, J. L., Dressler, A., Blakeslee, J. P., et al. 2001, ApJ, 546, 681
- Trakhtenbrot & Netzer (2010) Trakhtenbrot, B., & Netzer, H. 2010, MNRAS, 406, L35
- Trebitsch et al. (2018) Trebitsch, M., Volonteri, M., Dubois, Y., & Madau, P. 2018, MNRAS, 478, 5620
- Tully et al. (2009) Tully, R. B., Rizzi, L., Shaya, E. J., et al. 2009, AJ, 138, 323
- Tully et al. (2008) Tully, R. B., Shaya, E. J., Karachentsev, I. D., et al. 2008, ApJ, 676, 184
- van Dokkum (2001) van Dokkum, P. G. 2001, PASP, 113, 1420
- Veilleux et al. (2009) Veilleux, S., Kim, D.-C., Rupke, D. S. N., et al. 2009, ApJ, 701, 587
- Vogelsberger et al. (2014) Vogelsberger, M., Genel, S., Springel, V., et al. 2014, MNRAS, 444, 1518
- Volonteri & Natarajan (2009) Volonteri, M., & Natarajan, P. 2009, MNRAS, 400, 1911
- Volonteri et al. (2016) Volonteri, M., Dubois, P., Pichon, C., & Devriendt, J. 2016, MNRAS, 460, 2979
- Wandel (2002) Wandel, A. 2002, ApJ, 565, 762
- Worthey (1994) Worthey, G. 1994, ApJS, 95, 107
- Zibetti et al. (2009) Zibetti, S., Charlot, S., & Rix, H.-W. 2009, MNRAS, 400, 1181
|Object||RA||Dec||Date||Exp Time||Obs Setup|
|Mrk 335||0.0258||350.0||WIYN WHIRC H|
|2040.0||ACS HRC F550M|
|Mrk 1501||0.0893||1400.0||WIYN WHIRC H|
|2236.0||WFC3 UVIS2 F547M|
|PG 0026+129||0.1420||2559.8||NICMOS NIC2 F160W|
|Mrk 590||0.0264||1500.0||WIYN WHIRC H|
|1020.0||ACS HRC F550M|
|3C 120||0.0330||200 .0||WIYN WHIRC H|
|1020.0||ACS HRC F550M|
|Akn 120||0.0327||1000.0||WIYN WHIRC H|
|2040.0||ACS HRC F550M|
|Mrk 6||0.0188||720.0||WIYN WHIRC H|
|2620.0||WFC3 UVIS2 F547M|
|Mrk 79||0.0222||4140.0||WIYN WHIRC H|
|2040.0||ACS HRC F550M|
|PG 0844+349||0.0640||2559.8||NICMOS NIC2 F160W|
|1020.0||ACS HRC F550M|
|Mrk 110||0.0353||3500.0||WIYN WHIRC H|
|1020.0||ACS HRC F550M|
|NGC 3227||0.0039||1470.0||WIYN WHIRC H|
|2250.0||WFC3 UVIS2 F547M|
|NGC 3516||0.0088||1625.0||WIYN WHIRC H|
|2660.0||WFC3 UVIS2 F547M|
|SBS 1116+583A||0.0279||6300.0||WIYN WHIRC H|
|2510.0||WFC3 UVIS2 F547M|
|Arp 151||0.0211||3780.0||WIYN WHIRC H|
|2450.0||WFC3 UVIS2 F547M|
|Mrk 1310||0.0196||4500.0||WIYN WHIRC H|
|2240.0||WFC3 UVIS2 F547M|
|NGC 4051||0.0023||3060.0||WIYN WHIRC H|
|2340.0||WFC3 UVIS2 F547M|
|NGC 4151||0.0033||1005.0||WIYN WHIRC H|
|2310.0||WFC3 UVIS2 F547M|
|Mrk 202||0.0210||4800.0||WIYN WHIRC H|
|2510.0||WFC3 UVIS2 F547M|
|NGC 4253||0.0129||1700.0||WIYN WHIRC H|
|2270.0||WFC3 UVIS2 F547M|
|PG 1226+023||0.1583||2250.0||WIYN WHIRC H|
|2040.0||ACS HRC F550M|
|PG 1229+204||0.0630||2559.8||NICMOS NIC2 F160W|
|2040.0||ACS HRC F550M|
|NGC 4593||0.0090||960.0||WIYN WHIRC H|
|2240.0||WFC3 UVIS2 F547M|
|NGC 4748||0.0146||3600.0||WIYN WHIRC H|
|2250.0||WFC3 UVIS2 F547M|
|PG 1307+085||0.1550||2559.8||NICMOS NIC2 F160W|
|Mrk 279||0.0305||2300.0||WIYN WHIRC H|
|1020.0||ACS HRC F550M|
|PG 1411+442||0.0896||2559.8||NICMOS NIC2 F160W|
|2040.0||ACS HRC F550M|
|PG 1426+015||0.0866||2559.8||NICMOS NIC2 F160W|
|Mrk 817||0.0315||3520.0||WIYN WHIRC H|
|1020.0||ACS HRC F550M|
|PG 1613+658||0.1290||1260.0||WIYN WHIRC H|
|2040.0||ACS HRC F550M|
|PG 1617+175||0.1124||2559.8||NICMOS NIC2 F160W|
|PG 1700+518||0.2920||2559.8||NICMOS NIC2 F160W|
|2040.0||ACS HRC F550M|
|3C 390.3||0.0561||3060.0||WIYN WHIRC H|
|1020.0||ACS HRC F550M|
|Zw 229-015||0.0279||2000.0||WIYN WHIRC H|
|2320.0||WFC3 UVIS2 F547M|
|NGC 6814||0.0052||1200.0||WIYN WHIRC H|
|2240.0||WFC3 UVIS2 F547M|
|Mrk 509||0.0344||385.0||WIYN WHIRC H|
|PG 2130+099||0.0630||1800.0||WIYN WHIRC H|
|1020.0||ACS HRC F550M|
|NGC 7469||0.0163||300.0||WIYN WHIRC H|
|2240.0||WFC3 UVIS2 F547M|
|(mag)||(mag)||(arcsec)||( E of N)|
|Mrk 1501||inner disk|
|dust lane - inner|
|ring - inner|
|ring - outer|
|ring - inner|
|ring - outer|