Since the discovery of active galactic nuclei (AGN) and their subclasses, a unification scheme of AGN has been long sought. Orientation-based unified models predict that some of the diversity within AGN subclasses can be explained by the different viewing angles of the observer. Several orientation categorisations have been suggested, but a widely applicable measure has yet to be found. Using the properties of the ultraviolet and optical broad emission lines of quasars, in particular the velocity offsets and line widths of high-ionisation CIV and low-ionisation MgII lines, a correlation has been measured. It is postulated that this correlation is due to the viewing angle of the observer. Comparison with other orientation tracers shows consistency with this interpretation. Using a simulation of a wide angle disk-wind model for the broad emission line region, we successfully replicate the observed correlation with inclination. Future more detailed modelling will not only enable improved accuracy in the determination of the viewing angle to individual AGN, but will also substantially increase our understanding of the emitting regions of AGN.

Microlensing by stars within distant galaxies acting as strong gravitational lenses of multiply-imaged quasars, provides a unique and direct measurement of the internal structure of the lensed quasar on nano-arcsecond scales. The measurement relies on the temporal variation of high-magnification caustic crossings which vary on timescales of days to years. Multiwavelength observations provide information from distinct emission regions in the quasar. Through monitoring of these strong gravitational lenses, a full tomographic view can emerge with Astronomical-Unit scale resolution. Work to date has demonstrated the potential of this technique in about a dozen systems. In the 2020s there will be orders of magnitude more systems to work with. Monitoring of lens systems for caustic-crossing events to enable triggering of multi-platform, multi-wavelength observations in the 2020s will fulfill the potential of quasar microlensing as a unique and comprehensive probe of active black hole structure and dynamics.

Optical and X-ray observations of strongly gravitationally lensed quasars (especially when four separate images of the quasar are produced) determine not only the amount of matter in the lensing galaxy but also how much is in a smooth component and how much is composed of compact masses (e.g., stars, stellar remnants, primordial black holes, CDM sub-halos, and planets). Future optical surveys will discover hundreds to thousands of quadruply lensed quasars, and sensitive X-ray observations will unambiguously determine the ratio of smooth to clumpy matter at specific locations in the lensing galaxies and calibrate the stellar mass fundamental plane, providing a determination of the stellar $M/L$. A modest observing program with a sensitive, sub-arcsecond X-ray imager, combined with the planned optical observations, can make those determinations for a large number (hundreds) of the lensing galaxies, which will span a redshift range of $\sim$$0.25<z<1.5$

The stellar mass binary black hole (sBBH) mergers presently detected by LIGO may originate wholly or in part from binary black hole mergers embedded in disks of gas around supermassive black holes. Determining the contribution of these active galactic nucleus (AGN) disks to the sBBH merger rate enables us to uniquely measure important parameters of AGN disks, including their typical density, aspect ratio, and lifetime, thereby putting unique limits on an important element of galaxy formation. For the first time, gravitational waves will allow us to reveal the properties of the hidden interior of AGN disks, while electromagnetic radiation (EM) probes the disk photosphere. The localization of sBBH merger events from LIGO is generally insufficient for association with a single EM counterpart. However, the contribution to the LIGO event rate from rare source types (such as AGNs) can be determined on a statistical basis. To determine the contribution to the sBBH rate from AGNs in the next decade requires: {\em 1) a complete galaxy catalog for the LIGO search volume, 2) strategic multi-wavelength EM follow-up of LIGO events and 3) significant advances in theoretical understanding of AGN disks and the behavior of objects embedded within them.}

