This white paper summarizes major scientific challenges and opportunities in understanding magnetic reconnection and related explosive phenomena as a fundamental plasma process.

Nonthermal relativistic plasmas are ubiquitous in astrophysical systems like pulsar wind nebulae and active galactic nuclei, as inferred from their emission spectra. The underlying nonthermal particle acceleration (NTPA) processes have traditionally been modeled with a Fokker-Planck (FP) diffusion-advection equation in momentum space. In this paper, we directly test the FP framework in ab-initio kinetic simulations of driven magnetized turbulence in relativistic pair plasma. By statistically analyzing the motion of tracked particles, we demonstrate the diffusive nature of NTPA and measure the FP energy diffusion ($D$) and advection ($A$) coefficients as functions of particle energy $\gamma m_e c^2$. We find that $D(\gamma)$ scales as $\gamma^2$ in the high-energy nonthermal tail, in line with 2nd-order Fermi acceleration theory, but has a much weaker scaling at lower energies. We also find that $A$ is not negligible and reduces NTPA by tending to pull particles towards the peak of the particle energy distribution. This study provides strong support for the FP picture of turbulent NTPA, thereby enhancing our understanding of space and astrophysical plasmas.

Rapid gamma-ray flares pose an astrophysical puzzle, requiring mechanisms both to accelerate energetic particles and to produce fast observed variability. These dual requirements may be satisfied by collisionless relativistic magnetic reconnection. On the one hand, relativistic reconnection can energize gamma-ray emitting electrons. On the other, as previous kinetic simulations have shown, the reconnection acceleration mechanism preferentially focuses high-energy particles -- and their emitted photons -- into beams, which may create rapid blips in flux as they cross a telescope's line of sight. Using a series of 2D pair-plasma particle-in-cell simulations, we explicitly demonstrate the critical role played by radiative cooling in mediating the observable signatures of this `kinetic beaming' effect. Only in our efficiently cooled simulations do we measure kinetic beaming beyond one light crossing time of the reconnection layer. We find a correlation between the cooling strength and the photon energy range across which persistent kinetic beaming occurs: stronger cooling coincides with a wider range of beamed photon energies. We also apply our results to rapid gamma-ray flares in flat-spectrum radio quasars, suggesting that a paradigm of radiatively efficient kinetic beaming constrains relevant emission models. In particular, beaming-produced variability may be more easily realized in two-zone (e.g. spine-sheath) setups, with Compton seed photons originating in the jet itself, rather than in one-zone external Compton scenarios.

High-energy astrophysical systems frequently contain collisionless relativistic plasmas that are heated by turbulent cascades and cooled by emission of radiation. Understanding the nature of this radiative turbulence is a frontier of extreme plasma astrophysics. In this paper, we use particle-in-cell simulations to study the effects of external inverse Compton radiation on turbulence driven in an optically thin, relativistic pair plasma. We focus on the statistical steady state (where injected energy is balanced by radiated energy) and perform a parameter scan spanning from low magnetization to high magnetization ($0.04 \lesssim \sigma \lesssim 11$). We demonstrate that the global particle energy distributions are quasi-thermal in all simulations, with only a modest population of nonthermal energetic particles (extending the tail by a factor of $\sim 2$). This indicates that nonthermal particle acceleration (observed in similar non-radiative simulations) is quenched by strong radiative cooling. The quasi-thermal energy distributions are well fit by analytic models in which stochastic particle acceleration (due to, e.g., second-order Fermi mechanism or gyroresonant interactions) is balanced by the radiation reaction force. Despite the efficient thermalization of the plasma, nonthermal energetic particles do make a conspicuous appearance in the anisotropy of the global momentum distribution as highly variable, intermittent beams (for high magnetization cases). The beamed high-energy particles are spatially coincident with intermittent current sheets, suggesting that localized magnetic reconnection may be a mechanism for kinetic beaming. This beaming phenomenon may explain rapid flares observed in various astrophysical systems (such as blazar jets, the Crab nebula, and Sagittarius A*).

We use global magnetohydrodynamic simulations to study the influence of net vertical magnetic fields on the structure of geometrically thin ($H/r \approx 0.05$) accretion disks in the Newtonian limit. We consider initial mid-plane gas to magnetic pressure ratios $\beta_0 = 1000,\, 300$ and $100$, spanning the transition between weakly and strongly magnetized accretion regimes. We find that magnetic pressure is important for the disks' vertical structure in all three cases, with accretion occurring at $z/R\approx 0.2$ in the two most strongly magnetized models. The disk midplane shows outflow rather than accretion. Accretion through the surface layers is driven mainly by stress due to coherent large scale magnetic field rather than by turbulent stress. Equivalent viscosity parameters measured from our simulations show similar dependencies on initial $\beta_0$ to those seen in shearing box simulations, though the disk midplane is not magnetic pressure dominated even for the strongest magnetic field case. Winds are present but are not the dominant driver of disk evolution. Over the (limited) duration of our simulations, we find evidence that the net flux attains a quasi-steady state at levels that can stably maintain a strongly magnetized disk. We suggest that geometrically thin accretion disks in observed systems may commonly exist in a magnetically "elevated" state, characterized by non-zero but modest vertical magnetic fluxes, with potentially important implications for disk phenomenology in X-ray binaries (XRBs) and active galactic nuclei (AGN).

This is a science white paper submitted to the Astro-2020 and Plasma-2020 Decadal Surveys. The paper describes the present status and emerging opportunities in Extreme Plasma Astrophysics -- a study of astrophysically-relevant plasma processes taking place under extreme conditions that necessitate taking into account relativistic, radiation, and QED effects.

This white paper highlights compact object and fundamental physics science opportunities afforded by high-throughput broadband (0.1-60 keV) X-ray polarization observations. X-ray polarimetry gives new observables with geometric information about stellar remnants which are many orders of magnitude too small for direct imaging. The X-ray polarimetric data also reveal details about the emission mechanisms and the structure of the magnetic fields in and around the most extreme objects in the Universe. Whereas the Imaging X-ray Polarimetry Explorer (IXPE) to be launched in 2021 will obtain first results for bright objects, a follow-up mission could be one order of magnitude more sensitive and would be able to use a broader bandpass to perform physics type experiments for representative samples of sources.

We present the Spectroscopic Time-Resolving Observatory for Broadband Energy X-rays (STROBE-X), a probe-class mission concept selected for study by NASA. It combines huge collecting area, high throughput, broad energy coverage, and excellent spectral and temporal resolution in a single facility. STROBE-X offers an enormous increase in sensitivity for X-ray spectral timing, extending these techniques to extragalactic targets for the first time. It is also an agile mission capable of rapid response to transient events, making it an essential X-ray partner facility in the era of time-domain, multi-wavelength, and multi-messenger astronomy. Optimized for study of the most extreme conditions found in the Universe, its key science objectives include: (1) Robustly measuring mass and spin and mapping inner accretion flows across the black hole mass spectrum, from compact stars to intermediate-mass objects to active galactic nuclei. (2) Mapping out the full mass-radius relation of neutron stars using an ensemble of nearly two dozen rotation-powered pulsars and accreting neutron stars, and hence measuring the equation of state for ultradense matter over a much wider range of densities than explored by NICER. (3) Identifying and studying X-ray counterparts (in the post-Swift era) for multiwavelength and multi-messenger transients in the dynamic sky through cross-correlation with gravitational wave interferometers, neutrino observatories, and high-cadence time-domain surveys in other electromagnetic bands. (4) Continuously surveying the dynamic X-ray sky with a large duty cycle and high time resolution to characterize the behavior of X-ray sources over an unprecedentedly vast range of time scales. STROBE-X's formidable capabilities will also enable a broad portfolio of additional science.

Electron and ion energization (i.e., heating and nonthermal acceleration) is a fundamental, but poorly understood, outcome of plasma turbulence. In this work, we present new results on this topic from particle-in-cell simulations of driven turbulence in collisionless, relativistic electron-ion plasma. We focus on temperatures such that ions (protons) are sub-relativistic and electrons are ultra-relativistic, a regime relevant for high-energy astrophysical systems such as hot accretion flows onto black holes. We find that ions tend to be preferentially heated, gaining up to an order of magnitude more energy than electrons, and propose a simple empirical formula to describe the electron-ion energy partition as a function of the ratio of electron-to-ion gyroradii (which in turn is a function of initial temperatures and plasma beta). We also find that while efficient nonthermal particle acceleration occurs for both species in the ultra-relativistic regime, nonthermal electron populations are diminished with decreasing temperature whereas nonthermal ion populations are essentially unchanged. These results have implications for modeling and interpreting observations of hot accretion flows.

We carry out an extensive linear stability analysis of magnetized cylindrical jets in a global framework. Foregoing the commonly invoked force-free limit, we focus on the small-scale, internal instabilities triggered in regions of the jet dominated by a toroidal magnetic field, with a weak vertical field and finite thermal pressure gradient. Such regions are likely to occur far from the jet source and boundaries, and are potential sites of magnetic energy dissipation that is essential to explain the particle acceleration and radiation observed from astrophysical jets. We validate the local stability analysis of Begelman by verifying that the eigenfunctions of the most unstable modes are radially localized. This finding allows us to propose a generic stability criterion in the presence of a weak vertical field. A stronger vertical field with a radial gradient complicates the stability criterion, due to the competition between the destabilizing thermal pressure gradient and stabilizing magnetic pressure gradients. Nevertheless, we argue that the jet interiors generically should be subject to rapidly growing, small-scale instabilities, capable of producing current sheets that lead to dissipation. We identify some new instabilities, not predicted by the local analysis, which are sensitive to the background radial profiles but have smaller growth rates than the local instabilities, and discuss the relevance of our work to the findings of recent numerical jet simulations.

Rapid, large amplitude variability at optical to X-ray wavelengths is now seen in an increasing number of Seyfert galaxies and luminous quasars. The variations imply a global change in accretion power, but are too rapid to be communicated by inflow through a standard thin accretion disc. Such discs are long known to have difficulty explaining the observed optical/UV emission from active galactic nuclei. Here we show that alternative models developed to explain these observations have larger scale heights and shorter inflow times. Accretion discs supported by magnetic pressure in particular are geometrically thick at all luminosities, with inflow times as short as the observed few year timescales in extreme variability events to date. Future time-resolved, multi-wavelength observations can distinguish between inflow through a geometrically thick disc as proposed here, and alternative scenarios of extreme reprocessing of a central source or instability-driven limit cycles.

We have modeled direct collapse of a primordial gas within dark matter halos in the presence of radiative transfer, in high-resolution zoom-in simulations in a cosmological framework, down to the formation of the photosphere and the central object. Radiative transfer has been implemented in the flux-limited diffusion (FLD) approximation. Adiabatic models were run for comparison. We find that (a) the FLD flow forms an irregular central structure and does not exhibit fragmentation, contrary to adiabatic flow which forms a thick disk, driving a pair of spiral shocks, subject to Kelvin-Helmholtz shear instability forming fragments; (b) the growing central core in the FLD flow quickly reaches ~10 Mo and a highly variable luminosity of 10^{38}-10^{39} erg/s, comparable to the Eddington luminosity. It experiences massive recurrent outflows driven by radiation force and thermal pressure gradients, which mix with the accretion flow and transfer the angular momentum outwards; and (c) the interplay between these processes and a massive accretion, results in photosphere at ~10 AU. We conclude that in the FLD model (1) the central object exhibits dynamically insignificant rotation and slower than adiabatic temperature rise with density; (2) does not experience fragmentation leading to star formation, thus promoting the fast track formation of a supermassive black hole (SMBH) seed; (3) inclusion of radiation force leads to outflows, resulting in the mass accumulation within the central 10^{-3} pc, which is ~100 times larger than characteristic scale of star formation. The inclusion of radiative transfer reveals complex early stages of formation and growth of the central structure in the direct collapse scenario of SMBH seed formation.

We apply collisionless particle-in-cell simulations of relativistic pair plasmas to explore whether driven turbulence is a viable high-energy astrophysical particle accelerator. We characterize nonthermal particle distributions for varying system sizes up to $L/2\pi\rho_{e0} = 163$, where $L/2\pi$ is the driving scale and $\rho_{e0}$ is the initial characteristic Larmor radius. We show that turbulent particle acceleration produces power-law energy distributions that, when compared at a fixed number of large-scale dynamical times, slowly steepen with increasing system size. We demonstrate, however, that convergence is obtained by comparing the distributions at different times that increase with system size (approximately logarithmically). We suggest that the system-size dependence arises from the time required for particles to reach the highest accessible energies via Fermi acceleration. The converged power-law index of the energy distribution, $\alpha \approx 3.0$ for magnetization $\sigma = 3/8$, makes turbulence a possible explanation for nonthermal spectra observed in systems such as the Crab nebula.

Direct collapse within dark matter (DM) halos is a promising path to form supermassive black hole (SMBH) seeds at high redshifts. The outer part of this collapse remains optically thin, and has been studied intensively using numerical simulations. However, the innermost region of the collapse is expected to become optically thick and requires us to follow the radiation field in order to understand its subsequent evolution. So far, the adiabatic approximation has been used exclusively for this purpose. We apply radiative transfer in the flux-limited diffusion (FLD) approximation to solve the evolution of coupled gas and radiation, for isolated halos. For direct collapse within isolated DM halos, we find that (1) the photosphere forms at ~10^{-6} pc and rapidly expands outward. (2) A central core forms, with a mass of ~1 Mo, supported by thermal gas pressure gradients and rotation. (3) Growing thermal gas and radiation pressure gradients dissolve it. (4) This process is associated with a strong anisotropic outflow, and another core forms nearby and grows rapidly. (5) Typical radiation luminosity emerging from the photosphere encompassing these cores is ~5 x 10^{37}-5 x 10^{38} erg/s, of order the Eddington luminosity. (6) Two variability timescales are associated with this process: a long one, which is related to the accretion flow within the central ~10^{-4}-10^{-3} pc, and ~0.1 yr, which is related to radiation diffusion. (7) Adiabatic models have been run for comparison and their evolution differs profoundly from that of the FLD models, by forming a central geometrically-thick disk. Overall, an adiabatic equation of state is not a good approximation to the advanced stage of direct collapse, mainly because the radiation in the FLD is capable of escaping due to anisotropy in the optical depth and associated gradients.

Magnetic reconnection in relativistic collisionless plasmas can accelerate particles and power high-energy emission in various astrophysical systems. Whereas most previous studies focused on relativistic reconnection in pair plasmas, less attention has been paid to electron-ion plasma reconnection, expected in black hole accretion flows and relativistic jets. We report a comprehensive particle-in-cell numerical investigation of reconnection in an electron-ion plasma, spanning a wide range of ambient ion magnetizations $\sigma_i$, from the semirelativistic regime (ultrarelativistic electrons but nonrelativistic ions, 0.001<<$\sigma_i$<<1) to the fully relativistic regime (both species are ultrarelativistic, $\sigma_i$>>1). We investigate how the reconnection rate, electron and ion plasma flows, electric and magnetic field structures, electron/ion energy partitioning, and nonthermal particle acceleration depend on $\sigma_i$. Our key findings are: (1) the reconnection rate is about 0.1 of the Alfvenic rate across all regimes; (2) electrons can form concentrated moderately relativistic outflows even in the semirelativistic, small-$\sigma_i$ regime; (3) while the released magnetic energy is partitioned equally between electrons and ions in the ultrarelativistic limit, the electron energy fraction declines gradually with decreased $\sigma_i$ and asymptotes to about 0.25 in the semirelativistic regime; (4) reconnection leads to efficient nonthermal electron acceleration with a $\sigma_i$-dependent power-law index, $p(\sigma_i) \simeq $const$+0.7 {\sigma_i}^{-1/2}$. These findings are important for understanding black hole systems and lend support to semirelativistic reconnection models for powering nonthermal emission in blazar jets, offering a natural explanation for the spectral indices observed in these systems.

We examine the properties of strongly magnetized accretion discs in a global framework, with particular focus on the evolution of magnetohydrodynamic instabilities such as the magnetorotational instability (MRI). Work by Pessah and Psaltis showed that MRI is stabilized beyond a critical toroidal field in compressible, differentially rotating flows and, also, reported the appearance of two new instabilities beyond this field. Their results stemmed from considering geometric curvature effects due to the suprathermal background toroidal field, which had been previously ignored in weak-field studies. However, their calculations were performed under the local approximation, which poses the danger of introducing spurious behavior due to the introduction of global geometric terms in an otherwise local framework. In order to avoid this, we perform a global eigenvalue analysis of the linearized MHD equations in cylindrical geometry. We confirm that MRI indeed tends to be highly suppressed when the background toroidal field attains the Pessah-Psaltis limit. We also observe the appearance of two new instabilities that emerge in the presence of highly suprathermal toroidal fields. These results were additionally verified using numerical simulations in PLUTO. There are, however, certain differences between the the local and global results, especially in the vertical wavenumber occupancies of the various instabilities, which we discuss in detail. We also study the global eigenfunctions of the most unstable modes in the suprathermal regime, which are inaccessible in the local analysis. Overall, our findings emphasize the necessity of a global treatment for accurately modeling strongly magnetized accretion discs.

We describe results from particle-in-cell simulations of driven turbulence in collisionless, magnetized, relativistic pair plasma. This physical regime provides a simple setting for investigating the basic properties of kinetic turbulence and is relevant for high-energy astrophysical systems such as pulsar wind nebulae and astrophysical jets. In this paper, we investigate the statistics of turbulent fluctuations in simulations on lattices of up to $1024^3$ cells and containing up to $2 \times 10^{11}$ particles. Due to the absence of a cooling mechanism in our simulations, turbulent energy dissipation reduces the magnetization parameter to order unity within a few dynamical times, causing turbulent motions to become sub-relativistic. In the developed stage, our results agree with predictions from magnetohydrodynamic turbulence phenomenology at inertial-range scales, including a power-law magnetic energy spectrum with index near $-5/3$, scale-dependent anisotropy of fluctuations described by critical balance, log-normal distributions for particle density and internal energy density (related by a $4/3$ adiabatic index, as predicted for an ultra-relativistic ideal gas), and the presence of intermittency. We also present possible signatures of a kinetic cascade by measuring power-law spectra for the magnetic, electric, and density fluctuations at sub-Larmor scales.

We present results from particle-in-cell simulations of driven turbulence in magnetized, collisionless, and relativistic pair plasma. We find that fluctuations are consistent with the classical $k_\perp^{-5/3}$ magnetic energy spectrum at fluid scales and a steeper $k_\perp^{-4}$ spectrum at sub-Larmor scales, where $k_\perp$ is the wavevector perpendicular to the mean field. We demonstrate the development of a non-thermal, power-law particle energy distribution, $f(E) \sim E^{-\alpha}$, with index $\alpha$ that decreases with increasing magnetization and increases with increasing system size (relative to the characteristic Larmor radius). Our simulations indicate that turbulence can be a viable source of energetic particles in high-energy astrophysical systems, such as pulsar wind nebulae, if scalings asymptotically become insensitive to the system size.

We investigate the pre-disruption gravitational dynamics and post-disruption hydrodynamics of the tidal disruption of stars by supermassive black hole (SMBH) binaries. We focus on binaries with relatively low mass primaries ($10^6M_{\odot}$), moderate mass ratios, and separations with reasonably long gravitational wave inspiral times (tens of Myr). First, we generate a large ensemble (between 1 and 10 million) of restricted three-body integrations to quantify the statistical properties of tidal disruptions by circular SMBH binaries of initially-unbound stars. Compared to the reference case of a disruption by a single SMBH, the binary potential induces significant variance into the specific energy and angular momentum of the star at the point of disruption. Second, we use Newtonian numerical hydrodynamics to study the detailed evolution of the fallback debris from 120 disruptions randomly selected from the three-body ensemble (excluding only the most deeply penetrating encounters). We find that the overall morphology of the debris is greatly altered by the presence of the second black hole, and the accretion rate histories display a wide range of behaviors, including order of magnitude dips and excesses relative to control simulations that include only one black hole. Complex evolution persists, in some cases, for many orbital periods of the binary. We find evidence for power in the accretion curves on timescales related to the binary orbital period, though there is no exact periodicity. We discuss our results in the context of future wide-field surveys, and comment on the prospects of identifying and characterizing the subset of events occurring in nuclei with binary SMBHs.

The power mechanism and accretion geometry for low-power FR1 radio galaxies is poorly understood in comparison to Seyfert galaxies and QSOs. In this paper we use the diagnostic power of the Lya recombination line observed using the Cosmic Origins Spectrograph (COS) aboard the Hubble Space Telescope (HST) to investigate the accretion flows in three well-known, nearby FR1s: M87, NGC4696, and HydraA. The Ly$\alpha$ emission line's luminosity, velocity structure and the limited knowledge of its spatial extent provided by COS are used to assess conditions within a few parsecs of the super-massive black hole (SMBH) in these radio-mode AGN. We observe strong Ly$\alpha$ emission in all three objects with similar total luminosity to that seen in BL Lacertae objects. M87 shows a complicated emission line profile in Lya which varies spatially across the COS aperture and possibly temporally over several epochs of observation. In both NGC 4696 and M87, the Ly$\alpha$ luminosities $\sim 10^{40}$ ergs/s are closely consistent with the observed strength of the ionizing continuum in Case B recombination theory and with the assumption of near unity covering factor. It is possible that the Ly$\alpha$ emitting clouds are ionized largely by beamed radiation associated with the jets. Long-slit UV spectroscopy can be used to test this hypothesis. Hydra A and the several BL Lac objects studied in this and previous papers have Lya luminosities larger than M87 but their extrapolated, non-thermal continua are so luminous that they over-predict the observed strength of Ly$\alpha$, a clear indicator of relativistic beaming in our direction. Given their substantial space density ($\sim 4\times10^{-3} Mpc^{-3}$) the unbeamed Lyman continuum radiation of FR1s may make a substantial minority contribution (~10%) to the local UV background if all FR1s are similar to M87 in ionizing flux level.

