Accretion disk dynamo as the trigger for X-ray binary state transitions

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.

Viscous boundary layers of radiation-dominated, relativistic jets. II. The free-streaming jet model

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.

Viscous boundary layers of radiation-dominated, relativistic jets. I. The two-stream model

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.

Entropy Production in Relativistic Jet Boundary Layers

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.

Accreting Black Holes

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.

Reconciling models of luminous blazars with magnetic fluxes determined by radio core shift measurements

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.

The general relativistic equations of radiation hydrodynamics in the viscous limit

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.

Constraining the Location of Gamma-Ray Flares in Luminous Blazars

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.

Recollimation Boundary Layers in Relativistic Jets

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.

Delayed X-ray emission from fallback in compact-object mergers

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.