We present new HST WFPC3 imaging of four gravitationally lensed quasars: MG 0414+0534; RXJ 0911+0551; B 1422+231; WFI J2026-4536. In three of these systems we detect wavelength-dependent microlensing, which we use to place constraints on the sizes and temperature profiles of the accretion discs in each quasar. Accretion disc radius is assumed to vary with wavelength according to the power-law relationship $r\propto \lambda^p$, equivalent to a radial temperature profile of $T\propto r^{-1/p}$. The goal of this work is to search for deviations from standard thin disc theory, which predicts that radius goes as wavelength to the power $p=4/3$. We find a wide range of power-law indices, from $p=1.4^{+0.5}_{-0.4}$ in B 1422+231 to $p=2.3^{+0.5}_{-0.4}$ in WFI J2026-4536. The measured value of $p$ appears to correlate with the strength of the wavelength-dependent microlensing. We explore this issue with mock simulations using a fixed accretion disc with $p=1.5$, and find that cases where wavelength-dependent microlensing is small tend to under-estimate the value of $p$. This casts doubt on previous ensemble single-epoch measurements which have favoured low values using samples of lensed quasars that display only moderate chromatic effects. Using only our systems with strong chromatic microlensing we prefer $p>4/3$, corresponding to shallower temperature profiles than expected from standard thin disc theory.

A key to understanding quasar unification paradigms is the emission properties of broad absorption line quasars (BALQs). The fact that only a small fraction of quasar spectra exhibit deep absorption troughs blueward of the broad permitted emission lines provides a crucial clue to the structure of quasar emitting regions. To learn whether it is possible to discriminate between the BALQ and non-BALQ populations given the observed spectral properties of a quasar, we employ two approaches: one based on statistical methods and the other supervised machine learning classification, applied to quasar samples from the Sloan Digital Sky Survey. The features explored include continuum and emission line properties, in particular the absolute magnitude, redshift, spectral index, line width, asymmetry, strength, and relative velocity offsets of high-ionisation CIV $\lambda 1549$ and low-ionisation MgII $\lambda 2798$ lines. We consider a complete population of quasars, and assume that the statistical distributions of properties represent all angles where the quasar is viewed without obscuration. The distributions of the BALQ and non-BALQ sample properties show few significant differences. None of the observed continuum and emission line features are capable of differentiating between the two samples. Most published narrow disk-wind models are inconsistent with these observations, and an alternative disk-wind model is proposed. The key feature of the proposed model is a disk-wind filling a wide opening angle with multiple radial streams of dense clumps.

Black hole mergers detectable with LIGO can occur in active galactic nucleus (AGN) disks. Here we parameterize the merger rates, the mass spectrum and the spin spectrum of black holes (BH) in AGN disks. The predicted merger rate spans $\sim 10^{-4}-10^{4} \rm{Gpc}^{-1} \rm{yr}^{-1}$, so upper limits from LIGO ($<212\rm{Gpc}^{-1}\rm{yr}^{-1}$) already constrain it. The predicted mass spectrum has the form of a broken power-law consisting of a pre-existing BH powerlaw mass spectrum and a harder powerlaw mass spectrum resulting from mergers. The predicted spin spectrum is multi-peaked with the evolution of retrograde spin BH in the gas disk playing a key role. We outline the large uncertainties in each of these LIGO observables for this channel and we discuss ways in which they can be constrained in the future.

We assess the contribution of dynamical hardening by direct three-body scattering interactions to the rate of stellar-mass black hole binary (BHB) mergers in galactic nuclei. We derive an analytic model for the single-binary encounter rate in a nucleus with spherical and disk components hosting a super-massive black hole (SMBH). We determine the total number of encounters $N_{\rm GW}$ needed to harden a BHB to the point that inspiral due to gravitational wave emission occurs before the next three-body scattering event. This is done independently for both the spherical and disk components. Using a Monte Carlo approach, we refine our calculations for $N_{\rm GW}$ to include gravitational wave emission between scattering events. For astrophysically plausible models we find that typically $N_{\rm GW} \lesssim$ 10. We find two separate regimes for the efficient dynamical hardening of BHBs: (1) spherical star clusters with high central densities, low velocity dispersions and no significant Keplerian component; and (2) migration traps in disks around SMBHs lacking any significant spherical stellar component in the vicinity of the migration trap, which is expected due to effective orbital inclination reduction of any spherical population by the disk. We also find a weak correlation between the ratio of the second-order velocity moment to velocity dispersion in galactic nuclei and the rate of BHB mergers, where this ratio is a proxy for the ratio between the rotation- and dispersion-supported components. Because disks enforce planar interactions that are efficient in hardening BHBs, particularly in migration traps, they have high merger rates that can contribute significantly to the rate of BHB mergers detected by the advanced Laser Interferometer Gravitational-Wave Observatory.