We propose that the accretion disks fueling active galactic nuclei are supported vertically against gravity by a strong toroidal ($\phi-$direction) magnetic field that develops naturally as the result of an accretion disk dynamo. The magnetic pressure elevates most of the gas carrying the accretion flow at $R$ to large heights $z > 0.1 R$ and low densities, while leaving a thin dense layer containing most of the mass --- but contributing very little accretion --- around the equator. We show that such a disk model leads naturally to the formation of a broad emission line region through thermal instability. Extrapolating to larger radii, we demonstrate that local gravitational instability and associated star formation are strongly suppressed compared to standard disk models for AGN, although star formation in the equatorial zone is predicted for sufficiently high mass supply rates. This new class of accretion disk models thus appears capable of resolving two longstanding puzzles in the theory of AGN fueling: the formation of broad emission line regions and the suppression of fragmentation thought to inhibit accretion at the required rates. We show that the disk of stars that formed in the Galactic Center a few million years ago could have resulted from an episode of magnetically elevated accretion at $> 0.1$ of the Eddington limit.

The rate at which matter flows into a galactic nucleus during early phases of galaxy evolution can sometimes exceed the Eddington limit of the growing central black hole by several orders of magnitude. We discuss the necessary conditions for the black hole to actually accrete this matter at such a high rate, and consider the observational appearance and detectability of a hyperaccreting black hole. In order to be accreted at a hyper-Eddington rate, the infalling gas must have a sufficiently low angular momentum. Although most of the gas is accreted, a significant fraction accumulates in an optically thick envelope with luminosity $\sim L_{\rm Edd}$, probably pierced by jets of much higher power. If $\dot M > 10^3 \dot M_{\rm Edd}$, the envelope spectrum resembles a blackbody with a temperature of a few thousand K, but for lower (but still hyper-Eddington) accretion rates the spectrum becomes a very dilute and hard Wien spectrum. We consider the likelihood of various regimes of hyperaccretion, and discuss its possible observational signatures.

We have investigated the complex multiwavelength evolution of GRO J1655-40 during the rise of its 2005 outburst. We detected two hard X-ray flares, the first one during the transition from the soft state to the ultra-soft state, and the second one in the ultra-soft state. The first X-ray flare coincided with an optically thin radio flare. We also observed a hint of increased radio emission during the second X-ray flare. To explain the hard flares without invoking a secondary emission component, we fit the entire data set with the eqpair model. This single, hybrid Comptonization model sufficiently fits the data even during the hard X-ray flares if we allow reflection fractions greater than unity. In this case, the hard X-ray flares correspond to a Comptonizing corona dominated by non-thermal electrons. The fits also require absorption features in the soft and ultra-soft state which are likely due to a wind. In this work we show that the wind and the optically thin radio flare co-exist. Finally, we have also investigated the radio to optical spectral energy distribution, tracking the radio spectral evolution through the quenching of the compact jet and rise of the optically thin flare, and interpreted all data using state transition models.

We present a method for analyzing the interaction between radiation and matter in regions of intense, relativistic shear that can arise in many astrophysical situations. We show that there is a simple velocity profile that should be manifested in regions of large shear that have "lost memory" of their boundary conditions, and we use this self-similar velocity profile to construct the surface of last scattering, or $\tau \simeq 1$ surface, as viewed from any comoving point within the flow. We demonstrate that a simple treatment of scattering from this $\tau \simeq 1$ surface exactly conserves photon number, and derive the rate at which the radiation field is heated due to the shear present in the flow. The components of the comoving radiation energy-momentum tensor are calculated, and we show that they have relatively simple, approximate forms that interpolate between the viscous (small shear) and streaming (large shear) limits. We put our expression for the energy-momentum tensor in a covariant form that does not depend on the explicit velocity profile within the fluid and, therefore, represents a natural means for analyzing general, radiation-dominated, relativistic shear flows.

Motivated by indirect observational evidence for strongly magnetized accretion discs around black holes, and the novel theoretical properties of such solutions, we investigate how a strong magnetization state can develop and persist. To this end, we perform local simulations of accretion discs with an initially purely toroidal magnetic field of equipartition strength. We demonstrate that discs with zero net vertical magnetic flux and realistic boundary conditions cannot sustain a strong toroidal field. However, a magnetic pressure-dominated disc can form from an initial configuration with a sufficient amount of net vertical flux and realistic boundary conditions. Our results suggest that poloidal flux is a necessary prerequisite for the sustainability of strongly magnetized accretion discs.

A tidal disruption event (TDE) -- when a star is destroyed by the immense gravitational field of a supermassive black hole -- transforms a star into a stream of tidally-shredded debris. The properties of this debris ultimately determine the observable signatures of TDEs. Here we derive a simple, self-similar solution for the velocity profile of the debris streams produced from TDEs, and show that this solution agrees extremely well with numerical results. Using this self-similar solution, we calculate an analytic, approximate expression for the radial density profile of the stream. We show that there is a critical adiabatic index that varies as a function of position along the stream above (below) which the stream is unstable (stable) to gravitational fragmentation. We also calculate the impact of heating and cooling on this stability criterion.

We have developed a quasi-analytical model for the production of radiation in strong-line blazars, assuming a spine-sheath jet structure. The model allows us to study how the spine and sheath spectral components depend on parameters describing the geometrical and physical structure of "the blazar zone". We show that typical broad-band spectra of strong-line blazars can be reproduced by assuming the magnetization parameter to be of order unity and reconnection to be the dominant dissipation mechanism. Furthermore, we demonstrate that the spine-sheath model can explain why gamma-ray variations are often observed to have much larger amplitudes than the corresponding optical variations. The model is also less demanding of jet power than one-zone models, and can reproduce the basic features of extreme gamma-ray events.

Strongly magnetized accretion discs around black holes have attractive features that may explain enigmatic aspects of X-ray binary behaviour. The structure and evolution of these discs are governed by a dynamo-like mechanism, which channels part of the accretion power liberated by the magnetorotational instability (MRI) into an ordered toroidal magnetic field. To study dynamo activity, we performed three-dimensional, stratified, isothermal, ideal magnetohydrodynamic shearing box simulations. The strength of the self-sustained toroidal magnetic field depends on the net vertical magnetic flux, which we vary across almost the entire range over which the MRI is linearly unstable. We quantify disc structure and dynamo properties as a function of the initial ratio of mid-plane gas pressure to vertical magnetic field pressure, $\beta_0^{\rm mid} = p_{\rm gas} / p_B$. For $10^5 \geq \beta_0^{\rm mid} \geq 10$ the effective $\alpha$-viscosity parameter scales as a power-law. Dynamo activity persists up to and including $\beta_0^{\rm mid} = 10^2$, at which point the entire vertical column of the disc is magnetic pressure-dominated. Still stronger fields result in a highly inhomogeneous disc structure, with large density fluctuations. We show that the turbulent steady state $\beta^{\rm mid}$ in our simulations is well-matched by the analytic model of Begelman et al. (2015) describing the creation and buoyant escape of toroidal field, while the vertical structure of the disc can be broadly reproduced using this model. Finally, we discuss the implications of our results for observed properties of X-ray binaries.

Using two-dimensional particle-in-cell simulations, we characterize the energy spectra of particles accelerated by relativistic magnetic reconnection (without guide field) in collisionless electron-positron plasmas, for a wide range of upstream magnetizations $\sigma$ and system sizes $L$. The particle spectra are well-represented by a power law $\gamma^{-\alpha}$, with a combination of exponential and super-exponential high-energy cutoffs, proportional to $\sigma$ and $L$, respectively. For large $L$ and $\sigma$, the power-law index $\alpha$ approaches about 1.2.

We use cosmological adaptive mesh refinement (AMR) code Enzo zoom-in simulations to study the long term evolution of the collapsing gas within dark matter (DM) halos at high redshifts. This direct collapse process is a leading candidate for rapid formation of supermassive black hole (SMBH) seeds at high z. To circumvent the Courant condition at small radii, we have used the sink particle method, and focus on the evolution on scales ~0.01-10 pc. The collapse proceeds in two stages, with the secondary runaway happening within the central 10 pc, and with no detected fragmentation. The sink particles form when the collapsing gas requires additional refinement of the grid size at the highest refinement level. Their mass never exceeds ~10^3 Mo, with the sole exception of the central seed which grows dramatically to ~ 2 x 10^6 Mo in ~2 Myr, confirming the feasibility of this path to the SMBH. The time variability of angular momentum axis in the accreted gas results in the formation of two misaligned disks --- a small inner disk, and a more massive, outer disk which is inclined by ~45^o to the inner disk. The self-gravity of these disks is heavily diluted --- both disks lie within the Roche limit of the central seed. While the inner disk is geometrically thin and weakly asymmetric, the outer disk flares due to turbulent motions as a result of the massive inflow along a pair of penetrating filaments. The geometry of inflow via filaments determines the dominant and secondary Fourier modes in this disk --- these modes have a non-self-gravitational origin. We do not confirm that m=1 is a principal mode that drives the inflow in the presence of a central massive object. While the positions of the disks depend on the scale chosen to break the self-similar collapse, the overall configuration appears to be generic, and is expected to form when the central seed becomes sufficiently massive.

We investigate the distribution of particle acceleration sites, independently of the actual acceleration mechanism, during plasmoid-dominated, relativistic collisionless magnetic reconnection by analyzing the results of a particle-in-cell numerical simulation. The simulation is initiated with Harris-type current layers in pair plasma with no guide magnetic field, negligible radiative losses, no initial perturbation, and using periodic boundary conditions. We find that the plasmoids develop a robust internal structure, with colder dense cores and hotter outer shells, that is recovered after each plasmoid merger on a dynamical time scale. We use spacetime diagrams of the reconnection layers to probe the evolution of plasmoids, and in this context we investigate the individual particle histories for a representative sample of energetic electrons. We distinguish three classes of particle acceleration sites associated with (1) magnetic X-points, (2) regions between merging plasmoids, and (3) the trailing edges of accelerating plasmoids. We evaluate the contribution of each class of acceleration sites to the final energy distribution of energetic electrons -- magnetic X-points dominate at moderate energies, and the regions between merging plasmoids dominate at higher energies. We also identify the dominant acceleration scenarios, in order of decreasing importance -- (1) single acceleration between merging plasmoids, (2) single acceleration at a magnetic X-point, and (3) acceleration at a magnetic X-point followed by acceleration in a plasmoid. Particle acceleration is absent only in the vicinity of stationary plasmoids. The effect of magnetic mirrors due to plasmoid contraction does not appear to be significant in relativistic reconnection.

A tidal disruption event, which occurs when a star is destroyed by the gravitational field of a supermassive black hole, produces a stream of debris, the evolution of which ultimately determines the observational properties of the event. Here we show that a post-periapsis caustic -- a location where the locus of gas parcels comprising the stream would collapse into a two-dimensional surface if they evolved solely in the gravitational field of the hole -- occurs when the pericenter distance of the star is on the order of the tidal radius of the hole. It is demonstrated that this "pancake" induces significant density perturbations in the debris stream, and, for stiffer equations of state (adiabatic index $\gamma \gtrsim 5/3$), these fluctuations are sufficient to gravitationally destabilize the stream, resulting in its fragmentation into bound clumps. The results of our findings are discussed in the context of the observational properties of tidal disruption events.

Magnetohydrodynamic accretion disk simulations suggest that much of the energy liberated by the magnetorotational instability (MRI) can be channeled into large-scale toroidal magnetic fields through dynamo action. Under certain conditions, this field can dominate over gas and radiation pressure in providing vertical support against gravity, even close to the midplane. Using a simple model for the creation of this field, its buoyant rise, and its coupling to the gas, we show how disks could be driven into this magnetically dominated state and deduce the resulting vertical pressure and density profiles. Applying an established criterion for MRI to operate in the presence of a toroidal field, we show that magnetically supported disks can have two distinct MRI-active regions, separated by a "dead zone" where local MRI is suppressed, but where magnetic energy continues to flow upward from the dynamo region below. We suggest that the relative strengths of the MRI zones, and the local poloidal flux, determine the spectral states of X-ray binaries. Specifically, "intermediate" and "hard" accretion states occur when MRI is triggered in the hot, upper zone of the corona, while disks in "soft" states do not develop the upper MRI zone. We discuss the conditions under which various transitions should take place and speculate on the relationship of dynamo activity to the various types of quasi-periodic oscillations that sometimes appear in the hard spectral components. The model also explains why luminous accretion disks in the "soft" state show no signs of the thermal/viscous instability predicted by standard alpha models.

We analyze the interaction of a radiation-dominated jet and its surroundings using the equations of radiation hydrodynamics in the viscous limit. In a previous paper we considered the two-stream scenario, which treats the jet and its surroundings as distinct media interacting through radiation viscous forces. Here we present an alternative boundary layer model, known as the free-streaming jet model -- where a narrow stream of fluid is injected into a static medium -- and present solutions where the flow is ultrarelativistic and the boundary layer is dominated by radiation. It is shown that these jets entrain material from their surroundings and that their cores have a lower density of scatterers and a harder spectrum of photons, leading to observational consequences for lines of sight that look "down the barrel of the jet." These jetted outflow models may be applicable to the jets produced during long gamma-ray bursts and super-Eddington phases of tidal disruption events.

Using the relativistic equations of radiation hydrodynamics in the viscous limit, we analyze the boundary layers that develop between radiation-dominated jets and their environments. In this paper we present the solution for the self-similar, 2-D, plane-parallel two-stream problem, wherein the jet and the ambient medium are considered to be separate, interacting fluids, and we compare our results to those of previous authors. (In a companion paper we investigate an alternative scenario, known as the free-streaming jet model.) Consistent with past findings, we show that the boundary layer that develops between the jet and its surroundings creates a region of low-density material. These models may be applicable to sources such as super-Eddington tidal disruption events and long gamma-ray bursts.

We study the early stage of the formation of seed supermassive black holes via direct collapse in dark matter (DM) halos, in the cosmological context. We perform high-resolution zoom-in simulations of such collapse at high-$z$. Using the adaptive mesh refinement code ENZO, we resolve the formation and growth of a DM halo, until its virial temperature reaches $\sim 10^4$K, atomic cooling turns on, and collapse ensues. We demonstrate that direct collapse proceeds in two stages, although they are not well separated. The first stage is triggered by the onset of atomic cooling, and leads to rapidly increasing accretion rate with radius, from $\dot M\sim 0.1\,M_\odot {\rm yr^{-1}}$ at the halo virial radius to few $M_\odot \,{\rm yr^{-1}}$, around the scale radius $R_{\rm s}\sim 30$pc of the NFW DM density profile. The second stage of the collapse commences when the gas density takes precedence over the DM density. This is associated with the gas decoupling from the DM gravitational potential. The ensuing collapse approximates that of an isothermal sphere with $\dot M ( r )\sim $const. We confirm that the gas loses its angular momentum through non-axisymmetric perturbations and gravitational torques, to overcome the centrifugal barrier. During the course of the collapse, this angular momentum transfer process happens on nearly all spatial scales, and the angular momentum vector of the gas varies with position and time. Collapsing gas also exhibits supersonic turbulent motions which suppress gas fragmentation, and are characterized by density PDF consisting of a lognormal part and a high-density power law tail.

This is a White Paper in support of the mission concept of the Large Observatory for X-ray Timing (LOFT), proposed as a medium-sized ESA mission. We discuss the potential of LOFT for the study of the physics of accretion and ejection around compact objects. For a summary, we refer to the paper.

We present an analysis of the general relativistic Boltzmann equation for radiation, appropriate to the case where particles and photons interact through Thomson scattering, and derive the radiation energy-momentum tensor in the diffusion limit, with viscous terms included. Contrary to relativistic generalizations of the viscous stress tensor that appear in the literature, we find that the stress tensor should contain a correction to the comoving energy density proportional to the divergence of the four-velocity, as well as a finite bulk viscosity. These modifications are consistent with the framework of radiation hydrodynamics in the limit of large optical depth, and do not depend on thermodynamic arguments such as the assignment of a temperature to the zeroth-order photon distribution. We perform a perturbation analysis on our equations and demonstrate that, as long as the wave numbers do not probe scales smaller than the mean free path of the radiation, the viscosity contributes only decaying, i.e., stable, corrections to the dispersion relations. The astrophysical applications of our equations, including jets launched from super-Eddington tidal disruption events and those from collapsars, are discussed and will be considered further in future papers.

Hot relativistic jets, passing through a background medium with a pressure gradient p \propto r^{-\eta} where 2 < \eta <= 8/3, develop a shocked boundary layer containing a significant fraction of the jet power. In previous work, we developed a self-similar description of the boundary layer assuming isentropic flow, but we found that such models respect global energy conservation only for the special case \eta = 8/3. Here we demonstrate that models with \eta < 8/3 can be made self-consistent if we relax the assumption of constant specific entropy. Instead, the entropy must increase with increasing r along the boundary layer, presumably due to multiple shocks driven into the flow as it gradually collimates. The increase in specific entropy slows the acceleration rate of the flow and provides a source of internal energy that could be channeled into radiation. We suggest that this process may be important for determining the radiative characteristics of tidal disruption events and gamma-ray bursts from collapsars.

I outline the theory of accretion onto black holes, and its application to observed phenomena such as X-ray binaries, active galactic nuclei, tidal disruption events, and gamma-ray bursts. The dynamics as well as radiative signatures of black hole accretion depend on interactions between the relatively simple black-hole spacetime and complex radiation, plasma and magnetohydrodynamical processes in the surrounding gas. I will show how transient accretion processes could provide clues to these interactions. Larger global magnetohydrodynamic simulations as well as simulations incorporating plasma microphysics and full radiation hydrodynamics will be needed to unravel some of the current mysteries of black hole accretion.

Estimates of magnetic field strength in relativistic jets of active galactic nuclei (AGN), obtained by measuring the frequency-dependent radio core location, imply that the total magnetic fluxes in those jets are consistent with the predictions of the magnetically-arrested disk (MAD) scenario of jet formation. On the other hand, the magnetic field strength determines the luminosity of the synchrotron radiation, which forms the low-energy bump of the observed blazar spectral energy distribution (SED). The SEDs of the most powerful blazars are strongly dominated by the high-energy bump, which is most likely due to the external radiation Compton (ERC) mechanism. This high Compton dominance may be difficult to reconcile with the MAD scenario, unless 1) the geometry of external radiation sources (broad-line region, hot-dust torus) is quasi-spherical rather than flat, or 2) most gamma-ray radiation is produced in jet regions of low magnetization, e.g., in magnetic reconnection layers or in fast jet spines.

[abridged] Here we present the analysis of multifrequency data gathered for the FRII radio galaxy PKS B1358-113, hosted in the brightest cluster galaxy of Abell 1836. The galaxy harbors one of the most massive black holes known to date and our analysis of the optical data reveals that this black hole is only weakly active. Based on new Chandra and XMM-Newton X-ray observations and archival radio data we derive the preferred range for the jet kinetic luminosity $\sim (0.5-3) \times 10^{45}$ erg s$^{-1}$. This is above the values implied by various scaling relations proposed for radio sources in galaxy clusters, being instead very close to the maximum jet power allowed for the given accretion rate. We constrain the radio source lifetime as $\sim 40-70$ Myrs, and the total amount of deposited jet energy $\sim (2-8) \times 10^{60}$\,ergs. The detailed analysis of the X-ray data provides indication for the presence of a bow-shock driven by the expanding radio lobes into the Abell 1836 cluster environment, with the corresponding Mach number $\sim 2-4$. This, together with the recently growing evidence that powerful FRII radio galaxies may not be uncommon in the centers of clusters at higher redshifts, supports the idea that jet-induced shock heating may indeed play an important role in shaping the properties of clusters, galaxy groups, and galaxies in formation. We speculate on a possible bias against detecting jet-driven shocks in poorer environments, resulting from an inefficient electron heating at the shock front, combined with a relatively long electron-ion equilibration timescale.

Locating the gamma-ray emission sites in blazar jets is a long-standing and highly controversial issue. We investigate jointly several constraints on the distance scale r and Lorentz factor Gamma of the gamma-ray emitting regions in luminous blazars (primarily flat spectrum radio quasars, FSRQs). Working in the framework of one-zone external radiation Comptonization (ERC) models, we perform a parameter space study for several representative cases of actual gamma-ray flares in their multiwavelength context. We find a particularly useful combination of three constraints: from an upper limit on the collimation parameter Gamma*theta <~ 1, from an upper limit on the synchrotron self-Compton (SSC) luminosity L_SSC <~ L_X, and from an upper limit on the efficient cooling photon energy E_cool,obs <~ 100 MeV. These three constraints are particularly strong for sources with low accretion disk luminosity L_d. The commonly used intrinsic pair-production opacity constraint on Gamma is usually much weaker than the SSC constraint. The SSC and cooling constraints provide a robust lower limit on the collimation parameter Gamma*theta >~ 0.1 - 0.7. Typical values of r corresponding to moderate values of Gamma ~ 20 are in the range 0.1 - 1 pc, and are determined primarily by the observed variability time scale t_var,obs. Alternative scenarios motivated by the observed gamma-ray/mm connection, in which gamma-ray flares of t_var,obs ~ a few days are located at r ~ 10 pc, are in conflict with both the SSC and cooling constraints. Moreover, we use a simple light travel time argument to point out that the gamma-ray/mm connection does not provide a significant constraint on the location of gamma-ray flares. We argue that spine-sheath models of the jet structure do not offer a plausible alternative to external radiation fields at large distances, however, an extended broad-line region is an idea worth exploring.

The absence of thermal instability in the high/soft state of black hole X-ray binaries, in disagreement with the standard thin disk theory, is a long-standing riddle for theoretical astronomers. We have tried to resolve this question by studying the thermal stability of a thin disk with magnetically driven winds in the $\dot{M}- \Sigma$ plane. It is found that disk winds can greatly decrease the disk temperature and thus help the disk become more stable at a given accretion rate. The critical accretion rate $\dot{M}_{\rm crit}$ corresponding to the thermal instability threshold is increased significantly in the presence of disk winds. For $\alpha=0.01$ and $B_{\rm \phi}=10B_{\rm _p}$, the disk is quite stable even for a very weak initial poloidal magnetic field [$\beta_{\rm p,0}\sim 2000, \beta_{\rm p}=(P_{\rm {gas}}+P_{\rm rad})/(B_{\rm p}^2/8\pi)$]. But when $B_{\rm \phi}=B_{\rm _p}$ or $B_{\rm \phi}=0.1B_{\rm _p}$, a somewhat stronger (but still weak) field ($\beta_{\rm p,0}\sim 200$ or $\beta_{\rm p,0}\sim 20$) is required to make the disk stable. Nevertheless, despite the great increase of $\dot{M}_{\rm crit}$, the luminosity threshold corresponding to instability remains almost constant or even decreases slowly with increasing $\dot{M}_{\rm crit}$ due to the decrease of gas temperature. The advection and diffusion timescales of the large-scale magnetic field threading the disk are also investigated in this work. We find that the advection timescale can be smaller than the diffusion timescale in a disk with winds, because the disk winds take away most of the gravitational energy released in the disk, resulting in the decrease of the magnetic diffusivity $\eta$ and the increase of the diffusion timescale.

