We study the structure of accretion disks around supermassive black holes in the radial range $30\sim 100$ gravitational radii, using a three dimensional radiation magneto-hydrodynamic simulation. For typical conditions in this region of Active Galactic Nuclei (AGN), the Rosseland mean opacity is expected to be larger than the electron scattering value. We show that the iron opacity bump causes the disk to be convective unstable. Turbulence generated by convection puffs up the disk due to additional turbulent pressure support and enhances the local angular momentum transport. This also results in strong fluctuations in surface density and heating of the disk. The opacity drops with increasing temperature and convection is suppressed. The disk cools down and the whole process repeats again. This causes strong oscillations of the disk scale height and luminosity variations by more than a factor of $\approx 3-6$ over a few years' timescale. Since the iron opacity bump will move to different locations of the disk for black holes with different masses and accretion rates, we suggest that this is a physical mechanism that can explain the variability of AGN with a wide range of amplitudes over a time scale of years to decades.

We analyze two 3D general-relativistic magnetohydrodynamic accretion simulations in the context of how they would manifest in Event Horizon Telescope (EHT) observations of supermassive black holes. The two simulations differ only in whether the initial angular momentum of the plasma is aligned with the rapid (a = 0.9) spin of the black hole. Both have low net magnetic flux. Ray tracing is employed to generate resolved images of the synchrotron emission. When using parameters appropriate for Sgr A* and assuming a viewing angle aligned with the black hole spin, we find the most prominent difference is that the central shadow in the image is noticeably eccentric in tilted models, with the ring of emission relatively unchanged. Applying this procedure to M87 with a viewing angle based on the large-scale jet, we find that adding tilt increases the angular size of the ring for fixed black hole mass and distance, while at the same time increasing the number of bright spots in the image. Our findings illustrate observable features that can distinguish tilted from aligned flows. They also show that tilted models can be viable for M87, and that not accounting for tilt can bias inferences of physical parameters. Future modeling of horizon-scale observations should account for potential angular momentum misalignment, which is likely generic at the low accretion rates appropriate for EHT targets.

Inverse Compton scattering dominates the high energy part of the spectra of neutron star (NS) low mass X-ray binaries (LMXBs). It has been proposed that inverse Compton scattering also drives the radiative properties of kilohertz quasi periodic oscillations (kHz QPOs). In this work, we construct a model that predicts the energy dependence of the rms amplitude and time lag of the kHz QPOs. Using this model, we fit the rms amplitude and time lag energy spectra of the lower kHz QPO in the NS LMXB 4U 1636-53 over 11 frequency intervals of the QPO and report three important findings: (i) A medium that extends 1-8 km above the NS surface is required to fit the data; this medium can be sustained by the balance between gravity and radiation pressure, without forcing any equilibrium condition. (ii) We predict a time delay between the oscillating NS temperature, due to feedback, and the oscillating electron temperature of the medium which, with the help of phase resolved spectroscopy, can be used as a probe of the geometry and the feedback mechanism. (iii) We show that the observed variability as a function of QPO frequency is mainly driven by the oscillating electron temperature of the medium. This provides strong evidence that the Comptonising medium in LMXBs significantly affects, if not completely drives, the radiative properties of the lower kHz QPOs regardless of the nature of the dynamical mechanism that produces the QPO frequencies.

We use global three dimensional radiation magneto-hydrodynamic simulations to study the properties of inner regions of accretion disks around a 5\times 10^8 solar mass black hole with mass accretion rates reaching 7% and 20% of the Eddington value. This region of the disk is supported by magnetic pressure with surface density significantly smaller than the values predicted by the standard thin disk model but with a much larger disk scale height. The disks do not show any sign of thermal instability over many thermal time scales. More than half of the accretion is driven by radiation viscosity in the optically thin corona region for the lower accretion rate case, while accretion in the optically thick part of the disk is driven by the Maxwell and Reynolds stresses from MRI turbulence. Coronae with gas temperatures > 10^8 K are generated only in the inner \approx 10 gravitational radii in both simulations, being more compact in the higher accretion rate case. In contrast to the thin disk model, surface density increases with increasing mass accretion rate, which causes less dissipation in the optically thin region and a relatively weaker corona. The simulation results may explain the formation of X-ray coronae in Active Galactic Nuclei (AGNs), the compact size of such coronae, and the observed trend of optical to X-ray luminosity with Eddington ratio for many AGNs.

We conduct a systematic study of the properties of tilted accretion flows around spinning black holes, covering a range of tilt angles and black hole spins, using the general-relativistic magnetohydrodynamics code Athena++. The same initial magnetized torus is evolved around black holes with spins ranging from 0 to 0.9, with inclinations ranging from 0 degrees to 24 degrees. The tilted disks quickly reach a warped and twisted shape that rigidly precesses about the black hole spin axis with deformations in shape large enough to hinder the application of linear bending wave theory. Magnetized polar outflows form, oriented along the disk rotation axes. At sufficiently high inclinations a pair of standing shocks develops in the disks. These shocks dramatically affect the flow at small radii, driving angular momentum transport. At high spins they redirect material more effectively than they heat it, reducing the dissipation rate relative to the mass accretion rate and lowering the radiative efficiency of the flow.

(Abridged) Stars more massive than $20-30M_{\odot}$ are so luminous that the radiation force on the cooler, more opaque outer layers can balance or exceed the force of gravity. These near or super-Eddington outer envelopes represent a long standing challenge for calculating the evolution of massive stars in one dimension, a situation that limits our understanding of the stellar progenitors of some of the most exciting and energetic explosions in the universe. In particular, the proximity to the Eddington limit has been the suspected cause for the variability, large mass loss rate and giant eruptions of an enigmatic class of massive stars: the luminous blue variables (LBVs). When in quiescence, LBVs are usually found on the hot ($T_{eff} \approx 2 - 4 \times 10^4$ K) S Dor instability strip. While in outburst, most LBVs stay on the cold S Dor instability strip with a $T_{eff} \approx 9000$ K. Here we show that physically realistic three dimensional global radiation hydrodynamic simulations of radiation dominated massive stars with the largest supercomputers in the world naturally reproduce many observed properties of LBVs, specifically their location in the Hertzsprung-Russell (HR) diagram and their episodic mass loss with rates of $10^{-7}-10^{-5} M_{\odot}/yr$. The helium opacity peak is found to play an important role to determine these properties, which is not realized in the traditional one dimensional models of massive stars. The simulations also predict that convection causes irregular envelope oscillations and 10-30% brightness variations on a typical timescale of a few days. The variability is more prominent in our models that are on the cool part of the S Dor instability. These calculations pave the way to a quantitative understanding of the structure, stability and the dominant mode of mass loss of massive stars.

Warm Comptonization models for the soft X-ray excess in AGN do not self-consistently explain the relationship between the Comptonizing medium and the underlying accretion disc. Because of this, they cannot directly connect the fitted Comptonization temperatures and optical depths to accretion disc parameters. Since bulk velocities exceed thermal velocities in highly radiation pressure dominated discs, in these systems bulk Comptonization by turbulence may provide a physical basis in the disc itself for warm Comptonization models. We model the dependence of bulk Comptonization on fundamental accretion disc parameters, such as mass, luminosity, radius, spin, inner boundary condition, and $\alpha$. In addition to constraining warm Comptonization models, our model can help distinguish contributions from bulk Comptonization to the soft X-ray excess from those due to other physical mechanisms, such as absorption and reflection. By linking the time variability of bulk Comptonization to fluctuations in the disc vertical structure due to magnetorotational instability (MRI) turbulence, our results show that observations of the soft X-ray excess can be used to study disc turbulence in the radiation pressure dominated regime. Because our model connects bulk Comptonization to one dimensional vertical structure temperature profiles in a physically intuitive way, it will be useful for understanding this effect in future simulations run in new regimes.

