Tidal disruption events involve numerous physical processes (fluid dynamics, magnetohydrodynamics, radiation transport, self-gravity, general relativistic dynamics) in highly nonlinear ways, and, because TDEs are transients by definition, frequently in non-equilibrium states. For these reasons, numerical solution of the relevant equations can be an essential tool for studying these events. In this chapter, we present a summary of the key problems of the field for which simulations offer the greatest promise and identify the capabilities required to make progress on them. We then discuss what has been---and what cannot be---done with existing numerical methods. We close with an overview of what methods now under development may do to expand our ability to understand these events.

The black hole of an active galactic nucleus is encircled by an accretion disk. The surface density of the disk is always too low to affect the tidal disruption of a star, but it can be high enough that a vigorous interaction results when the debris stream returns to pericenter and punches through the disk. Shocks excited in the disk dissipate the kinetic energy of the disk interior to the impact point and expedite inflow toward the black hole. Radiatively efficient disks with luminosity $\gtrsim10^{-3}$ Eddington have high enough surface density that the initial stream-disk interaction leads to energy dissipation at a super-Eddington rate. Because of the rapid inflow, only part of this dissipated energy emerges as radiation, while the rest is advected into the black hole. Dissipation, inflow, and cooling balance to keep the bolometric luminosity at an Eddington-level plateau that lasts tens of days. After the plateau, the luminosity decreases in proportion to the disk surface density, with a power-law index between $-3$ and $-2$ at earlier times, and possibly a steeper index at later times.

A fraction of tidal disruption events (TDEs) occur in active galactic nuclei (AGNs) whose black holes possess accretion disks; these TDEs can be confused with common AGN flares. The disruption itself is unaffected by the disk, but the evolution of the bound debris stream is modified by its collision with the disk when it returns to pericenter. The outcome of the collision is largely determined by the ratio of the stream mass current to the azimuthal mass current of the disk rotating underneath the stream footprint, which in turn depends on the mass and luminosity of the AGN. To characterize TDEs in AGNs, we simulated a suite of stream--disk collisions with various mass current ratios. The collision excites shocks in the disk, leading to inflow and energy dissipation orders of magnitude above Eddington; however, much of the radiation is trapped in the inflow and advected into the black hole, so the actual bolometric luminosity may be closer to Eddington. The emergent spectrum may not be thermal, TDE-like, or AGN-like. The rapid inflow causes the disk interior to the impact point to be depleted within a fraction of the mass return time. If the stream is heavy enough to penetrate the disk, part of the outgoing material eventually hits the disk again, dissipating its kinetic energy in the second collision; another part becomes unbound, emitting synchrotron radiation as it shocks with surrounding gas.

This paper introduces a series of papers presenting a quantitative theory for the tidal disruption of main sequence stars by supermassive black holes. Using fully general relativistic hydrodynamics simulations and MESA-model initial conditions, we explore the pericenter-dependence of tidal disruption properties for eight stellar masses ($0.15 \leq M_*/M_\odot \leq 10$) and six black hole masses ($10^5 \leq M_{BH}/M_\odot \leq 5 \times 10^7$). We present here the results most relevant to observations. The effects of internal stellar structure and relativity decouple for both the disruption cross section and the characteristic energy width of the debris. Moreover, the full disruption cross section is almost independent of $M_*$ for $M_*/M_\odot \lesssim 3$. Independent of $M_*$, relativistic effects increase the critical pericenter distance for full disruptions by up to a factor $\sim 3$ relative to the Newtonian prediction. The probability of a direct capture is also independent of $M_*$; at $M_{BH}/M_\odot \simeq 5 \times 10^6$ this probability is equal to that of a complete disruption. The width of the debris energy distribution $\Delta E$ can differ from the standard estimate by factors from 0.35 to 2, depending on $M_*$ and $M_{BH}$, implying a corresponding change in the characteristic mass-return timescale. The "frozen-in approximation" is inconsistent with $\Delta E$, and mass-loss continues over a long span of time. We provide analytic forms, suitable for use in both event rate estimates and parameter inference, to describe all these trends. For partial disruptions, we find a nearly-universal relation between the star's angular momentum and the fraction of $M_*$ remaining. Within the "empty loss-cone" regime, partial disruptions must precede full disruptions. These partial disruptions can drastically affect the rate and appearance of subsequent total disruptions.

This is the second in a series of papers presenting the results of fully general relativistic simulations of stellar tidal disruptions in which the stars' initial states are realistic main-sequence models. In the first paper (Paper I), we gave an overview of this program and discussed the principal observational implications of our work. Here we describe our calculational method and provide details about the outcomes of full disruptions, focusing on the stellar mass dependence of the outcomes for a black hole of mass $10^{6}\rm{M}_{\odot}$. We consider eight different stellar masses, from $0.15~{\rm M}_\odot$ to $10~{\rm M}_\odot$. We find that, relative to the traditional order-of-magnitude estimate $r_{\rm t}$, the physical tidal radius of low-mass stars ($M_{\star} \lesssim 0.7~ {\rm M}_\odot$) is larger by tens of percent, while for high-mass stars ($M_{\star} \gtrsim1~ {\rm M}_\odot$) it is smaller by a factor 2--2.5. The traditional estimate of the range of energies found in the debris is $\approx 1.4\times$ too large for low-mass stars, but is a factor $\sim 2$ too small for high-mass stars; in addition, the energy distribution for high-mass stars has significant wings. For all stars undergoing tidal encounters, we find that mass-loss continues for many stellar vibration times because the black hole's tidal gravity competes with the instantaneous stellar gravity at the star's surface until the star has reached a distance from the black hole $\sim O(10)r_{\rm t}$.

In this paper, the third in this series, we continue our study of tidal disruption events of main-sequence stars by a non-spinning $10^{6}~\rm{M}_\odot$ supermassive black hole. Here we focus on the stellar mass dependence of the outcomes of partial disruptions. As the encounter becomes weaker, the debris mass is increasingly concentrated near the outer edges of the energy distribution. As a result, the mass fallback rate can deviate substantially from a $t^{-5/3}$ power-law, becoming more like a single peak with a tail declining as $t^{-p}$ with $p\simeq2-5$. Surviving remnants are spun-up in the prograde direction and are hotter than main sequence stars of the same mass. Their specific orbital energy is $\simeq10^{-3}\times$ that of the debris, but of either sign with respect to the black hole potential, while their specific angular momentum is close to that of the original star. Even for strong encounters, remnants have speeds at infinity relative to the black hole potential $\lesssim 300$ km s$^{-1}$, so they are unable to travel far out into the galactic bulge. The remnants most deeply bound to the black hole go through a second tidal disruption event upon their first return to pericenter; if they have not thermally relaxed, they will be completely disrupted.

Using a suite of fully relativistic hydrodynamic simulations applied to main-sequence stars with realistic internal density profiles, we examine full and partial tidal disruptions across a wide range of black hole mass ($10^{5}\leq M_{\rm BH}/\mathrm{M}_{\odot}\leq 5\times 10^{7}$) and stellar mass ($0.3 \leq M_{\star} /\mathrm{M}_{\odot}\leq 3$) as larger $M_{\rm BH}$ leads to stronger relativistic effects. For fixed $M_{\star}$, as $M_{\rm BH}$ increases, the ratio of the maximum pericenter distance yielding full disruptions ($\mathcal{R}_{\rm t}$) to its Newtonian prediction rises rapidly, becoming triple the Newtonian value for $M_{\rm BH} = 5\times10^{7}~{\rm M}_\odot$, while the ratio of the energy width of the stellar debris for full disruptions to the Newtonian prediction decreases steeply, resulting in a factor of two correction at $M_{\rm BH} = 5 \times 10^7~{\rm M}_\odot$. We find that for partial disruptions, the fractional remnant mass for a given ratio of the pericenter to $\mathcal{R}_{\rm t}$ is higher for larger $M_{\rm BH}$. These results have several implications. As $M_{\rm BH}$ increases above $\sim 10^7~{\rm M}_\odot$, the cross section for complete disruptions is suppressed by competition with direct capture. However, the cross section ratio for partial to complete disruptions depends only weakly on $M_{\rm BH}$. The relativistic correction to the debris energy width delays the time of peak mass-return rate and diminishes the magnitude of the peak return rate. For $M_{\rm BH} \gtrsim 10^7~{\rm M}_\odot$, the $M_{\rm BH}$-dependence of the full disruption cross section and the peak mass-return rate and time is influenced more by relativistic effects than by Newtonian dynamics.

Tidal disruption events (TDEs), events in which a star passes very close to a supermassive black hole, are generally imagined as leading either to the star's complete disruption or to its passage directly into the black hole. In the former case it is widely believed that in all cases the bound portion of the debris quickly "circularizes" due to relativistic apsidal precession, i.e., forms a compact accretion disk, and emits a flare of standardized lightcurve and spectrum. We show here that TDEs are more diverse and can be grouped into several distinct categories on the basis of stellar pericenter distance $r_p$; we calculate the relative frequency of these categories. In particular, because rapid circularization requires $r_p \lesssim 10r_g$ ($r_g \equiv GM_{\rm BH}/c^2$), it can happen in only a minority of total disruptions, $\lesssim 1/4$ when the black hole has mass $M_{\rm BH} = 10^6 M_\odot$. For larger pericenter distances, $10 < r_p/r_g < 27$ (for $M_{\rm BH}=10^6M_\odot$), main sequence stars are completely disrupted, but the bound debris orbits are highly eccentric and possess semimajor axes $\sim 100\times$ the scale of the expected compact disk. Partial disruptions with fractional mass-loss $\gtrsim 10\%$ should occur with a rate similar to that of total disruptions; for fractional mass-loss $\gtrsim 50\%$, the rate is $\approx 1/3$ as large. Partial disruptions---which must precede total disruptions when the stars' angular momenta evolve in the "empty loss-cone" regime---change the orbital energy by factors $\gtrsim O(1)$. Remnants of partial disruptions are in general far from thermal equilibrium. Depending on the orbital energy of the remnant and conditions within the stellar cluster surrounding the SMBH, it may return after hundreds or thousands of years and be fully disrupted, or it may rejoin the stellar cluster.

In this paper, the third in this series, we continue our study of tidal disruption events of main-sequence stars by a non-spinning $10^{6}~\rm{M}_\odot$ supermassive black hole. Here we focus on the outcomes of partial disruptions. As the encounter becomes weaker, the debris mass is increasingly concentrated near the outer edges of the energy distribution. As a result, the mass fallback rate can deviate substantially from a $t^{-5/3}$ power-law, becoming more like a single peak with a tail declining as $t^{-p}$ with $p\simeq2-5$. Surviving remnants are spun-up in the prograde direction and are hotter than MS stars of the same mass. Their specific orbital energy is $\simeq10^{-3}\times$ that of the debris (but of either sign with respect to the black hole potential) while their specific angular momentum is close to that of the original star. Even for strong encounters, remnants have speeds at infinity relative to the black hole potential $\lesssim 300$ km s$^{-1}$, so they are unable to travel far out into the galactic bulge. Remnants bound to the black hole can possibly go through a second tidal disruption event.

A star is tidally disrupted by a supermassive black hole when their separation is shorter than the "tidal radius". This quantity is often estimated on an order-of-magnitude basis without reference to the star's internal structure. Using MESA models for main sequence stars and fully general relativistic dynamics, we find the physical tidal radius for complete disruption $\cal{R}_t$ for a $10^6M_\odot$ black hole (BH). We find that across a factor $\sim20$ in stellar mass $M_*$, i.e., $0.15M_{\odot}\leq M_*\leq3M_\odot$, $\cal{R}_t\sim27\times$(BH's gravitational radius). When comparing $\cal{R}_t$ with the commonly used order-of-magnitude estimate $r_t$, we find that $\cal{R}_t\sim1.05-1.45r_t$ for $0.15M_\odot\leq M_*\leq0.5M_\odot$, but between $0.5 M_\odot$ and $1 M_\odot$, $\cal{R}_t$ drops to $\sim 0.45r_t$, and it remains at this value up to $10 M_\odot$. The near-constancy of $\cal{R}_t$ implies a weaker dependence of the full disruption rate on $M_*$ than when predicted with $r_t$. The characteristic energy width of the debris $\Delta E$ ranges from $\sim1.2\Delta\cal{E}$ for low-mass stars to $\sim 0.35\Delta\cal{E}$ for higher-mass stars, where $\Delta\cal{E}=GM_{\rm BH}R_*/\cal{R}_t^{2}$. We present analytic fits for the $M_*$ dependence of $\cal{R}_t$ and $\Delta E$; these fits lead to analytic expressions for the time of peak mass fallback rate and the maximal mass fallback rate. Our results also bear on the fraction of events leading to fast or slow circularization, as well as on the character of the tidal event occurring when the remnant of a partial disruption returns to the black hole. Using a semi-analytic model, we show that $\cal{R}_t$ is primarily determined by the star's central density rather than its mean density. For high-mass stars, the full disruption rate is roughly 1/4 the partial disruption rate, while this ratio is close to unity for low-mass stars.

This is the second in a series of papers presenting the results of fully general relativistic simulations of stellar tidal disruptions in which the stars' initial states are realistic main-sequence models. We consider eight different stellar masses, from $0.15~{\rm M}_\odot$ to $10~{\rm M}_\odot$. In the first paper (Ryu et al. 2019a), we gave an overview of this program and discussed the principal observational implications of our work. Here we describe our calculational method and provide details about the outcomes of full disruptions. We find that, relative to the traditional order-of-magnitude estimate $r_{\rm t}$, the physical tidal radius of low-mass stars is larger by tens of percent, while for high-mass stars ($M_{\star} \gtrsim1~ {\rm M}_\odot$) it is smaller by a factor $2-2.5$. The traditional estimate of the range of energies found in the debris is approximately accurate for low-mass stars, but is a factor $\sim 2$ too small for high-mass stars; in addition, the energy distribution for high-mass stars has significant wings. For all stars undergoing tidal encounters, we find that mass-loss continues for a long time because the ${\it instantaneous}$ tidal radius, the distance out to which the black hole's tidal gravity competes with the instantaneous stellar gravity at the star's surface, stays comparable to the distance to the black hole until the star has reached $O(10)~r_{\rm t}$. These findings indicate significant failings in the popular "frozen-in" approximation.

When a star gets too close to a supermassive black hole, it is torn apart by the tidal forces. Roughly half of the stellar mass becomes unbound and flies away at tremendous velocities - around $10^4$ km/s. In this work we explore the idea that the shock produced by the interaction of the unbound debris with the ambient medium gives rise to the synchrotron radio emission observed in several TDEs. We use a moving mesh numerical simulation to study the evolution of the unbound debris and the bow shock around it. We find that as the periapse distance of the star decreases, the outflow becomes faster and wider. A tidal disruption event whose periapse distance is a factor of 7 smaller than the tidal radius can account for the radio emission observed in ASASSN-14li. This model also allows us to obtain a more accurate estimate for the gas density around the centre of the host galaxy of ASASSN-14li.

Distinguishing the photon output of an accreting supermassive black hole binary system from that of a single supermassive black hole accreting at the same rate is intrinsically difficult because the majority of the light emerges from near the ISCOs of the black holes. However, there are two possible signals that can distinctively mark binaries, both arising from the gap formed in circumbinary accretion flows inside approximately twice the binary separation. One of these is a "notch" cut into the thermal spectra of these systems in the IR/optical/UV, the other a periodically-varying excess hard X-ray luminosity whose period is of order the binary orbital period. Using data from detailed galaxy evolution simulations, we estimate the distribution function in mass, mass ratio, and accretion rate for accreting supermassive black hole binaries as a function of redshift and then transform this distribution function into predicted source counts for these two potential signals. At flux levels >~10^{-13}~erg~cm^{-2}~s^{-1}, there may be ~O(10^2) such systems in the sky, mostly in the redshift range ~0.5 < z < ~1. Roughly 10% should have periods short enough (<~5~yr) to detect the X-ray modulation; this is also the period range accessible to PTA observations.

In this paper we continue the first ever study of magnetized mini-disks coupled to circumbinary accretion in a supermassive binary black hole (SMBBH) approaching merger reported in Bowen et al. 2018. We extend this simulation from 3 to 12 binary orbital periods. We find that relativistic SMBBH accretion acts as a resonant cavity, where quasi-periodic oscillations tied to the the frequency at which the black hole's orbital phase matches a non-linear $m=1$ density feature, or ``lump'', in the circumbinary accretion disk permeate the system. The rate of mass accretion onto each of the mini-disks around the black holes is modulated at the beat frequency between the binary frequency and the lump's mean orbital frequency, i.e., $\Omega_{\rm beat} = \Omega_{\rm bin} - \bar{\Omega}_{\rm lump}$, while the total mass accretion rate of this equal-mass binary is modulated at two different frequencies, $\gtrsim \bar{\Omega}_{\rm lump}$ and $\approx 2 \Omega_{\rm beat}$. The instantaneous rotation rate of the lump itself is also modulated at two frequencies close to the modulation frequencies of the total accretion rate, $\bar{\Omega}_{\rm lump}$ and $2 \Omega_{\rm beat}$. Because of the compact nature of the mini-disks in SMBBHs approaching merger, the inflow times within the mini-disks are comparable to the period on which their mass-supply varies, so that their masses---and the accretion rates they supply to their black holes---are strongly modulated at the same frequency. In essence, the azimuthal symmetry of the circumbinary disk is broken by the dynamics of orbits near a binary, and this $m=1$ asymmetry then drives quasi-periodic variation throughout the system, including both accretion and disk-feeding. In SMBBHs approaching merger, such time variability could introduce distinctive, increasingly rapid, fluctuations in their electromagnetic emission.

The grand challenges of contemporary fundamental physics---dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem---all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horizons, singularities and ergoregions. The hitherto invisible landscape of the gravitational Universe is being unveiled before our eyes: the historical direct detection of gravitational waves by the LIGO-Virgo collaboration marks the dawn of a new era of scientific exploration. Gravitational-wave astronomy will allow us to test models of black hole formation, growth and evolution, as well as models of gravitational-wave generation and propagation. It will provide evidence for event horizons and ergoregions, test the theory of General Relativity itself, and may reveal the existence of new fundamental fields. The synthesis of these results has the potential to radically reshape our understanding of the cosmos and of the laws of Nature. The purpose of this work is to present a concise, yet comprehensive overview of the state of the art in the relevant fields of research, summarize important open problems, and lay out a roadmap for future progress.

We describe results from a new technique for the prediction of complete, self-consistent X-ray spectra from three-dimensional General Relativistic magnetohydrodynamic (GRMHD) simulations of black hole accretion flows. Density and cooling rate data from a HARM3D GRMHD simulation are processed by both an improved version of the Monte Carlo radiation transport code PANDURATA (in the corona) and the Feautrier solver PTRANSX (in the disk), with XSTAR subroutines. The codes are run in a sequential but iterative fashion to achieve globally energy-conserving and self-consistent radiation fields, temperature maps, and photoionization equilibria. The output is the X-ray spectrum as seen by a distant observer. For the example cases we consider here---a non-rotating $10 M_\odot$ black hole with solar abundances, accreting at 0.01, 0.03, 0.1, or 0.3 Eddington---we find spectra resembling actual observations of stellar-mass black holes in the soft or steep power-law state: broad thermal peaks (at 1-3 keV), steep power-laws extending to high energy ($\Gamma$ = 2.7-4.5), and prominent, asymmetric Fe K$\alpha$ emission lines with equivalent widths in the range 40-400 eV (larger EW at lower accretion rates). By starting with simulation data, we obviate the need for parameterized descriptions of the accretion flow geometry---no a priori specification of the corona's shape or flux, or the disk temperature or density, etc., are needed. Instead, we apply the relevant physical principles to simulation output using appropriate numerical techniques; this procedure allows us to calculate inclination-dependent spectra after choosing only a small number of physically meaningful parameters: black hole mass and spin, accretion rate, and elemental abundances.