We suggest that the hysteretic cycle of black hole state transitions arises from two established properties of accretion disks: the increase in turbulent stress in disks threaded by a net magnetic field and the ability of thick (but not thin) disks to advect such a field radially. During quiescence, magnetic field loops are generated by the magnetorotational instability at the interface between the inner hot flow and outer thin disk. Vertical flux is advected into and accumulates stochastically within the inner flow, where it stimulates the turbulence so that $\alpha \sim 1$. The transition to a geometrically thin inner disk occurs when $L \sim \alpha^2 L_{\rm Edd} \sim L_{\rm Edd}$, and the first "thin" disk to form is itself moderately thick, strongly magnetized, and able to advect field inward. These properties favor episodic jet production. As the accretion rate declines magnetic flux escapes, $\alpha$ decreases to $\alpha \sim 0.01 - 0.1$, and a hot inner flow is not re-established until $L \ll L_{\rm Edd}$. We discuss possible observational consequences of our scenario.

The discovery of rapid synchrotron gamma-ray flares above 100 MeV from the Crab Nebula has attracted new interest in alternative particle acceleration mechanisms in pulsar wind nebulae. Diffuse shock-acceleration fails to explain the flares because particle acceleration and emission occur during a single or even sub-Larmor timescale. In this regime, the synchrotron energy losses induce a drag force on the particle motion that balances the electric acceleration and prevents the emission of synchrotron radiation above 160 MeV. Previous analytical studies and 2D particle-in-cell (PIC) simulations indicate that relativistic reconnection is a viable mechanism to circumvent the above difficulties. The reconnection electric field localized at X-points linearly accelerates particles with little radiative energy losses. In this paper, we check whether this mechanism survives in 3D, using a set of large PIC simulations with radiation reaction force and with a guide field. In agreement with earlier works, we find that the relativistic drift kink instability deforms and then disrupts the layer, resulting in significant plasma heating but few non-thermal particles. A moderate guide field stabilizes the layer and enables particle acceleration. We report that 3D magnetic reconnection can accelerate particles above the standard radiation reaction limit, although the effect is less pronounced than in 2D with no guide field. We confirm that the highest energy particles form compact bunches within magnetic flux ropes, and a beam tightly confined within the reconnection layer, which could result in the observed Crab flares when, by chance, the beam crosses our line of sight.

The Crab Nebula was formed after the collapse of a massive star about a thousand years ago, leaving behind a pulsar that inflates a bubble of ultra-relativistic electron-positron pairs permeated with magnetic field. The observation of brief but bright flares of energetic gamma rays suggests that pairs are accelerated to PeV energies within a few days; such rapid acceleration cannot be driven by shocks. Here, it is argued that the flares may be the smoking gun of magnetic dissipation in the Nebula. Using 2D and 3D particle-in-cell simulations, it is shown that the observations are consistent with relativistic magnetic reconnection, where pairs are subject to strong radiative cooling. The Crab flares may highlight the importance of relativistic magnetic reconnection in astrophysical sources.

After the destruction of the star during a tidal disruption event (TDE), the cataclysmic encounter between a star and the supermassive black hole (SMBH) of a galaxy, approximately half of the original stellar debris falls back onto the hole at a rate that can initially exceed the Eddington limit by orders of magnitude. We argue that the angular momentum of this matter is too low to allow it to attain a disk-like configuration with accretion proceeding at a mildly super-Eddington rate, the excess energy being carried away by a combination of radiative losses and radially distributed winds. Instead, we propose that the infalling gas traps accretion energy until it inflates into a weakly-bound, quasi-spherical structure with gas extending nearly to the poles. We study the structure and evolution of such "Zero-Bernoulli accretion" flows (ZEBRAs) as a model for the super-Eddington phase of TDEs. We argue that such flows cannot stop extremely super-Eddington accretion from occurring, and that once the envelope is maximally inflated, any excess accretion energy escapes through the poles in the form of powerful jets. We compare the predictions of our model to Swift J1644+57, the putative super-Eddington TDE, and show that it can qualitatively reproduce some of its observed features. Similar models, including self-gravity, could be applicable to gamma-ray bursts from collapsars and the growth of supermassive black hole seeds inside quasi-stars.

We investigate the contribution of the beamed jet component to the high energy emission in young and compact extragalactic radio sources, focusing for the first time on the gamma-ray band. We derive predictions on the gamma-ray luminosities associated with the relativistic jet assuming a leptonic radiative model. The high energy emission is produced via Compton scattering by the relativistic electrons in a spherical region at the considered scales ($\lesssim$10 kpc). Simulations show a wide range of gamma-ray luminosities, with intensities up to $\sim10^{46}-10^{48}$ erg s$^{-1}$ depending on the assumed jet parameters. We find a highly linear relation between the simulated X-ray and gamma-ray luminosities that can be used to select candidates for a gamma-ray detection. We compare the simulated luminosity distributions in the radio, X-ray and gamma-ray regimes with observations for the largest sample of X-ray detected young radio quasars. Our analysis of $\sim$4-year Fermi Large Area Telescope (LAT) data does not give any statistically significant detection. However, the majority of the model-predicted gamma-ray fluxes for the sample are near or below the current Fermi-LAT flux threshold and compatible with the derived upper limits. Our study gives constraints on the minimum jet power ($L_{jet,kin}/L_{disk}>0.01$), on a potential jet contribution to the X-ray emission in the most compact sources ($\lesssim1$ kpc) and on the particles to magnetic field energy density ratio in broad agreement with equipartition assumption.

We perform a suite of two- and three-dimensional magnetohydrodynamic (MHD) simulations with the Athena code of the non-driven Kelvin-Helmholtz instability in the subsonic, weak magnetic field limit. Focusing the analysis on the non-linear turbulent regime, we quantify energy transfer on a scale-by-scale basis and identify the physical mechanisms responsible for energy exchange by developing the diagnostic known as spectral energy transfer function analysis. At late times when the fluid is in a state of MHD turbulence, magnetic tension mediates the dominant mode of energy injection into the magnetic reservoir, whereby turbulent fluid motions twist and stretch the magnetic field lines. This generated magnetic energy turbulently cascades to smaller scales, while being exchanged backwards and forwards with the kinetic energy reservoir, until finally being dissipated. Incorporating explicit dissipation pushes the dissipation scale to larger scales than if the dissipation were entirely numerical. For scales larger than the dissipation scale, we show that the physics of energy transfer in decaying MHD turbulence is robust to numerical effects.

We use numerical simulations to explore whether direct collapse can lead to the formation of SMBH seeds at high-z. We follow the evolution of gas within DM halos of 2 x 10^8 Mo and 1 kpc. We adopt cosmological density profiles and j-distributions. Our goal is to understand how the collapsing flow overcomes the centrifugal barrier and whether it is subject to fragmentation. We find that the collapse leads either to a central runaway or to off-center fragmentation. A disk-like configuration is formed inside the centrifugal barrier. For more cuspy DM distribution, the gas collapses more and experiences a bar-like perturbation and a central runaway. We have followed this inflow down to ~10^{-4} pc. The flow remains isothermal and the specific angular momentum is efficiently transferred by gravitational torques in a cascade of nested bars. This cascade supports a self-similar, disk-like collapse. In the collapsing phase, virial supersonic turbulence develops and fragmentation is damped. For larger initial DM cores the timescales become longer. In models with more organized initial rotation, a torus forms and appears to be supported by turbulent motions. The evolution depends on the competition between two timescales, corresponding to the onset of the central runaway and off-center fragmentation. For less organized rotation, the torus is greatly weakened, the central accretion timescale is shortened, and off-center fragmentation is suppressed --- triggering the central runaway even in previously `stable' models. The resulting SMBH masses lie in the range 2 x 10^4 Mo - 2 x 10^6 Mo, much higher than for Population III remnants. We argue that the above upper limit appears to be more realistic mass. Corollaries of this model include a possible correlation between SMBH and DM halo masses, and similarity between the SMBH and halo mass functions, at time of formation.

Since the discovery of kiloparsec-scale X-ray emission from quasar jets, the physical processes responsible for their high-energy emission have been poorly defined. A number of mechanisms are under active debate, including synchrotron radiation, inverse-Comptonized CMB (IC/CMB) emission, and other Comptonization processes. In a number of cases, the optical and X-ray emission of jet regions are inked by a single spectral component, and in those, high- resolution multi-band imaging and polarimetry can be combined to yield a powerful diagnostic of jet emission processes. Here we report on deep imaging photometry of the jet of PKS 1136$-$135 obtained with the {\it Hubble Space Telescope.} We find that several knots are highly polarized in the optical, with fractional polarization $\Pi>30%$. When combined with the broadband spectral shape observed in these regions, this is very difficult to explain via IC/CMB models, unless the scattering particles are at the lowest-energy tip of the electron energy distribution, with Lorentz factor $\gamma \sim 1$, and the jet is also very highly beamed ($\delta \geq 20$) and viewed within a few degrees of the line of sight. We discuss both the IC/CMB and synchrotron interpretation of the X-ray emission in the light of this new evidence, presenting new models of the spectral energy distribution and also the matter content of this jet. The high polarizations do not completely rule out the possibility of IC/CMB optical-to-X-ray emission in this jet, but they do strongly disfavor the model. We discuss the implications of this finding, and also the prospects for future work.

It is generally accepted that astrophysical sources cannot emit synchrotron radiation above 160 MeV in their rest frame. This limit is given by the balance between the accelerating electric force and the radiation reaction force acting on the electrons. The discovery of synchrotron gamma-ray flares in the Crab Nebula, well above this limit, challenges this classical picture of particle acceleration. To overcome this limit, particles must accelerate in a region of high electric field and low magnetic field. This is possible only with a non-ideal magnetohydrodynamic process, like magnetic reconnection. We present the first numerical evidence of particle acceleration beyond the synchrotron burnoff limit, using a set of 2D particle-in-cell simulations of ultra-relativistic pair plasma reconnection. We use a new code, Zeltron, that includes self-consistently the radiation reaction force in the equation of motion of the particles. We demonstrate that the most energetic particles move back and forth across the reconnection layer, following relativistic Speiser orbits. These particles then radiate >160 MeV synchrotron radiation rapidly, within a fraction of a full gyration, after they exit the layer. Our analysis shows that the high-energy synchrotron flux is highly variable in time because of the strong anisotropy and inhomogeneity of the energetic particles. We discover a robust positive correlation between the flux and the cut-off energy of the emitted radiation, mimicking the effect of relativistic Doppler amplification. A strong guide field quenches the emission of >160 MeV synchrotron radiation. Our results are consistent with the observed properties of the Crab flares, supporting the reconnection scenario.

Constraining the accretion flow geometry of black hole binaries in outburst is complicated by the inability of simplified multi-colour disc models to distinguish between changes in the inner disc radius and alterations to the emergent spectrum, parameterised by the phenomenological colour correction factor, f_col. We analyse Rossi X-ray Timing Explorer observations of the low mass Galactic black hole X-ray binary, GX 339-4, taken over seven epochs when the source was experiencing a state transition. The accretion disc component is isolated using a pipeline resulting in robust detections for disc luminosities, 0.001 < L_disc / L_Edd < 0.5. Assuming that the inner disc remains situated at the innermost stable circular orbit over the course of a state transition, we measure the relative degree of change in f_col required to explain the spectral evolution of the disc component. A variable f_col that increases by a factor of ~ 2.0 - 3.5 as the source transitions from the high/soft state to the low/hard state can adequately explain the observed disc spectral evolution. For the observations dominated by a disc component, the familiar scaling between the disc luminosity and effective temperature, L_disc ~ T_eff^4, is observed; however, significant deviations from this relation appear when GX 339-4 is in the hard intermediate and low/hard states. Allowing for an evolving f_col between spectral states, the L_disc-T_eff^4 law is recovered over the full range of disc luminosities, although this depends heavily on the physically conceivable range of f_col. We demonstrate that physically reasonable changes in f_col provide a viable description for multiple state transitions of a black hole binary without invoking radial motion of the inner accretion disc.

We present the results of special relativistic, adaptive mesh refinement, 3D simulations of gamma-ray burst jets expanding inside a realistic stellar progenitor. Our simulations confirm that relativistic jets can propagate and break out of the progenitor star while remaining relativistic. This result is independent of the resolution, even though the amount of turbulence and variability observed in the simulations is greater at higher resolutions. We find that the propagation of the jet head inside the progenitor star is slightly faster in 3D simulations compared to 2D ones at the same resolution. This behavior seems to be due to the fact that the jet head in 3D simulations can wobble around the jet axis, finding the spot of least resistance to proceed. Most of the jet properties, such as density, pressure, and Lorentz factor, are only marginally affected by the dimensionality of the simulations and therefore results from 2D simulations can be considered reliable.

We argue that the magnetic flux threading the black hole, rather than black hole spin or Eddington ratio, is the dominant factor in launching powerful jets and thus determining the radio loudness of active galactic nuclei (AGN). Most AGN are radio quiet because the thin accretion disks that feed them are inefficient in depositing magnetic flux close to the black hole. Flux accumulation is more likely to occur during a hot accretion (or thick disk) phase, and we argue that radio-loud quasars and strong emission-line radio galaxies occur only when a massive, cold accretion event follows an episode of hot accretion. Such an event might be triggered by the merger of a giant elliptical galaxy with a disk galaxy. This picture supports the idea that flux accumulation can lead to the formation of a so-called magnetically choked accretion flow (MCAF). The large observed range in radio loudness reflects not only the magnitude of the flux pressed against the black hole, but also the decrease in UV flux from the disk, due to its disruption by the "magnetosphere" associated with the accumulated flux. While the strongest jets result from the secular accumulation of flux, moderate jet activity can also be triggered by fluctuations in the magnetic flux deposited by turbulent, hot inner regions of otherwise thin accretion disks, or by the dissipation of turbulent fields in accretion disk coronae. These processes could be responsible for jet production in Seyferts and low-luminosity AGN, as well as jets associated with X-ray binaries.

We present the results of a set of numerical simulations of long-duration gamma-ray burst jets associated with massive, compact stellar progenitors. The simulations extend to large radii and allow us to locate the region in which the peak frequency of the advected radiation is set before the radiation is released at the photosphere. Light curves and spectra are calculated for different viewing angles as well as different progenitor structures and jet properties. We find that the radiation released at the photosphere of matter-dominated jets is able to reproduce the observed Amati and energy-Lorentz factor correlations. Our simulations also predict a correlation between the burst energy and the radiative efficiency of the prompt phase, consistent with observations.

We study the collimation of a highly magnetized jet by a surrounding cocoon that forms as a result of the interaction of the jet with the external medium. We show that in regions where the jet is well confined by the cocoon, current-driven instabilities should develop over timescales shorter than the expansion time of the jet's head. We speculate that these instabilities would give rise to complete magnetic field destruction, whereby the jet undergoes a transition from high-to-low sigma above the collimation zone. Using this assumption, we construct a self-consistent model for the evolution of the jet-cocoon system in an ambient medium of arbitrary density profile. We apply the model to jet breakout in long GRBs, and show that the jet is highly collimated inside the envelope of the progenitor star, and is likely to remain confined well after breakout. We speculate that this strong confinement may provide a channel for magnetic field conversion in GRB outflows, whereby the hot, low-sigma jet section thereby produced is the source of the photospheric emission observed in many bursts.

We study the collimation of relativistic magnetohydrodynamic jets by the pressure of an ambient medium, in the limit where the jet interior loses causal contact with its surroundings. This follows up a hydrodynamic study in a previous paper, adding the effects of a toroidal magnetic field threading the jet. As the ultrarelativistic jet encounters an ambient medium with a pressure profile with a radial scaling of p ~ r^-eta where 2<eta<4, it loses causal contact with its surroundings and forms a boundary layer with a large pressure gradient. By constructing self-similar solutions to the fluid equations within this boundary layer, we examine the structure of this layer as a function of the external pressure profile. We show that the boundary layer always becomes magnetically dominated far from the source, and that in the magnetic limit, physical self-similar solutions are admitted in which the total pressure within the layer decreases linearly with distance from the contact discontinuity inward. These solutions suggest a `hollow cone' behavior of the jet, with the boundary layer thickness prescribed by the value of eta. In contrast to the hydrodynamical case, however, the boundary layer contains an asymptotically vanishing fraction of the jet energy flux.

We perform a local (short-wavelength) linear stability analysis of an axisymmetric column of magnetized plasma with a nearly toroidal magnetic field and a smooth poloidal velocity shear by perturbing the equations of relativistic magnetohydrodynamics. We identify two types of unstable modes, which we call 'exponential' and 'overstable', respectively. The exponential modes are present in the static equilibria and their growth rates decrease with increasing velocity shear. The overstable modes are driven by the effects of velocity shear and dominate the exponential modes for sufficiently high shear values. We argue that these local instabilities can provide an important energy dissipation mechanism in astrophysical relativistic jets. Strong co-moving velocity shear arises naturally in the magnetic acceleration mechanism, therefore it may play a crucial role in converting Poynting-flux-dominated jets into matter-dominated jets, regulating the global acceleration and collimation processes, and producing the observed emission of blazars and gamma-ray bursts.

We study theoretical implications of a rapid Very-High-Energy (VHE) flare detected by MAGIC in the Flat-Spectrum Radio Quasar PKS 1222+216. The minimum distance from the jet origin at which this flare could be produced is 0.5 pc. A moderate Doppler factor of the VHE source, D_{VHE} ~ 20, is allowed by all opacity constraints. The concurrent High-Energy (HE) emission observed by Fermi provides estimates of the total jet power and the jet magnetic field strength. Energetic constraints for the VHE flare are extremely tight: for an isotropic particle distribution they require a huge co-moving energy density in the emitting region and a very efficient radiative process. We disfavor hadronic processes due to their low radiative efficiency, as well as the synchrotron scenario recently proposed for the case of HE flares in the Crab Nebula, since the parameters needed to overcome the radiative losses are quite extreme. The VHE emission can be explained by the Synchrotron Self-Compton (SSC) mechanism for D_{VHE} ~ 20 or by the External Radiation Compton (ERC) mechanism involving the infrared radiation of the dusty torus for D_{VHE} ~ 50. After discussing several alternative scenarios, we propose that the extreme energy density constraint can be satisfied when the emission comes from highly anisotropic short-lived bunches of particles formed by the kinetic beaming mechanism in magnetic reconnection sites. By focusing the emitting particles into very narrow beams, this mechanism allows one to relax the causality constraint on the source size, decreasing the required energy density by 4 orders of magnitude.

We report on the first study of the angular distribution of energetic particles and radiation generated in relativistic collisionless electron-positron pair plasma reconnection, using two-dimensional particle-in-cell simulations. We discover a strong anisotropy of the particles accelerated by reconnection and the associated strong beaming of their radiation. The focusing of particles and radiation increases with their energy; in this sense, this "kinetic beaming" effect differs fundamentally from the relativistic Doppler beaming usually invoked in high-energy astrophysics, in which all photons are focused and boosted achromatically. We also present, for the first time, the modeling of the synchrotron emission as seen by an external observer during the reconnection process. The expected lightcurves comprise several bright symmetric sub-flares emitted by the energetic beam of particles sweeping across the line of sight intermittently, and exhibit super-fast time variability as short as about one tenth of the system light-crossing time. The concentration of the energetic particles into compact regions inside magnetic islands and particle anisotropy explain the rapid variability. This radiative signature of reconnection can account for the brightness and variability of the gamma-ray flares in the Crab Nebula and in blazars.

We present a set of numerical simulations of stellar explosions induced by relativistic jets emanating from a central engine sitting at the center of compact, dying stars. We explore a wide range of durations of the central engine activity, two candidate stellar progenitors, and two possible values of the total energy release. We find that even if the jets are narrowly collimated, their interaction with the star unbinds the stellar material, producing a stellar explosion. We also find that the outcome of the explosion can be very different depending on the duration of the engine activity. Only the longest-lasting engines result in successful gamma-ray bursts. Engines that power jets only for a short time result in relativistic supernova explosions, akin to observed engine-driven SNe such as SN2009bb. Engines with intermediate durations produce weak gamma-ray bursts, with properties similar to nearby bursts such as GRB 980425. Finally, we find that the engines with the shortest durations, if they exist in nature, produce stellar explosions that lack sizable amounts of relativistic ejecta and are therefore dynamically indistinguishable from ordinary core-collapse supernovae.

We propose that the growth of supermassive black holes is associated mainly with brief episodes of highly super-Eddington infall of gas ("hyperaccretion"). This gas is not swallowed in real time, but forms an envelope of matter around the black hole that can be swallowed gradually, over a much longer timescale. However, only a small fraction of the black hole mass can be stored in the envelope at any one time. We argue that any infalling matter above a few per cent of the hole's mass is ejected as a result of the plunge in opacity at temperatures below a few thousand degrees K, corresponding to the Hayashi track. The speed of ejection of this matter, compared to the velocity dispersion (sigma) of the host galaxy's core, determines whether the ejected matter is lost forever or returns eventually to rejoin the envelope, from which it can be ultimately accreted. The threshold between matter recycling and permanent loss defines a relationship between the maximum black hole mass and sigma that resembles the empirical M_BH-sigma relation.

We study the collimation of relativistic hydrodynamic jets by the pressure of an ambient medium in the limit where the jet interior has lost causal contact with its surroundings. For a jet with an ultrarelativistic equation of state and external pressure that decreases as a power of spherical radius, p \propto r^(-eta), the jet interior will lose causal contact when eta > 2. However, the outer layers of the jet gradually collimate toward the jet axis as long as eta < 4, leading to the formation of a shocked boundary layer. Assuming that pressure-matching across the shock front determines the shape of the shock, we study the resulting structure of the jet in two ways: first by assuming that the pressure remains constant across the shocked boundary layer and looking for solutions to the shock jump equations, and then by constructing self-similar boundary-layer solutions that allow for a pressure gradient across the shocked layer. We demonstrate that the constant-pressure solutions can be characterized by four initial parameters that determine the jet shape and whether the shock closes to the axis. We show that self-similar solutions for the boundary layer can be constructed that exhibit a monotonic decrease in pressure across the boundary layer from the contact discontinuity to the shock front, and that the addition of this pressure gradient in our initial model generally causes the shock front to move outwards, creating a thinner boundary layer and decreasing the tendency of the shock to close. We discuss trends based on the value of the pressure power-law index eta.

