We present a large compilation of reddening estimates from broad-line Balmer decrements for active galactic nuclei (AGNs) with measured X-ray column densities. The median reddening is E(B-V) ~ 0.77 +/- 0.10 for type-1 to type-1.9 AGNs with reported column densities. This is notably higher than the median reddening of AGNs from the SDSS. We attribute this to the selection bias of the SDSS towards blue AGNs. For other AGNs we find evidence of a publication bias against reporting low column densities. We find a significant correlation between N_H and E(B-V) but with a large scatter of +/- 1 dex (without accounting for uncertainties in measurements). On average the X-ray columns are consistent with those predicted from E(B-V) for a solar neighbourhood dust-to-gas ratio. We argue that the large scatter of column densities and reddenings can be explained by X-ray column-density variability. For AGNs with detectable broad-line regions (BLRs) that have undergone significant changes of Seyfert type (changing-look AGNs) we do not find any statistically significant differences in N_H or E(B-V) compared to standard type-1 to type-1.9 AGNs. There is no evidence for any type-1 AGNs being Compton-thick. We also analyze type-2 AGNs and find no significant correlation between N_H and narrow-line region reddening. The median column density of LINERs is 22.68 +/- 0.75 compared to a column density of 22.90 +/- 0.28 for type-2s. We find the majority of low column density type-2s are LINERs, but N_H is probably underestimated because of scattered X-ray light diluting the photoelectric edge.

In the standard AGN reverberation-mapping model variations in broad-line region (BLR) fluxes are predicted from optical continuum variability (taken as a proxy for variations in the ionizing continuum) convolved with an impulse function that depends on the geometry of the emitting and reprocessing regions. However, it has long been known that BLR variability can deviate from these predictions. We analyze extensive long-term H$\beta$ and continuum monitoring of NGC 5548 and a large sample of high-quality H$\beta$ light curves of other AGNs to investigate the frequency and characteristics of anomalous responses of the BLR. Anomalies are very common and probably occur in every object. Onsets can be on a timescale only slightly longer than the light-crossing time and durations are of the order of the characteristic timescale of variability of the optical continuum to several times longer. Anomalies are larger when NGC 5548 is in a low state, but otherwise there is no correlation with continuum variability. There is abundant evidence for the optical continuum of AGNs varying independently of the higher-energy continua and this is sufficient to explain the anomalous responses. There are good reasons for believing that the frequent lack of correlation between different spectral regions is due to anisotropic and non-axisymmetric emission. Rapid changes in line profiles and velocity-dependent lags are consistent with this. Motion of compact absorbing clouds across the line of sight is another possible cause of anomalies. The prevalence of anomalies needs to be taken into account when planning reverberation-mapping campaigns.

Two major challenges to unification schemes for active galactic nuclei (AGN) are the existence of Narrow-Line Seyfert 1s (NLS1s) and the existence of changing-look (CL) AGNs. AGNs can drastically change their spectral appearance in the optical (changing their Seyfert type) and/or in the X-ray region. We illustrate the CL phenomenon with our multi-wavelength monitoring of NGC 2617 and discuss its properties compared with NLS1s. There are few examples of CL NLS1s and the changes are mostly only in the X-ray region. It has been proposed that some of these could be cases of a tidal-disruption events (TDE) or supernova events. If BLRs have a flat geometry and NLS1s are seen face-on then we have to see CL cases only if the orientation of the BLR changes as a result of a TDE or a close encounter of a star without a TDE. If NLS1s include both high Eddington accretion rate and low-inclination AGNs then a significant fraction of NLS1s could be obscured and would not be identified as NLS1s. CL cases might happen more in such objects if dust sublimation occurs following a strong increase in the optical luminosity.

We present results of a study of the correlation between the infrared (JHKL) and optical(B) fluxes of the nucleus of the Seyfert galaxy NGC 4151 for the years 2010 - 2015 using our own data (partially published) in combination with published data of Roberts and Rumstey (2012), Guo et al. (2014) and Schnulle et al. (2013, 2015).We find similar lags for each of the HKL passbands relative to the optical of 37 +-3 days. The lags are the same to within the accuracy of measurement. We do not confirm a significant decrease in the lag for HKL in 2013-2014 previously reported by Schnulle et al. (2015), but we find that the lag of the short-lag component of J increased. We discuss our results within the framework of the standard model, where the variable infrared radiation is mainly due to the thermal re-emission of short-wave radiation by dust clouds close to a variable central source. There is also some contribution to the IR emission from the accretion disk, and this contribution increases with decreasing wavelength. The variability in J and K is not entirely simultaneous, which may be due to the differing contributions of the radiation from the accretion disk in these bands. The absence of strong wavelength-dependent changes in infrared lag across the HKL passbands can be explained by having the dust clouds during 2010-2015 be located beyond the sublimation radius. The relative wavelength independence of the infrared lags is also consistent with the hollow bi-conical outflow model of Oknyansky et al. (2015).

We report the discovery of non-stellar hydrogen Balmer and metastable helium absorption lines accompanying a transient, high-velocity (0.05$c$) broad absorption line (BAL) system in the optical spectra of the tidal disruption event (TDE) AT2018zr ($z=0.071$). In the HST UV spectra, absorption of high- and low-ionization lines are also present at this velocity, making AT2018zr resemble a low-ionization broad absorption line (LoBAL) QSO. We conclude that these transient absorption features are more likely to arise in fast outflows produced by the TDE than absorbed by the unbound debris. In accordance with the outflow picture, we are able to reproduce the flat-topped H$\alpha$ emission in a spherically expanding medium, without invoking the typical prescription of an elliptical disk. We also report the appearance of narrow ($\sim$1000~km~s$^{-1}$) NIII$\lambda$4640, HeII$\lambda4686$, H$\alpha$, and H$\beta$, emission in the late-time optical spectra of AT2018zr, which may be a result of UV continuum hardening at late time as observed by Swift. Including AT2018zr, we find a high association rate (3 out of 4) of BALs in the UV spectra of TDEs. This suggests that outflows may be ubiquitous among TDEs and may be less sensitive to viewing angle effects compared to QSO outflows.

The profiles of the broad emission lines of active galactic nuclei (AGNs) and the time delays in their response to changes in the ionizing continuum ("lags") give information about the structure and kinematics of the inner regions of AGNs. Line profiles are also our main way of estimating the masses of the supermassive black holes (SMBHs). However, the profiles often show ill-understood, asymmetric structure and velocity-dependent lags vary with time. Here we show that partial obscuration of the broad-line region (BLR) by outflowing, compact, dusty clumps produces asymmetries and velocity-dependent lags similar to those observed. Our model explains previously inexplicable changes in the ratios of the hydrogen lines with time and velocity, the lack of correlation of changes in line profiles with variability of the central engine, the velocity dependence of lags, and the change of lags with time. We propose that changes on timescales longer than the light-crossing time do not come from dynamical changes in the BLR, but are a natural result of the effect of outflowing dusty clumps driven by radiation pressure acting on the dust. The motion of these clumps offers an explanation of long-term changes in polarization. The effects of the dust complicate the study of the structure and kinematics of the BLR and the search for sub-parsec SMBH binaries. Partial obscuration of the accretion disc can also provide the local fluctuations in luminosity that can explain sizes deduced from microlensing.

We present a Hubble Space Telescope (HST) survey of extended [O III] {\lambda}5007 emission for a sample of 12 nearby (z < 0.12), luminous Type 2 quasars (QSO2s), which we use to measure the extent and kinematics of their AGN-ionized gas. We find the size of the observed [O III] regions scale with luminosity in comparison to nearby, less luminous Seyfert galaxies and radially outflowing kinematics to exist in all targets. We report an average maximum outflow radius of $\sim$600 pc, with gas continuing to be kinematically influenced by the central AGN out to an average radius of $\sim$1130 pc. These findings question the effectiveness of AGN being capable of clearing material from their host bulge in the nearby universe and suggest that disruption of gas by AGN activity may prevent star formation without requiring evacuation. Additionally, we find a dichotomy in our targets when comparing [O III] radial extent and nuclear FWHM, where QSO2s with compact [O III] morphologies typically possess broader nuclear emission-lines.