We present the results of local, vertically stratified, radiation magnetohydrodynamic shearing box simulations of magnetorotational instability (MRI) turbulence for a (hydrogen poor) composition applicable to accretion disks in AM CVn type systems. Many of these accreting white dwarf systems are helium analogues of dwarf novae (DNe). We utilize frequency-integrated opacity and equation of state tables appropriate for this regime to accurately portray the relevant thermodynamics. We find bistability of thermal equilibria in the effective temperature, surface mass density plane typically associated with disk instabilities. Along this equilibrium curve (i.e. the S-curve) we find that the stress to thermal pressure ratio $\alpha$ varied with peak values of $\sim 0.15$ near the tip of the upper branch. Similar to DNe, we found enhancement of $\alpha$ near the tip of the upper branch caused by convection; this increase in $\alpha$ occurred despite our choice of zero net vertical magnetic flux. Two notable differences we find between DN and AM CVn accretion disk simulations are that AM CVn disks are capable of exhibiting persistent convection in outburst, and ideal MHD is valid throughout quiescence for AM CVns. In contrast, DNe simulations only show intermittent convection, and non-ideal MHD effects are likely important in quiescence. By combining our previous work with these new results, we also find that convective enhancement of the MRI is anticorrelated with mean molecular weight.

Vertically stratified shearing box simulations of magnetorotational turbulence commonly exhibit a so-called butterfly diagram of quasi-periodic azimuthal field reversals. However, in the presence of hydrodynamic convection, field reversals no longer occur. Instead, the azimuthal field strength fluctuates quasi-periodically while maintaining the same polarity, which can either be symmetric or antisymmetric about the disc midplane. Using data from the simulations of Hirose et al. (2014), we demonstrate that the lack of field reversals in the presence of convection is due to hydrodynamic mixing of magnetic field from the more strongly magnetized upper layers into the midplane, which then annihilate field reversals that are starting there. Our convective simulations differ in several respects from those reported in previous work by others, in which stronger magnetization likely plays a more important role than convection.

Bulk velocities exceed thermal velocities for sufficiently radiation pressure dominated accretion flows. We model the contribution of bulk Comptonization to the soft X-ray excess in AGN. Bulk Comptonization is due to both turbulence and the background shear. We calculate spectra both taking into account and not taking into account bulk velocities using scaled data from radiation magnetohydrodynamic (MHD) shearing box simulations. We characterize our results with temperatures and optical depths to make contact with other warm Comptonization models of the soft excess. We chose our fiducial mass, $M = 2 \times 10^6 M_{\odot}$, and accretion rate, $L/L_{\rm Edd} = 2.5$, to correspond to those fit to the super-Eddington narrow line Seyfert 1 (NLS1) RE1034+396. The temperatures, optical depths, and Compton $y$ parameters we find broadly agree with those fit to RE1034+396. The effect of bulk Comptonization is to shift the Wien tail to higher energy and lower the gas temperature, broadening the spectrum. Observations of the soft excess in NLS1s can constrain the properties of disc turbulence if the bulk Comptonization contribution can be separated out from contributions from other physical effects, such as reflection and absorption.

We use three dimensional radiation magneto-hydrodynamic simulations to study the effects of magnetic fields on the energy transport and structure of radiation pressure dominated main sequence massive star envelopes at the region of the iron opacity peak. We focus on the regime where the local thermal timescale is shorter than the dynamical timescale, corresponding to inefficient convective energy transport. We begin with initially weak magnetic fields relative to the thermal pressure, from 100-1000G in differing geometries. The unstable density inversion amplifies the magnetic field, increasing the magnetic energy density to values close to equipartition with the turbulent kinetic energy density. By providing pressure support, the magnetic field's presence significantly increases the density fluctuations in the turbulent envelope, thereby enhancing the radiative energy transport by allowing photons to diffuse out through low density regions. Magnetic buoyancy brings small scale magnetic fields to the photosphere and increases the vertical energy transport with the energy advection velocity proportional to the Alfv\'en velocity, although in all cases we study photon diffusion still dominates the energy transport. The increased radiative and advective energy transport causes the stellar envelope to shrink by several scale heights. We also find larger turbulent velocity fluctuations compared to the purely hydrodynamic case, reaching $\approx$ 100 km/s at the stellar photosphere. The photosphere also shows vertical oscillations with similar averaged velocities and periods of a few hours. The increased turbulent velocity and oscillations will have strong impacts on the line broadening and periodic signals in massive stars.

The phenomenological Disc Instability Model has been successful in reproducing the observed light curves of dwarf nova outbursts by invoking an enhanced Shakura-Sunyaev $\alpha$ parameter $\sim0.1-0.2$ in outburst compared to a low value $\sim0.01$ in quiescence. Recent thermodynamically consistent simulations of magnetorotational (MRI) turbulence with appropriate opacities and equation of state for dwarf nova accretion discs have found that thermal convection enhances $\alpha$ in discs in outburst, but only near the hydrogen ionization transition. At higher temperatures, convection no longer exists and $\alpha$ returns to the low value comparable to that in quiescence. In order to check whether this enhancement near the hydrogen ionization transition is sufficient to reproduce observed light curves, we incorporate this MRI-based variation in $\alpha$ into the Disc Instability Model, as well as simulation-based models of turbulent dissipation and convective transport. These MRI-based models can successfully reproduce observed outburst and quiescence durations, as well as outburst amplitudes, albeit with different parameters from the standard Disc Instability Models. The MRI-based model lightcurves exhibit reflares in the decay from outburst, which are not generally observed in dwarf novae. However, we highlight the problematic aspects of the quiescence physics in the Disc Instability Model and MRI simulations that are responsible for this behavior.

Motivated by recent studies showing an apparent correlation between the high-frequency quasi-periodic oscillations (QPOs) and the low-frequency, type-C QPO in low-mass, black hole X-ray binaries (LMXBs), we explore a model that explains all three QPOs in terms of an oscillating, precessing hot flow in the truncated-disk geometry. Our model favors attributing the two high-frequency QPOs, often occurring in a near 3:2 frequency ratio, to the breathing and vertical epicyclic frequency modes of the hot, thick flow, although we can not rule out the Keplerian and m=-1 radial epicyclic modes. In either case, the type-C QPO is attributed to precession. The correlation of the QPOs comes from the fact that all three frequencies are associated with the same geometrical structure. While the exact QPO frequencies are sensitive to the black hole mass and spin, their evolution over the course of an outburst is mainly tied to the truncation radius between the geometrically thin, optically thick disk and the inner, hot flow. We show that, in the case of the LMXB GRO J1655-40, this model can explain the one simultaneous observation of all three QPOs and that an extrapolation of the model appears to match lower frequency observations where only two of the three components are seen. Thus, this model may be able to unify multiple QPO observations using the properties of a single, simple, geometrical model.

Radiation pressure dominated accretion discs around compact objects may have turbulent velocities that greatly exceed the electron thermal velocities within the disc. Bulk Comptonization by the turbulence may therefore dominate over thermal Comptonization in determining the emergent spectrum. Bulk Comptonization by divergenceless turbulence is due to radiation viscous dissipation only. It can be treated as thermal Comptonization by solving the Kompaneets equation with an equivalent "wave" temperature, which is a weighted sum over the power present at each scale in the turbulent cascade. Bulk Comptonization by turbulence with non-zero divergence is due to both pressure work and radiation viscous dissipation. Pressure work has negligible effect on photon spectra in the limit of optically thin turbulence, and in this limit radiation viscous dissipation alone can be treated as thermal Comptonization with a temperature equivalent to the full turbulent power. In the limit of extremely optically thick turbulence, radiation viscous dissipation is suppressed, and the evolution of local photon spectra can be understood in terms of compression and expansion of the strongly coupled photon and gas fluids. We discuss the consequences of these effects for self-consistently resolving and interpreting turbulent Comptonization in spectral calculations in radiation MHD simulations of high luminosity accretion flows.