We present the results of a numerical investigation of current-driven instability in magnetized jets. Utilizing the well-tested, relativistic magnetohydrodynamic code Athena, we construct an ensemble of local, co-moving plasma columns in which initial radial force balance is achieved through various combinations of magnetic, pressure, and rotational forces. We then examine the resulting flow morphologies and energetics to determine the degree to which these systems become disrupted, the amount of kinetic energy amplification attained, and the non-linear saturation behaviors. Our most significant finding is that the details of initial force balance have a pronounced effect on the resulting flow morphology. Models in which the initial magnetic field is force-free deform, but do not become disrupted. Systems that achieve initial equilibrium by balancing pressure gradients and/or rotation against magnetic forces, however, tend to shred, mix, and develop turbulence. In all cases, the linear growth of current-driven instabilities is well-represented by analytic models. CDI-driven kinetic energy amplification is slower and saturates at a lower value in force-free models than in those that feature pressure gradients and/or rotation. In rotating columns, we find that magnetized regions undergoing rotational shear are driven toward equipartition between kinetic and magnetic energies. We show that these results are applicable for a large variety of physical parameters, but we caution that algorithmic decisions (such as choice of Riemann solver) can affect the evolution of these systems more than physically motivated parameters.

High resolution spectra of Galactic Black Holes (GBH) reveal the presence of highly ionised absorbers. In one GBH, accreting close to the Eddington limit for more than a decade, a powerful accretion disc wind is observed to be present in softer X-ray states and it has been suggested that it can carry away enough mass and energy to quench the radio jet. Here we report that these winds, which may have mass outflow rates of the order of the inner accretion rate or higher, are an ubiquitous component of the jet-free soft states of all GBH. We furthermore demonstrate that these winds have an equatorial geometry with opening angles of few tens of degrees, and so are only observed in sources in which the disc is inclined at a large angle to the line of sight. The decrease in Fe XXV / Fe XXVI line ratio with Compton temperature, observed in the soft state, suggests a link between higher wind ionisation and harder spectral shapes. Although the physical interaction between the wind, accretion flow and jet is still not fully understood, the mass flux and power of these winds, and their presence ubiquitously during the soft X-ray states suggests they are fundamental components of the accretion phenomenon.

The gamma-ray space telescopes AGILE and Fermi detected short and bright synchrotron gamma-ray flares at photon energies above 100 MeV in the Crab Nebula. This discovery suggests that electron-positron pairs in the nebula are accelerated to PeV energies in a milliGauss magnetic field, which is difficult to explain with classical models of particle acceleration and pulsar wind nebulae. We investigate whether particle acceleration in a magnetic reconnection layer can account for the puzzling properties of the flares. We numerically integrate relativistic test-particle orbits in the vicinity of the layer, including the radiation reaction force, and using analytical expressions for the large-scale electromagnetic fields. As they get accelerated by the reconnection electric field, the particles are focused deep inside the current layer where the magnetic field is small. The electrons suffer less from synchrotron losses and are accelerated to extremely high energies. Population studies show that, at the end of the layer, the particle distribution piles up at the maximum energy given by the electric potential drop and is focused into a thin fan beam. Applying this model to the Crab Nebula, we find that the emerging synchrotron emission spectrum peaks above 100 MeV and is close to the spectral shape of a single electron. The flare inverse Compton emission is negligible and no detectable emission is expected at other wavelengths. This mechanism provides a plausible explanation for the gamma-ray flares in the Crab Nebula and could be at work in other astrophysical objects such as relativistic jets in active galactic nuclei.

We reformulate the adiabatic inflow-outflow (ADIOS) model for radiatively inefficient accretion flows, treating the inflow and outflow zones on an equal footing. For purely adiabatic flows (i.e., with no radiative losses), we show that the mass flux in each zone must satisfy Mdot ~ R^n with n=1, in contrast to previous work in which 0<n< 1 is a free parameter but in rough agreement with numerical simulations. We also demonstrate that the resulting two-zone ADIOS models are not dynamically self-consistent without the introduction of an energy source close in to the central regions of the flow; we identify this with the energy liberated by accretion. We explore the parameter space of non-radiative models and show that both powerful winds and gentle breezes are possible. When small radiative losses (with fixed efficiency) are included, any centrally injected energy flux is radiated away and the system reverts to a power-law behavior with n < 1, where n falls in a small range determined by the fractional level of radiative losses. We also present an ADIOS model for hypercritical (super-Eddington) disk accretion, in which the radiative losses are closely related to the flow geometry. We suggest that hyperaccretion can lead to either winds or breezes.

The recent discovery of day-long gamma-ray flares in the Crab Nebula, presumed to be synchrotron emission by PeV (10^{15} eV) electrons in milligauss magnetic fields, presents a strong challenge to particle acceleration models. The observed photon energies exceed the upper limit (~100 MeV) obtained by balancing the acceleration rate and synchrotron radiation losses under standard conditions where the electric field is smaller than the magnetic field. We argue that a linear electric accelerator, operating at magnetic reconnection sites, is able to circumvent this difficulty. Sufficiently energetic electrons have gyroradii so large that their motion is insensitive to small-scale turbulent structures in the reconnection layer and is controlled only by large-scale fields. We show that such particles are guided into the reconnection layer by the reversing magnetic field as they are accelerated by the reconnection electric field. As these electrons become confined within the current sheet, they experience a decreasing perpendicular magnetic field that may drop below the accelerating electric field. This enables them to reach higher energies before suffering radiation losses and hence to emit synchrotron radiation in excess of the 100 MeV limit, providing a natural resolution to the Crab gamma-ray flare paradox.

We present the results of a numerical investigation of the spectra and light curves of the emission from the photospheres of long-duration gamma-ray burst jets. We confirm that the photospheric emission has high efficiency and we show that the efficiency increases slightly with the off-axis angle. We show that the peak frequency of the observed spectrum is proportional to the square root of the photosphere's luminosity, in agreement with the Amati relation. However, a quantitative comparison reveals that the thermal peak frequency is too small for the corresponding total luminosity. As a consequence, the radiation must be out of thermal equilibrium with the baryons in order to reproduce the observations. Finally, we show that the spectrum integrated over the emitting surface is virtually indistinguishable from a Planck law, and therefore an additional mechanism has to be identified to explain the non-thermal behavior of the observed spectra at both high and low frequencies.

We discuss strong gravitational lensing of gravitational waves from merging of massive black hole binaries in the context of the LISA mission. Detection of multiple events would provide invaluable information on competing theories of gravity, evolution and formation of structures and, with complementary observations, constraints on H_0 and other cosmological parameters. Most of the optical depth for lensing is provided by intervening massive galactic halos, for which wave optics effects are negligible. Probabilities to observe multiple events are sizable for a broad range of formation histories. For the most optimistic models, up to 4 multiple events with a signal to noise ratio >= 8 are expected in a 5-year mission. Chances are significant even for conservative models with either light (<= 60%) or heavy (<= 40%) seeds. Due to lensing amplification, some intrinsically too faint signals are brought over threshold (<= 2 per year).

We present a numerical simulation of a gamma-ray burst jet from a long-lasting engine in the core of a 16 solar mass Wolf-Rayet star. The engine is kept active for 6000 s with a luminosity that decays in time as a power-law with index -5/3. Even though there is no short time-scale variability in the injected engine luminosity, we find that the jet's kinetic luminosity outside the progenitor star is characterized by fluctuations with relatively short time scale. We analyze the temporal characteristics of those fluctuations and we find that they are consistent with the properties of observed flares in X-ray afterglows. The peak to continuum flux ratio of the flares in the simulation is consistent with some, but not all, the observed flares. We propose that propagation instabilities, rather than variability in the engine luminosity, are responsible for the X-ray flares with moderate contrast. Strong flares such as the one detected in GRB 050502B, instead, cannot be reproduced by this model and require strong variability in the engine activity.

We present the results of hydrodynamical simulations of gamma-ray burst jets propagating through their stellar progenitor material and subsequently through the surrounding circumstellar medium. We consider both jets that are injected with constant properties in the center of the star and jets injected with a variable luminosity. We show that the variability properties of the jet outside the star are a combination of the variability injected by the engine and the variability caused by the jet propagation through the star. Comparing power spectra for the two cases shows that the variability injected by the engine is preserved even if the jet is heavily shocked inside the star. Such shocking produces additional variability at long time scales, of order several seconds. Our findings suggest that the broad pulses of several seconds duration typically observed in gamma-ray bursts are due to the interaction of the jet with the progenitor, while the short-timescale variability, characterized by fluctuations on time scales of milliseconds, has to be injected at the base of the jet. Studying the properties of the fast variability in GRBs may therefore provide clues to the nature of the inner engine and the mechanisms of energy extraction from it.

We study the spectrum of high frequency radiation emerging from mildly dissipative photospheres of long-duration gamma-ray burst outflows. Building on the results of recent numerical investigations, we assume that electrons are heated impulsively to mildly relativistic energies by either shocks or magnetic dissipation at Thomson optical depths of several and subsequently cool by inverse Compton, scattering off the thermal photons of the photosphere. We show that even in the absence of magnetic field and non-thermal leptons, inverse Compton scattering produces power-law tails that extend from the peak of the thermal radiation, at several hundred keV, to several tens of MeV, and possibly up to GeV energies. The slope of the high-frequency power-law is predicted to vary substantially during a single burst, and the model can easily account for the diversity of high-frequency spectra observed by BATSE. Our model works in baryonic as well as in magnetically dominated outflows, as long as the magnetic field component is not overwhelmingly dominant.

We construct a numerical model of emission from minijets, localized flows driven by magnetic reconnection inside Poynting-flux-dominated jets proposed to explain the ultrafast variability of blazars. The geometrical structure of the model consists of two wedge-like regions of relativistically flowing gas, separated by a stationary shock. The dynamics is based on solutions of relativistic magnetic reconnection with a guide field from Lyubarsky (2005). Electron distributions in each region are chosen to the match the pressure and density of the local plasma. Synchrotron emission from both regions is used to calculate Compton scattering, Compton drag and photon-photon opacity effects, with exact treatment of anisotropy and the Klein-Nishina regime. Radiative effects on plasma are taken into account, including the dependence of pressure on electron radiative losses and adiabatic heating of the flow decelerating under Compton drag. The results are applied to the July 2006 flare in the BL Lac object PKS 2155-304, with the aim of matching TeV flux measurements by H.E.S.S. with models that satisfy the variability constraints, while keeping X-ray emission below simultaneous Chandra observations. We find that models of isolated minijets with a significant guide field overproduce X-ray emission, and that we must take into account the radiative interaction of oppositely-oriented minijets in order to achieve a high enough dominance by Comptonized TeV radiation. We argue that such interactions are likely to occur in a jet where there is substantial internal reconnection, producing a large number of misaligned minijets. Finally, we show that large jet magnetizations are indeed required to satisfy all observational constraints and that the effective Lorentz factor of the minijet plasma has to be larger than 50, in agreement with earlier one-zone estimates.