Optical and near-infrared photometry, optical spectroscopy, and soft X-ray and UV monitoring of the changing-look active galactic nucleus NGC 2617 show that it continues to have the appearance of a type-1 Seyfert galaxy. An optical light curve for 2010$-$2017 indicates that the change of type probably occurred between 2010 October and 2012 February and was not related to the brightening in 2013. In 2016 and 2017 NGC 2617 brightened again to a level of activity close to that in 2013 April. However, in 2017 from the end of the March to end of July 2017 it was in very low level and starting to change back to a Seyfert 1.8. We find variations in all passbands and in both the intensities and profiles of the broad Balmer lines. A new displaced emission peak has appeared in H$\beta$. X-ray variations are well correlated with UV$-$optical variability and possibly lead by $\sim$ 2$-$3 d. The $K$ band lags the $J$ band by about 21.5 $\pm$ 2.5 d and lags the combined $B + J$ bands by $\sim$ 25 d. $J$ lags $B$ by $\sim$ 3 d. This could be because $J$-band variability arises predominantly from the outer part of the accretion disc, while $K$-band variability is dominated by thermal re-emission by dust. We propose that spectral-type changes are a result of increasing central luminosity causing sublimation of the innermost dust in the hollow bi-conical outflow. We briefly discuss various other possible reasons that might explain the dramatic changes in NGC 2617.

Optical and near-infrared photometry, optical spectroscopy, and soft X-ray and UV monitoring of the changing look active galactic nucleus NGC 2617 show that it continues to have the appearance of a type-1 Seyfert galaxy. An optical light curve for 2010-2016 indicates that the change of type probably occurred between 2010 October and 2012 February and was not related to the brightening in 2013. In 2016 NGC 2617 brightened again to a level of activity close to that in 2013 April. We find variations in all passbands and in both the intensities and profiles of the broad Balmer lines. A new displaced emission peak has appeared in H$\beta$. X-ray variations are well correlated with UV-optical variability and possibly lead by $\sim$ 2-3 d. The $K$ band lags the $J$ band by about 21.5 $\pm$ 2.5 d. and lags the combined $B+J$ filters by $\sim$ 25 d. $J$ lags $B$ by about 3 d. This could be because $J$-band variability arises from the outer part of the accretion disc, while $K$-band variability comes from thermal re-emission by dust. We propose that spectral-type changes are a result of increasing central luminosity causing sublimation of the innermost dust in the hollow biconical outflow. We briefly discuss various other possible reasons that might explain the dramatic changes in NGC 2617.

Low-redshift active galactic nuclei (AGNs) with extremely blue optical spectral indices are shown to have a mean, velocity-averaged, broad-line H$\alpha$/H$\beta$ ratio of $\thickapprox 2.72 \pm 0.04$, consistent with a Baker-Menzel Case B value. Comparison of a wide range of properties of the very bluest AGNs with those of a luminosity-matched subset of the Dong et al. blue AGN sample indicates that the only difference is the internal reddening. Ultraviolet fluxes are brighter for the bluest AGNs by an amount consistent with the flat AGN reddening curve of Gaskell et al.(2004). The lack of a significant difference in the GALEX (FUV--NUV) colour index strongly rules out a steep SMC-like reddening curve and also argues against an intrinsically harder spectrum for the bluest AGNs. For very blue AGNs the Ly$\alpha$/H$\beta$ ratio is also consistent with being the Case B value. The Case B ratios provide strong support for the self-shielded broad-line model of Gaskell, Klimek & Nazarova. It is proposed that the greatly enhanced Ly$\alpha$/H$\beta$ ratio at very high velocities is a consequence of continuum fluorescence in the Lyman lines (Case C). Reddenings of AGNs mean that the far-UV luminosity is often underestimated by up to an order of magnitude. This is a major factor causing the discrepancies between measured accretion disc sizes and the predictions of simple accretion disc theory. Dust covering fractions for most AGNs are lower than has been estimated. The total mass in lower mass supermassive black holes must be greater than hitherto estimated.

We use mid-IR to UV observations to derive a mean attenuation curve out to the rest-frame extreme ultraviolet (EUV) for "BAL dust" -- the dust causing the additional extinction of active galactic nuclei (AGNs) with broad absorption lines (BALQSOs). In contrast to the normal, relatively flat, mean AGN attenuation curve, BAL dust is well fit by a steeply rising, SMC-like curve. We confirm the shape of the theoretical Weingartner & Draine SMC curve out to 700 \AA, but the drop in attenuation at still shorter wavelengths is less than predicted. The identical attenuation curve for low-ionization BALQSOs (LoBALs) does not support them being a "break out" phase in the life of AGNs. Although attenuation in the optical due to BAL dust is low ($E(B-V) \sim 0.03 - 0.05$), the attenuation rises to one magnitude in the EUV because of the steep extinction curve. Here the dust optical depth is at the optimum value for radiative acceleration of dusty gas. Because the spectral energy distribution of AGNs peaks in the EUV where the optical depth is highest, the force on the dust dominates the acceleration of BAL gas. For LoBALs we get a negative attenuation curve in the optical. This is naturally explained if there is more light scattered into our line of sight in LoBALs compared with non-BALQSOs. We suggest that this and partial covering are causes when attenuation curves appear to be steeper in the UV that an SMC curve.

We present an analysis of 43 years (1972 to 2015) of spectroscopic observations of the Seyfert 1 galaxy NGC 5548. This includes 12 years of new unpublished observations (2003 to 2015). We compiled about 1600 H$\beta$ spectra and analyzed the long-term spectral variations of the 5100 \AA\ continuum and the H$\beta$ line. Our analysis is based on standard procedures including the Lomb-Scargle method, which is known to be rather limited to such heterogeneous data sets, and new method developed specifically for this project that is more robust and reveals a $\sim$5700 day periodicity in the continuum light curve, the H$\beta$ light curve and the radial velocity curve of the red wing of the H$\beta$ line. The data are consistent with orbital motion inside the broad emission line region of the source. We discuss several possible mechanisms that can explain this periodicity, including orbiting dusty and dust-free clouds, a binary black hole system, tidal disruption events, and the effect of an orbiting star periodically passing through an accretion disc.

We use the Balmer decrements of the broad-line regions (BLRs) and narrow-line regions (NLRs) of active galactic nuclei (AGNs) as reddening indicators to investigate the location of the dust for four samples of AGNs with reliable estimates of the NLR contribution to the Balmer lines. Intercomparison of the NLR and BLR Balmer decrements indicates that the reddening of the NLR sets a lower limit to the reddening of the BLR. Almost no objects have high NLR reddening but low BLR reddening. The reddening of the BLR is often substantially greater than the reddening of the NLR. The BLR reddening is correlated with the equivalent widths of [O III] lines and the intensity of the [O III] lines relative to broad H\beta. We find these relationships to be consistent with the predictions of a simple model where the additional dust reddening the BLR is interior to the NLR. We thus conclude that the dust causing the additional reddening of the accretion disc and BLR is mostly located at a smaller radius than the NLR.

We show that, contrary to simple predictions, most AGNs show at best only a small increase of lags in the J, H, K, and L bands with increasing wavelength. We suggest that a possible cause of this near simultaneity of the variability from the near-IR to the mid-IR is that the hot dust is in a hollow bi-conical outflow of which we only see the near side. Although most AGNs show near simultaneity of IR variability, there was at least one epoch when NGC 4151 showed the sharply increasing IR lag with the increase of the wavelength. This behaviour might also be present in GQ Comae. We discuss these results briefly. The relative wavelength independence of IR lags simplifies the use of IR lags for estimating cosmological parameters.

The high-ionization lines of the broad-line region (BLR) of thermal active galactic nuclei (AGNs) show blueshifts of a few hundred km/s to several thousand km/sec with respect to the low-ionization lines. This has long been thought to be due to the high-ionization lines of the BLR arising in a wind of which the far side of the outflow is blocked from our view by the accretion disc. Evidence for and against the disc-wind model is discussed. The biggest problem for the model is that velocity-resolved reverberation mapping repeatedly fails to show the expected kinematic signature of outflow of the BLR. The disc-wind model also cannot readily reproduce the red side of the line profiles of high-ionization lines. The rapidly falling density in an outflow makes it difficult to obtain high equivalent widths. We point out a number of major problems with associating the BLR with the outflows producing broad absorption lines. An explanation which avoids all these problems and satisfies the constraints of both the line profiles and velocity-resolved reverberation-mapping is a model in which the blueshifting is due to scattering off material spiraling inwards with an inflow velocity of half the velocity of the blueshifting. We discuss how recent reverberation mapping results are consistent with the scattering-plus-inflow model but do not support a disc-wind model. We propose that the anti-correlation of the apparent redshifting of H$\beta$ with the blueshifting of CIV is a consequence of contamination of the red wings of H$\beta$ by the broad wings of [O III].