We present new spectroscopy of the z=3.62 gravitationally lensed quasar B1422+117 from the Gemini North GMOS integral field spectrograph. We observe significant differential magnifications between the broad emission lines and the continuum, as well as across the velocity structure of the Lyman-alpha line. We take advantage of this differential microlensing to algebraically decompose the quasar spectrum into the absorbed broad emission line and absorbed continuum components. We use the latter to derive the intrinsic Ly-alpha forest absorption spectrum. The proximity effect is clearly detected, with a proximity zone edge of 8.6-17.3 Mpc from the quasar, implying (perhaps intermittent) activity over at least 28 Myrs. The Ly-alpha line profile exhibits a blue excess that is inconsistent with a symmetric fit to the unabsorbed red side. This has important implications for the use of this fitting technique in estimating the absorbed blue Ly-alpha wings of Gunn-Peterson trough quasars.

The structure and kinematics of the broad line region (BLR) in quasars are still not well established. One popular BLR model is the disk-wind model that offers a geometric unification of a quasar based on the angle of viewing. We construct a simple kinematical disk-wind model with a narrow outflowing wind angle. The model is combined with radiative transfer in the Sobolev, or high velocity, limit. We examine how angle of viewing affects the observed characteristics of the emission line, especially the line widths and velocity offsets. The line profiles exhibit distinct properties depending on the orientation, wind opening angle, and region of the wind where the emission arises. At low inclination angle (close to face-on), we find the shape of the emission line is asymmetric with narrow width and significantly blueshifted. As the inclination angle increases (close to edge-on), the line profile becomes more symmetric, broader, and less blueshifted. Additionally, lines that arise close to the base of the disk wind, near the accretion disk, tend to be broad and symmetric. The relative increase in blueshift of the emission line with increasing wind vertical distance is larger for polar winds compared with equatorial winds. By considering the optical thickness of the wind, single-peaked line profiles are recovered for the intermediate and equatorial outflowing wind. The model is able to reproduce a faster response in either the red and blue sides of the line profile found in reverberation mapping studies. A quicker response in the red side is achieved in the model with a polar wind and intermediate wind opening angle at low viewing angle. The blue side response is faster for an equatorial wind seen at high inclination.

The Large Synoptic Survey Telescope is designed to provide an unprecedented optical imaging dataset that will support investigations of our Solar System, Galaxy and Universe, across half the sky and over ten years of repeated observation. However, exactly how the LSST observations will be taken (the observing strategy or "cadence") is not yet finalized. In this dynamically-evolving community white paper, we explore how the detailed performance of the anticipated science investigations is expected to depend on small changes to the LSST observing strategy. Using realistic simulations of the LSST schedule and observation properties, we design and compute diagnostic metrics and Figures of Merit that provide quantitative evaluations of different observing strategies, analyzing their impact on a wide range of proposed science projects. This is work in progress: we are using this white paper to communicate to each other the relative merits of the observing strategy choices that could be made, in an effort to maximize the scientific value of the survey. The investigation of some science cases leads to suggestions for new strategies that could be simulated and potentially adopted. Notably, we find motivation for exploring departures from a spatially uniform annual tiling of the sky: focusing instead on different parts of the survey area in different years in a "rolling cadence" is likely to have significant benefits for a number of time domain and moving object astronomy projects. The communal assembly of a suite of quantified and homogeneously coded metrics is the vital first step towards an automated, systematic, science-based assessment of any given cadence simulation, that will enable the scheduling of the LSST to be as well-informed as possible.