We explore the cosmic evolution of massive black hole (MBH) seeds forming within 'quasistars' (QSs), accreting black holes embedded within massive hydrostatic gaseous envelopes. These structures could form if the infall of gas into the center of a halo exceeds about 1 solar mass per year. We use a merger-tree approach to estimate the rate at which QSs might form as a function of redshift, and the statistical properties of the resulting QS and seed black hole populations. We relate the triggering of runaway infall to major mergers of gas-rich galaxies, and to a threshold for global gravitational instability, which we link to the angular momentum of the host. This is the main parameter of our models. Once infall is triggered, its rate is determined by the halo potential; the properties of the resulting QS and seed black hole depend on this rate. After the epoch of QSs, we model the growth of MBHs within their hosts in a merger-driven accretion scenario. We compare MBH seeds grown inside quasistars to a seed model that derives from the remnants of the first metal-free stars, and also study the case in which both channels of MBH formation operate simultaneously. We find that a limited range of QS/MBH formation efficiencies exists that allows one to reproduce observational constraints. Our models match the density of z = 6 quasars, the cumulative mass density accreted onto MBHs (according to Soltan's argument), and the current mass density of MBHs. The mass function of QSs peaks at mass ~ 1e6 solar masses, and we calculate the number counts for the JWST in the 2-10 micron band. We find that JWST could detect up to several QSs per field at z ~ 5 - 10.

We present 2D high-resolution hydrodynamic simulations of the relativistic outflows of long-duration gamma-ray burst progenitors. We analyze the properties of the outflows at wide off-axis angles, produced by the expansion of the hot cocoon that surrounds the jet inside the progenitor star. We find that the cocoon emission at wide angles may have properties similar to those of the subclass of short-duration gamma-ray bursts with persistent X-ray emission. We compute the predicted duration distribution, redshift distribution, and afterglow brightness and we find that they are all in agreement with the observed properties of short GRBs with persistent emission. We suggest that a SN component, the properties of the host galaxies, and late afterglow observations can be used as a crucial test to verify this model.

In a dynamical-radiative model we recently developed to describe the physics of compact, GHz-Peaked-Spectrum (GPS) sources, the relativistic jets propagate across the inner, kpc-sized region of the host galaxy, while the electron population of the expanding lobes evolves and emits synchrotron and inverse-Compton (IC) radiation. Interstellar-medium gas clouds engulfed by the expanding lobes, and photoionized by the active nucleus, are responsible for the radio spectral turnover through free-free absorption (FFA) of the synchrotron photons. The model provides a description of the evolution of the spectral energy distribution (SED) of GPS sources with their expansion, predicting significant and complex high-energy emission, from the X-ray to the gamma-ray frequency domain. Here, we test this model with the broad-band SEDs of a sample of eleven X-ray emitting GPS galaxies with Compact-Symmetric-Object (CSO) morphology, and show that: (i) the shape of the radio continuum at frequencies lower than the spectral turnover is indeed well accounted for by the FFA mechanism; (ii) the observed X-ray spectra can be interpreted as non-thermal radiation produced via IC scattering of the local radiation fields off the lobe particles, providing a viable alternative to the thermal, accretion-disk dominated scenario. We also show that the relation between the hydrogen column densities derived from the X-ray (N_H) and radio (N_HI) data of the sources is suggestive of a positive correlation, which, if confirmed by future observations, would provide further support to our scenario of high-energy emitting lobes.

The jet of the radio galaxy M87 is misaligned, resulting in a Doppler factor delta~1 for emission of plasma moving parallel to the jet. This makes the observed fast TeV flares on timescales of t_v~5R_g/c harder to understand as emission from the jet. In previous work, we have proposed a jets-in-a-jet model for the ultra-fast TeV flares with t_v<<R_g/c seen in Mrk 501 and PKS 2155-304. Here, we show that about half of the minijets beam their emission outside the jet cone. Minijets emitting off the jet axis result in rapidly evolving TeV (and maybe lower energy) flares that can be observed in nearby radio galaxies. The TeV flaring from M87 fits well into this picture, if M87 is a misaligned blazar.

Supermassive stars, with masses greater than a million solar masses, are possible progenitors of supermassive black holes in galactic nuclei. Because of their short nuclear burning timescales, such objects can be formed only when matter is able to accumulate at a rate exceeding ~ 1 solar mass/yr. Here we revisit the structure and evolution of rotationally-stabilized supermassive stars, taking into account their continuous accumulation of mass and their thermal relaxation. We show that the outer layers of supermassive stars are not thermally relaxed during much of the star's main sequence lifetime. As a result, they do not resemble n=3 polytropes, as assumed in previous literature, but rather consist of convective (polytropic) cores surrounded by convectively stable envelopes that contain most of the mass. We compute the structures of these envelopes, in which the specific entropy is proportional to the enclosed mass M(R) to the 2/3-power. By matching the envelope solutions to convective cores, we calculate the core mass as a function of time. We estimate the initial black hole masses formed as a result of core-collapse, and their subsequent growth via accretion from the bloated envelopes ("quasistars") that result. The seed black holes formed in this way could have typical masses in the range ~ 10^4-10^5 solar masses, considerably larger than the remnants thought to be left by the demise of Population III stars. Supermassive black holes therefore could have been seeded during an epoch of rapid infall considerably later than the era of Pop III star formation.

The emission mechanisms in extragalactic jets include synchrotron and various inverse-Compton processes. At low (radio through infrared) energies, it is widely agreed that synchrotron emission dominates in both low-power (FR I) and high-power (FR II and quasar) jets, because of the power-law nature of the spectra observed and high polarizations. However, at higher energies, the emission mechanism for high-power jets at kpc scales is hotly debated. Two mechanisms have been proposed: either inverse-Compton of cosmic microwave background photons or synchrotron emission from a second, high-energy population of electrons. Here we discuss optical polarimetry as a method for diagnosing the mechanism for the high-energy emission in quasar jets, as well as revealing the jet's three-dimensional energetic and magnetic field structure. We then discuss high-energy emission mechanisms for powerful jets in the light of the HST polarimetry of PKS 1136-135.

Rapid inflows associated with early galaxy formation lead to the accumulation of self-gravitating gas in the centers of proto-galaxies. Such gas accumulations are prone to non-axisymmetric instabilities, as in the well-known Maclaurin sequence of rotating ellipsoids, which are accompanied by a catastrophic loss of angular momentum (J). Self-gravitating gas is also intuitively associated with star formation. However, recent simulations of the infall process display highly turbulent continuous flows. We propose that J-transfer, which enables the inflow, also suppresses fragmentation. Inefficient J loss by the gas leads to decay of turbulence, triggering global instabilities and renewed turbulence driving. Flow regulated in this way is stable against fragmentation, whilst staying close to the instability threshold for bar formation -- thick self-gravitating disks are prone to global instabilities before they become unstable locally. On smaller scales, the fraction of gravitationally unstable matter swept up by shocks in such a flow is a small and decreasing function of the Mach number. We conclude counterintuitively that gas able to cool down to a small fraction of its virial temperature will not fragment as it collapses. This provides a venue for supermassive black holes to form via direct infall, without the intermediary stage of forming a star cluster. Some black holes could have formed or grown in massive halos at low redshifts. Thus the fragmentation is intimately related to J redistribution within the system: it is less dependent on the molecular and metal cooling but is conditioned by the ability of the flow to develop virial, supersonic turbulence.

We study the polarization properties of the prompt emission of gamma-ray bursts produced by fragmented fireballs. Such fireballs, known in the literature under various names, are made by the superposition of many individual blobs, each of which produces a spike in the light curve. The differences between pulses are due both to the intrinsic diversity in the blobs' properties and to their orientation with respect to the line of sight. We show that the peak flux and the polarization of each pulse are connected through the orientation of the blob that produces the pulse, while the position angle fluctuates randomly from one pulse to the next and is constant within an individual pulse. The most polarized pulses are those with approximately one tenth the peak flux of the brightest pulse. These conclusions do not depend on the assumed radiation mechanism nor on the energy and Lorentz factor of the blobs. We compare the prediction of this model to a simulated set of observations, showing that a limited sample of GRBs with time-resolved polarization measurements would provide a crucial test for this model. We finally show that a hint of the predicted correlation may have been already observed.

We numerically analyze the evolution of a long-duration gamma-ray burst jet as it leaves the progenitor star and propagates to the photospheric radius, where radiation can be released. We find that the interaction of the relativistic material with the progenitor star has influences well beyond the stellar surface. Tangential collimation shocks are observed throughout the jet evolution, out to about 100 stellar radii, which is the whole range of our simulation. We find that the jet is internally hot at the photospheric radius and we compute the photospheric emission. The photosphere is a very efficient radiator, capable of converting more than half of the total energy of the jet into radiation. We show that bright photospheres are a common feature of jets born inside massive progenitor stars and that this effect can explain the high radiative efficiency observed in long-duration bursts.

The fast TeV variability of the blazars Mrk 501 and PKS 2155--304 implies a compact emitting region that moves with a bulk Lorentz factor of Gamma_{em}~100 toward the observer. The Lorentz factor is clearly in excess of the jet Lorentz factors Gamma_j\simless 10 measured on sub-pc scales in these sources. We propose that the TeV emission originates from compact emitting regions that move relativistically {\it within} a jet of bulk Gamma_j~10. This can be physically realized in a Poynting flux-dominated jet. We show that if a large fraction of the luminosity of the jet is prone to magnetic dissipation through reconnection, then material outflowing from the reconnection regions can efficiently power the observed TeV flares through synchrotron-self-Compton emission. The model predicts simultaneous far UV/soft X-ray flares.

The study of the broad-band emission of GHz-Peaked-Spectrum (GPS) radio galaxies is a powerful tool to investigate the physical processes taking place in the central, kpc-sized region of their active hosts, where the jets propagate and the lobes expand, interacting with the surrounding interstellar medium (ISM). We recently developed a new dynamical-radiative model to describe the evolution of the GPS phenomenon (Stawarz et al. 2008): as the relativistic jets propagate through the ISM, gradually engulfing narrow-line emitting gas clouds along their way, the electron population of the expanding lobes evolves, emitting synchrotron light, as well as inverse-Compton radiation via up-scattering of the photon fields from the host galaxy and its active nucleus. The model, which successfully reproduces the key features of the GPS radio sources as a class, provides a description of the evolution of their spectral energy distribution (SED) with the lobes' expansion, predicting significant and complex X-ray to gamma-ray emission. We apply here the model to the broad-band SED's of a sample of known, X-ray emitting GPS galaxies, and show that: (i) the free-free absorption mechanism enables us to reproduce the radio continuum at frequencies below the turnover; (ii) the lobes' non-thermal, inverse-Compton emission can account for the observed X-ray spectra, providing a viable alternative to the thermal, accretion-dominated scenario. We also show that, in our sample, the relationship between the X-ray and radio hydrogen column densitities, N_H and N_HI, is suggestive of a positive correlation, which, if confirmed, would support the scenario of high-energy emitting lobes.

When double neutron star or neutron star-black hole binaries merge, the final remnant may comprise a central solar-mass black hole surrounded by a 0.01-0.1 solar masses torus. The subsequent evolution of this disc may be responsible for short gamma-ray bursts (SGRBs). A comparable amount of mass is ejected into eccentric orbits and will eventually fall back to the merger site after approximately 0.01 seconds. In this Paper, we investigate analytically the fate of the fallback matter, which may provide a luminous signal long after the disc is exhausted. We find that matter in the eccentric tail returns at a super-Eddington rate and is eventually (0.1 sec) unable to cool via neutrino emission and accrete all the way to the black hole. Therefore, contrary to previous claims, our analysis suggests that fallback matter is not an efficient source of late time accretion power and is unlikely to cause the late flaring activity observed in SGRB afterglows. The fallback matter rather forms a radiation-driven wind or a bound atmosphere. In both cases, the emitting plasma is very opaque and photons are released with a degraded energy in the X-ray band. We therefore suggest that compact binary mergers could be followed by an "X-ray renaissance", as late as several days to weeks after the merger. This might be observed by the next generation of X-ray detectors.

[ABRIDGED] We study supermassive black hole binary mergers driven by angular momentum loss to small-scale gas discs. Such binaries form after major galaxy mergers, but their fate is unclear since hardening through stellar scattering becomes very inefficient at sub-parsec distances. Gas discs may dominate binary dynamics on these scales, and promote mergers. Using numerical simulations, we investigate the evolution of the orbits of binaries embedded within geometrically thin gas discs. Our simulations directly resolve angular momentum transport within the disc, which at the radii of interest is likely dominated by disc self-gravity. We show that the binary decays at a rate which is in good agreement with analytical estimates, while the eccentricity grows. Saturation of eccentricity growth is not observed up to values e > 0.35. Accretion onto the black holes is variable, and is roughly modulated by the binary orbital frequency. Scaling our results, we analytically estimate the maximum rate of binary decay that is possible without fragmentation occuring within the surrounding gas disc, and compare that rate to an estimate of the stellar dynamical hardening rate. For binary masses in the range 10^5 < M/Msun < 10^8 we find that decay due to gas discs may dominate for separations below a ~ 0.01--0.1 pc, where the disc is optically thick. The minimum merger time scale is shorter than the Hubble time for M < 10^7 Msun. This implies that gas discs could commonly attend relatively low mass black hole mergers, and that a significant population of binaries might exist at separations of a few hundredths of a pc, where the combined decay rate is slowest. For more massive binaries, we suggest that scattering of stars formed within a fragmenting gas disc could act as a significant additional sink of binary angular momentum.

We compute the average luminosity of X-ray flares as a function of time, for a sample of 10 long-duration gamma-ray burst afterglows. The mean luminosity, averaged over a timescale longer than the duration of the individual flares, declines as a power-law in time with index ~-1.5. We elaborate on the properties of the central engine that can produce such a decline. Assuming that the engine is an accreting compact object, and for a standard conversion factor between accretion rate and jet luminosity, the switch between a neutrino-cooled thin disk and a non-cooled thick disk takes place at the transition from the prompt to the flaring phase. We discuss the implications of this coincidence under different scenarios for the powering of the GRB outflow. We also show that the interaction of the outflow with the envelope of the progenitor star cannot produce flares out of a continuous relativistic flow, and conclude that it is the dynamics of the disk or the jet-launching mechanism that generates an intrinsically unsteady outflow on timescales much longer than the dynamical timescale of the system. This is consistent with the fact that X-ray flares are observed in short-duration GRBs as well as in long-duration ones.

Flares overlaid on the smooth power-law decay of Swift X-ray afterglows are rather common, appearing in roughly half the observed light curves. They are a manifestation of the late time activity of the inner engine, since their temporal evolution is too fast to be linked to activity taking place in the external shock blastwave. In this paper we show that the energy emitted in the form of flares decreases with time as a power-law. We discuss several possibilities in which the flares can be powered and the source of the observed variability. We show that late time accretion from a disk can provide the necessary energy input in both classes of short duration and long duration gamma-ray bursts.

We study the structure and evolution of "quasistars", accreting black holes embedded within massive hydrostatic gaseous envelopes. These configurations may model the early growth of supermassive black hole seeds. The accretion rate onto the black hole adjusts so that the luminosity carried by the convective envelope equals the Eddington limit for the total mass. This greatly exceeds the Eddington limit for the black hole mass alone, leading to rapid growth of the black hole. We use analytic models and numerical stellar structure calculations to study the structure and evolution of quasistars. We derive analytically the scaling of the photospheric temperature with the black hole mass and envelope mass, and show that it decreases with time as the black hole mass increases. Once the photospheric temperature becomes lower than 10,000 K, the photospheric opacity drops precipitously and the photospheric temperature hits a limiting value, analogous to the Hayashi track for red giants and protostars, below which no hydrostatic solution for the convective envelope exists. For metal-free (Population III) opacities this limiting temperature is approximately 4000 K. After a quasistar reaches this limiting temperature, the envelope is rapidly dispersed by radiation pressure. We find that black hole seeds with masses between 1000 and 10000 solar masses could form via this mechanism in less than a few Myr.

We consider axisymmetric relativistic jets with a toroidal magnetic field and an ultrarelativistic equation of state, with the goal of studying the lateral structure of jets whose pressure is matched to the pressure of the medium through which they propagate. We find all self-similar steady-state solutions of the relativistic MHD equations for this setup. One of the solutions is the case of a parabolic jet being accelerated by the pressure gradient as it propagates through a medium with pressure declining as p(z)\propto z^{-2}. As the jet material expands due to internal pressure gradients, it runs into the ambient medium resulting in a pile-up of material along the jet boundary, while the magnetic field acts to produce a magnetic pinch along the axis of the jet. Such jets can be in a lateral pressure equilibrium only if their opening angle \theta_j at distance z is smaller than about 1/\gamma, where \gamma is the characteristic bulk Lorentz-factor at this distance; otherwise, different parts of the jet cannot maintain causal contact. We construct maps of optically thin synchrotron emission from our self-similar models. We suggest that the boundary pile-up may be the reason for the limb-brightening of the sub-parsec jet of M87. We find that if the synchrotron emissivity falls with the distance from the jet axis, the polarization fraction rises toward the edge, as seen in 3C273 and Mkn501. Projection effects and the emissivity pattern of the jet have a strong effect on the observed polarization signal, so the interpretation of the polarization data in terms of the geometry of magnetic fields is rather uncertain. For example, jets with toroidal magnetic fields display the `spine-sheath' polarization angle pattern seen in some BL Lac objects.

Here we discuss evolution and broad-band emission of compact (< kpc) lobes in young radio sources. We propose a simple dynamical description for these objects, consisting of a relativistic jet propagating into a uniform gaseous medium in the central parts of an elliptical host. In the framework of the proposed model, we follow the evolution of ultrarelativistic electrons injected from a terminal hotspot of a jet to expanding lobes, taking into account their adiabatic energy losses as well as radiative cooling. This allows us to discuss the broad-band lobe emission of young radio sources. In particular, we argue that the observed spectral turnover in the radio synchrotron spectra of these objects cannot originate from the synchrotron self-absorption process but is most likely due to free-free absorption effects connected with neutral clouds of interstellar medium engulfed by the expanding lobes and photoionized by active centers. We also find a relatively strong and complex high-energy emission component produced by inverse-Compton up-scattering of various surrounding photon fields by the lobes' electrons. We argue that such high energy radiation is strong enough to account for several observed properties of GHz-peaked-spectrum (GPS) radio galaxies at UV and X-ray frequencies. In addition, this emission is expected to extend up to GeV (or possibly even TeV) photon energies and can thus be probed by several modern gamma-ray instruments. In particular, we suggest that GPS radio galaxies should constitute a relatively numerous class of extragalactic sources detected by GLAST.

We consider the effects of eccentricity on the fragmentation of gravitationally unstable accretion disks, using numerical hydrodynamics. We find that eccentricity does not affect the overall stability of the disk against fragmentation, but significantly alters the manner in which such fragments accrete gas. Variable tidal forces around an eccentric orbit slow the accretion process, and suppress the formation of weakly-bound clumps. The "stellar" mass function resulting from the fragmentation of an eccentric disk is found to have a significantly higher characteristic mass than that from a corresponding circular disk. We discuss our results in terms of the disk(s) of massive stars at ~0.1pc from the Galactic Center, and find that the fragmentation of an eccentric accretion disk, due to gravitational instability, is a viable mechanism for the formation of these systems.

We discuss the implications of rapid (few-minute) variability in the TeV flux of blazars, which has been observed recently with the HESS and MAGIC telescopes. The variability timescales seen in PKS 2155-304 and Mrk 501 are much shorter than inferred light-crossing times at the black hole horizon, suggesting that the variability involves enhanced emission in a small region within an outflowing jet. The enhancement could be triggered by dissipation in part of the black hole's magnetosphere at the base of the outflow, or else by instabilities in the jet itself. By considering the energetics of the observed flares, along with the requirement that TeV photons escape without producing pairs, we deduce that the bulk Lorentz factors in the jets must be >50. The distance of the emission region from the central black hole is less well-constrained. We discuss possible consequences for multi-wavelength observations.

Rapid infall of gas in the nuclei of galaxies could lead to the formation of black holes by direct collapse, without first forming stars. Black holes formed in this way would have initial masses of a few solar masses, but would be embedded in massive envelopes that would allow them to grow at a highly super-Eddington rate. Thus, seed black holes as large as 10^3-10^4 solar masses could form very rapidly. I will sketch the basic physics of the direct collapse process and the properties of the accreting envelopes.

Most cool core clusters of galaxies possess active galactic nuclei (AGN) in their centers. These AGN inflate buoyant bubbles containing non-thermal radio emitting particles. If such bubbles efficiently confine cosmic rays (CR) then this could explain ``radio ghosts'' seen far from cluster centers. We simulate the diffusion of cosmic rays from buoyant bubbles inflated by AGN. Our simulations include the effects of the anisotropic particle diffusion introduced by magnetic fields. Our models are consistent with the X-ray morphology of AGN bubbles, with disruption being suppressed by the magnetic draping effect. We conclude that for such magnetic field topologies, a substantial fraction of cosmic rays can be confined inside the bubbles on buoyant rise timescales even when the parallel diffusivity coefficient is very large. For isotropic diffusion at a comparable level, cosmic rays would leak out of the bubbles too rapidly to be consistent with radio observations. Thus, the long confinement times associated with the magnetic suppression of CR diffusion can explain the presence of radio ghosts. We show that the partial escape of cosmic rays is mostly confined to the wake of the rising bubbles, and speculate that this effect could: (1) account for the excitation of the H$\alpha$ filaments trailing behind the bubbles in the Perseus cluster, (2) inject entropy into the metal enriched material being lifted by the bubbles and, thus, help to displace it permanently from the cluster center and (3) produce observable $\gamma$-rays via the interaction of the diffusing cosmic rays with the thermal intracluster medium (ICM).

We study the long-term evolution of relativistic jets in collapsars and examine the effects of viewing angle on the subsequent gamma ray bursts. We carry out a series of high-resolution simulations of a jet propagating through a stellar envelope in 2D cylindrical coordinates using the FLASH relativistic hydrodynamics module. For the first time, simulations are carried out using an adaptive mesh that allows for a large dynamic range inside the star while still being efficient enough to follow the evolution of the jet long after it breaks out from the star. Our simulations allow us to single out three phases in the jet evolution: a precursor phase in which relativistic material turbulently shed from the head of the jet first emerges from the star, a shocked jet phase where a fully shocked jet of material is emerging, and an unshocked jet phase where the jet consists of a free-streaming, unshocked core surrounded by a thin boundary layer of shocked jet material. The appearance of these phases will be different to observers at different angles. The precursor has a wide opening angle and would be visible far off axis. The shocked phase has a relatively narrow opening angle that is constant in time. During the unshocked jet phase the opening angle increases logarithmically with time. As a consequence, some observers see prolonged dead times of emission even for constant properties of the jet injected in the stellar core. We also present an analytic model that is able to reproduce the overall properties of the jet and its evolution. We finally discuss the observational implications of our results, emphasizing the possible ways to test progenitor models through the effects of jet propagation in the star. In an appendix, we present 1D and 2D tests of the FLASH relativistic hydrodynamics module.

We present the first results from a 500 ksec Chandra ACIS-I observation of M87. At soft energies (0.5-1.0 keV), we detect filamentary structures associated with the eastern and southwestern X-ray and radio arms. Many filaments are spatially resolved with widths of ~300 pc. This filamentary structure is particularly striking in the eastern arm where we suggest the filaments are outer edges of a series of plasma-filled, buoyant bubbles whose ages differ by ~6 x 10^6 years. These X-ray structures may be influenced by magnetic filamentation. At hard energies (3.5-7.5 keV), we detect a nearly circular ring of outer radius 2.8' (13 kpc) which provides an unambiguous signature of a weak shock, driven by an outburst from the SMBH. The density rise in the shock is ~1.3 (Mach number, M~1.2). The observed spectral hardening in the ring corresponds to a temperature rise T_shock / T_0 ~ 1.2, or M~1.2, in agreement with the Mach number derived independently from the gas density. Thus, for the first time, we detect gas temperature and density jumps associated with a classical shock in the atmosphere around a supermassive black hole. We also detect two additional surface brightness edges and pressure enhancements at radii of ~0.6' and ~1'. The ~0.6' feature may be over-pressurized thermal gas surrounding the relativistic plasma in the radio cocoon, the ``piston'', produced by the current episode of AGN activity. The over-pressurized gas is surrounded by a cool gas shell. The ~1' feature may be an additional weak shock from a secondary outburst. In an earlier episode, the ``piston'' was responsible for driving the 2.8' shock.

It is now generally accepted that long gamma-ray bursts are associated with the final evolutionary stages of massive stars. As a consequence, their jets must propagate through the stellar progenitor and break out on their surface, before they can reach the photospheric radius and produce the gamma-ray photons. We investigate the role of the progenitor star in shaping the jet properties. We show that even a jet powered by a steady engine can develop a rich phenomenology at the stellar surface. We present special-relativistic simulations and compare the results to analytic considerations. We show that the jet is complex in the time as well as in the angular domain, so that observers located along different lines of sight detect significantly different bursts.

Simulations and analytic arguments suggest that the turbulence driven by magnetorotational instability (MRI) in accretion discs can amplify the toroidal (azimuthal) component of the magnetic field to a point at which magnetic pressure exceeds the combined gas + radiation pressure in the disc. Arguing from the recent analysis by Pessah and Psaltis, and other MRI results in the literature, we conjecture that the limiting field strength for a thin disc is such that the Alfven speed roughly equals the geometric mean of the Keplerian speed and the gas sound speed. We examine the properties of such magnetically-dominated discs, and show that they resolve a number of outstanding problems in accretion disc theory. The discs would be thicker than standard (Shakura-Sunyaev) discs at the same radius and accretion rate, and would tend to have higher colour temperatures. If they transport angular momentum according to an alpha-prescription, they would be stable against the thermal and viscous instabilities that are found in standard disc models. In discs fuelling active galactic nuclei, magnetic pressure support could also alleviate the restriction on accretion rate imposed by disc self-gravity.

The association of at least some long gamma-ray bursts with type Ic supernova explosions has been established beyond reasonable doubt. Theoretically, the challenge is to explain the presence of a light hyper-relativistic flow propagating through a massive stellar core without losing those properties. We discuss the role of the jet-star interaction in shaping the properties of the outflow emerging on the surface of the star. We show that the nature of the inner engine is hidden from the observer for most of the evolution, well beyond the time of the jet breakout on the stellar surface. The discussion is based on analytical considerations as well as high resolution numerical simulations. Finally, the observational consequences of the scenario are addressed in light of the present capabilities.

We consider the dynamical evolution of a disk of stars orbiting a central black hole. In particular, we focus on the effect of the stellar mass function on the evolution of the disk, using both analytic arguments and numerical simulations. We apply our model to the ring of massive stars at ~0.1pc from the Galactic Center, assuming that the stars formed in a cold, circular disk, and find that our model requires the presence of a significant population of massive (>100Msun) stars in order to explain the the observed eccentricities of 0.2-0.3. Moreover, in order to limit the damping of the heavier stars' eccentricities, we also require fewer low-mass stars than expected from a Salpeter mass function, giving strong evidence for a significantly ``top-heavy'' mass function in the rings of stars seen near to the Galactic Center. We also note that the maximum possible eccentricities attainable from circular initial conditions at ages of <10Myr are around 0.4-0.5, and suggest that any rings of stars found with higher eccentricities were probably not formed from circular disks.

The existence of the radius of marginal stability means that accretion flows around black holes invariably undergo a transition from a MHD turbulent disk-like flow to an inward plunging flow. We argue that the plunging inflow can greatly enhance the trapping of large scale magnetic field on the black hole, and therefore may increase the importance of the Blandford-Znajek (BZ) effect relative to previous estimates that ignore the plunge region. We support this hypothesis by constructing and analyzing a toy-model of the dragging and trapping of a large scale field by a black hole disk, revealing a strong dependence of this effect on the effective magnetic Prandtl number of the MHD turbulent disk. Furthermore, we show that the enhancement of the BZ effect depends on the geometric thickness of the accretion disk. This may be, at least in part, the physical underpinnings of the empirical relation between the inferred geometric thickness of a black hole disk and the presence of a radio jet.

We describe a mechanism by which supermassive black holes can form directly in the nuclei of protogalaxies, without the need for seed black holes left over from early star formation. Self-gravitating gas in dark matter halos can lose angular momentum rapidly via runaway, global dynamical instabilities, the so-called "bars within bars" mechanism. This leads to the rapid buildup of a dense, self-gravitating core supported by gas pressure - surrounded by a radiation pressure-dominated envelope - which gradually contracts and is compressed further by subsequent infall. These conditions lead to such high temperatures in the central region that the gas cools catastrophically by thermal neutrino emission, leading to the formation and rapid growth of a central black hole. We estimate the initial mass and growth rate of the black hole for typical conditions in metal-free halos with T_vir ~ 10^4 K, which are the most likely to be susceptible to runaway infall. The initial black hole should have a mass of <~20 solar masses, but in principle could grow at a super-Eddington rate until it reaches ~ 10^4-10^6 solar masses. Rapid growth may be limited by feedback from the accretion process and/or disruption of the mass supply by star formation or halo mergers. Even if super-Eddington growth stops at \~10^3-10^4 solar masses, this process would give black holes ample time to attain quasar-size masses by a redshift of 6, and could also provide the seeds for all supermassive black holes seen in the present universe.

The periodic precession (162--day) and nodding (6.3--day) motions of the jets in SS433 are driven in the outer regions of the disc, whereas the jets themselves, being relativistic, are launched near the black hole at the disc centre. Given that the nutation period is comparable to the dynamical timescale in the outer regions of the disc, it seems unlikely that these periods can be communicated efficiently to the disc centre. Here we propose that the massive outflow observed in SS433 is launched at large radii in the disc, about 1/10 of the outer disc edge, and that it is this outflow which responds to the oscillations of the outer disc and determines the direction of the jets. The massive outflow is launched at large radius because the mass transfer rate is hyper-Eddington. This implies not only that the total luminosity of SS433 exceeds $\le$ by a considerable factor, but also that the radiative output is collimated along the outflow. We thus suggest that SS433 is an ultraluminous X-ray source (ULX) viewed `from the side'. We also suggest that the obscured {\it INTEGRAL} sources may be SS433--like objects, but with slightly lower mass transfer rates.

As soon as it was realized that long GRBs lie at cosmological distances, attempts have been made to use them as cosmological probes. Besides their use as lighthouses, a task that presents mainly the technological challenge of a rapid deep high resolution follow-up, researchers attempted to find the Holy Grail: a way to create a standard candle from GRB observables. We discuss here the attempts and the discovery of the Ghirlanda correlation, to date the best method to standardize the GRB candle. Together with discussing the promises of this method, we will underline the open issues, the required calibrations and how to understand them and keep them under control. Even though GRB cosmology is a field in its infancy, ongoing work and studies will clarify soon if and how GRBs will be able to keep up to the promises.

We consider the effects of "photon bubble" shock trains on the vertical structure of radiation pressure-dominated accretion disks. These density inhomogeneities are expected to develop spontaneously in radiation-dominated accretion disks where magnetic pressure exceeds gas pressure, even in the presence of magnetorotational instability. They increase the rate at which radiation escapes from the disk, and may allow disks to exceed the Eddington limit by a substantial factor. We first generalize the theory of photon bubbles to include the effects of finite optical depths and radiation damping. Modifications to the diffusion law at low optical depth tend to fill in the low-density regions of photon bubbles, while radiation damping inhibits the formation of photon bubbles at large radii, small accretion rates, and small heights above the equatorial plane. Accretion disks dominated by photon bubble transport may reach luminosities of 10 to >100 times the Eddington limit (L_E), depending on the mass of the central object, while remaining geometrically thin. However, photon bubble-dominated disks with alpha-viscosity are subject to the same thermal and viscous instabilities that plague standard radiation pressure-dominated disks, suggesting that they may be intrinsically unsteady. Photon bubbles can lead to a "core-halo" vertical disk structure. In super-Eddington disks the halo forms the base of a wind, which carries away substantial energy and mass, but not enough to prevent the luminosity from exceeding L_E. Photon bubble-dominated disks may have smaller color corrections than standard accretion disks of the same luminosity. They remain viable contenders for some ultraluminous X-ray sources and may play a role in the rapid growth of supermassive black holes at high redshift.

We study the early afterglows of gamma-ray bursts produced by geometrically thick fireballs, following the development of the external shock as energy is continually supplied to the shocked material. We study the dependence of the early afterglow slope on the luminosity history of the central engine. The resulting light curves are modeled with power-law functions and the importance of a correct choice of the reference time t_0 is investigated. We find that deviations from a simple power-law are observed only if a large majority of the energy is released at late times. The light curve in this case can be described as a simple power-law if the reference time is set to be close to the end of the burst. We applied our analysis to the cases of GRB 050219a and GRB 050315. We show that the early steep decay of the afterglow cannot result from the interaction of the fireball with the ambient medium. We conclude that the early X-ray afterglow emission is associated with the prompt phase and we derive limits on the radius at which the prompt radiation is produced.

If quasar jets are accelerated by magnetic fields but terminate as matter dominated, where and how does the transition occur between the Poynting-dominated and matter-dominated regimes? To address this question, we study constraints which are imposed on the jet structure by observations at different spatial scales. We demonstrate that observational data are consistent with a scenario where the acceleration of a jet occurs within 10^{3-4} R_g. In this picture, the non-thermal flares -- important defining attributes of the blazar phenomenon - are produced by strong shocks formed in the region where the jet inertia becomes dominated by matter. Such shocks may be formed due to collisions between the portions of a jet accelerated to different velocities, and the acceleration differentiation is very likely to be related to global MHD instabilities.

We derive an analytic model for nonlinear "photon bubble" wave trains driven by buoyancy forces in magnetized, radiation pressure-dominated atmospheres. Continuous, periodic wave solutions exist when radiative diffusion is slow compared to the dynamical timescale of the atmosphere. We identify these waves with the saturation of a linear instability discovered by Arons - therefore, these wave trains should develop spontaneously. The buoyancy-driven waves are physically distinct from photon bubbles in the presence of rapid diffusion, which evolve into trains of gas pressure-dominated shocks as they become nonlinear. Like the gas pressure-driven shock trains, buoyancy-driven photon bubbles can exhibit very large density contrasts, which greatly enhance the flow of radiation through the atmosphere. However, steady-state solutions for buoyancy-driven photon bubbles exist only when an extra source of radiation is added to the energy equation, in the form of a flux divergence. We argue that this term is required to compensate for the radiation flux lost via the bubbles, which increases with height. We speculate that an atmosphere subject to buoyancy-driven photon bubbles, but lacking this compensating energy source, would lose pressure support and collapse on a timescale much shorter than the radiative diffusion time in the equivalent homogeneous atmosphere.

We argue that the velocity dispersions and masses of galactic bulges and spheroids are byproducts of the feedback that regulates rapid black hole growth in protogalaxies. We suggest that the feedback energy liberated by accretion must pass through the accreting material, in an energy-conserving flux close-in and a momentum-conserving flux further out. If the inflowing gas dominates the gravitational potential outside the Bondi radius, feedback from Eddington-limited accretion drives the density profile of the gas to that of a singular isothermal sphere. We find that the velocity dispersion associated with the isothermal potential, sigma, increases with time as the black hole mass M grows, in such a way that M is proportional to sigma^4. The coefficient of this proportionality depends on the radius at which the flow switches from energy conserving to momentum conserving, and gives the observed M-sigma relation if the transition occurs at ~100 Schwarzschild radii. We associate this transition with radiative cooling and show that bremsstrahlung, strongly boosted by inverse Compton scattering in a two-temperature (T_p >> T_e) plasma, leads to a transition at the desired radius. According to this picture, bulge masses M_b are insensitive to the virial masses of their dark matter haloes, but correlate linearly with black hole mass. Our analytic model also explains the M_b-sigma (Faber-Jackson) relation as a relic of black hole accretion. The model naturally explains why the M-sigma relation has less scatter than either the M-M_b (Magorrian) or the Faber-Jackson relation. It suggests that the M-sigma relation could extend down to very low velocity dispersions, and predicts that the relation should not evolve with redshift.

We consider the time-dependent evolution of a relativistic jet following its breakout through the surface of a massive compact star, as envisaged in the collapsar model of gamma-ray bursts. At breakout, the jet is tightly collimated by the pressure of its hot cocoon, which is created as the jet traverses the star. After breakout, the cocoon pressure drops and the jet expands toward its natural opening angle. We show that the evolving opening angle of the jet produces a stratification of the total energy with the off-axis angle, resulting in a universal morphology. The angular structure is largely independent of the initial beam pattern and depends only on the luminosity of the central engine. With the minimal assumption of a constant energy release we reproduce the theta^{-2} profile required to explain observations of afterglows.

After many years of speculation, recent observations have confirmed the association of gamma-ray bursts with core-collapse supernova explosions from massive stars. This association carries with it important consequences. The burst relativistic jet has to propagate through the cold dense stellar material before it reaches the transparency radius and the burst photons are produced. This propagation is likely to affect the initial properties of the jet, shaping it and changing its energy composition. The variability injected at the base of the jet is also likely to be erased by the jet-star interaction. Despite this, GRBs seem to have remarkably predictable properties once the radiative phase sets in, as emphasized by the recent discovery of several tight correlation between spectral, geometric and energetic properties of the jet. In this contribution we discuss the jet interaction with the star, emphasizing its time-dependent properties and the resulting energy distribution. We finally emphasize the surprising predictability of jet and radiation properties outside the star and underline its implication for standardizing the GRB candle.

We consider the production of electron cyclotron maser emission by low-density, highly magnetized plasmas in relativistic jets. The population inversion required to drive cyclotron maser instability could occur in localized, transient sites where hydromagnetic instabilities, shocks, and/or turbulence lead to magnetic mirroring along current-carrying flux tubes. The maser is pumped as electrons are accelerated by the parallel electric field that develops as a result of the mirror. We estimate the maximum brightness temperatures that can be obtained in a single maser site and in an array of many masers operating simultaneously, under conditions likely to apply in blazar jets. Synchrotron absorption, by relativistic electrons within the jet, presents the largest obstacle to the escape of the maser radiation, and may render most of it invisible. However, we argue that a high brightness temperature could be produced in a thin boundary layer outside the synchrotron photosphere, perhaps in the shear layer along the wall of the jet. Induced Compton scattering provides additional constraints on the maximum brightness temperature of a masing jet. We suggest that recent observations of diffractive scintillation in the blazar J1819+3845, indicating intrinsic brightness temperatures greater than 10^{14} K at 5 GHz, may be explained in terms of cyclotron maser emission. High brightness temperature maser emission from blazar jets may extend to frequencies as high as ~100 GHz, with the maximum possible T_B scaling roughly as 1/frequency. Less massive relativistic jet sources, such as microquasars, are even better candidates for producing cyclotron maser emission, primarily in the infrared and optical bands.

The formation of relativistic astrophysical jets is presumably mediated by magnetic fields threading accretion disks and central, rapidly rotating objects. As it is accelerated by magnetic stresses, the jet's kinetic energy flux grows at the expense of its Poynting flux. However, it is unclear how efficient is the conversion from magnetic to kinetic energy and whether there are any observational signatures of this process. We address this issue in the context of jets in quasars. Using data from all spatial scales, we demonstrate that in these objects the conversion from Poynting-flux-dominated to matter-dominated jets is very likely to take place closer to the black hole than the region where most of the Doppler boosted radiation observed in blazars is produced. We briefly discuss the possibility that blazar activity can be induced by global MHD instabilities, e.g., via the production of localized velocity gradients that lead to dissipative events such as shocks or magnetic reconnection, where acceleration of relativistic particles and production of non-thermal flares is taking place.

We propose that the knots of the large-scale jets in powerful radio sources represent moving and separate portions of the jet matter, with the excess kinetic power. This can explain many morphological properties of quasar jets - like high knot-to-interknot brightness contrasts, frequency-independent knot profiles and almost universal extents of the knot regions - independently of the exact emission mechanisms responsible for producing broad-band jet emission. We briefly discuss a possible connection between this scenario and the idea of highly modulated/intermittent activity of the jet engine.

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.

A general formalism is developed for calculating the luminosity function and the expected number $N$ of observed GRBs above a peak photon flux $S$ for any GRB jet structure. This new formalism directly provides the true GRB rate without the need for a `correction factor'. We apply it to the uniform jet (UJ) and universal structured jet (USJ) models for the structure of GRB jets and perform fits to the observed log(N)-log(S) distribution from the GUSBAD catalog which contains 2204 BATSE bursts. A core angle $\theta_c$ and an outer edge at $\theta_{max}$ are introduced for the structured jet, and a finite range of half-opening angles $\theta_{min}\leq\theta_j\leq\theta_{max}$ is assumed for the uniform jets. The efficiency $\epsilon_\gamma$ for producing gamma-rays, and the energy per solid angle $\epsilon$ in the jet are allowed to vary with $\theta_j$ (the viewing angle $\theta_{obs}$) in the UJ (USJ) model, $\epsilon_\gamma\propto\theta^{-b}$ and $\epsilon\propto\theta^{-a}$. We find that a single power-law luminosity function provides a good fit to the data. Such a luminosity function arises naturally in the USJ model, while in the UJ model it implies a power-law probability distribution for $\theta_j$, $P(\theta_j)\propto\theta_j^{-q}$. The value of $q$ cannot be directly determined from the fit to the observed log(N)-log(S) distribution, and an additional assumption on the value of $a$ or $b$ is required. Alternatively, an independent estimate of the true GRB rate would enable one to determine $a$, $b$ and $q$. The implied values of $\theta_c$ (or $\theta_{min}$) and $\theta_{max}$ are close to the current observational limits. The true GRB rate for the USJ model is found to be $R_{GRB}(z=0)=0.86^{+0.14}_{-0.05} Gpc^{-3} yr^{-1}$.

We present three-dimensional simulations of viscous dissipation of AGN induced gas motions and waves in clusters of galaxies. These simulations are motivated by recent detections of ripples in the Perseus and Virgo clusters. Although the sound waves generated by buoyant bubbles decay with distance from the cluster center, we show that these waves can contribute substantially to offsetting the radiative cooling at distances significantly exceeding the bubble size. The energy flux of the waves declines more steeply with radius than the inverse-square law predicted by energy conservation, implying that dissipation plays an important role in tapping the wave energy. We show that such dispersing sound waves/weak shocks are detectable as ripples on unsharp-masked X-ray cluster maps, and point out that the interfaces between the intracluster medium and old bubbles are also clearly detectable in unsharp-masked X-ray maps. This opens up the possibility of detecting fossil bubbles that are difficult to detect in radio emission. This mode of heating is consistent with other observational constraints, such as the presence of cool rims around the bubbles and the absence of strong shocks. Thus, the mechanism offers a way of heating clusters in a spatially distributed and gentle fashion. We also discuss the energy transfer between the central AGN and the surrounding medium. In our numerical experiments, we find that roughly 65 per cent of the energy injected by the AGN is transferred to the intracluster medium and approximately 25 percent of the injected energy is dissipated by viscous effects and contributes to heating of the gas. The overall transfer of heat from the AGN to the gas is comparable to the radiative cooling losses. The simulations were performed with the FLASH adaptive mesh refinement code.

Recent X-ray observations of clusters of galaxies have shown that the entropy of the intracluster medium (ICM), even at radii as large as half the virial radius, is higher than that expected from gravitational processes alone. This is thought to be the result of nongravitational processes influencing the physical state of the ICM. In this paper, we investigate whether heating by a central AGN can explain the distribution of excess entropy as a function of radius. The AGN is assumed to inject buoyant bubbles into the ICM, which heat the ambient medium by doing pdV work as they rise and expand. Several authors have suggested that this "effervescent heating" mechanism could allow the central regions of clusters to avoid the ``cooling catastrophe''. Here we study the effect of effervescent heating at large radii. Our calculations show that such a heating mechanism is able to solve the entropy problem. The only free parameters of the model are the time-averaged luminosity and the AGN lifetime. The results are mainly sensitive to the total energy injected into the cluster. Our model predicts that the total energy injected by AGN should be roughly proportional to the cluster mass. The expected correlation is consistent with a linear relation between the mass of the central black hole(s) and the mass of the cluster, which is reminiscent of the Magorrian relation between the black hole and bulge mass.

We argue that the observed `accretion rate gap' - between black holes in radio-loud active galactic nuclei (AGN) accreting at close to the Eddington limit and those accreting at considerably lower rates - can be explained in terms of the adiabatic inflow-outflow (ADIOS) scenario for radiatively inefficient accretion. Whenever the accretion rate falls below a threshold value (corresponding to a luminosity L_crit) that depends on the viscosity parameter, alpha, the inner region of the accretion disk - extending from the marginally stable orbit to ~ 1000 Schwarzschild radii - is susceptible to becoming hot and radiatively inefficient. If this happens, the disk luminosity decreases by a factor of ~100, as most of the matter originally destined to be swallowed is instead expelled in a wind. According to our conjecture, accretion flows onto black holes never radiate steadily in the range ~ 0.01 L_crit < L < L_crit, hence the inferred accretion rate gap. We expect the gap to exist also for black holes in X-ray binaries, where it may be responsible for state transitions and the luminosity fluctuations associated with X-ray nova outbursts.

We simulate the effects of viscous dissipation of waves that are generated by AGN activity in clusters of galaxies. We demonstrate that the amount of viscous heating associated with the dissipation of these waves can offset radiative cooling rates in cooling flow clusters of galaxies. This heating mechanism leads to spatially distributed and approximately symmetrical dissipation. The heating waves reach a given distance from the cluster center on a timescale shorter than the cooling time. This means that this heating mechanism has the potential of quenching cooling flows in a quasi-stable fashion. Moreover, the heating is gentle as no strong shocks are present in the simulations. We first investigated whether a single continuous episode of AGN activity can lead to adequate dissipation to balance cooling rates. These simulations demonstrated that, whereas secondary waves generated by the interaction of the rising bubble with the intracluster medium are clearly present, viscous heating associated with the dissipation of these waves is insufficient to balance radiative cooling. It is only when the central source is intermittent that the viscous dissipation of waves associated with subsequent episodes of activity can offset cooling. This suggests that the ripples observed in the Perseus cluster can be interpreted as being due to the AGN duty cycle, i.e., they trace AGN activity history. The simulations were performed using the PPM adaptive mesh refinement code FLASH in two dimensions.

Jets and winds are significant channels for energy loss from accreting black holes. These outflows mechanically heat their surroundings, through shocks as well as gentler forms of heating. We discuss recent efforts to understand the nature and distribution of mechanical heating by central AGNs in clusters of galaxies, using numerical simulations and analytic models. Specifically, we will discuss whether the relatively gentle `effervescent heating' mechanism can compensate for radiative losses in the central regions of clusters, and account for the excess entropy observed at larger radii.

We discuss morphological properties of the large-scale jets in powerful radio sources, which are now observed at radio, optical and X-ray frequencies, in order to determine the origin of their X-ray radiation and the nature of the particle acceleration responsible for their multiwavelength emission. We show that modeling knots of the these objects as stationary regions of energy dissipation within the uniform and continuous jet flow leads to several problems. Such problems are especially pronounced if the observed X-ray emission is due to inverse-Compton scattering of the CMB radiation. We explore another possibility, namely that the knots represent moving and separate portions of the jet matter, with excess kinetic power. We suggest a possible connection between this scenario and the idea of intermittent/highly modulated jet activity. The proposed model can explain some morphological properties of quasar jets -- like high knot-to-interknot brightness contrasts, frequency-independent knot profiles and almost universal extents of the knot regions -- independently of the exact emission mechanism responsible for producing the X-rays. In this context, we consider different possibilities for the production of such X-ray radiation, and discuss the related issue of particle acceleration. We conclude that the appropriate process cannot be designated for certain yet. We suggest, however, that X-ray observations of the jet in 3C 120 (in which intermittent jet activity is quite distinct) already seem to support the synchrotron origin of the X-ray emission, at least in this object.

(abridged) We present a detailed spectral analysis of the first observation of the Seyfert 1 galaxy MCG-6-30-15 by the European Photon Imaging Camera on board the XMM-Newton observatory, together with contemporaneous data from the Proportional Counter Array on the Rossi X-ray Timing Explorer. Confirming our previously published result, we find that the presence of extremely broadened reflection features from an ionized relativistic accretion disk is required even when one employs the latest X-ray reflection models and includes the effect of complex absorption. The extremely broadened reflection features are also present if the primary continuum is modeled with a thermal Comptonisation spectrum rather than a simple power-law continuum. With this fact established, we examine these data using a relativistic smearing function corresponding to a ``generalized thin accretion disk'' model. We find strong evidence for torquing of the central parts of the accretion disk (presumably through magnetic interactions with the plunging region of the disk and/or the rotating black hole itself). We also perform a study of spectral variability within our observation. We find that the disk reflection features maintain roughly a constant equivalent width with respect to the observed continuum, as predicted by simple reflection models. Taken together with other studies of MCG-6-30-15 that find disk features to possess constant intensity at higher flux states, we suggest that the flux of disk features undergoes a saturation once the source emerges from a Deep Minimum state.

We present results from axisymmetric, time-dependent magnetohydrodynamic (MHD) simulations of the collapsar model for gamma-ray bursts. Our main conclusion is that, within the collapsar model, MHD effects alone are able to launch, accelerate and sustain a strong polar outflow. We also find that the outflow is Poynting flux-dominated, and note that this provides favorable initial conditions for the subsequent production of a baryon-poor fireball.

We present results from axisymmetric, time-dependent magnetohydrodynamic (MHD) simulations of the collapsar model for gamma-ray bursts. We begin the simulations after the 1.7 MSUN iron core of a 25 MSUN presupernova star has collapsed and study the ensuing accretion of the 7 MSUN helium envelope onto the central black hole formed by the collapsed iron core. We consider a spherically symmetric progenitor model, but with spherical symmetry broken by the introduction of a small, latitude-dependent angular momentum and a weak radial magnetic field. Our MHD simulations include a realistic equation of state, neutrino cooling, photodisintegration of helium, and resistive heating. Our main conclusion is that, within the collapsar model, MHD effects alone are able to launch, accelerate and sustain a strong polar outflow. We also find that the outflow is Poynting flux-dominated, and note that this provides favorable initial conditions for the subsequent production of a baryon-poor fireball.

We propose a gamma-ray burst scenario involving relativistic jets dominated by Poynting flux with alternating toroidal magnetic fields. Such a structure may arise naturally if the jet is formed and powered by the accretion flow in the core of a collapsar. We conjecture that the polarity of the toroidal magnetic field changes randomly due to hydromagnetic turbulence driven by the magneto-rotational instability (MRI), with the typical reversal time determined by the time scale for amplifying magnetic fields up to a dynamically important level. Poynting flux-dominated jets with reversing B-fields provide a natural and efficient way to dissipate energy via the reconnection process. Gamma-rays are produced at the spatially separated reconnection sites. In this scenario, the emergent synchrotron radiation can be highly polarized and can form both smooth and spiky light curves. We note the possibility that cold and dense filaments can form in the reconnection zones as the result of thermal instability. One could then explain the production of very hard X-ray spectra, as due to bound-free absorption of the synchrotron radiation.

When gas accretes onto a black hole, at a rate either much less than or much greater than the Eddington rate, it is likely to do so in an "adiabatic" or radiatively inefficient manner. Under fluid (as opposed to MHD) conditions, the disk should become convective and evolve toward a state of marginal instability. The resulting disk structure is "gyrentropic," with convection proceeding along common surfaces of constant angular momentum, Bernoulli function and entropy, called "gyrentropes." We present a family of two-dimensional, self-similar models which describes the time-averaged disk structure. We then suppose that there is a self-similar, Newtonian torque and that the Prandtl number is large. This torque drives inflow and meridional circulation and the resulting flow is computed. Convective transport will become ineffectual near the disk surface. It is conjectured that this will lead to a large increase of entropy across a "thermal front" which we identify as the effective disk surface and the base of an outflow. The conservation of mass, momentum and energy across this thermal front permits a matching of the disk models to self-similar outflow solutions. We then demonstrate that self-similar disk solutions can be matched smoothly onto relativistic flows at small radius and thin disks at large radius. This model of adiabatic accretion is contrasted with some alternative models that have been discussed recently. The disk models developed in this paper should be useful for interpreting numerical, fluid dynamical simulations. Related principles to those described here may govern the behaviour of astrophysically relevant, magnetohydrodynamic disk models.

We report on the second phase of our study of slightly rotating accretion flows onto black holes. We consider magnetohydrodynamical (MHD) accretion flows with a spherically symmetric density distribution at the outer boundary, but with spherical symmetry broken by the introduction of a small, latitude-dependent angular momentum and a weak radial magnetic field. We study accretion flows by means of numerical 2D, axisymmetric, MHD simulations with and without resistive heating. Our main result is that the properties of the accretion flow depend mostly on an equatorial accretion torus which is made of the material that has too much angular momentum to be accreted directly. The torus accretes, however, because of the transport of angular momentum due to the magnetorotational instability (MRI). Initially, accretion is dominated by the polar funnel, as in the hydrodynamic inviscid case, where material has zero or very low angular momentum. At the later phase of the evolution, the torus thickens towards the poles and develops a corona or an outflow or both. Consequently, the mass accretion through the funnel is stopped. The accretion of rotating gas through the torus is significantly reduced compared to the accretion of non-rotating gas (i.e., the Bondi rate). It is also much smaller than the accretion rate in the inviscid, weakly rotating case.Our results do not change if we switch on or off resistive heating. Overall our simulations are very similar to those presented by Stone, Pringle, Hawley and Balbus despite different initial and outer boundary conditions. Thus, we confirm that MRI is very robust and controls the nature of radiatively inefficient accretion flows.

Accreting black holes can release enormous amounts of energy to their surroundings, in various forms. Such feedback may profoundly influence a black hole's environment. After briefly reviewing the possible types of feedback, I focus on the injection of kinetic energy through jets and powerful winds. The effects of these outflows may be especially apparent in the heating of the X-ray--emitting atmospheres that pervade clusters of galaxies. Analogous heating effects, during the epoch of galaxy formation, could regulate the growth of supermassive black holes.

The masses of supermassive black holes (SBHs) show correlations with bulge properties in disk and elliptical galaxies. We study the formation of galactic structure within flat-core triaxial haloes and show that these correlations can be understood within the framework of a baryonic component modifying the orbital structure in the underlying potential. In particular, we find that terminal properties of bulges and their central SBHs are constrained by the destruction of box orbits in the harmonic cores of dark haloes and the emergence of progressively less eccentric loop orbits there. SBH masses, M_SBH, should exhibit a tighter correlation with bulge velocity dispersions, sigma_B, than with bulge masses, M_B, in accord with observations, if there is a significant scatter in the M_H-sigma_H relation for the halo. In the context of this model the observed M_SBH-sigma_B relation implies that haloes should exhibit a Faber-Jackson type relationship between their masses and velocity dispersions. The most important prediction of our model is that halo properties determine the bulge and SBH parameters. The model also has important implications for galactic morphology and the process of disk formation.

Radiation dominated accretion disks are likely to be subject to the ``photon bubble'' instability, which may lead to strong density inhomogeneities on scales much shorter than the disk scale height. Such disks -- and magnetized, radiation-dominated atmospheres in general -- could radiate well above the Eddington limit without being disrupted. When density contrasts become large over distances of order the photon mean free path, radiative transfer cannot be described adequately using either the standard diffusion approximation or existing prescriptions for flux-limited diffusion. Using analytical and Monte Carlo techniques, we consider the effects of strong density gradients deep within radiation- and scattering-dominated atmospheres. We find that radiation viscosity -- i.e., the off-diagonal elements of the radiation stress tensor -- has an important effect on radiative transfer under such conditions. We compare analytical and numerical results in the specific case of a plane-parallel density wave structure and calculate Eddington enhancement factors due to the porosity of the atmosphere. Our results can be applied to the study of dynamical coupling between radiation forces and density inhomogeneities in radiation dominated accretion disks in two or three dimensions.

We report on the first phase of our study of slightly rotating accretion flows onto black holes. We consider inviscid accretion flows with a spherically symmetric density distribution at the outer boundary, but with spherical symmetry broken by the introduction of a small, latitude-dependent angular momentum. We study accretion flows by means of numerical 2D, axisymmetric, hydrodynamical simulations. Our main result is that the properties of the accretion flow do not depend as much on the outer boundary conditions (i.e., the amount as well as distribution of the angular momentum) as on the geometry of the non-accreting matter. The material that has too much angular momentum to be accreted forms a thick torus near the equator. Consequently, the geometry of the polar region, where material is accreted (the funnel), and the mass accretion rate through it are constrained by the size and shape of the torus. Our results show one way in which the mass accretion rate of slightly rotating gas can be significantly reduced compared to the accretion of non-rotating gas (i.e., the Bondi rate), and set the stage for calculations that will take into account the transport of angular momentum and energy.

There is mounting observational evidence from Chandra for strong interaction between keV gas and AGN in cooling flows. It is now widely accepted that the temperatures of cluster cores are maintained at a level of 1 keV and that the mass deposition rates are lower than earlier ROSAT/Einstein values. Recent theoretical results suggest that thermal conduction can be very efficient even in magnetized plasmas. Motivated by these discoveries, we consider a ``double heating model'' which incorporates the effects of simultaneous heating by both the central AGN and thermal conduction from the hot outer layers of clusters. Using hydrodynamical simulations, we demonstrate that there exists a family of solutions that does not suffer from the cooling catastrophe. In these cases, clusters relax to a stable final state, which is characterized by minimum temperatures of order 1 keV and density and temperature profiles consistent with observations. Moreover, the accretion rates are much reduced, thereby reducing the need for excessive mass deposition rates required by the standard cooling flow models.

The time- and ensemble-averaged mechanical energy outputs of radio galaxies may be large enough to offset much of the cooling inferred from X-ray observations of galaxy clusters. But does this heating actually counterbalance the cooling, diminishing cooling flows or quenching them altogether? I will argue that energy injection by radio galaxies may be important even in clusters where no active source is present, due to the likely intermittency of the jets. If the energy injected by radio galaxies percolates through the intracluster medium without excessive mixing, it could stabilize the atomic cooling responsible for X-ray emission.

During a previous outburst phase, the neutron star 4U 2129+47 exhibited evidence for a spatially extended corona via broad, partial X-ray eclipses occurring periodically on the 5.24 hr orbit. Since 1983, however, 4U 2129+47 has entered a quiescent state several orders of magnitude fainter. We have performed a 37 ksec Chandra observation of 4U 2129+47 to determine whether an extended coronal structure also exists in quiescence. Total eclipses are found, and the rapidity of the eclipse ingress and egress are used to place upper limits on the size of the X-ray source. The spectrum is comprised of a soft component plus a fairly hard power law tail. The former is seen to be sinusoidally modulated over the orbital period in a manner consistent with neutral column variations, possibly due to the interaction of an accretion stream with a (small) disk about the neutron star. We fit realistic atmosphere models to the soft spectra, and comment on the consistency of the spectra with cooling neutron star models. It has been suggested that the 4U 2129+47 system is a hierarchical triple, with the outer body being an F star. We use differential astrometry to show that the X-ray point source and F star are spatially coincident to within 0.1". We further compare newly determined upper limits for the extrinsic neutral column to the reddening of the F star. Finally, we discuss how future X-ray observations can further constrain models of quiescent neutron star emission, as well as directly verify the triple hypothesis.

Radiation pressure-dominated accretion disks are predicted to exhibit strong density inhomogeneities on scales much smaller than the disk scale height, due to the nonlinear development of photon bubble instability. Radiation would escape from such a "leaky" disk at a rate higher than that predicted by standard accretion disk theory. The disk scale height is then smaller than that of a similar disk without small-scale inhomogeneities, and the disk can remain geometrically thin even as the flux approaches and exceeds the Eddington limit. An idealized one-zone model for disks with radiation-driven inhomogeneities suggests that the escaping flux could exceed L_E by a factor of up to ~10-100, depending on the mass of the central object. Such luminous disks would develop strong mass loss, but the resulting decrease in accretion rate would not necessarily prevent the luminosity from exceeding L_E. We suggest that the observed "ultraluminous X-ray sources" are actually thin, super-Eddington accretion disks orbiting stellar-mass black holes, and need not indicate the existence of a class of intermediate-mass black holes.

We present a numerical investigation of dead, or relic, radio galaxies and the environmental impact that radio galaxy activity has on the host galaxy or galaxy cluster. We perform axisymmetric hydrodynamical calculations of light, supersonic, back-to-back jets propagating in a beta-model galaxy/cluster atmosphere. We then shut down the jet activity and let the resulting structure evolve passively. The dead source undergoes an initial phase of pressure driven expansion until it achieves pressure equilibrium with its surroundings. Thereafter, buoyancy forces drive the evolution and lead to the formation of two oppositely directed plumes that float high into the galaxy/cluster atmosphere. These plumes entrain a significant amount of low entropy material from the galaxy/cluster core and lift it high into the atmosphere. An important result is that a large fraction (at least half) of the energy injected by the jet activity is thermalized in the ISM/ICM core. The whole ISM/ICM atmosphere inflates in order to regain hydrostatic equilibrium. This inflation is mediated by an approximately spherical disturbance which propagates into the atmosphere at the sound speed. The fact that such a large fraction of the injected energy is thermalized suggests that radio galaxies may have an important role in the overall energy budget of rich ISM/ICM atmospheres. In particular, they may act as a strong and highly time-dependent source of negative feedback for galaxy/cluster cooling flows.

We present Chandra observations of the galaxy cluster A4059. We find strong evidence that the FR-I radio galaxy PKS 2354-35 at the center of A4059 is inflating cavities with radii ~20 kpc in the intracluster medium, similar to the situation seen in Perseus A and Hydra A. We also find evidence for interaction between the intracluster medium and PKS 2354-35 on small scales in the very center of the cluster. Arguments are presented suggesting that this radio galaxy has faded significantly in radio power (possibly from an FR-II state) over the past 10^8 yrs.

The transfer of polarized radiation in stochastic synchrotron sources is explored by means of analytic treatment and Monte Carlo simulations. We argue that the main mechanism responsible for the circular polarization properties of compact synchrotron sources is likely to be Faraday conversion and that, contrary to common expectation, a significant rate of Faraday rotation does not necessarily imply strong depolarization. The long-term persistence of the sign of circular polarization, observed in many sources, is most likely due to a small net magnetic flux generated in the central engine, carried along the jet axis and superimposed on a highly turbulent magnetic field. We show that the mean levels of circular and linear polarizations depend on the number of field reversals along the line of sight and that the gradient in Faraday rotation across turbulent regions can lead to ``correlation depolarization''. Our model is potentially applicable to a wide range of synchrotron sources. In particular, we demonstrate how our model can naturally explain the excess of circular over linear polarization in the Galactic Center and the nearby spiral galaxy M81 and discuss its application to the quasar 3C 279, the intraday variable blazar PKS 1519-273 and the X-ray binary SS 433.

We present XMM-Newton European Photon Imaging Camera (EPIC) observations of the bright Seyfert 1 galaxy MCG-6-30-15, focusing on the broad Fe K$\alpha$ line at ~6keV and the associated reflection continuum, which is believed to originate from the inner accretion disk. We find these reflection features to be extremely broad and red-shifted, indicating its origin from the very most central regions of the accretion disk. It seems likely that we have caught this source in the ``deep minimum'' state first observed by Iwasawa et al. (1996). The implied central concentration of X-ray illumination is difficult to understand in any pure accretion disk model. We suggest that we are witnessing the extraction and dissipation of rotational energy from a spinning black hole by magnetic fields connecting the black hole or plunging region to the disk.

Theories of high energy radiation production in quasar jets can be verified by studies of both time-averaged spectra and variability patterns. While the former has been explored extensively, the latter is in its infancy. In this paper, we study the production of short-term flares in the shock-in-jet model. We examine how the flares' profiles depend on such parameters as shock/dissipation lifetime, electron-injection time profile, adiabaticity, and half-opening angle of the jet. In particular, we demonstrate the large difference between flare profiles produced in the radiative and adiabatic regimes. We apply our model to the $\sim$day timescale flares observed in optically violently variable (OVV) quasars, checking whether the external-radiation-Compton (ERC) model for $\gamma$-ray flares at energies $> 30$ MeV (EGRET range) can be reconciled with the flares observed at lower energies. Specifically, we show that the strict correlation between X-ray and $\gamma$-ray flares strongly supports the dominance of the synchrotron self-Compton mechanism in the X-ray band. We also derive conditions that must be satisfied by the ERC model in order to explain a lag of the $\gamma$-ray peak behind the optical one, as claimed to be observed in PKS 1406-076. Finally, we predict that in ERC models where the MeV peak is related to the break in electron distribution due to inefficient cooling of electrons below a certain energy, the flares should decay significantly more slowly in the soft $\gamma$-ray band than at energies greater than 30 MeV.

We show that magnetized, radiation dominated atmospheres can support steady state patterns of density inhomogeneity that enable them to radiate at far above the Eddington limit, without suffering mass loss. The inhomogeneities consist of periodic shock fronts bounding narrow, high-density regions, interspersed with much broader regions of low density. The radiation flux avoids the regions of high density, which are therefore weighed down by gravity, while gas in the low-density regions is slammed upward into the shock fronts by radiation force. As the wave pattern moves through the atmosphere, each parcel of matter alternately experiences upward and downward forces, which balance on average. Magnetic tension shares the competing forces between regions of different densities, preventing the atmosphere from blowing apart. We calculate the density structure and phase speed of the wave pattern, and relate these to the wavelength, the density contrast, and the factor by which the net radiation flux exceeds the Eddington limit. In principle, this factor can be as large as the ratio of magnetic pressure to mean gas pressure, or the ratio of radiation pressure to gas pressure, whichever is smaller. Although the magnetic pressure must be large compared to the mean gas pressure in order to support a large density contrast, it need not be large compared to the radiation pressure. These highly inhomogeneous flows could represent the nonlinear development of the "photon bubble" instability discovered by Gammie. We briefly discuss the applicability of these solutions to astrophysical systems.

Motivated by hydrodynamic simulations, we discuss the X-ray appearance of radio galaxies embedded in the intracluster medium (ICM) of a galaxy cluster. We distinguish three regimes. In the early life of a powerful source, the entire radio cocoon is expanding supersonically and hence drives a strong shock into the ICM. Eventually, the sides of the cocoon become subsonic and the ICM is disturbed by the sonic booms of the jet's working surface. In both of these regimes, X-ray observations would find an X-ray shell. In the strong shock regime, this shell will be hot and relatively thin. However, in the weak shock (sonic-boom) regime, the shell will be approximately the same temperature as the undisturbed ICM. If a cooling flow is present, the observed shell may even be cooler than the undisturbed ICM due to the lifting of cooler material into the shell from the inner (cooler) regions of the cluster. In the third and final regime, the cocoon has collapsed and no well-defined X-ray shell will be seen. We discuss ways of estimating the power and age of the source once its regime of behavior has been determined.

We present results from 170ksec long RXTE observations of LMC X-1 and LMC X-3, taken in 1996 December, where their spectra can be described by a disc black body plus an additional soft (Gamma~2.8) high-energy power-law (detected up to 50keV in LMC X-3). These observations, as well as archival ASCA observations, constrain any narrow Fe line present in the spectra to have an equivalent width <90eV, broad lines (~150eV EW, \sigma ~ 1keV) are permitted. We also study the variability of LMC X-1. Its X-ray power spectral density (PSD) is approximately f^{-1} between 10^{-3} and 0.3Hz with a rms variability of ~7%. Above 5keV the PSD shows evidence of a break at f > 0.2Hz, possibly indicating an outer disc radius of ~1000GM/c^2 in this likely wind-fed system. Furthermore, the coherence function between variability in the > 5keV band and variablity in the lower energy bands is extremely low. We discuss the implications of these observations for the mechanisms.

We present the analysis of an RXTE monitoring campaign of the canonical soft state black hole candidates LMC X-1 and LMC X-3. In contrast to LMC X-1, which does not exhibit any periodic spectral changes, we find that LMC X-3 exhibits stron spectral variability on time scales of days to weeks. For typical RXTE ASM count rates, the luminosity variations of LMC X-3 are due to changes of the phenomenological disk blackbody temperature. During episodes of especially low luminosity (ASM count rates < 0.6 counts/sec), kT strongly decreases and the power law significantly hardens to a photon index of ~1.8. These changes are consistent with state changes of LMC X-3 from the soft state to the canonical hard state of galactic black hole candidates. We argue that the long term variability of LMC X-3 might be due to a wind-driven limit cycle such as discussed by Shields et al. (1986)

We explore the possibility that extragalactic radio jets might be accelerated by highly disorganized magnetic fields that are strong enough to dominate the dynamics until the terminal Lorentz factor is reached. Following the twin-exhaust model by Blandford & Rees (1974), the collimation under this scenario is provided by the stratified thermal pressure from an external medium. The acceleration efficiency then depends on the pressure gradient of that medium. In order for this mechanism to work there must be continuous tangling of the magnetic field, changing the magnetic equation of state away from pure flux freezing (otherwise conversion of Poynting flux to kinetic energy flux is suppressed). This is a complementary approach to models in which the plasma is accelerated by large scale ordered fields. We include a simple prescription for magnetic dissipation, which leads to tradeoffs among conversion of magnetic energy into bulk kinetic energy, random particle energy, and radiation. We present analytic dynamical solutions of such jets, assess the effects of radiation drag, and comment on observational issues, such as the predicted polarization and synchrotron brightness. Finally, we try to make the connection to observed radio galaxies and gamma-ray bursts.

We consider possible evolutionary models for SS 433. We assume that common-envelope evolution is avoided if radiation pressure is able to expel most of a super-Eddington accretion flow from a region smaller than the accretor's Roche lobe. This condition is satisfied, at least initially, for largely radiative donors with masses in the range 4-12 solar masses. For donors more massive than about 5 solar masses, moderate mass ratios q = M_2/M_1 > 1 are indicated, thus tending to favor black-hole accretors. For lower mass donors, evolutionary considerations do not distinguish between a neutron star or black hole accretor. In all cases the mass transfer (and mass loss) rates are much larger than the likely mass-loss rate in the precessing jets. Almost all of the transferred mass is expelled at radii considerably larger than the jet acceleration region, producing the "stationary" H-alpha line, the infrared luminosity, and accounting for the low X-ray luminosity.

The black hole candidate LMC X-3 varies by a factor of four in the soft X-rays on a timescale of 200 or 100days. We have monitored LMC X-3 with RXTE in three to four week intervals starting in December 1996, obtaining a large observational database that sheds light on the nature of the long term X-ray variability in this source. In this paper we present the results of this monitoring campaign, focusing on evidence of recurring hard states in this canonical soft state black hole candidate.

We propose that gamma ray bursts (GRBs) are produced by a shower of heavy blobs running into circumstellar material at highly relativistic speeds. The gamma ray emission is produced in the shocks these bullets drive into the surrounding medium. The short term variability seen in GRBs is set by the slowing-down time of the bullets while the overall duration of the burst is set by the lifetime of the central engine. A requirement of this model is that the ambient medium be dense, consistent with a strong stellar wind. The efficiency of the burst can be relatively high.

Two-dimensional (axially symmetric) numerical hydrodynamical calculations of accretion flows which cannot cool through emission of radiation are presented. The calculations begin from an equilibrium configuration consisting of a thick torus with constant specific angular momentum. Accretion is induced by the addition of a small anomalous azimuthal shear stress which is characterized by a function \nu. We study the flows generated as the amplitude and form of \nu are varied. A spherical polar grid which spans more than two orders of magnitude in radius is used to resolve the flow over a wide range of spatial scales. We find that convection in the inner regions produces significant outward mass motions that carry away both the energy liberated by, and a large fraction of the mass participating in, the accretion flow. Although the instantaneous structure of the flow is complex and dominated by convective eddies, long time averages of the dynamical variables show remarkable correspondence to certain steady-state solutions. Near the equatorial plane, the radial profiles of the time-averaged variables are power-laws with an index that depends on the radial scaling of the shear stress. We find that regardless of the adiabatic index of the gas, or the form or magnitude of the shear stress, the mass inflow rate is a strongly increasing function of radius, and is everywhere nearly exactly balanced by mass outflow. The net mass accretion rate through the disc is only a fraction of the rate at which mass is supplied to the inflow at large radii, and is given by the local, viscous accretion rate associated with the flow properties near the central object.

Recent work on Cygnus X-2 suggests that neutron-star or black-hole binaries survive highly super-Eddington mass transfer rates without undergoing common-envelope evolution. We suggest here that the accretion flows in such cases are radiation pressure-dominated versions of the "ADIOS" picture proposed by Blandford and Begelman (1999), in which almost all the mass is expelled from large radii in the accretion disk. We estimate the maximum radius from which mass loss is likely to occur, and show that common-envelope evolution is probably avoided in any binary in which a main-sequence donor transfers mass on a thermal timescale to a neutron star or black hole, even though the mass transfer rate may reach values of 0.001 solar masses per year. This conclusion probably applies also to donors expanding across the Hertzsprung gap, provided that their envelopes are radiative. SS433 may be an example of a system in this state.

We present hard X-ray observations of the nearby radio galaxy M87 and the core of the Virgo cluster using the Rossi X-ray Timing Explorer. These are the first hard X-ray observations of M87 not affected by contamination from the nearby Seyfert 2 galaxy NGC4388. Thermal emission from Virgo's intracluster medium is clearly detected and has a spectrum indicative of kT=2.5keV plasma with approximately 25% cosmic abundances. No non-thermal (power-law) emission from M87 is detected in the hard X-ray band, with fluctuations in the Cosmic X-ray Background being the limiting factor. Combining with ROSAT data, we infer that the X-ray spectrum of the M87 core and jet must be steep (Gamma_core>1.90$ and Gamma_jet>1.75), and we discuss the implications of this result. In particular, these results are consistent with M87 being a mis-aligned BL-Lac object.

We have recently shown that a `sphere+disk' geometry Compton corona model provides a good description of Rossi X-ray Timing Explorer (RXTE) observations of the hard/low state of Cygnus X--1. Separately, we have analyzed the temporal data provided by RXTE. In this paper we consider the implications of this timing analysis for our best-fit `sphere+disk' Comptonization models. We focus our attention on the observed Fourier frequency-dependent time delays between hard and soft photons. We consider whether the observed time delays are: created in the disk but are merely reprocessed by the corona; created by differences between the hard and soft photon diffusion times in coronae with extremely large radii; or are due to `propagation' of disturbances through the corona. We find that the time delays are most likely created directly within the corona; however, it is currently uncertain which specific model is the most likely explanation. Models that posit a large coronal radius [or equivalently, a large Advection Dominated Accretion Flow (ADAF) region] do not fully address all the details of the observed spectrum. The Compton corona models that do address the full spectrum do not contain dynamical information. We show, however, that simple phenomenological propagation models for the observed time delays for these latter models imply extremely slow characteristic propagation speeds within the coronal region.

We analyze HST and ground based spectra of the brightest BALQSO in the UV: PG~0946+301. A detailed study of the absorption troughs as a function of velocity is presented, facilitated by the use of a new algorithm to solve for the optical depth as a function of velocity for multiplet lines. We find convincing evidence for saturation in parts of the troughs. This supports our previous assertion that saturation is common in BALs and therefore cast doubts on claims for very high metallicity in BAL flows. Due to the importance of BAL saturation we also discuss its evidence in other objects. In PG~0946+301 large differences in ionization as a function of velocity are detected and our findings supports the hypothesis that the line of sight intersects a number of flow components that combine to give the appearance of the whole trough. Based on the optical depth profiles, we develop a geometrical-kinematical model for the flow. We have positively identified 16 ions of 8 elements (H~I, C~III, C~IV, N~III, N~IV, N~V, O~III, O~IV, O~V, O~VI, Ne~V, Ne~VIII, P~V, Si~IV, S~V, S~VI) and have a probable identifications of Mg~X and S~IV. Unlike earlier analysis of IUE data, we find no evidence for BALs arising from excited ionic states in the HST spectrum of PG~0946+301.

We present results from a 20 ksec RXTE observation of the black hole candidate Cyg X-1. We apply self-consistent accretion disk corona models to these hard state data and show that Comptonization in a spherical corona irradiated by soft photons from an exterior cold disk is able to successfully model the spectrum. We also present the power spectrum, the coherence function, and the time lags for lightcurves from four energy bands. By modeling the high-resolution lightcurves with stochastic linear state space models, we show that the rapid hard state variability of Cyg X-1 can be explained with a single timescale.

Of all known persistent stellar-mass black hole candidates, only LMC X-1 and LMC X-3 consistently show spectra that are dominated by a soft, thermal component. We present results from long (170ksec) Rossi X-ray Timing Explorer (RXTE) observations of LMC X-1 and LMC X-3 made in 1996 December. The spectra can be described by a multicolor disk blackbody plus an additional high-energy power-law. Even though the spectra are very soft (Gamma is about 2.5), RXTE detected a significant signal from LMC X-3 up to energies of 50keV, the hardest energy at which the object was ever detected. Focusing on LMC X-3, we present results from the first year of an ongoing monitoring campaign with RXTE which started in 1997 January. We show that the appearance of the object changes considerably over its ~200d long cycle. This variability can either be explained by periodic changes in the mass transfer rate or by a precessing accretion disk analogous to Her X-1.

I propose that the moving `wisps' near the center of the Crab Nebula result from nonlinear Kelvin-Helmholtz instabilities in the equatorial plane of the shocked pulsar wind. Recent observations suggest that the wisps trace out circular wavefronts in this plane, expanding radially at speeds approximately less than c/3. Instabilities could develop if there is sufficient velocity shear between a faster-moving equatorial zone and a slower moving shocked pulsar wind at higher latitudes. The development of shear could be related to the existence of a neutral sheet -- with weak magnetic field -- in the equatorial zone, and could also be related to a recent suggestion by Begelman that the magnetic field in the Crab pulsar wind is much stronger than had been thought. I show that plausible conditions could lead to the growth of instabilities at the radii and speeds observed, and that their nonlinear development could lead to the appearance of sharp wisplike features.

The relativistically broad X-ray iron line seen in many AGN spectra is thought to originate from the central regions of the putative black hole accretion disk. Both the line profile and strength will vary in response to rapid variability of the primary X-ray continuum source. The temporal response of the line contains information on the accretion disk structure, the X-ray source geometry, and the spin of the black hole. Since the X-ray source will have a size comparable to the fluorescing region of the accretion disk, the general reverberation problem is not invertible. However, progress can be made since, empirically, AGN light curves are seen to undergo dramatic short timescale variability which presumably corresponds to the creation of a single new active region within the distributed X-ray source. The iron line response to these individual events can be described using linear transfer theory. We consider the line response to the activation/flaring of a new X-ray emitting region. Most of our detailed calculations are performed for the case of an X-ray source on the symmetry axis and at some height above the disk plane around a Kerr black hole. We also present preliminary calculations for off-axis flares. We suggest ways in which future, high-throughput X-ray observatories such as XMM and the Constellation X-ray Mission may use these reverberation signatures to probe both the mass and spin of AGN black holes, as well as the X-ray source geometry.

Gas supplied conservatively to a black hole at rates well below the Eddington rate may not be able to radiate effectively and the net energy flux, including the energy transported by the viscous torque, is likely to be close to zero at all radii. This has the consequence that the gas accretes with positive energy so that it may escape. Accordingly, we propose that only a small fraction of the gas supplied actually falls onto the black hole and that the binding energy it releases is transported radially outward by the torque so as to drive away the remainder in the form of a wind. This is a generalization of and an alternative to an "ADAF" solution. Some observational implications and possible ways to distinguish these two types of flow are briefly discussed.

We present timing analysis for a Rossi X-ray Timing Explorer observation of Cygnus X-1 in its hard/low state. This was the first RXTE observation of Cyg X-1 taken after it transited back to this state from its soft/high state. RXTE's large effective area, superior timing capabilities, and ability to obtain long, uninterrupted observations have allowed us to obtain measurements of the power spectral density (PSD), coherence function, and Fourier time lags to a decade lower in frequency and half a decade higher in frequency than typically was achieved with previous instruments. Notable aspects of our observations include a weak 0.005 Hz feature in the PSD coincident with a coherence recovery; a `hardening' of the high-frequency PSD with increasing energy; a broad frequency range measurement of the coherence function, revealing rollovers from unity coherence at both low and high frequency; and an accurate determination of the Fourier time lags over two and a half decades in frequency. As has been noted in previous similar observations, the time delay is approximately proportional to f^(-0.7), and at a fixed Fourier frequency the time delay of the hard X-rays compared to the softest energy channel tends to increase logarithmically with energy. Curiously, the 0.01-0.2 Hz coherence between the highest and lowest energy bands is actually slightly greater than the coherence between the second highest and lowest energy bands. We carefully describe all of the analysis techniques used in this paper, and we make comparisons of the data to general theoretical expectations. In a companion paper, we make specific comparisons to a Compton corona model that we have successfully used to describe the energy spectral data from this observation.

Recent work by Pringle and by Maloney, Begelman & Pringle has shown that geometrically thin, optically thick, accretion disks are unstable to warping driven by radiation torque from the central source. In this paper we generalize the study of radiation-driven warping to include general power-law surface density distributions, $\Sigma\propto R^{-\delta}$. We consider the range $\delta=3/2$ (isothermal disks) to $\delta=-3/2$, which corresponds to a radiation-pressure-supported disk; this spans the range of surface density distributions likely to be found in real astrophysical disks. There is a critical minimum size for unstable disks. The critical radius and the steady-state precession rate depend only weakly on $\delta$. The case $\delta=1$ divides the solutions into two qualitatively different regimes. Nonlinear effects must be important if the warp extends to the disk inner edge for $\delta \ge 1$, but for $\delta < 1$ nonlinearity will be important only if the warp amplitude is large at the origin. The effects of shadowing of the central source by the warp will thus be very different in the two regimes of $\delta.$ In real accretion disks the outer boundary condition is likely to be different from the zero-crossing condition that we have assumed. In accretion disks around massive black holes in active galactic nuclei, the disk will probably become optically thin before the outer disk boundary is reached, while in X-ray binaries, there will be an outer disk region (outside the circularization radius) in which the inflow velocity is zero but angular momentum is still transported. We show that in both these cases the solutions are similar to the zero-crossing eigenfunctions.

We present the results of the analysis of the broad-band spectrum of Cygnus X-1 from 3.0 to 200 keV, using data from a 10 ksec observation by the Rossi X-ray Timing Explorer. The spectrum can be well described phenomenologically by an exponentially cut-off power law with a photon index Gamma = 1.45 +/- 0.02 (a value considerably harder than typically found), e-folding energy E_fold = 162 +/- 9 keV, plus a deviation from a power law that formally can be modeled as a thermal blackbody with temperature kT_bb = 1.2 +/1 0.2 keV. Although the 3 - 30 keV portion of the spectrum can be fit with a reflected power law with Gamma = 1.81 +/- 0.01 and covering fraction f = 0.35 +/- 0.02, the quality of the fit is significantly reduced when the HEXTE data in the 30 - 100 keV range is included, as there is no observed hardening in the power law within this energy range. As a physical description of this system, we apply the accretion disc corona models of Dove, Wilms & Begelman (1997) --- where the temperature of the corona is determined self-consistently. A spherical corona with a total optical depth tau = 1.6 +/- 0.1 and an average temperature kT_c = 87 +/- 5 keV, surrounded by an exterior cold disc, does provide a good description of the data (reduced chi-squared = 1.55). These models deviate from the data by up to 7% in the 5 - 10 keV range, and we discuss possible reasons for these discrepancies. However, considering how successfully the spherical corona reproduces the 10 - 200 keV data, such ``photon-starved'' coronal geometries seem very promising for explaining the accretion processes of Cygnus X-1.

We present a simple model for an evolving radio galaxy, as it expands into the hot, X-ray emitting interstellar medium (ISM) of its host galaxy or similarly into the hot intracluster medium (ICM) of its host cluster. We solve the governing equations numerically on a grid of model parameters in order to present simple analytical tools for X-ray observations of the shocked shell that is pushed out by the cocoon and the associated cavity in the cluster emission. We apply these tools to the well known example of Perseus A to show that its time-averaged kinetic luminosity probably exceeds 10^46 ergs/sec, much larger than the estimated current power. We show how future observations can be used to extract useful source parameters such as the average kinetic power and the source age, and discuss detectability of sources at various stages of their lives.

Many of the powerful radio galaxies observed at high redshift are very small, presumably because they are very young. A simple model, which treats the radio source as a bubble expanding into a radially stratified medium, is capable of explaining the main features of their size and luminosity evolution. In particular, the model predicts strong negative luminosity evolution with increasing size, thus accounting for the large number of short-lived small sources in flux-limited samples. A variety of observational characteristics can be understood within the framework of the bubble model. I address three of these: 1) the number versus size statistics, which may provide a clue that jets in radio galaxies are intermittent, with "on-off" cycles lasting only ten to one hundred thousand yr; 2) the rarity of spectral steepening due to synchrotron cooling at frequencies of a few GHz, and its implications for source geometry and assumptions about equipartition; and 3) the prevalence of gigahertz-peaked spectra among small sources. I argue that all of these phenomena are primarily the consequence of the youth and compactness of the radio source, and depend only weakly on environment. There is no need to postulate an ISM with properties any different from that of an ordinary nearby galaxy.