A growing body of evidence suggests that part of, if not all, scattering regions of active galactic nuclei (AGNs) are clumpy. Hence. in this paper, we run radiative transfer models in the optical/UV for a variety of AGN reprocessing regions with different distributions of clumpy scattering media. We use the latest version of the Monte Carlo code STOKES presented in the first two papers of this series to model AGN reprocessing regions of increasing morphological complexity. We replace previously uniform-density media with up to thousands of constant-density clumps. We couple a continuum source to fragmented equatorial scattering regions, polar outflows, and toroidal, obscuring dust regions and investigate a wide range of geometries. We also consider different levels of fragmentation in each scattering region to evaluate importance of fragmentation for the net polarization of the AGN. We find that, in comparison with uniform-density models, equatorial distributions of gas and dust clouds result in grayer spectra, and show a decrease of the net polarization percentage at all lines of sight. The resulting polarization position angle depends on the morphology of the clumpy structure, with extended tori favoring parallel polarization while compact tori produce orthogonal polarization position angles. In the case of polar scattering regions, fragmentation increases the net polarization unless the cloud filling factor is small. A complete set of AGN models constructed from the individual, fragmented regions is investigated. Our modeling shows that the introduction of fragmented dusty tori significantly alters the resulting net polarization of an AGN. Comparison of our models to polarization observations of large AGN samples greatly favors geometrically compact clumpy tori over extended ones.

We investigate the correlation between infrared (JHKL) and optical (B) fluxes of the variable nucleus of the Seyfert galaxy NGC 4151 using partially published data for the last 6 years (2008-2013.). Here we are using the same data as in Oknyansky et al. (2014), but include also optical (B) data from Guo et al. We find that the lag of flux in all the infrared bands is the same, 40 +- 6 days, to within the measurement accuracy. Variability in the J and K bands is not quite simultaneous, perhaps due to the differing contributions of the accretion disk in these bands. The lag found for the K band compared with the B band is not significantly different from earlier values obtained for the period 2000-2007. However, finding approximately the same lags in all IR bands for 2008-2013 differs from previous results at earlier epochs when the lag increased with increasing wavelength. Examples of almost the same lag in different IR bands are known for some other active nuclei. In the case of NGC 4151 it appears that the relative lags between the IR bands may be different in different years. The available data, unfortunately, do not allow us to investigate a possible change in the lags during the test interval. We discuss our results in the framework of the standard model where the variable infrared radiation is mainly due to thermal re-emission from the part of the dusty torus closest to the central source. There is also a contribution of some IR emission from the accretion disk, and this contribution increases with decreasing wavelength. Some cosmological applications of obtained results are discussed.

The STOKES Monte Carlo radiative transfer code has been extended to model the velocity dependence of the polarization of emission lines. We use STOKES to present improved modelling of the velocity-dependent polarization of broad emission lines in active galactic nuclei. We confirm that off-axis continuum emission can produce observed velocity dependencies of both the degree and position angle of polarization. The characteristic features are a dip in the percentage polarization and an S-shaped swing in the position angle of the polarization across the line profile. Some differences between our STOKES results and previous modelling of polarization due to off-axis emission are noted. In particular we find that the presence of an offset between the maximum in line flux and the dip in the percentage of polarization or the central velocity of the swing in position angle does not necessarily imply that the scattering material is moving radially. Our model is an alternative scenario to the equatorial scattering disk described by Smith et al. (2005). We discuss strategies to discriminate between both interpretations and to constrain their relative contributions to the observed velocity-resolved line and polarization.

The tidal disruption event (TDE) PS1-10jh lacked strong Balmer lines but showed strong, broad, He II emission both before maximum light and for at least 8 months thereafter. Gezari et al. interpreted this as evidence for the disruption of a rare hydrogen-deficient star. However, Guillochon et al. have argued instead that the disrupted star was a normal main-sequence star and that the strength of the He II emission compared with the Balmer lines is a result the emission being similar to the broad-line region (BLR) of an AGN, but lacking the outer, lower-ionization BLR gas. We show that the profile of He II 4686 in PS1-10jh is similar to the blueshifted profiles of high-ionization lines in AGNs. We find an He II 4686/Halpha ratio for PS1-10jh of ~ 3.7. We show that both the high-velocity gas of the inner BLR of normal AGNs and the spectra of type II-P supernovae right after shock break out also produce very high He II 4686/Halpha ratios. A high He II 4686/Halpha ratio can thus be produced with a solar H/He abundance ratio. We demonstrate from photoionization modelling that the estimated He II 4686/Halpha ratio can be produced with a BLR truncated before the He++ Stroemgren length if the density is ~ 10^{11} cm^{-3}. The similarity of the He II 4686 emission in PS1-10jh to the emission from the inner BLRs of AGNs supports the idea that the emission after a TDE event is similar to that of normal AGNs.

Velocity-resolved reverberation mapping suggests that the broad-line regions (BLRs) of AGNs can have significant net inflow. We use the STOKES radiative transfer code to show that electron and Rayleigh scattering off the BLR and torus naturally explains the blueshifted profiles of high-ionization lines and the ionization dependence of the blueshifts. This result is insensitive to the geometry of the scattering region. If correct, this model resolves the long-standing conflict between the absence of outflow implied by velocity-resolved reverberation mapping and the need for outflow if the blueshifting is the result of obscuration. The accretion rate implied by the inflow is sufficient to power the AGN. We suggest that the BLR is part of the outer accretion disk and that similar MHD processes are operating. In the scattering model the blueshifting is proportional to the accretion rate so high-accretion-rate AGNs will show greater high-ionization line blueshifts as is observed. Scattering can lead to systematically too high black hole mass estimates from the C IV line. We note many similarities between narrow-line region (NLR) and BLR blueshiftings, and suggest that NLR blueshiftings have a similar explanation. Our model explains the higher blueshifts of broad absorption line QSOs if they are more highly inclined. Rayleigh scattering from the BLR and torus could be more important in the UV than electron scattering for predominantly neutral material around AGNs. The importance of Rayleigh scattering versus electron scattering can be assessed by comparing line profiles at different wavelengths arising from the same emission-line region.

The anti-correlations between the equivalent widths of emission lines and the continuum luminosity in AGNs, known as the Baldwin effect are well established for broad lines, but are less well studied for narrow lines. In this paper we explore the Baldwin effect of narrow emission lines over a wide range of ionization levels and critical densities using a large sample of broad-line, radio-quiet AGNs taken from Sloan Digital Sky Survey (SDSS) Data Release 4. These type1 AGNs span three orders of magnitude in continuum luminosity. We show that most narrow lines show a similar Baldwin effect slope of about -0.2 while the significant deviations of the slopes for [NII] 6583, [OII] 3727, [NeV] 3425, and the narrow component of Ha can be explained by the influence of metallicity, star-formation contamination and possibly by difference in the shape of the UV-optical continuum. The slopes do not show any correlation with either the ionization potential or the critical density. We show that a combination of 50% variations in continuum near 5100A and a log-normal distribution of observed luminosity can naturally reproduce a constant Baldwin effect slope of -0.2 for all narrow lines. The variations of the continuum could be due to variability, intrinsic anisotropic emission, or an inclination effect.

We model the spectropolarimetric signature resulting from the radiative coupling between the innermost parts of active galactic nuclei (AGNs). We use a new public version of STOKES, a Monte Carlo radiative transfer code presented in the first paper of this series. The code has been significantly improved for computational speed and polarization imaging has been implemented. We couple continuum sources with equatorial scattering regions, polar outflows, and toroidal obscuring dust and we study the resulting polarization. We compute a grid of thermal AGN models for different half-opening angles of the torus and polar winds. We also consider a range of optical depths for equatorial and polar electron scattering and investigate how the model geometry influences the type-1/type-2 polarization dichotomy for thermal AGNs. We put new constrains on the inflowing medium within the inner walls of the torus. The inflow should be confined to the common equatorial plane of the torus and the accretion disc and have a radial optical depth of 1 < tau < 3. Our modeling of type-1 AGNs also indicates that the torus is more likely to have a large (~ 60{\deg}) half-opening angle. Polarization perpendicular to the axis of the torus may arise at a type-1 viewing angle for a torus half-opening angle of 30{\deg}- 45{\deg} or polar outflows with an optical depth near unity. Our modeling suggests that most Seyfert-2 AGN must have a half-opening angle > 60{\deg} to match the level of perpendicular polarization expected. If outflows are collimated by the torus inner walls, they must not be optically thick (tau < 1) in order to preserve the polarization dichotomy. The wind's optical depth is found not to play a critical role for the degree of polarization of type-2 thermal AGNs but it has a significant impact on the type-1/type-2 polarization dichotomy when the optical depth exceeds tau = 0.3.