We re-examine the possibility of ultrahigh energy cosmic rays being accelerated in nearby dormant quasars. We particularize our study to heavy nuclei to accommodate the spectrum and nuclear composition recently reported by the Pierre Auger Collaboration. Particle acceleration is driven by the Blandford-Znajek mechanism, which wires the dormant spinning black holes as Faraday unipolar dynamos. We demonstrate that energy losses are dominated by photonuclear interactions on the ambient photon fields. We argue that the local dark fossils of the past quasar activity can be classified on the basis of how source parameters (mass of the central engine and photon background surrounding the accelerator) impact the photonuclear interaction. In this classification it is possible to distinguish two unequivocal type of sources: those in which nuclei are completely photodisintegrated before escaping the acceleration region and those in which photopion production is the major energy damping mechanism. We further argue that the secondary nucleons from the photodisintegrated nuclei (which have a steep spectral index at injection) can populate the energy region below "the ankle" feature in the cosmic ray spectrum, whereas heavy and medium mass nuclei (with a harder spectral index) populate the energy region beyond "the ankle", all the way to the high energy end of the spectrum. In addition, we show that five potential quasar remnants from our cosmic backyard correlate with the hot-spot observed by the Telescope Array.

We present new integral field spectroscopy of the gravitationally lensed broad absorption line (BAL) quasar H1413+117, covering the ultraviolet to visible rest-frame spectral range. We observe strong microlensing signatures in lensed image D, and we use this microlensing to simultaneously constrain both the broad emission and broad absorption line gas. By modeling the lens system over the range of probable lensing galaxy redshifts and using on a new argument based on the wavelength-independence of the broad line lensing magnifications, we determine that there is no significant broad line emission from smaller than ~20 light days. We also perform spectral decomposition to derive the intrinsic broad emission line (BEL) and continuum spectrum, subject to BAL absorption. We also reconstruct the intrinsic BAL absorption profile, whose features allow us to constrain outflow kinematics in the context of a disk-wind model. We find a very sharp, blueshifted onset of absorption of 1,500 km/s in both C IV and N V that may correspond to an inner edge of a disk-wind's radial outflow. The lower ionization Si IV and Al III have higher-velocity absorption onsets, consistent with a decreasing ionization parameter with radius in an accelerating outflow. There is evidence of strong absorption in the BEL component which indicates a high covering factor for absorption over two orders of magnitude in outflow radius.

We study a combined sample of 264 star-forming, 51 composite, and 73 active galaxies using optical spectra from SDSS and mid-infrared (mid-IR) spectra from the Spitzer Infrared Spectrograph. We examine optical and mid-IR spectroscopic diagnostics that probe the amount of star formation and relative energetic contributions from star formation and an active galactic nucleus (AGN). Overall we find good agreement between optical and mid-IR diagnostics. Misclassifications of galaxies based on the SDSS spectra are rare despite the presence of dust obscuration. The luminosity of the [NeII] 12.8 \mu m emission-line is well correlated with the star formation rate (SFR) measured from the SDSS spectra, and this holds for the star forming, composite, and AGN-dominated systems. AGN show a clear excess of [NeIII] 15.6 \mu m emission relative to star forming and composite systems. We find good qualitative agreement between various parameters that probe the relative contributions of the AGN and star formation, including: the mid-IR spectral slope, the ratio of the [NeV] 14.3 \mu m to [NeII] \mu m 12.8 fluxes, the equivalent widths of the 7.7, 11.3, and 17 $\mu m$ PAH features, and the optical "D" parameter which measures the distance a source lies from the locus of star forming galaxies in the optical BPT emission-line diagnostic diagram. We also consider the behavior of the three individual PAH features by examining how their flux ratios depend upon the degree of AGN-dominance. We find that the PAH 11.3 \mu m feature is significantly suppressed in the most AGN-dominated systems.