There have been a number of recent spectral models that have been successful in reproducing the observed X-ray spectra of galactic black hole candidates (GBHC). However, there still exists controversy over such issues as: what are the sources of hard radiation, what is the system's geometry, is the accretion efficient or inefficient, etc. A potentially powerful tool for distinguishing among these possibilities, made possible by the Rossi X-ray Timing Explorer (RXTE), is the variability data, especially the observed phase lags and variability coherence. These data, in conjunction with spectral modeling, have the potential of determining physical sizes of the system, as well as placing strong constraints on both Compton corona and advection models. As an example, we present RXTE variability data of Cygnus X-1

The radiation-driven warping instability discovered by Pringle holds considerable promise as the mechanism responsible for producing warped, precessing accretion disks in X-ray binaries. This instability is an inherently global mode of the disk, thereby avoiding the difficulties with earlier models for the precession. Here we follow up earlier work to study the linear behavior of the instability in the specific context of a binary system. We treat the influence of the companion as an orbit-averaged, quadrupole torque on the disk. The presence of this external torque allows the existence of solutions in which the direction of precession of the warp is retrograde with respect to disk rotation, in addition to the prograde solutions which exist in the absence of external torques.

Simultaneous observations of blazars in different spectral regimes imply that the amplitude of variability depends on the observational band. Both in the low energy spectral component and in the high energy spectral component, the amplitude increases with the photon energy. We show that such behavior can be explained in terms of a two-component model, where the spectra observed during flares are superpositions of spectra arising from at least two distinct, spatially separated sources. Using the two-component model for blazar variability we derive constraints on jet physics in quasars. Particularly, we demonstrate that if production of X-rays during rapid flares is dominated by Comptonization of external radiation, then the upper limits for the total power of jets imply that the jet plasma is strongly dominated by pairs.