Observations of the optical polarization of NGC 4151 in 1997-2003 show variations of an order of magnitude in the polarized flux while the polarization position angle remains constant. The amplitude of variability of the polarized flux is comparable to the amplitude of variability of the total U-band flux, except that the polarized flux follows the total flux with a lag of 8 +/- 3 days. The time lag and the constancy of the position angle strongly favor a scattering origin for the variable polarization rather than a non-thermal synchrotron origin. The orientation of the position angle of the polarized flux (parallel to the radio axis) and the size of the lag imply that the polarization arises from electron scattering in a flattened region within the low-ionization component of the broad-line-region. Polarization from dust scattering in the equatorial torus is ruled out as the source of the lag in polarized flux because it would produce a larger lag and, unless the half-opening angle of the torus is more than 53 degrees, the polarization would be perpendicular to the radio axis. We note a long-term change in the percentage of polarization at similar total flux levels and this could be due either to changing non-axisymmetry in the optical continuum emission, or a change in the number of scatterers on a timescale of years.

A simple refinement is proposed to the Dibai method for determining black hole masses in type-1 thermal AGNs. Comparisons with reverberation mapping black hole masses and host galaxy bulge properties suggest that the method is accurate to +/- 0.15 dex. Contrary to what was thought when the black hole mass - stellar velocity dispersion ("M - sigma") relationship was first discovered, it does not have a lower dispersion than the black hole mass - bulge luminosity ("M - L") relationship. The dispersion in the M - L relationship for AGNs decreases strongly with increasing black hole mass or bulge luminosity. This is naturally explained as a consequence of the black hole - bulge relationships being the result of averaging due to mergers. Simulations show that the decrease in dispersion in the M - L relationship with increasing mass is in qualitative agreement with being driven by mergers. The large scatter in AGN black hole masses at lower masses rules out significant AGN feedback. A non-causal origin of the correlations between black holes and bulges explains the frequent lack of supermassive black holes in late-type galaxies, and the lack of correlation of black hole mass with pseudo-bulges.

The general picture of how thermal AGNs work has become clearer in recent years but major observational puzzles threaten to undermine this picture. These puzzles include AGNs with extremely asymmetric emission line profiles, inconsistent multi-wavelength variability, rapid apparent changes in the sizes of emitting regions and in the direction of gas flow, a curious insensitivity of gas in some narrow velocity ranges to changes in the ionizing continuum, and differing dependencies of polarization on gas velocity. It is proposed that all these puzzles can readily be explained by off-axis variability, and that there is no need to invoke exotic explanations such as binary supermassive black holes or recoiling black holes.

We use homogeneous samples of radio-quiet Seyfert 1 galaxies and QSOs selected from the Sloan Digital Sky Survey to investigate the connection between the velocity shift and the equivalent width (EW) of the [OIII] 5007 emission line, and their correlations with physical parameters of active galactic nuclei (AGNs). We find a significant and negative correlation between the EW of the core component, EW(core), and the blueshift of either the core (the peak), the wing, or the total profile of [OIII] emission; it is fairly strong for the blueshift of the total profile particularly. However, both quantities (EW and velocity shift) generally have only weak, if any, correlations with fundamental AGN parameters such as the nuclear continuum luminosity at 5100 L_{5100}, black hole mass (M_{BH}), and the Eddington ratio (L/L_{Edd}); these correlations include the classical Baldwin effect of EW(core), an inverse Baldwin effect of EW(wing), and the relationship between velocity shifts and \lratio. Our findings suggest that both the large object-to-object variation in the strength of [OIII] emission and the blueshift--EW(core) connection are not governed primarily by fundamental AGN parameters such as L_{5100}, M_{BH} and L/L_{Edd}. We propose that the ISM conditions of the host galaxies play a major role instead in the diversity of the [OIII] properties in active galaxies. This suggests that the use of[OIII] 5007 luminosity as proxy of AGN luminosity does not depend strongly on the above-mentioned fundamental AGN parameters.

It is shown that broad-line region (BLR) line profiles ranging from the classic "logarithmic" profile to double-peaked, disk-like profiles are readily explained by the distribution of BLR gas proposed by Gaskell, Klimek, & Nazarova (2007; GKN) without any need to invoke fundamental differences in the AGNs other than differing viewing angles. It is argued that the highly-variable thermal energy generation in AGNs originates off axis in regions that cannot be axially symmetric. This off-axis model readily explains the varying degrees of temporal correlation found in multi-wavelength variability studies, the strong, variable asymmetry of BLR line profiles, the varying time delays in the response of the BLR to different continuum events, how narrow velocity ranges of line profiles will often appear to respond differently or not at all to continuum variability, complex changes in the Balmer decrement with velocity, inconsistent and variable inflow/outflow signatures found in velocity-resolved reverberation mapping, the diversity of velocity-dependent polarizations observed, and polarization variability. The fundamentally non-axisymmetric nature of AGN continuum variability severely limits what can be learned from reverberation mapping. In particular, high-fidelity reverberation mapping is not possible. There will be systematic orientation-dependent errors in black hole mass determinations. The effects of off-axis emission will mask subtle signatures of possible close supermassive black hole binaries. Some tests of the off-axis-variability model are proposed.

We present the final results from a high sampling rate, multi-month, spectrophotometric reverberation mapping campaign undertaken to obtain either new or improved Hbeta reverberation lag measurements for several relatively low-luminosity AGNs. We have reliably measured thetime delay between variations in the continuum and Hbeta emission line in six local Seyfert 1 galaxies. These measurements are used to calculate the mass of the supermassive black hole at the center of each of these AGNs. We place our results in context to the most current calibration of the broad-line region (BLR) R-L relationship, where our results remove outliers and reduce the scatter at the low-luminosity end of this relationship. We also present velocity-resolved Hbeta time delay measurements for our complete sample, though the clearest velocity-resolved kinematic signatures have already been published.

There is a strong decrease in scatter in the black hole mass versus bulge luminosity relationship with increasing luminosity and very little scatter for the most luminous galaxies. It is shown that this is a natural consequence of the substantial initial dispersion in the ratio of black hole mass to total stellar mass and of subsequent galaxy growth through hierarchical mergers. "Fine-tuning" through feedback between black hole growth and bulge growth is neither necessary nor desirable.

A new empirical formulae is given for estimating the masses of black holes in AGNs from the H beta velocity dispersion and the continuum luminosity at 5100 Angstroms. It is calibrated to reverberation-mapping and stellar-dynamical estimates of black hole masses. The resulting mass estimates are as accurate as reverberation-mapping and stellar-dynamical estimates. The new mass estimates show that there is very little scatter in the M_{bh} - L_{bulge} relationship for high-luminosity galaxies, and that the scatter increases substantially in lower-mass galaxies.

The evidence from velocity-resolved reverberation mapping showing a net infall of the broad-line region (BLR) of AGNs is reviewed. Different lines in many objects at different epochs give a consistent picture of BLR motions. The motions are dominated by virialized motion (rotation plus turbulence) with significant net inflow. The BLR mass influx is sufficient to power the AGN. The increasing blueshifting of lines with increasing ionization potential is a consequence of scattering off infalling material. The high blueshiftings of the UV lines in Narrow-Line Seyfert 1s are due to enhanced BLR inflow rates rather than strong winds. Seemingly conflicting cases of apparent outflow reverberation mapping signatures are a result of the breakdown of the axial-symmetry assumption in reverberation mapping. There are several plausible causes of this breakdown: high-energy variability tends to be intrinsically anisotropic, regions of variability are necessarily located off-axis, and X-ray observations reveal major changes in line-of-sight column densities close to the black hole. Results from reverberation mapping campaigns dominated by a single event need to be treated with caution.

I review progress made in understanding the nature of the broad-line region (BLR) of active galactic nuclei (AGNs) and the role BLRs play in the AGN phenomenon. The high equivalent widths of the lines imply a high BLR covering factor, and the absence of clear evidence for absorption by the BLR means that the BLR has a flattened distribution and that we always view it near pole-on. The BLR gas is strongly self-shielding near the equatorial plane. Velocity-resolved reverberation mapping has long strongly excluded significant outflow of the BLR and shows instead that the predominant motions are Keplerian with large turbulence and a significant net inflow. The rotation and turbulence are consistent with the inferred geometry. The blueshifting of high-ionization lines is a consequence of scattering off inflowing material rather than the result of an outflowing wind. The rate of inflow of the BLR is sufficient to provide the accretion rate needed to power the AGN. Because the motions of the BLR are gravitationally dominated, and the BLR structure is very similar in most AGNs, consistent black hole masses can be determined. The good correlation between these estimates and masses predicted from the bulge luminosities of host galaxies provides strong support for the similarity of AGN continuum shapes and the correctness of the BLR picture presented. It is concluded that although many mysteries remain about the details of how AGNs work, a general overall picture of the torus and BLR is becoming clear.