We perform three dimensional radiation hydrodynamic simulations of the structure and dynamics of radiation dominated envelopes of massive stars at the location of the iron opacity peak. One dimensional hydrostatic calculations predict an unstable density inversion at this location, whereas our simulations reveal a complex interplay of convective and radiative transport whose behavior depends on the ratio of the photon diffusion time to the dynamical time. The latter is set by the ratio of the optical depth per pressure scale height, $\tau_0$, to $\tau_c=c/c_g$, where $c_g \approx$ 50 km/s is the isothermal sound speed in the gas alone. When $\tau_0 \gg \tau_c$, convection reduces the radiation acceleration and removes the density inversion. The turbulent energy transport in the simulations agrees with mixing length theory and provides its first numerical calibration in the radiation dominated regime. When $\tau_0 \ll \tau_c$, convection becomes inefficient and the turbulent energy transport is negligible. The turbulent velocities exceed $c_g$, driving shocks and large density fluctuations that allow photons to preferentially diffuse out through low-density regions. However, the effective radiation acceleration is still larger than the gravitational acceleration so that the time average density profile contains a modest density inversion. In addition, the simulated envelope undergoes large-scale oscillations with periods of a few hours. The turbulent velocity field may affect the broadening of spectral lines and therefore stellar rotation measurements in massive stars, while the time variable outer atmosphere could lead to variations in their mass loss and stellar radius.

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.

We propose a new model of the steep power law state of luminous black hole X-ray binaries. The model uses the fact that at high luminosities, the inner radii of radiation pressure dominated accretion discs are expected to i) become effectively optically thin and ii) produce significant luminosities. The gas temperature therefore rises sharply inwards, producing local saturated Compton spectra with rapidly increasing peak energy. These spectra sum together to form a steep power law tail to the spectrum. A given photon energy on this tail corresponds to a narrow range in radius, so that local vertical oscillations of the disc naturally produce high quality high-frequency quasi-periodic oscillations (HFQPOs) in the hard X-ray band. The two lowest order modes have a robust frequency ratio of ~1.53. This model explains the appearance of steep power law spectra and HFQPOs at high luminosity, the 3:2 HFQPO frequency ratios, and their association with the power law spectral component. We predict an increase in QPO quality factor when the power spectrum is restricted to a narrower photon energy band, and an increase in HFQPO frequency at higher X-ray energies or lower luminosities. Future X-ray telescopes could detect additional HFQPOs from higher order modes. We demonstrate how this model could be used to measure black hole spin from HFQPOs, and qualitatively estimate the spin of GRO J1655-40 as a/M ~ 0.4 - 0.7.

Simulations of hot, pressure supported, tilted black hole accretion flows, in which the angular momentum of the flow is misaligned with the black hole spin axis, can exhibit two non-axisymmetric shock structures in the inner regions of the flow. We analyze the strength and significance of these shock structures in simulations with tilt angles of 10 and 15 degrees. By integrating fluid trajectories in the simulations through the shocks, and tracking the variations of fluid quantities along these trajectories, we show that these shocks are strong, with substantial compression ratios, in contrast to earlier claims. However, they are only moderately relativistic. We also show that the two density enhancements resembling flow streams in their shape are in fact merely post-shock compressions, as fluid trajectories cut across, rather than flow along, them. The dissipation associated with the shocks is a substantial fraction ($\simeq3-12$ percent) of the rest mass energy advected into the hole, and therefore comparable to the dissipation expected from turbulence. The shocks should therefore make order unity changes in the observed properties of black hole accretion flows that are tilted.

We present spectral calculations of non-LTE accretion disk models appropriate for high luminosity stellar mass black hole X-ray binary systems. We first use a dissipation profile based on scaling the results of shearing box simulations of Hirose et al. (2009) to a range of annuli parameters. We simultaneously scale the effective temperature, orbital frequency and surface density with luminosity and radius according to the standard \alpha-model (Shakura & Sunyaev, 1973). This naturally brings increased dissipation to the disk surface layers (around the photospheres) at small radii and high luminosities. We find that the local spectrum transitions directly from a modified black body to a saturated Compton scattering spectrum as we increase the effective temperature and orbital frequency while decreasing midplane surface density. Next, we construct annuli models based on the parameters of a L/L_Edd=0.8 disk orbiting a 6.62 solar mass black hole using two modified dissipation profiles that explicitly put more dissipation per unit mass near the disk surface. The new dissipation profiles are qualitatively similar to the one found by Hirose et al. (2009), but produce strong near power-law spectral tails. Our models also include physically motivated magnetic acceleration support based once again on scaling the Hirose et al. (2009) results. We present three full-disk spectra each based on one of the dissipation prescriptions. Our most aggressive dissipation profile results in a disk spectrum that is in approximate quantitative agreement with certain observations of the steep power law (SPL) spectral states from some black hole X-ray binaries.

I provide a broad overview of the basic theoretical paradigms of black hole accretion flows. Models that make contact with observations continue to be mostly based on the four decade old alpha stress prescription of Shakura & Sunyaev (1973), and I discuss the properties of both radiatively efficient and inefficient models, including their local properties, their expected stability to secular perturbations, and how they might be tied together in global flow geometries. The alpha stress is a prescription for turbulence, for which the only existing plausible candidate is that which develops from the magnetorotational instability (MRI). I therefore also review what is currently known about the local properties of such turbulence, and the physical issues that have been elucidated and that remain uncertain that are relevant for the various alpha-based black hole accretion flow models.

We study the spatial and temporal behavior of fluid in fully three-dimensional, general relativistic, magnetohydrodynamical simulations of both tilted and untilted black hole accretion flows. We uncover characteristically greater variability in tilted simulations at frequencies similar to those predicted by the formalism of trapped modes, but ultimately conclude that its spatial structure is inconsistent with a modal interpretation. We find instead that previously identified, transient, over-dense clumps orbiting on roughly Keplerian trajectories appear generically in our global simulations, independent of tilt. Associated with these fluctuations are acoustic spiral waves interior to the orbits of the clumps. We show that the two nonaxisymmetric standing shock structures that exist in the inner regions of these tilted flows effectively amplify the variability caused by these spiral waves to markedly higher levels than in untilted flows, which lack standing shocks. Our identification of clumps, spirals, and spiral-shock interactions in these fully general relativistic, magnetohydrodynamical simulations suggests that these features may be important dynamical elements in models which incorporate tilt as a way to explain the observed variability in black hole accretion flows.

We study local radiation magnetohydrodynamic instabilities in static, optically thick, vertically stratified media with constant flux mean opacity. We include the effects of vertical gradients in a horizontal background magnetic field. Assuming rapid radiative diffusion, we use the zero gas pressure limit as an entry point for investigating the coupling between the photon bubble instability and the Parker instability. Apart from factors that depend on wavenumber orientation, the Parker instability exists for wavelengths longer than a characteristic wavelength lambda_{tran}, while photon bubbles exist for wavelengths shorter than lambda_{tran}. The growth rate in the Parker regime is independent of the orientation of the horizontal component of the wavenumber when radiative diffusion is rapid, but the range of Parker-like wavenumbers is extended if there exists strong horizontal shear between field lines (i.e. horizontal wavenumber perpendicular to the magnetic field). Finite gas pressure introduces an additional short wavelength limit to the Parker-like behavior, and also limits the growth rate of the photon bubble instability to a constant value at short wavelengths. We also consider the effects of differential rotation with accretion disk applications in mind. Our results may explain why photon bubbles have not yet been observed in recent stratified shearing box accretion disk simulations. Photon bubbles may physically exist in simulations with high radiation to gas pressure ratios, but higher spatial resolution will be needed to resolve the asymptotically growing unstable wavelengths.

(Abridged) Narrow Line Seyfert 1 (NLS1) galaxies have low mass black holes and mass accretion rates close to (or exceeding) Eddington, so a standard blackbody accretion disc should peak in the EUV. However, the lack of true absorption opacity in the disc means that the emission is better approximated by a colour temperature corrected blackbody, and this colour temperature correction is large enough ($\sim 2.4$) that the bare disc emission from a zero spin black hole can extend into the soft X-ray bandpass. Part of the soft X-ray excess seen in these objects must be intrinsic emission from the disc unless the vertical structure is very different to that predicted. However, the soft excess is much broader than predicted by a bare disc spectrum, indicating some Compton upscattering by cool, optically thick material. We associate this with the disc itself, so it must ultimately be powered by mass accretion. We build an energetically self consistent model assuming that the emission thermalises at large radii, but that at smaller radii the gravitational energy is split between powering optically thick Comptonised disc emission (forming the soft X-ray excess) and an optically thin corona above the disc (forming the tail to higher energies). We show examples of this model fit to the extreme NLS1 REJ1034+396, and to the much lower Eddington fraction Broad Line Seyfert 1 PG1048+231. We use these to guide our fits and interpretations of three template spectra made from co-adding multiple sources to track out a sequence of AGN spectra as a function of $L/L_{Edd}$. The new model is publically available within the {\sc xspec} spectral fitting package.

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.