We present a series of simulations in both pure hydrodynamics (HD) and magnetohydrodynamics (MHD) exploring the degree to which alignment of disks subjected to external precessional torques (e.g., as in the `Bardeen-Petterson' effect) is dependent upon the disk sound speed c_s. Across the range of sound speeds examined, we find that the influence of the sound speed can be encapsulated in a simple "lumped-parameter" model proposed by Sorathia et al. (2013a). In this model, alignment fronts propagate outward at a speed ~0.2 rOmega_precess(r), where Omega_precess is the local test-particle precession frequency. Meanwhile, transonic radial motions transport angular momentum both inward and outward at a rate that may be described roughly in terms of an orientation diffusion model with diffusion coefficient ~2c_s^2/Omega, for local orbital frequency Omega. The competition between the two leads, in isothermal disks, to a stationary position for the alignment front at a radius proportional to c_s^(-4/5). For alignment to happen at all, the disk must either be turbulent due to the magnetorotational instability in MHD, or, in HD, it must be cool enough for the bending waves driven by disk warp to be nonlinear at their launch point. Contrary to long-standing predictions, warp propagation in MHD disks is diffusive independent of the parameter c_s/(alpha v_orb$, for orbital speed v_orb and ratio of stress to pressure of alpha. In purely HD disks, i.e., those with no internal stresses other than bulk viscosity, warmer disks align weakly or not at all; cooler disks align qualitatively similarly to MHD disks.

Spectropolarimetry is a powerful technique that has provided critical support for the geometric unification model of local active galactic nuclei. In this paper, we present optical (rest-frame UV) Keck spectropolarimetry of five luminous obscured (Type 2) and extremely red quasars (ERQs) at z~2.5. Three objects reach polarization fractions of >10% in the continuum. We propose a model in which dust scattering is the dominant scattering and polarization mechanism in our targets, though electron scattering cannot be completely excluded. Emission lines are polarized at a lower level than is the continuum. This suggests that the emission-line region exists on similar spatial scales as the scattering region. In three objects we detect an intriguing 90 degree swing in the polarization position angle as a function of line-of-sight velocity in the emission lines of Ly-alpha, CIV and NV. We interpret this phenomenon in the framework of a geometric model with an equatorial dusty scattering region in which the material is outflowing at several thousand km/sec. Emission lines may also be scattered by dust or resonantly. This model explains several salient features of observations by scattering on scales of a few tens of pc. Our observations provide a tantalizing view of the inner region geometry and kinematics of high-redshift obscured and extremely red quasars. Our data and modeling lend strong support for toroidal obscuration and powerful outflows on the scales of the UV emission-line region, in addition to the larger scale outflows inferred previously from the optical emission-line kinematics.

We present the first fully relativistic prediction of the electromagnetic emission from the surrounding gas of a supermassive binary black hole system approaching merger. Using a ray-tracing code to post-process data from a general relativistic 3-d MHD simulation, we generate images and spectra, and analyze the viewing angle dependence of the light emitted. When the accretion rate is relatively high, the circumbinary disk, accretion streams, and mini-disks combine to emit light in the UV/EUV bands. We posit a thermal Compton hard X-ray spectrum for coronal emission; at high accretion rates, it is almost entirely produced in the mini-disks, but at lower accretion rates it is the primary radiation mechanism in the mini-disks and accretion streams as well. Due to relativistic beaming and gravitational lensing, the angular distribution of the power radiated is strongly anisotropic, especially near the equatorial plane.

We present a "multipatch" infrastructure for numerical simulation of fluid problems in which sub-regions require different gridscales, different grid geometries, different physical equations, or different reference frames. Its key element is a sophisticated client-router-server framework for efficiently linking processors supporting different regions ("patches") that must exchange boundary data. This infrastructure may be used with a wide variety of fluid dynamics codes; the only requirement is that their primary dependent variables be the same in all patches, e.g., fluid mass density, internal energy density, and velocity. Its structure can accommodate either Newtonian or relativistic dynamics. The overhead imposed by this system is both problem- and computer cluster architecture-dependent. Compared to a conventional simulation using the same number of cells and processors, the increase in runtime can be anywhere from negligible to a factor of a few; however, one of the infrastructure's advantages is that it can lead to a very large reduction in the total number of zone-updates.

Eccentric disks arise in such astrophysical contexts as tidal disruption events, but it is unknown whether the magnetorotational instability (MRI), which powers accretion in circular disks, operates in eccentric disks as well. We examine the linear evolution of unstratified, incompressible MRI in an eccentric disk orbiting a point mass. We consider vertical modes of wavenumber $k$ on a background flow with uniform eccentricity $e$ and vertical Alfv\'en speed $v_\mathrm A$ along an orbit with mean motion $n$. We find two mode families, one with dominant magnetic components, the other with dominant velocity components; the former is unstable at $(1-e)^3f^2\lesssim3$, where $f\equiv kv_\mathrm A/n$, the latter at $e\gtrsim0.8$. For $f^2\lesssim3$, MRI behaves much like in circular disks, but the growth per orbit declines slowly with increasing $e$; for $f^2\gtrsim3$, modes grow by parametric amplification, which is resonant for $0<e\ll1$. MRI growth and the attendant angular momentum and energy transport happen chiefly near pericenter, where orbital shear dominates magnetic tension.

We present the first magnetohydrodynamic simulation in which a circumbinary disk around a relativistic binary black hole feeds mass to individual accretion disks ("mini-disks") around each black hole. Mass flow through the accretion streams linking the circumbinary disk to the mini-disks is modulated quasi-periodically by the streams' interaction with a nonlinear $m=1$ density feature, or "lump", at the inner edge of the circumbinary disk: the stream supplying each mini-disk comes into phase with the lump at a frequency $0.74$ times the binary orbital frequency. Because the binary is relativistic, the tidal truncation radii of the mini-disks are not much larger than their innermost stable circular orbits; consequently, the mini-disks' inflow times are shorter than the conventional estimate and are comparable to the stream modulation period. As a result, the mini-disks are always in inflow disequilibrium, with their masses and spiral density wave structures responding to the stream's quasi-periodic modulation. The fluctuations in each mini-disk's mass are so large that as much as $75\%$ of the total mini-disk mass can be contained within a single mini-disk. Such quasi-periodic modulation of the mini-disk structure may introduce distinctive time-dependent features in the binary's electromagnetic emission.

Near-Eddington radiation from active galactic nuclei (AGNs) has significant dynamical influence on the surrounding dusty gas, plausibly furnishing AGNs with geometrically thick obscuration. We investigate this paradigm with radiative magnetohydrodynamics simulations. The simulations solve the magnetohydrodynamics equations simultaneously with the infrared (IR) and ultraviolet (UV) radiative transfer (RT) equations; no approximate closure is used for RT. We find that our torus, when given a suitable sub-Keplerian angular momentum profile, spontaneously evolves toward a state in which its opening angle, density distribution, and flow pattern change only slowly. This "steady" state lasts for as long as there is gas resupply toward the inner edge. The torus is best described as a mid-plane inflow and a high-latitude outflow. The outflow is launched from the torus inner edge by UV radiation and expands in solid angle as it ascends; IR radiation continues to drive the wide-angle outflow outside the central hole. The dusty outflow obscures the central source in soft X-rays, the IR, and the UV over three quarters of solid angle, and each decade in column density covers roughly equal solid angle around the central source; these obscuration properties are similar to what observations imply.

Reverberation observations have uncovered an Fe K{\alpha} fluorescence line in the tidal disruption event (TDE) Swift J1644+57 (Kara et al. 2016). The discovery paper used the lag spectrum to argue that the X-ray continuum source was located very close to the blackhole (~30 gravitational radii) and moved sub-relativistically. We reanalyze the lag spectrum, pointing out that dilution effects cause it to indicate a geometric scale an order of magnitude larger than previously inferred. If the X-ray continuum is produced by a relativistic jet, as suggested by rapid variability, high luminosity and hard spectrum, this larger scale predicts an Fe ionization state consistent with efficient K{\alpha} production. Moreover, the momentum of the jet radiation impinging on the surrounding accretion flow on this larger scale accelerates a layer of gas to speeds ~0.1-0.2c, consistent with the blueshifted line profile. Implications of our results on the global picture of jetted TDEs are discussed. A power-law {\gamma}/X-ray spectrum may be produced by external UV-optical photons being repetitively inverse-Compton scattered by cold electrons in the jet, although our model for the K{\alpha} reverberation does not depend on the jet radiation mechanism (magnetic reconnection in a Poynting jet is still a viable mechanism). The non-relativistic wind driven by jet radiation may explain the late-time radio rebrightening in Swift J1644+57. This energy injection may also cause the thermal UV-optical emission from jetted TDEs to be systematically brighter than in non-jetted ones.

We present the first exploration of relativistic gas dynamics in the immediate vicinity of binary black holes as the system inspirals close to merger in the gravitational radiation-driven regime. We focus on 2D hydrodynamical studies of comparable-mass, non-spinning systems. Relativistic effects alter the dynamics of gas in this environment in several ways. Because the gravitational potential between the two black holes becomes shallower than in the Newtonian regime, the mini-disks stretch toward the L1 point and the amount of gas passing back and forth between the mini-disks increases sharply with decreasing binary separation. This "sloshing" is quasi-periodically modulated at $2$ and $2.75$ times the binary orbital frequency, corresponding to timescales of hours to days for supermassive binary black holes. In addition, relativistic effects add an $m=1$ component to the tidally-driven spiral waves in the disks that are purely $m=2$ in Newtonian gravity; this component becomes dominant when the separation is $\lesssim 100$ gravitational radii. Both the sloshing and the spiral waves have the potential to create distinctive radiation features that may uniquely mark supermassive binary black holes in the relativistic regime.

Models for tidal disruption events (TDEs) in which a supermassive black hole disrupts a star commonly assume that the highly eccentric streams of bound stellar debris promptly form a circular accretion disk at the pericenter scale. However, recent numerical simulations (Shiokawa et al., 2015) demonstrated that dissipation via hydrodynamical shocks is insufficient to circularize debris, and the flow retains its initial semi-major axis scale throughout the first ~10 orbits of the event. The bolometric peak luminosity of most TDE candidates, a few x 10^{44} erg/s, implies that we observe only ~1% of the energy expected from radiatively efficient accretion. Motivated by these results, (Piran et al., 2015) suggested that the observed optical TDE emission is powered by shocks at the apocenter between freshly infalling material and earlier arriving matter. This model explains the small radiated energy, the low temperature, and the large radius implied by the observations as well as the t^{-5/3} light curve. However the question of the system's low {bolometric} efficiency remains unanswered. We suggest that the high orbital energy and low angular momentum of the flow make it possible for magnetic stresses to reduce the matter's already small angular momentum to the point at which it can fall ballistically into the SMBH before circularization. As a result, the efficiency is only ~1--10% of a standard accretion disk's efficiency. Thus, the intrinsically high eccentricity of the tidal debris naturally explains why most TDE candidates are fainter than expected.

When a circumbinary disk surrounds a binary whose secondary's mass is at least $\sim 10^{-2}\times$ the primary's mass, a nearly empty cavity with radius a few times the binary separation is carved out of the disk. Narrow streams of material pass from the inner edge of the circumbinary disk into the domain of the binary itself, where they eventually join onto the small disks orbiting the members of the binary. Using data from 3-d MHD simulations of this process, we determine the luminosity of these streams; it is mostly due to weak laminar shocks, and is in general only a few percent of the luminosity of adjacent regions of either the circumbinary disk or the "mini-disks". This luminosity therefore hardly affects the deficit in the thermal continuum predicted on the basis of a perfectly dark gap region.

Stars that pass within the Roche radius of a supermassive black hole will be tidally disrupted, yielding a sudden injection of gas close to the black hole horizon which produces an electromagnetic flare. A few dozen of these flares have been discovered in recent years, but current observations provide poor constraints on the bolometric luminosity and total accreted mass of these events. Using images from the Wide-field Infrared Survey Explorer (WISE), we have discovered transient 3.4 micron emission from several previously known tidal disruption flares. The observations can be explained by dust heated to its sublimation temperature due to the intense radiation of the tidal flare. From the break in the infrared light curve we infer that this hot dust is located ~0.1 pc from the supermassive black hole. Since the dust has been heated by absorbing UV and (potentially) soft X-ray photons of the flare, the reprocessing light curve yields an estimate of the bolometric flare luminosity. For the flare PTF-09ge, we infer that the most likely value of the luminosity integrated over frequencies at which dust can absorb photons is $8\times 10^{44}$ erg/s, with a factor of 3 uncertainty due to the unknown temperature of the dust. This bolometric luminosity is a factor ~10 larger than the observed black body luminosity. Our work is the first to probe dust in the nuclei of non-active galaxies on sub-parsec scales. The observed infrared luminosity implies a covering factor ~1% for the nuclear dust in the host galaxies.

Substantial evidence points to dusty, geometrically thick tori obscuring the central engines of active galactic nuclei (AGNs), but so far no mechanism satisfactorily explains why cool dust in the torus remains in a puffy geometry. Near-Eddington infrared (IR) and ultraviolet (UV) luminosities coupled with high dust opacities at these frequencies suggest that radiation pressure on dust can play a significant role in shaping the torus. To explore the possible effects of radiation pressure, we perform three-dimensional radiative hydrodynamics simulations of an initially smooth torus. Our code solves the hydrodynamics equations, the time-dependent multi-angle group IR radiative transfer (RT) equation, and the time-independent UV RT equation. We find a highly dynamic situation. IR radiation is anisotropic, leaving primarily through the central hole. The torus inner surface exhibits a break in axisymmetry under the influence of radiation and differential rotation; clumping follows. In addition, UV radiation pressure on dust launches a strong wind along the inner surface; when scaled to realistic AGN parameters, this outflow travels at $\sim 5000 (M/10^7 M_\odot)^{1/4} [L_\mathrm{UV}/(0.1 L_\mathrm E)]^{1/4} \mathrm{km}\,\mathrm s^{-1}$ and carries $\sim 0.1 (M/10^7 M_\odot)^{3/4} [L_\mathrm{UV}/(0.1 L_\mathrm E)]^{3/4} M_\odot\,\mathrm{yr}^{-1}$, where $M$, $L_\mathrm{UV}$, and $L_\mathrm E$ are the mass, UV luminosity, and Eddington luminosity of the central object respectively.

ASASSN-14li is a recently-discovered tidal disruption event with an exceptionally rich data-set: spectra and lightcurves in soft X-rays, UV, optical, and radio. To understand its emission properties in all these bands, we have extended our model for post-tidal disruption accretion and photon production to estimate both soft X-ray radiation produced by the prompt accretion phase and synchrotron emission associated with the bow shock driven through an external medium by the unbound tidal debris, as well as optical and UV light. We find that fiducial values of the stellar mass ($1 M_\odot$) and black hole mass ($10^{6.5} M_{\odot}$) yield: quantitative agreement with the optical/UV luminosity, lightcurve, and color temperature; approximate agreement with the somewhat uncertain soft X-ray spectrum and lightcurve; and quantitative agreement with the radio luminosity, spectrum and lightcurve. Equipartition analysis of the radio data implies that the radio-emitting region expands with a constant speed, and its magnitude is comparable to the speed expected for the unbound stellar ejecta. Both facts provide strong support to our model. We find that the disruption event took place in mid-September 2014. Two independent parameters, the magnitude and logarithmic radial gradient of the ambient gas density near the black hole, must be fit to the data to explain the radio emission; their inferred values are comparable to those found near both Sgr A* and the TDE candidate Swift J1644.

We present first results from a new technique for the prediction of Fe K$\alpha$ profiles directly from general relativistic magnetohydrodynamic (GRMHD) simulations. Data from a GRMHD simulation are processed by a Monte Carlo global radiation transport code, which determines the X-ray flux irradiating the disk surface and the coronal electron temperature self-consistently. With that irradiating flux and the disk's density structure drawn from the simulation, we determine the reprocessed Fe K$\alpha$ emission from photoionization equilibrium and solution of the radiation transfer equation. We produce maps of the surface brightness of Fe K$\alpha$ emission over the disk surface, which---for our example of a $10 M_\odot$, Schwarzschild black hole accreting at $1\%$ the Eddington value---rises steeply one gravitational radius outside the radius of the innermost stable circular orbit and then falls $\propto r^{-2}$ at larger radii. We explain these features of the Fe K$\alpha$ radial surface brightness profile as consequences of the disk's ionization structure and an extended coronal geometry, respectively. We also present the corresponding Fe K$\alpha$ line profiles as would be seen by distant observers at several inclinations. Both the shapes of the line profiles and the equivalent widths of our predicted K$\alpha$ lines are qualitatively similar to those typically observed from accreting black holes. Most importantly, this work represents a direct link between theory and observation: in a fully self-consistent way, we produce observable results---iron fluorescence line profiles---from the theory of black hole accretion with almost no phenomenological assumptions.

We present the results of a new series of global 3D relativistic magneto-hydrodynamic (MHD) simulations of thin accretion disks around spinning black holes. The disks have aspect ratios of $H/R\sim 0.05$ and spin parameters $a/M=0, 0.5, 0.9$, and $0.99$. Using the ray-tracing code Pandurata, we generate broad-band thermal spectra and polarization signatures from the MHD simulations. We find that the simulated spectra can be well fit with a simple, universal emissivity profile that better reproduces the behavior of the emission from the inner disk, compared to traditional analyses carried out using a Novikov-Thorne thin disk model. Lastly, we show how spectropolarization observations can be used to convincingly break the spin-inclination degeneracy well-known to the continuum fitting method of measuring black hole spin.

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

When an accretion disk surrounds a binary rotating in the same sense, the binary exerts strong torques on the gas. Analytic work in the 1D approximation indicated that these torques sharply diminish or even eliminate accretion from the disk onto the binary. However, recent 2D and 3D simulational work has shown at most modest diminution. We present new MHD simulations demonstrating that for binaries with mass ratios of 1 and 0.1 there is essentially no difference between the accretion rate at large radius in the disk and the accretion rate onto the binary. To resolve the discrepancy with earlier analytic estimates, we identify the small subset of gas trajectories traveling from the inner edge of the disk to the binary and show how the full accretion rate is concentrated onto them.

Using only physical mechanisms, i.e., 3D MHD with no phenomenological viscosity, we have simulated the dynamics of a moderately thin accretion disk subject to torques whose radial scaling mimics those produced by lowest-order post-Newtonian gravitomagnetism. In this simulation, we have shown how, in the presence of MHD turbulence, a time-steady transition can be achieved between an inner disk region aligned with the equatorial plane of the central mass's spin and an outer region orbiting in a different plane. The position of the equilibrium orientation transition is determined by a balance between gravitomagnetic torque and warp-induced inward mixing of misaligned angular momentum from the outer disk. If the mixing is interpreted in terms of diffusive transport, the implied diffusion coefficient is ~(0.6--0.8)c_s^2/Omega for sound speed c_s and orbital frequency Omega. This calibration permits estimation of the orientation transition's equilibrium location given the central mass, its spin parameter, and the disk's surface density and scaleheight profiles. However, the alignment front overshoots before settling into an equilibrium, signaling that a diffusive model does not fully represent the time-dependent properties of alignment fronts under these conditions. Because the precessional torque on the disk at the alignment front is always comparable to the rate at which misaligned angular momentum is brought inward to the front by warp-driven radial motions, no break forms between the inner and outer portions of the disk in our simulation. Our results also raise questions about the applicability to MHD warped disks of the traditional distinction between "bending wave" and "diffusive" regimes.

A tidal disruption event (TDE) takes place when a star passes near enough to a massive black hole to be disrupted. About half the star's matter is given elliptical trajectories with large apocenter distances, the other half is unbound. To "circularize", i.e., to form an accretion flow, the bound matter must lose a significant amount of energy, with the actual amount depending on the characteristic scale of the flow measured in units of the black hole's gravitational radius ($\sim 10^{51} (R/1000R_g)^{-1}$~erg). Recent numerical simulations \citep{Shiokawa+2015} have revealed that the circularization scale is close to the scale of the most-bound initial orbits, $\sim 10^3 M_{BH,6.5}^{-2/3} R_g \sim 10^{15} M_{BH,6.5}^{1/3}$~cm from the black hole, and the corresponding circularization energy dissipation rate is $\sim 10^{44} M_{BH,6.5}^{-1/6}$~erg/s. We suggest that the energy liberated during circularization, rather then energy liberated by accretion onto the black hole, powers the observed optical TDE candidates. The observed rise times, luminosities, temperatures, emission radii, and line widths seen in these TDEs \citep[e.g.][]{Arcavi+2014} are all more readily explained in terms of heating associated with circularization than in terms of accretion.

We consider the evolution of a supermassive black hole binary (SMBHB) surrounded by a retrograde accretion disk. Assuming the disk is exactly in the binary plane and transfers energy and angular momentum to the binary via direct gas accretion, we calculate the time evolution of the binary's semi-major axis $a$ and eccentricity $e$. Because the gas is predominantly transferred when the binary is at apocenter, we find the eccentricity grows rapidly while maintaining constant $a(1+e)$. After accreting only a fraction of the secondary's mass, the eccentricity grows to nearly unity; from then on, gravitational wave emission dominates the evolution, preserving constant $a(1-e)$. The high-eccentricity waveforms redistribute the peak gravitational wave power from the nHz to $\mu$Hz bands, substantially affecting the signal that might be detected with pulsar timing arrays. We also estimate the torque coupling binaries of arbitrary eccentricity with obliquely aligned circumbinary disks. If the outer edge of the disk is not an extremely large multiple of the binary separation, retrograde accretion can drive the binary into the gravitational wave-dominated state before these torques align the binary with the angular momentum of the mass supply.

We study how the matter dispersed when a supermassive black hole tidally disrupts a star joins an accretion flow. Combining a relativistic hydrodynamic simulation of the stellar disruption with a relativistic hydrodynamics simulation of the tidal debris motion, we track such a system until ~80% of the stellar mass bound to the black hole has settled into an accretion flow. Shocks near the stellar pericenter and also near the apocenter of the most tightly-bound debris dissipate orbital energy, but only enough to make the characteristic radius comparable to the semi-major axis of the most-bound material, not the tidal radius as previously thought. The outer shocks are caused by post-Newtonian effects, both on the stellar orbit during its disruption and on the tidal forces. Accumulation of mass into the accretion flow is non-monotonic and slow, requiring ~3--10x the orbital period of the most tightly-bound tidal streams, while the inflow time for most of the mass may be comparable to or longer than the mass accumulation time. Deflection by shocks does, however, remove enough angular momentum and energy from some mass for it to move inward even before most of the mass is accumulated into the accretion flow. Although the accretion rate rises sharply and then decays roughly as a power-law, its maximum is ~0.1x the previous expectation, and the duration of the peak is ~5x longer than previously predicted. The geometric mean of the black hole mass and stellar mass inferred from a measured event timescale is therefore ~0.2x the value given by classical theory.

We present the results of local, vertically stratified, radiation MHD shearing box simulations of MRI turbulence appropriate for the hydrogen ionizing regime of dwarf nova and soft X-ray transient outbursts. We incorporate the frequency-integrated opacities and equation of state for this regime, but neglect non-ideal MHD effects and surface irradiation, and do not impose net vertical magnetic flux. We find two stable thermal equilibrium tracks in the effective temperature versus surface mass density plane, in qualitative agreement with the S-curve picture of the standard disk instability model. We find that the large opacity at temperatures near $10^4$K, a corollary of the hydrogen ionization transition, triggers strong, intermittent thermal convection on the upper stable branch. This convection strengthens the magnetic turbulent dynamo and greatly enhances the time-averaged value of the stress to thermal pressure ratio $\alpha$, possibly by generating vertical magnetic field that may seed the axisymmetric magnetorotational instability, and by increasing cooling so that the pressure does not rise in proportion to the turbulent dissipation. These enhanced stress to pressure ratios may alleviate the order of magnitude discrepancy between the $\alpha$-values observationally inferred in the outburst state and those that have been measured from previous local numerical simulations of magnetorotational turbulence that lack net vertical magnetic flux.

Observations indicate that most massive galaxies contain a supermassive black hole, and theoretical studies suggest that when such galaxies have a major merger, the central black holes will form a binary and eventually coalesce. Here we discuss two spectral signatures of such binaries that may help distinguish them from ordinary AGN. These signatures are expected when the mass ratio between the holes is not extreme and the system is fed by a circumbinary disk. One such signature is a notch in the thermal continuum that has been predicted by other authors; we point out that it should be accompanied by a spectral revival at shorter wavelengths and also discuss its dependence on binary properties such as mass, mass ratio, and separation. In particular, we note that the wavelength $\lambda_n$ at which the notch occurs depends on these three parameters in such a way as to make the number of systems displaying these notches $\propto \lambda_n^{16/3}$; longer wavelength searches are therefore strongly favored. A second signature, first discussed here, is hard X-ray emission with a Wien-like spectrum at a characteristic temperature $\sim 100$ keV produced by Compton cooling of the shock generated when streams from the circumbinary disk hit the accretion disks around the individual black holes. We investigate the observability of both signatures. The hard X-ray signal may be particularly valuable as it can provide an indicator of black hole merger a few decades in advance of the event.

When matter orbits around a central mass obliquely with respect to the mass's spin axis, the Lense-Thirring effect causes it to precess at a rate declining sharply with radius. Ever since the work of Bardeen & Petterson (1975), it has been expected that when a fluid fills an orbiting disk, the orbital angular momentum at small radii should then align with the mass's spin. Nearly all previous work has studied this alignment under the assumption that a phenomenological "viscosity" isotropically degrades fluid shears in accretion disks, even though it is now understood that internal stress in flat disks is due to anisotropic MHD turbulence. In this paper we report a pair of matched simulations, one in MHD and one in pure (non-viscous) HD in order to clarify the specific mechanisms of alignment. As in the previous work, we find that disk warps induce radial flows that mix angular momentum of different orientation; however, we also show that the speeds of these flows are generically transonic and are only very weakly influenced by internal stresses other than pressure. In particular, MHD turbulence does not act in a manner consistent with an isotropic viscosity. When MHD effects are present, the disk aligns, first at small radii and then at large; alignment is only partial in the HD case. We identify the specific angular momentum transport mechanisms causing alignment and show how MHD effects permit them to operate more efficiently. Lastly, we relate the speed at which an alignment front propagates outward (in the MHD case) to the rate at which Lense-Thirring torques deliver angular momentum at smaller radii.

Recent studies of accretion onto supermassive black hole binaries suggest that much, perhaps most, of the matter eventually accretes onto one hole or the other. If so, then for binaries whose inspiral from ~1 pc to 0.001 - 0.01 pc is driven by interaction with external gas, both the binary orbital axis and the individual black hole spins can be reoriented by angular momentum exchange with this gas. Here we show that, unless the binary mass ratio is far from unity, the spins of the individual holes align with the binary orbital axis in a time few-100 times shorter than the binary orbital axis aligns with the angular momentum direction of the incoming circumbinary gas; the spin of the secondary aligns more rapidly than that of the primary by a factor ~(m_1/m_2)^{1/2}>1. Thus the binary acts as a stabilizing agent, so that for gas-driven systems, the black hole spins are highly likely to be aligned (or counteraligned if retrograde accretion is common) with each other and with the binary orbital axis. This alignment can significantly reduce the recoil speed resulting from subsequent black hole merger.

Global disk simulations provide a powerful tool for investigating accretion and the underlying magnetohydrodynamic turbulence driven by the magneto-rotational instability (MRI). Using them to predict accurately quantities such as stress, accretion rate, and surface brightness profile requires that purely numerical effects, arising from both resolution and algorithm, be understood and controlled. We use the flux-conservative Athena code to conduct a series of experiments on disks having a variety of magnetic topologies to determine what constitutes adequate resolution. We develop and apply several resolution metrics: Qz and Qphi, the ratio of the grid zone size to the characteristic MRI wavelength, alpha_mag, the ratio of the Maxwell stress to the magnetic pressure, and the ratio of radial to toroidal magnetic field energy. For the initial conditions considered here, adequate resolution is characterized by Qz > 15, Qphi > 20, alpha_mag = 0.45, and a field energy ratio of 0.2. These values are associated with > 35 zones per scaleheight, a result consistent with shearing box simulations. Numerical algorithm is also important. Use of the HLLE flux solver or second-order interpolation can significantly degrade the effective resolution compared to the HLLD flux solver and third-order interpolation. Resolution at this standard can be achieved only with large numbers of grid zones, arranged in a fashion that matches the symmetries of the problem and the scientific goals of the simulation.

We present the results of a new global radiation transport code coupled to a general relativistic magneto-hydrodynamic simulation of an accreting, non-rotating black hole. For the first time, we are able to explain from first principles in a self-consistent way all the components seen in the X-ray spectra of stellar-mass black holes, including a thermal peak and all the features associated with strong hard X-ray emission: a power-law extending to high energies, a Compton reflection hump, and a broad iron line. Varying only the mass accretion rate, we are able to reproduce a wide range of X-ray states seen in most galactic black hole sources. The temperature in the corona is T_e ~ 10 keV in a boundary layer near the disk and rises smoothly to T_e >~ 100 keV in low-density regions far above the disk. Even as the disk's reflection edge varies from the horizon out to ~ 6M as the accretion rate decreases, we find that the shape of the Fe K\alpha line is remarkably constant. This is because photons emitted from the plunging region are strongly beamed into the horizon and never reach the observer. We have also carried out a basic timing analysis of the spectra and find that the fractional variability increases with photon energy and viewer inclination angle, consistent with the coronal hot spot model for X-ray fluctuations.

In this document, we describe the scientific potential of blazar observations with a X-ray polarimetry mission like GEMS (Gravity and Extreme Magnetism SMEX). We describe five blazar science investigations that such a mission would enable: (i) the structure and the role of magnetic fields in AGN jets, (ii) analysis of the polarization of the synchrotron X-ray emission from AGN jets, (iii) discrimination between synchrotron self-Compton and external Compton models for blazars with inverse Compton emission in the X-ray band, (iv) a precision study of the polarization properties of the X-ray emission from Cen-A, (v) tests of Lorentz Invariance based on X-ray polarimetric observations of blazars. We conclude with a discussion of a straw man observation program and recommended accompanying multiwavelength observations.

Orbiting disks may exhibit bends due to a misalignment between the angular momentum of the inner and outer regions of the disk. We begin a systematic simulational inquiry into the physics of warped disks with the simplest case: the relaxation of an unforced warp under pure fluid dynamics, i.e. with no internal stresses other than Reynolds stress. We focus on the nonlinear regime in which the bend rate is large compared to the disk aspect ratio. When warps are nonlinear, strong radial pressure gradients drive transonic radial motions along the disk's top and bottom surfaces that efficiently mix angular momentum. The resulting nonlinear decay rate of the warp increases with the warp rate and the warp width, but, at least in the parameter regime studied here, is independent of the sound speed. The characteristic magnitude of the associated angular momentum fluxes likewise increases with both the local warp rate and the radial range over which the warp extends; it also increases with increasing sound speed, but more slowly than linearly. The angular momentum fluxes respond to the warp rate after a delay that scales with the square-root of the time for sound waves to cross the radial extent of the warp. These behaviors are at variance with a number of the assumptions commonly used in analytic models to describe linear warp dynamics.

We present a new code for radiation transport around Kerr black holes, including arbitrary emission and absorption terms, as well as electron scattering and polarization. The code is particularly useful for analyzing accretion flows made up of optically thick disks and optically thin coronae. We give a detailed description of the methods employed in the code, and also present results from a number of numerical tests to assess its accuracy and convergence.

We examine the expected X-ray polarization properties of neutron-star X-ray sources of various types, e.g., accretion and rotation powered pulsars, magnetars, and low-mass X-ray binaries. We summarize the model calculations leading to these expected properties. We describe how a comparison of these with their observed properties, as inferred from GEMS data, will probe the essential dynamical, electromagnetic, plasma, and emission processes in neutron-star binaries, discriminate between models of these processes, and constrain model parameters. An exciting goal is the first observational demonstration in this context of the existence of vacuum resonance, a fundamental quantum electrodynamical phenomenon first described in the 1930s.

We present here a summary of the scientific goals behind the Gravity and Extreme Magnetism SMEX (GEMS) X-ray polarimetry mission's black hole (BH) observing program. The primary targets can be divided into two classes: stellar-mass galactic BHs in accreting binaries, and super-massive BHs in the centers of active galactic nuclei (AGN). The stellar-mass BHs can in turn be divided into various X-ray spectral states: thermal-dominant (disk), hard (radio jet), and steep power-law (hot corona). These different spectral states are thought to be generated by different accretion geometries and emission mechanisms. X-ray polarization is an ideal tool for probing the geometry around these BHs and revealing the specific properties of the accreting gas.

A hyperaccretion disk formed around a stellar mass black hole is a plausible model for the central engine that powers gamma-ray bursts (GRBs). If the central black hole rotates and a poloidal magnetic field threads its horizon, a powerful relativistic jet may be driven by a process resembling the Blandford-Znajek mechanism. We estimate the luminosity of such a jet assuming that the poloidal magnetic field strength is comparable to the inner accretion disk pressure. We show that the jet efficiency attains its maximal value when the accretion flow is cooled via optically-thin neutrino emission. The jet luminosity is much larger than the energy deposition through neutrino-antineutrino annihilation provided that the black hole is spinning rapidly enough. When the accretion rate onto a rapidly spinning black hole is large enough (> 0.003-0.01M_sun/sec), the predicted jet luminosity is sufficient to drive a GRB.

We explore the temporal structure of tidal disruption events pointing out the corresponding transitions in the lightcurves of the thermal accretion disk and of the jet emerging from such events. The hydrodynamic time scale of the disrupted star is the minimal time scale of building up the accretion disk and the jet and it sets a limit on the rise time. This suggest that Swift J1644+57, that shows several flares with a rise time as short as a few hundred seconds could not have arisen from a tidal disruption of a main sequence star whose hydrodynamic time is a few hours. The disrupted object must have been a white dwarf. A second important time scale is the Eddington time in which the accretion rate changes form super to sub Eddington. It is possible that such a transition was observed in the light curve of Swift J2058+05. If correct this provides intersting constraints on the parameters of the system.

In this white paper, we discuss the concept of a next-generation X-ray mission called BEST (Black hole Evolution and Space Time). The mission concept uses a 3000 square centimeter effective area mirror (at 6 keV) to achieve unprecedented sensitivities for hard X-ray imaging spectrometry (5-70 keV) and for broadband X-ray polarimetry (2-70 keV). BEST can make substantial contributions to our understanding of the inner workings of accreting black holes, our knowledge about the fabric of extremely curved spacetime, and the evolution of supermassive black holes. BEST will allow for time resolved studies of accretion disks. With a more than seven times larger mirror area and a seven times wider bandpass than GEMS, BEST will take X-ray polarimetry to a new level: it will probe the time variability of the X-ray polarization from stellar mass and supermassive black holes, and it will measure the polarization properties in 30 independent energy bins. These capabilities will allow BEST to conduct tests of accretion disk models and the underlying spacetimes. With three times larger mirror area and ten times better angular resolution than NuSTAR, BEST will be able to make deep field observations with a more than 15 times better sensitivity than NuSTAR. The mission will be able to trace the evolution of obscured and unobscured black holes in the redshift range from zero to six, covering the most important epoch of supermassive black hole growth. The hard X-ray sensitivity of BEST will enable a deep census of non-thermal particle populations. BEST will give us insights into AGN feedback by measuring the particle luminosity injected by AGNs into the interstellar medium (ISM) of their hosts, and will map the emission from particles accelerated at large scale structure shocks. Finally, BEST has the potential to constrain the equation of state of neutron stars (NS).

As 2 black holes bound to each other in a close binary approach merger their inspiral time becomes shorter than the characteristic inflow time of surrounding orbiting matter. Using an innovative technique in which we represent the changing spacetime in the region occupied by the orbiting matter with a 2.5PN approximation and the binary orbital evolution with 3.5PN, we have simulated the MHD evolution of a circumbinary disk surrounding an equal-mass non-spinning binary. Prior to the beginning of the inspiral, the structure of the circumbinary disk is predicted well by extrapolation from Newtonian results. The binary opens a low-density gap whose radius is roughly two binary separations, and matter piles up at the outer edge of this gap as inflow is retarded by torques exerted by the binary; nonetheless, the accretion rate is diminished relative to its value at larger radius by only about a factor of 2. During inspiral, the inner edge of the disk at first moves inward in coordination with the shrinking binary, but as the orbital evolution accelerates, the rate at which the inner edge moves toward smaller radii falls behind the rate of binary compression. In this stage, the rate of angular momentum transfer from the binary to the disk slows substantially, but the net accretion rate decreases by only 10-20%. When the binary separation is tens of gravitational radii, the rest-mass efficiency of disk radiation is a few percent, suggesting that supermassive binary black holes in galactic nuclei could be very luminous at this stage of their evolution. If the luminosity were optically thin, it would be modulated at a frequency that is a beat between the orbital frequency of the disk's surface density maximum and the binary orbital frequency. However, a disk with sufficient surface density to be luminous should also be optically thick; as a result, the periodic modulation may be suppressed.

We present the first three-dimensional magnetohydrodynamic (MHD) simulations of a circumbinary disk surrounding an equal mass binary. The binary maintains a fixed circular orbit of separation $a$. As in previous hydrodynamical simulations, strong torques by the binary can maintain a gap of radius $\simeq 2a$. Streams curve inward from $r \simeq 2a$ toward the binary; some of their mass passes through the inner boundary, while the remainder swings back out to the disk. However, we also find that near its inner edge the disk develops both a strong $m=1$ asymmetry and growing orbital eccentricity. Because the MHD stresses introduce more matter into the gap, the total torque per unit disk mass is $\simeq 14$ times larger than found previously. The inner boundary accretion rate per unit disk mass is $\simeq 40$ times greater than found from previous hydrodynamical calculations. The implied binary shrinkage rate is determined by a balance between the rate at which the binary gains angular momentum by accretion and loses it by gravitational torque. The large accretion rate brings these two rates nearly into balance, but in net, we find that $\dot a/a < 0$, and its magnitude is about 2.7 times larger than predicted by the earlier hydrodynamic simulations. If the binary comprises two massive black holes, the accretion rate may be great enough for one or both to be AGN, with consequences for the physical state of the gas both in the disk body and in its inner gap.

Tidal disruption of main sequence stars by black holes has generally been thought to lead to a signal dominated by UV emission. If, however, the black hole spins rapidly and the poloidal magnetic field intensity on the black hole horizon is comparable to the inner accretion disk pressure, a powerful jet may form whose luminosity can easily exceed the thermal UV luminosity. When the jet beam points at Earth, its non-thermal luminosity can dominate the emitted spectrum. The thermal and non-thermal components decay differently with time. In particular, the thermal emission should remain roughly constant for a significant time after the period of maximum accretion, beginning to diminish only after a delay, whereas after the peak accretion rate, the non-thermal jet emission decays, but then reaches a plateau. Both transitions are tied to a characteristic timescale $t_{\rm Edd}$ at which the accretion rate falls below Eddington. Making use of this timescale in a new parameter-inference formalism for tidal disruption events with significant emission from a jet, we analyze the recent flare source Swift J2058. It is consistent with an event in which a main sequence solar-type staris disrupted by a black hole of mass $\sim 4 \times 10^7 M_{\odot}$. The beginning of the flat phase in the non-thermal emission from this source can possibly be seen in the late-time lightcurve. Optical photometry over the first $\simeq 40$ d of this flare is also consistent with this picture, but is only weakly constraining because the bolometric correction is very uncertain. We suggest that future searches for main sequence tidal disruptions use methods sensitive to jet radiation as well as to thermal UV radiation.

We propose that the remarkable object Swift J1644+57, in which multiple recurring hard X-ray flares were seen over a span of several days, is a system in which a white dwarf was tidally disrupted by an intermediate mass black hole. Disruption of a white dwarf rather than a main sequence star offers a number of advantages in understanding the multiple, and short, timescales seen in the light curve of this system. In particular, the short internal dynamical timescale of a white dwarf offers a more natural way of understanding the short rise times ($\sim 100$ s) observed. The relatively long intervals between flares ($\sim 5 \times 10^4$ s) may also be readily understood as the period between successive pericenter passages of the remnant white dwarf. In addition, the expected jet power is larger when a white dwarf is disrupted. If this model is correct, the black hole responsible must have mass $\lesssim 10^5 M_{\odot}$.

Recent general relativistic magneto-hydrodynamic (MHD) simulations of accretion onto black holes have shown that, contrary to the basic assumptions of the Novikov-Thorne model, there can be substantial magnetic stress throughout the plunging region. Additional dissipation and radiation can therefore be expected. We use data from a particularly well-resolved simulation of accretion onto a non-spinning black hole to compute both the radiative efficiency of such a flow and its spectrum if all emitted light is radiated with a thermal spectrum whose temperature matches the local effective temperature. This disk is geometrically thin enough (H/r ~= 0.06) that little heat is retained in the flow. In terms of light reaching infinity (i.e., after allowance for all relativistic effects and for photon capture by the black hole), we find that the radiative efficiency is at least ~=6-10% greater than predicted by the Novikov-Thorne model (complete radiation of all heat might yield another ~6%). We also find that the spectrum more closely resembles the Novikov-Thorne prediction for a/M ~= 0.2--0.3 than for the correct value, a/M=0. As a result, if the spin of a non-spinning black hole is inferred by model-fitting to a Novikov-Thorne model with known black hole mass, distance, and inclination, the inferred a/M is too large by ~= 0.2--0.3.

Discretized numerical simulations are a powerful tool for investigation of nonlinear MHD turbulence in accretion disks. However, confidence in their quantitative predictions requires a demonstration that further refinement of the spatial gridscale would not result in any significant change. This has yet to be accomplished, particularly for global disk simulations. In this paper, we combine data from previously published stratified shearing box simulations and new global disk simulations to calibrate several quantitative diagnostics by which one can estimate progress toward numerical convergence. Using these diagnostics, we find that the established criterion for an adequate numerical description of linear growth of the magneto-rotational instability (the number of cells across a wavelength of the fastest-growing vertical wavenumber mode) can be extended to a criterion for adequate description of nonlinear MHD disk turbulence, but the standard required is more stringent. We also find that azimuthal resolution, which has not often been extensively examined in previous studies, can significantly affect the evolution of the poloidal magnetic field. We further analyze the comparative resolution requirements of a small sample of initial magnetic field geometries; not surprisingly, more complicated initial field geometries require higher spatial resolution. Otherwise, they tend to evolve to qualitatively similar states if evolved for sufficient time. Applying our quantitative resolution criteria to a sample of previously published global simulations, we find that, with perhaps a single exception, they are significantly under-resolved, and therefore underestimate the magnetic turbulence and resulting stress levels throughout the accretion flow.

Standard models of radiation supported accretion disks generally assume that diffusive radiation flux is solely responsible for vertical heat transport. This requires that heat must be generated at a critical rate per unit volume if the disk is to be in hydrostatic and thermal equilibrium. This raises the question of how heat is generated and how energy is transported in MHD turbulence. By analysis of a number of radiation/MHD stratified shearing-box simulations, we show that the divergence of the diffusive radiation flux is indeed capped at the critical rate, but deep inside the disk, substantial vertical energy flux is also carried by advection of radiation. Work done by radiation pressure is a significant part of the energy budget, and much of this work is dissipated later through damping by radiative diffusion. We show how this damping can be measured in the simulations, and identify its physical origins. Radiative damping accounts for as much as tens of percent of the total dissipation, and is the only realistic physical mechanism for dissipation of turbulence that can actually be resolved in numerical simulations of accretion disks. Buoyancy associated with dynamo-driven, highly magnetized, nearly-isobaric nonlinear slow magnetosonic fluctuations is responsible for the radiation advection flux, and also explains the persistent periodic magnetic upwelling seen at all values of the radiation to gas pressure ratio. The intimate connection between radiation advection and magnetic buoyancy is the first example we know of in astrophysics in which a dynamo has direct impact on the global energetics of a system.

Type 2 quasars are luminous active galactic nuclei (AGN) whose central regions are obscured by large amounts of gas and dust. In this paper, we present a catalog of type 2 quasars from the Sloan Digital Sky Survey (SDSS), selected based on their optical emission lines. The catalog contains 887 objects with redshifts z < 0.83; this is six times larger than the previous version and is by far the largest sample of type 2 quasars in the literature. We derive the [OIII]5008 luminosity function for 10^8.3 Lsun < L[OIII] < 10^10 Lsun (corresponding to intrinsic luminosities up to M[2400A]-28 mag or bolometric luminosities up to 4x10^47 erg/sec). This luminosity function provides strong lower limits to the actual space density of obscured quasars, due to our selection criteria, the details of the spectroscopic target selection, as well as other effects. We derive the equivalent luminosity function for the complete sample of type 1 (unobscured) quasars; then, we determine the ratio of type 2/type 1 quasar number densities. Our best data constrain this ratio to be at least 1.5:1 for 10^8.3 Lsun < L[OIII] < 10^9.5 Lsun at z < 0.3, and at least 1.2:1 for L[OIII]=10^10 Lsun at 0.3 < z < 0.83. Type 2 quasars are at least as abundant as type 1 quasars in the relatively nearby Universe (z < 0.8) for the highest luminosities.

We explore the parameter dependence of inner disk stress in black hole accretion by contrasting the results of a number of simulations, all employing 3-d general relativistic MHD in a Schwarzschild spacetime. Five of these simulations were performed with the intrinsically conservative code HARM3D, which allows careful regulation of the disk aspect ratio, H/R; our simulations span a range in H/R from 0.06 to 0.17. We contrast these simulations with two previously reported simulations in a Schwarzschild spacetime in order to investigate possible dependence of the inner disk stress on magnetic topology. In all cases, much care was devoted to technical issues: ensuring adequate resolution and azimuthal extent, and averaging only over those time-periods when the accretion flow is in approximate inflow equilibrium. We find that the time-averaged radial-dependence of fluid-frame electromagnetic stress is almost completely independent of both disk thickness and poloidal magnetic topology. It rises smoothly inward at all radii (exhibiting no feature associated with the ISCO) until just outside the event horizon, where the stress plummets to zero. Reynolds stress can also be significant near the ISCO and in the plunging region; the magnitude of this stress, however, depends on both disk thickness and magnetic topology. The two stresses combine to make the net angular momentum accreted per unit rest-mass 7-15% less than the angular momentum of the ISCO.

We study the properties of the turbulence driven by the magnetorotational instability (MRI) in a stratified shearing box with outflow boundary conditions and an equation of state determined by self-consistent dissipation and radiation losses. A series of simulations with increasing resolution are performed within a fixed computational box. We achieve numerical convergence with respect to radial and azimuthal resolution. As vertical resolution is improved, the ratio of stress to pressure increases only slowly, but the absolute levels of both the stress and the pressure increase noticeably. The vertical correlation length of the magnetic field within the core of the disk is highly time-dependent, but averaged over time, it is $\simeq 3$ times larger than found in previous unstratifed simulations. This correlation length decreases slowly as vertical resolution increases. We suggest that this greater degree of correlation in stratified than unstratified disks, even near the midplane, is due to field buoyancy. We further show that the undulatory Parker instability drives the magnetic field upwelling at several scaleheights from the midplane that is characteristic of stratified MHD-turbulent disks.

We present new calculations of X-ray polarization from accreting black holes (BHs), using a Monte-Carlo ray-tracing code in full general relativity. In our model, an optically thick disk in the BH equatorial plane produces thermal seed photons with polarization oriented parallel to the disk surface. These seed photons are then inverse-Compton scattered through a hot (but thermal) corona, producing a hard X-ray power-law spectrum. We consider three different models for the corona geometry: a wedge ``sandwich'' with aspect ratio H/R and vertically-integrated optical depth tau_0 constant throughout the disk; an inhomogeneous ``clumpy'' corona with a finite number of hot clouds distributed randomly above the disk within a wedge geometry; and a spherical corona of uniform density, centered on the BH and surrounded by a truncated thermal disk with inner radius R_edge. In all cases we find a characteristic transition from horizontal polarization at low energies to vertical polarization above the thermal peak; the vertical direction is defined as the projection of the BH spin axis on the plane of the sky. We show how the details of the spectropolarization signal can be used to distinguish between these models and infer various properties of the corona and BH. Although the bulk of this paper focuses on stellar-mass BHs, we also consider the effects of coronal scattering on the X-ray polarization signal from supermassive BHs in active galactic nuclei.

Although recent work in numerical relativity has made tremendous strides in quantifying the gravitational wave luminosity of black hole mergers, very little is known about the electromagnetic luminosity that might occur in immediate conjunction with these events. We show that whenever the heat deposited in the gas near a pair of merging black holes is proportional to its total mass, and the surface density of the gas in the immediate vicinity is greater than the (quite small) amount necessary to make it optically thick, the characteristic scale of the luminosity emitted in direct association with the merger is the Eddington luminosity independent of the gas mass. The duration of the photon signal is proportional to the gas mass, and is generally rather longer than the merger event. At somewhat larger distances, dissipation associated with realigning the gas orbits to the new spin orientation of the black hole can supplement dissipation of the energy gained from orbital adjustment to the mass lost in gravitational radiation; these two heat sources can combine to augment the electromagnetic radiation over longer timescales.

Explorers have made breakthroughs in many fields of astrophysics. The science from both these missions contributed to three Nobel Prizes - Giacconi (2002), Mather, and Smoot (2006). Explorers have: marked the definitive beginning of precision cosmology, discovered that short gamma-ray bursts are caused by compact star mergers and have measured metalicity to redshifts z>6. NASA Explorers do cutting-edge science that cannot be done by facility-class instruments. The Explorer program provides a rapid response to changing science and technology, to enable cutting-edge science at moderate cost. Explorers also enable innovation, and engage & train scientists, managers and engineers, adding human capital to NASA and the nation. The astrophysics Explorer launch rate now being achieved is 1 per 3 years, and budget projections are in the $150M/year range for the next five years. A newly Vigorous Explorer Program should be created to: 1. Reach the long-stated goal of annual astrophysics launches; 2. Find additional launch options for Explorers and actively encourage cost savings in launchers and spacecraft, such as new commercial vehicles and innovative partnerships. 3. Mitigate risk via stronger technical development and sub-orbital programs, and through longer, more thorough, Phase A programs, potentially reducing the need for a 30% contingency; 4. Strive to protect the funding for missions that have reached Phase B, to prevent significant launch slips and cancellations, with a goal of 4 to 5 years from Phase B to launch; 5. Review the project management procedures and requirements to seek cost reductions, including the risk management strategy and the review and reporting process; 6. Review and possibly modify the cost caps for all Explorer classes to optimize scientific returns per dollar. [ABRIDGED]

We report on a global, three-dimensional GRMHD simulation of an accretion torus embedded in a large scale vertical magnetic field orbiting a Schwarzschild black hole. This simulation investigates how a large scale vertical field evolves within a turbulent accretion disk and whether global magnetic field configurations suitable for launching jets and winds can develop. We find that a "coronal mechanism" of magnetic flux motion, which operates largely outside the disk body, dominates global flux evolution. In this mechanism, magnetic stresses driven by orbital shear create large-scale half-loops of magnetic field that stretch radially inward and then reconnect, leading to discontinuous jumps in the location of magnetic flux. In contrast, little or no flux is brought in directly by accretion within the disk itself. The coronal mechanism establishes a dipole magnetic field in the evacuated funnel around the orbital axis with a field intensity regulated by a combination of the magnetic and gas pressures in the inner disk. These results prompt a reevaluation of previous descriptions of magnetic flux motion associated with accretion. Local pictures are undercut by the intrinsically global character of magnetic flux. Formulations in terms of an "effective viscosity" competing with an "effective resistivity" are undermined by the nonlinearity of of the magnetic dynamics and the fact that the same turbulence driving mass motion (traditionally identified as "viscosity") can alter magnetic topology.

Magnetic fields connecting the immediate environs of rotating black holes to large distances appear to be the most promising mechanism for launching relativistic jets, an idea first developed by Blandford and Znajek in the mid-1970s. To enable an understanding of this process, we provide a brief introduction to dynamics and electromagnetism in the spacetime near black holes. We then present a brief summary of the classical Blandford-Znajek mechanism and its conceptual foundations. Recently, it has become possible to study these effects in much greater detail using numerical simulation. After discussing which aspects of the problem can be handled well by numerical means and which aspects remain beyond the grasp of such methods, we summarize their results so far. Simulations have confirmed that processes akin to the classical Blandford-Znajek mechanism can launch powerful electromagnetically-dominated jets, and have shown how the jet luminosity can be related to black hole spin and concurrent accretion rate. However, they have also shown that the luminosity and variability of jets can depend strongly on magnetic field geometry. We close with a discussion of several important open questions.

We present the results of a series of radiation-MHD simulations of a local patch of an accretion disk, with fixed vertical gravity profile but with different surface mass densities and a broad range of radiation to gas pressure ratios. Each simulation achieves a thermal equilibrium that lasts for many cooling times. After averaging over times long compared to a cooling time, we find that the vertically integrated stress is approximately proportional to the vertically-averaged total thermal (gas plus radiation) pressure. We map out--for the first time on the basis of explicit physics--the thermal equilibrium relation between stress and surface density: the stress decreases (increases) with increasing surface mass density when the simulation is radiation (gas) pressure dominated. The dependence of stress on surface mass density in the radiation pressure dominated regime suggests the possibility of a Lightman-Eardley inflow instability, but global simulations or shearing box simulations with much wider radial boxes will be necessary to confirm this and determine its nonlinear behavior.

We present the results of one and three-dimensional radiative transfer calculations of polarized spectra emerging from snapshots of radiation magnetohydrodynamical simulations of the local vertical structure of black hole accretion disks. The simulations cover a wide range of physical regimes relevant for the high/soft state of black hole X-ray binaries. We constrain the uncertainties in theoretical spectral color correction factors due to the presence of magnetic support of the disk surface layers and strong density inhomogeneities. For the radiation dominated simulation, magnetic support increases the color correction factor by about ten percent, but this is largely compensated by a ten percent softening due to inhomogeneities. We also compute the effects of inhomogeneities and Faraday rotation on the resulting polarization. Magnetic fields in the simulations are just strong enough to produce significant Faraday depolarization near the spectral peak of the radiation field. X-ray polarimetry may therefore be a valuable diagnostic of accretion disk magnetic fields, being able to directly test simulations of magnetorotational turbulence.

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

On the basis of data from an energy-conserving 3D general relativistic MHD simulation, we predict the statistical character of variability in the coronal luminosity from accreting black holes. When the inner boundary of the corona is defined to be the electron scattering photosphere, its location depends only on the mass accretion rate in Eddington units (\dot{M}). Nearly independent of viewing angle and \dot{M}, the power spectrum over the range of frequencies from approximately the orbital frequency at the innermost stable circular orbit (ISCO) to ~100 times lower is well approximated by a power-law with index -2, crudely consistent with the observed power spectra of hard X-ray fluctuations in AGN and the hard states of Galactic binary black holes. The underlying physical driver for variability in the light curve is variations in the accretion rate caused by the chaotic character of MHD turbulence, but the power spectrum of the coronal light output is significantly steeper. Part of this contrast is due to the fact that the mass accretion rate can be significantly modulated by radial epicyclic motions that do not result in dissipation, and therefore do not drive luminosity fluctuations. The other part of this contrast is due to the inward decrease of the characteristic inflow time, which leads to decreasing radial coherence length with increasing fluctuation frequency.

We present new calculations of X-ray polarization from black hole accretion disks in the thermally-dominated state, using a Monte-Carlo ray-tracing code in full general relativity. In contrast to many previously published studies, our approach allows us to include returning radiation that is deflected by the strong-field gravity of the BH and scatters off of the disk before reaching a distant observer. Although carrying a relatively small fraction of the total observed flux, the scattered radiation tends to be highly polarized and in a direction perpendicular to the direct radiation. We show how these new features of the polarization spectra may be developed into a powerful tool for measuring black hole spin and probing the gas flow in the innermost disk.

We present new calculations of X-ray polarization from black hole (BH) accretion disks in the thermally-dominated state, using a Monte-Carlo ray-tracing code in full general relativity. In contrast to many previously published studies, our approach allows us to include returning radiation that is deflected by the strong-field gravity of the BH and scatters off of the disk before reaching a distant observer. Although carrying a relatively small fraction of the total observed flux, the scattered radiation tends to be highly polarized and in a direction perpendicular to the direct radiation. For moderately large spin parameters (a/M >~ 0.9), this scattered returning radiation dominates the polarization signal at energies above the thermal peak, giving a net rotation in the polarization angle of 90 deg. We show how these new features of the polarization spectra from BHs in the thermal state may be developed into a powerful tool for measuring BH spin and probing the gas flow in the innermost disk. In addition to determining the emission profile, polarization observations can be used to constrain other properties of the system such as BH mass, inclination, and distance. New instruments currently under development should be able to exploit this tool in the near future.

SDSSJ092712.65+294344.0 was identified by the SDSS as a quasar, but has the unusual property of having two emission line systems offset by 2650 km/s. One of these contains the usual combination of broad and narrow lines, the other only narrow lines. In the first paper commenting on this system (Komossa et al. 2008), it was interpreted as a galaxy in which a pair of black holes had merged, imparting a several thousand km/s recoil to the new, larger black hole. In two other papers (Bogdanovic, Eracleous & Sigurdsson 2008; Dotti et al. 2008), it was interpreted as a small-separation binary black hole. We propose a new interpretation: that this system is a more distant analog of NGC1275, a large and a small galaxy interacting near the center of a rich cluster.

When the accretion rate is more than a small fraction of Eddington, the inner regions of accretion disks around black holes are expected to be radiation-dominated. However, in the alpha-model, these regions are also expected to be thermally unstable. In this paper, we report two 3-d radiation MHD simulations of a vertically-stratified shearing box in which the ratio of radiation to gas pressure is ~ 10, and yet no thermal runaway occurs over a timespan ~ 40 cooling times. Where the time-averaged dissipation rate is greater than the critical dissipation rate that creates hydrostatic equilibrium by diffusive radiation flux, the time-averaged radiation flux is held to the critical value, with the excess dissipated energy transported by radiative advection. Although the stress and total pressure are well-correlated as predicted by the alpha-model, we show that stress fluctuations precede pressure fluctuations, contrary to the usual supposition that the pressure controls the saturation level of the magnetic energy. This fact explains the thermal stability. Using a simple toy-model, we show that independently-generated magnetic fluctuations can drive radiation pressure fluctuations, creating a correlation between the two while maintaining thermal stability.

Numerical simulation of magnetohydrodynamic (MHD) turbulence makes it possible to study accretion dynamics in detail. However, special effort is required to connect inflow dynamics (dependent largely on angular momentum transport) to radiation (dependent largely on thermodynamics and photon diffusion). To this end we extend the flux-conservative, general relativistic MHD code HARM from axisymmetry to full 3D. The use of an energy conserving algorithm allows the energy dissipated in the course of relativistic accretion to be captured as heat. The inclusion of a simple optically thin cooling function permits explicit control of the simulated disk's geometric thickness as well as a direct calculation of both the amplitude and location of the radiative cooling associated with the accretion stresses. Fully relativistic ray-tracing is used to compute the luminosity received by distant observers. For a disk with aspect ratio H/r ~ 0.1 accreting onto a black hole with spin parameter a/M = 0.9, we find that there is significant dissipation beyond that predicted by the classical Novikov-Thorne model. However, much of it occurs deep in the potential, where photon capture and gravitational redshifting can strongly limit the net photon energy escaping to infinity. In addition, with these parameters and this radiation model, significant thermal and magnetic energy remains with the gas and is accreted by the black hole. In our model, the net luminosi ty reaching infinity is 6% greater than the Novikov-Thorne prediction. If the accreted thermal energy were wholly radiated, the total luminosity of the accretion flow would be ~20% greater than the Novikov-Thorne value.

Type 2 quasars are luminous Active Galactic Nuclei whose central engines are seen through large amounts of gas and dust. We present Spitzer spectra of twelve type 2 quasars selected on the basis of their optical emission line properties. Within this sample, we find a surprising diversity of spectra, from those that are featureless to those showing strong PAH emission, deep silicate absorption at 10 micron, hydrocarbon absorption, high-ionization emission lines and H_2 rotational emission lines. About half of the objects in the sample are likely Compton-thick, including the two with the deepest Si absorption. The median star-formation luminosity of the objects in our sample measured from the strength of the PAH features is 5x10^11 L_sun, much higher than for field galaxies or for any other AGN sample, but similar to other samples of type 2 quasars. This suggests an evolutionary link between obscured quasars and peak star formation activity in the host galaxy. Despite the high level of star formation, the bolometric output is dominated by the quasar in all cases. For a given strength of 10 micron Si absorption, ULIRGs are significantly colder than are type 2 quasars (their F_nu[14.5 micron]/F_nu[27.5 micron] ratio is 0.5 dex lower), perhaps reflecting different obscuration geometries in these sources. We find that the appearance of the 10 micron feature (i.e., whether it shows in emission or in absorption) is well-correlated with the optical classification in type 1 and type 2 quasars, contrary to some models of clumpy obscuration. Furthermore, this correlation is significantly stronger in quasars (L_bol>10^45 erg/s) than it is in Seyfert galaxies (L_bol<<10^45 erg/s).

General Relativistic (GR) Magnetohydrodynamic (MHD) simulations of black hole accretion find significant magnetic stresses near and inside the innermost stable circular orbit (ISCO), suggesting that such flows could radiate in a manner noticeably different from the prediction of the standard model, which assumes that there are no stresses in that region. We provide estimates of how phenomenologically interesting parameters like the ``radiation edge", the innermost ring of the disc from which substantial thermal radiation escapes to infinity, may be altered by stresses near the ISCO. These estimates are based on data from a large number of three-dimensional GRMHD simulations combined with GR ray-tracing. For slowly spinning black holes ($a/M<0.9$), the radiation edge lies well inside where the standard model predicts, particularly when the system is viewed at high inclination. For more rapidly spinning black holes, the contrast is smaller. At fixed total luminosity, the characteristic temperature of the accretion flow increases between a factor of $1.2-2.4$ over that predicted by the standard model, whilst at fixed mass accretion rate, there is a corresponding enhancement of the accretion luminosity which may be anywhere from tens of percent to order unity. When all these considerations are combined, we find that, for fixed black hole mass, luminosity, and inclination angle, our uncertainty in the characteristic temperature of the radiation reaching distant observers due to uncertainty in dissipation profile (around a factor of 3) is {\it greater} than the uncertainty due to a complete lack of knowledge of the black hole's spin (around a factor of 2) and furthermore that spin estimates based on the stress-free inner boundary condition provide an upper limit to $a/M$.

We model the spectra and light curves of circumbinary accretion disks during the time after the central black holes merge. The most immediate effect of this merger is the dissipation of energy in the outer regions of the disk due to the gravitational wave energy and linear momentum flux released at merger. This has the effect of perturbing the gas in the disk, which then radiates the dissipated energy over a cooling timescale, giving a characteristic infrared signal for tens of thousands of years when the total black hole mass is M~10^8 M_sun. On the basis of a simple cosmological merger model in which a typical supermassive black hole undergoes a few major mergers during its lifetime, we predict that ~10^4-10^5 of these IR sources should be observable today and discuss the possibility of identifying them with multi-wavelength surveys such as SWIRE/XMM-LSS/XBootes and COSMOS.

We have chosen 6 Ultra-Luminous X-ray sources from the {\it XMM-Newton} archive whose spectra have high signal-to-noise and can be fitted solely with a disk model without requiring any power-law component. To estimate systematic errors in the inferred parameters, we fit every spectrum to two different disk models, one based on local blackbody emission (KERRBB) and one based on detailed atmosphere modelling (BHSPEC). Both incorporate full general relativistic treatment of the disk surface brightness profile, photon Doppler shifts, and photon trajectories. We found in every case that they give almost identical fits and similar acceptable parameters. The best-fit value of the most interesting parameter, the mass of the central object, is between 23 and 73 M$_\sun$ in 5 of the 6 examples. In every case, the best-fit inclination angle and mass are correlated, in the sense that large mass corresponds to high inclination. Even after allowing for this degeneracy, we find that, with $\gtrsim 99.9%$ formal statistical confidence, 3 of the 6 objects have mass $\gtrsim 25 M_\sun$; for the other 3, these data are consistent with a wide range of masses. A mass greater than several hundred $M_\sun$ is unlikely for the 3 best-constrained objects. These fits also suggest comparatively rapid black hole spin in the 3 objects whose masses are relatively well-determined, but our estimate of the spin is subject to significant systematic error having to do with uncertainty in the underlying surface brightness profile.

The dynamics and structure of toroidal obscuration around AGN remain uncertain and controversial. In this paper we extend earlier work on the dynamical role of infrared radiation pressure by adding the effects of two kinds of distributed heating: Compton-heating due to hard X-rays from the nucleus and local starlight heating. We find numerical solutions to the axisymmetric hydrostatic equilibrium, energy balance, and photon diffusion equations including these effects. Within the regime of typical parameters, the two different sources of additional heating have very similar effects: the density profile within the torus becomes shallower both radially and vertically, but for plausible heating rates, there is only minor change (relative to the source-free case) in the distribution of column density with solid angle. The most interesting consequence of distributed heating is that it selects out a relatively narrow range of parameters permitting an equilibrium, particularly $(L/L_E)/\tau_T$. We discuss the implications of both the narrowness of the permitted range and its approximate coincidence with the range inferred from observations.

Recent work has suggested that the fraction of obscured AGN declines with increasing luminosity, but it has been difficult to quantify this trend. Here, we attempt to measure this fraction as a function of luminosity by studying the ratio of mid-infrared to intrinsic nuclear bolometric luminosity in unobscured AGN. Because the mid-infrared is created by dust reprocessing of shorter wavelength nuclear light, this ratio is a diagnostic of f_obsc, the fraction of solid angle around the nucleus covered by obscuring matter. In order to eliminate possible redshift-dependences while also achieving a large dynamic range in luminosity, we have collected archival 24 micron MIPS photometry from objects with z~1 in the Sloan Digital Sky Survey (SDSS), the Great Observatories Origins Deep Survey (GOODS) and the Cosmic Evolution Survey (COSMOS). To measure the bolometric luminosity for each object, we used archival optical data supplemented by GALEX data. We find that the mean ratio of 24 microns to bolometric luminosity decreases by a factor of ~3 in the L_bol=10^44-3x10^47 ergs s^-1 range, but there is also a large scatter at constant L_bol. Using radiation transfer solutions for model geometries, we show how the IR/bolometric ratio relates to f_obsc and compare these values with those obtained obtained from samples of X-ray selected AGN. Although we find approximate agreement, our method indicates somewhat higher values of f_obsc, particularly in the middle range of luminosities, suggesting that there may be a significant number of heavily obscured AGN missed by X-ray surveys.

Because the magneto-rotational instability is capable of exponentially amplifying weak preexisting magnetic fields, it might be hoped that the character of the magnetic field in accretion disks is independent of the nature of the seed field. However, the divergence-free nature of magnetic fields in highly conducting fluids ensures that their large-scale topology is preserved, no matter how greatly the field intensity is changed. By performing global two-dimensional and three-dimensional general relativistic magnetohydrodynamic disk simulations with several different topologies for the initial magnetic field, we explore the degree to which the character of the flows around black holes depends on the initial topology. We find that while the qualitative properties of the accretion flow are nearly independent of field topology, jet-launching is very sensitive to it: a sense of vertical field consistent for at least an inner disk inflow time is essential to the support of strong jets.

With the fundamental stress mechanism of accretion disks identified--correlated MHD turbulence driven by the magneto-rotational instability--it has become possible to make numerical simulations of accretion disk dynamics based on well-understood physics. A sampling of results from both Newtonian 3-d shearing box and general relativistic global disk MHD simulations is reported. Among other things, these simulations have shown that: contrary to long-held assumptions, stress is continuous through the marginally stable and plunging regions around black holes, so that rotating black holes can give substantial amounts of angular momentum electromagnetically to surrounding matter; the upper layers of accretion disks are primarily supported by magnetic pressure, potentially leading to interesting departures from local black-body emitted spectra; and initially local magnetic fields in accretion flows can, in some cases, spontaneously generate large-scale fields that connect rotating black holes to infinity and mediate strong relativistic jets.

We explore the thermodynamic and global structural properties of a local patch of an accretion disk whose parameters were chosen so that radiation pressure and gas pressure would be comparable in magnitude. Heating, radiative transport, and cooling are computed self-consistently with the structure by solving the equations of radiation MHD in the shearing-box approximation. Using a fully 3-d and energy-conserving code, we compute the structure and energy balance of this disk segment over a span of more than forty cooling times. As is also true when gas pressure dominates, the disk's upper atmosphere is magnetically-supported. However, unlike the gas-dominated case, no steady-state is reached; instead, the total (i.e., radiation plus gas) energy content fluctuates by factors of 3--4 over timescales of several tens of orbits, with no secular trend. Because the radiation pressure varies much more than the gas pressure, the ratio of radiation pressure to gas pressure varies over the approximate range 0.5--2. The volume-integrated dissipation rate generally increases with increasing total energy, but the mean trend is somewhat slower than linear, and the instantaneous dissipation rate is often a factor of two larger or smaller than the mean for that total energy level. Locally, the dissipation rate per unit volume scales approximately in proportion to the current density; the time-average dissipation rate per unit mass is proportional to m^{-1/2}, where m is the horizontally-averaged mass column density to the nearer of the top or bottom surface. As in our earlier study of a gas-dominated shearing-box, we find that energy transport is completely dominated by radiative diffusion, with Poynting flux carrying less than 1% of the energy lost from the box.

We have employed a 3-d energy-conserving radiation MHD code to simulate the vertical structure and thermodynamics of a shearing box whose parameters were chosen so that the radiation and gas pressures would be comparable. The upper layers of this disk segment are magnetically-dominated, creating conditions appropriate for both photon bubble and Parker instabilities. We find little evidence for photon bubbles, even though the simulation has enough spatial resolution to see them and their predicted growth rates are high. On the other hand, there is strong evidence for Parker instabilities, and they appear to dominate the evolution of the magnetically supported surface layers. The disk photosphere is complex, with large density inhomogeneities at both the scattering and effective (thermalization) photospheres of the evolving horizontally-averaged structure. Both the dominant magnetic support and the inhomogeneities are likely to have strong effects on the spectrum and polarization of thermal photons emerging from the disk atmosphere. The inhomogeneities are also large enough to affect models of reflection spectra from the atmospheres of accretion disks.

Explicit 2-d axisymmetric solutions are found to the hydrostatic equilibrium, energy balance, and photon diffusion equations within obscuring tori around active galactic nuclei. These solutions demonstrate that infrared radiation pressure can support geometrically thick structures in AGN environments subject to certain constraints: the bolometric luminosity must be roughly 0.03--1 times the Eddington luminosity; and the Compton optical depth of matter in the equatorial plane should be order unity, with a tolerance of about an order of magnitude up or down. Both of these constraints are at least roughly consistent with observations. In addition, angular momentum must be redistributed so that the fractional rotational support against gravity rises from the inner edge of the torus to the outer in a manner specific to the detailed shape of the gravitational potential. This model also predicts that the column densities observed in obscured AGN should range from about 10^{22} to about 10^{24} cm^{-2}.

The inward flow of matter through accretion disks is driven by MHD turbulence. Global general relativistic MHD simulations shed quantitative light on this process, revealing a number of aspects of accretion previously unrecognized. Among them are strong stresses in the marginally stable and plunging regions near the black hole and electromagnetically-dominated conical relativistic jets that can form spontaneously from the accretion flow. The energy release associated with both of these effects can significantly augment the classical energy-release estimates based on purely hydrodynamic models.

Type II quasars are luminous Active Galactic Nuclei whose centers are obscured by large amounts of gas and dust. In this paper we present 3-band HST images of nine type II quasars with redshifts 0.2 < z < 0.4 selected from the Sloan Digital Sky Survey based on their emission line properties. The intrinsic luminosities of these AGN are estimated to be -24 > M_B > -26, but optical obscuration allows their host galaxies to be studied unencumbered by bright nuclei. Each object has been imaged in three continuum filters (`UV', `blue' and `yellow') placed between the strong emission lines. The spectacular, high quality images reveal a wealth of details about the structure of the host galaxies and their environments. Six of the nine galaxies in the sample are ellipticals with de Vaucouleurs light profiles, one object has a well-defined disk component and the remaining two have marginal disks. Stellar populations of type II quasar hosts are more luminous (by a median of 0.3-0.7 mag, depending on the wavelength) and bluer (by about 0.4 mag) than are M* galaxies at the same redshift. When smooth fits to stellar light are subtracted from the images, we find both positive and negative residuals that become more prominent toward shorter wavelengths. We argue that the negative residuals are due to kpc-scale dust obscuration, while most positive residuals are due to the light from the nucleus scattered off interstellar material in the host galaxy. Scattered light makes a significant contribution to the broad band continuum emission and can be the dominant component of the extended emission in the UV in extreme cases.

We use a relativistic ray-tracing code to calculate the light curves observed from a global general relativistic magneto-hydrodynamic simulation of an accretion flow onto a Schwarzschild black hole. We apply three basic emission models to sample different properties of the time-dependent accretion disk. With one of these models, which assumes thermal blackbody emission and free-free absorption, we can predict qualitative features of the high-frequency power spectrum from stellar-mass black holes in the "Thermal Dominant" state. The simulated power spectrum is characterized by a power law of index Gamma ~ 3 and total rms fractional variance of <~ 2% above 10 Hz. For each emission model, we find that the variability amplitude should increase with increasing inclination angle. On the basis of a newly-developed formalism for quantifying the significance of quasi-periodic oscillations (QPOs) in simulation data, we find that these simulations are able to identify any such features with (rms/mean) amplitudes >~ 1 % near the orbital frequency at the inner-most stable orbit. Initial results indicate the existence of transient QPO peaks with frequency ratios of nearly 2:3 at a 99.9% confidence limit, but they are not generic features because at any given time they are seen only from certain observer directions. Additionally, we present detailed analysis of the azimuthal structure of the accretion disk and the evolution of density perturbations in the inner disk. These "hot spot" structures appear to be roughly self-similar over a range of disk radii, with a single characteristic size \delta\phi=25 deg and \delta r/r=0.3, and typical lifetimes T_l ~ 0.3 T_orb.

We analyze observations obtained with the Chandra X-ray Observatory of bright Compton thick active galactic nuclei (AGNs), those with column densities in excess of 1.5 x 10^{24} cm^{-2} along the lines of sight. We therefore view the powerful central engines only indirectly, even at X-ray energies. Using high spatial resolution and considering only galaxies that do not contain circumnuclear starbursts, we reveal the variety of emission AGNs alone may produce. Approximately 1% of the continuum's intrinsic flux is detected in reflection in each case. The only hard X-ray feature is the prominent Fe K alpha fluorescence line, with equivalent width greater than 1 keV in all sources. The Fe line luminosity provides the best X-ray indicator of the unseen intrinsic AGN luminosity. In detail, the morphologies of the extended soft X-ray emission and optical line emission are similar, and line emission dominates the soft X-ray spectra. Thus, we attribute the soft X-ray emission to material that the central engines photoionize. Because the resulting spectra are complex and do not reveal the AGNs directly, crude analysis techniques such as hardness ratios would mis-classify these galaxies as hosts of intrinsically weak, unabsorbed AGNs and would fail to identify the luminous, absorbed nuclei that are present. We demonstrate that a three-band X-ray diagnostic can correctly classify Compton thick AGNs, even when significant soft X-ray line emission is present. The active nuclei produce most of the galaxies' total observed emission over a broad spectral range, and much of their light emerges at far-infrared wavelengths. Stellar contamination of the infrared emission can be severe, however, making long-wavelength data alone unreliable indicators of the buried AGN luminosity.

Simulations of black hole accretion have shown that magnetic stresses are present near and inside the innermost stable circular orbit (ISCO). This finding suggests that such flows may be more luminous than predicted by the standard accretion disk model. Here we apply a prescription for heat dissipation within the simulated accretion flows to estimate their implied radiative efficiency. We assume that dissipation is proportional to the current density squared, and find that the resulting azimuthally-averaged and shell-integrated radial profile is well-matched to the radial heat dissipation profile of the standard disk model for the region outside the ISCO, particularly when it is adjusted to account for additional stress at the ISCO. In contrast to the standard model, however, the dissipation profile derived from the current density continues past the ISCO and through the plunging region. The total predicted dissipation rate is between $\simeq 30%$ and $\simeq 100%$ greater than that predicted by the standard model, depending on the black hole spin. Most of the additional dissipation takes place just outside the ISCO. To predict luminosities, we assume instantaneous radiation and zero optical depth, but allow for photon capture. The net radiative efficiency seen by a distant observer is increased relative to the standard model by $\simeq 25%$--80%, with the largest fractional increase for intermediate black hole spins because the increase in dissipation from enhanced stress that occurs for rapid spin is partially offset by the increased likelihood that the additional photons will be captured by the hole.

We investigate the behavior of the magneto-rotational instability in the limit of extremely weak magnetic field, i.e., as the ratio of ion cyclotron frequency to orbital frequency (X) becomes small. Considered only in terms of cold two-fluid theory, instability persists to arbitrarily small values of X, and the maximum growth rate is of order the orbital frequency except for the range m_e/m_i < |X| < 1, where it can be rather smaller. In this range, field aligned with rotation (X > 0) produces slower growth than anti-aligned field (X < 0). The maximum growth rate is generally achieved at smaller and smaller wavelengths as |X| diminishes. When |X| < m_e/m_i, new unstable "electromagnetic-rotational" modes appear that do not depend on the equilibrium magnetic field. Because the most rapidly-growing modes have extremely short wavelengths when |X| is small, they are often subject to viscous or resistive damping, which can result in suppressing all but the longest wavelengths, for which growth is much slower. We find that this sort of damping is likely to curtail severely the frequently-invoked mechanism for cosmological magnetic field growth in which a magnetic field seeded by the Biermann battery is then amplified by the magneto-rotational instability. On the other hand, the small |X| case may introduce interesting effects in weakly-ionized disks in which dust grains carry most of the electric charge.

Type II quasars are luminous Active Galactic Nuclei whose centers are obscured by large amounts of gas and dust. In this contribution we present 3-band HST images of nine type II quasars with redshifts 0.25<z<0.4 selected from the Sloan Digital Sky Survey based on their emission line properties. The intrinsic luminosities of these quasars are thought to be in the range -24>M_B>-26, but optical obscuration implies that host galaxies can be studied unencumbered by bright nuclei. Each object has been imaged in three filters (`red', `green' and `blue') placed between the strong emission lines. The spectacular, high quality images reveal a wealth of details about the structure of the host galaxies and their environments. Most galaxies in the sample are ellipticals, but strong deviations from de Vaucouleurs profiles are found, especially in the blue band. We argue that most of these deviations are due to the light from the nucleus scattered off interstellar material in the host galaxy. This scattered component can make a significant contribution to the broad-band flux and complicates the analysis of the colors of the stellar populations in the host galaxy. This extended component can be difficult to notice in unobscured luminous quasars and may bias the results of host galaxy studies.

Stellar atmosphere models of ionized accretion disks have generally neglected the contribution of magnetic fields to the vertical hydrostatic support, although magnetic fields are widely believed to play a critical role in the transport of angular momentum. Simulations of magnetorotational turbulence in a vertically stratified shearing box geometry show that magnetic pressure support can be dominant in the upper layers of the disk. We present calculations of accretion disk spectra that include this magnetic pressure support, as well as a vertical dissipation profile based on simulation. Magnetic pressure support generically produces a more vertically extended disk atmosphere with a larger density scale height. This acts to harden the spectrum compared to models that neglect magnetic pressure support. We estimate the significance of this effect on disk-integrated spectra by calculating an illustrative disk model for a stellar mass black hole, assuming that similar magnetic pressure support exists at all radii.

We compute a series of three-dimensional general relativistic magnetohydrodynamic simulations of accretion flows in the Kerr metric to investigate the properties of the unbound outflows that result. The overall strength of these outflows increases sharply with increasing black hole rotation rate, but a number of generic features are found in all cases. The mass in the outflow is concentrated in a hollow cone whose opening angle is largely determined by the effective potential for matter orbiting with angular momentum comparable to that of the innermost stable circular orbit. The dominant force accelerating the matter outward comes from the pressure of the accretion disk's corona. The principal element that shapes the outflow is therefore the centrifugal barrier preventing accreting matter from coming close to the rotation axis. Inside the centrifugal barrier, the cone contains very little matter and is dominated by electromagnetic fields that rotate at a rate tied closely to the rotation of the black hole. These fields carry an outward-going Poynting flux whose immediate energy source is the rotating spacetime of the Kerr black hole. When the spin parameter a/M of the black hole exceeds ~0.9, the energy carried to infinity by these outflows can be comparable to the nominal radiative efficiency predicted in the Novikov-Thorne model. Similarly, the expelled angular momentum can be comparable to that accreted by the black hole. Both the inner electromagnetic part and the outer matter part can contribute in significant fashion to the energy and angular momentum of the outflow.

(shortened) We calculate the vertical structure of a local patch of an accretion disk in which heating by dissipation of MRI-driven MHD turbulence is balanced by radiative cooling. Heating, radiative transport, and cooling are computed self-consistently with the structure by solving the equations of radiation MHD in the shearing-box approximation. Using a fully 3-d and energy-conserving code, we compute the structure of this disk segment over a span of more than five cooling times. After a brief relaxation period, a statistically steady-state develops. Measuring height above the midplane in units of the scale-height H predicted by a Shakura-Sunyaev model, we find that the disk atmosphere stretches upward, with the photosphere rising to about 7H, in contrast to the approximately 3H predicted by conventional analytic models. This more extended structure, as well as fluctuations in the height of the photosphere, may lead to departures from Planckian form in the emergent spectra. Dissipation is distributed across the region within roughly 3H of the midplane, but is very weak at greater altitudes. Because fluctuations in the dissipation are particularly strong away from the midplane, the emergent radiation flux can track dissipation fluctuations with a lag that is only 0.1--0.2 times the mean cooling time of the disk. Long timescale asymmetries in the dissipation distribution can also cause significant asymmetry in the flux emerging from the top and bottom surfaces of the disk. Radiative diffusion dominates Poynting flux in the vertical energy flow throughout the disk.

We are carrying out sensitive X-ray observations with Chandra and XMM of type II quasars selected from the Sloan Digital Sky Survey based on their optical emission line properties. We present observations of four objects at redshifts 0.4 < z < 0.8 and an analysis of the archival data for four additional objects in the same redshift range. Six of the eight were detected in X-rays; five of them have sufficient signal to derive spectral information. All of the detected sources have intrinsic luminosities L(2-10 keV) > 5 x 10^43 erg s^-1. The five with sufficient counts for spectral fitting show evidence for significant absorption (N_H >~ a few x 10^22 cm^-2). At least three of the objects likely have N_H > 10^23 cm^-2; some may be Compton-thick (N_H > 10^{24} cm^-2). In the five objects for which we could fit spectra, the slopes tend to be significantly flatter than is typically observed in AGN; it is possible that this is due either to reprocessing of the nuclear emission or to a line of sight that passes through patchy absorption.

We report results from {\it XMM-Newton} observations of thirteen X-ray bright BL Lacertae objects, selected from the {\it Einstein} Slew Survey sample. The survey was designed to look for evidence of departures of the X-ray spectra from a simple power law shape (i.e., curvature and/or line features), and to find objects worthy of deeper study. Our data are generally well fit by power-law models, with three cases having hard ($\Gamma<2; dN/dE \propto E^{-\Gamma}$) spectra that indicate synchrotron peaks at $E \gsim 5$ keV. Previous data had suggested a presence of absorption features in the X-ray spectra of some BL Lacs. In contrast, none of these spectra show convincing examples of line features, either in absorption or emission, suggesting that such features are rare amongst BL Lacs, or, more likely, artifacts caused by instrumental effects. We find significant evidence for intrinsic curvature (steepening by $d\Gamma / d({\rm log} E) = 0.4 \pm 0.15$) in fourteen of the seventeen X-ray spectra. This cannot be explained satisfactorily via excess absorption, since the curvature is essentially constant from $0.5-6$ keV, an observation which is inconsistent with the modest amounts of absorption that would be required. We use the {\it XMM-Newton} Optical Monitor data with concurrent radio monitoring to derive broadband spectral energy distributions and peak frequency estimates. From these we examine models of synchrotron emission and model the spectral curvature we see as the result of episodic particle acceleration.

We have calculated the structures and the emergent spectra of stationary, geometrically thin accretion disks around 100 and 1000 M_sun black holes in both the Schwarzschild and extreme Kerr metrics. Equations of radiative transfer, hydrostatic equilibrium, energy balance, ionization equilibrium, and statistical equilibrium are solved simultaneously and consistently. The six most astrophysically abundant elements (H, He, C, N, O, and Fe) are included, as well as energy transfer by Comptonization. The observed spectrum as a function of viewing angle is computed incorporating all general relativistic effects. We find that, in contrast with the predictions of the commonly-used multi-color disk (MCD) model, opacity associated with photoionization of heavy elements can significantly alter the spectrum near its peak. These ionization edges can create spectral breaks visible in the spectra of slowly-spinning black holes viewed from almost all angles and in the spectra of rapidly-spinning black holes seen approximately pole-on. For fixed mass and accretion rate relative to Eddington, both the black hole spin and the viewing angle can significantly shift the observed peak energy of the spectrum, particularly for rapid spin viewed obliquely or edge-on. We present a detailed test of the approximations made in various forms of the MCD model. Linear limb-darkening is confirmed to be a reasonable approximation for the integrated flux, but not for many specific frequencies of interest.

The emission line luminosity function of active galactic nuclei (AGN) is measured from about 3000 AGN included in the main galaxy sample of the Sloan Digital Sky Survey within a redshift range of $0<z<0.15$. The $\Ha$ and [OIII]$\lambda 5007$ luminosity functions for Seyferts cover luminosity range of $10^{5-9}$$L_\odot$ in H$\alpha$ and the shapes are well fit by broken power laws, without a turnover at fainter nuclear luminosities. Assuming a universal conversion from emission line strength to continuum luminosity, the inferred B band magnitude luminosity function is comparable both to the AGN luminosity function of previous studies and to the low redshift quasar luminosity function derived from the 2dF redshift survey. The inferred AGN number density is approximately 1/5 of all galaxies and about $6\times 10^{-3}$ of the total light of galaxies in the $r$-band comes from the nuclear activity. The numbers of Seyfert 1s and Seyfert 2s are comparable at low luminosity, while at high luminosity, Seyfert 1s outnumber Seyfert 2s by a factor of 2-4. In making the luminosity function measurements, we assumed that the nuclear luminosity is independent of the host galaxy luminosity, an assumption we test {\it a posteriori}, and show to be consistent with the data. Given the relationship between black hole mass and host galaxy bulge luminosity, the lack of correlation between nuclear and host luminosity suggests that the main variable that determines the AGN luminosity is the Eddington ratio, not the black hole mass. This appears to be different from luminous quasars, which are most likely to be shining near the Eddington limit.

We have compiled a large sample of low-redshift active galactic nuclei (AGN) identified via their emission line characteristics from the spectroscopic data of the Sloan Digital Sky Survey. Since emission lines are often contaminated by stellar absorption lines, we developed an objective and efficient method of subtracting the stellar continuum from every galaxy spectrum before making emission line measurements. The distribution of the measured H$\alpha$ Full Width at Half Maxima values of emission line galaxies is strongly bimodal, with two populations separated at about 1,200km s$^{-1}$. This feature provides a natural separation between narrow-line and broad-line AGN. The narrow-line AGN are identified using standard emission line ratio diagnostic diagrams. 1,317 broad-line and 3,074 narrow-line AGN are identified from about 100,000 galaxy spectra selected over 1151 square degrees. This sample is used in a companion paper to determine the emission-line luminosity function of AGN.

We have conducted spectropolarimetry of 12 type II (obscured) quasar candidates selected from the spectroscopic database of the Sloan Digital Sky Survey based on their emission line properties. Polarization was detected in all objects, with nine being highly polarized (> 3%) and with polarization reaching as high as 17% in two objects. Broad lines were detected in the polarized spectra of five objects. These observations prove beyond a reasonable doubt that the objects in our sample are indeed type II quasars, in that they harbor luminous UV-excess AGNs in their centers and that the direct view to the AGN is highly obscured. For three of the objects in this paper, we have obtained HST images in three bands. The HST observations, combined with the spectropolarimetry data, imply that scattering off material outside the obscuration plane is the dominant polarization mechanism. In all three objects the sizes of scattering regions are a few kpc. For one object, the extent of the scattering region, coupled with the characteristics of the polarized spectrum, argue strongly that dust scattering rather than electron scattering dominates the polarized light. Our observations are well-described by the basic orientation-based unification model of toroidal obscuration and off-plane scattering, implying that the model can be extended to include at least some high-luminosity AGNs.

This paper continues the analysis of a set of general relativistic 3D MHD simulations of accreting tori in the Kerr metric with different black hole spins. We focus on bound matter inside the initial pressure maximum, where the time-averaged motion of gas is inward and an accretion disk forms. We use the flows of mass, angular momentum, and energy in order to understand dynamics in this region. The sharp reduction in accretion rate with increasing black hole spin reported in Paper I of this series is explained by a strongly spin-dependent outward flux of angular momentum conveyed electromagnetically; when a/M > 0.9, this flux can be comparable to the inward angular momentum flux carried by the matter. In all cases, there is outward electromagnetic angular momentum flux throughout the flow; in other words, contrary to the assertions of traditional accretion disk theory, there is in general no "stress edge", no surface within which the stress is zero. The retardation of accretion in the inner disk by electromagnetic torques also alters the radial distribution of surface density, an effect that may have consequences for observable properties such as Compton reflection. The net accreted angular momentum is sufficiently depressed by electromagnetic effects that in the most rapidly-spinning black holes mass growth can lead to spindown. Spinning black holes also lose energy by Poynting flux; this rate is also a strongly increasing function of black hole spin, rising to 10% or more of the rest-mass accretion rate at very high spin. As the black hole spins faster, the path of the Poynting flux changes from being predominantly within the accretion disk to predominantly within the funnel outflow.

The origin and nature of Ultra-Luminous X-ray sources (ULXs) is a contentious and controversial topic. There are ongoing debates about the masses of the objects responsible, their sources of mass for accretion, and their relation to stellar populations in galaxies. A new picture of these objects is proposed in which they are intermediate-mass black holes confined to the disks of their host galaxies and accreting from the interstellar medium. They are then preferentially found in or near molecular clouds. This model correctly predicts the shape of the observed luminosity function and requires only a very small fraction of the baryonic mass of a galaxy to be in the form of intermediate-mass black holes. Because the X-rays they produce strongly heat nearby interstellar gas and because they move relatively rapidly in and out of dense regions, ULXs are predicted to have brief episodes of high luminosity, perhaps of order 10^5 yr in duration, but they may recur many times.

We have carried out fully relativistic numerical simulations of accretion disks in the Kerr metric. In this paper we focus on the unbound outflows that emerge self-consistently from the accretion flow. These outflows are found in the axial funnel region and consist of two components: a hot, fast, tenuous outflow in the axial funnel proper, and a colder, slower, denser jet along the funnel wall. Although a rotating black hole is not required to produce these unbound outflows, their strength is enhanced by black hole spin. The funnel-wall jet is excluded from the axial funnel due to elevated angular momentum, and is also pressure-confined by a magnetized corona. The tenuous funnel outflow accounts for a significant fraction of the energy transported to large distances in the higher-spin simulations. We compare the outflows observed in our simulations with those seen in other simulations.

Type II quasars are luminous AGNs whose central engines and broad-line regions are obscured by intervening material; such objects only recently have been discovered in appreciable numbers. We study the multiwavelength properties of 291 type II AGN candidates (0.3 < z < 0.8) selected based on their optical emission line properties from the spectroscopic database of the Sloan Digital Sky Survey. This sample includes about 150 objects luminous enough to be classified as type II quasars. We matched the sample to the FIRST (20 cm), IRAS (12-100 micron), 2MASS (JHK_S) and RASS (0.1-2.4 keV) surveys. Roughly 10% of optically selected type II AGN candidates are radio-loud, comparable to the AGN population as a whole. About 40 objects are detected by IRAS at 60 micron and/or 100 micron, and the inferred mid/far-IR luminosities lie in the range nu L_nu=10^45-3x10^46 erg/sec. Average IR-to-[OIII]5007 ratios of objects in our sample are consistent with those of other AGNs. Objects from our sample are ten times less likely to have soft X-ray counterparts in RASS than type I AGNs with the same redshifts and [OIII]5007 luminosities. The few type II AGN candidates from our sample that are detected by RASS have harder X-ray spectra than those of type I AGNs. The multiwavelength properties of the type II AGN candidates from our sample are consistent with their interpretation as powerful obscured AGNs.

We calculate the efficiency of iron K line emission and iron K absorption in photoionized models using a new set of atomic data. These data are more comprehensive than those previously applied to the modeling of iron K lines from photoionized gases, and allow us to systematically examine the behavior of the properties of line emission and absorption as a function of the ionization parameter, density and column density of model constant density clouds. We show that, for example, the net fluorescence yield for the highly charged ions is sensitive to the level population distribution produced by photoionization, and these yields are generally smaller than those predicted assuming the population is according to statistical weight. We demonstrate that the effects of the many strongly damped resonances below the K ionization thresholds conspire to smear the edge, thereby potentially affecting the astrophysical interpretation of absorption features in the 7-9 keV energy band. We show that the centroid of the ensemble of K$\alpha$ lines, the K$\beta$ energy, and the ratio of the K$\alpha_1$ to K$\alpha_2$ components are all diagnostics of the ionization parameter of our model slabs

We report results from {\it XMM-Newton} observations of thirteen X-ray bright BL Lacertae objects, selected from the {\it Einstein} Slew Survey sample. The spectra are generally well fit by power-law models, with four objects having hard ($\alpha<1; F_\nu \propto \nu^{-\alpha}$) spectra that indicates synchrotron peaks at $>5$ keV. None of our spectra show line features, indicating that soft X-ray absorption ``notches'' must be rare amongst BL Lacs, rather than common or ubiquitous as had previously been asserted. We find significant curvature in most of the spectra. This curvature is almost certainly intrinsic, as it appears nearly constant from 0.5 to 6 keV, an observation which is inconsistent with the small columns seen in these sources.

A selection of results from the general relativistic MHD accretion simulations described in the previous talk are presented. We find that the magnetic field strength increases sharply with decreasing radius and is also enhanced near rapidly-spinning black holes. The greater magnetic field strength associated with rapid black hole rotation leads to a large outward electromagnetic angular momentum flux that substantially reduces both the mean accretion rate and the net accreted angular momentum per unit rest-mass. This electromagnetic stress strongly violates the traditional guess that the accretion stress vanishes at and inside the marginally stable orbit. Possible observational consequences include a constraint on the maximum spin of black holes, enhancement to the radiative efficiency, and concentration of fluorescent Fe Kalpha to the innermost part of the accretion disk.

We present a detailed analysis of the magnetic field structure found in general relativistic 3D MHD simulations of accreting tori in the Kerr metric with different black hole spins. Among the properties analyzed are the field strength as a function of position and black hole spin, the shapes of field lines, the degree to which they connect different regions, and their degree of tangling. We investigate prior speculations about the structure of the magnetic fields and discuss how frequently certain configurations are seen in the simulations. We also analyze the distribution of current density, with a view toward identifying possible locations for magnetic energy dissipation.

Type II quasars are suggested by unification models of Active Galactic Nuclei as luminous analogs of Seyfert 2 galaxies. A number of different methods have been used hitherto to discover type II quasars, but only a handful have been found. We selected about 150 type II quasar candidates from the SDSS spectroscopic database; these objects have strong narrow high-ionization emission lines. We describe the selection procedure and estimate intrinsic luminosities of these optically obscured AGN. A campaign to perform a multi-wavelength follow-up of the sample is now underway.

Type II quasars are the long-sought luminous analogs of type II (narrow emission line) Seyfert galaxies, suggested by unification models of active galactic nuclei (AGN) and postulated to account for an appreciable fraction of the cosmic hard X-ray background. We present a sample of 291 type II AGN at redshifts 0.3<Z<0.83 from the spectroscopic data of the Sloan Digital Sky Survey. These objects have narrow (FWHM<2000 km/s), high equivalent width emission lines with high-ionization line ratios. We describe the selection procedure and discuss the optical properties of the sample. About 50% of the objects have [OIII] \lambda_air 5007 line luminosities in the range 3\times 10^8-10^10 L_Sun, comparable to those of luminous (-27<M_B<-23) quasars; this, along with other evidence, suggests that the objects in the luminous subsample are type II quasars.

This is the first in a series of papers that investigate the properties of accretion flows in the Kerr metric through three-dimensional, general relativistic magnetohydrodynamic simulations of tori with a near-Keplerian initial angular velocity profile. We study four models with increasing black hole spin, from a/M=0 to 0.998, for which the structural parameters of the initial tori are maintained nearly constant. The subsequent accretion flows arise self-consistently from stresses and turbulence created by the magnetorotational instability. We investigate the overall evolution and the late-time global structure in the resulting non-radiative accretion flows, including the magnetic fields within the disks, the properties of the flow in the plunging region, and the flux of conserved quantities into the black hole. Independent of black hole spin, the global structure is described in terms of five regions: the main disk body, the coronal envelope, the inner disk, consisting of an inner torus and plunging region, an evacuated axial funnel, and a bi-conical outflow confined to the corona-funnel boundary. We find evidence for lower accretion rates, stronger funnel-wall outflows, and increased stress in the near hole region with increasing black hole spin.

We examine the small-scale dynamics of black hole accretion disks in which radiation pressure exceeds gas pressure. Local patches of disk are modeled by numerically integrating the equations of radiation MHD in the flux-limited diffusion approximation. The shearing-box approximation is used, and the vertical component of gravity is neglected. Magneto-rotational instability (MRI) leads to turbulence in which accretion stresses are due primarily to magnetic torques. When radiation is locked to gas over the length and time scales of fluctuations in the turbulence, the accretion stress, density contrast, and dissipation differ little from those in the corresponding calculations with radiation replaced by extra gas pressure. However, when radiation diffuses each orbit a distance that is comparable to the RMS vertical wavelength of the MRI, radiation pressure is less effective in resisting squeezing. Large density fluctuations occur, and radiation damping of compressive motions converts PdV work into photon energy. The accretion stress in calculations having a net vertical magnetic field is found to be independent of opacity over the range explored, and approximately proportional to the square of the net field. In calculations with zero net magnetic flux, the accretion stress depends on the portion of the total pressure that is effective in resisting compression. The stress is lower when radiation diffuses rapidly with respect to the gas. We show that radiation-supported Shakura-Sunyaev disks accreting via internal magnetic stresses are likely in their interiors to have radiation marginally coupled to turbulent gas motions.

We present a catalog of 224 broad absorption line quasars (BALQSOs) from the Sloan Digital Sky Survey's Early Data Release Quasar Catalog, including a relatively complete and homogeneous subsample of 131 BALQSOs. Since the identification of BALQSOs is subject to considerable systematic uncertainties, we attempt to create a complete sample of SDSS BALQSOs by combining the results of two automated selection algorithms and a by-eye classification scheme. One of these automated algorithms finds broad absorption line troughs by comparison with a composite quasar spectrum. We present the details of this algorithm and compare this method to that which uses a power-law fit to the continuum. The BALQSOs in our sample are further classified as high-ionization BALQSOs (HiBALs), low-ionization BALQSOs (LoBALs), and BALQSOs with excited iron absorption features (FeLoBALs); composite spectra of each type are presented. We further present a study of the properties of the BALQSOs in terms of the balnicity distribution, which rises with decreasing balnicity. This distribution of balnicities suggests that the fraction of quasars with intrinsic outflows may be significantly underestimated.

A new broad absorption line quasar (BAL) sample is derived from the first data released by the Sloan Digital Sky Survey. With 116 objects, it is the largest BAL sample yet assembled. Over the redshift range 1.8 < z < 3.8, the crude fraction with broad absorption in the CIV line is roughly 15%. This fraction may be subject to small selection-efficiency adjustments. There are also hints of redshift-dependence in the BAL fraction. The sample is large enough to permit the first estimate of the distribution of "balnicity index": subject to certain arbitrary parameters in the definition of this quantity, it is very broad, with (roughly) equal numbers of objects per logarithmic interval of balnicity. BAL quasars are also found to be redder on average than non-BAL quasars. The fraction of radio-loud BAL quasars is (weakly) consistent with the fraction of radio-loud ordinary quasars.

(abridged) Based on observations of the Seyfert nucleus in NGC1068 with ASCA, RXTE and BeppoSAX, we report the discovery of a flare (increase in flux by a factor of ~1.6) in the 6.7 keV Fe K line component between observations obtained 4 months apart, with no significant change in the other (6.21, 6.4, and 6.97 keV) Fe K_alpha line components. During this time, the continuum flux decreased by ~20%. The RXTE spectrum requires an Fe K absorption edge near 8.6 keV (Fe XXIII - XXV). The spectral data indicate that the 2-10 keV continuum emission is dominated (~2/3 of the luminosity) by reflection from a previously unidentified region of warm, ionized gas located <~ 0.2 pc from the AGN. The remaining ~1/3 of the observed X-ray emission is reflected from optically thick, neutral gas. The inferred properties of the warm reflector (WR) are: size (diameter) <~0.2 pc, gas density n >~ 10^{5.5} /cm3, ionization parameter xi approx 10^{3.5} erg cm/s, and covering fraction 0.003 (L_0/10^{43.5} erg/s)^{-1} < (Omega/4 pi) < 0.024 (L_0/10^{43.5})^{-1}, where L_0 is the intrinsic 2-10 keV X-ray luminosity of the AGN. We suggest that the WR gas is the source of the (variable) 6.7 keV Fe line emission, and the 6.97 keV Fe line emission. The 6.7 keV line flare is assumed to be due to an increase in the emissivity of the WR gas from a decrease (by 20-30%) in L_0. The properties of the WR are most consistent with an intrinsically X-ray weak AGN with L_0 approx 10^{43.0} erg/s. The optical and UV emission that scatters from the WR into our line of sight is required to suffer strong extinction, which can be reconciled if the line-of-sight skims the outer surface of the torus. Thermal bremsstrahlung radio emission from the WR may be detectable in VLBA radio maps of the NGC 1068 nucleus.

We analyze X-ray spectra of heavily obscured (N_H > 10^{24} cm^{-2}) active galaxies obtained with Chandra, concentrating on the iron K alpha fluorescence line. We measure very large equivalent widths in most cases, up to 5 keV in the most extreme example. The geometry of an obscuring torus of material near the active galactic nucleus (AGN) determines the Fe emission, which we model as a function of torus opening angle, viewing angle, and optical depth. The starburst/AGN composite galaxies in this sample require small opening angles. Starburst/AGN composite galaxies in general therefore present few direct lines of sight to their central engines. These composite galaxies are common, and their large covering fractions and heavy obscuration effectively hide their intrinsically bright X-ray continua. While few distant obscured AGNs have been identified, we propose to exploit their signature large Fe K alpha equivalent widths to find more examples in X-ray surveys.

Although soft X-ray absorption features in AGN were discovered almost ten years ago, the nature and location of the gas creating them has remained controversial. However, by making use of the newly-available high-resolution spectra provided by XMM-Newton and Chandra, we should be able to make substantial advances. The first such spectra indicate that multiple ionization states coexist in the absorber; this is a natural consequence of photoionization physics. Photoionized evaporation in the presence of a copious mass source locks the ratio of ionizing intensity to pressure to a critical value. A broad range of temperatures can all coexist in equilibrium for this value of the ratio of ionizing intensity to pressure. Consequently, the flow is expected to be strongly inhomogeneous in temperature. The inferred distance of this material from the source of ionizing radiation depends on how much matter exists at the highest-obtainable temperature. This distance can be measured by monitoring how ionic column densities respond to changes in the ionizing continuum on timescales of days to years.

What is meant by the "inner edge" of an accretion disk around a black hole depends on the property that defines the edge. We discuss four such definitions using data from recent high-resolution numerical simulations. These are: the "turbulence edge", where flux-freezing becomes more important than turbulence in determining the magnetic field structure; the "stress edge", where plunging matter loses dynamical contact with the outer accretion flow; the "reflection edge", the smallest radius capable of producing significant X-ray reflection features; and the "radiation edge", the innermost place from which significant luminosity emerges. All these edges are dependent on the accretion rate and are non-axisymmetric and time-variable. Although all are generally located in the vicinity of the marginally stable orbit, significant displacements can occur, and data interpretations placing the disk edge precisely at this point can be misleading. If observations are to be used successfully as diagnostics of accretion in strong gravity, the models used to interpret them must take careful account of these distinctions.

We present new VLT spectroscopic observations of the most distant quasar known, SDSS J1030+0524 at z=6.28 which was recently discovered by the Sloan Digital Sky Survey. We confirm the presence of a complete Gunn-Peterson trough caused by neutral hydrogen in the intergalactic medium. There is no detectable flux over the wavelength range from 8450 to 8710 A. We set an improved limit on the drop of the flux level blueward of the Ly alpha line: a factor of > 200. Below 8450 A the spectrum shows a rise in flux, with a large fraction (> 60 %) of the total emission produced by few narrow features of transmitted flux. We discuss the obvious proximity effect around this quasar, with the presence of transmitted flux with many absorption features in a region of about 23h^{-1} comoving Mpc. If assuming the surrounding medium were completely neutral, the size of this region would imply a quasar lifetime of \~1.3x10^7 years. We also present near-IR spectroscopy of both SDSS J1030+0524 and of SDSS J1306+05, the second most distant quasar known at redshift 6.0. We combine measurements of the CIV line and limits on the HeII emission with the NV line measurements from the optical spectra to derive line ratios, and by implication the abundances of these early quasar environments. The results are indistinguishable from those of lower redshift quasars and indicate little or no evolution in the abundances from z ~ 6 to z ~ 2. The line ratios suggest supersolar metallicities, implying that the first stars around the quasars must have formed at least a few hundreds of Myrs prior to the observation, i.e. at redshift higher than 8.

We present the results of UV imaging polarimetry of the Seyfert 2 galaxy Mrk 477 taken by the Faint Object Camera onboard the Hubble Space Telescope (HST). From a previous HST UV image (lambda ~ 2180A), Mrk 477 has been known to have a pointlike bright UV hotspot in the central region, peculiar among nearby Seyfert 2 galaxies. There are also claims of UV/optical variability, unusual for a Seyfert 2 galaxy. Our data show that there is an off-nuclear scattering region ~ 0."6 (~ 500 pc) NE from the hotspot. The data, after the subtraction of the instrumental effect due to this bright hotspot region, might indicate that the scattered light is also detected in the central 0."2 radius region and is extended to a very wide angle. The hotspot location is consistent with the symmetry center of the PA pattern, which represents the location of the hidden nucleus, but our data do not provide a strong upper limit to the distance between the symmetry center and the hotspot. We have obtained high spatial resolution color map of the continuum which shows that the nuclear spiral arm of 0."4 scale (~ 300pc) is significantly bluer than the off-nuclear mirror and the hotspot region. The nature of the hotspot is briefly discussed.

We present UV imaging polarimetry data of the Seyfert 2 galaxy Mrk 3 taken by the Hubble Space Telescope. The polarized flux is found to be extended to ~1 kpc from the nucleus, and the position angles of polarization are centrosymmetric, confirming that the polarization is caused by scattering. We determine the location of the hidden nucleus as the center of this centrosymmetric pattern. From the polarization images taken in two broad bands, we have obtained the color distribution of the polarized flux. Some regions have blue polarized flux, consistent with optically-thin dust scattering, but some bright knots have a color similar to that of Seyfert 1 nucleus. Also, the recent Chandra X-ray observation suggests that the ratio of scattered UV flux to scattered X-ray flux is rather similar to the intrinsic UV/X-ray ratio in a Seyfert 1 nucleus, if the observed extended X-ray continuum is scattered light. While the scattered X-ray would be essentially from electron scattering, the UV slope and UV/X-ray ratio both being similar to Seyfert 1's would lead to two possibilities as to the nature of the UV scatterers. One is that the UV may also be scattered by electrons, in which case the scattering gas is somehow dust-free. The other is that the UV is scattered by dust grains, but the wavelength-independent UV scattering with low efficiency indicated by the UV slope and UV/X-ray ratio would suggest that the grains reside in UV-opaque clouds, or the dust might be mainly composed of large grains and lacks small-grain population.

New three dimensional magnetohydrodynamic simulations of accretion disk dynamics in a pseudo-Newtonian Paczynski-Wiita potential are presented. These have finer resolution in the inner disk than any previously reported. Finer resolution leads to increased magnetic field strength, greater accretion rate, and greater fluctuations in the accretion rate. One simulation begins with a purely poloidal magnetic field, the other with a purely toroidal field. Compared to the poloidal initial field simulation, a purely toroidal initial field takes longer to reach saturation of the magnetorotational instability and produces less turbulence and weaker magnetic field energies. For both initial field configurations, magnetic stresses continue across the marginally stable orbit; measured in units corresponding to the Shakura-Sunyaev alpha parameter, the stress grows from ~0.1 in the disk body to as much as ~10 deep in the plunging region. Matter passing the inner boundary of the simulation has ~10% greater binding energy and ~10% smaller angular momentum than it did at the marginally stable orbit. Both the mass accretion rate and the integrated stress fluctuate widely on a broad range of timescales.

We fit general relativistic, geometrically thin accretion disk models with non-LTE atmospheres to near simultaneous multiwavelength data of 3C~273, extending from the optical to the far ultraviolet. Our model fits show no flux discontinuity associated with a hydrogen Lyman edge, but they do exhibit a spectral break which qualitatively resembles that seen in the data. This break arises from relativistic smearing of Lyman emission edges which are produced locally at tens of gravitational radii in the disk. We discuss the possible effects of metal line blanketing on the model spectra, as well as the substantial Comptonization required to explain the observed soft X-ray excess. Our best fit accretion disk model underpredicts the near ultraviolet emission in this source, and also has an optical spectrum which is too red. We discuss some of the remaining physical uncertainties, and suggest in particular that an extension of our models to the slim disk regime and/or including nonzero magnetic torques across the innermost stable circular orbit may help resolve these discrepancies.

Although soft X-ray absorption features in AGN were discovered almost ten years ago, the nature and location of the gas creating them has remained controversial. Guided by the results of recent high-resolution X-ray spectroscopy, we argue that these features are created in a multi-temperature wind whose source of matter is photoionized evaporation from the inner edge of the obscuring torus often found surrounding AGN. Photoionized evaporation in the presence of a copious mass source locks the ratio of ionizing intensity to pressure to a critical value. However, a broad range of temperatures can all coexist in equilibrium for this value of the ratio of ionizing intensity to pressure. Consequently, the flow is expected to be strongly inhomogeneous in temperature. The inferred distance of this material from the source of ionizing radiation depends on how much matter exists at the highest-obtainable temperature. This distance can be measured by monitoring how ionic column densities respond to changes in the ionizing continuum on timescales of days to years.

The Sloan Digital Sky Survey has discovered a population of broad absorption line quasars with various extreme properties. Many show absorption from metastable states of FeII with varying excitations; several objects are almost completely absorbed bluewards of MgII; at least one shows stronger absorption from FeIII than FeII, indicating temperatures T>35000 K in the absorbing region; and one object even seems to have broad H-beta absorption. Many of these extreme BALs are also heavily reddened, though `normal' BALs (particularly LoBALs) from SDSS also show evidence for internal reddening.

We present 13 Broad Absorption Line (BAL) quasars, including 12 new objects, which were identified in the Sloan Digital Sky Survey (SDSS) and matched within 2'' to sources in the FIRST radio survey catalog. The surface density of this sample of radio-detected BAL quasars is 4.5 +- 1.2 per 100 deg^2, i.e. approximately 4 times larger than previously found by the shallower FIRST Bright Quasar Survey (FBQS). A majority of these radio-detected BAL quasars are moderately radio-loud objects. The fraction of BAL quasars in the entire radio quasar sample, 4.8 +- 1.3 %, is comparable to the fraction of BAL quasars among the SDSS optical quasar sample (ignoring selection effects). We estimate that the true fraction of BAL quasars (mostly HiBALs) in the radio sample is 9.2 +- 2.6 % once selection effects are accounted for. We caution that the absorption troughs of 4 of the 13 radio-detected quasars considered do not strictly satisfy the standard BALnicity criterion. One or possibly two of the new radio-detected BAL quasars are of the rare ``iron LoBAL'' type. BAL quasars are generally redder than the median SDSS quasar at the same redshift.

We show that if all gamma-ray bursts emit X-rays in a way similar to those observed by BeppoSax, much of the extinction along the line of sight in the host galaxy of the burst can be destroyed. Two mechanisms are principally responsible for dust destruction: grain heating and grain charging. The latter, which can lead to electrostatic stresses greater than the tensile strength of the grains, is often the more important. Grains may regularly be destroyed at distances as large as 100 parsecs. This dust destruction can permit us to see the UV/optical afterglow even when the burst is embedded deep within a highly-obscured region. Because the destruction rate depends on grain composition and size, it may be possible to observe the amount and wavelength-dependence of extinction change during the course of the burst and first few minutes of the afterglow. It may also be possible to detect interstellar absorption lines in the afterglow spectrum that would not exist but for the return of heavy elements to the gas phase.

We extend our models of the vertical structure and emergent radiation field of accretion disks around supermassive black holes described in previous papers of this series. Our models now include both a self-consistent treatment of Compton scattering and the effects of continuum opacities of the most important metal species (C, N, O, Ne, Mg, Si, S, Ar, Ca, Fe, Ni). With these new effects incorporated, we compute the predicted spectrum from black holes accreting at nearly the Eddington luminosity (L/L_Edd = 0.3) and central masses of 10^6, 10^7, and 10^8 M_sun. We also consider two values of the Shakura-Sunyaev alpha parameter, 0.1 and 0.01. Although it has little effect when M > 10^8 M_sun, Comptonization grows in importance as the central mass decreases and the central temperature rises. It generally produces an increase in temperature with height in the uppermost layers of hot atmospheres. Compared to models with coherent electron scattering, Comptonized models have enhanced EUV/soft X-ray emission, but they also have a more sharply declining spectrum at very high frequencies. Comptonization also smears the hydrogen and the He II Lyman edges. The effects of metals on the overall spectral energy distribution are smaller than the effects of Comptonization for these parameters. Compared to pure hydrogen-helium models, models with metal continuum opacities have reduced flux in the high frequency tail, except at the highest frequencies, where the flux is very low. Metal photoionization edges are not present in the overall disk-integrated model spectra. In addition to our new grid of models, we also present a simple analytic prescription for the vertical temperature structure of the disk in the presence of Comptonization, and show under what conditions a hot outer layer (a corona) is formed.

The standard equilibrium for radiation-dominated accretion disks has long been known to be viscously, thermally, and convectively unstable, but the nonlinear development of these instabilities---hence the actual state of such disks---has not yet been identified. By performing local two-dimensional hydrodynamic simulations of disks, we demonstrate that convective motions can release heat sufficiently rapidly as to substantially alter the vertical structure of the disk. If the dissipation rate within a vertical column is proportional to its mass, the disk settles into a new configuration thinner by a factor of two than the standard radiation-supported equilibrium. If, on the other hand, the vertically-integrated dissipation rate is proportional to the vertically-integrated total pressure, the disk is subject to the well-known thermal instability. Convection, however, biases the development of this instability toward collapse. The end result of such a collapse is a gas pressure-dominated equilibrium at the original column density.

For the past twenty-five years, nearly all analyses of accretion disk dynamics have assumed that stress inside the disk is locally proportional to pressure (the "alpha-model") and that this stress goes to zero at the marginally stable orbit. Recently, it has been demonstrated that MHD turbulence accounts for the bulk of internal disk stress. In contradiction with the traditional view, the stress from this MHD turbulence does not diminish near the marginally stable orbit, and the ratio of magnetic stress to pressure rises sharply there. Examples of the consequences include: an increase in accretion efficiency that may also be time- and circumstance-dependent; a decrease in the rate of black hole spin-up by accretion; and generation of disk luminosity fluctuations. Preliminary results from numerical simulations lend support to analytic estimates of these effects.

Mid-infrared images of the Seyfert 1 galaxy Mrk 279 obtained with the ISO satellite are presented together with the results of a one-year monitoring campaign of the 2.5-11.7 micron spectrum. Contemporaneous optical photometric and spectrophotometric observations are also presented. The galaxy appears as a point-like source at the resolution of the ISOCAM instrument, 4-5". The 2.5-11.7 micron average spectrum of the nucleus in Mrk 279 shows a strong power law continuum with a spectral index alpha = -0.80+/-0.05 and weak PAH emission features. The Mrk 279 spectral energy distribution shows a mid-IR bump, which extends from 2 to 15-20 micron . The mid-IR bump is consistent with thermal emission from dust grains at a distance of >= 100 light-days. No significant variations of the mid-IR flux have been detected during our observing campaign, consistent with the relatively low amplitude (~10 % rms) of the optical variability during the campaign. The time delay for the Hbeta line emission in response to the optical continuum variations is 16.7 +5.3/5.6 days, consistent with previous measurements.

The mass of the central black hole in many active galactic nuclei has been estimated on the basis of the assumption that the dynamics of the broad emission line gas are dominated by the gravity of the black hole. The most commonly-employed method is to estimate a characteristic size-scale $r_*$ from reverberation mapping experiments and combine it with a characteristic velocity $v_*$ taken from the line profiles; the inferred mass is then estimated by $r_* v_*^2/G$. We critically discuss the evidence supporting the assumption of gravitational dynamics and find that the arguments are still inconclusive. We then explore the range of possible systematic error if the assumption of gravitational dynamics is granted. Inclination relative to a flattened system may cause a systematic underestimate of the central mass by a factor $\sim (h/r)^2$, where $h/r$ is the aspect ratio of the flattening. The coupled effects of a broad radial emissivity distribution, an unknown angular radiation pattern of line emission, and sub-optimal sampling in the reverberation experiment can cause additional systematic errors as large as a factor of 3 or more in either direction.

We present a detailed three dimensional magnetohydrodynamic (MHD) simulation describing the inner region of a disk accreting onto a black hole. To avoid the technical complications of general relativity, the dynamics are treated in Newtonian fashion using the pseudo-Newtonian Pacz\'ynski-Wiita potential. The disk evolves due to angular momentum transport which is produced naturally from MHD turbulence generated by the magnetorotational instability. We find that the resulting stress is continuous across the marginally stable orbit, in contradiction with the widely-held assumption that the stress should go to zero there. As a consequence, the specific angular momentum of the matter accreted into the hole is smaller than the specific angular momentum at the marginally stable orbit. The disk exhibits large fluctuations in almost every quantity, both spatially and temporally. In particular, the ratio of stress to pressure (the local analog of the Shakura-Sunyaev $\alpha$ parameter) exhibits both systematic gradients and large fluctuations; from $\sim 10^{-2}$ in the disk midplane at large radius, it rises to $\sim 10$ both at a few gas density scaleheights above the plane at large radius, and near the midplane well inside the plunging region. Driven in part by large-amplitude waves excited near the marginally stable orbit, both the mass accretion rate and the integrated stress exhibit large fluctuations whose Fourier power spectra are smooth "red" power-laws stretching over several orders of magnitude in timescale.

Numerous variations have been proposed on the original suggestion by Blandford and Znajek that magnetic fields could be used to extract rotational energy from black holes. A new categorization of these variations is proposed so that they may be considered in a systematic way. "Black hole spindown" is defined more precisely, distinguishing decrease in the spin parameter a/M from decreases in angular momentum a and rotational kinetic energy, M - M_i. Several key physical questions are raised: Can the "stretched horizon" of a black hole communicate with the outside world? Do accretion disks bring any net magnetic flux to the black holes at their centers? Is the magnetic field adjacent to a black hole force-free everywhere?

We study radiation-hydrodynamical normal modes of radiation-supported accretion disks in the WKB limit. It has long been known that in the large optical depth limit the standard equilibrium is unstable to convection. We study how the growth rate depends on location within the disk, optical depth, disk rotation, and the way in which the local dissipation rate depends on density and pressure. The greatest growth rates are found near the disk surface. Rotation stabilizes vertical wavevectors, so that growing modes tend to have nearly-horizontal wavevectors. Over the likely range of optical depths, the linear growth rate for convective instability has only a weak dependence on disk opacity. Perturbations to the dissipation have little effect on convective mode growth rates, but can cause growth of radiation sound waves.

During 1998 April 13-16, NGC 3516 was monitored almost continuously with HST for 10.3 hr in the UV and 2.8 d in the optical, and simultaneous RXTE and ASCA monitoring covered the same period. The X-rays were strongly variable with the soft (0.5-2 keV) showing stronger variations (~65% peak-to-peak) than the hard (2-10 keV; ~50% peak-to-peak). The optical continuum showed much smaller but highly significant variations: a slow ~2.5% rise followed by a faster ~3.5% decline. The short UV observation did not show significant variability. The soft and hard X-ray light curves were strongly correlated with no significant lag. Likewise, the optical continuum bands (3590 and 5510 A) were also strongly correlated with no measurable lag above limits of <0.15 d. However no significant correlation or simple relationship could be found for the optical and X-ray light curves. These results appear difficult to reconcile with previous reports of correlations between X-ray and optical variations and of measurable lags within the optical band for some other Seyfert 1s. These results also present serious problems for "reprocessing" models in which the X-ray source heats a stratified accretion disk which then reemits in the optical/ultraviolet: the synchronous variations within the optical would suggest that the emitting region is <0.3 lt-d across, while the lack of correlation between X-ray and optical variations would indicate, in the context of this model, that any reprocessing region must be >1 lt-d in size. It may be possible to resolve this conflict by invoking anisotropic emission or special geometry, but the most natural explanation appears to be that the bulk of the optical luminosity is generated by some other mechanism than reprocessing.

We have constructed a grid of non-LTE disk models for a wide range of black hole mass and mass accretion rate, for several values of viscosity parameter alpha, and for two extreme values of the black hole spin: the maximum-rotation Kerr black hole, and the Schwarzschild (non-rotating) black hole. Our procedure calculates self-consistently the vertical structure of all disk annuli together with the radiation field, without any approximations imposed on the optical thickness of the disk, and without any ad hoc approximations to the behavior of the radiation intensity. The total spectrum of a disk is computed by summing the spectra of the individual annuli, taking into account the general relativistic transfer function. The grid covers nine values of the black hole mass between M = 1/8 and 32 billion solar masses with a two-fold increase of mass for each subsequent value; and eleven values of the mass accretion rate, each a power of 2 times 1 solar mass/year. The highest value of the accretion rate corresponds to 0.3 Eddington. We show the vertical structure of individual annuli within the set of accretion disk models, along with their local emergent flux, and discuss the internal physical self-consistency of the models. We then present the full disk-integrated spectra, and discuss a number of observationally interesting properties of the models, such as optical/ultraviolet colors, the behavior of the hydrogen Lyman limit region, polarization, and number of ionizing photons. Our calculations are far from definitive in terms of the input physics, but generally we find that our models exhibit rather red optical/UV colors. Flux discontinuities in the region of the hydrogen Lyman limit are only present in cool, low luminosity models, while hotter models exhibit blueshifted changes in spectral slope.

Magnetic connections to the plunging region can exert stresses on the inner edge of an accretion disk around a black hole. We recompute the relativistic corrections to the thin-disk dynamics equations when these stresses take the form of a time-steady torque on the inner edge of the disk. The additional dissipation associated with these stresses is concentrated relatively close outside the marginally stable orbit, scaling as r to the -7/2 at large radius. As a result of these additional stresses: spin-up of the central black hole is retarded; the maximum spin-equilibrium accretion efficiency is 36%, and occurs at a/M=0.94; the disk spectrum is extended toward higher frequencies; line profiles (such as Fe K-alpha) are broadened if the line emissivity scales with local flux; limb-brightening, especially at the higher frequencies, is enhanced; and the returning radiation fraction is substantially increased, up to 58%. This last effect creates possible explanations for both synchronized continuum fluctuations in AGN, and polarization rises shortward of the Lyman edge in quasars. We show that no matter what additional stresses occur, when a/M < 0.36, the second law of black hole dynamics sets an absolute upper bound on the accretion efficiency.

We show how analysis of a quasar high-magnification microlensing event may be used to construct a map of the frequency-dependent surface brightness of the quasar accretion disk. The same procedure also allows determination of the disk inclination angle, the black hole mass (modulo the caustic velocity), and possibly the black hole spin. This method depends on the validity of one assumption: that the optical and ultraviolet continuum of the quasar is produced on the surface of an azimuthally symmetric, flat equatorial disk, whose gas follows prograde circular orbits in a Kerr spacetime (and plunges inside the marginally stable orbit). Given this assumption, we advocate using a variant of first-order linear regularization to invert multi-frequency microlensing lightcurves to obtain the disk surface brightness as a function of radius and frequency. The other parameters can be found by minimizing chi-square in a fashion consistent with the regularized solution for the surface brightness. We present simulations for a disk model appropriate to the Einstein Cross quasar, an object uniquely well-suited to this approach. These simulations confirm that the surface brightness can be reconstructed quite well near its peak, and that there are no systematic errors in determining the other model parameters. We also discuss the observational requirements for successful implementation of this technique.

Qualitative arguments are presented to demonstrate that the energy density of magnetic fields in matter accreting onto a black hole inside the marginally stable orbit is automatically comparable to the rest-mass energy density of the accretion flow. Several consequences follow: magnetic effects must be dynamically significant, but cannot be so strong as to dominate; outward energy transport in Alfven waves may alter the effective efficiency of energy liberation; and vertical magnetic stresses in this region may contribute to "coronal" activity.

As part of an extensive multi-wavelength monitoring campaign, the International Ultraviolet Explorer satellite was used to observe the broad-line radio galaxy 3C 390.3 during the period 1994 December 31 to 1996 March 5. Spectra were obtained every 6-10 days. The UV continuum varied by a factor of 7 through the campaign, while the broad emission-lines varied by factors of 2-5. Unlike previously monitored Seyfert 1 galaxies, in which the X-ray continuum generally varies with a larger amplitude than the UV, in 3C 390.3 the UV continuum light-curve is similar in both amplitude and shape to the X-ray light-curve observed by ROSAT. The UV broad emission-line variability lags that of the UV continuum by 35-70 days for Ly-alpha and CIV 1549; values larger than those found for Seyfert 1 galaxies of comparable UV luminosity. These lags are also larger than those found for the Balmer lines in 3C 390.3 over the same period. The red and blue wings of CIV and Ly-alpha vary in phase, suggesting that radial motion does not dominate the kinematics of the UV line-emitting gas. Comparison with archival data provides evidence for velocity-dependent changes in the Ly-alpha and CIV line profiles, indicating evolution in the detailed properties and/or distribution of the broad-line emitting gas.

We present and discuss HST (WFPC2 and FOC) images and UV GHRS spectra plus ground-based near UV through to near IR spectra of three Seyfert 2 nuclei (NGC 7130, NGC 5135 and IC 3639). These galaxies, together to Mrk 477, were selected from a bigger sample that comprises the 20 brightest Seyfert 2 nuclei, with the goal to study the origin of the UV-optical-near IR featureless continuum in Seyfert 2 nuclei. These four galaxies have bolometric luminosities, as computed with the four IRAS bands, of 10^11 Lsol. They are close enough to be resolved with HST the nuclear zone. This makes these Seyfert 2 galaxies benchmarks to study the Starburst-AGN connection in more distant galaxies. The data provide direct evidence of the existence of a central nuclear starburst that dominates the UV light, and that seem to be responsible for the origin of the so called featureless continuum. These starbursts are dusty and compact. They have sizes (from less than 100 pc to a few hundred pc) much smaller and closer to the nucleus than that seen in the prototype Seyfert 2 galaxy NGC 1068. The bolometric luminosity of these starbursts is similar to the estimated bolometric luminosities of their obscured Seyfert 1 nuclei, and thus they contribute in the same amount to the overall energetics of these galaxies.

MHD turbulence is generally believed to have two important functions in accretion disks: it transports angular momentum outward, and the energy in its shortest wavelength modes is dissipated into the heat that the disks radiate. In this paper we examine a pair of mechanisms which may play an important role in regulating the amplitude and spectrum of this turbulence: photon diffusion and viscosity. We demonstrate that in radiation pressure-dominated disks, photon damping of compressive MHD waves is so rapid that it likely dominates all other dissipation mechanisms.

Accretion disks around black holes in which the shear stress is proportional to the total pressure, the accretion rate is more than a small fraction of Eddington, and the matter is distributed smoothly are both thermally and viscously unstable in their inner portions. The nonlinear endstate of these instabilities is uncertain. Here a new inhomogeneous equilibrium is proposed which is both thermally and viscously stable. In this equilibrium the majority of the mass is in dense clumps, while a minority reaches temperatures $\sim 10^9$ K. The requirements of dynamical and thermal equilibrium completely determine the parameters of this system, and these are found to be in good agreement with the parameters derived from observations of accreting black holes, both in active galactic nuclei and in stellar binary systems.

Spectropolarimetry of Broad Absorption Line Quasars (BALQSO's) has demonstrated that the geometry of the absorbing material is far from spherically symmetric. Calculations of accretion disk spectra suggest that the intrinsic radiation pattern of quasars is also anisotropic. Because the quasar counts distribution is very steep, if the orientations of these two anisotropies are correlated, optical flux-limited samples would be very strongly biassed with respect to discovery of BALQSO's. In particular, currently favored models suggest that BALQSO's may be much more common than in the samples. If so, the intrinsic covering fraction of absorbing matter must be relatively large. Numerous consequences follow, including prediction of a new population of hard X-ray sources (as may have already been seen by ASCA), and a possible explanation of the anomalously strong NV~1240 lines seen in many high redshift quasars.

A critical review is presented of the current status of our understanding of the the spectral output of thermally-radiating accretion disks around black holes, with special attention to those around supermassive black holes in active galactic nuclei. Uncertainties in disk physics (e.g. the nature of the inter-ring stresses, the distribution of dissipation with altitude) contribute to uncertainties in the output. So, too, do inadequacies in calculations (e.g. incomplete treatment of non-LTE effects, neglect of transitions in elements heavier than He, difficulty of coupling Comptonization to genuine transfer). Existing models do not yet adequately explain the absence of Lyman edge features, or the observed flux in the EUV.

We observed the Seyfert 1 galaxy NGC 4151 on eleven occasions at 1-2 day intervals using the Berkeley spectrometer during the ORFEUS-SPAS II mission in 1996 November. The mean spectrum covers 912-1220 A at ~0.3 A resolution with a total exposure of 15,658 seconds. The mean flux at 1000 A was 4.7e-13 erg/cm^2/s/A. We identify the neutral hydrogen absorption with a number of components that correspond to the velocity distribution of \ion{H}{1} seen in our own Galaxy as well as features identified in the CIV 1549 absorption profile by Weymann et al. The main component of neutral hydrogen in NGC 4151 has a total column density of log N_HI = 18.7 +/- 1.5 cm^{-2} for a Doppler parameter b=250 +/- 50 km/s, and it covers 84 +/- 6% of the source. This is consistent with previous results obtained with the Hopkins Ultraviolet Telescope. Other intrinsic far-UV absorption features are not resolved, but the CIII* 1176 absorption line has a significantly higher blueshift relative to NGC 4151 than the CIII 977 resonance line. This implies that the highest velocity region of the outflowing gas has the highest density. Variations in the equivalent width of the CIII* 1176 absorption line anticorrelate with continuum variations on timescales of days. For an ionization timescale <1 day, we set an upper limit of 25 pc on the distance of the absorbing gas from the central source. The OVI 1034 and HeII 1085 emission lines also vary on timescales of 1-2 days, but their response to the continuum variations is complex. For some continuum variations they show no response, while for others the response is instantaneous to the limit of our sampling interval.

We report the detection with the Rossi X-ray Timing Explorer of a Compton reflection signature in the Seyfert galaxy MCG-5-23-16. RXTE also resolves the Fe K-alpha fluorescence line with FWHM ~48,000 km s^{-1}. This measurement provides the first independent confirmation of ASCA detections in Seyfert galaxies of broad Fe K-alpha lines that are thought to be the signature of emission from the inner regions of an accretion disk orbiting a black hole. Under the assumption that reflection arises from an isotropic source located above a neutral accretion disk, and using a theoretical model that accounts for the dependence of the reflected spectrum on inclination angle, we derive a 90% confidence range for the disk inclination of i = 50 to 81 degrees. The large inclination is consistent with that expected from the unified model for MCG-5-23-16 based on its Seyfert 1.9 classification. If we assume that the high-energy cutoff in the incident spectrum lies at energies larger than a few hundred keV, then the equivalent width of the Fe K-alpha line is much larger than predicted for the amount of reflection. This implies either an enhanced iron abundance, a covering factor of reflecting material > 0.5, or a cutoff in the incident spectrum at energies between ~60 and ~200 keV.

We present an improved calculation of the vertical structure and ultraviolet spectrum of a dissipative accretion disk in an AGN. We calculate model spectra in which the viscous stress is proportional to the total pressure, the gas pressure only and the geometric mean of the radiation and gas pressures (cf. Laor & Netzer 1989: LN89). As a result of a more complete treatment of absorptive opacity, we find greater overall spectral curvature than did LN89, as well as larger amplitudes in both the Lyman and HeII photoionization edges. The local black body approximation is not a good description of the near UV spectrum. With relativistic corrections (appropriate to non-rotating black holes) included, we find that the near UV spectrum hardens with increasing m-dot / m_8 (m-dot is the accretion rate in Eddington units, m_8 the black hole mass in units of 10^8 M_Sun). The near UV spectrum is consistent with observations if m-dot/ m_8 \sim 10^{-3}, but disks this cold would have large, and unobserved, absorption features at the Lyman edge. The edge amplitude is reduced when m-dot/m_8 is larger, but then the near-UV slope is too hard to match observations. We conclude that models in which conventional disks orbit non-rotating black holes do not adequately explain UV continuum production in AGN.

Accreting black holes radiate in one of several spectral states, switching from one to another for reasons that are as yet not understood. Using the best studied example, Cyg X-1, we identify the geometry and physical conditions characterizing these states. In particular, we show that in the hard state most of the accretion energy is dissipated in a corona-like structure which fills the inner few tens of gravitational radii around the black hole and has Compton optical depth of order unity. In this state, an optically thick accretion disc extends out to greater distance, but penetrates only a short way into the coronal region. In the soft state, the optically thick disc moves inward and receives the majority of the dissipated energy, while the "corona" becomes optically thin and extends around much of the inner disc. The mass accretion rate in both states is $~10^{-8} M_{\odot}$ yr$^{-1}$.

Using simple analytical estimates we show how the physical parameters characterizing different spectral states of the galactic black hole candidates can be determined using spectral data presently available.

We calculate the flux and spectrum of \gamma-rays emitted by a two-temperature advection-dominated accretion flow (ADAF) around a black hole. The \gamma-rays are from the decay of neutral pions produced through proton-proton collisions. We discuss both thermal and power-law distributions of proton energies and show that the \gamma-ray spectra in the two cases are very different. We apply the calculations to the \gamma-ray source, 2EG J1746-2852, detected by EGRET from the direction of the Galactic Center. We show that the flux and spectrum of this source are consistent with emission from an ADAF around the supermassive accreting black hole Sgr A^* if the proton distribution is a power-law. The model uses accretion parameters within the range made likely by other considerations. If this model is correct, it provides evidence for the presence of a two temperature plasma in Sgr A^*, and predicts \gamma-ray fluxes from other accreting black holes which could be observed with more sensitive detectors.

We present the results of ultraviolet spectral monitoring of the Seyfert galaxy NGC 4151 performed by the Hopkins Ultraviolet Telescope during the Astro-2 mission. We show how the response of the Lyman continuum optical depth to variations in the continuum flux may be used to derive estimates of the distance of the absorbing gas from the central object. In this case, it seems that the absorbing gas is located $\leq 50$~pc away from the central source.

We study the effects of thermal conduction in a hot, active corona above an accretion disk. We assume that all of the dissipative heating takes place in the corona. We find that the importance of conduction decreases with increases in the local dissipative compactness of the corona, l_{diss,loc}, and increases with increasing abundance of electron-positron pairs. For l_{diss,loc} < 1, a significant fraction of the energy released in the corona may be carried away by the conductive flux, leading to formation of a relatively hot transition layer below the base of the corona. Comptonization of disk radiation in such a layer may account for the presence of soft X-ray excesses in the spectra emitted by disk-corona systems.

Motivated by recent work indicating that the UV continuum in AGN may be produced by reradiation of energy absorbed from X-rays irradiating an accretion disk, we present a calculation of the vertical structures and ultraviolet spectra of X-ray irradiated accretion disks around massive non-rotating black holes. After finding the radial dependence of vertically-integrated quantities for these disks, we solve the equations of hydrostatic equilibrium, energy balance, and frequency-dependent radiation transfer as functions of altitude. To solve the last set of equations, we use a variable Eddington factor method. We include electron scattering, free-free, and HI, HeI, and HeII bound-free opacities and the corresponding continuum cooling processes. While the overall spectral shape predicted by X-ray irradiation may be compatible with observations, the Lyman edge emission feature it predicts is not. This finding raises questions for many otherwise plausible models in which X-ray irradiation plays a major role.

A global view is given of how emission lines have been, and may in the future, be used to enhance our understanding of AGN. Lines from the microwave to the X-ray bands all contribute. Although we have a deep understanding of the physical processes by which line photons are generated, when the circumstances of line emission are complicated, models become too unreliable to provide strong inferences about the rest of the AGN system. At present, the lines which appear the most promising for helping answer the most important questions about AGN are the 22~GHz H_2 O rotational transition and the 6.4 keV Fe Kalpha fluorescence.

We observed a total of 13 low-redshift AGN with the Hopkins Ultraviolet Telescope (HUT) during the course of the Astro-1 and Astro-2 missions. Of these, 4 show intrinsic Lyman limit absorption---NGC 1068, NGC 4151, NGC 3516 and NGC 3227. All galaxies with optically thick Lyman limits have extended narrow-line regions with bi-polar morphologies. All AGN with no absorption are Seyfert 1's with compact NLR's. These observations support geometrical shadowing as the means for collimating the ionizing radiation in unified models of AGN, most likely in a photoionized atmosphere above the obscuring torus.

It is generally supposed that when the "compactness" $l \equiv L\sigma_T/(r m_e c^3)$ in photons above the pair-production threshold is large, few $\gamma$-rays can escape. We demonstrate that even when $l \gg 1$, if the high energy and low energy photons are produced in geometrically-separated regions, many of the $\gamma$-rays can, in fact, escape. Pair-production along a thin surface separating the two sources creates enough Compton optical depth to deflect most of the low energy photons away from the high energy ones. Those few low-energy photons which penetrate the shielding surface are reduced in opacity by advection to large distance and small density, by relativistic beaming along the inner edge of the surface, and by Compton upscattering to higher energies inside the surface. The pairs in this surface flow outward relativistically, forming a structure resembling a pair-dominated mildly relativistic jet.

We obtained X-ray spectra of the Seyfert 1 galaxy NGC~3516 in March 1995 using ASCA. Simultaneous far-UV observations were obtained with HUT on the Astro-2 shuttle mission. The ASCA spectrum shows a lightly absorbed power law of energy index 0.78. The low energy absorbing column is significantly less than previously seen. Prominent O~vii and O~viii absorption edges are visible, but, consistent with the much lower total absorbing column, no Fe K absorption edge is detectable. A weak, narrow Fe~K$\alpha$ emission line from cold material is present as well as a broad Fe~K$\alpha$ line. These features are similar to those reported in other Seyfert 1 galaxies. A single warm absorber model provides only an imperfect description of the low energy absorption. In addition to a highly ionized absorber with ionization parameter $U = 1.66$ and a total column density of $1.4 \times 10^{22}~\rm cm^{-2}$, adding a lower ionization absorber with $U = 0.32$ and a total column of $6.9 \times 10^{21}~\rm cm^{-2}$ significantly improves the fit. The contribution of resonant line scattering to our warm absorber models limits the Doppler parameter to $< 160~\rm km~s^{-1}$ at 90\% confidence. Turbulence at the sound speed of the photoionized gas provides the best fit. None of the warm absorber models fit to the X-ray spectrum can match the observed equivalent widths of all the UV absorption lines. Accounting for the X-ray and UV absorption simultaneously requires an absorbing region with a broad range of ionization parameters and column densities.

We observed the Seyfert 1 galaxy NGC 3516 twice during the flight of Astro-2 using the Hopkins Ultraviolet Telescope in March 1995. Simultaneous X-ray observations were performed with ASCA. Our far-ultraviolet spectra cover the spectral range 820-1840 A with a resolution of 2-4 A. No significant variations were found between the two observations. The total spectrum shows a red continuum, $f_\nu \sim \nu^{-1.89}$, with an observed flux of $\rm 2.2 \times 10^{-14}~erg~cm^{-2}~s^{-1}~\AA^{-1}$ at 1450 A, slightly above the historical mean. Intrinsic absorption in Lyman $\beta$ is visible as well as absorption from O~vi 1032,1038, N~v 1239,1243, Si~iv 1394,1403, and C~iv 1548,1551. The UV absorption lines are far weaker than is usual for NGC~3516, and also lie closer to the emission line redshift rather than showing the blueshift typical of these lines when they are strong. The neutral hydrogen absorption, however, is blueshifted by $400~\rm km~s^{-1}$ relative to the systemic velocity, and it is opaque at the Lyman limit. The sharpness of the cutoff indicates a low effective Doppler parameter, $b < \rm 20~km~s^{-1}$. For $b = \rm 10~km~s^{-1}$ the derived intrinsic column is $\rm 3.5 \times 10^{17}~cm^{-2}$. As in NGC~4151, a single warm absorber cannot produce the strong absorption visible over the wide range of observed ionization states. Matching both the UV and X-ray absorption simultaneously requires absorbers spanning a range of $10^3$ in both ionization parameter and column density.

We derive both total flux and velocity-resolved response functions for the CIV 1549 emission line from the data obtained in the 1993 NGC 5548 monitoring campaign. These response functions imply: 1.) the emission region stretches from inside 1 lt-d to outside 10 lt-d, and is probably better described as round than flat; 2.) the velocity field is dominated by a red/blue symmetric component (e.g. 2--d or 3--d random motions, or rotation in a disk) but there is also {\it significant} radial infall. Quantitative modelling indicates that the random speeds are typically a few times as large as the radial speed. However, no simple model gives a completely acceptable fit to the data. These inferences rule out numerous simple and otherwise plausible models for broad line region dynamics, including outflowing winds, radial free-fall, rotation in a disk, or collisionless orbital motion.