It has been proposed that when the peaks of the broad emission lines in active galactic nuclei (AGNs) are significantly blueshifted or redshifted from the systemic velocity of the host galaxy, this could be a consequence of orbital motion of a supermassive blackhole binary (SMB) (Gaskell 1983). The AGN J1536+0441 (=SDSS J153636.22+044127.0) has recently been proposed as an example of this phenomenon (Boroson & Lauer 2009). It is proposed here instead that 1536+044 is an example of line emission from a disc. If this is correct, the lack of clear optical spectral evidence for close SMBs is significant and argues either that the merging of close SMBs is much faster than has generally been hitherto thought, or if the approach is slow, that when the separation of the binary is comparable to the size of the torus and broad-line region, the feeding of the black holes is disrupted.

We present results from a multi-month reverberation mapping campaign undertaken primarily at MDM Observatory with supporting observations from around the world. We measure broad line region (BLR) radii and black hole masses for six objects. A velocity-resolved analysis of the H_beta response shows the presence of diverse kinematic signatures in the BLR.

A detailed analysis of the data from a high sampling rate, multi-month reverberation mapping campaign, undertaken primarily at MDM Observatory with supporting observations from telescopes around the world, reveals that the Hbeta emission region within the broad line regions (BLRs) of several nearby AGNs exhibit a variety of kinematic behaviors. While the primary goal of this campaign was to obtain either new or improved Hbeta reverberation lag measurements for several relatively low luminosity AGNs (presented in a separate work), we were also able to unambiguously reconstruct velocity-resolved reverberation signals from a subset of our targets. Through high cadence spectroscopic monitoring of the optical continuum and broad Hbeta emission line variations observed in the nuclear regions of NGC 3227, NGC 3516, and NGC 5548, we clearly see evidence for outflowing, infalling, and virialized BLR gas motions, respectively.

We present the results of extensive multi-frequency monitoring of the radio galaxy 3C 120 between 2002 and 2007 at X-ray, optical, and radio wave bands, as well as imaging with the Very Long Baseline Array (VLBA). Over the 5 yr of observation, significant dips in the X-ray light curve are followed by ejections of bright superluminal knots in the VLBA images. Consistent with this, the X-ray flux and 37 GHz flux are anti-correlated with X-ray leading the radio variations. This implies that, in this radio galaxy, the radiative state of accretion disk plus corona system, where the X-rays are produced, has a direct effect on the events in the jet, where the radio emission originates. The X-ray power spectral density of 3C 120 shows a break, with steeper slope at shorter timescale and the break timescale is commensurate with the mass of the central black hole based on observations of Seyfert galaxies and black hole X-ray binaries. These findings provide support for the paradigm that black hole X-ray binaries and active galactic nuclei are fundamentally similar systems, with characteristic time and size scales linearly proportional to the mass of the central black hole. The X-ray and optical variations are strongly correlated in 3C 120, which implies that the optical emission in this object arises from the same general region as the X-rays, i.e., in the accretion disk-corona system. We numerically model multi-wavelength light curves of 3C 120 from such a system with the optical-UV emission produced in the disk and the X-rays generated by scattering of thermal photons by hot electrons in the corona. From the comparison of the temporal properties of the model light curves to that of the observed variability, we constrain the physical size of the corona and the distances of the emitting regions from the central BH.

An improved transformation of the full width at half maxima (FWHM) of the [O III] 5007 line in AGNs to the stellar velocity dispersion, sigma, of the host galaxy is given. This significantly reduces the systematic errors in using the [O III] FWHM as a proxy for sigma. AGN black hole masses, M, estimated using the Dibai single epoch spectrum method, are combined with the new estimates of sigma to give a revised AGN M-sigma relationship extending up to high masses. This shows that the masses of the most massive black holes are systematically higher than predicted by extrapolation of M \propto sigma^4 to high masses. This supports recent suggestions that stellar dynamical masses of the most massive black holes have been systematically underestimated. The steepening of the M-sigma relationship is consistent with the absence of very high sigma galaxies in the local universe and with the curvature of the Faber-Jackson relationship. There appears to be significantly less intrinsic scatter in the M-sigma relationship for galaxies with black hole masses > 10^9 solar masses. It is speculated that this is connected with the core elliptical versus extra-light elliptical dichotomy. The low scatter in the high end of the M-sigma relationship implies that the transformation proposed here and the Dibai method are good indicators of stellar velocity dispersion and mass respectively. There is no evidence for evolution of the M-sigma relationship over time.

We present the first results from a high sampling rate, multi-month reverberation mapping campaign undertaken primarily at MDM Observatory with supporting observations from telescopes around the world. The primary goal of this campaign was to obtain either new or improved Hbeta reverberation lag measurements for several relatively low luminosity AGNs. We feature results for NGC 4051 here because, until now, this object has been a significant outlier from AGN scaling relationships, e.g., it was previously a ~2-3sigma outlier on the relationship between the broad-line region (BLR) radius and the optical continuum luminosity - the R_BLR-L relationship. Our new measurements of the lag time between variations in the continuum and Hbeta emission line made from spectroscopic monitoring of NGC 4051 lead to a measured BLR radius of R_BLR = 1.87 (+0.54 -0.50) light days and black hole mass of M_BH = 1.73 (+0.55 -0.52) x 10^6 M_sun. This radius is consistent with that expected from the R_BLR-L relationship, based on the present luminosity of NGC 4051 and the most current calibration of the relation by Bentz et al. (2009a). We also present a preliminary look at velocity-resolved Hbeta light curves and time delay measurements, although we are unable to reconstruct an unambiguous velocity-resolved reverberation signal.

Two effects have substantially increased the scatter in the AGN black hole mass - host galaxy bulge luminosity relationship derived from SDSS spectra. The first is that at a fixed black hole mass, the SDSS spectrum depends strongly on redshift because an SDSS fiber sees a larger fraction of the total light of more distant galaxies. The second is that at a given redshift, the fraction of host-galaxy light in the fiber increases with decreasing galaxy luminosity. We illustrate the latter effect using the Kormendy et al. (2009) light profiles of Virgo ellipticals. With allowance for the two effects, we obtain a black hole mass - bulge luminosity relationship for AGNs which has a scatter of only +/- 0.23 dex in mass. This is less than the scatter found for inactive galaxies, and is consistent with the measuring errors. We show that there is a correspondingly tight linear relationship between the fraction of host galaxy light in AGN spectra and the Eddington ratio. This linearity implies that at a given black hole mass, the host luminosities of high-accretion-rate AGNs (NLS1s) and low-accretion-rate AGNs are similar. The relationship between the fraction of host galaxy light and the Eddington ratio provides a simple means of estimating the fraction of host galaxy light in AGN spectra. This means that the real amplitude of variability of low-accretion-rate AGNs is increased relative to NLS1s.

Some progress in understanding AGN variability is reviewed. Reprocessing of X-ray radiation to produce significant amounts of longer-wavelength continua seems to be ruled out. In some objects where there has been correlated X-ray and optical variability, the amplitude of the optical variability has exceeded the amplitude of X-ray variability. We suggest that accelerated particles striking material could be linking X-ray and optical variability (as in activity in the solar chromosphere). Beaming effects could be significant in all types of AGN. The diversity in optical/X-ray relationships at different times in the same object, and between different objects, might be explained by changes in geometry and directions of motion relative to our line of sight. Linear shot-noise models of the variability are ruled out; instead there must be large-scale organization of variability. Variability occurs on light-crossing timescales rather than viscous timescales and this probably rules out the standard Shakura-Sunyaev accretion disk. Radio-loud and radio-quiet AGNs have similar continuum shapes and similar variability properties. This suggests similar continuum origins and variability mechanisms. Despite their extreme X-ray variability, narrow-line Seyfert 1s (NLS1s) do not show extreme optical variability.

The first set of supermassive black hole mass estimates, published from 1980 to 1984 by E. A. Dibai, are shown to be in excellent agreement with recent reverberation-mapping estimates. Comparison of the masses of 17 AGNs covering a mass range from about 10^6 to 10^9 solar masses shows that the Dibai mass estimates agree with reverberation-mapping mass estimates to significantly better than 0.3 dex and were, on average, only 0.14 dex (~ 40%) systematically lower than masses obtained from reverberation mapping. This surprising agreement with the results of over a quarter of a century ago has important implication for the structure and kinematics of AGNs and implies that AGNs are very similar. Our results give strong support to the use of the single-epoch-spectrum (Dibai) method for investigating the co-evolution of supermassive black holes and their host galaxies.