We analyze the time dependence of fluid variables in general relativistic, magnetohydrodynamic simulations of accretion flows onto a black hole with dimensionless spin parameter a/M=0.9. We consider both the case where the angular momentum of the accretion material is aligned with the black hole spin axis (an untilted flow) and where it is misaligned by 15 degrees (a tilted flow). In comparison to the untilted simulation, the tilted simulation exhibits a clear excess of inertial variability, that is, variability at frequencies below the local radial epicyclic frequency. We further study the radial structure of this inertial-like power by focusing on a radially extended band at 118 (M/10Msol)^-1 Hz found in each of the three analyzed fluid variables. The three dimensional density structure at this frequency suggests that the power is a composite oscillation whose dominant components are an over dense clump corotating with the background flow, a low order inertial wave, and a low order inertial-acoustic wave. Our results provide preliminary confirmation of earlier suggestions that disk tilt can be an important excitation mechanism for inertial waves.

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.

A crucial difficulty in understanding the nature of the putative accretion disk in AGNs is that some of its key intrinsic spectral signatures cannot be observed directly. The strong emissions from the broad-line region (BLR) and the obscuring torus, which are generally yet to be spatially resolved, essentially 'bury' such signatures. Here we argue that we can actually isolate the disk emission spectrum by using optical and near-infrared polarization of quasars and uncover the important spectral signatures. In these quasars, the polarization is considered to originate from electron scattering interior to the BLR, so that the polarized flux shows the disk spectrum with all the emissions from the BLR and torus eliminated. The polarized flux observations have now revealed a Balmer edge feature in absorption and a blue near-infrared spectral shape consistent with a specific and robust theoretical prediction. These results critically verify the long-standing picture of an optically thick and locally heated disk in AGNs.

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.

Quasars are thought to be powered by supermassive black holes accreting surrounding gas. Central to this picture is a putative accretion disk which is believed to be the source of the majority of the radiative output. It is well known, however, that the most extensively studied disk model -- an optically thick disk which is heated locally by the dissipation of gravitational binding energy -- is apparently contradicted by observations in a few major respects. In particular, the model predicts a specific blue spectral shape asymptotically from the visible to the near-infrared, but this is not generally seen in the visible wavelength region where the disk spectrum is observable. A crucial difficulty was that, toward the infrared, the disk spectrum starts to be hidden under strong hot dust emission from much larger but hitherto unresolved scales, and thus has essentially been impossible to observe. Here we report observations of polarized light interior to the dust-emiting region that enable us to uncover this near-infrared disk spectrum in several quasars. The revealed spectra show that the near-infrared disk spectrum is indeed as blue as predicted. This indicates that, at least for the outer near-infrared-emitting radii, the standard picture of the locally heated disk is approximately correct. The model problems at shorter wavelengths should then be directed toward a better understanding of the inner parts of the revealed disk. The newly uncovered disk emission at large radii, with more future measurements, will also shed totally new light on the unanswered critical question of how and where the disk ends.

This work presents a detailed analysis of the overall flow structure and unique features of the inner region of the tilted disk simulations described in Fragile et al. (2007). The primary new feature identified in the main disk body is a latitude-dependent radial epicyclic motion driven by pressure gradients attributable to the gravitomagnetic warping of the disk. The induced motion of the gas is coherent over the scale of the entire disk and is fast enough that it could be observable in features such as relativistic iron lines. The eccentricity of the associated fluid element trajectories increases with decreasing radius, leading to a crowding of orbit trajectories near their apocenters. This results in a local density enhancement akin to a compression. These compressions are sufficiently strong to produce a pair of weak shocks in the vicinity of the black hole. These shocks are roughly aligned with the line-of-nodes between the black-hole symmetry plane and disk midplane, with one shock above the line-of-nodes on one side of the black hole and the other below on the opposite side. These shocks enhance angular momentum transport and energy dissipation near the hole, forcing some material to plunge toward the black hole from well outside the innermost stable circular orbit. A new, extended simulation, which was evolved for more than a full disk precession period, allows us to confirm that these shocks and the previously identified ``plunging streams'' precess with the disk in such a way as to remain aligned relative to the line-of-nodes, as expected based on our physical understanding of these phenomena. Such a precessing structure would likely present a strong quasi-periodic signal.

We measure long (2200-4000 ang) and short (1450-2200 ang) wavelength spectral slopes \alpha (F_\nu proportional to \nu^\alpha) for quasar spectra from the Sloan Digital Sky Survey. The long and short wavelength slopes are computed from 3646 and 2706 quasars with redshifts in the z=0.76-1.26 and z=1.67-2.07 ranges, respectively. We calculate mean slopes after binning the data by monochromatic luminosity at 2200 ang and virial mass estimates based on measurements of the MgII line width and 3000 ang continuum luminosity. We find little evidence for mass dependent variations in the mean slopes, but a significant luminosity dependent trend in the near UV spectral slopes is observed with larger (bluer) slopes at higher luminosities. The far UV slopes show no clear variation with luminosity and are generally lower (redder) than the near UV slopes at comparable luminosities, suggesting a slightly concave quasar continuum shape. We compare these results with Monte Carlo distributions of slopes computed from models of thin accretion disks, accounting for uncertainties in the mass estimates. The model slopes produce mass dependent trends which are larger than observed, though this conclusion is sensitive to the assumed uncertainties in the mass estimates. The model slopes are also generally bluer than observed, and we argue that reddening by dust intrinsic to the source or host galaxy may account for much of the discrepancy.

We study epicyclic oscillations of fluid tori around black holes (in the Paczynski-Wiita potential) and derive exact analytic expressions for their radial and vertical eigenfrequencies nu_r and nu_z to second-order accuracy in the width of the torus. We prove that pressure effects make the eigenfrequencies smaller than those for free particles. However, the particular ratio nu_z/nu_r=3/2, which is important for the theory of high-frequency quasi-periodic oscillations (QPOs), occurs when the fluid tori epicyclic frequencies nu_r and nu_z are about 15% higher than the ones corresponding to free particles. Our results therefore suggest that previous estimates of black hole spins from QPOs have produced values that are too high.

This paper presents a continuation of our efforts to numerically study accretion disks that are misaligned (tilted) with respect to the rotation axis of a Kerr black hole. Here we present results of a global numerical simulation which fully incorporates the effects of the black hole spacetime as well as magnetorotational turbulence that is the primary source of angular momentum transport in the flow. This simulation shows dramatic differences from comparable simulations of untilted disks. Accretion onto the hole occurs predominantly through two opposing plunging streams that start from high latitudes with respect to both the black-hole and disk midplanes. This is due to the aspherical nature of the gravitational spacetime around the rotating black hole. These plunging streams start from a larger radius than would be expected for an untilted disk. In this regard the tilted black hole effectively acts like an untilted black hole of lesser spin. Throughout the duration of the simulation, the main body of the disk remains tilted with respect to the symmetry plane of the black hole; thus there is no indication of a Bardeen-Petterson effect in the disk at large. The torque of the black hole instead principally causes a global precession of the main disk body. In this simulation the precession has a frequency of $3 (M_\odot/M)$ Hz, a value consistent with many observed low-frequency quasi-periodic oscillations. However, this value is strongly dependent on the size of the disk, so this frequency may be expected to vary over a large range.

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.

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.

I briefly review the theoretical models of radiatively efficient, geometrically thin and optically thick accretion disk spectra that currently exist for AGN. I then discuss three recent observational developments that have real potential to teach us about the physics of these flows. Finally, I present results on the most recent, thermodynamically consistent simulations of magnetorotational turbulence and discuss what these simulations are suggesting about the vertical structure of accretion disks.