The mass of the central black hole in the active galaxy NGC 4258 (M106) has been measured to be $M=3.6\times10^7\Msun$ (Miyoshi et al. 1995). The Eddington luminosity corresponding to this mass is $L_E=4.5\times10^{45}$ erg s$^{-1}$. By contrast the X-ray luminosity of the nucleus of NGC 4258 between $2-10$ keV is $(4\pm 1)\times10^{40}~{\rm erg\,s^{-1}}$ while the optical/UV luminosity is less than $1.5\times10^{42} ~{\rm erg\,s^{-1}}$. The luminosity of NGC 4258 is therefore extremely sub-Eddington, $L\sim10^{-5}L_E$ in X-rays and $L\sim3\times10^{-4} L_E$ even if we take the maximum optical/UV luminosity. Assuming the usual accretion efficiency of 0.1 would imply accretion rates orders of magnitude lower than in Seyfert galaxies and quasars. We show that the properties of the AGN in NGC 4258 can be explained by an accretion flow in the form of a very hot, optically-thin plasma which advects most of the viscously generated thermal energy into the central black hole and radiates only a small fraction of the energy. In this case the accretion rate in Eddington units could be as high as $\sim 0.16\alpha$, where $\alpha$ is the standard viscosity parameter; and the size of the hot disk should be larger than $\sim 10$ times the Schwarzschild radius. We compare the predictions of this model with observations and discuss its application to other low luminosity AGN.