We present the results of concurrent X-ray and optical monitoring of the Seyfert 1 galaxy Mrk 79 over a period of more than five years. We find that on short to medium time-scales (days to a few tens of days) the 2-10 keV X-ray and optical u and V band fluxes are significantly correlated, with a delay between the bands consistent with zero days. We show that most of these variations may be well reproduced by a model where the short-term optical variations originate from reprocessing of X-rays by an optically thick accretion disc. The optical light curves, however, also display long time-scale variations over thousands of days, which are not present in the X-ray light curve. These optical variations must originate from an independent variability mechanism and we show that they can be produced by variations in the (geometrically) thin disc accretion rate as well as by varying reprocessed fractions through changes in the location of the X-ray corona.

Gaskell et al. (2007) introduced polarization reverberation mapping as a new technique to explore the structure of active galactic nuclei. We present modeling results for the time-dependent polarization signal expected from scattering inside a centrally illuminated spheroid. Such a model setup describes a larger corona surrounding the compact source of an active nucleus. Time-delays between the polarized and the total flux are computed and related to the geometry of the cloud and the viewing angle. When including the in-flow dynamics of the cloud, it is possible to constrain its optical depth and velocity, which enables estimations of the mass inflow rate.

Energetics considerations imply that the broad-line region (BLR) has a high covering factor. The absence of absorption from the BLR means that the BLR has to have a flattened distribution and be seen through a polar hole. The BLR is the inward extension of the torus and they have similar geometries and covering factors. Reconciling velocity-resolved reverberation mapping, spectropolarimetry, and the increasing blueshifting of BLR lines with decreasing distance from the centre, implies that the BLR has a significant inflow component. This inflow provides the mass inflow rate needed to power the AGN. We suggest that the mechanism producing the outward transport of angular momentum necessary for the net inflow of the BLR is the magneto-rotational instability, and that the BLR and outer accretion disc are one and the same.

We present colour transformations for the conversion of the {\em 2MASS} photometric system to the Johnson-Cousins $UBVRI$ system and further into the {\em SDSS} $ugriz$ system. We have taken {\em SDSS} $gri$ magnitudes of stars measured with the 2.5-m telescope from $SDSS$ Data Release 5 (DR5), and $BVRI$ and $JHK_{s}$ magnitudes from Stetson's catalogue and \citet{Cu03}, respectively. We matched thousands of stars in the three photometric systems by their coordinates and obtained a homogeneous sample of 825 stars by the following constraints, which are not used in previous transformations: 1) the data are de-reddened, 2) giants are omitted, and 3) the sample stars selected are of the highest quality. We give metallicity, population type, and transformations dependent on two colours. The transformations provide absolute magnitude and distance determinations which can be used in space density evaluations at short distances where some or all of the {\em SDSS} $ugriz$ magnitudes are saturated. The combination of these densities with those evaluated at larger distances using {\em SDSS} $ugriz$ photometry will supply accurate Galactic model parameters, particularly the local space densities for each population.

Theory and observations of the dominant thermal continuum emission in AGNs are examined. After correction for reddening, the steady state AGN optical--UV spectral energy distributions (SEDs) are very similar. The SEDs are dominated energetically by the big blue bump (BBB), but this bump never shows the nu^{+1/3} spectrum predicted for a standard thin accretion disk with a r^{-0.75} radial temperature gradient. Instead, the observed optical-UV SED implies a temperature gradient of r^{-0.57} independent of the thickness of the disk. This means that there is some flow of heat outwards in the disk. The disk is large and the region emitting the optical continuum is as large as the inner broad-line region (BLR). Because optical variability is seen in all AGNs on the light-crossing time of the BLR, variations must propagate at close to the speed of light, rather than on dynamical timescales. This argues that the energy-generation mechanism is electromagnetic rather that hydrodynamic. Since the velocities are near the speed of light, there can be significant local anisotropy in the emission. The large rapid variations of the BBB imply that the magnetohydrodynamic energy generation is fundamentally unstable. Because of the inevitable radial temperature gradient in the accreting material, different spectral regions come predominantly from different radii, and variations in different spectral regions correspond to variability at different radii. This explains the frequently observed independence of X-ray and optical variations, cases of variability at lower energies leading variability at higher energies, and rapid changes in emission-line reverberation lags. Some observational tests of the local variability hypothesis are proposed.

We consider in detail the spectral energy distribution (SED) and multi-wavelength variability of NGC5548. Comparison with the SEDs of other AGNs implies that the internal reddening of NGC5548 is E(B-V) = 0.17 mag. The extinction curve is consistent with the mean curve of other AGNs found by Gaskell & Benker, but inconsistent with an SMC-type reddening curve. Because most IR emission originates exterior to the broad-line region (BLR), the SED seen by the inner BLR is different from that seen by the outer BLR and from the earth. The most likely BLR covering factor is ~ 40% and it is not possible to get an overall BLR covering factor of less than 20%. This requires that the BLR is not spherically symmetric and that we are viewing through a hole. Line-continuum variability transfer functions are consistent with this geometry. The covering factor and geometry imply that near the equatorial plane the BLR covering approaches 100%. The spectrum seen by the outer regions of the BLR and by the torus is thus modified by the absorption in the inner BLR. This shielding solves the problem of observed BLR ionization stratification being much greater than implied by photoionization models. The BLR obscuration also removes the problem of the torus covering factor being greater than the BLR covering factor, and gives consistency with the observed fraction of obscured AGNs. The flux reduction at the torus also reduces the problem of AGN dust-reverberation lags giving sizes smaller than the dust-sublimation radii.

We present intrinsic extinction curves for 14 AGNs. The AGNs have reddenings, E(B-V), of up to 0.36 mag. The majority (13 out of 14) of the extinction curves are not steep in the UV. Of the seven best determined extinction curves, five have extinction curves that are as flat as the standard Galactic curve in the optical and near UV, but flatter in the far UV, and without the 2175 Angstrom feature. One AGN, B3 0754+394, has a steep SMC-like extinction curve, and another, Mrk 304, has an LMC-like extinction curve, including a probable 2175 Ang. bump. The remaining seven, lower-quality, extinction curves have overall shapes that are consistent with an LMC-like shape or a flatter shape. Two have possible 2175 features, and one might be identical to the Galactic curve. The flatter curves that predominate in our best determined extinction curves are not as flat as the Gaskell et al. (2004) extinction curve for radio-loud AGNs. This suggests that the previous radio-loud extinction curve might be slightly too flat in the range 4 < 1/lambda < 6.5 inverse microns because of luminosity-dependent reddening biases in the composite spectra, but further investigation is needed. We present a parameterized average AGN extinction curve. Observed variations in the continuum properties of the AGNs are inconsistent with intrinsic object-to-object variations because observed differences are least in the far UV where changes in the accretion disk spectrum should be greatest. We suggest that the steepening of AGN spectra around Lyman alpha is the result of a small amount of SMC-like dust (E(B-V) ~ 0.03). We find the largest object-to-object differences in spectral shape to be in the Fe II emission of the "small blue bump".

We discuss the use of Sloan Digital Sky Survey (SDSS) ugriz point-spread function (PSF) photometry for setting the zero points of UBVRI CCD images. From a comparison with the Landolt (1992) standards and our own photometry we find that there is a fairly abrupt change in B, V, R, & I zero points around g, r, i ~ 14.5, and in the U zero point at u ~ 16. These changes correspond to where there is significant interpolation due to saturation in the SDSS PSF fluxes. There also seems to be another, much smaller systematic effect for stars with g, r > 19.5. The latter effect is consistent with a small Malmquist bias. Because of the difficulties with PSF fluxes of brighter stars, we recommend that comparisons of ugriz and UBVRI photometry should only be made for unsaturated stars with g, r and i in the range 14.5 - 19.5, and u in the range 16 - 19.5. We give a prescription for setting the UBVRI zero points for CCD images, and general equations for transforming from ugriz to UBVRI.

We have used the STOKES radiative transfer code, to model polarization induced by dust scattering in the polar regions of Active Galactic Nuclei (AGN). We discuss the wavelength-dependence of the spectral intensity and polarization over the optical/UV range at different viewing angles for two different types of dust: a Galactic dust model, and a dust model inferred from extinction properties of AGN. The STOKES code and documentation are freely available at http://www.stokes-program.info/.