The Spitzer-SDSS-GALEX Spectroscopic Survey (SSGSS) provides a new sample of 101 star-forming galaxies at z < 0.2 with unprecedented multi-wavelength coverage. New mid- to far-infrared spectroscopy from the Spitzer Space Telescope is added to a rich suite of previous imaging and spectroscopy, including ROSAT, Galaxy Evolution Explorer, Sloan Digital Sky Survey, Two Micron All Sky Survey, and Spitzer/SWIRE. Sample selection ensures an even coverage of the full range of normal galaxy properties, spanning two orders of magnitude in stellar mass, color, and dust attenuation. In this paper we present the SSGSS data set, describe the science drivers, and detail the sample selection, observations, data reduction, and quality assessment. Also in this paper, we compare the shape of the thermal continuum and the degree of silicate absorption of these typical, star-forming galaxies to those of starburst galaxies. We investigate the link between star formation rate, infrared luminosity, and total polycyclic aromatic hydrocarbon luminosity, with a view to calibrating the latter for spectral energy distribution models in photometric samples and at high redshift. Last, we take advantage of the 5-40 micron spectroscopic and far-infrared photometric coverage of this sample to perform detailed fitting of the Draine et al. dust models, and investigate the link between dust mass and star formation history and active galactic nucleus properties.

The detailed workings of the central engines of powerful quasars remain a mystery. This is primarily due to the fact that, at their cosmological distances, the inner regions of these quasars are spatially unresolvable. Reverberation mapping is now beginning to unlock the physics of the Broad Emission Line Region (BELR) in nearby, low-luminosity quasars, however it is still unknown whether this gas is dominated by virial motion, by outflows, or infall. The challenge is greater for more distant, powerful sources due to the very long response time of the BELR to changes in the continuum. We present a new technique for probing the kinematic properties of the BELR and accretion disk of high-z quasars using differential microlensing, and show how substantial information can be gained through a single observation of a strongly-lensed quasar using integral field spectroscopy. We apply this technique to GMOS IFU observations of the multiply-imaged quasar Q2237+0305, and find that the observed microlensing signature in the CIII] broad emission line favours gravitationally-dominated dynamics over an accelerating outflow.

We investigate the use of MIR PAH bands, continuum and emission lines as probes of star-formation and AGN activity in a sample of 100 'normal' and local (z~0.1) emission-line galaxies. The MIR spectra were obtained with the Spitzer Space Telescope Infrared Spectrograph (IRS) as part of the Spitzer-SDSS-GALEX Spectroscopic Survey (SSGSS) which includes multi-wavelength photometry from the UV to the FIR and optical spectroscopy. The continuum and features were extracted using PAHFIT (Smith et al. 2007), a decomposition code which we find to yield PAH equivalent widths up to ~30 times larger than the commonly used spline methods. Despite the lack of extreme objects in our sample (such as strong AGNs, low metallicity galaxies or ULIRGs), we find significant variations in PAH, continuum and emission line properties and systematic trends between these MIR properties and optically derived physical properties such as age, metallicity and radiation field hardness. We revisit the diagnostic diagram relating PAH equivalent widths and [Ne II]12.8micrometers/[O IV]25.9micrometers line ratios and find it to be in much better agreement with the standard optical star-formation/AGN classification than when spline decompositions are used, while also potentially revealing obscured AGNs. The luminosity of individual PAH components, of the continuum, and with poorer statistics, of the neon emission lines and molecular hydrogen lines, are found to be tightly correlated to the total IR luminosity, making individual MIR components good gauges of the total dust emission in SF galaxies. Like the total IR luminosity, these individual components can be used to estimate dust attenuation in the UV and in Halpha lines based on energy balance arguments. We also propose average scaling relations between these components and dust corrected, Halpha derived star-formation rates.