We posit that accreting compact objects, including stellar mass black holes and neutron stars as well as supermassive black holes, may undergo extended periods of accretion during which the angular momentum of the disk at large scales is misaligned with that of the compact object. In such a scenario, Lense-Thirring precession caused by the rotating compact object can dramatically affect the disk. In this presentation we describe results from a three-dimensional relativistic magnetohydrodynamic simulation of an MRI turbulent disk accreting onto a tilted rapidly rotating black hole. For this case, the disk does not achieve the commonly described Bardeen-Petterson configuration; rather, it remains nearly planar, undergoing a slow global precession. Accretion from the disk onto the hole occurs predominantly through two opposing plunging streams that start from high latitudes with respect to both the black-hole and disk midplanes. This is a consequence of the non-sphericity of the gravitational spacetime of the black hole.

One primary difficulty in understanding the nature of the putative accretion disk in the central engine of AGNs is that some of its key intrinsic spectral signatures are buried under the emissions from the surrounding regions, i.e. the broad line region (BLR) and the obscuring torus. We argue here that these signatures can be revealed by using optical and near-IR polarization. At least in some quasars, the polarization is seen only in the continuum and is not shared by emission lines. In this case, the polarized flux is considered to show the intrinsic spectrum interior to the BLR, removing off the emissions from the BLR and torus. We have used this polarization to reveal the Balmer-edge feature and near-IR spectral shape of the central engine, both of which are important for testing the fundamental aspects of the models.

Turbulent Comptonization, a potentially important damping and radiation mechanism in relativistic accretion flows, is discussed. Particular emphasis is placed on the physical basis, relative importance, and thermodynamics of turbulent Comptonization. The effects of metal-absorption opacity on the spectral component resulting from turbulent Comptonization is considered as well.

We present results from spectral modeling of three black hole X-ray binaries: LMC X-3, GRO J1655-40, and XTE J1550-564. Using a sample of disk dominated observations, we fit the data with a range of spectral models that includes a simple multitemperature blackbody (DISKBB), a relativistic accretion disk model based on color-corrected blackbodies (KERRBB), and a relativistic model based on non-LTE atmosphere models within an alpha prescription (BHSPEC). BHSPEC provides the best fit for a BeppoSAX observation of LMC X-3, which has the broadest energy coverage of our sample. It also provides the best fit for multiple epochs of Rossi X-ray Timing Explorer (RXTE) data in this source, except at the very highest luminosity (L/L_{Edd} >~ 0.7), where additional physics must be coming into play. BHSPEC is also the best-fit model for multi-epoch RXTE observations of GRO J1655-40 and XTE J1550-564, although the best-fit inclination of the inner disk differs from the binary inclination. All our fits prefer alpha=0.01 to alpha=0.1, in apparent disagreement with the large stresses inferred from the rapid rise times observed in outbursts of these two sources. In all three sources our fits imply moderate black hole spins (a/M ~ 0.1 - 0.8), but this is sensitive to the reliability of independent measurements of these system parameters and to the physical assumptions which underly our spectral models.

Quasi-periodic oscillations (QPOs) in the X-ray lightcurves of accreting neutron star and black hole binaries have been widely interpreted as being due to standing wave modes in accretion disks. These disks are thought to be highly turbulent due to the magnetorotational instability (MRI). We study wave excitation by MRI turbulence in the shearing box geometry. We demonstrate that axisymmetric sound waves and radial epicyclic motions driven by MRI turbulence give rise to narrow, distinct peaks in the temporal power spectrum. Inertial waves, on the other hand, do not give rise to distinct peaks which rise significantly above the continuum noise spectrum set by MRI turbulence, even when the fluid motions are projected onto the eigenfunctions of the modes. This is a serious problem for QPO models based on inertial waves.

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.