We introduce a new, publicly available Monte Carlo radiative transfer code, STOKES, which has been developed to model polarization induced by scattering off free electrons and dust grains. It can be used in a wide range of astrophysical applications. Here, we apply it to model the polarization produced by the equatorial obscuring and scattering tori assumed to exist in active galactic nuclei (AGNs). We present optical/UV modeling of dusty tori with a curved inner shape and for two different dust types: one composition reproduces extinction properties of our Galaxy, and the other is derived from composite quasar spectra. The polarization spectra enable us to clearly distinguish between the two dust compositions. The STOKES code and its documentation can be freely downloaded from http://www.stokes-program.info/.

[abridged] We investigate the effects of various AGN scattering region geometries on the polarized flux. We introduce a new, publicly available Monte Carlo radiative transfer code, Stokes, which models polarization induced by scattering off free electrons and dust grains. We find that the shape of the funnel of the dusty torus has a significant impact on the polarization efficiency. A compact torus with a steep inner surface scatters more light toward type-2 viewing angles than a large torus of the same half-opening angle, theta0. For theta0 < 53 deg, the scattered light is polarized perpendicularly to the symmetry axis, whilst for theta0 > 60 deg it is polarized parallel to the symmetry axis. In between these intervals the orientation of the polarization depends on the viewing angle. The degree of polarization ranges between 0% and 20% and is wavelength-independent for a large range of theta0. Observed wavelength-independent optical and near-UV polarization thus does not necessarily imply electron scattering. For polar dust, scattering spectra are reddened for type-1 viewing angles, and made bluer for type-2 viewing angles. Polar electron-scattering cones are very efficient polarizers at type-2 viewing angles, whilst the polarized flux of the torus is weak. We predict that the net polarization of Seyfert-2 galaxies decreases with luminosity, and conclude that the degree of polarization should be correlated with the relative strength of the thermal IR flux. We find that a flattened, equatorial, electron-scattering disk, of relatively low optical depth, reproduces type-1 polarization. This is insensitive to the exact geometry, but the observed polarization requires a limited range of optical depth.

I present new comparisons of AGN optical, UV, and X-ray variations. These reveal complex relationships between the different passbands that can change with time in a given object. While there is evidence from several objects that X-ray and optical activity levels are correlated on long timescales, variations on shorter timescales can occur independently. It is proposed that the combination of correlated and uncorrelated short-timescale variability is a consequence of anisotropic high-energy emission. It is also argued that the correlation between X-ray and optical variability on long timescales must be due to a common underlying factor and not to reprocessing of X-ray radiation.

We report on the relative optical variability of the three brightest nearby quasars, 3C 273, PDS 456, and PHL 1811. All three have comparable absolute magnitudes, but PDS 456 and PHL 1811 are radio quiet. PDS 456 is a broad-line object, but PHL 1811 could be classified as a high-luminosity Narrow-Line Seyfert 1 (NLS1). Both of the radio-quiet quasars show significant variability on a timescale of a few days. The seasonal rms V-band variability amplitudes of 3C 273 and PDS 456 are indistinguishable, and the seasonal rms variability amplitude of PHL 1811 was only exceeded by 3C 273 once in 30 years of monitoring. We find no evidence that the optical variability of 3C 273 is greater than or more rapid than the variability of the comparably-bright, radio-quiet quasars. This suggests that not only do radio-loud and radio-quiet AGNs have similar spectral energy distributions, but that the variability mechanisms are also similar. The optical variability of 3C 273 is not dominated by a "blazer" component.

A model of wavelength-dependent lags in optical continuum variability of AGNs is proposed which avoids the problems of the popular ``lamp-post'' model. Rather than being due to reprocessing of high-energy radiation from a hypothetical source above the accretion disk, the wavelength-dependent delays observed from the B to I bands are instead due to contamination of an intrinsically coherently variable continuum with the Wien tail of the thermal emission from the hot dust in the surrounding torus. The new model correctly gives the size, wavelength dependence, and luminosity dependence of the lags, and quantitatively predicts observed color hysteresis. The model also explains how the measured delays vary with epoch of observation. There must also be contamination by scattered light and this can be detected by a lag in the polarized flux.

We present the first reliable determination of chemical abundances in an AGN outflow. The abundances are extracted from the deep and simultaneous FUSE and HST/STIS observations of Mrk 279. This data set is exceptional for its high signal-to-noise, unblended doublet troughs and little Galactic absorption contamination. These attributes allow us to solve for the velocity-dependent covering fraction, and therefore obtain reliable column densities for many ionic species. For the first time we have enough such column densities to simultaneously determine the ionization equilibrium and abundances in the flow. Our analysis uses the full spectral information embedded in these high-resolution data. Slicing a given trough into many independent outflow elements yields the extra constraints needed for a physically meaningful abundances determination. We find that relative to solar the abundances in the Mrk 279 outflow are (linear scaling): carbon 2.2+/-0.7, nitrogen 3.5+/-1.1 and oxygen 1.6+/-0.8. Our UV-based photoionization and abundances results are in good agreement with the independent analysis of the simultaneous Mrk 279 X-ray spectra. This is the best agreement between the UV and X-ray analyses of the same outflow to date.

We present low-resolution ultraviolet spectra of 14 low redshift (z<0.8) quasars observed with HST/STIS as part of a Snap project to understand the relationship between quasar outflows and luminosity. By design, all observations cover the CIV emission line. Nine of the quasars are from the Hamburg-ESO catalog, three are from the Palomar-Green catalog, and one is from the Parkes catalog. The sample contains a few interesting quasars including two broad absorption line (BAL) quasars (HE0143-3535, HE0436-2614), one quasar with a mini-BAL (HE1105-0746), and one quasar with associated narrow absorption (HE0409-5004). These BAL quasars are among the brightest known (though not the most luminous) since they lie at z<0.8. We compare the properties of these BAL quasars to the z<0.5 Palomar-Green and z>1.4 Large Bright Quasar samples. By design, our objects sample luminosities in between these two surveys, and our four absorbed objects are consistent with the v ~ L^0.62 relation derived by Laor & Brandt (2002). Another quasar, HE0441-2826, contains extremely weak emission lines and our spectrum is consistent with a simple power-law continuum. The quasar is radio-loud, but has a steep spectral index and a lobe-dominated morphology, which argues against it being a blazar. The unusual spectrum of this quasar resembles the spectra of the quasars PG1407+265, SDSSJ1136+0242, and PKS1004+13 for which several possible explanations have been entertained.

We present a Chandra-LETGS observation of the Seyfert 1 galaxy Mrk 279. This observation was carried out simultaneously with HST-STIS and FUSE, in the context of a multiwavelength study of this source. The Chandra pointings were spread over ten days for a total exposure time of ~360 ks. The spectrum of Mrk279 shows evidence of broad emission features, especially at the wavelength of the OVII triplet. We quantitatively explore the possibility that this emission is produced in the broad line region (BLR). We modeled the broad UV emission lines seen in the FUSE and HST-STIS spectra following the ``locally optimally emitting cloud" approach. We find that the X-ray lines luminosity derived from the best fit BLR model can match the X-ray features, suggesting that the gas producing the UV lines is sufficient to account also for the X-ray emission. The spectrum is absorbed by ionized gas whose total column density is ~5x10^{20} cm^{-2}. The absorption spectrum can be modeled by two distinct gas components (log xi ~ 0.47 and 2.49, respectively) both showing a significant outflow velocity. However, the data allow also the presence of intermediate ionization components. The distribution of the column densities of such extra components as a function of the ionization parameter is not consistent with a continuous, power law-like, absorber, suggesting a complex structure for the gas outflow for Mrk 279 (abridged).

We present an analysis of the intrinsic UV absorption in the Seyfert 1 galaxy Mrk 279 based on simultaneous long observations with the Hubble Space Telescope (41 ks) and the Far Ultraviolet Spectroscopic Explorer (91 ks). To extract the line-of-sight covering factors and ionic column densities, we separately fit two groups of absorption lines: the Lyman series and the CNO lithium-like doublets. For the CNO doublets we assume that all three ions share the same covering factors. The fitting method applied here overcomes some limitations of the traditional method using individual doublet pairs; it allows for the treatment of more complex, physically realistic scenarios for the absorption-emission geometry and eliminates systematic errors that we show are introduced by spectral noise. We derive velocity-dependent solutions based on two models of geometrical covering -- a single covering factor for all background emission sources, and separate covering factors for the continuum and emission lines. Although both models give good statistical fits to the observed absorption, we favor the model with two covering factors because: (a) the best-fit covering factors for both emission sources are similar for the independent Lyman series and CNO doublet fits; (b) the fits are consistent with full coverage of the continuum source and partial coverage of the emission lines by the absorbers, as expected from the relative sizes of the nuclear emission components; and (c) it provides a natural explanation for variability in the Ly$\alpha$ absorption detected in an earlier epoch. We also explore physical and geometrical constraints on the outflow from these results.