Jets and pulsar-fed supernova remnants (plerions) tend to develop highly organized toroidal magnetic field. Such a field structure could explain the polarization properties of some jets, and contribute to their lateral confinement. A toroidal field geometry is also central to models for the Crab Nebula - the archetypal plerion - and leads to the deduction that the Crab pulsar's wind must have a weak magnetic field. Yet this `Z-pinch' field configuration is well known to be locally unstable, even when the magnetic field is weak and/or boundary conditions slow or suppress global modes. Thus, the magnetic field structures imputed to the interiors of jets and plerions are unlikely to persist. To demonstrate this, I present a local analysis of Z-pinch instabilities for relativistic fluids in the ideal MHD limit. Kink instabilities dominate, destroying the concentric field structure and probably driving the system toward a more chaotic state in which the mean field strength is independent of radius (and in which resistive dissipation of the field may be enhanced). I estimate the timescales over which the field structure is likely to be rearranged and relate these to distances along relativistic jets and radii from the central pulsar in a plerion. I conclude that a concentric toroidal field is unlikely to exist well outside the Crab pulsar's wind termination shock. There is thus no dynamical reason to conclude that the magnetic energy flux carried by the pulsar wind is much weaker than the kinetic energy flux. Abandoning this inference would resolve a long-standing puzzle in pulsar wind theory.