We investigate the use of mid-infrared PAH bands, continuum and emission lines as probes of star-formation and AGN activity in a sample of 100 `normal' and local (z~0.1) galaxies. The MIR spectra were obtained with the Spitzer IRS as part of the Spitzer-SDSS-GALEX Spectroscopic Survey (SSGSS) which includes multi-wavelength photometry from the UV to the FIR and optical spectroscopy. The spectra were decomposed using PAHFIT (Smith et al. 2007), which we find to yield PAH equivalent widths (EW) up to ~30 times larger than the commonly used spline methods. Based on correlations between PAH, continuum and emission line properties and optically derived physical properties (gas phase metallicity, radiation field hardness), we revisit the diagnostic diagram relating PAH EWs and [NeII]/[OIV] and find it more efficient as distinguishing weak AGNs from star-forming galaxies than when spline decompositions are used. The luminosity of individual MIR component (PAH, continuum, Ne and molecular hydrogen lines) are found to be tightly correlated to the total IR luminosity and can be used to estimate dust attenuation in the UV and in Ha lines based on energy balance arguments.

We present the analysis of the Polycyclic Aromatic Hydrocarbon (PAH) spectra of a sample of 92 typical star forming galaxies at 0.03 < z < 0.2 observed with the Spitzer IRS. We compare the relative strengths of PAH emission features with SDSS optical diagnostics to probe the relationship between PAH grain properties and star formation and AGN activity. Short-to-long wavelength PAH ratios, and in particular the 7.7-to-11.3 micron feature ratio, are strongly correlated with the star formation diagnostics D_n(4000) and H-alpha equivalent width, increasing with younger stellar populations. This ratio also shows a significant difference between active and non-active galaxies, with the active galaxies exhibiting weaker 7.7 micron emission. A hard radiation field as measured by [OIII]/H-beta and [NeIII]_15.6/[NeII]_12.8 affects PAH ratios differently depending on whether this field results from starburst activity or an AGN. Our results are consistent with a picture in which larger PAH molecules grow more efficiently in richer media and in which smaller PAH molecules are preferentially destroyed by AGN.

We present data from Integral Field Spectroscopy for 3 supercompact UV-Luminous Galaxies (ScUVLGs). As nearby (z~0.2), compact (R_50~1-2 kpc), bright Paschen-alpha sources, with unusually high star formation rates (SFR=3-100 M_sun/yr), ScUVLGs are an ideal population for studying detailed kinematics and dynamics in actively star-forming galaxies. In addition, ScUVLGs appear to be excellent analogs to high redshift Lyman Break Galaxies (LBGs) and our results may offer additional insight into the dynamics of LBGs. Previous work by our team has shown that the morphologies of these galaxies exhibit tidal features and companions, and in this study we find that the dynamics of ScUVLGs are dominated by disturbed kinematics of the emission line gas-- suggestive that these galaxies have undergone recent feedback, interactions or mergers. While 2 of the 3 galaxies do display rotation, v/sigma < 1 -- suggesting dispersion dominated kinematics rather than smooth rotation. We also simulate how these observations would appear at z~2. Lower resolution and loss of low surface brightness features causes some apparent discrepancies between the low-z (observed) and high-z (simulated) interpretations and quantitatively gives different values for v/sigma, yet simulations of these low-z analogs manage to detect the brightest regions well and resemble actual high-z observations of LBGs.

We present spatially resolved spectroscopic images of the gravitationally lensed QSO 2237+0305 taken with the GMOS Integral Field Unit (IFU) on the Gemini North telescope. These observations have the best spatial resolution of any IFU observations of this object to date and include the redshifted CIII] and MgII QSO broad lines. Unlike Mediavilla et al. 1998, we find no evidence for an arc of resolved broad line emission in either the CIII] or MgII lines. We calculate the image flux ratios of both the integrated emission lines and the surrounding continua. The flux ratios of the CIII] and MgII emission lines are consistent with each other but differ substantially from their corresponding continuum flux ratios and the radio/mid-IR flux ratios previously published. We argue that the broad emission line region must be microlensed and the CIII] and MgII emission regions must be approximately the same size and co-located along the line-of-sight. Assuming a simple model for the broad emission line region and the continuum region, we show the size of the CIII] / MgII broad line region is ~0.06 h_{70}^{1/2} pc and the continuum region is <= 0.02 h_{70}^{1/2} pc.

We investigate the luminosity evolution of the host galaxies of radio-loud AGN through Hubble Space Telescope imaging of 72 BL Lac objects, including new STIS imaging of nine z > 0.6 BL Lacs. With their intrinsically low accretion rates and their strongly beamed jets, BL Lacs provide a unique opportunity to probe host galaxy evolution independent of the biases and ambiguities implicit in quasar studies. We find that the host galaxies of BL Lacs evolve strongly, consistent with passive evolution from a period of active star formation in the range 0.5 <~ z <~ 2.5, and inconsistent with either passive evolution from a high formation redshift or a non-evolving population. This evolution is broadly consistent with that observed in the hosts of other radio-loud AGN, and inconsistent with the flatter luminosity evolution of quiescent early types and radio-quiet hosts. This indicates that active star formation, and hence galaxy interactions, are associated with the formation for radio-loud AGN, and that these host galaxies preferentially accrete less material after their formation epoch than galaxies without powerful radio jets. We discuss possible explanations for the link between merger history and the incidence of a radio jet.

H1517+656 is an unusual source, even for a BL Lac object. It is one of the most luminous BL Lacs known, with extreme emission properties at radio, optical, and X-ray frequencies. Furthermore, in our WFPC2 snapshot survey we discovered a series of peculiar arcs describing a 2.4 arcsec radius ring surrounding the source. This paper describes follow-up observations with additional WFPC2 bands and the STIS longpass filter, which have revealed this structure to be the remnants of a very recent galaxy merger. Population synthesis modelling has shown that regions of the arcs have stellar populations with age < 20 Myrs. Additionally, the circularity of the arcs indicates that the plane of the collision and hence accretion is very close to the plane of the sky. Given that BL Lac jets are thought to be aligned with the line of sight, this observation may provide a direct link between the transfer of angular momentum in an interaction and the generation of a radio source.

We have studied the host galaxies of a sample of radio-loud AGN spanning more than four decades in the energy output of the nucleus. The core sample includes 40 low-power sources (BL Lac objects) and 22 high-power sources (radio-loud quasars) spanning the redshift range z~0.15 to z~0.5, all imaged with the high spatial resolution of HST. All of the sources are found to lie in luminous elliptical galaxies, which follow the Kormendy relation for normal ellipticals. A very shallow trend is detected between nuclear brightness (corrected for beaming) and host galaxy luminosity. Black hole masses are estimated for the entire sample, using both the bulge luminosity--black hole mass and the velocity dispersion--black hole mass relations for local galaxies. The latter involves a new method, using the host galaxy morphological parameters, mu_e and r_e, to infer the velocity dispersion, sigma, via the fundamental plane correlation. Both methods indicate that the entire sample of radio-loud AGN are powered by very massive central black holes, with M_{black hole} ~ 10^8 to 10^10 M_{sun}$. Eddington ratios range from L/L_{Eddington} ~ 2 x 10^-4 to ~1, with the high- power sources having higher Eddington ratios than the low-power sources. Overall, radio-loud AGN appear to span a very large range in accretion efficiency, which is all but independent of the mass of the host galaxy.

We compare the host galaxies properties of BL Lac objects with those of radio loud quasars (RLQs) imaged by the WFPC2 on board of HST. The considered objects (z<0.5) are always well resolved and their host galaxies satisfactorily modelled by ellipticals. After homogeneous treatment of the data we find RLQs hosts are systematically more luminous (by ~0.7 mag) with respect to the hosts of BL Lacs. Using the M_BH - L_bulge relation, derived for nearby elliptical galaxies, we have evaluated the central black hole masses of our sample of active galaxies. These data are discussed in conjunction with the nuclear luminosity and the Eddington ratio.

The HST NICMOS Camera 2 was used for H-band imaging of 12 BL Lacertae objects taken from the larger sample observed with the WFPC2 in the R band (Urry et al. 2000; Scarpa et al. 2000). Ten of the 12 BL Lacs are clearly resolved, and the detected host galaxies are large, bright ellipticals with average H-band absolute magnitude M=-26.2+-0.45 mag and effective radius 10+-5 kpc. The rest-frame integrated color of the host galaxies is on average R-H=2.3+-0.3, consistent with the value for both radio galaxies and normal, non-active elliptical galaxies, and indicating the dominant stellar population is old. The host galaxies tend to be bluer in their outer regions than in their cores, with average color gradient Delta(R-H)/Delta(log r)=-0.2 mag, again consistent with results for normal non-active elliptical galaxies. The infrared Kormendy relation, derived for the first time for BL Lac host galaxies, is m(e) = 3.8*log(R)+14.8 (where m(e) is the surface brightness at the effective radius R), fully in agreement with the relation for normal ellipticals. The close similarity between BL Lac host galaxies and normal ellipticals suggests the active nucleus has surprisingly little effect on the host galaxy. This supports a picture in which all elliptical galaxies harbor black holes which can be actively accreting for some fraction of their lifetime.

We have used the HST WFPC2 camera to survey 132 BL Lac objects comprising seven complete radio-, X-ray-, and optically-selected samples. We obtained useful images for 110 targets spanning the redshift range 0 < z < 1.3. In two thirds of the BL Lac images, host galaxies are detected, including nearly all for z < 0.5 (58 of 63). The highest redshift host galaxy detected is in a BL Lac object at z=0.664. In 58 of the 72 resolved host galaxies, a de Vaucouleurs profile is significantly preferred, at >99% confidence, over a pure exponential disk; the two fits are comparable in the remaining 14 cases. These results limit the number of disk systems to at most 8% of BL Lacs (at 99% confidence), and are consistent with all BL~Lac host galaxies being ellipticals. The detected host galaxies are luminous ellipticals with a median absolute K-corrected magnitude of M_R= -23.7 +- 0.6 mag, at least one magnitude brighter than M* and comparable to brightest cluster galaxies. The galaxy morphologies are generally smooth and undisturbed, with small or negligible ellipticities (<0.2). There is no correlation between host galaxy and observed nuclear magnitude or estimated jet power corrected for beaming. If black hole mass is correlated linearly with bulge mass in general, this implies a large range in Eddington ratio. Present data strongly support the unification picture with FR I galaxies constituting the bulk of the parent population of BL Lac objects.

We analyze images of BL Lacertae objects obtained with the {\it Hubble Space Telescope} WFPC2 and the F814W filter. The nine objects cover a redshift range of 0.19 to 0.997. The relatively deep images are sufficient to detect galaxies at least one magnitude below M$^{*}_{I}$ (--21.4) and in most cases to three magnitudes below M$^{*}$. Galaxy enhancement over the average background is found around four out of the nine objects. Results for some cases are confirmed by ground-based imaging. In the other cases, the redshifts of the target BL Lac objects may be incorrect or they are truly isolated. These findings reinforce the idea that on average, BL Lac objects are found in regions of above average galaxy density. However, isolated objects apparently can host BL Lac nuclei too, a result that has implications for the processes that trigger/fuel the nuclear activity.

We present HST observations of seven unusual objects from the HST ``snapshot survey'' of BL Lac objects, of which four are gravitational lens candidates. In three cases a double point sources is observed: 0033+595, with 1.58 arcsec separation, and 0502+675 and 1440+122, each with $\sim 0.3$ arcsec separation. The last two also show one or more galaxies, which could be either host or lensing galaxies. If any are confirmed as lenses, these BL Lac objects are excellent candidates for measuring H$_0$ via gravitational time delay because of their characteristic rapid, high amplitude variability. An additional advantage is that, like other blazars, they are likely superluminal radio sources, in which case the source plane is mapped out over a period of years, providing strong additional constraints on the lensing mass distribution. The fourth gravitational lens candidate is 1517+656, which is surrounded by three arclets forming an almost perfect ring of radius 2.4 arcsec. If this is indeed an Einstein ring, it is most likely a background source gravitationally lensed by the BL Lac object host galaxy and possibly a surrounding group or cluster. In the extreme case that all four candidates are true lenses, the derived frequency of gravitational lensing in this BL Lac sample would be an order of magnitude higher than in comparable quasar samples. We also report on three other remarkable BL Lac objects: 0138-097, which is surrounded by a large number of close companion galaxies; 0806+524, whose host galaxy contains an uncommon arc-like structure; and 1959+650, which is hosted by a gas rich elliptical galaxy with a prominent dust lane of $\sim 5\times 10^5 M_\odot$.