We present extinction curves derived from the broad emission lines and continua of large samples of both radio-loud and radio-quiet AGNs. The curves are significantly flatter in the UV than are curves for the local ISM. The reddening curves for the radio-quiet LBQS quasars are slightly steeper than those of the radio-loud quasars in the UV, probably because of additional reddening by dust further out in the host galaxies of the former. The UV extinction curves for the radio-loud AGNs are very flat. This is explicable with slight modifications to standard MRN dust models: there is a relative lack of small grains in the nuclear dust. Our continuum and broad-emission line reddening curves agree in both shape and amplitude, confirming that the continuum shape is indeed profoundly affected by reddening for all but the bluest AGNs. With correction by our generic extinction curve, all of the radio-loud AGNs have continuous optical-UV spectra consistent with a single shape. We show that radio-quiet AGNs have very similar intrinsic UV to optical shape over orders of magnitude in luminosity. We also argue that radio-loud and radio-quiet AGNs probably share the same underlying continuum shape and that most of the systematic differences between their observed continuum shapes are due to higher nuclear reddening in radio-selected AGNs, and additional reddening from dust further out in the host galaxies in radio-quiet AGNs. Our conclusions have important implications for the modelling of quasar continua and the analysis of quasar demographics.

We present results of a broad-band photometric study of the optical variability of six Narrow-Line Seyfert 1 (NLS1) galaxies observed at 172 epochs. We searched for microvariability on 33 nights. Strong evidence for microvariability was found only for our lowest luminosity object, NGC 4051, on one night. Weaker evidence suggests such variability on a few other nights for two other objects, but the data are not as convincing. Intra-night variability in NLS1s is thus concluded to be rare and of low amplitude. We give illustrations of how variable image quality can produce spurious variability. We find that for well-studied non-NLS1s there is a spread in the amplitude of seasonal variability (i.e., in some years an AGN is more variable than in others). We find that the means of the variability amplitudes of non-NLS1s over several seasons vary from object to object (i.e., some AGNs are, on average, more variable than others). NLS1s also show a spread in seasonal variabilities. The best-studied NLS1, Ark 564, shows a range of amplitudes of variability from season to season that is comparable to the range found in BLS1s, and in one season Ark 564 was as variable as the most variable non-NLS1. The seasonal amplitudes of variability for NLS1s are mostly in the lower half of the range of non-BLS1 seasonal amplitudes, but the absence of a suitable control sample makes a precise comparison difficult. However, on long timescales (weeks to years), NLS1s as a class are not more variable than non-NLS1s. The extreme variability seen in the X-rays was not seen in the optical. This has consequences for the models of AGNs in general as well as NLS1s in particular.

We present results from a study of high- and low-ionization emission line ratios as a function of projected velocity for a sample of eight active galactic nuclei (AGNs). Our results are based on analysis of high signal-to-noise optical and Hubble Space Telescope (HST) UV spectra. Comparing the emission line ratios to those predicted by photoionization models indicates that the physical conditions responsible for the high-ionization emission lines are consistent with a wind, whereas those of the low-ionization lines are consistent with a virialized disk.

A combined analysis of the profiles of the main broad quasar emission lines in both Hubble Space Telescope and optical spectra shows that while the profiles of the strong UV lines are quite similar, there is frequently a strong increase in the Ly-alpha/H-alpha ratio in the high-velocity gas. We show that the suggestion that the high velocity gas is optically-thin presents many problems. We show that the relative strengths of the high velocity wings arise naturally in an optically-thick BLR component. An optically-thick model successfully explains the equivalent widths of the lines, the Ly-alpha/H-alpha ratios and flatter Balmer decrements in the line wings, the strengths of CIII] and the lambda 1400 blend, and the strong variability of high-velocity, high-ionization lines (especially HeII and HeI).

We present the results of a two-year long optical monitoring program of the narrow-line Seyfert 1 galaxy Akn 564. The majority of this monitoring project was also covered by X-ray observations (RXTE) and for a period of ~50 days, we observed the galaxy in UV (HST) and X-rays (RXTE & ASCA) simultaneously with the ground-based observations. Rapid and large-amplitude variations seen in the X-ray band, on a daily and hourly time-scale, were not detected at optical and UV wavelengths, which in turn exhibited much lower variability either on short (one day) or long (several months) time-scales. The only significant optical variations can be described as two 2--4 day events with ~10% flux variations. We detect no significant optical line variations and thus cannot infer a reverberation size for the broad-line region. Similarly, the large X-ray variations seem to vanish when the light curve is smoothed over a period of 30 days. The UV continuum follows the X-rays with a lag of ~0.4 days, and the optical band lags the UV band by ~2 days. No significant correlation was found between the entire X-ray dataset and the optical band. Focusing on a 20-day interval around the strongest optical event we detect a significant X-ray--optical correlation with similar events seen in the UV and X-rays. Our data are consistent with reprocessing models on the grounds of the energy emitted in this single event. However, several large X-ray flares produced no corresponding optical emission.

Results of a ground-based optical monitoring campaign on NGC5548 in June 1998 are presented. The broad-band fluxes (U,B,V), and the spectrophotometric optical continuum flux F_lambda(5100 A) monotonically decreased in flux while the broad-band R and I fluxes and the integrated emission-line fluxes of Halpha and Hbeta remained constant to within 5%. On June 22, a short continuum flare was detected in the broad band fluxes. It had an amplitude of about ~18% and it lasted only ~90 min. The broad band fluxes and the optical continuum F_lambda(5100 A) appear to vary simultaneously with the EUV variations. No reliable delay was detected for the broad optical emission lines in response to the EUVE variations. Narrow Hbeta emission features predicted as a signature of an accretion disk were not detected during this campaign. However, there is marginal evidence for a faint feature at lambda = 4962 A with FWHM=~6 A redshifted by Delta v = 1100 km/s with respect to Hbeta_narrow.

The evidence is reviewed that the Broad Line Region (BLR) probably has two distinct components located at about the same distance from the central black hole. One component, BLR II, is optically-thick, low-ionization emission at least some of which arises from a disc and the other, BLR I, is probably optically-thin emission from a more spherically symmetric halo or atmosphere. The high Fe II/H-beta ratios seen in Narrow-Line Seyfert 1 galaxies (NLS1s) are not due to strong Fe II emission, as is commonly thought, but to unusually weak Balmer emission, probably caused by higher densities. NLS1s probably differ from non-NLS1s because of the higher density of gas near the black hole. This produces a higher accretion rate, a denser BLR, and a view of the central regions that is more face-on.

Some quasars show Doppler shifted broad emission line peaks. I give new statistics of the occurrence of these peaks and show that, while the most spectacular cases are in quasars with strong radio jets inclined to the line of sight, they are also almost as common in radio-quiet quasars. Theories of the origin of the peaks are reviewed and it is argued that the displaced peaks are most likely produced by the supermassive binary model. The separations of the peaks in the 3C 390.3-type objects are consistent with orientation-dependent "unified models" of quasar activity. If the supermassive binary model is correct, all members of "the jet set" (astrophysical objects showing jets) could be binaries.

Analysis of spectra of the quasar 3C 390.3 covering a period of over 20 yr shows that the blueshifted peak of H beta has been changing its radial velocity at an almost constant rate during this time. The radial velocity has increased by over 1500 km/s. The lower limit to the period of radial velocity changes is 210 yr. Although very long periods cannot be excluded by the radial velocity curve alone, other considerations suggest that the period is about 300 yr. If the radial velocity changes are due to orbital motion, the radius of the orbit is about 0.3 pc and the total mass of the system is about 7 * 10**9 solar masses. In the binary black-hole model the masses of the two holes are about 2.2 * 10**9 solar masses. and about 4.4 * 10**9 solar masses. A possible third peak claimed to be present in some 1974-75 spectra is shown to be an instrumental artifact. The narrowness of the displaced peaks in 3C 390.3 objects requires that the broad line region cloud motions are not governed solely by gravity. The black hole masses derived by Koratkar & Gaskell (1991) need to be increased significantly. This probably lowers accretion efficiences to less than 1\% of the Eddington limit.