We propose that the intensity changes and spectral evolution along the M87 jet can be explained by adiabatic changes to the particle momentum distribution function and the magnetic field. This is supported by the lack of any significant variation in the radio-to-optical spectral index along the jet and the moderate changes in radio brightness. Assuming a simple scaling law between magnetic field and density, we use the deprojection of a 2 cm VLA intensity map by Sparks, Biretta, & Macchetto (1996) to predict the spectral evolution along the jet. We derive limits for the magnetic field and the total pressure by comparing our results with the spatially resolved fit to spectral data by Neumann, Meisenheimer, & Roeser (1997) of a model spectrum that cuts off at approx 10^15 Hz. To explain the weakness of synchrotron cooling along the jet, the magnetic field strength must lie below the equipartition value. Although the inferred pressure in the limit of nonrelativistic bulk flow lies far above the estimated pressure of the interstellar matter in the center of M87, bulk Lorentz factors Gamma_jet in the range of 3 - 5 and inclination angles approx less than 25 deg lead to pressure estimates close to the ISM pressure. The average best fit magnetic fields we derive fall in the range of 20 - 40 microG, departing from equipartition by a factor approx 1.5 - 5. This model is consistent with the proposal by Bicknell & Begelman (1996) that the knots in the M87 jet are weak, oblique shocks. First-order Fermi acceleration will then have a minimal effect on the slope of the radio-- to-optical spectrum while possibly accounting for the X-ray spectrum.

We suggest that extragalactic radio sources are intermittent on timescales of 10^4-10^5 yr. Using a simple spherical model of a cocoon/shock system, it is found that inactive sources fade rapidly in radio luminosity but the shock in the ambient medium continues to expand supersonically, thereby keeping the whole source structure intact during the inactive phases. The fading of inactive sources, and the effect of the intermittency on the expansion velocity, can readily explain the observed over-abundance of small radio sources. In particular, the plateau in the observed distribution of sizes found by O'Dea & Baum (1997) can be interpreted as being due to intermittency. The model predicts that very young sources will be particularly radio luminous, once the effects of absorption have been accounted for. Furthermore, it predicts the existence of a significant number of faint `coasting' sources. These might be detectable in deep, low-frequency radio maps, or via the X-ray and optical emission line properties of the shock front.

The fluorescent iron Ka line is a powerful observational probe of the inner regions of black holes accretion disks. Previous studies have assumed that only material outside the radius of marginal stability can contribute to the observed line emission. Here, we show that fluorescence by material inside the radius of marginal stability, which is in the process of spiralling towards the event horizon, can have a observable influence on the iron line profile and equivalent width. For concreteness, we consider the case of a geometrically thin accretion disk, around a Schwarzschild black hole, in which fluorescence is excited by an X-ray source placed at some height above the disk and on the axis of the disk. Fully relativistic line profiles are presented for various source heights and efficiencies. It is found that the extra line flux generally emerges in the extreme red wing of the iron line, due to the large gravitational redshift experienced by photons from the region within the radius of marginal stability. We apply our models to the variable iron line seen in the ASCA spectrum of the Seyfert nucleus MCG-6-30-15. It is found that the change in the line profile, equivalent width, and continuum normalization, can be well explained as being due to a change in the height of the source above the disk. We discuss the implications of these results for distinguishing rapidly-rotating black holes from slowly rotating holes using iron line diagnostics.

We apply our self-consistent accretion disk corona (ADC) model, with two different geometries, to the broad-band X-ray spectrum of the black hole candidate Cygnus X-1. As shown in a companion paper (Dove, Wilms, and Begelman), models where the Comptonizing medium is a slab surrounding the cold accretion disk cannot have a temperature higher than about 120 keV for optical depths greater than 0.2, resulting in spectra that are much softer than the observed 10-30 keV spectrum of Cyg X-1. In addition, the slab geometry models predict a substantial ``soft excess'' at low energies, a feature not observed for Cyg X-1, and Fe K\alpha fluorescence lines that are stronger than observed. Previous Comptonization models in the literature invoke a slab geometry with the optical depth \tau_T \gta 0.3 and the coronal temperature T_c \sim 150 keV, but they are not self-consistent. Therefore, ADC models with a slab geometry are not appropriate for explaining the X-ray spectrum of Cyg X-1. Models with a spherical corona and an exterior disk, however, predict much higher self-consistent coronal temperatures than the slab geometry models. The higher coronal temperatures are due to the lower amount of reprocessing of coronal radiation in the accretion disk, giving rise to a lower Compton cooling rate. Therefore, for the sphere+disk geometry, the predicted spectrum can be hard enough to describe the observed X-ray continuum of Cyg X-1 while predicting Fe fluorescence lines having an equivalent width of \sim 40 eV. Our best-fit parameter values for the sphere+disk geometry are \tau_T \approx 1.5 and T_c \approx 90 keV.

We present the properties of accretion disk corona (ADC) models, where the radiation field, the temperature, and the total opacity of the corona are determined self-consistently. We use a non-linear Monte Carlo code to perform the calculations. As an example, we discuss models where the corona is situated above and below a cold accretion disk with a plane-parallel (slab) geometry, similar to the model of Haardt and Maraschi. By Comptonizing the soft radiation emitted by the accretion disk, the corona is responsible for producing the high-energy component of the escaping radiation. Our models include the reprocessing of radiation in the accretion disk. Here, the photons either are Compton reflected or photo-absorbed, giving rise to fluorescent line emission and thermal emission. The self-consistent coronal temperature is determined by balancing heating (due to viscous energy dissipation) with Compton cooling, determined using the fully relativistic, angle-dependent cross-sections. The total opacity is found by balancing pair productions with annihilations. We find that, for a disk temperature kT_bb \lta 200 eV, these coronae are unable to have a self-consistent temperature higher than \sim 120 keV if the total optical depth is \gta 0.2, regardless of the compactness parameter of the corona and the seed opacity. This limitation corresponds to the angle-averaged spectrum of escaping radiation having a photon index \gta 1.8 within the 5 keV - 30 keV band. Finally, all models that have reprocessing features also predict a large thermal excess at lower energies. These constraints make explaining the X-ray spectra of persistent black hole candidates with ADC models very problematic.

The Rossi X-ray Timing Explorer (RXTE) has made feasible for the first time the search for high-frequency (~ 100 Hz) periodic features in black hole candidate (BHC) systems. Such a feature, with a 67 Hz frequency, recently has been discovered in the BHC GRS 1915+105 (Morgan, Remillard, & Greiner). This feature is weak (rms variability ~0.3%-1.6%), stable in frequency (to within ~2 Hz) despite appreciable luminosity fluctuations, and narrow (quality factor Q ~ 20). Several of these properties are what one expects for a ``diskoseismic'' g-mode in an accretion disk about a 10.6 M_sun (nonrotating) - 36.3 M_sun (maximally rotating) black hole (if we are observing the fundamental mode frequency). We explore this possibility by considering the expected luminosity modulation, as well as possible excitation and growth mechanisms---including turbulent excitation, damping, and ``negative'' radiation damping. We conclude that a diskoseismic interpretation of the observations is viable.

A geometrically thin, optically thick, warped accretion disk with a central source of luminosity is subject to non-axisymmetric forces due to radiation pressure; the resulting torque acts to modify the warp. In a recent paper, \cite{pri96} used a local analysis to show that initially planar accretion disks are unstable to warping driven by radiation torque. Here we extend this work with a global analysis of the stable and unstable modes. We confirm Pringle's conclusion that thin centrally-illuminated accretion disks are generically unstable to warping via this mechanism; we discuss the time-evolution and likely steady-state of such systems and show specifically that this mechanism can explain the warping of the disk of water masers in NGC 4258 and the 164-day precession period of the accretion disk in SS 433. Radiation-driven warping and precession provides a robust mechanism for producing warped, precessing accretion disks in active galactic nuclei and X-ray binary systems.

A geometrically thin, optically thick, warped accretion disk with a central source of luminosity is subject to non-axisymmetric forces due to radiation pressure; the resulting torque acts to modify the warp. Initially planar accretion disks are unstable to warping driven by radiation torque, as shown in a local analysis by Pringle (1996) and a global analysis of the stable and unstable modes by Maloney, Begelman, & Pringle (1996). In general, the warp also precesses. We discuss the nature of this instability, and its possible implications for accretion disks in X-ray binaries and active galactic nuclei. Specifically, we argue that this effect provides a plausible explanation for the misalignment and precession of the accretion disks in X-ray binaries such as SS 433 and Her X-1; the same mechanism explains why the maser disk in NGC 4258 is warped.

We show that the sharp cutoff in the hard X-ray spectrum of NGC 4151, unusual for Seyfert 1 galaxies, can be reconciled with the average Seyfert 1 spectrum if we assume that the central source is completely hidden from our line of sight by the thick part of the accretion disk or by the broad emission line clouds. The observed X-ray radiation is produced by scattering of the Seyfert 1-type spectrum in the higher, cooler parts of the accretion disk corona, or in a wind. A sharp cutoff appears as a result of the Compton recoil effect. This model naturally explains a discrepancy regarding the inclination of the central source, inferred to be low (face-on) from observations of the iron $K\alpha$ emission line, but inferred to be high on the basis of optical and UV observations.

For the first time, detailed radiative transfer calculations of Comptonized X-ray and gamma-ray radiation in a hot pair plasma above a cold accretion disk are performed using two independent codes and methods. The simulations include both energy and pair balance as well as reprocessing of the X- and gamma-rays by the cold disk. We study both plane-parallel coronae as well as active dissipation regions having shapes of hemispheres and pill boxes located on the disk surface. It is shown, contrary to earlier claims, that plane-parallel coronae in pair balance have difficulties in selfconsistently reproducing the ranges of 2-20 keV spectral slopes, high energy cutoffs, and compactnesses inferred from observations of type 1 Seyfert galaxies. Instead, the observations are consistent with the X-rays coming from a number of individual active regions located on the surface of the disk. A number of effects such as anisotropic Compton scattering, the reflection hump, feedback to the soft photon source by reprocessing, and an active region in pair equilibrium all conspire to produce the observed ranges of X-ray slopes, high energy cutoffs, and compactnesses. The spread in spectral X-ray slopes can be due to a spread in the properties of the active regions such as their compactnesses and their elevations above the disk surface. Simplified models invoking isotropic Comptonization in spherical clouds are no longer sufficient when interpreting the data.

If gamma-ray bursts originate in our Galaxy, they probably involve violent disturbances in the magnetospheres of neutron stars. Any event of this kind is likely to trigger the sudden expulsion of magnetic flux and plasma at relativistic speed. We show that such ejecta would be braked by the interstellar medium (ISM), and that a gamma-ray flash with duration $\sim 0.1 - 100 $ s may result from this interaction. We suggest that the stronger,``smooth"class of bursts, recently identified from analyses of BATSE data, could result from the deceleration of such a blast wave by the ISM. The radiative efficiency,and hence the detectability, of such bursts would depend on the density of the circumstellar ISM. Therefore, even if the neutron star ``triggers" were uniformly distributed in space (at least within 1--2 kpc of the Sun), the bserved locations of bursts would correlate with regions of above-average ISM density.