When an AGN is obscured, the warm reflecting gas nearby can be seen by a combination of bremsstrahlung, intrinsic line emission, and reflection of the nuclear continuum, both by electron scattering and by resonance line scattering. Resonance lines, due both to intrinsic emission and scattering, are particularly prominent in the soft X-ray band. When our line of sight to the nucleus is not obscured, the dominant effect is absorption. In the soft X-ray band, ionization edges of highly ionized species and resonance lines contribute comparably to the opacity; in the ultraviolet, the gas is almost transparent except for a small number of resonance lines. We identify the "warm absorbers" seen in many AGN X-ray spectra with this gas, but argue that most of the UV absorption lines seen must be due to a small amount of more weakly ionized gas which is embedded in the main body of the warm, reflecting gas. Because the ionization equilibration timescales of some ions may be as long as the variability timescales in AGN, the ionic abundances indicated by the transmission spectra may not be well-described by ionization equilibrium.

The hard X-ray spectra of accreting black holes systems are generally well-fit by thermal Comptonization models with temperatures $\sim 100$ keV. We demonstrate why, over many orders of magnitude in heating rate and seed photon supply, hot plasmas radiate primarily by inverse Compton scattering, and find equilibrium temperatures within a factor of a few of 100 keV. We also determine quantitatively the (wide) bounds on heating rate and seed photon supply for which this statement is true. Plasmas in thermal balance in this regime obey two simple scaling laws, one relating the product of temperature and optical depth to the ratio of seed photon luminosity to plasma heating rate $l_s/l_h$, the other relating the spectral index of the output power-law to $l_s/l_h$. Because $\alpha$ is almost independent of everything but $l_s/l_h$, the observed power law index may be used to estimate $l_s/l_h$. In both AGN and stellar black holes, the mean value estimated this way is $l_s/l_h \sim 0.1$. As a corollary, $\Theta \tau_T$ must be $\simeq 0.1$ -- 0.2, depending on plasma geometry.