Accretion flows with pressure gradients permit the existence of standing waves which may be responsible for observed quasi-periodic oscillations (QPO's) in X-ray binaries. We present a comprehensive treatment of the linear modes of a hydrodynamic, non-self-gravitating, polytropic slender torus, with arbitrary specific angular momentum distribution, orbiting in an arbitrary axisymmetric spacetime with reflection symmetry. We discuss the physical nature of the modes, present general analytic expressions and illustrations for those which are low order, and show that they can be excited in numerical simulations of relativistic tori. The mode oscillation spectrum simplifies dramatically for near Keplerian angular momentum distributions, which appear to be generic in global simulations of the magnetorotational instability. We discuss our results in light of observations of high frequency QPO's, and point out the existence of a new pair of modes which can be in an approximate 3:2 ratio for arbitrary black hole spins and angular momentum distributions, provided the torus is radiation pressure dominated. This mode pair consists of the axisymmetric vertical epicyclic mode and the lowest order axisymmetric breathing mode.

Fluids in external gravity may oscillate with frequencies characteristic of the epicyclic motions of test particles. We explicitly demonstrate that global oscillations of a slender, perfect fluid torus around a Kerr black hole admit incompressible vertical and radial epicyclic modes. Our results may be directly relevant to one of the most puzzling astrophysical phenomena -- high (hundreds of hertz) frequency quasiperiodic oscillations (QPOs) detected in X-ray fluxes from several black hole sources. Such QPOs are pairs of stable frequencies in the 3/2 ratio. It seems that they originate a few gravitational radii away from the black hole and thus observations of them have the potential to become an accurate probe of super-strong gravity.

The near-infrared shape of the big blue bump component in quasar spectra has been essentially unknown. It usually cannot be observed directly, due to the strong hot dust emission which dominates quasar spectra longward of ~1micron. However this is quite an important part of the spectrum theoretically. At least bare disk models provide quite a robust prediction for the overall continuum shape in the near-infrared. Self-gravity should become important in the outer, near-infrared emitting regions of the putative disk, possibly leaving a signature of disk truncation in the near-infrared. We propose here that this important part of the spectrum can be revealed for the first time by observing polarized flux from normal quasars. At least in some polarized quasars, the emission lines are all unpolarized and so the polarized flux should originate interior to the broad line region, and therefore also interior to the dust emitting region. This can then be used to eliminate the dust emission. We present the results of near-infrared polarimetry for such three quasars (Ton202, 4C37.43, B2 1208+32). The data for Ton202 have the highest S/N, and the near-infrared polarized flux in this case is measured to have quite a blue shape, nu^+0.42+-0.29 in F_nu, intriguingly consistent with the simple multi-temperature black body, bare disk prediction of nu^+1/3. All these data, although still with quite low S/N for the other two objects, demonstrate the unique potential of the technique with future better data. We also present similar data for other quasars and radio galaxies, and briefly discuss the nature of the polarization.

We examine the effects of photon bubble instability in radiation-dominated accretion disks such as those found around black holes in active galactic nuclei and X-ray binary star systems. Two- and 3-D numerical radiation MHD calculations of small patches of disk are used. Modes with wavelengths shorter than the gas pressure scale height grow faster than the orbital frequency in the surface layers. The fastest growth rate observed is five times the orbital frequency and occurs on nearly-vertical magnetic fields. The spectrum of linear modes agrees with a WKB analysis indicating still faster growth at unresolved scales, with a maximum proportional to the gravity and inversely proportional to the gas sound speed. Disturbances reaching non-linear amplitudes steepen into trains of shocks similar to a 1-D periodic non-linear analytic solution. Variations in propagation speed result in merging of adjacent fronts, and over time the shock spacing and amplitude increase. Growth is limited by the strength of the field, and the structure is disrupted when the ram pressure exceeds the magnetic pressure. The largest horizontal density variations are similar to the ratio of magnetic to gas pressure, and in our calculations are more than 100. Under the conditions considered, radiation diffuses through the inhomogeneneous flow 5x faster than through the initial hydrostatic equilibrium, and the net cooling rate is several times greater than in a similar calculation without magnetic fields showing the effects of convection. These results indicate photon bubbles may be important in cooling radiation-dominated disks.

An explanation for the large random velocities of pulsars is presented. Like many other models, we propose that the momentum imparted to the star is given at birth. The ultimate source of energy is provided by the intense optically thick neutrino flux that is responsible for radiating the proto-neutron star's gravitational binding energy during the Kelvin-Helmholtz phase. The central feature of the kick mechanism is a radiative-driven magnetoacoustic instability, which we refer to as ``neutrino bubbles.'' Identical in nature to the photon bubble instability, the neutrino bubble instability requires the presence of an equilibrium radiative flux as well as a coherent steady background magnetic field. Over regions of large magnetic flux densities, the neutrino bubble instability is allowed to grow on dynamical timescales ~ 1ms, potentially leading to large luminosity enhancements and density fluctuations. Local luminosity enhancements, which preferentially occur over regions of strong magnetic field, lead to a net global asymmetry in the neutrino emission and the young neutron star is propelled in the direction opposite to these regions. For favorable values of magnetic field structure, size, and strength as well as neutrino bubble saturation amplitude, momentum kicks in excess of 1000 km/s can be achieved. Since the neutrino-powered kick is delivered over the duration of the Kelvin-Helmholtz time ~ a few seconds, one expects spin-kick alignment from this neutrino bubble powered model.

We investigate the ultraviolet-to-optical spectral energy distributions (SEDs) of 17 active galactic nuclei (AGNs) using quasi-simultaneous spectrophotometry spanning 900-9000 Angstrom (rest frame). We employ data from the Far Ultraviolet Spectroscopic Explorer (FUSE), the Hubble Space Telescope (HST), and the 2.1-meter telescope at Kitt Peak National Observatory (KPNO). Taking advantage of the short-wavelength coverage, we are able to study the so-called "big blue bump," the region where the energy output peaks, in detail. Most objects exhibit a spectral break around 1100 Angstrom. Although this result is formally associated with large uncertainty for some objects, there is strong evidence in the data that the far-ultraviolet spectral region is below the extrapolation of the near-ultraviolet-optical slope, indicating a spectral break around 1100 Angstrom. We compare the behavior of our sample to those of non-LTE thin-disk models covering a range in black-hole mass, Eddington ratio, disk inclination, and other parameters. The distribution of ultraviolet-optical spectral indices redward of the break, and far-ultraviolet indices shortward of the break, are in rough agreement with the models. However, we do not see a correlation between the far-ultraviolet spectral index and the black hole mass, as seen in some accretion disk models. We argue that the observed spectral break is intrinsic to AGNs, although intrinsic reddening as well as Comptonization can strongly affect the far-ultraviolet spectral index. We make our data available online in digital format.

We present calculations of non-LTE, relativistic accretion disk models applicable to the high/soft state of black hole X-ray binaries. We include the effects of thermal Comptonization and bound-free and free-free opacities of all abundant ion species. We present spectra calculated for a variety of accretion rates, black hole spin parameters, disk inclinations, and stress prescriptions. We also consider nonzero inner torques on the disk, and explore different vertical dissipation profiles, including some which are motivated by recent radiation MHD simulations of magnetorotational turbulence. Bound-free metal opacity generally produces significantly less spectral hardening than previous models which only considered Compton scattering and free-free opacity. It also tends to keep the effective photosphere near the surface, resulting in spectra which are remarkably independent of the stress prescription and vertical dissipation profile, provided little dissipation occurs above the effective photosphere. We provide detailed comparisons between our models and the widely used multicolor disk model. Frequency dependent discrepancies exist that may affect the parameters of other spectral components when this simpler disk model is used to fit modern X-ray data. For a given source, our models predict that the luminosity in the high/soft state should approximately scale with the fourth power of the empirically inferred maximum temperature, but with a slight hardening at high luminosities. This is in good agreement with observations. (abridged)

Polarimetric study in the UV and optical has been one of the keys to reveal the structure and nature of active galactic nuclei (AGN). Combined with the HST's high spatial resolution, it has directly confirmed the predicted scattering geometry of ~100 pc scale and accurately located the hidden nuclear position in some nearby active galaxies. Recently, we are using optical spectropolarimetry to reveal the nature of the accretion flow in the central engine of quasars. This is to use polarized light to de-contaminate the spectrum and investigate the Balmer edge spectral feature, which is otherwise buried under the strong emission from the outer region surrounding the central engine.

In our previous paper, we have reported the detection of a Balmer edge absorption feature in the polarized flux of one quasar (Ton 202). We have now found similar Balmer edge features in the polarized flux of four more quasars (4C09.72, 3C95, B2 1208+32, 3C323.1), and possibly a few more, out of 14 newly observed with the VLT and Keck telescopes. In addition, we also re-observed Ton 202, but we did not detect such a dramatic feature, apparently due to polarization variability (the two observations are one-year apart). The polarization measurements of some quasars are affected by an interstellar polarization in our Galaxy, but the measurements have been corrected for this effect reasonably well. Since the broad emission lines are essentially unpolarized and the polarization is confined only to the continuum in the five quasars including Ton 202 in both epochs, the polarized flux is considered to originate interior to the broad emission line region. The Balmer edge feature seen in the polarized flux is most simply interpreted as an intrinsic spectral feature of the quasar UV/optical continuum, or the ``Big Blue Bump'' emission. In this case, the edge feature seen in absorption indeed indicates the thermal and optically-thick nature of the continuum emitted. However, we also discuss other possible interpretations.

Spectroscopic studies indicate that gas in the photospheres of young O stars moves at speeds up to the sound speed. We show, using two-dimensional radiation MHD calculations and results from a local linear analysis, that the motions may be due to photon bubble instability if young O stars have magnetic fields.

The ultraviolet/optical continuum of quasars is thought to be from an accretion flow around a supermassive black hole, and it dominates the radiative output of quasars. However, the nature of this emission, often called the Big Blue Bump, has not been well understood. Robust evidence for its thermal nature would be a continuum opacity edge feature intrinsic to this component, but this has not been clearly confirmed despite the predictions by many models. We are now developing and exploiting a new method to detect the Balmer edge of hydrogen opacity. The method overcomes for certain objects the major obstacle of the heavy contamination from emissions outside the nucleus by taking the polarized flux spectrum. If our interpretation of the polarized flux is correct, our data show that the Big Blue Bump emission has a Balmer edge in absorption, indicating that the emission is indeed thermal, and the emitter is optically thick.

Numerical calculations and linear theory of radiation magnetohydrodynamic flows indicate that the photon bubble and magnetorotational instability (MRI) may produce large density inhomogeneities in radiation pressure supported media. We study the effects of the photon bubble instability on accretion disk spectra using 2-D Monte Carlo (MC) and 1-D Feautrier radiative transfer calculations on a snapshot of a 2-D numerical simulation domain. We find an enhancement in the thermalization of the MC spectra over that of the Feautrier calculation. In the inner-most regions of these disks, the turbulent magnetic pressure may greatly exceed that of the gas. It is then possible for bulk turbulent Alfvenic motions driven by the MRI to exceed the thermal velocity making turbulent Comptonization the dominant radiative process. We estimate the spectral distortion due to turbulent Comptonization utilizing a 1-D MC calculation.

Numerical simulations of the interiors of radiation dominated accretion disks show that significant density inhomogeneities can be generated in the gas. Here, we present the first results of our study on X-ray reflection spectra from such heterogeneous density structures. We consider two cases: first, we produce a number of toy models where a sharp increase or decrease in density of variable width is placed at different depths in a uniform slab. Comparing the resulting reflection spectra to those from an unaltered slab shows that the inhomogeneity can affect the emission features, in particular the Fe K and O VIII Lyalpha lines. The magnitude of any differences depends on both the parameters of the density change and the ionizing power of the illuminating radiation, but the inhomogeneity is required to be within ~2 Thomson depths of the surface to cause an effect. However, only relatively small variations in density (on the order of a few) are necessary for significant changes in the reflection features to be possible. Our second test was to compute reflection spectra from the density structure predicted by a simulation of the non-linear outcome of the photon bubble instability. The resulting spectra also exhibited differences from the constant density models, caused primarily by a strong 6.7 keV iron line. Nevertheless, constant density models can provide a good fit to simulated spectra, albeit with a low reflection fraction, between 2 and 10 keV. Below 2 keV, differences in the predicted soft X-ray line emission result in very poor fits with a constant density ionized disk model. The results indicate that density inhomogeneities may further complicate the relationship between the Fe K equivalent width and the X-ray continuum.

We model the spectral energy distribution of the type 1 Seyfert galaxies, fitting data from simultaneous optical, UV, and X-ray monitoring observations. We assume a geometry consisting of a hot central Comptonizing region surrounded by a thin accretion disk. The properties of the disk and the hot central region are determined by the feedback occurring between the hot Comptonizing region and thermal reprocessing in the disk that, along with viscous dissipation, provides the seed photons for the Comptonization process. The constraints imposed upon this model by existing multiwavelength data allow us to derive limits on the central black hole mass, the accretion rate, and the radius of the transition region between the thin outer disk and the geometrically thick, hot inner region in these objects.

Thermal models for the quasar Big Blue Bump generally lead to bound-free continuum features, which may be in absorption or emission. Searches for Lyman edges attributable to quasar atmospheres (in particular accretion disk atmospheres) have been ambiguous at best, but various relativistic and non-LTE effects may make them hard to detect. The Balmer edge features tentatively predicted by such models might be easier to detect since they'd form farther out in the potential wells. These can be sought in certain cases using spectropolarimetry to remove the effects of atomic emission (i.e. the Small Blue Bump) from the spectrum. We do find apparent Balmer edges in absorption in several quasars using this method! Although the features we see are believable and apparently common, more data and more modeling are needed to verify that they are best interpreted as atmospheric Balmer edges.

In the inner-most regions of radiation pressure supported accretion disks, the turbulent magnetic pressure may greatly exceed that of the gas. If this is the case, it is possible for bulk Alfvenic motions driven by the magnetorotational instability (MRI) to surpass the electron thermal velocity. Bulk rather than thermal Comptonization may then be the dominant radiative process which mediates gravitational energy release. For sufficiently large turbulent stresses, we show that turbulent Comptonization produces a significant contribution to the far-UV and X-ray emission of black hole accretion disks. The existence of this spectral component provides a means of obtaining direct observational constraints on the nature of the turbulence itself. We describe how this component may affect the spectral energy distributions and variability properties of X-ray binaries and active galactic nuclei.

We have found for the first time a Balmer edge feature in the polarized flux spectrum of a quasar (Ton 202). The edge feature is seen as a discontinuity in the slope, rather than as a discontinuity in the absolute flux. Since the polarized flux contains essentially no broad emission lines, it is considered to arise interior to the broad emission line region, showing the spectrum with all the emissions outside the nucleus scraped off and removed. Therefore, the polarized flux spectrum is likely to reveal features intrinsic to the Big Blue Bump emission. In this case, the existence of the Balmer edge feature, seen in absorption in the shorter wavelength side, indicates that the Big Blue Bump is indeed thermal and optically thick.

We examine the local conditions for radiative damping and driving of short wavelength, propagating hydrodynamic and magnetohydrodynamic (MHD) waves in static, optically thick, stratified equilibria. We show that so-called strange modes in stellar oscillation theory and magnetic photon bubbles are intimately related and are both fundamentally driven by the background radiation flux acting on compressible waves. We identify the necessary criteria for unstable driving of these waves, and show that this driving can exist in both gas and radiation pressure dominated media, as well as pure Thomson scattering media in the MHD case. The equilibrium flux acting on opacity fluctuations can drive both hydrodynamic acoustic waves and magnetosonic waves unstable. In addition, magnetosonic waves can be driven unstable by a combination of the equilibrium flux acting on density fluctuations and changes in the background radiation pressure along fluid displacements. We briefly describe the conditions under which these instabilities might be manifested in both main sequence stellar envelopes and accretion disks.

We have found for the first time a Balmer edge feature in the Big Blue Bump emission of a quasar. The feature is seen in the polarized flux spectrum of the quasar, where all the emissions from outside the nucleus are scraped off and removed. The existence of the Balmer-edge absorption feature directly indicates that the Big Blue Bump is indeed thermal and optically-thick.

We model the spectral energy distribution of the type 1 Seyfert galaxy NGC 5548, fitting data from simultaneous optical, UV, and X-ray monitoring observations. We assume a geometry consisting of a hot central Comptonizing region surrounded by a thin accretion disk. The properties of the disk and the hot central region are determined by the feedback occurring between the hot Comptonizing region and thermal reprocessing in the disk that, along with viscous dissipation, provides the seed photons for the Comptonization process. The constraints imposed upon this model by the multiwavelength data allow us to derive limits on the central black hole mass, $M \la 2 \times 10^7$ \Msun, the accretion rate, $\Mdot \la 2.5\times 10^5$ M$_\odot^2$yr$^{-1}/M$, and the radius of the transition region between the thin outer disk and the geometrically thick, hot inner region, $\sim 2$--5$\times 10^{14}$ cm.

These lectures provide an overview of the theory of accretion disks with application to bright sources containing black holes. I focus on the fundamental physics of these flows, stressing modern developments and outstanding questions wherever possible. After a review of standard Shakura-Sunyaev based models and their problems and uncertainties, I describe the basic principles that determine the overall spectral energy distribution produced by the flow. I then describe the physics of angular momentum transport in black hole accretion disks, stressing the important role of magnetic fields. Finally, I discuss the physics of radiation magnetohydrodynamics and how it might affect the overall flow structure in the innermost regions near the black hole.

We consider the dynamical evolution of bound, hierarchical triples of supermassive black holes that might be formed in the nuclei of galaxies undergoing sequential mergers. The tidal force of the outer black hole on the inner binary produces eccentricity oscillations through the Kozai mechanism, and this can substantially reduce the gravitational wave merger time of the inner binary. We numerically calculate the merger time for a wide range of initial conditions and black hole mass ratios, including the effects of octupole interactions in the triple as well as general relativistic periastron precession in the inner binary. The semimajor axes and the mutual inclination of the inner and outer binaries are the most important factors affecting the merger time. We find that for a random distribution of inclination angles and approximately equal mass black holes, it is possible to reduce the merger time of a near circular inner binary by more than a factor of ten in over fifty percent of all cases. We estimate that a typical exterior quadrupole moment from surrounding matter in the galaxy may also be sufficient to excite eccentricity oscillations in supermassive black hole binaries, and also accelerate black hole mergers.

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.

We present an application of the thermal Comptonization/disk reprocessing model recently proposed by Zdziarski, Lubi\'nski, & Smith. We show that the absence of strong optical variations in the presence of strong concurrent X-ray variations, similar to those found by HST/RXTE monitoring observations of NGC 3516, can be explained by changing the geometry of the Comptonizing plasma rather than the accretion disk itself. The total X-ray luminosity of the Comptonizing plasma must decrease as its spatial extent increases. In contrast, the disk inner radius must be roughly fixed in order not to produce optical/ultraviolet color variations stronger than observed. By including emission due to internal viscous dissipation in the disk, we can roughly match the optical and X-ray flux levels and variability amplitudes seen from NGC 3516 during the HST/RXTE campaign.

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.

We present a general linear dispersion relation which describes the coupled behavior of magnetorotational, photon bubble, and convective instabilities in weakly magnetized, differentially rotating accretion disks. We presume the accretion disks to be geometrically thin and supported vertically by radiation pressure. We fully incorporate the effects of a nonzero radiative diffusion length on the linear modes. In an equilibrium with purely vertical magnetic field, the vertical magnetorotational modes are completely unaffected by compressibility, stratification, and radiative diffusion. However, in the presence of azimuthal fields, which are expected in differentially rotating flows, the growth rate of all magnetorotational modes can be reduced substantially below the orbital frequency. This occurs if diffusion destroys radiation sound waves on the length scale of the instability, and the magnetic energy density of the azimuthal component exceeds the non-radiative thermal energy density. While sluggish in this case, the magnetorotational instability still persists and will still tap the free energy of the differential rotation. Photon bubble instabilities are generically present in radiation pressure dominated flows where diffusion is present. We show that their growth rates are limited to a maximum value which is reached at short wavelengths where the modes may be viewed as unstable slow magnetosonic waves. We also find that vertical radiation pressure destabilizes upward propagating fast waves, and that Alfv\'en waves can be unstable. These instabilities typically have smaller growth rates than the photon bubble/slow modes. We discuss how all these modes behave in various regimes of interest, and speculate how they may affect the dynamics of real accretion disk flows.

Recent spectropolarimetry of QSOs has revealed a surprising rise in polarization in the Lyman continuum of several objects. We discuss several recent attempts to interpret this feature, including the role of Lyman limit absorption in PG 1222+228. We present new theoretical results involving electron scattering in a hot corona or wind above an accretion disk, and polarization resulting from the relativistic returning radiation. These mechanisms can lead to potentially observable polarization at short wavelengths, but neither has quantitative success in fitting the observed Lyman continuum polarization rises. A renewed capability for ultraviolet spectropolarimetry from space is urgently needed to provide further clues to the nature of this phenomenon.

Using the Long Wavelength Spectrometer on Keck I, we have imaged the gravitationally lensed radio quiet quasi-stellar object (QSO) 2237+0305 at 8.9 and 11.7 microns for the first time. The mid-infrared flux ratios are inconsistent with the optical flux ratios, but agree with the radio flux ratios and with some published gravitational lens models. These flux ratios indicate that the infrared emission is not affected by microlensing, which rules out the synchrotron emission model. The infrared emission is likely produced by hot dust extended on a length scale of more than 0.03 parsecs. The spectral energy distribution further implies a narrow range of dust temperatures, suggesting that the dust may be located in a shell extending between ~1 to 3 parsecs from the nucleus, and intercepting about half of the QSO luminosity.

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.

We present simultaneous observations by EUVE, ASCA, and RXTE of the type~1 Seyfert galaxy NGC 5548. These data indicate that variations in the EUV emission (at $\sim 0.2 $keV) appear to lead similar modulations in higher energy ($\ga 1 $keV) X-rays by $\sim$10--30 ks. This is contrary to popular models which attribute the correlated variability of the EUV, UV and optical emission in type~1 Seyferts to reprocessing of higher energy radiation. This behavior instead suggests that the variability of the optical through EUV emission is an important driver for the variability of the harder X-rays which are likely produced by thermal Comptonization. We also investigate the spectral characteristics of the fluorescent iron K$\alpha$ line and Compton reflection emission. In contrast to prior measurements of these spectral features, we find that the iron K$\alpha$ line has a relatively small equivalent width ($W_{K\alpha} \sim 100 $eV) and that the reflection component is consistent with a covering factor which is significantly less than unity ($\Omega/2\pi \sim 0.4$--0.5). Notably, although the 2--10 keV X-ray flux varies by $\sim \pm 25$% and the derived reflection fraction appears to be constant throughout our observations, the flux in the Fe~K$\alpha$ line is also constant. This behavior is difficult to reconcile in the context of standard Compton reflection models.

Accretion disk models still do not provide a satisfactory explanation of the optical/ultraviolet continuum observed in Seyferts and quasars. Substantial theoretical progress has been made in understanding one aspect of the problem: the dearth of spectral features at the Lyman limit. Promising solutions have also been proposed to explain the surprising observations of large polarization in the Lyman continuum observed in some sources. I review the recent progress in this field, and try to point out future research directions which would be fruitful in trying to obtain a complete, self-consistent model of the continuum emitting regions.

We have calculated the spectrum and polarization of a standard thin accretion disk with parameters appropriate for a bright quasar. This model improves upon previous work by including ultraviolet metal line opacities, assumed for now to be in LTE. Though not yet fully self-consistent, our calculations demonstrate that metal lines can change the spectral slope, reduce the polarization, and reduce the Lyman edge feature in accretion disk spectra. Some observational differences between quasar spectra and accretion disk models might be reconciled with the inclusion of metal lines.

We present results of radiative transfer calculations of Compton scattering of polarized radiation in a corona above an accretion disk. The corona has a plane-parallel geometry and is assumed to be homogeneous with a fixed temperature. We employ a Feautrier method to calculate the radiative transfer, incorporating the full relativistic, polarized Compton redistribution matrix. This allows us to explore a broad range of parameter space, in particular large optical depths and/or low corona temperatures, for which the iterative scattering method is very inefficient. Compared to Thomson scattering, Compton scattering generally produces smaller polarizations for featureless spectra. Distinct frequency-dependent patterns in polarization are produced near lines and edges, which can be used as diagnostics for the corona. We discuss the physics which gives rise to these results. Compton scattering can smear out Lyman edges in total flux, and at the same time produce large increases in polarization. This might help explain the steep rises in polarization observed near the Lyman limit in some quasars.

We discuss a model of the central source in Seyfert 1 galaxy NGC 5548. The model assumes a three phase disk structure consisting of a cold outer disk, a hot central disk constituting a Comptonizing X/gamma source, and an intermediate unstable and complex phase emitting a soft excess component. The model qualitatively explains broad-band spectrum and variability behavior assuming that the soft excess contributes significantly to the continuum emission and drives variability by geometrical changes of the intermediate disk zone.

We analyze a composite broad-band optical/UV/Xgamma-ray spectrum of the Seyfert 1 galaxy NGC 5548. The spectrum consists of an average of simultaneous optical/IUE/Ginga observations accompanied by ROSAT and GRO/OSSE data from non-simultaneous observations. We show that the broad-band continuum is inconsistent with simple disk models extending to the soft X-rays. Instead, the soft-excess is well described by optically thick, low temperature, thermal Comptonization which may dominate the entire big blue bump. This might explain the observed tight UV/soft X-ray variability correlation and absence of a Lyman edge in this object. However, the plasma parameters inferred by the spectrum need stratification in optical depth and/or temperature to prevent physical inconsistency. The optical/UV/soft X-ray component contributes about half of the total source flux. The spectral variations of the soft-excess are consistent with that of the UV and argue that the components are closely related. The overall pattern of spectral variability suggests variations of the source geometry, and shows the optical/UV/soft X-ray component to be harder when brighter, while the hard X-ray component is softer when brighter.

Equipartition magnetic fields can dramatically affect the polarization of radiation emerging from accretion disk atmospheres in active galactic nuclei. We extend our previous work on this subject by exploring the interaction between Faraday rotation and absorption opacity in local, plane-parallel atmospheres with parameters appropriate for accretion discs. Faraday rotation in pure scattering atmospheres acts to depolarize the radiation field by rotating the polarization planes of photons after last scattering. Absorption opacity in an unmagnetized atmosphere can increase or decrease the polarization compared to the pure scattering case, depending on the thermal source function gradient. Combining both Faraday rotation and absorption opacity, we find the following results. If absorption opacity is much larger than scattering opacity throughout the atmosphere, then Faraday rotation generally has only a small effect on the emerging polarization because of the small electron column density along a photon mean free path. However, if the absorption opacity is not too large and it acts alone to increase the polarization, then the effects of Faraday rotation can be enhanced over those in a pure scattering atmosphere. Finally, while Faraday rotation often depolarizes the radiation field, it can in some cases increase the polarization when the thermal source function does not rise too steeply with optical depth. We confirm the correctness of the Silant'ev (1979) analytic calculation of the high magnetic field limit of the pure scattering atmosphere, which we incorrectly disputed in our previous paper.

We examine the possibility that gamma-ray bursts arise from sources in the Oort comet cloud, basing most of our arguments on accepted models for the formation and spatial distribution of the cloud. We identify three severe problems with such models: (1) There is no known mechanism for producing bursts that can explain the observed burst rate and energetics without violating other observational constraints. (2) The bright source counts cannot be reconciled with standard models for the phase-space distribution of objects in the Oort cloud. (3) The observed isotropy of the available burst data is inconsistent with the expected angular distribution of sources in the Oort cloud. We therefore assert that Oort cloud models of gamma-ray bursts are extremely implausible.