We present a code for modelling the ionization conditions of optically thin astrophysical gas structures. Given the gas hydrogen density, equilibrium temperature, elemental abundances, and the ionizing spectrum, the code solves the equilibrium ionization fractions and number densities for all ions from hydrogen to zinc. The included processes are photoionization, Auger ionization, direct collisional ionization, excitation auto-ionization, charge exchange ionization, two-body radiative recombination, dielectronic recombination, and charge exchange recombination. The ionizing spectrum can be generalized to include the ultraviolet background (UVB) and/or Starburst99 stellar populations of various masses, ages, metallicities, and distances. The ultimate goal with the code is to provide fast computation of the ionization conditions of gas in N-body + hydrodynamics cosmological simulations, in particular adaptive mesh refinement codes, in order to facilitate absorption line analysis of the simulated gas for comparison with observations. To this end, we developed a method to linearize the rate equations and efficiently solve the rate matrix with a minimum number of iterations. Comparing the code to Cloudy 13.03 (Ferland), we find that the hydrogen and helium ionization fractions and the metal species ionization corrections are in excellent agreement. We discuss the science drivers and plans for further development of the ionization code to a full radiative hydrodynamic routine that can be employed for processing the simulations. A stand-alone version of the code has been made publicly available (http://astronomy.nmsu.edu/cwc/Software/Ioncode/).

We present the first observation of a galaxy (z=0.2) that exhibits metal-line absorption back-illuminated by the galaxy ("down-the-barrel") and transversely by a background quasar at a projected distance of 58 kpc. Both absorption systems, traced by MgII, are blueshifted relative to the galaxy systemic velocity. The quasar sight-line, which resides almost directly along the projected minor axis of the galaxy, probes MgI and MgII absorption obtained from Keck/LRIS and Lya, SiII and SiIII absorption obtained from HST/COS. For the first time, we combine two independent models used to quantify the outflow properties for down-the-barrel and transverse absorption. We find that the modeled down-the-barrel deprojected outflow velocities range between $V_{dtb}=45-255$ km/s. The transverse bi-conical outflow model, assuming constant-velocity flows perpendicular to the disk, requires wind velocities $V_{outflow}=40-80$ km/s to reproduce the transverse MgII absorption kinematics, which is consistent with the range of $V_{dtb}$. The galaxy has a metallicity, derived from H$\alpha$ and NII, of $[{\rm O/H}]=-0.21\pm0.08$, whereas the transverse absorption has $[{\rm X/H}]=-1.12\pm0.02$. The galaxy star-formation rate is constrained between $4.6-15$ M$_{\odot}$/yr while the estimated outflow rate ranges between $1.6-4.2$ M$_{\odot}$/yr and yields a wind loading factor ranging between $0.1-0.9$. The galaxy and gas metallicities, the galaxy-quasar sight-line geometry, and the down-the-barrel and transverse modeled outflow velocities collectively suggest that the transverse gas originates from ongoing outflowing material from the galaxy. The $\sim$1 dex decrease in metallicity from the base of the outflow to the outer halo suggests metal dilution of the gas by the time it reached 58 kpc.

Any successful model of galaxy formation needs to explain the low rate of star formation in the small progenitors of today's galaxies. This inefficiency is necessary for reproducing the low stellar-to-virial mass fractions, suggested by current abundance matching models. A possible driver of this low efficiency is the radiation pressure exerted by ionizing photons from massive stars. The effect of radiation pressure in cosmological, zoom-in galaxy formation simulations is modeled as a non-thermal pressure that acts only in dense and optically thick star-forming regions. We also include photoionization and photoheating by massive stars. The full photoionization of hydrogen reduces the radiative cooling in the $10^{4-4.5}$ K regime. The main effect of radiation pressure is to regulate and limit the high values of gas density and the amount of gas available for star formation. This maintains a low star formation rate of $\sim 1 \ {\rm M_\odot} \ {\rm yr}^{-1}$ in halos with masses about $10^{11} \ {M_\odot}$ at $z\simeq3$. Infrared trapping and photoionization/photoheating processes are secondary effects in this mass range. The galaxies residing in these low-mass halos contain only $\sim0.6\%$ of the total virial mass in stars, roughly consistent with abundance matching. Radiative feedback maintains an extended galaxy with a rising circular velocity profile.

We have identified 469 MgII doublet systems having W_r >= 0.02 {\AA} in 252 Keck/HIRES and UVES/VLT quasar spectra over the redshift range 0.1 < z < 2.6. Using the largest sample yet of 188 weak MgII systems (0.02 {\AA} <= W_r < 0.3 {\AA}), we calculate their absorber redshift path density, dN/dz. We find clear evidence of evolution, with dN/dz peaking at z ~ 1.2, and that the product of the absorber number density and cross section decreases linearly with increasing redshift; weak MgII absorbers seem to vanish above z ~ 2.7. If the absorbers are ionized by the UV background, we estimate number densities of 10^6 - 10^9 per Mpc^3 for spherical geometries and 10^2 - 10^5 per Mpc^3 for more sheetlike geometries. We also find that dN/dz toward intrinsically faint versus bright quasars differs significantly for weak and strong (W_r >= 1.0 {\AA}) absorbers. For weak absorption, dN/dz toward bright quasars is ~ 25% higher than toward faint quasars (10 sigma at low redshift, 0.4 <= z <= 1.4, and 4 sigma at high redshift, 1.4 < z <= 2.34). For strong absorption the trend reverses, with dN/dz toward faint quasars being ~ 20% higher than toward bright quasars (also 10 sigma at low redshift and 4 sigma at high redshift). We explore scenarios in which beam size is proportional to quasar luminosity and varies with absorber and quasar redshifts. These do not explain dN/dz's dependence on quasar luminosity.

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.

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.

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.

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.

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 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 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.

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 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 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.