Long-term observations have shown that black hole X-ray binaries exhibit strong, aperiodic variability on time-scales of a few milliseconds to seconds. The observed light curves display various characteristic features like a log-normal distribution of flux and a linear rms-flux relation, which indicate that the underlying variability process is stochastic in nature. It is also thought to be intrinsic to accretion. This variability has been modelled as inward propagating fluctuations of mass accretion rate, although the physical process driving the fluctuations remains puzzling. In this work, we analyse five exceptionally long duration general relativistic magnetohydrodynamic (GRMHD) simulations of optically thin, geometrically thick, black hole accretion flows to look for hints of propagating fluctuations in the simulation data. We find that the accretion profiles from these simulations do show evidence for inward propagating fluctuations below the viscous frequency by featuring strong radial coherence and positive time lags when comparing smaller to larger radii, although these time lags are generally shorter than the viscous time-scale and frequency independent. Our simulations also support the notion that the fluctuations in $\dot{M}$ build up in a multiplicative manner, as the simulations exhibit linear rms-mass flux relations, as well as log-normal distributions of their mass fluxes. When combining the mass fluxes from the simulations with an assumed emissivity profile, we additionally find broad agreement with observed power spectra and time lags, including a recovery of the frequency dependency of the time lags.

Realistic synthetic observations of theoretical source models are essential for our understanding of real observational data. In using synthetic data, one can verify the extent to which source parameters can be recovered and evaluate how various data corruption effects can be calibrated. These studies are important when proposing observations of new sources, in the characterization of the capabilities of new or upgraded instruments, and when verifying model-based theoretical predictions in a comparison with observational data. We present the SYnthetic Measurement creator for long Baseline Arrays (SYMBA), a novel synthetic data generation pipeline for Very Long Baseline Interferometry (VLBI) observations. SYMBA takes into account several realistic atmospheric, instrumental, and calibration effects. We used SYMBA to create synthetic observations for the Event Horizon Telescope (EHT), a mm VLBI array, which has recently captured the first image of a black hole shadow. After testing SYMBA with simple source and corruption models, we study the importance of including all corruption and calibration effects. Based on two example general relativistic magnetohydrodynamics (GRMHD) model images of M87, we performed case studies to assess the attainable image quality with the current and future EHT array for different weather conditions. The results show that the effects of atmospheric and instrumental corruptions on the measured visibilities are significant. Despite these effects, we demonstrate how the overall structure of the input models can be recovered robustly after performing calibration steps. With the planned addition of new stations to the EHT array, images could be reconstructed with higher angular resolution and dynamic range. In our case study, these improvements allowed for a distinction between a thermal and a non-thermal GRMHD model based on salient features in reconstructed images.

The Event Horizon Telescope image of the supermassive black hole in the galaxy M87 is dominated by a bright, unresolved ring. General relativity predicts that embedded within this image lies a thin "photon ring," which is composed of an infinite sequence of self-similar subrings that are indexed by the number of photon orbits around the black hole. The subrings approach the edge of the black hole "shadow," becoming exponentially narrower but weaker with increasing orbit number, with seemingly negligible contributions from high order subrings. Here, we show that these subrings produce strong and universal signatures on long interferometric baselines. These signatures offer the possibility of precise measurements of black hole mass and spin, as well as tests of general relativity, using only a sparse interferometric array.

We explore a simple spherical model of optically thin accretion on a Schwarzschild black hole, and study the properties of the image as seen by a distant observer. We show that a dark circular region in the center --- a shadow --- is always present. The outer edge of the shadow is located at the photon ring radius $b_{\rm ph} \equiv \sqrt{27}r_g$, where $r_g=GM/c^2$ is the gravitational radius of the accreting mass $M$. The location of the shadow edge is independent of the inner radius at which the accreting gas stops radiating. The size of the observed shadow is thus a signature of the spacetime geometry and it is hardly influenced by accretion details. We briefly discuss the relevance of these results for the Event Horizon Telescope image of the supermassive black hole in M87.

We report on a new class of solutions of black hole accretion disks that we have found through three-dimensional, global, radiative magnetohydrodynamic simulations in general relativity. It combines features of the canonical thin, slim and thick disk models but differs in crucial respects from each of them. We expect these new solutions to provide a more realistic description of black hole disks than the slim disk model. We are presenting a disk solution for a non-spinning black hole at a sub-Eddington mass accretion rate, $\dot M=0.6\,\dot M_{\rm Edd}$. By the density scale-height measure the disk appears to be thin, having a high density core near the equatorial plane of height $h_{\rho} \sim 0.1 \,r$, but most of the inflow occurs through a highly advective, turbulent, optically thick, Keplerian region that sandwiches the core and has a substantial geometrical thickness comparable to the radius, $H \sim r$. The accreting fluid is supported above the midplane in large part by the magnetic field, with the gas and radiation to magnetic pressure ratio $\beta \sim 1$, this makes the disk thermally stable, even though the radiation pressure strongly dominates over gas pressure. A significant part of the radiation emerging from the disk is captured by the black hole, so the disk is less luminous than a thin disk would be at the same accretion rate.

Recent developments in compact object astrophysics, especially the discovery of merging neutron stars by LIGO, the imaging of the black hole in M87 by the Event Horizon Telescope (EHT) and high precision astrometry of the Galactic Center at close to the event horizon scale by the GRAVITY experiment motivate the development of numerical source models that solve the equations of general relativistic magnetohydrodynamics (GRMHD). Here we compare GRMHD solutions for the evolution of a magnetized accretion flow where turbulence is promoted by the magnetorotational instability from a set of nine GRMHD codes: Athena++, BHAC, Cosmos++, ECHO, H-AMR, iharm3D, HARM-Noble, IllinoisGRMHD and KORAL. Agreement between the codes improves as resolution increases, as measured by a consistently applied, specially developed set of code performance metrics. We conclude that the community of GRMHD codes is mature, capable, and consistent on these test problems.

The joint detection of the gravitational wave (GW) GW170817 and its electromagnetic (EM) counterparts GRB170817A and kilonova AT 2017gfo has triggered extensive study of the EM emission of binary neutron star mergers. A parameter which is common to and plays a key role in both the GW and the EM analyses is the viewing angle of the binary's orbit. If a binary is viewed from different angles, the amount of GW energy changes (implying that orientation and distance are correlated) and the EM signatures can vary, depending on the structure of the emission. Information about the viewing angle of the binary orbital plane is therefore crucial to the interpretation of both the GW and the EM data, and can potentially be extracted from either side. In the first part of this study, we present a systematic analysis of how well the viewing angle of binary neutron stars can be measured from the GW data. We show that if the sky position and the redshift of the binary can be identified via the EM counterpart and an associated host galaxy, then for 50$\%$ of the systems the viewing angle can be constrained to $\leq 7^{\circ}$ uncertainty from the GW data, independent of electromagnetic emission models. On the other hand, if no redshift measurement is available, the measurement of the viewing angle with GW alone is not informative, unless the true viewing angle is close to $90^{\circ}$. This holds true even if the sky position is measured independently. Then, we consider the case where some constraints on the viewing angle can be placed from the EM data itself. We show that the EM measurements can then be used in the analysis of GW data to improve the precision of the luminosity distance, and hence of the Hubble constant, by a factor of 2 to 3.

In the absence of frequent binary collisions to isotropize the plasma, the fulfillment of the magnetohydrodynamic (MHD) Rankine-Hugoniot jump conditions by collisionless shocks is not trivial. In particular, the presence of an external magnetic field can allow for stable anisotropies, implying some departures from the isotropic MHD jumps. The functional dependence of such anisotropies in terms of the field is yet to be determined. By hypothesizing a kinetic history of the plasma through the shock front, we recently devised a theory of the downstream anisotropy, hence of the density jump, in terms of the field strength for a parallel shock [J. Plasma Phys. (2018), vol. 84, 905840604]. Here, we extend the analysis to the case of a perpendicular shock. We still find that the field reduces the density jump, but the effect is less pronounced than in the parallel case.

We present two-temperature, radiative general relativistic magnetohydrodynamic simulations of Magnetically Arrested Discs (MAD) that launch powerful relativistic jets. The mass accretion rates of our simulations are scaled to match the luminosity of the accretion flow around the supermassive black hole in M87. We consider two sub-grid prescriptions for electron heating: one based on a Landau-damped turbulent cascade, and the other based on heating from trans-relativistic magnetic reconnection. The simulations produce jets with power on the order of the observed value for M87. Both simulations produce spectra that are consistent with observations of M87 in the radio, millimetre, and submillimetre. Furthermore, the predicted image core-shifts in both models at frequencies between 15 GHz and 86 GHz are consistent with observations. At 43 and 86 GHz, both simulations produce wide opening angle jets consistent with VLBI images. Both models produce 230~GHz images with distinct black hole shadows that are resolvable by the Event Horizon Telescope (EHT), although at a viewing angle of 17 degrees, the 230 GHz images are too large to match EHT observations from 2009 and 2012. The 230 GHz images from the simulations are dynamic on time-scales of months to years, suggesting that repeated EHT observations may be able to detect the motion of rotating magnetic fields at the event horizon.

The Lunar Ultraviolet Cosmic Imager (LUCI) is a near-ultraviolet (NUV) telescope with all-spherical mirrors, designed and built to fly as a scientific payload on a lunar mission with Team Indus - the original Indian entry to the Google Lunar X-Prize. Observations from the Moon provide a unique opportunity of a stable platform with an unobstructed view of the space at all wavelengths due to the absence of atmosphere and ionosphere. LUCI is an 80 mm aperture telescope, with a field of view of 27.6'x 20.4' and a spatial resolution of 5'', will scan the sky in the NUV (200-320 nm) domain to look for transient sources. We describe here the assembly, alignment, and calibration of the complete instrument. LUCI is now in storage in a class 1000 clean room and will be delivered to our flight partner in readiness for flight.

The plasma in low-luminosity accretion flows, such as the one around the black hole at the center of M87 or Sgr A* at our Galactic Center, is expected to be collisioness and two-temperature, with protons hotter than electrons. Here, particle heating is expected to be controlled by magnetic reconnection in the transrelativistic regime $\sigma_{w}\sim 0.1$-$1$, where the magnetization $\sigma_{w}$ is the ratio of magnetic energy density to plasma enthalpy density. By means of large-scale 2D particle-in-cell simulations, we explore for a fiducial $\sigma_w=0.1$ how the dissipated magnetic energy gets partitioned between electrons and protons, as a function of $\beta_{\rm i}$ (the ratio of proton thermal pressure to magnetic pressure) and of the strength of a guide field $B_{\rm g}$ perpendicular to the reversing field $B_0$. At low $\beta_{\rm i}\;(\lesssim 0.1)$, we find that the fraction of initial magnetic energy per particle converted into electron irreversible heat is nearly independent of $B_{\rm g}/B_0$, whereas protons get heated much less with increasing $B_{\rm g}/B_0$. As a result, for large $B_{\rm g} /B_{0}$, electrons receive the overwhelming majority of irreversible particle heating (${\sim}93\%$ for $B_{\rm g} /B_{0}=6$). This is significantly different than the antiparallel case $B_{\rm g}/B_0=0$, in which electron irreversible heating accounts for only ${\sim}18\%$ of the total particle heating. At $\beta_{\rm i} \sim 2$, when both species start already relativistically hot (for our fiducial $\sigma_w=0.1$), electrons and protons each receive ${\sim}50\%$ of the irreversible particle heating, regardless of the guide field strength. Our results provide important insights into the plasma physics of electron and proton heating in hot accretion flows around supermassive black holes.

Recent polarization measurements of fast radio bursts (FRBs) provide new insights on these enigmatic sources. We show that the nearly 100% linear polarization and small variation of the polarization position angles (PAs) of multiple bursts from the same source suggest that the radiation is produced near the surface of a strongly magnetized neutron star. As the emitted radiation travels through the magnetosphere, the electric vector of the X-mode wave adiabatically rotates and stays perpendicular to the local magnetic field direction. The PA freezes at a radius where the plasma density becomes too small to be able to turn the electric vector. At the freeze-out radius, the electric field is perpendicular to the magnetic dipole moment of the neutron star projected in the plane of the sky, independent of the radiation mechanism or the orientation of the magnetic field in the emission region. We discuss a number of predictions of the model. The variation of PAs from repeating FRBs should follow the rotational period of the underlying neutron star (but the burst occurrence may not be periodic). Measuring this period will provide crucial support for the neutron star nature of the progenitors of FRBs. For FRB 121102, the small range of PA variation means that the magnetic inclination angle is less than about 20 degrees and that the observer's line of sight is outside the magnetic inclination cone. Other repeating FRBs may have a different range of PA variation from that of FRB 121102, depending on the magnetic inclination and the observer's viewing angle.

We report results from general relativistic radiation MHD (GRRMHD) simulations of a super-Eddington black hole (BH) accretion disk formed as a result of a tidal disruption event (TDE). We consider the fiducial case of a solar mass star on a mildly penetrating orbit disrupted by a supermassive BH of mass $10^6 \, M_\odot$, and consider the epoch of peak fall back rate. We post-process the simulation data to compute viewing angle dependent spectra. We perform a parameter study of the dynamics of the accretion disk as a function of BH spin and magnetic flux, and compute model spectra as a function of the viewing angle of the observer. We also consider detection limits based on the model spectra. We find that an accretion disk with a relatively weak magnetic field around the BH (so-called SANE regime of accretion) does not launch a relativistic jet, whether or not the BH is rotating. Such models reasonably reproduce several observational properties of non-jetted TDEs. The same is also true for a non-rotating BH with a strong magnetic field (MAD regime). One of our simulations has a rapidly rotating BH (spin parameter 0.9) as well as a MAD accretion disk. This model launches a powerful relativistic jet, which is powered by the BH spin energy. It reproduces the high energy emission and jet structure of the jetted TDE Swift J1644+57 surprisingly well. Jetted TDEs may thus correspond to the subset of TDE systems that have both a rapidly spinning BH and MAD accretion.

The Event Horizon Telescope (EHT), a global 230 GHz VLBI array, achieves an angular resolution of ${\approx} 20\,\mu{\rm as}$, sufficient to resolve the supermassive black hole Sagittarius A* (Sgr A*). This resolution may soon enable measurement of the black hole "shadow" size and asymmetry, predicted to be ${\approx}50\,\mu$as and ${\lesssim} 3\,\mu$as, respectively. Measurements that depart from these values could indicate a violation of the "no-hair theorem." However, refractive scattering by the turbulent ionized interstellar medium distorts the image of Sgr A*, affecting its apparent size and asymmetry. In this paper, we present a general analytic approach to quantify the expected image wander, distortion, and asymmetry from refractive scattering. If the turbulence in the scattering material of Sgr A* is close to Kolmogorov, we estimate the mean refractive image wander, distortion, and asymmetry to be 0.53 $\mu$as, 0.72 $\mu$as, and 0.52 $\mu$as at 230 GHz. However, alternative scattering models with flatter power spectra can yield larger values, up to 2.1 $\mu$as, 6.3 $\mu$as, and 5.0 $\mu$as, respectively. We demonstrate that these effects can be reduced by averaging images over multiple observations. For a small number of observations, the effects of scattering can be comparable to or greater than those from black hole spin, and they determine a fundamental limit for testing general relativity via images of Sgr A*.

We compare shadows cast by Schwarzschild black holes with those produced by two classes of naked singularities that result from gravitational collapse of spherically symmetric matter. The latter models consist of an interior naked singularity spacetime restricted to radii $r\leq R_b$, matched to Schwarzschild spacetime outside the boundary radius $R_b$. While a black hole always has a photon sphere and always casts a shadow, we find that the naked singularity models have photon spheres only if a certain parameter $M_0$ that characterizes these models satisfies $M_0\geq 2/3$, or equivalently, if $R_b\leq 3M$, where $M$ is the total mass of the object. Such models do produce shadows. However, models with $M_0<2/3$ (or $R_b>3M$) have no photon sphere and do not produce a shadow. Instead, they produce an interesting `full-moon' image. These results imply that the presence of a shadow does not by itself prove that a compact object is necessarily a black hole. The object could be a naked singularity with $M_0\geq 2/3$, and we will need other observational clues to distinguish the two possibilities. On the other hand, the presence of a full-moon image would certainly rule out a black hole and might suggest a naked singularity with $M_0<2/3$. It would be worthwhile to generalize the present study, which is restricted to spherically symmetric models, to rotating black holes and naked singularities.

Collisionless shocks follow the Rankine-Hugoniot jump conditions to a good approximation. However, for a shock propagating parallel to a magnetic field, magnetohydrodynamics states that the shock properties are independent of the field strength, whereas recent Particle-in-Cell simulations reveal a significant departure from magnetohydrodynamics behavior for such shocks in the collisionless regime. This departure is found to be caused by a field-driven anisotropy in the downstream pressure, but the functional dependence of this anisotropy on the field strength is yet to be determined. Here, we present a non-relativistic model of the plasma evolution through the shock front, allowing for a derivation of the downstream anisotropy in terms of the field strength. Our scenario assumes double adiabatic evolution of a pair plasma through the shock front. As a result, the perpendicular temperature is conserved. If the resulting downstream is firehose stable, then the plasma remains in this state. If unstable, it migrates towards the firehose stability threshold. In both cases, the conservation equations, together with the relevant hypothesis made on the temperature, allows a full determination of the downstream anisotropy in terms of the field strength.

Sagittarius A* (Sgr A*) is the variable radio, near-infrared (NIR), and X-ray source associated with accretion onto the Galactic center black hole. We present an analysis of the most comprehensive NIR variability dataset of Sgr A* to date: eight 24-hour epochs of continuous monitoring of Sgr A* at 4.5 $\mu$m with the IRAC instrument on the Spitzer Space Telescope, 93 epochs of 2.18 $\mu$m data from Naos Conica at the Very Large Telescope, and 30 epochs of 2.12 $\mu$m data from the NIRC2 camera at the Keck Observatory, in total 94,929 measurements. A new approximate Bayesian computation method for fitting the first-order structure function extracts information beyond current Fast Fourier Transformation (FFT) methods of power spectral density (PSD) estimation. With a combined fit of the data of all three observatories, the characteristic coherence timescale of Sgr A* is $\tau_{b} = 243^{+82}_{-57}$ minutes ($90\%$ credible interval). The PSD has no detectable features on timescales down to 8.5 minutes ($95\%$ credible level), which is the ISCO orbital frequency for a dimensionless spin parameter $a = 0.92$. One light curve measured simultaneously at 2.12 and 4.5 $\mu$m during a low flux-density phase gave a spectral index $\alpha_s = 1.6 \pm 0.1$ ($F_\nu \propto \nu^{-\alpha_s}$). This value implies that the Sgr A* NIR color becomes bluer during higher flux-density phases. The probability densities of flux densities of the combined datasets are best fit by log-normal distributions. Based on these distributions, the Sgr A* spectral energy distribution is consistent with synchrotron radiation from a non-thermal electron population from below 20 GHz through the NIR.

Radio images of the Galactic Center supermassive black hole, Sagittarius A* (Sgr A*), are dominated by interstellar scattering. Previous studies of Sgr A* have adopted an anisotropic Gaussian model for both the intrinsic source and the scattering, and they have extrapolated the scattering using a purely $\lambda^2$ scaling to estimate intrinsic properties. However, physically motivated source and scattering models break all three of these assumptions. They also predict that refractive scattering effects will be significant, which have been ignored in standard model fitting procedures. We analyze radio observations of Sgr A* using a physically motivated scattering model, and we develop a prescription to incorporate refractive scattering uncertainties when model fitting. We show that an anisotropic Gaussian scattering kernel is an excellent approximation for Sgr A* at wavelengths longer than 1cm, with an angular size of $(1.380 \pm 0.013) \lambda_{\rm cm}^2\,{\rm mas}$ along the major axis, $(0.703 \pm 0.013) \lambda_{\rm cm}^2\,{\rm mas}$ along the minor axis, and a position angle of $81.9^\circ \pm 0.2^\circ$. We estimate that the turbulent dissipation scale is at least $600\,{\rm km}$, with tentative support for $r_{\rm in} = 800 \pm 200\,{\rm km}$, suggesting that the ion Larmor radius defines the dissipation scale. We find that the power-law index for density fluctuations in the scattering material is $\beta < 3.47$, shallower than expected for a Kolmogorov spectrum ($\beta=11/3$), and we estimate $\beta = 3.38^{+0.08}_{-0.04}$ in the case of $r_{\rm in} = 800\,{\rm km}$. We find that the intrinsic structure of Sgr A* is nearly isotropic over wavelengths from 1.3mm to 1.3cm, with a size that is roughly proportional to wavelength. We discuss implications for models of Sgr A*, for theories of interstellar turbulence, and for imaging Sgr A* with the Event Horizon Telescope.

The accretion flow around the Galactic Center black hole Sagittarius A* (Sgr A*) is expected to have an electron temperature that is distinct from the ion temperature, due to weak Coulomb coupling in the low-density plasma. We present four two-temperature general relativistic radiative magnetohydrodynamic (GRRMHD) simulations of Sgr A* performed with the code KORAL. These simulations use different electron heating prescriptions, motivated by different models of the underlying plasma microphysics. We compare the Landau-damped turbulent cascade model used in previous work with a new prescription we introduce based on the results of particle-in-cell simulations of magnetic reconnection. With the turbulent heating model, electrons are preferentially heated in the polar outflow, whereas with the reconnection model electrons are heated by nearly the same fraction everywhere in the accretion flow. The spectra of the two models are similar around the submillimetre synchrotron peak, but the models heated by magnetic reconnection produce variability more consistent with the level observed from Sgr A*. All models produce 230~GHz images with distinct black hole shadows which are consistent with the image size measured by the Event Horizon Telescope, but only the turbulent heating produces an anisotropic `disc-jet' structure where the image is dominated by a polar outflow or jet at frequencies below the synchrotron peak. None of our models can reproduce the observed radio spectral slope, the large near-infrared and X-ray flares, or the near-infrared spectral index, all of which suggest non-thermal electrons are needed to fully explain the emission from Sgr A*.

Scattering in the ionized interstellar medium is commonly observed to be anisotropic, with theories of magnetohydrodynamic (MHD) turbulence explaining the anisotropy through a preferred magnetic field direction throughout the scattering regions. In particular, the line of sight to the Galactic Center supermassive black hole, Sgr A*, exhibits strong and anisotropic scattering, which dominates its observed size at wavelengths of a few millimeters and longer. Therefore, inferences of the intrinsic structure of \sgra\ at these wavelengths are sensitive to the assumed scattering model. In addition, extrapolations of the scattering model from long wavelengths, at which its parameters are usually estimated, to 1.3 mm, where the Event Horizon Telescope (EHT) seeks to image Sgr A* on Schwarzschild-radius scales, are also sensitive to the assumed scattering model. Past studies of Sgr A* have relied on simple Gaussian models for the scattering kernel that effectively presume an inner scale of turbulence far greater than the diffractive scale; this assumption is likely violated for Sgr A* at 1.3 mm. We develop a physically motivated model for anisotropic scattering, using a simplified model for MHD turbulence with a finite inner scale and a wandering transverse magnetic field direction. We explore several explicit analytic models for this wandering and derive the expected observational properties --- scatter broadening and refractive scintillation --- for each. For expected values of the inner scale, the scattering kernel for all models is markedly non-Gaussian at 1.3 mm but is straightforward to calculate and depends only weakly on the assumed model for the wandering of the magnetic field direction. On the other hand, in all models, the refractive substructure depends strongly on the wandering model and may be an important consideration in imaging Sgr A* with the EHT.

Interferometric imaging now achieves angular resolutions as fine as 10 microarcsec, probing scales that are inaccessible to single telescopes. Traditional synthesis imaging methods require calibrated visibilities; however, interferometric calibration is challenging, especially at high frequencies. Nevertheless, most studies present only a single image of their data after a process of "self-calibration," an iterative procedure where the initial image and calibration assumptions can significantly influence the final image. We present a method for efficient interferometric imaging directly using only closure amplitudes and closure phases, which are immune to station-based calibration errors. Closure-only imaging provides results that are as non-committal as possible and allows for reconstructing an image independently from separate amplitude and phase self-calibration. While closure-only imaging eliminates some image information (e.g., the total image flux density and the image centroid), this information can be recovered through a small number of additional constraints. We demonstrate that closure-only imaging can produce high fidelity results, even for sparse arrays such as the Event Horizon Telescope, and that the resulting images are independent of the level of systematic amplitude error. We apply closure imaging to VLBA and ALMA data and show that it is capable of matching or exceeding the performance of traditional self-calibration and CLEAN for these data sets.

Recent X-ray observations of merger shocks in galaxy clusters have shown that the post-shock plasma is two-temperature, with the protons being hotter than the electrons. In this work, the second of a series, we investigate by means of two-dimensional particle-in-cell simulations the efficiency of electron irreversible heating in perpendicular low Mach number shocks. We consider values of plasma beta (ratio of thermal and magnetic pressures) in the range $4\lesssim \beta_{p0}\lesssim 32$ and sonic Mach number (ratio of shock speed to pre-shock sound speed) in the range $2\lesssim M_{s}\lesssim 5$, as appropriate for galaxy cluster shocks. As shown in Paper I, magnetic field amplification - induced by shock compression of the pre-shock field, or by strong proton cyclotron and mirror modes accompanying the relaxation of proton temperature anisotropy - can drive the electron temperature anisotropy beyond the threshold of the electron whistler instability. The growth of whistler waves breaks the electron adiabatic invariance, and allows for efficient entropy production. We find that the post-shock electron temperature $T_{e2}$ exceeds the adiabatic expectation $T_{e2,\rm ad}$ by an amount $(T_{e2}-T_{e2,\rm ad})/T_{e0}\simeq 0.044 \,M_s (M_s-1)$ (here, $T_{e0}$ is the pre-shock temperature), which depends only weakly on the plasma beta, over the range $4\lesssim \beta_{p0}\lesssim 32$ which we have explored, and on the proton-to-electron mass ratio (the coefficient of $\simeq 0.044$ is measured for our fiducial $m_i/m_e=49$, and we estimate that it will decrease to $\simeq 0.03$ for the realistic mass ratio). Our results have important implications for current and future observations of galaxy cluster shocks in the radio band (synchrotron emission and Sunyaev-Zel'dovich effect) and at X-ray frequencies.

Recent X-ray observations of merger shocks in galaxy clusters have shown that the post-shock plasma is two-temperature, with the protons hotter than the electrons. By means of two-dimensional particle-in-cell simulations, we study the physics of electron irreversible heating in perpendicular low Mach number shocks, for a representative case with sonic Mach number of 3 and plasma beta of 16. We find that two basic ingredients are needed for electron entropy production: (i) an electron temperature anisotropy, induced by field amplification coupled to adiabatic invariance; and (ii) a mechanism to break the electron adiabatic invariance itself. In shocks, field amplification occurs at two major sites: at the shock ramp, where density compression leads to an increase of the frozen-in field; and farther downstream, where the shock-driven proton temperature anisotropy generates strong proton cyclotron and mirror modes. The electron temperature anisotropy induced by field amplification exceeds the threshold of the electron whistler instability. The growth of whistler waves breaks the electron adiabatic invariance, and allows for efficient entropy production. We find that the electron heating efficiency displays only a weak dependence on mass ratio (less than 30 percent drop, as we increase the mass ratio from 49 up to 1600). We develop an analytical model of electron irreversible heating and show that it is in excellent agreement with our simulation results.

Hot collisionless accretion flows, such as the one in Sgr A$^{*}$ at our Galactic center, provide a unique setting for the investigation of magnetic reconnection. Here, protons are non-relativistic while electrons can be ultra-relativistic. By means of two-dimensional particle-in-cell simulations, we investigate electron and proton heating in the outflows of trans-relativistic reconnection (i.e., $\sigma_w\sim 0.1-1$, where the magnetization $\sigma_w$ is the ratio of magnetic energy density to enthalpy density). For both electrons and protons, we find that heating at high $\beta_{\rm i}$ (here, $\beta_{\rm i}$ is the ratio of proton thermal pressure to magnetic pressure) is dominated by adiabatic compression ('adiabatic heating'), while at low $\beta_{\rm i}$ it is accompanied by a genuine increase in entropy ('irreversible heating'). For our fiducial $\sigma_w=0.1$, the irreversible heating efficiency at $\beta_{\rm i}\lesssim 1$ is nearly independent of the electron-to-proton temperature ratio $T_{\rm e}/T_{\rm i}$ (which we vary from $0.1$ up to $1$), and it asymptotes to $\sim 2\%$ of the inflowing magnetic energy in the low-$\beta_{\rm i}$ limit. Protons are heated more efficiently than electrons at low and moderate $\beta_{\rm i}$ (by a factor of $\sim7$), whereas the electron and proton heating efficiencies become comparable at $\beta_{\rm i}\sim 2$ if $T_{\rm e}/T_{\rm i}=1$, when both species start already relativistically hot. We find comparable heating efficiencies between the two species also in the limit of relativistic reconnection ($\sigma_w\gtrsim 1$). Our results have important implications for the two-temperature nature of collisionless accretion flows, and may provide the sub-grid physics needed in general relativistic MHD simulations.

We use VLBA+VLA observations to measure the sizes of the scatter-broadened images of 6 of the most heavily scattered known pulsars: 3 within the Galactic Centre (GC) and 3 elsewhere in the inner Galactic plane. By combining the measured sizes with temporal pulse broadening data from the literature and using the thin-screen approximation, we locate the scattering medium along the line of sight to these 6 pulsars. At least two scattering screens are needed to explain the observations of the GC sample. We show that the screen inferred by previous observations of SGR J1745-2900 and Sgr A*, which must be located far from the GC, falls off in strength on scales < 0.2 degree. A second scattering component closer to (< 2 kpc) or even (tentatively) within (< 700 pc) the GC produces most or all of the temporal broadening observed in the other GC pulsars. Outside the GC, the scattering locations for all three pulsars are ~2 kpc from Earth, consistent with the distance of the Carina-Sagittarius or Scutum spiral arm. For each object the 3D scattering origin coincides with a known HII region (and in one case also a supernova remnant), suggesting that such objects preferentially cause the intense interstellar scattering seen towards the Galactic plane. We show that the HII regions should contribute > 25% of the total dispersion measure (DM) towards these pulsars, and calculate reduced DM distances. Those distances for other pulsars lying behind HII regions may be similarly overestimated.

Current simulations of hot accretion flows around black holes assume either a single-temperature gas or, at best, a two-temperature gas with thermal ions and electrons. However, processes like magnetic reconnection and shocks can accelerate electrons into a nonthermal distribution, which will not quickly thermalise at the very low densities found in many systems. Such nonthermal electrons have been invoked to explain the infrared and X-ray spectra and strong variability of Sagittarius A* (Sgr A*), the black hole at the Galactic Center. We present a method for self-consistent evolution of a nonthermal electron population in the GRMHD code KORAL. The electron distribution is tracked across Lorentz factor space and is evolved in space and time, in parallel with thermal electrons, thermal ions, and radiation. In the present study, for simplicity, energy injection into the nonthermal distribution is taken as a fixed fraction of the local electron viscous heating rate. Numerical results are presented for a model with a low mass accretion rate similar to that of Sgr A*. We find that the presence of a nonthermal population of electrons has negligible effect on the overall dynamics of the system. Due to our simple uniform particle injection prescription, the radiative power in the nonthermal simulation is enhanced at large radii. The energy distribution of the nonthermal electrons shows a synchrotron cooling break, with the break Lorentz factor varying with location and time, reflecting the complex interplay between the local viscous heating rate, magnetic field strength, and fluid velocity.

We implement the KORAL code to perform two sets of very long general relativistic radiation magnetohydrodynamic simulations of an axisymmetric optically thin magnetized flow around a non-rotating black hole: one with a new term in the electromagnetic field tensor due to the radiation pressure felt by the plasma electrons on the comoving frame of the electron-proton plasma, and one without. The source of the radiation is the accretion flow itself. Without the new term, the system evolves to a standard accretion flow due to the development of the magneto-rotational instability (MRI). With the new term, however, the system eventually evolves to a magnetically arrested state (MAD) in which a large scale jet-like magnetic field threads the black hole horizon. Our results confirm the secular action of the Cosmic Battery in accretion flows around astrophysical black holes.

According to magnetohydrodynamics (MHD), the encounter of two collisional magnetized plasmas at high velocity gives rise to shock waves. Investigations conducted so far have found that the same conclusion still holds in the case of collisionless plasmas. For the case of a flow-aligned field, MHD stipulates that the field and the fluid are disconnected, so that the shock produced is independent of the field. We present a violation of this MHD prediction when considering the encounter of two cold pair plasmas along a flow-aligned magnetic field. As the guiding magnetic field grows, isotropization is progressively suppressed, resulting in a strong influence of the field on the resulting structure. A micro-physics analysis allows to understand the mechanisms at work. Particle-in-cell simulations also support our conclusions and show that the results are not restricted to a strictly parallel field.

We have performed a set of simulations of expanding, spherically symmetric nebulae inflated by winds from accreting black holes in ultraluminous X-ray sources (ULXs). We implemented a realistic cooling function to account for free-free and bound-free cooling. For all model parameters we considered, the forward shock in the interstellar medium becomes radiative at a radius $\sim $ 100 pc. The emission is primarily in the optical and UV, and the radiative luminosity is about 50% of the total kinetic luminosity of the wind. In contrast, the reverse shock in the wind is adiabatic so long as the terminal outflow velocity of the wind $v_{\rm w} \sim 0.003c$. The shocked wind in these models radiates in X-rays, but with a luminosity of only $\sim 10^{35} \rm\,erg\,s^{-1}$. For wind velocities $v_{\rm w} \sim 0.001c$, the shocked wind becomes radiative, but it is no longer hot enough to produce X-rays. Instead it emits in optical and UV, and the radiative luminosity is comparable to 100% of the wind kinetic luminosity. We suggest that measuring the optical luminosities and putting limits on the X-ray and radio emission from shock-ionized ULX bubbles may help in estimating the mass outflow rate of the central accretion disk and the velocity of the outflow.

Many black hole (BH) candidates have been discovered in X-ray binaries and in the nuclei of galaxies. The prediction of Einstein's general relativity is that BHs have an event horizon --- a one-way membrane through which particles fall into the BH but cannot exit. However, except for the very few nearby supermassive BH candidates, our telescopes are unable to resolve and provide a direct proof of the event horizon. Here, we propose a novel observation that supports the existence of event horizons around supermassive BH candidates heavier than 10^{7.5} solar masses. Instead of an event horizon, if the BH candidate has a hard surface, when a star falls onto the surface, the shocked baryonic gas will form a radiation pressure supported envelope that shines at the Eddington luminosity for an extended period of time from months to years. We show that such emission has already been ruled out by the Pan-STARRS1 3pi survey if supermassive BH candidates have a hard surface at radius larger than 1 + 10^{-4.4} times the Schwarzschild radius. Future observations by LSST should be able to improve the limit to 1 + 10^{-6}.

We present general relativistic radiation MHD simulations of super-Eddington accretion on a $10M_\odot$ black hole. We consider a range of mass accretion rates, black hole spins, and magnetic field configurations. We compute the spectra and images of the models as a function of viewing angle, and compare them with the observed properties of ultraluminous X-ray sources (ULXs). The models easily produce apparent luminosities in excess of $10^{40}~{\rm erg\,s^{-1}}$ for pole-on observers. However, the angle-integrated radiative luminosities rarely exceed $2.5\times10^{39}~{\rm erg\,s^{-1}}$ even for mass accretion rates of tens of Eddington. The systems are thus radiatively inefficient, though they are energetically efficient when the energy output in winds and jets is also counted. The simulated models reproduce the main empirical types of spectra --- disk-like, supersoft, soft, hard --- observed in ULXs. The magnetic field configuration, whether MAD ("magnetically arrested disk") or SANE ("standard and normal evolution"), has a strong effect on the results. In SANE models, the X-ray spectral hardness is almost independent of accretion rate, but decreases steeply with increasing inclination. MAD models with non-spinning black holes produce significantly softer spectra at higher values of $\dot{M}$, even at low inclinations. MAD models with rapidly spinning black holes are quite different from all other models. They are radiatively efficient (efficiency factor $\sim10-20\%$), super-efficient when the mechanical energy output is also included ($70\%$), and produce hard blazar-like spectra. In all models, the emission shows strong geometrical beaming, which disagrees with the more isotropic illumination favored by observations of ULX bubbles.

We present a numerical method which evolves a two-temperature, magnetized, radiative, accretion flow around a black hole, within the framework of general relativistic radiation magnetohydrodynamics. As implemented in the code KORAL, the gas consists of two sub-components -- ions and electrons -- which share the same dynamics but experience independent, relativistically consistent, thermodynamical evolution. The electrons and ions are heated independently according to a standard prescription from the literature for magnetohydrodynamical turbulent dissipation. Energy exchange between the particle species via Coulomb collisions is included. In addition, electrons gain and lose energy and momentum by absorbing and emitting synchrotron and bremsstrahlung radiation, and through Compton scattering. All evolution equations are handled within a fully covariant framework in the relativistic fixed-metric spacetime of the black hole. Numerical results are presented for five models of low luminosity black hole accretion. In the case of a model with a mass accretion rate $\dot{M}\sim10^{-8} \dot M_{\rm Edd}$, we find that radiation has a negligible effect on either the dynamics or the thermodynamics of the accreting gas. In contrast, a model with a larger $\dot{M}\sim 4\times 10^{-4} \dot M_{\rm Edd}$ behaves very differently. The accreting gas is much cooler and the flow is geometrically less thick, though it is not quite a thin accretion disk.

We present an extension to the general relativistic radiation magnetohydrodynamic code HARMRAD to account for emission and absorption by thermal cyclo-synchrotron, double Compton, bremsstrahlung, low-temperature OPAL opacities as well as Thomson and Compton scattering. We approximate the radiation field as a Bose-Einstein distribution and evolve it using the radiation number-energy-momentum conservation equations in order to track photon hardening. We perform various simulations to study how these extensions affect the radiative properties of magnetically-arrested disks accreting at Eddington to super-Eddington rates. We find that double Compton dominates bremsstrahlung in the disk within a radius of $r\sim 15r_g$ (gravitational radii) at a hundred times the Eddington accretion rate, and within smaller radii at lower accretion rates. Double Compton and cyclo-synchrotron regulate radiation and gas temperatures in the corona, while cyclo-synchrotron regulates temperatures in the jet. Interestingly, as the accretion rate drops to Eddington, an optically thin corona develops whose gas temperature of $T\sim 10^9$K is $\sim 100$ times higher than the disk's black body temperature. Our results show the importance of double Compton and synchrotron in super-Eddington disks, magnetized coronae, and jets.

Images of the linear polarization of synchrotron radiation around Active Galactic Nuclei (AGN) identify their projected magnetic field lines and provide key data for understanding the physics of accretion and outflow from supermassive black holes. The highest resolution polarimetric images of AGN are produced with Very Long Baseline Interferometry (VLBI). Because VLBI incompletely samples the Fourier transform of the source image, any image reconstruction that fills in unmeasured spatial frequencies will not be unique and reconstruction algorithms are required. In this paper, we explore extensions of the Maximum Entropy Method (MEM) to linear polarimetric VLBI imaging. In contrast to previous work, our polarimetric MEM algorithm combines a Stokes I imager that uses only bispectrum measurements that are immune to atmospheric phase corruption with a joint Stokes Q and U imager that operates on robust polarimetric ratios. We demonstrate the effectiveness of our technique on 7- and 3-mm wavelength quasar observations from the VLBA and simulated 1.3-mm Event Horizon Telescope observations of Sgr A* and M87. Consistent with past studies, we find that polarimetric MEM can produce superior resolution compared to the standard CLEAN algorithm when imaging smooth and compact source distributions. As an imaging framework, MEM is highly adaptable, allowing a range of constraints on polarization structure. Polarimetric MEM is thus an attractive choice for image reconstruction with the EHT.

At super-Eddington rates accretion flows onto black holes have been described as slim (aspect ratio $H/R \lesssim 1$) or thick (H/R >1) discs, also known as tori or (Polish) doughnuts. The relation between the two descriptions has never been established, but it was commonly believed that at sufficiently high accretion rates slim discs inflate, becoming thick. We wish to establish under what conditions slim accretion flows become thick. We use analytical equations, numerical 1+1 schemes, and numerical radiative MHD codes to describe and compare various accretion flow models at very high accretion rates.We find that the dominant effect of advection at high accretion rates precludes slim discs becoming thick. At super-Eddington rates accretion flows around black holes can always be considered slim rather than thick.

Newly recognized effects of refractive scattering in the ionized interstellar medium have broad implications for very long baseline interferometry (VLBI) at extreme angular resolutions. Building upon work by Blandford & Narayan (1985), we present a simplified, geometrical optics framework, which enables rapid, semi-analytic estimates of refractive scattering effects. We show that these estimates exactly reproduce previous results based on a more rigorous statistical formulation. We then derive new expressions for the scattering-induced fluctuations of VLBI observables such as closure phase, and we demonstrate how to calculate the fluctuations for arbitrary quantities of interest using a Monte Carlo technique.

We study energy flows in geometrically thick accretion discs, both optically thick and thin, using general relativistic, three-dimensional simulations of black hole accretion flows. We find that for non-rotating black holes the efficiency of the total feedback from thick accretion discs is $3\%$ - roughly half of the thin disc efficiency. This amount of energy is ultimately distributed between outflow and radiation, the latter scaling weakly with the accretion rate for super-critical accretion rates, and returned to the interstellar medium. Accretion on to rotating black holes is more efficient because of the additional extraction of rotational energy. However, the jet component is collimated and likely to interact only weakly with the environment, whereas the outflow and radiation components cover a wide solid angle.

We construct models of static, spherically symmetric shells supported by the radiation flux of a luminous neutron star in the Schwarzschild metric. The atmospheres are disconnected from the star and levitate above its surface. Gas pressure and density inversion appear in the inner region of these atmospheres, which is a purely relativistic phenomenon. We account for the scattering opacity dependence on temperature and utilize the relativistic M1 closure scheme for the radiation tensor, hence allowing for a fully GR-consistent treatment of the photon flux and radiation tensor anisotropy. In this way we are able to address atmospheres of both large and moderate/low optical depths with the same set of equations. We discuss properties of the levitating atmospheres and find that they may indeed be optically thick, with the distance between star surface and the photosphere expanding as luminosity increases. These results may be relevant for the photosphereric radius expansion X-ray bursts.

We propose here that the well-known black hole paradoxes such as the information loss and teleological nature of the event horizon are restricted to a particular idealized case, which is the homogeneous dust collapse model. In this case, the event horizon, which defines the boundary of the black hole, forms initially, and the singularity in the interior of the black hole at a later time. We show that, in contrast, gravitational collapse from physically more realistic initial conditions typically leads to the scenario in which the event horizon and space-time singularity form simultaneously. We point out that this apparently simple modification can mitigate the causality and teleological paradoxes, and also lends support to two recently suggested solutions to the information paradox, namely, the `firewall' and `classical chaos' proposals.

Earth-space interferometry with RadioAstron provides the highest direct angular resolution ever achieved in astronomy at any wavelength. RadioAstron detections of the classic quasar 3C273 on interferometric baselines up to 171,000 km suggest brightness temperatures exceeding expected limits from the "inverse-Compton catastrophe" by two orders of magnitude. We show that at 18 cm, these estimates most probably arise from refractive substructure introduced by scattering in the interstellar medium. We use the scattering properties to estimate an intrinsic brightness temperature of 7*10^12 K, which is consistent with expected theoretical limits, but which is ~15 times lower than estimates that neglect substructure. At 6.2 cm, the substructure influences the measured values appreciably but gives an estimated brightness temperature that is comparable to models that do not account for the substructure. At 1.35 cm, the substructure does not affect the extremely high inferred brightness temperatures, in excess of 10^13 K. We also demonstrate that for a source having a Gaussian surface brightness profile, a single long-baseline estimate of refractive substructure determines an absolute minimum brightness temperature, if the scattering properties along a given line of sight are known, and that this minimum accurately approximates the apparent brightness temperature over a wide range of total flux densities.

In general relativity static gaseous atmospheres may be in hydrostatic balance in the absence of a supporting stellar surface, provided that the luminosity is close to the Eddington value. We construct analytic models of optically thin, spherically symmetric shells supported by the radiation pressure of a luminous central body in the Schwarzschild metric. Opacity is assumed to be dominated by Thomson scattering. The inner parts of the atmospheres, where the luminosity locally has supercritical values, are characterized by a density and pressure inversion. The atmospheres are convectively and Rayleigh-Taylor stable, and there is no outflow of gas.

Near a black hole, differential rotation of a magnetized accretion disk is thought to produce an instability that amplifies weak magnetic fields, driving accretion and outflow. These magnetic fields would naturally give rise to the observed synchrotron emission in galaxy cores and to the formation of relativistic jets, but no observations to date have been able to resolve the expected horizon-scale magnetic-field structure. We report interferometric observations at 1.3-millimeter wavelength that spatially resolve the linearly polarized emission from the Galactic Center supermassive black hole, Sagittarius A*. We have found evidence for partially ordered fields near the event horizon, on scales of ~6 Schwarzschild radii, and we have detected and localized the intra-hour variability associated with these fields.

This paper describes the Polarization Spectroscopic Telescope Array (PolSTAR), a mission proposed to NASA's 2014 Small Explorer (SMEX) announcement of opportunity. PolSTAR measures the linear polarization of 3-50 keV (requirement; goal: 2.5-70 keV) X-rays probing the behavior of matter, radiation and the very fabric of spacetime under the extreme conditions close to the event horizons of black holes, as well as in and around magnetars and neutron stars. The PolSTAR design is based on the technology developed for the Nuclear Spectroscopic Telescope Array (NuSTAR) mission launched in June 2012. In particular, it uses the same X-ray optics, extendable telescope boom, optical bench, and CdZnTe detectors as NuSTAR. The mission has the sensitivity to measure ~1% linear polarization fractions for X-ray sources with fluxes down to ~5 mCrab. This paper describes the PolSTAR design as well as the science drivers and the potential science return.

We describe HEROIC, an upgraded version of the relativistic radiative post-processor code HERO described in a previous paper, but which now Includes Comptonization. HEROIC models Comptonization via the Kompaneets equation, using a quadratic approximation for the source function in the short characteristics radiation solver. It employs a simple form of accelerated lambda iteration to handle regions of high scattering opacity. In addition to solving for the radiation field, HEROIC also solves for the gas temperature by applying the condition of radiative equilibrium. We present benchmarks and tests of the Comptonization module in HEROIC with simple 1D and 3D scattering problems. We also test the ability of the code to handle various relativistic effects using model atmospheres and accretion flows in a black hole space-time. We present two applications of HEROIC to general relativistic MHD simulations of accretion discs. One application is to a thin accretion disc around a black hole. We find that the gas below the photosphere in the multi-dimensional HEROIC solution is nearly isothermal, quite different from previous solutions based on 1D plane parallel atmospheres. The second application is to a geometrically thick radiation-dominated accretion disc accreting at 11 times the Eddington rate. The multi-dimensional HEROIC solution shows that, for observers who are on axis and look down the polar funnel, the isotropic equivalent luminosity could be more than ten times the Eddington limit, even though the spectrum might still look thermal and show no signs of relativistic beaming.

We present a set of four three-dimensional, general relativistic, radiation MHD simulations of black hole accretion at super-critical mass accretion rates, $\dot{M} > \dot{M}_{\rm Edd}$. We use these simulations to study how disk properties are modified when we vary the black hole mass, the black hole spin, or the mass accretion rate. In the case of a non-rotating black hole, we find that the total efficiency is of order $3\%\dot M c^2$, approximately a factor of two less than the efficiency of a standard thin accretion disk. The radiation flux in the funnel along the axis is highly super-Eddington, but only a small fraction of the energy released by accretion escapes in this region. The bulk of the $3\%\dot M c^2$ of energy emerges farther out in the disk, either in the form of photospheric emission or as a wind. In the case of a black hole with a spin parameter of 0.7, we find a larger efficiency of about $8\%\dot M c^2$. By comparing the relative importance of advective and diffusive radiation transport, we show that photon trapping is effective near the equatorial plane. However, near the disk surface, vertical transport of radiation by diffusion dominates. We compare the properties of our fiducial three-dimensional run with those of an equivalent two-dimensional axisymmetric model with a mean-field dynamo. The latter simulation runs nearly 100 times faster than the three-dimensional simulation, and gives very similar results for time-averaged properties of the accretion flow.

We introduce a new method for treating Comptonization in computational fluid dynamics. By construction, this method conserves the number of photons. Whereas the traditional "blackbody Comptonization" approach assumes that the radiation is locally a perfect blackbody and therefore uses a single parameter, the radiation temperature, to describe the radiation, the new "photon-conserving Comptonization" approach treats the photon gas as a Bose-Einstein fluid and keeps track of both the radiation temperature and the photon number density. We have implemented photon-conserving Comptonization in the general relativistic radiation magnetohydrodynamical code KORAL and we describe its impact on simulations of mildly super-critical black hole accretion disks. We find that blackbody Comptonization underestimates the gas and radiation temperature by up to a factor of two compared to photon-conserving Comptonization. This discrepancy could be serious when computing spectra. The photon-conserving simulation indicates that the spectral color correction factor of the escaping radiation in the funnel region of the disk could be as large as 5.

HERO (Hybrid Evaluator for Radiative Objects) is a 3D general relativistic radiative transfer code which has been tailored to the problem of analyzing radiation from simulations of relativistic accretion discs around black holes. HERO is designed to be used as a postprocessor. Given some fixed fluid structure for the disc (i.e. density and velocity as a function of position from a hydrodynamics or magnetohydrodynamics simulation), the code obtains a self-consistent solution for the radiation field and for the gas temperatures using the condition of radiative equilibrium. The novel aspect of HERO is that it combines two techniques: 1) a short characteristics (SC) solver that quickly converges to a self consistent disc temperature and radiation field, with 2) a long characteristics (LC) solver that provides a more accurate solution for the radiation near the photosphere and in the optically thin regions. By combining these two techniques, we gain both the computational speed of SC and the high accuracy of LC. We present tests of HERO on a range of 1D, 2D and 3D problems in flat space and show that the results agree well with both analytical and benchmark solutions. We also test the ability of the code to handle relativistic problems in curved space. Finally, we discuss the important topic of ray-defects, a major limitation of the SC method, and describe our strategy for minimizing the induced error.

We explore the variability properties of long, high cadence GRMHD simulations across the electromagnetic spectrum using an efficient, GPU-based radiative transfer algorithm. We focus on both disk- and jet-dominated simulations with parameters that successfully reproduce the time-averaged spectral properties of Sgr A* and the size of its image at 1.3mm. We find that the disk-dominated models produce short timescale variability with amplitudes and power spectra that closely resemble those inferred observationally. In contrast, jet-dominated models generate only slow variability, at lower flux levels. Neither set of models show any X-ray flares, which most likely indicate that additional physics, such as particle acceleration mechanisms, need to be incorporated into the GRMHD simulations to account for them. The disk-dominated models show strong, short-lived mm/IR flares, with short (<~ 1hr) time lags between the mm and IR wavelengths, that arise from strong-field gravitational lensing of magnetic flux tubes near the horizon. Such events provide a natural explanation for the observed IR flares with no X-ray counterparts.

The 6 billion solar mass supermassive black hole at the center of the giant elliptical galaxy M87 powers a relativistic jet. Observations at millimeter wavelengths with the Event Horizon Telescope have localized the emission from the base of this jet to angular scales comparable to the putative black hole horizon. The jet might be powered directly by an accretion disk or by electromagnetic extraction of the rotational energy of the black hole. However, even the latter mechanism requires a confining thick accretion disk to maintain the required magnetic flux near the black hole. Therefore, regardless of the jet mechanism, the observed jet power in M87 implies a certain minimum mass accretion rate. If the central compact object in M87 were not a black hole but had a surface, this accretion would result in considerable thermal near-infrared and optical emission from the surface. Current flux limits on the nucleus of M87 strongly constrain any such surface emission. This rules out the presence of a surface and thereby provides indirect evidence for an event horizon.

We describe a set of simulations of super-critical accretion onto a non-rotating supermassive BH. The accretion flow is radiation pressure dominated and takes the form of a geometrically thick disk with twin low-density funnels around the rotation axis. For accretion rates $\gtrsim 10 \dot M_{\rm Edd}$, there is sufficient gas in the funnel to make this region optically thick. Radiation from the disk first flows into the funnel, after which it accelerates the optically thick funnel gas along the axis. The resulting jet is baryon-loaded and has a terminal density-weighted velocity $\approx 0.3c$. Much of the radiative luminosity is converted into kinetic energy by the time the escaping gas becomes optically thin. For an observer viewing down the axis, the isotropic equivalent luminosity of total energy is as much as $10^{48}\,\rm erg\,s^{-1}$ for a $10^7 M_\odot$ BH accreting at $10^3$ Eddington. Therefore, energetically, the simulated jets are consistent with observations of the most powerful tidal disruption events, e.g., Swift J1644. The jet velocity is, however, too low to match the Lorentz factor $\gamma > 2$ inferred in J1644. There is no such conflict in the case of other tidal disruption events. Since favorably oriented observers see isotropic equivalent luminosities that are highly super-Eddington, the simulated models can explain observations of ultra-luminous X-ray sources, at least in terms of luminosity and energetics, without requiring intermediate mass black holes. Finally, since the simulated jets are baryon-loaded and have mildly relativistic velocities, they match well the jets observed in SS433. The latter are, however, more collimated than the simulated jets. This suggests that, even if magnetic fields are not important for acceleration, they may perhaps still play a role in confining the jet.

Previous MHD simulations have shown that wind must exist in black hole hot accretion flows. In this paper, we continue our study by investigating the detailed properties of wind, such as mass flux and poloidal speed, and the mechanism of wind production. For this aim, we make use of a three dimensional GRMHD simulation of hot accretion flows around a Schwarzschild black hole. The simulation is designed so that the magnetic flux is not accumulated significantly around the black hole. To distinguish real wind from turbulent outflows, we track the trajectories of the virtual Largrangian particles from simulation data. We find two types of real outflows, i.e., a quasi-relativistic jet close to the axis and a sub-relativistic wind subtending a much larger solid angle. Most of the wind originates from the surface layer of the accretion flow. The poloidal wind speed almost remains constant once they are produced, but the flux-weighted wind speed roughly follows $v_{\rm p, wind}(r)\approx 0.25 v_k(r)$. The mass flux of jet is much lower but the speed is much higher, $v_{\rm p,jet}\sim (0.3-0.4) c$. Consequently, both the energy and momentum fluxes of the wind are much larger than those of the jet. We find that the wind is produced and accelerated primarily by the combination of centrifugal force and magnetic pressure gradient, while the jet is mainly accelerated by magnetic pressure gradient. Finally, we find that the wind production efficiency $\epsilon_{\rm wind}\equiv\dot{E}_{\rm wind}/\dot{M}_{\rm BH}c^2\sim 1/1000$, in good agreement with the value required from large-scale galaxy simulations with AGN feedback.

In systems accreting well below the Eddington rate, the plasma in the innermost regions of the disk is collisionless and two-temperature, with the ions hotter than the electrons. Yet, whether a collisionless faster-than-Coulomb energy transfer mechanism exists in two-temperature accretion flows is still an open question. We study the physics of electron heating during the growth of ion velocity-space instabilities, by means of multi-dimensional particle-in-cell (PIC) simulations. A large-scale compression - embedded in a novel form of the PIC equations - continuously amplifies the field. This constantly drives a pressure anisotropy P_perp > P_parallel, due to the adiabatic invariance of the particle magnetic moments. We find that, for ion plasma beta values beta_i ~ 5-30 appropriate for the midplane of low-luminosity accretion flows, mirror modes dominate if the electron-to-proton temperature ratio is > 0.2, whereas if it is < 0.2 the ion cyclotron instability triggers the growth of strong Alfven-like waves, that pitch-angle scatter the ions to maintain marginal stability. We develop an analytical model of electron heating during the growth of the ion cyclotron instability, which we validate with PIC simulations. We find that for cold electrons (beta_e < m_e/m_i), the electron energy gain is controlled by the magnitude of the E-cross-B velocity induced by the ion cyclotron waves. This term is independent of the initial electron temperature, so it provides a solid energy floor even for electrons starting with extremely low temperatures. On the other hand, the electron energy gain for beta_e > m_e/m_i - governed by the conservation of the magnetic moment in the growing fields of the instability - is proportional to the initial electron temperature. Our results have implications for two-temperature accretion flows as well as the solar wind and intracluster plasmas. [abridged]

Electron acceleration to non-thermal energies is known to occur in low Mach number (M<5) shocks in galaxy clusters and solar flares, but the electron acceleration mechanism remains poorly understood. Using two-dimensional (2D) particle-in-cell (PIC) plasma simulations, we showed in Paper I that electrons are efficiently accelerated in low Mach number (M=3) quasi-perpendicular shocks via a Fermi-like process. The electrons bounce between the upstream region and the shock front, with each reflection at the shock resulting in energy gain via shock drift acceleration. The upstream scattering is provided by oblique magnetic waves, that are self-generated by the electrons escaping ahead of the shock. In the present work, we employ additional 2D PIC simulations to address the nature of the upstream oblique waves. We find that the waves are generated by the shock-reflected electrons via the firehose instability, which is driven by an anisotropy in the electron velocity distribution. We systematically explore how the efficiency of wave generation and of electron acceleration depend on the magnetic field obliquity, the flow magnetization (or equivalently, the plasma beta), and the upstream electron temperature. We find that the mechanism works for shocks with high plasma beta (>20) at nearly all magnetic field obliquities, and for electron temperatures in the range relevant for galaxy clusters. Our findings offer a natural solution to the conflict between the bright radio synchrotron emission observed from the outskirts of galaxy clusters and the low electron acceleration efficiency usually expected in low Mach number shocks.

Recent advances in general relativistic magnetohydrodynamic simulations have expanded and improved our understanding of the dynamics of black-hole accretion disks. However, current simulations do not capture the thermodynamics of electrons in the low density accreting plasma. This poses a significant challenge in predicting accretion flow images and spectra from first principles. Because of this, simplified emission models have often been used, with widely different configurations (e.g., disk- versus jet-dominated emission), and were able to account for the observed spectral properties of accreting black-holes. Exploring the large parameter space introduced by such models, however, requires significant computational power that exceeds conventional computational facilities. In this paper, we use GRay, a fast GPU-based ray-tracing algorithm, on the GPU cluster El Gato, to compute images and spectra for a set of six general relativistic magnetohydrodynamic simulations with different magnetic field configurations and black-hole spins. We also employ two different parametric models for the plasma thermodynamics in each of the simulations. We show that, if only the spectral properties of Sgr A* are used, all twelve models tested here can fit the spectra equally well. However, when combined with the measurement of the image size of the emission using the Event Horizon Telescope, current observations rule out all models with strong funnel emission, because the funnels are typically very extended. Our study shows that images of accretion flows with horizon-scale resolution offer a powerful tool in understanding accretion flows around black-holes and their thermodynamic properties.

Space astronomy in the last 40 years has largely been done from spacecraft in low Earth orbit (LEO) for which the technology is proven and delivery mechanisms are readily available. However, new opportunities are arising with the surge in commercial aerospace missions. We describe here one such possibility: deploying a small instrument on the Moon. This can be accomplished by flying onboard the Indian entry to the Google Lunar X PRIZE competition, Team Indus mission, which is expected to deliver a nearly 30 kgs of payloads to the Moon, with a rover as its primary payload. We propose to mount a wide-field far-UV (130--180 nm) imaging telescope as a payload on the Team Indus lander. Our baseline operation is a fixed zenith pointing but with the option of a mechanism to allow observations of different attitudes. Pointing towards intermediate ecliptic latitude (50 deg or above) ensures that the Sun is at least 40 deg off the line of sight at all times. In this position, the telescope can cover higher galactic latitudes as well as parts of Galactic plane. The scientific objectives of such a prospective are delineated and discussed.

The transient Swift J1644+57 is believed to have been produced by an unlucky star wandering too close to a supermassive black hole (BH) leading to a tidal disruption event. This unusual flare displayed highly super-Eddington X-ray emission which likely originated in a relativistic, collimated jet. This presents challenges to modern accretion and jet theory as upper limits of prior BH activity, which we obtain from the radio afterglow of this event, imply that both the pre-disruption BH and stellar magnetic fluxes fall many orders of magnitude short of what is required to power the observed X-ray luminosity. We argue that a pre-existing, "fossil" accretion disc can contain a sufficient reservoir of magnetic flux and that the stellar debris stream is capable of dragging this flux into the BH. To demonstrate this, we perform local, 3D magnetohydrodynamic simulations of the disc--stream interaction and demonstrate that the interface between the two is unstable to mixing. This mixing entrains a sufficient amount of fossil disc magnetic flux into the infalling stellar debris to power the jet. We argue that the interaction with the fossil disc can have a pronounced effect on the structure and dynamics of mass fallback and likely the resulting transient. Finally, we describe possible ramifications of these interactions on unresolved problems in tidal disruption dynamics, in particular, the efficiency of debris circularization, and effects of the disruption on the preexisting black hole system. Animations online: http://goo.gl/T84tLs

Observations of the black hole in the center of the Milky Way with the Event Horizon Telescope at 1.3 mm have revealed a size of the emitting region that is smaller than the size of the black-hole shadow. This can be reconciled with the spectral properties of the source, if the accretion flow is seen at a relatively high inclination (50-60 degrees). Such an inclination makes the angular momentum of the flow, and perhaps of the black hole, nearly aligned with the angular momenta of the orbits of stars that lie within 3 arcsec from the black hole. We discuss the implications of such an alignment for the properties of the black hole and of its accretion flow. We argue that future Event-Horizon-Telescope observations will not only refine the inclination of Sgr A* but also measure precisely its orientation on the plane of the sky.

The image of the emission surrounding the black hole in the center of the Milky Way is predicted to exhibit the imprint of general relativistic (GR) effects, including the existence of a shadow feature and a photon ring of diameter ~50 microarcseconds. Structure on these scales can be resolved by millimeter-wavelength very long baseline interferometry (VLBI). However, strong-field GR features of interest will be blurred at lambda >= 1.3 mm due to scattering by interstellar electrons. The scattering properties are well understood over most of the relevant range of baseline lengths, suggesting that the scattering may be (mostly) invertible. We simulate observations of a model image of Sgr A* and demonstrate that the effects of scattering can indeed be mitigated by correcting the visibilities before reconstructing the image. This technique is also applicable to Sgr A* at longer wavelengths.

Using the Very Long Baseline Array, we have measured a trigonometric parallax for the micro quasar GRS 1915+105, which contains a black hole and a K-giant companion. This yields a direct distance estimate of 8.6 (+2.0,-1.6) kpc and a revised estimate for the mass of the black hole of 12.4 (+2.0,-1.8) Msun. GRS 1915+105 is at about the same distance as some HII regions and water masers associated with high-mass star formation in the Sagittarius spiral arm of the Galaxy. The absolute proper motion of GRS 1915+105 is -3.19 +/- 0.03 mas/y and -6.24 +/- 0.05 mas/y toward the east and north, respectively, which corresponds to a modest peculiar speed of 22 +/-24 km/s at the parallax distance, suggesting that the binary did not receive a large velocity kick when the black hole formed. On one observational epoch, GRS 1915+105 displayed superluminal motion along the direction of its approaching jet. Considering previous observations of jet motions, the jet in GRS 1915+105 can be modeled with a jet inclination to the line of sight of 60 +/- 5 deg and a variable flow speed between 0.65c and 0.81c, which possibly indicates deceleration of the jet at distances from the black hole >2000 AU. Finally, using our measurements of distance and estimates of black hole mass and inclination, we provisionally confirm our earlier result that the black hole is spinning very rapidly.

Electron acceleration to non-thermal energies in low Mach number (M<5) shocks is revealed by radio and X-ray observations of galaxy clusters and solar flares, but the electron acceleration mechanism remains poorly understood. Diffusive shock acceleration, also known as first-order Fermi acceleration, cannot be directly invoked to explain the acceleration of electrons. Rather, an additional mechanism is required to pre-accelerate the electrons from thermal to supra-thermal energies, so they can then participate in the Fermi process. In this work, we use two- and three-dimensional particle-in-cell plasma simulations to study electron acceleration in low Mach number shocks. We focus on the particle energy spectra and the acceleration mechanism in a reference run with M=3 and a quasi-perpendicular pre-shock magnetic field. We find that about 15 percent of the electrons can be efficiently accelerated, forming a non-thermal power-law tail in the energy spectrum with a slope of p~2.4. Initially, thermal electrons are energized at the shock front via shock drift acceleration. The accelerated electrons are then reflected back upstream, where their interaction with the incoming flow generates magnetic waves. In turn, the waves scatter the electrons propagating upstream back toward the shock, for further energization via shock drift acceleration. In summary, the self-generated waves allow for repeated cycles of shock drift acceleration, similarly to a sustained Fermi-like process. This mechanism offers a natural solution to the conflict between the bright radio synchrotron emission observed from the outskirts of galaxy clusters and the low electron acceleration efficiency usually expected in low Mach number shocks.

Astronomers have discovered two populations of black holes: (i) stellar-mass black holes with masses in the range 5 to 30 solar masses, millions of which are present in each galaxy in the universe, and (ii) supermassive black holes with masses in the range 10^6 to 10^{10} solar masses, one each in the nucleus of every galaxy. There is strong circumstantial evidence that all these objects are true black holes with event horizons. The measured masses of supermassive black hole are strongly correlated with properties of their host galaxies, suggesting that these black holes, although extremely small in size, have a strong influence on the formation and evolution of entire galaxies. Spin parameters have recently been measured for a handful of black holes. Based on the data, there is an indication that the kinetic power of at least one class of relativistic jet ejected from accreting black holes may be correlated with black hole spin. If verified, it would suggest that these jets are powered by a generalized Penrose process mediated by magnetic fields.

We present a sub-grid model that emulates the magnetic dynamo operating in magnetized accretion disks. We have implemented this model in the general relativisic radiation magnetohydrodynamic (GRRMHD) code \koral, using results from local shearing sheet simulations of the magnetorotational instability to fix the parameters of the dynamo. With the inclusion of this dynamo, we are able to run 2D axisymmetric GRRMHD simulations of accretion disks for arbitrarily long times. The simulated disks exhibit sustained turbulence, with the poloidal and toroidal magnetic field components driven towards a state similar to that seen in 3D studies. Using this dynamo code, we present a set of long-duration global simulations of super-Eddington, optically-thick disks around non-spinning and spinning black holes. Super-Eddington disks around non-rotating black holes exhibit a surprisingly large efficiency, $\eta\approx0.04$, independent of the accretion rate, where we measure efficiency in terms of the total energy output, both radiation and mechanical, flowing out to infinity. Super-Eddington disks around spinning black holes are even more efficient, and appear to extract black hole rotational energy through a process similar to the Blandford-Znajek mechanism. All the simulated models are characterized by highly super-Eddington radiative fluxes collimated along the rotation axis. We also present a set of simulations that were designed to have Eddington or slightly sub-Eddington accretion rates ($\dot{M} \lesssim 2\dot M_{\rm Edd}$). None of these models reached a steady state. Instead, the disks collapsed as a result of runaway cooling, presumably because of a thermal instability.

The formation of purely baryonic globular clusters with no gravitationally bound dark matter is still a theoretical challenge. We show that these objects might form naturally whenever there is a relative stream velocity between baryons and dark matter. The stream velocity causes a phase shift between linear modes of baryonic and dark matter perturbations, which translates to a spatial offset between the two components when they collapse. For a 2sigma (3sigma) density fluctuation, baryonic clumps with masses in the range 1e5 - 2.5e6 Msun (1e5 - 4e6 Msun) collapse outside the virial radii of their counterpart dark matter halos. These objects could survive as long-lived dark matter-free objects and might conceivably become globular clusters. In addition, their dark matter counterparts, which were deprived of gas, might become dark satellite galaxies.

Collisionless shocks in plasmas play an important role in space physics (Earth's bow shock) and astrophysics (supernova remnants, relativistic jets, gamma-ray bursts, high energy cosmic rays). While the formation of a fluid shock through the steepening of a large amplitude sound wave has been understood for long, there is currently no detailed picture of the mechanism responsible for the formation of a collisionless shock. We unravel the physical mechanism at work and show that an electromagnetic Weibel shock always forms when two relativistic collisionless, initially unmagnetized, plasma shells encounter. The predicted shock formation time is in good agreement with 2D and 3D particle-in-cell simulations of counterstreaming pair plasmas. By predicting the shock formation time, experimental setups aiming at producing such shocks can be optimised to favourable conditions.

In Gou et al. (2011), we reported that the black hole primary in the X-ray binary Cygnus X-1 is a near-extreme Kerr black hole with a spin parameter a*>0.95(3{\sigma}). We confirm this result while setting a new and more stringent limit: a*>0.983 at the 3{\sigma}(99.7%) level of confidence. The earlier work, which was based on an analysis of all three useful spectra that were then available, was possibly biased by the presence in these spectra of a relatively strong Compton power-law component: The fraction of the thermal seed photons scattered into the power law was f_s=23-31%, while the upper limit for reliable application of the continuum-fitting method is f_s<25%. We have subsequently obtained six additional spectra of Cygnus X-1 suitable for the measurement of spin. Five of these spectra are of high quality with f_s in the range 10% to 19%, a regime where the continuum-fitting method has been shown to deliver reliable results. Individually, the six spectra give lower limits on the spin parameter that range from a*>0.95 to a*>0.98, allowing us to conservatively conclude that the spin of the black hole is a*>0.983 (3{\sigma}).

We show that, in principle, a slowly evolving gravitationally collapsing perfect fluid cloud can asymptotically settle to a static spherically symmetric equilibrium configuration with a naked singularity at the center. We consider one such asymptotic final configuration with a finite outer radius, and construct a toy model in which it is matched to a Schwarzschild exterior geometry. We examine the properties of circular orbits in this model. We then investigate observational signatures of a thermal accretion disk in this spacetime, comparing them with the signatures expected for a disk around a black hole of the same mass. Several notable differences emerge. A disk around the naked singularity is much more luminous than one around an equivalent black hole. Also, the disk around the naked singularity has a spectrum with a high frequency power law segment that carries a major fraction of the total luminosity. Thus, at least some naked singularities can, in principle, be distinguished observationally from black holes of the same mass. We discuss possible implications of these results.

Black hole binaries exhibit a wide range of variability phenomena, from large-scale state changes to broadband noise and quasi-periodic oscillations, but the physical nature of much of this variability is poorly understood. We examine the variability properties of three GRMHD simulations of thin accretion disks around black holes of varying spin, producing light curves and power spectra as would be seen by observers. We find that the simulated power spectra show a broad feature at high frequency, which increases in amplitude with the inclination of the observer. We show that this high-frequency feature is a product of the Doppler effect and that its location is a function of the mass and spin of the black hole. This Doppler feature demonstrates that power spectral properties of the accretion disk can be tied to, and potentially used to determine, physical properties of the black hole.

In the case involving particles the necessary and sufficient condition for the Penrose process to extract energy from a rotating black hole is absorption of particles with negative energies and angular momenta. No torque at the black-hole horizon occurs. In this article we consider the case of arbitrary fields or matter described by an unspecified, general energy-momentum tensor $T_{\mu \nu}$ and show that the necessary and sufficient condition for extraction of a black hole's rotational energy is analogous to that in the mechanical Penrose process: absorption of negative energy and negative angular momentum. We also show that a necessary condition for the Penrose process to occur is for the Noether current (the conserved energy-momentum density vector) to be spacelike or past directed (timelike or null) on some part of the horizon. In the particle case, our general criterion for the occurrence of a Penrose process reproduces the standard result. In the case of relativistic jet-producing "magnetically arrested disks" we show that the negative energy and angular-momentum absorption condition is obeyed when the Blandford-Znajek mechanism is at work, and hence the high energy extraction efficiency up to $\sim 300\%$ found in recent numerical simulations of such accretion flows results from tapping the black hole's rotational energy through the Penrose process. We show how black-hole rotational energy extraction works in this case by describing the Penrose process in terms of the Noether current.

Black hole accretion flows can be divided into two broad classes: cold and hot. Cold accretion flows, which consist of cool optically thick gas, are found at relatively high mass accretion rates. Prominent examples are the standard thin disk, which occurs at a fraction of the Eddington mass accretion rate, and the slim disk at super-Eddington rates. These accretion flows are responsible for luminous systems such as active galactic nuclei radiating at or close to the Eddington luminosity and black hole X-ray binaries in the soft state. Hot accretion flows, the topic of this review, are virially hot and optically thin. They occur at lower mass accretion rates, and are described by models such as the advection-dominated accretion flow and luminous hot accretion flow. Because of energy advection, the radiative efficiency of these flows is in general lower than that of a standard thin accretion disk. Moreover, the efficiency decreases with decreasing mass accretion rate. Observations show that hot accretion flows are associated with jets. In addition, theoretical arguments suggest that hot flows should produce strong winds. Hot accretion flows are believed to be present in low-luminosity active galactic nuclei and in black hole X-ray binaries in the hard and quiescent states. The prototype is Sgr A*, the ultra-low-luminosity supermassive black hole at our Galactic center. The jet, wind and radiation from a supermassive black hole with a hot accretion flow can interact with the external interstellar medium and modify the evolution of the host galaxy. Details of this "maintenance-mode feedback" could, in principle, be worked out through theoretical studies and numerical simulations of hot accretion flows.

Black hole (BH) accretion flows and jets are dynamic hot relativistic magnetized plasma flows whose radiative opacity can significantly affect flow structure and behavior. We describe a numerical scheme, tests, and an astrophysically relevant application using the M1 radiation closure within a new three-dimensional (3D) general relativistic (GR) radiation (R) magnetohydrodynamics (MHD) massively parallel code called HARMRAD. Our 3D GRRMHD simulation of super-Eddington accretion (about $20$ times Eddington) onto a rapidly rotating BH (dimensionless spin $j=0.9375$) shows sustained non-axisymmemtric disk turbulence, a persistent electromagnetic jet driven by the Blandford-Znajek effect, and a total radiative output consistently near the Eddington rate. The total accretion efficiency is of order $20\%$, the large-scale electromagnetic jet efficiency is of order $10\%$, and the total radiative efficiency that reaches large distances remains low at only order $1\%$. However, the radiation jet and the electromagnetic jet both emerge from a geometrically beamed polar region, with super-Eddington isotropic equivalent luminosities. Such simulations with HARMRAD can enlighten the role of BH spin vs.\ disks in launching jets, help determine the origin of spectral and temporal states in x-ray binaries, help understand how tidal disruption events (TDEs) work, provide an accurate horizon-scale flow structure for M87 and other active galactic nuclei (AGN), and isolate whether AGN feedback is driven by radiation or by an electromagnetic, thermal, or kinetic wind/jet. For example, the low radiative efficiency and weak BH spin-down rate from our simulation suggest that BH growth over cosmological times to billions of solar masses by redshifts of $z\sim 6-8$ is achievable even with rapidly rotating BHs and ten solar mass BH seeds.

Using long-duration general relativistic magnetohydrodynamic simulations of radiatively inefficient accretion discs, the energy, momentum and mass outflow rates from such systems are estimated. Outflows occur via two fairly distinct modes: a relativistic jet and a sub-relativistic wind. The jet power depends strongly on the black hole spin and on the magnetic flux at the horizon. Unless these are very small, the energy output in the jet dominates over that in the wind. For a rapidly spinning black hole accreting in the magnetically arrested limit, it is confirmed that jet power exceeds the total rate of accretion of rest mass energy. However, because of strong collimation, the jet probably does not have a significant feedback effect on its immediate surroundings. The power in the wind is more modest and shows a weaker dependence on black hole spin and magnetic flux. Nevertheless, because the wind subtends a large solid angle, it is expected to provide efficient feedback on a wide range of scales inside the host galaxy. Empirical formulae are obtained for the energy and momentum outflow rates in the jet and the wind.

A new general relativistic radiation magnetohydrodynamical code KORAL, is described, which employs the M1 scheme to close the radiation moment equations. The code has been successfully verified against a number of tests. Axisymmetric simulations of super-critical magnetized accretion on a non-rotating black hole (a=0.0) and a spinning black hole (a=0.9) are presented. The accretion rates in the two models are \dot M = 100-200 \dot M_Edd. These first general relativistic simulations of super-critical black hole accretion are potentially relevant to tidal disruption events and hyper-accreting supermassive black holes in the early universe. Both simulated models are optically and geometrically thick, and have funnels through which energy escapes in the form of relativistic gas, Poynting flux and radiative flux. The jet is significantly more powerful in the a=0.9 run. The net energy outflow rate in the two runs correspond to efficiencies of 5% (a=0) and 33% (a=0.9), as measured with respect to the mass accretion rate at the black hole. These efficiencies agree well with those measured in previous simulations of non-radiative geometrically thick disks. Furthermore, in the a=0.9 run, the outflow power appears to originate in the spinning black hole, suggesting that the associated physics is again similar in non-radiative and super-critical accretion flows. While the two simulations are efficient in terms of total energy outflow, both runs are radiatively inefficient. Their luminosities are only \sim 1-10 L_Edd, which corresponds to a radiative efficiency \sim 0.1%. Interestingly, most of the radiative luminosity emerges through the funnels, which subtend a very small solid angle. Therefore, measured in terms of a local radiative flux, the emitted radiation is highly super-Eddington.

General relativistic magnetohydrodynamic (GRMHD) simulations are providing influential models for black hole spin measurements, gamma ray bursts, and supermassive black hole feedback. Many of these simulations use the same initial condition: a rotating torus of fluid in hydrostatic equilibrium. A persistent concern is that simulation results sometimes depend on arbitrary features of the initial torus. For example, the Bernoulli parameter (which is related to outflows), appears to be controlled by the Bernoulli parameter of the initial torus. In this paper, we give a new equilibrium torus solution and describe two applications for the future. First, it can be used as a more physical initial condition for GRMHD simulations than earlier torus solutions. Second, it can be used in conjunction with earlier torus solutions to isolate the simulation results that depend on initial conditions. We assume axisymmetry, an ideal gas equation of state, constant entropy, and ignore self-gravity. We fix an angular momentum distribution and solve the relativistic Euler equations in the Kerr metric. The Bernoulli parameter, rotation rate, and geometrical thickness of the torus can be adjusted independently. Our torus tends to be more bound and have a larger radial extent than earlier torus solutions. While this paper was in preparation, several GRMHD simulations appeared based on our equilibrium torus. We believe it will continue to provide a more realistic starting point for future simulations.

Recently it has been observed that the scaling of jet power with black hole spin in galactic X-ray binaries is consistent with the predictions of the Blandford-Znajek (BZ) jet model. These observations motivate us to revisit the BZ model using general relativistic magnetohydrodynamic simulations of magnetized jets from accreting (h/r ~ 0.3), spinning (0 < a_* < 0.98) black holes. We have three main results. First, we quantify the discrepancies between the BZ jet power and our simulations: assuming maximum efficiency and uniform fields on the horizon leads to a ~10% overestimate of jet power, while ignoring the accretion disk leads to a further ~50% overestimate. Simply reducing the standard BZ jet power prediction by 60% gives a good fit to our simulation data. Our second result is to show that the membrane formulation of the BZ model correctly describes the physics underlying simulated jets: torques, dissipation, and electromagnetic fields on the horizon. This provides intuitive yet rigorous pictures for the black hole energy extraction process. Third, we compute the effective resistance of the load region and show that the load and the black hole achieve near perfect impedance matching. Taken together, these results increase our confidence in the BZ model as the correct description of jets observed from astrophysical black holes.

The standard thin accretion disc model predicts that discs around stellar mass black holes become radiation pressure dominated and thermally unstable once their luminosity exceeds L>0.02 L_Edd. Observationally, discs in the high/soft state of X-ray binaries show little variability in the range 0.01 L_Edd < L < 0.5 L_Edd, implying that these discs in nature are in fact quite stable. In an attempt to reconcile this conflict, we investigate one-zone disc models including turbulent and convective modes of vertical energy transport. We find both mixing mechanisms to have a stabilizing effect, leading to an increase in the L threshold up to which the disc is thermally stable. In the case of stellar mass black hole systems, convection alone leads to only a minor increase in this threshold, up to ~5 per cent of Eddington. However turbulent mixing has a much greater effect -- the threshold rises up to ~20 per cent Eddington under reasonable assumptions. In optimistic models with superefficient turbulent mixing, we even find solutions that are completely thermally stable for all accretion rates. Similar results are obtained for supermassive black holes, except that all critical accretion rates are a factor ~10 lower in Eddington ratio.

The spins of ten stellar black holes have been measured using the continuum-fitting method. These black holes are located in two distinct classes of X-ray binary systems, one that is persistently X-ray bright and another that is transient. Both the persistent and transient black holes remain for long periods in a state where their spectra are dominated by a thermal accretion disk component. The spin of a black hole of known mass and distance can be measured by fitting this thermal continuum spectrum to the thin-disk model of Novikov and Thorne; the key fit parameter is the radius of the inner edge of the black hole's accretion disk. Strong observational and theoretical evidence links the inner-disk radius to the radius of the innermost stable circular orbit, which is trivially related to the dimensionless spin parameter a_* of the black hole (|a_*| < 1). The ten spins that have so far been measured by this continuum-fitting method range widely from a_* \approx 0 to a_* > 0.95. The robustness of the method is demonstrated by the dozens or hundreds of independent and consistent measurements of spin that have been obtained for several black holes, and through careful consideration of many sources of systematic error. Among the results discussed is a dichotomy between the transient and persistent black holes; the latter have higher spins and larger masses. Also discussed is recently discovered evidence in the transient sources for a correlation between the power of ballistic jets and black hole spin.

Magnetic fields appear to be present in all galaxies and galaxy clusters. Recent measurements indicate that a weak magnetic field may be present even in the smooth low density intergalactic medium. One explanation for these observations is that a seed magnetic field was generated by some unknown mechanism early in the life of the Universe, and was later amplified by various dynamos in nonlinear objects like galaxies and clusters. We show that a primordial magnetic field is expected to be generated in the early Universe on purely linear scales through vorticity induced by scale-dependent temperature fluctuations or equivalently, a spatially varying speed of sound of the gas. Residual free electrons left over after recombination tap into this vorticity to generate magnetic field via the Biermann battery process. Although the battery operates even in the absence of any relative velocity between dark matter and gas at the time of recombination, the presence of such a relative velocity modifies the predicted spatial power spectrum of the magnetic field. At redshifts of order a few tens, we estimate a root mean square field strength of order 1e-25 -1e-24 G on comoving scales~10 kpc. This field, which is generated purely from linear perturbations, is expected to be amplified significantly after reionization, and to be further boosted by dynamo processes during nonlinear structure formation.

We perform detailed magnetohydrodynamic simulations of the gas cloud G2 interacting with the accretion flow around the Galactic Center black hole Sgr A*. We take as our initial conditions a steady-state, converged solution of the accretion flow obtained earlier using the general-relativistic magnetohydrodynamic code HARM. Using the observed parameters for the cloud's orbit, we compute the interaction of the cloud with the ambient gas and identify the shock structure that forms ahead of the cloud. We show that for many configurations, the cloud front crosses orbit pericenter 7 to 9 months earlier than the center-of-mass.

Collisionless shocks are ubiquitous in astrophysics and in the lab. Recent numerical simulations and experiments have shown how they can arise from the encounter of two collisionless plasma shells. When the shells interpenetrate, the overlapping region turns unstable, triggering the shock formation. As a first step towards a microscopic understanding of the process, we analyze here in detail the initial instability phase. On the one hand, 2D relativistic PIC simulations are performed where two symmetric initially cold pair plasmas collide. On the other hand, the instabilities at work are analyzed, as well as the field at saturation and the seed field which gets amplified. For mildly relativistic motions and onward, Weibel modes govern the linear phase. We derive an expression for the duration of the linear phase in good agreement with the simulations. This saturation time constitutes indeed a lower-bound for the shock formation time.

We calculate radio light curves produced by the bow shock that is likely to form in front of the G2 cloud when it penetrates the accretion disk of Sgr A*. The shock acceleration of the radio-emitting electrons is captured self-consistently by means of first-principles particle-in-cell simulations. We show that the radio luminosity is expected to reach maximum in early 2013, roughly a month after the bow shock crosses the orbit pericenter. We estimate the peak radio flux at 1.4 GHz to be 1.4 - 22 Jy depending on the assumed orbit orientation and parameters. We show that the most promising frequencies for radio observations are in the 0.1<nu<1 GHz range, for which the bow shock emission will be much stronger than the intrinsic radio flux for all the models considered.

It has for long been an article of faith among astrophysicists that black hole spin energy is responsible for powering the relativistic jets seen in accreting black holes. Two recent advances have strengthened the case. First, numerical general relativistic magnetohydrodynamic simulations of accreting spinning black holes show that relativistic jets form spontaneously. In at least some cases, there is unambiguous evidence that much of the jet energy comes from the black hole, not the disk. Second, spin parameters of a number of accreting stellar-mass black holes have been measured. For ballistic jets from these systems, it is found that the radio luminosity of the jet correlates with the spin of the black hole. This suggests a causal relationship between black hole spin and jet power, presumably due to a generalized Penrose process.

A numerical scheme is described for including radiation in multi-dimensional general-relativistic conservative fluid dynamics codes. In this method, a covariant form of the M1 closure scheme is used to close the radiation moments, and the radiative source terms are treated semi-implicitly in order to handle both optically thin and optically thick regimes. The scheme has been implemented in a conservative general relativistic radiation hydrodynamics code KORAL. The robustness of the code is demonstrated on a number of test problems, including radiative relativistic shock tubes, static radiation pressure supported atmosphere, shadows, beams of light in curved spacetime, and radiative Bondi accretion. The advantages of M1 closure relative to other approaches such as Eddington closure and flux-limited diffusion are discussed, and its limitations are also highlighted.

Using 5 GHz radio luminosity at light-curve maximum as a proxy for jet power and black-hole spin measurements obtained via the continuum-fitting method, Narayan & McClintock (2012) presented the first direct evidence for a relationship between jet power and black hole spin for four transient black-hole binaries. We test and confirm their empirical relationship using a fifth source, H1743-322, whose spin was recently measured. We show that this relationship is consistent with Fe-line spin measurements provided that the black hole spin axis is assumed to be aligned with the binary angular momentum axis. We also show that, during a major outburst of a black hole transient, the system reasonably approximates an X-ray standard candle. We further show, using the standard synchrotron bubble model, that the radio luminosity at light-curve maximum is a good proxy for jet kinetic energy. Thus, the observed tight correlation between radio power and black hole spin indicates a strong underlying link between mechanical jet power and black hole spin. Using the fitted correlation between radio power and spin for the above five calibration sources, we predict the spins of six other black holes in X-ray/radio transient systems with low-mass companions. Remarkably, these predicted spins are all relatively low, especially when compared to the high measured spins of black holes in persistent, wind-fed systems with massive companions.

We present a first attempt to model the narrow-line (NL) region of active galactic nuclei (AGN) in hydrodynamic simulations of galaxy mergers, using a novel physical prescription. This model is used to determine the origin of double-peaked NL (dNL) AGN in merging galaxies and their connection to supermassive black hole (SMBH) pairs, motivated by recent observations of such objects. We find that dNL AGN induced by the relative motion of SMBH pairs are a generic but short-lived feature of gaseous major mergers. dNL AGN should often be observed in late-stage mergers, during the kpc-scale phase of SMBH inspiral or soon after the SMBH merger. However, even within the kpc-scale phase, only a minority of dNL AGN are directly induced by SMBH motion; their lifetimes are typically a few Myr. Most double peaks arise from gas kinematics near the SMBH, although prior to the SMBH merger up to 80% of all dNL profiles may be influenced by SMBH motion via altered peak ratios or velocity offsets. The total lifetimes of dNL AGN depend strongly on viewing angle and on properties of the merging galaxies. Also, in a typical merger, at least 10-40% of the double peaks induced by SMBH motion have small projected separations, 0.1-1 kpc, such that dual peaks of stellar surface brightness are not easily resolved. Diffuse tidal features can indicate late-stage galaxy mergers, although they do not distinguish SMBH pairs from merged SMBHs. We show that dNL profiles with peak velocity splittings > 500 km s^-1 or with measurable overall velocity shifts are often associated with SMBH pairs. Our results support the notion that selection of dNL AGN is a promising method for identifying dual SMBH candidates, but demonstrate the critical importance of high-resolution, multi-wavelength follow-up observations, and the use of multiple lines of evidence, for confirming the dual nature of candidate SMBH pairs. (Abridged)

Almost all hydrodynamic accretion disk models parametrize viscosity with the dimensionless parameter alpha. There is no detailed model for alpha, so it is usually taken to be a constant. However, global simulations of magnetohydrodynamic disks find that alpha varies with distance from the central object. Also, Newtonian simulations tend to find smaller alpha's than general relativistic simulations. We seek a one-dimensional model for alpha that can reproduce these two observations. We are guided by data from six general relativistic magnetohydrodynamic accretion disk simulations. The variation of alpha in the inner, laminar regions of the flow results from stretching of mean magnetic field lines by the flow. The variation of alpha in the outer, turbulent regions results from the dependence of the magnetorotational instability on the dimensionless shear rate. We give a one-dimensional prescription for alpha(r) that captures these two effects and reproduces the radial variation of alpha observed in the simulations. For thin disks, the prescription simplifies to the formula alpha(r)=0.025[q(r)/1.5]^6, where the shear parameter, q(r), is an analytical function of radius in the Kerr metric. The coefficient and exponent are inferred from our simulations and will change as better simulation data becomes available. We conclude that the alpha-viscosity prescription can be extended to the radially varying alpha's observed in simulations. It is possible that Newtonian simulations find smaller alpha's than general relativistic simulations because the shear parameter is lower in Newtonian flows.

A dense ionized cloud of gas has been recently discovered to be moving directly toward the supermassive black hole, Sgr A*, at the Galactic Center. In June 2013, at the pericenter of its highly eccentric orbit, the cloud will be approximately 3100 Schwarzschild radii from the black hole and will move supersonically through the ambient hot gas with a velocity of v_p ~ 5400 km/s. A bow shock is likely to form in front of the cloud and could accelerate electrons to relativistic energies. We estimate via particle-in-cell simulations the energy distribution of the accelerated electrons and show that the non-thermal synchrotron emission from these electrons might exceed the quiescent radio emission from Sgr A* by a factor of several. The enhanced radio emission should be detectable at GHz and higher frequencies around the time of pericentric passage and in the following months. The bow shock emission is expected to be displaced from the quiescent radio emission of Sgr A* by ~33 mas. Interferometric observations could resolve potential changes in the radio image of Sgr A* at wavelengths < 6 cm.

We investigate the distribution of neutron star masses in different populations of binaries, employing Bayesian statistical techniques. In particular, we explore the differences in neutron star masses between sources that have experienced distinct evolutionary paths and accretion episodes. We find that the distribution of neutron star masses in non-recycled eclipsing high-mass binaries as well as of slow pulsars, which are all believed to be near their birth masses, has a mean of 1.28 M_solar and a dispersion of 0.24 M_solar. These values are consistent with expectations for neutron star formation in core-collapse supernovae. On the other hand, double neutron stars, which are also believed to be near their birth masses, have a much narrower mass distribution, peaking at 1.33 M_solar but with a dispersion of only 0.05 M_solar. Such a small dispersion cannot easily be understood and perhaps points to a particular and rare formation channel. The mass distribution of neutron stars that have been recycled has a mean of 1.48 M_solar and a dispersion of 0.2 M_solar, consistent with the expectation that they have experienced extended mass accretion episodes. The fact that only a very small fraction of recycled neutron stars in the inferred distribution have masses that exceed ~2 M_solar suggests that only a few of these neutron stars cross the mass threshold to form low mass black holes.

We present results from two long-duration GRMHD simulations of advection-dominated accretion around a non-spinning black hole. The first simulation was designed to avoid significant accumulation of magnetic flux around the black hole. This simulation was run for a time of 200,000GM/c^3 and achieved inflow equilibrium out to a radius \sim90GM/c^2. Even at this relatively large radius, the mass outflow rate \dot{M}_{out} is found to be only 60% of the net mass inflow rate \dot{M}_{BH} into the black hole. The second simulation was designed to achieve substantial magnetic flux accumulation around the black hole in a magnetically arrested disc. This simulation was run for a shorter time of 100,000GM/c^3. Nevertheless, because the mean radial velocity was several times larger than in the first simulation, it reached inflow equilibrium out to a radius \sim170GM/c^2. Here, \dot{M}_{out} becomes equal to \dot{M}_{BH} at r\sim 160GM/c^2. Since the mass outflow rates in the two simulations do not show robust convergence with time, it is likely that the true outflow rates are lower than our estimates. The effect of black hole spin on mass outflow remains to be explored. Neither simulation shows strong evidence for convection, though a complete analysis including the effect of magnetic fields is left for the future.

It is generally thought that the light coming from the inner plunging region of black hole accretion discs contributes negligibly to the disc's overall spectrum, i.e. the plunging fluid is swallowed by the black hole before it has time to radiate. In the standard disc model used to fit X-ray observations of accretion discs, the plunging region is assumed to be perfectly dark. However, numerical simulations that include the full physics of the magnetized flow predict that a small fraction of the disc's total luminosity emanates from the plunging region. We investigate the observational consequences of this neglected inner light. We compute radiative transfer based disc spectra that correspond to 3D general relativistic magnetohydrodynamic simulated discs (which produce light inside their plunging regions). In the context of black hole spin estimation, we find that the neglected inner light only has a modest effect (this bias is less than typical observational systematic errors). For rapidly spinning black holes, we find that the combined emission from the plunging region produces a weak power-law tail at high energies. This indicates that infalling matter is the origin for some of the `coronal' emission observed in the thermal dominant and steep power-law states of X-ray binaries.

We present global solutions of optically thin, two-temperature black hole accretion disks incorporating magnetic fields. We assume that the {\pi}{\phi}-component of the Maxwell stress is proportional to the total pressure, and prescribe the radial dependence of the magnetic flux advection rate in order to complete the set of basic equations. We obtained magnetically supported (low-{\beta}) disk solutions, whose luminosity exceeds the maximum luminosity for an advection-dominated accretion flow (ADAF), L > 0.4 {\alpha}^2 L_Edd, where L_Edd is the Eddington luminosity. The accretion flow is composed of the outer ADAF, a luminous hot accretion flow (LHAF) inside the transition layer from the outer ADAF to the low-{\beta} disk, the low-{\beta} disk, and the inner ADAF. The low-{\beta} disk region becomes wider as the mass-accretion rate increases further. In the low-{\beta} disk, the magnetic heating balances the radiative cooling, and the electron temperature decreases from ~ 10^9.5 K to ~ 10^8 K as the luminosity increases. These results are consistent with the anti-correlation between the energy cutoff in X-ray spectra (hence the electron temperature) and the luminosity when L > 0.1 L_Edd, observed in the bright/hard state during the bright hard-to-soft transitions of transient outbursts in galactic black hole candidates.

Recent advances in general relativistic magnetohydrodynamic modeling of jets offer unprecedented insights into the inner workings of accreting black holes that power the jets in active galactic nuclei (AGN) and other accretion systems. I will present the results of recent studies that determine spin-dependence of jet power and discuss the implications for the AGN radio loud/quiet dichotomy and recent observations of high jet power in a number of AGN.

We show that the 5-GHz radio flux of transient ballistic jets in black hole binaries correlates with the dimensionless black hole spin parameter a* estimated via the continuum-fitting method. The data suggest that jet power scales either as the square of a* or the square of the angular velocity of the horizon. This is the first direct evidence that jets may be powered by black hole spin energy. The observed correlation validates the continuum-fitting method of measuring spin. In addition, for those black holes that have well-sampled radio observations of ballistic jets, the correlation may be used to obtain rough estimates of their spins.

We develop here a new procedure within Einstein's theory of gravity to generate equilibrium configurations that result as the final state of gravitational collapse from regular initial conditions. As a simplification, we assume that the collapsing fluid is supported only by tangential pressure. We show that the equilibrium geometries generated by this method form a subset of static solutions to the Einstein equations, and that they can either be regular or develop a naked singularity at the center. When a singularity is present, there are key differences in the properties of stable circular orbits relative to those around a Schwarzschild black hole with the same mass. Therefore, if an accretion disk is present around such a naked singularity it could be observationally distinguished from a disk around a black hole.

We report a direct and accurate measurement of the distance to the X-ray binary Cygnus X-1, which contains the first black hole to be discovered. The distance of 1.86 (-0.11,+0.12) kpc was obtained from a trigonometric parallax measurement using the Very Long Baseline Array. The position measurements are also sensitive to the 5.6 d binary orbit and we determine the orbit to be clockwise on the sky. We also measured the proper motion of Cygnus X-1 which, when coupled to the distance and Doppler shift, gives the three-dimensional space motion of the system. When corrected for differential Galactic rotation, the non-circular (peculiar) motion of the binary is only about 21 km/s, indicating that the binary did not experience a large "kick" at formation.

Cygnus X-1 is a binary star system that is comprised of a black hole and a massive giant companion star in a tight orbit. Building on our accurate distance measurement reported in the preceding paper, we first determine the radius of the companion star, thereby constraining the scale of the binary system. To obtain a full dynamical model of the binary, we use an extensive collection of optical photometric and spectroscopic data taken from the literature. By using all of the available observational constraints, we show that the orbit is slightly eccentric (both the radial velocity and photometric data independently confirm this result) and that the companion star rotates roughly 1.4 times its pseudosynchronous value. We find a black hole mass of M =14.8\pm1.0 M_{\sun}, a companion mass of M_{opt}=19.2\pm1.9 M_{\sun}, and the angle of inclination of the orbital plane to our line of sight of i=27.1\pm0.8 deg.

The compact primary in the X-ray binary Cygnus X-1 was the first black hole to be established via dynamical observations. We have recently determined accurate values for its mass and distance, and for the orbital inclination angle of the binary. Building on these results, which are based on our favored (asynchronous) dynamical model, we have measured the radius of the inner edge of the black hole's accretion disk by fitting its thermal continuum spectrum to a fully relativistic model of a thin accretion disk. Assuming that the spin axis of the black hole is aligned with the orbital angular momentum vector, we have determined that Cygnus X-1 contains a near-extreme Kerr black hole with a spin parameter a/M>0.95 (3\sigma). For a less probable (synchronous) dynamical model, we find a/M>0.92 (3\sigma). In our analysis, we include the uncertainties in black hole mass, orbital inclination angle and distance, and we also include the uncertainty in the calibration of the absolute flux via the Crab. These four sources of uncertainty totally dominate the error budget. The uncertainties introduced by the thin-disk model we employ are particularly small in this case given the extreme spin of the black hole and the disk's low luminosity.

We describe global, 3D, time-dependent, non-radiative, general-relativistic, magnetohydrodynamic simulations of accreting black holes (BHs). The simulations are designed to transport a large amount of magnetic flux to the center, more than the accreting gas can force into the BH. The excess magnetic flux remains outside the BH, impedes accretion, and leads to a magnetically arrested disc. We find powerful outflows. For a BH with spin parameter a = 0.5, the efficiency with which the accretion system generates outflowing energy in jets and winds is eta ~ 30%. For a = 0.99, we find eta ~ 140%, which means that more energy flows out of the BH than flows in. The only way this can happen is by extracting spin energy from the BH. Thus the a = 0.99 simulation represents an unambiguous demonstration, within an astrophysically plausible scenario, of the extraction of net energy from a spinning BH via the Penrose-Blandford-Znajek mechanism. We suggest that magnetically arrested accretion might explain observations of active galactic nuclei with apparent eta ~ few x 100%.

We consider a model in which the ultra-relativistic jet in a gamma-ray burst (GRB) is cold and magnetically accelerated. We assume that the energy flux in the outflowing material is partially thermalized via internal shocks or a reverse shock, and we estimate the maximum amount of radiation that could be produced in such magnetized shocks. We compare this estimate with the available observational data on prompt gamma-ray emission in GRBs. We find that, even with highly optimistic assumptions, the magnetized jet model is radiatively too inefficient to be consistent with observations. One way out is to assume that much of the magnetic energy in the post-shock, or even pre-shock, jet material is converted to particle thermal energy by some unspecified process, and then radiated. This can increase the radiative efficiency sufficiently to fit observations. Alternatively, jet acceleration may be driven by thermal pressure rather than magnetic fields. In this case, which corresponds to the traditional fireball model, sufficient prompt GRB emission could be produced either from shocks at a large radius or from the jet photosphere closer to the center.

Rapid variability on a time scale much faster than the light-crossing time of the central supermassive black hole has been seen in TeV emission from the blazar PKS 2155-304. The most plausible explanation of this puzzling observation is that the radiating fluid in the relativistic jet is divided into a large number of sub-regions which move in random directions with relativistic speeds. We consider two versions of this "jets in a jet" model. In the first, the "subjets" model, stationary regions in the mean jet frame emit relativistic subjets that produce the observed radiation. The variability time scale is determined by the size of the sub-regions in the mean jet frame. This model, which is motivated by magnetic reconnection, has great difficulty explaining the observations in PKS 2155-304. In the alternate "turbulence" model, various sub-regions move relativistically in random directions and the variability time scale is determined by the size of these regions in their own comoving frames. This model fits the data much more comfortably. We consider collisions between TeV photons emitted from different sub-regions and find that, in both the subjets and turbulence models, the mean bulk Lorentz factor of the jet needs to be greater than 25 to avoid the pair catastrophe.

Angular momentum transport owing to hydrodynamic turbulent convection is studied using local three dimensional numerical simulations employing the shearing box approximation. We determine the turbulent viscosity from non-rotating runs over a range of values of the shear parameter and use a simple analytical model in order to extract the non-diffusive contribution (Lambda-effect) to the stress in runs where rotation is included. Our results suggest that the turbulent viscosity is of the order of the mixing length estimate and weakly affected by rotation. The Lambda-effect is non-zero and a factor of 2-4 smaller than the turbulent viscosity in the slow rotation regime. We demonstrate that for Keplerian shear, the angular momentum transport can change sign and be outward when the rotation period is greater than the turnover time, i.e. when the Coriolis number is below unity. This result seems to be relatively independent of the value of the Rayleigh number.

Active galactic nuclei (AGN), powered by long-term accretion onto central supermassive black holes, produce relativistic jets with lifetimes of greater than one million yr that preclude observations at birth. Transient accretion onto a supermassive black hole, for example through the tidal disruption of a stray star, may therefore offer a unique opportunity to observe and study the birth of a relativistic jet. On 2011 March 25, the Swift {\gamma}-ray satellite discovered an unusual transient source (Swift J164449.3+573451) potentially representing such an event. Here we present the discovery of a luminous radio transient associated with Swift J164449.3+573451, and an extensive set of observations spanning centimeter to millimeter wavelengths and covering the first month of evolution. These observations lead to a positional coincidence with the nucleus of an inactive galaxy, and provide direct evidence for a newly-formed relativistic outflow, launched by transient accretion onto a million solar mass black hole. While a relativistic outflow was not predicted in this scenario, we show that the tidal disruption of a star naturally explains the high-energy properties, radio luminosity, and the inferred rate of such events. The weaker beaming in the radio compared to {\gamma}-rays/X-rays, suggests that radio searches may uncover similar events out to redshifts of z ~ 6.

Slim-disk models describe accretion flows at high luminosities, while reducing to the standard thin disk form in the low luminosity limit. We have developed a new spectral model, slimbb, within the framework of XSPEC, which describes fully relativistic slim-disk accretion and includes photon ray-tracing that starts from the disk photosphere, rather than the equatorial plane. We demonstrate the features of this model by applying it to RXTE spectra of the persistent black-hole X-ray binary LMC X-3. LMC X-3 has the virtues of exhibiting large intensity variations while maintaining itself in soft spectral states which are well described using accretion-disk models, making it an ideal candidate to test the aptness of slimbb. Our results demonstrate consistency between the low-luminosity (thin-disk) and high luminosity (slim-disk) regimes. We also show that X-ray continuum-fitting in the high accretion rate regime can powerfully test black-hole accretion disk models.

A supermassive black hole in the nucleus of an elliptical galaxy at the centre of a cool-core group or cluster of galaxies is immersed in hot gas. Bondi accretion should occur at a rate determined by the properties of the gas at the Bondi radius and the mass of the black hole. X-ray observations of massive nearby elliptical galaxies, including M87 in the Virgo cluster, indicate a Bondi accretion rate Mdot which roughly matches the total kinetic power of the jets, suggesting that there is a tight coupling between the jet power and the mass accretion rate. While the Bondi model considers non-rotating gas, it is likely that the external gas has some angular momentum, which previous studies have shown could decrease the accretion rate drastically. We investigate here the possibility that viscosity acts at all radii to transport angular momentum outward so that the accretion inflow proceeds rapidly and steadily. The situation corresponds to a giant Advection Dominated Accretion Flow (ADAF) which extends from beyond the Bondi radius down to the black hole. We find solutions of the ADAF equations in which the gas accretes at just a factor of a few less than Mdot. These solutions assume that the atmosphere beyond the Bondi radius rotates with a sub-Keplerian velocity and that the viscosity parameter is large, alpha~0.1. The infall time of the ADAF solutions is no more than a few times the free-fall time. Thus the accretion rate at the black hole is closely coupled to the surrounding gas, enabling tight feedback to occur. We show that jet powers of a few per cent of Mdot c^2 are expected if either a fraction of the accretion power is channeled into the jet or the black hole spin energy is tapped by a strong magnetic field pressed against the black hole by the pressure of the accretion flow.(Truncated)

With an inferred bolometric luminosity exceeding 10^42 erg/s, HLX-1 in ESO 243-49 is the most luminous of ultraluminous X-ray sources and provides one of the strongest cases for the existence of intermediate mass black holes. We obtain good fits to disk-dominated observations of the source with BHSPEC, a fully relativistic black hole accretion disk spectral model. Due to degeneracies in the model arising from the lack of independent constraints on inclination and black hole spin, there is a factor of 100 uncertainty in the best-fit black hole mass M. Nevertheless, spectral fitting of XMM-Newton observations provides robust lower and upper limits with 3000 Msun < M < 3 x 10^5 Msun, at 90% confidence, placing HLX-1 firmly in the intermediate-mass regime. The lower bound on M is entirely determined by matching the shape and peak energy of the thermal component in the spectrum. This bound is consistent with (but independent of) arguments based solely on the Eddington limit. Joint spectral modelling of the XMM-Newton data with more luminous Swift and Chandra observations increases the lower bound to 6000 Msun, but this tighter constraint is not independent of the Eddington limit. The upper bound on M is sensitive to the maximum allowed inclination i, and is reduced to M < 10^5 Msun if we limit i < 75 deg.

The X-ray spectra of accretion discs of eight stellar-mass black holes have been analyzed to date using the thermal continuum fitting method, and the spectral fits have been used to estimate the spin parameters of the black holes. However, the underlying model used in this method of estimating spin is the general relativistic thin-disc model of Novikov & Thorne, which is only valid for razor-thin discs. We therefore expect errors in the measured values of spin due to inadequacies in the theoretical model. We investigate this issue by computing spectra of numerically calculated models of thin accretion discs around black holes, obtained via three-dimensional general relativistic magnetohydrodynamic (GRMHD) simulations. We apply the continuum fitting method to these computed spectra to estimate the black hole spins and check how closely the values match the actual spin used in the GRMHD simulations. We find that the error in the dimensionless spin parameter is up to about 0.2 for a non-spinning black hole, depending on the inclination. For black holes with spins of 0.7, 0.9 and 0.98, the errors are up to about 0.1, 0.03 and 0.01 respectively. These errors are comparable to or smaller than those arising from current levels of observational uncertainty. Furthermore, we estimate that the GRMHD simulated discs from which these error estimates are obtained correspond to effective disc luminosities of about 0.4-0.7 Eddington, and that the errors will be smaller for discs with luminosities of 0.3 Eddington or less, which are used in the continuum-fitting method. We thus conclude that use of the Novikov-Thorne thin-disc model does not presently limit the accuracy of the continuum-fitting method of measuring black hole spin.

A typical galaxy is thought to contain tens of millions of stellar-mass black holes, the collapsed remnants of once massive stars, and a single nuclear supermassive black hole. Both classes of black holes accrete gas from their environments. The accreting gas forms a flattened orbiting structure known as an accretion disk. During the past several years, it has become possible to obtain measurements of the spins of the two classes of black holes by modeling the X-ray emission from their accretion disks. Two methods are employed, both of which depend upon identifying the inner radius of the accretion disk with the innermost stable circular orbit (ISCO), whose radius depends only on the mass and spin of the black hole. In the Fe K method, which applies to both classes of black holes, one models the profile of the relativistically-broadened iron line with a special focus on the gravitationally redshifted red wing of the line. In the continuum-fitting method, which has so far only been applied to stellar-mass black holes, one models the thermal X-ray continuum spectrum of the accretion disk. We discuss both methods, with a strong emphasis on the continuum-fitting method and its application to stellar-mass black holes. Spin results for eight stellar-mass black holes are summarized. These data are used to argue that the high spins of at least some of these black holes are natal, and that the presence or absence of relativistic jets in accreting black holes is not entirely determined by the spin of the black hole.

We use dynamical mass measurements of 16 black holes in transient low-mass X-ray binaries to infer the stellar black hole mass distribution in the parent population. We find that the observations are best described by a narrow mass distribution at 7.8 +/- 1.2 Msolar. We identify a selection effect related to the choice of targets for optical follow-ups that results in a flux-limited sample. We demonstrate, however, that this selection effect does not introduce a bias in the observed distribution and cannot explain the absence of black holes in the 2-5 solar mass range. On the high mass end, we argue that the rapid decline in the inferred distribution may be the result of the particular evolutionary channel followed by low-mass X-ray binaries. This is consistent with the presence of high-mass black holes in the persistent, high-mass X-ray binary sources. If the paucity of low-mass black holes is caused by a sudden decrease of the supernova explosion energy with increasing progenitor mass, this would have observable implications for ongoing transient surveys that target core-collapse supernovae. Our results also have significant implications for the calculation of event rates from the coalescence of black hole binaries for gravitational wave detectors.

We measure the spin of XTE J1550-564 in two ways: by modelling the thermal continuum spectrum of the accretion disc, and independently by modeling the broad red wing of the reflection fluorescence Fe-K line. We find that the spin measurements conducted independently using both leading methods are in agreement with one another. For the continuum-fitting analysis, we use a data sample consisting of several dozen RXTE spectra, and for the Fe-K analysis, we use a pair of ASCA spectra from a single epoch. Our spin estimate for the black hole primary using the continuum-fitting method is -0.11 < a* < 0.71 (90 per cent confidence), with a most likely spin of a* = 0.34. In obtaining this result, we have thoroughly explored the dependence of the spin value on a wide range of model-dependent systematic errors and observational errors; our precision is limited by uncertainties in the distance and orbital inclination of the system. For the Fe-line method, our estimate of spin is a* = 0.55(+0.15,-0.22). Combining these results, we conclude that the spin of this black hole is moderate, a* = 0.49(+0.13,-0.20), which suggests that the jet of this microquasar is powered largely by its accretion disc rather than by the spin energy of the black hole.

During its yearlong outburst in 1975--76, the transient source A0620--00 reached an intensity of 50 Crab, an all-time record for any X-ray binary. The source has been quiescent since. We recently determined accurate values for the black hole mass, orbital inclination angle and distance. Building on these results, we have measured the radius of the inner edge of the accretion disk around the black hole primary by fitting its thermal continuum spectrum to our version of the relativistic Novikov-Thorne thin-disk model. We have thereby estimated the spin of the black hole. Although our spin estimate depends on a single high quality spectrum, which was obtained in 1975 by OSO-8, we are confident of our result because of the consistent values of the inner-disk radius that we have obtained for hundreds of observations of other sources: H1743-322, XTE J1550-564, and notably LMC X-3. We have determined the dimensionless spin parameter of the black hole to be a*=0.12+/-0.19, with a*<0.49 and a*>-0.59 at the 3 sigma level of confidence. This result takes into account all sources of observational and model-parameter uncertainties. Despite the low spin, the intensity and properties of the radio counterpart, both in outburst and quiescence, attest to the presence of a strong jet. If jets are driven by black hole spin, then current models indicate that jet power should be a steeply increasing function of a*. Consequently, the low spin of A0620--00 suggests that its jet may be disk-driven.

The black-hole binary system LMC X-3 has been observed by virtually every X-ray mission since the inception of X-ray astronomy. Among the persistent sources, LMC X-3 is uniquely both habitually soft and highly variable. Using a fully relativistic accretion-disk model, we analyze hundreds of spectra collected during eight X-ray missions that span 26 years. For a selected sample of 391 RXTE spectra we find that to within ~2 percent the inner radius of the accretion disk is constant over time and unaffected by source variability. Even considering an ensemble of eight X-ray missions, we find consistent values of the radius to within ~4-6 percent. Our results provide strong evidence for the existence of a fixed inner-disk radius. The only reasonable inference is that this radius is closely associated with the general relativistic innermost stable circular orbit (ISCO). Our findings establish a firm foundation for the measurement of black hole spin.

The standard general relativistic model of a razor-thin accretion disk around a black hole, developed by Novikov & Thorne (NT) in 1973, assumes the shear stress vanishes at the radius of the innermost stable circular orbit (ISCO) and that, outside the ISCO, the shear stress is produced by an effective turbulent viscosity. However, astrophysical accretion disks are not razor-thin, it is uncertain whether the shear stress necessarily vanishes at the ISCO, and the magnetic field, which is thought to drive turbulence in disks, may contain large-scale structures that do not behave like a simple local scalar viscosity. We describe three-dimensional general relativistic magnetohydrodynamic simulations of accretion disks around black holes with a range of spin parameters, and we use the simulations to assess the validity of the NT model. Our fiducial initial magnetic field consists of multiple (alternating polarity) poloidal field loops whose shape is roughly isotropic in the disk in order to match the isotropic turbulence expected in the poloidal plane. For a thin disk with an aspect ratio |h/r| ~ 0.07 around a non-spinning black hole, we find a decrease in the accreted specific angular momentum of 2.9% relative to the NT model and an excess luminosity from inside the ISCO of 3.5%. The deviations in the case of spinning black holes are of the same order. In addition, the deviations decrease with decreasing |h/r|. We therefore conclude that magnetized thin accretion disks in x-ray binaries in the thermal/high-soft spectral state ought to be well-described by the NT model, especially at luminosities below 30% of Eddington where we expect a very small disk thickness |h/r| <~ 0.05. We also discuss how the stress and the luminosity inside the ISCO depend on the assumed initial magnetic field geometry. (abridged)

In prior work, {\it Chandra} and Gemini-North observations of the eclipsing X-ray binary M33 X-7 have yielded measurements of the mass of its black hole primary and the system's orbital inclination angle of unprecedented accuracy. Likewise, the distance to the binary is known to a few percent. In an analysis based on these precise results, fifteen {\it Chandra} and {\it XMM-Newton} X-ray spectra, and our fully relativistic accretion disk model, we find that the dimensionless spin parameter of the black hole primary is $a_* = 0.77 \pm 0.05$. The quoted 1-$\sigma$ error includes all sources of observational uncertainty. Four {\it Chandra} spectra of the highest quality, which were obtained over a span of several years, all lead to the same estimate of spin to within statistical errors (2%), and this estimate is confirmed by 11 spectra of lower quality. There are two remaining uncertainties: (1) the validity of the relativistic model used to analyze the observations, which is being addressed in ongoing theoretical work; and (2) our assumption that the black hole spin is approximately aligned with the angular momentum vector of the binary, which can be addressed by a future X-ray polarimetry mission.

Achromatic breaks in afterglow light curves of gamma-ray bursts (GRBs) arise naturally if the product of the jet's Lorentz factor \gamma and opening angle \Theta_j satisfies (\gamma \Theta_j) >> 1 at the onset of the afterglow phase, i.e., soon after the conclusion of the prompt emission. Magnetohydrodynamic (MHD) simulations of collimated GRB jets generally give (\gamma \Theta_j) <~ 1, suggesting that MHD models may be inconsistent with jet breaks. We work within the collapsar paradigm and use axisymmetric relativistic MHD simulations to explore the effect of a finite stellar envelope on the structure of the jet. Our idealized models treat the jet-envelope interface as a collimating rigid wall, which opens up outside the star to mimic loss of collimation. We find that the onset of deconfinement causes a burst of acceleration accompanied by a slight increase in the opening angle. In our fiducial model with a stellar radius equal to 10^4.5 times that of the central compact object, the jet achieves an asymptotic Lorentz factor \gamma ~ 500 far outside the star and an asymptotic opening angle \Theta_j ~ 0.04 rad ~ 2 deg, giving (\gamma \Theta_j) ~ 20. These values are consistent with observations of typical long-duration GRBs, and explain the occurrence of jet breaks. We provide approximate analytic solutions that describe the numerical results well.

Radio loud active galactic nuclei (AGN) are on average 1000 times brighter in the radio band compared to radio quiet AGN. We investigate whether this radio loud/quiet dichotomy can be due to differences in the spin of the central black holes that power the radio-emitting jets. Using general relativistic magnetohydrodynamic simulations, we construct steady state axisymmetric numerical models for a wide range of black hole spins (dimensionless spin parameter 0.1 <= a <= 0.9999 and a variety of jet geometries. We assume that the total magnetic flux through the black hole horizon at radius r_H(a) is held constant. If the black hole is surrounded by a thin accretion disk, we find that the total black hole power output depends approximately quadratically on the angular frequency of the hole, P \propto \Omega_H^2 \propto (a/r_H)^2. We conclude that, in this scenario, differences in the black hole spin can produce power variations of only a few tens at most. However, if the disk is thick such that the jet subtends a narrow solid angle around the polar axis, then the power dependence becomes much steeper, P \propto \Omega_H^4 or even \propto \Omega_H^6. Power variations of 1000 are then possible for realistic black hole spin distributions. We derive an analytic solution that accurately reproduces the steeper scaling of jet power with \Omega_H, and we provide a numerical fitting formula that reproduces all our simulation results. We discuss other physical effects that might contribute to the observed radio loud/quiet dichotomy of AGN.

The connection between collimation and acceleration of magnetized relativistic jets is discussed. The focus is on recent numerical simulations which shed light on some longstanding problems.

We use the Novikov-Thorne thin disk model to fit the thermal continuum X-ray spectra of black hole X-ray binaries, and thereby extract the dimensionless spin parameter a* = a/M of the black hole as a parameter of the fit. We summarize the results obtained to date for six systems and describe work in progress on additional systems. We also describe recent methodological advances, our current efforts to make our analysis software fully available to others, and our theoretical efforts to validate the Novikov-Thorne model.

All prior work on measuring the spins of stellar-mass black holes via the X-ray continuum-fitting method has relied on the use of weakly-Comptonized spectra obtained in the thermal dominant state. Using a self-consistent Comptonization model, we show that one can analyze spectra that exhibit strong power-law components and obtain values of the inner disk radius, and hence spin, that are consistent with those obtained in the thermal dominant state. Specifically, we analyze many RXTE spectra of two black hole transients, H1743-322 and XTE J1550-564, and we demonstrate that the radius of the inner edge of the accretion disk remains constant to within a few percent as the strength of the Comptonized component increases by an order of magnitude, i.e., as the fraction of the thermal seed photons that are scattered approaches 25%. We conclude that the continuum-fitting method can be applied to a much wider body of data than previously thought possible, and to sources that have never been observed to enter the thermal dominant state (e.g., Cyg X-1).

The first extragalactic X-ray binary, LMC X-1, was discovered in 1969. In the 1980s, its compact primary was established as the fourth dynamical black-hole candidate. Recently, we published accurate values for the mass of the black hole and the orbital inclination angle of the binary system. Building on these results, we have analyzed 53 X-ray spectra obtained by RXTE and, using a selected sample of 18 of these spectra, we have determined the dimensionless spin parameter of the black hole to be a* = 0.92(-0.07,+0.05). This result takes into account all sources of observational and model-parameter uncertainties. The standard deviation around the mean value of a* for these 18 X-ray spectra, which were obtained over a span of several years, is only 0.02. When we consider our complete sample of 53 RXTE spectra, we find a somewhat higher value of the spin parameter and a larger standard deviation. Finally, we show that our results based on RXTE data are confirmed by our analyses of selected X-ray spectra obtained by the XMM-Newton, BeppoSAX and Ginga missions.

Unconfined relativistic outflows from rotating, magnetized compact objects are often well-modeled by assuming the field geometry is approximately a split-monopole at large radii. Earlier work has indicated that such an unconfined flow has an inefficient conversion of magnetic energy to kinetic energy. This has led to the conclusion that ideal magnetohydrodynamical (MHD) processes fail to explain observations of, e.g., the Crab pulsar wind at large radii where energy conversion appears efficient. In addition, as a model for astrophysical jets, the monopole field geometry has been abandoned in favor of externally confined jets since the latter appeared to be generically more efficient jet accelerators. We perform time-dependent axisymmetric relativistic MHD simulations in order to find steady state solutions for a wind from a compact object endowed with a monopole field geometry. Our simulations follow the outflow for 10 orders of magnitude in distance from the compact object, which is large enough to study both the initial "acceleration zone" of the magnetized wind as well as the asymptotic "coasting zone." We obtain the surprising result that acceleration is actually {\it efficient} in the polar region, which develops a jet despite not being confined by an external medium. Our models contain jets that have sufficient energy to account for moderately energetic long and short gamma-ray burst (GRB) events (~10^{51}--10^{52} erg), collimate into narrow opening angles (opening half-angle \theta_j \approx 0.03 rad), become matter-dominated at large radii (electromagnetic energy flux per unit matter energy flux \sigma<1), and move at ultrarelativistic Lorentz factors (\gamma_j ~ 200 for our fiducial model). (abridged)

We explore the connection between different classes of active galactic nuclei (AGNs) and the evolution of their host galaxies, by deriving host galaxy properties, clustering, and Eddington ratios of AGNs selected in the radio, X-ray, and infrared. We study a sample of 585 AGNs at 0.25 < z < 0.8 using redshifts from the AGN and Galaxy Evolution Survey (AGES) and data in the radio (WSRT 1.4 GHz), X-rays (Chandra XBootes), and mid-IR (IRAC Shallow Survey). The radio, X-ray, and IR AGN samples show modest overlap, indicating that to the flux limits of the survey, they represent largely distinct classes of AGNs. We derive host galaxy colors and luminosities, as well as Eddington ratios (lambda), for obscured or optically faint AGNs. We also measure the two-point cross-correlation between AGNs and galaxies on scales of 0.3-10 h^-1 Mpc, and derive typical dark matter halo masses. We find that: (1) radio AGNs are mainly found in luminous red galaxies, are strongly clustered (with M_halo ~ 3x10^13 h^-1 M_sun), and have very low lambda <~ 10^-3; (2) X-ray-selected AGNs are preferentially found in galaxies in the "green valley" of color-magnitude space and are clustered similarly to typical AGES galaxies (M_halo ~ 10^13 h^-1 M_sun), with 10^-3 <~ lambda <~ 1; (3) IR AGNs reside in slightly bluer, less luminous galaxies than X-ray AGNs, are weakly clustered (M_halo <~ 10^12 h^-1 M_sun), and have lambda > 10^-2. We interpret these results in terms of a simple model of AGN and galaxy evolution, whereby a "quasar" phase and the growth of the stellar bulge occurs when a galaxy's dark matter halo reaches a critical mass between ~10^12 and 10^13 M_sun. Subsequently, star formation ceases and AGN accretion shifts from radiatively efficient (optical- and IR- bright) to radiatively inefficient (optically-faint, radio-bright) modes.

We present a dynamical model of the high mass X-ray binary LMC X-1 based on high-resolution optical spectroscopy and extensive optical and near-infrared photometry. From our new optical data we find an orbital period of P=3.90917 +/- 0.00005 days. We present a refined analysis of the All Sky Monitor data from RXTE and find an X-ray period of P=3.9094 +/- 0.0008 days, which is consistent with the optical period. A simple model of Thomson scattering in the stellar wind can account for the modulation seen in the X-ray light curves. The V-K color of the star (1.17 +/- 0.05) implies A_V = 2.28 +/- 0.06, which is much larger than previously assumed. For the secondary star, we measure a radius of R_2 = 17.0 +/- 0.8 solar radii and a projected rotational velocity of V_rot*sin(i) = 129.9 +/- 2.2 km/s. Using these measured properties to constrain the dynamical model, we find an inclination of i = 36.38 +/- 1.92 deg, a secondary star mass of M_2 = 31.79 +/- 3.48 solar masses, and a black hole mass of 10.91 +/- 1.41 solar masses. The present location of the secondary star in a temperature-luminosity diagram is consistent with that of a star with an initial mass of 35 solar masses that is 5 Myr past the zero-age main sequence. The star nearly fills its Roche lobe (~90% or more), and owing to the rapid change in radius with time in its present evolutionary state, it will encounter its Roche lobe and begin rapid and possibly unstable mass transfer on a timescale of a few hundred thousand years.

Black hole event horizons, causally separating the external universe from compact regions of spacetime, are one of the most exotic predictions of General Relativity (GR). Until recently, their compact size has prevented efforts to study them directly. Here we show that recent millimeter and infrared observations of Sagittarius A* (Sgr A*), the supermassive black hole at the center of the Milky Way, all but requires the existence of a horizon. Specifically, we show that these observations limit the luminosity of any putative visible compact emitting region to below 0.4% of Sgr A*'s accretion luminosity. Equivalently, this requires the efficiency of converting the gravitational binding energy liberated during accretion into radiation and kinetic outflows to be greater than 99.6%, considerably larger than those implicated in Sgr A*, and therefore inconsistent with the existence of such a visible region. Finally, since we are able to frame this argument entirely in terms of observable quantities, our results apply to all geometric theories of gravity that admit stationary solutions, including the commonly discussed f(R) class of theories.

We present an empirical model of Comptonization for fitting the spectra of X-ray binaries. This model, "simpl", has been developed as a package implemented in XSPEC. With only two free parameters, simpl is competitive as the simplest empirical model of Compton scattering. Unlike other empirical models, such as the standard power-law model, simpl incorporates the basic physics of Compton scattering of soft photons by energetic coronal electrons. Using a simulated spectrum, we demonstrate that simpl closely matches the behavior of physical Comptonization models which consider the effects of optical depth, coronal electron temperature, and geometry. We present fits to RXTE spectra of the black-hole transient H1743-322 and a BeppoSAX spectrum of LMC X-3 using both simpl and the standard power-law model. A comparison of the results shows that simpl gives equally good fits and a comparable spectral index, while eliminating the troublesome divergence of the standard power-law model at low energies. Importantly, simpl is completely flexible and can be used self-consistently with any seed spectrum of photons. We show that simpl - unlike the standard power law - teamed up with diskbb (the standard model of disk accretion) gives results for the inner-disk radius that are unaffected by strong Comptonization, a result of great importance for the determination of black hole spin via the continuum-fitting method.

We consider a two-parameter family of cylindrical force-free equilibria, modeled to match numerical simulations of relativistic force-free jets. We study the linear stability of these equilibria, assuming a rigid impenetrable wall at the outer cylindrical radius R_j. We find that equilibria in which the Lorentz factor \gamma(R) increases monotonically with increasing radius R are stable. On the other hand, equilibria in which \gamma(R) reaches a maximum value at an intermediate radius and then declines to a smaller value \gamma_j at R_j are unstable. The most rapidly growing mode is an m=1 kink instability which has a growth rate ~ (0.4 / \gamma_j) (c/R_j). The e-folding length of the equivalent convected instability is ~2.5 \gamma_j R_j. For a typical jet with an opening angle \theta_j ~ few / \gamma_j, the mode amplitude grows weakly with increasing distance from the base of the jet, much slower than one might expect from a naive application of the Kruskal-Shafranov stability criterion.

We show that the excellent optical and gamma-ray data available for GRB 080319B rule out the internal shock model for the prompt emission. The data instead point to a model in which the observed radiation was produced close to the deceleration radius ($\sim10^{17}$ cm) by a turbulent source with random Lorentz factors $\sim10$ in the comoving frame. The optical radiation was produced by synchrotron emission from relativistic electrons, and the gamma-rays by inverse Compton scattering of the synchrotron photons. The gamma-ray emission originated both in eddies and in an inter-eddy medium, whereas the optical radiation was mostly from the latter. Therefore, the gamma-ray emission was highly variable whereas the optical was much less variable. The model explains all the observed features in the prompt optical and gamma-ray data of GRB 080319B. We are unable to determine with confidence whether the energy of the explosion was carried outward primarily by particles (kinetic energy) or magnetic fields. Consequently, we cannot tell whether the turbulent medium was located in the reverse shock (we can rule out the forward shock) or in a Poynting-dominated jet.

A popular paradigm to explain the rapid temporal variability observed in gamma-ray burst (GRB) lightcurves is the internal shock model. We propose an alternative model in which the radiating fluid in the GRB shell is relativistically turbulent with a typical eddy Lorentz factor $\gamma_t$. In this model, all pulses in the gamma-ray lightcurve are produced at roughly the same distance $R$ from the center of the explosion. The burst duration is $\sim R/c\Gamma^2$, where $\Gamma$ is the bulk Lorentz factor of the expanding shell, and the duration of individual pulses in the lightcurve is $\sim R/c\Gamma^2\gamma_t^2$. The model naturally produces highly variable lightcurves with $\sim\gamma_t^2$ individual pulses. Even though the model assumes highly inhomogeneous conditions, nevertheless the efficiency for converting jet energy to radiation is high.

Spin parameters of stellar-mass black holes in X-ray binaries are currently being estimated by fitting the X-ray continuum spectra of their accretion disk emission. For this method, it is necessary to know the inclination of the X-ray-producing inner region of the disk. Since the inner disk is expected to be oriented perpendicular to the spin axis of the hole, the usual practice is to assume that the black hole spin is aligned with the orbital angular momentum vector of the binary, and to estimate the inclination of the latter from ellipsoidal modulations in the light curve of the secondary star. We show that the inclination of the disk can be inferred directly if we have both spectral and polarization information on the disk radiation. The predicted degree of polarization varies from 0% to 5% as the disk inclination changes from face-on to edge-on. With current X-ray polarimetric techniques the polarization degree of a typical bright X-ray binary could be measured to an accuracy of 0.1% by observing the source for about 10 days. Such a measurement would constrain the disk inclination to within a degree or two and would significantly improve the reliability of black hole spin estimates. In addition, it would provide new information on the tilt between the black hole spin axis and the orbital rotation axis of the binary, which would constrain any velocity kicks experienced by stellar-mass black holes during their formation.

We describe three-dimensional general relativistic magnetohydrodynamic simulations of a geometrically thin accretion disk around a non-spinning black hole. The disk has a thickness $h/r\sim0.05-0.1$ over the radial range $(2-20)GM/c^2$. In steady state, the specific angular momentum profile of the inflowing magnetized gas deviates by less than 2% from that of the standard thin disk model of Novikov & Thorne (1973). Also, the magnetic torque at the radius of the innermost stable circular orbit (ISCO) is only $\sim2%$ of the inward flux of angular momentum at this radius. Both results indicate that magnetic coupling across the ISCO is relatively unimportant for geometrically thin disks.

We determine some basic properties of stars that produce spectacular gamma-ray bursts at the end of their life. We assume that accretion of the outer portion of the stellar core by a central black hole fuels the prompt emission, and that fall-back and accretion of the stellar envelope later produces the plateau in the X-ray light curve seen in some bursts. Using X-ray data for three bursts we estimate the radius of the stellar core to be ~ 1-3 x 10^10 cm, and that of the stellar envelope to be ~ 1-2 x 10^11 cm. The density profile in the envelope is fairly shallow, with \rho ~ r^-2. The rotation speeds of the core and envelope are ~ 0.05 and ~ 0.2 of the local Keplerian speed, respectively.

We calculate the rate of in-fall of stellar matter on an accretion disk during the collapse of a rapidly rotating massive star, and estimate the luminosity of the relativistic jet that results from accretion on to the central black hole. We find that the jet luminosity remains high for about 100 seconds, at a level comparable to the typical luminosity observed in gamma-ray bursts (GRBs). The luminosity then decreases rapidly with time for about 10^3 seconds, roughly as ~ t^-3; the duration depends on the size and rotation speed of the stellar core. The rapid decrease of the jet power explains the steeply declining X-ray flux observed at the end of most long duration GRBs. A X-ray plateau is also produced by continued fall-back of matter -- either from an extended stellar envelope or from material that failed to escape with the supernova ejecta. In a few GRBs, the X-ray light curve is observed to drop suddenly at the end of the plateau phase, while in others the decline is ~ t^-1 - t ^-2. These features arise naturally in the accretion model depending on the radius and mean specific angular momentum of the stellar envelope. The accretion model thus provides a coherent explanation for the diverse and puzzling features observed in the early X-ray light curves of GRBs. (Abridged)

Long-duration gamma-ray bursts (GRBs) require an engine capable of driving a jet of plasma to ultrarelativistic bulk Lorentz factors of up to several hundred and into narrow opening angles of a few degrees. We use global axisymmetric stationary solutions of magnetically-dominated (force-free) ultrarelativistic jets to test whether the popular magnetic-driving paradigm can generate the required Lorentz factors and opening angles. Our global solutions are obtained via time-dependent relativistic ideal magnetodynamical numerical simulations which follow the jet from the central engine to beyond six orders of magnitude in radius. Our model is primarily motivated by the collapsar model, in which a jet is produced by a spinning black hole or neutron star and then propagates through a massive stellar envelope. We find that the size of the presupernova progenitor star and the radial profile of pressure inside the star determine the terminal Lorentz factor and opening angle of the jet. At the radius where the jet breaks out of the star, our well-motivated fiducial model generates a Lorentz factor $\gamma\sim 400$ and a half-opening angle $\theta_j\sim 2^\circ$, consistent with observations of many long-duration GRBs. Other models with slightly different parameters give $\gamma$ in the range 100 to 5000 and $\theta_j$ from $0.1^\circ$ to $10^\circ$, thus reproducing the range of properties inferred for GRB jets. A potentially observable feature of some of our solutions is that the maximum Poynting flux in the jet is found at $\theta \sim \theta_j$ with the jet power concentrated in a hollow cone, while the maximum in the Lorentz factor occurs at an angle $\theta$ substantially smaller than $\theta_j$ also in a hollow cone. [abridged]

In a previous paper we demonstrated that non-radial hydrodynamic oscillations of a thermally-supported (Bonnor-Ebert) sphere embedded in a low-density, high-temperature medium persist for many periods. The predicted column density variations and molecular spectral line profiles are similar to those observed in the Bok globule B68 suggesting that the motions in some starless cores may be oscillating perturbations on a thermally supported equilibrium structure. Such oscillations can produce molecular line maps which mimic rotation, collapse or expansion, and thus could make determining the dynamical state from such observations alone difficult. However, while B68 is embedded in a very hot, low-density medium, many starless cores are not, having interior/exterior density contrasts closer to unity. In this paper we investigate the oscillation damping rate as a function of the exterior density. For concreteness we use the same interior model employed in Broderick et al. (2007), with varying models for the exterior gas. We also develop a simple analytical formalism, based upon the linear perturbation analysis of the oscillations, which predicts the contribution to the damping rates due to the excitation of sound waves in the external medium. We find that the damping rate of oscillations on globules in dense molecular environments is always many periods, corresponding to hundreds of thousands of years, and persisting over the inferred lifetimes of the globules.

We place constraints on the dynamics of the Local Group (LG) by comparing the dipole of the Cosmic Microwave Background (CMB) with the peculiar velocity induced by the 2MRS galaxy sample. The analysis is limited by the lack of surveyed galaxies behind the Zone of Avoidance (ZoA). We therefore allow for a component of the LG velocity due to unknown mass concentrations behind the ZoA, as well as for an unknown transverse velocity of the Milky Way relative to the Andromeda galaxy. We infer extra motion along the direction of the Galactic center (where Galactic confusion and dust obscuration peaks) at the 95% significance level. With a future survey of the ZoA it might be possible to constrain the transverse velocity of the Milky Way relative to Andromeda.

As the luminosity of an accreting black hole drops to a few percent of Eddington, the spectrum switches from the familiar soft state to a hard state that is well-described by a distended and tenuous advection-dominated accretion flow (ADAF). An ADAF is a poor radiator, and the ion temperature can approach 10^{12} K near the center, although the electrons are cooler, with their temperature typically capped at ~10^{9-11} K. The foundational papers predicted that the large thermal energy in an ADAF would drive strong winds and jets, as later observed and also confirmed in computer simulations. Of chief interest, however, is the accreting gas that races inward. It carries the bulk of the accretion energy as stored thermal energy, which vanishes without a trace as the gas passes through the hole's event horizon. One thus expects black holes in the ADAF regime to be unusually faint. Indeed, this is confirmed by a comparison of accreting stellar-mass black holes and neutron stars, which reside in very similar transient X-ray binary systems. The black holes are on average observed to be fainter by a factor of ~100-1000. The natural explanation is that a neutron star must radiate the advected thermal energy from its surface, whereas a black hole can hide the energy behind its event horizon. The case for an event horizon in Sagittarius A*, which is immune to caveats on jet outflows and is furthermore independent of the ADAF model, is especially compelling. These two lines of evidence for event horizons are impervious to counterarguments that invoke strong gravity or exotic stars.

When we model black hole accretion sources such as active galactic nuclei and black hole X-ray binaries as advection-dominated accretion flows (ADAFs), it is neccesary to use the global solution to the equations rather than the simpler self-similar solution, since the latter is inaccurate in the region near the black hole where most of the radiation is emitted. However, technically, it is a difficult task to calculate the global solution because of the transonic nature of the flow, which makes it a two-point boundary value problem. In this paper we propose a simplified approach for calculating the global ADAF solution. We replace the radial momentum equation by a simple algebraic relation between the angular velocity of the gas and the Keplerian angular velocity, while keeping all other equations unchanged. It is then easy to solve the differential energy equations to obtain an approximate global solution. By adjusting the free parameters, we find that for almost any accretion rate and for $\alpha=0.1-0.3$ we can get good simplified global solutions. The predicted spectra from the approximate solutions are very close to the spectra obtained from the true global solutions.

We calculate the very high energy (sub-GeV to TeV) inverse Compton emission of GRB afterglows. We argue that this emission provides a powerful test of the currently accepted afterglow model. We focus on two processes: synchrotron self-Compton (SSC) emission within the afterglow blast wave, and external inverse Compton (EIC) emission which occurs when flare photons (produced by an internal process) pass through the blast wave. We show that if our current interpretations of the Swift XRT data are correct, there should be a canonical high energy afterglow emission light curve. Our predictions can be tested with high energy observatories such as GLAST, Whipple, H.E.S.S. and MAGIC. Under favorable conditions we expect afterglow detections in all these detectors.

We consider a simple Newtonian model of a steady accretion disk around a black hole. The model is based on height-integrated hydrodynamic equations, alpha-viscosity, and a pseudo-Newtonian potential that results in an innermost stable circular orbit (ISCO) that closely approximates the one predicted by GR. We find that the hydrodynamic models exhibit increasing deviations from the standard disk model of Shakura & Sunyaev as disk thickness H/R or the value of alpha increases. The latter is an analytical model in which the viscous torque is assumed to vanish at the ISCO. We consider the implications of the results for attempts to estimate black hole spin by using the standard disk model to fit continuum spectra of black hole accretion disks. We find that the error in the spin estimate is quite modest so long as H/R < 0.1 and alpha < 0.2. At worst the error in the estimated value of the spin parameter is 0.1 for a non-spinning black hole; the error is much less for a rapidly spinning hole. We also consider the density and disk thickness contrast between the gas in the disk and that inside the ISCO. The contrast needs to be large if black hole spin is to be successfully estimated by fitting the relativistically-broadened X-ray line profile of fluorescent iron emission from reflection off an accretion disk. In our hydrodynamic models, the contrast in density and thickness is low when H/R>0.1, sugesting that the iron line technique may be most reliable in extemely thin disks. We caution that these results have been obtained with a viscous hydrodynamic model and need to be confirmed with MHD simulations of radiatively cooled thin disks.

We have used the Novikov-Thorne thin disk model to fit the continuum X-ray spectra of three transient black hole X-ray binaries in the thermal state. From the fits we estimate the dimensionless spin parameters of the black holes to be: 4U 1543-47, a* = a/M = 0.7-0.85; GRO J1655-40, a* = 0.65-0.8; GRS 1915+105, a* = 0.98-1. We plan to expand the sample of spin estimates to about a dozen over the next several years. Some unresolved theoretical issues are briefly discussed.

In the absence of an active dynamo, purely poloidal magnetic field configurations are unstable to large-scale dynamical perturbations, and decay via reconnection on an Alfvenic timescale. Nevertheless, a number of classes of dynamo-free stars do exhibit significant, long-lived, surface magnetic fields. Numerical simulations suggest that the large-scale poloidal field in these systems is stabilized by a toroidal component of the field in the stellar interior. Using the principle of conservation of total helicity, we develop a variational principle for computing the structure of the magnetic field inside a conducting sphere surrounded by an insulating vacuum. We show that, for a fixed total helicity, the minimum energy state corresponds to a force-free configuration. We find a simple class of axisymmetric solutions, parametrized by angular and radial quantum numbers. However, these solutions have a discontinuity in the toroidal magnetic field at the stellar surface which will exert a toroidal stress on the surface of the star. We then describe two other classes of solutions, the standard spheromak solutions and ones with fixed surface magnetic fields, the latter being relevant for neutron stars with rigid crusts. We discuss the implications of our results for the structure of neutron star magnetic fields, the decay of fields, and the origin of variability and outbursts in magnetars.

This Supplementary Information provides details about the spectral extraction (crowding issues and the removal of nebular lines), a discussion about the distance to M33, a model for the O-star wind and the measurement of the true photospheric X-ray eclipse width Theta, and details about ellipsoidal modelling. It also contains five related figures, one related table, and additional references.

Stellar-mass black holes are discovered in X-ray emitting binary systems, where their mass can be determined from the dynamics of their companion stars. Models of stellar evolution have difficulty producing black holes in close binaries with masses >10 solar masses, which is consistent with the fact that the most massive stellar black holes known so all have masses within 1 sigma of 10 solar masses. Here we report a mass of 15.65 +/- 1.45 solar masses for the black hole in the recently discovered system M33 X-7, which is located in the nearby galaxy Messier 33 (M33) and is the only known black hole that is in an eclipsing binary. In order to produce such a massive black hole, the progenitor star must have retained much of its outer envelope until after helium fusion in the core was completed. On the other hand, in order for the black hole to be in its present 3.45 day orbit about its 70.0 +/- 6.9 solar mass companion, there must have been a ``common envelope'' phase of evolution in which a significant amount of mass was lost from the system. We find the common envelope phase could not have occured in M33 X-7 unless the amount of mass lost from the progenitor during its evolution was an order of magnitude less than what is usually assumed in evolutionary models of massive stars.

In a previous paper, we modeled the oscillations of a thermally-supported (Bonnor-Ebert) sphere as non-radial, linear perturbations following a standard analysis developed for stellar pulsations. The predicted column density variations and molecular spectral line profiles are similar to those observed in the Bok globule B68 suggesting that the motions in some starless cores may be oscillating perturbations on a thermally supported equilibrium structure. However, the linear analysis is unable to address several questions, among them the stability, and lifetime of the perturbations. In this paper we simulate the oscillations using a three-dimensional numerical hydrodynamic code. We find that the oscillations are damped predominantly by non-linear mode-coupling, and the damping time scale is typically many oscillation periods, corresponding to a few million years, and persisting over the inferred lifetime of gobules.

We consider photon emission from a supercritical accretion disk in which photons in the inner regions are trapped and advected in towards the center. Such disks are believed to be present in failed supernovae and gamma-ray bursts, and possibly also in ultra-luminous X-ray sources. We show that the luminosity from a supercritical accretion disk is proportional to the logarithm of the mass accretion rate when the vertical profile of the matter density is exponential. We obtain analytical results also for other density profiles, but the dependence of the luminosity on the accretion rate is no longer logarithmic.

We describe a program that we have embarked on to estimate the spins of stellar-mass black holes in X-ray binaries. We fit the continuum X-ray spectrum of the radiation from the accretion disk using the standard thin disk model, and extract the dimensionless spin parameter a* = a/M of the black hole as a parameter of the fit. We have obtained results on three systems, 4U 1543-47 (a* = 0.7-0.85), GRO J1655-40 (0.65-0.8), and GRS 1915+105 (0.98-1), and have nearly completed analysis of two additional systems. We anticipate expanding the sample of spin estimates to about a dozen over the next several years.

The Rossi X-ray Timing Explorer has detected nearly coherent oscillations in the tails of type I X-ray bursts from 17 low-mass X-ray binaries. The oscillations are thought to be generated by brightness fluctuations associated with a surface mode on the rotating neutron star. The mechanism that drives the modes is, however, not understood, since the burning layer is stable to thermal perturbations. We show here via a linear perturbation analysis that, even under conditions when pure thermal perturbations are stable, nonradial surface modes may still be unstable by the epsilon mechanism. Specifically, we find that, if helium-burning reactions supply a reasonable fraction of the outgoing flux during burst decay, nonradial surface modes will grow in time. On the other hand, the same modes are likely to be stable in the presence of hydrogen burning via the rp-process. The results naturally explain why oscillations in the decay phase of type I X-ray bursts are detected only from short-duration bursts.

Active galactic nuclei, x-ray binaries, pulsars, and gamma-ray bursts are all believed to be powered by compact objects surrounded by relativistic plasma flows driving phenomena such as accretion, winds, and jets. These flows are often accurately modelled by the relativistic magnetohydrodynamics (MHD) approximation. Time-dependent numerical MHD simulations have proven to be especially insightful, but one regime that remains difficult to simulate is when the energy scales (kinetic, thermal, magnetic) within the plasma become disparate. We develop a numerical scheme that significantly improves the accuracy and robustness of the solution in this regime. We use a modified form of the WENO method to construct a finite-volume general relativistic hydrodynamics code called WHAM that converges at fifth order. We avoid (1) field-by-field decomposition by adaptively reducing down to 2-point stencils near discontinuities for a more accurate treatment of shocks, and (2) excessive reduction to low order stencils, as in the standard WENO formalism, by maintaining high order accuracy in smooth monotonic flows. Our scheme performs the proper surface integral of the fluxes, converts cell averaged conserved quantities to point conserved quantities before performing the reconstruction step, and correctly averages all source terms. We demonstrate that the scheme is robust in strong shocks, very accurate in smooth flows, and maintains accuracy even when the energy scales in the flow are highly disparate.

We use the two-zone model of Cooper & Narayan to study the onset and time evolution of hydrogen-triggered type I X-ray bursts on accreting neutron stars. At the lowest accretion rates, thermally unstable hydrogen burning ignites helium as well and produces a mixed hydrogen and helium burst. For somewhat higher accretion rates, thermally unstable hydrogen burning does not ignite helium and thus triggers only a weak hydrogen flash. The peak luminosities of weak hydrogen flashes are typically much lower than the accretion luminosity. These results are in accord with previous theoretical work. We find that a series of weak hydrogen flashes generates a massive layer of helium that eventually ignites in an energetic pure helium flash. Although previously conjectured, this is the first time such bursting behavior has been actually demonstrated in a theoretical model. For yet higher accretion rates, hydrogen burning is thermally stable and thus steadily generates a layer of helium that ultimately ignites in a pure helium flash. We find that, for a narrow range of accretion rates between the mixed hydrogen and helium burst and weak hydrogen flash regimes, unstable hydrogen burning ignites helium only after a short series of weak hydrogen flashes has generated a sufficiently deep layer of helium. These bursts have fluences that are intermediate between those of normal mixed hydrogen and helium bursts and energetic pure helium flashes.

The gravastar model, which postulates a strongly correlated thin shell of anisotropic matter surrounding a region of anti-de Sitter space, has been proposed as an alternative to black holes. We discuss constraints that present-day observations of well-known black hole candidates place on this model. We focus upon two black hole candidates known to have extraordinarily low luminosities: the supermassive black hole in the Galactic Center, Sagittarius A*, and the stellar-mass black hole, XTE J1118+480. We find that the length scale for modifications of the type discussed in Chapline et al. (2003) must be sub-Planckian.

We investigate the latitude at which type I X-ray bursts are ignited on rapidly rotating accreting neutron stars. We find that, for a wide range of accretion rates, ignition occurs preferentially at the equator, in accord with the work of Spitkovsky et al. However, for a range of accretion rates below the critical rate above which bursts cease, ignition occurs preferentially at higher latitudes. The range of accretion rates over which nonequatorial ignition occurs is an increasing function of the neutron star spin frequency. These findings have significant implications for thermonuclear flame propagation, and they may explain why oscillations during the burst rise are detected predominantly when the accretion rate is high. They also support the suggestion of Bhattacharyya & Strohmayer that non-photospheric radius expansion double-peaked bursts and the unusual harmonic content of oscillations during the rise of some bursts result from ignition at or near a rotational pole.

In paper I, we showed that time-dependent general relativistic magnetohydrodynamic (GRMHD) numerical models of accretion disks, although being highly turbulent, have surprisingly simple electromagnetic properties. In particular, the toroidal current density in the disk takes the form $dI_\phi/dr \propto r^{-5/4}$. Guided by this simplicity, we use a time-dependent general relativistic force-free electrodynamics (GRFFE) code to study an idealized problem in which the accretion disk is replaced by an infinitely thin rotating equatorial current sheet. We consider both an $r^{-5/4}$ current profile and an $r^{-1}$ profile, the latter corresponding to the paraboloidal model of Blandford & Znajek (1977). The force-free magnetosphere we obtain with the $r^{-5/4}$ current sheet matches remarkably well to the Poynting-dominated jet seen in GRMHD numerical models. By comparing to the non-rotating force-free model studied in paper I, rotation is seen to lead to mild decollimation of the jet suggesting that hoop-stress forces nearly cancel centrifugal forces. In order to study the process that generates the corona and disk wind and destroys the ordered field in the corona in GRMHD numerical models, the force-free field with the $r^{-5/4}$ current distribution is embedded in an accretion disk and followed in a GRMHD simulation. Reconnection and magnetic stresses contribute to a magnetized, thermal wind without the aid of an ordered field threading the disk.

General relativistic numerical simulations of magnetized accretion flows around black holes show a disordered electromagnetic structure in the disk and corona and a highly relativistic, Poynting-dominated funnel jet in the polar regions. The polar jet is nearly consistent with the stationary paraboloidal Blandford-Znajek model of an organized field threading the polar regions of a rotating black hole. How can a disordered accretion disk and corona lead to an ordered jet? We show that the polar jet is associated with a strikingly simple angular-integrated toroidal current distribution $dI_\phi/dr \propto r^{-5/4}$, where $I_\phi(r)$ is the toroidal current enclosed inside radius $r$. We demonstrate that the poloidal magnetic field in the simulated jet agrees well with the force-free field solution for a non-rotating thin disk with an $r^{-5/4}$ toroidal current, suggesting rotation leads to negligible self-collimation. We find that the polar field is confined/collimated by the corona. The electromagnetic field in the disk also scales as $r^{-5/4}$, which is consistent with some Newtonian accretion models that assume rough equipartition between magnetic and gas pressure. However, the agreement is accidental since toward the black hole the magnetic pressure increases faster than the gas pressure. This field dominance near the black hole is associated with magnetic stresses that imply a large effective viscosity parameter $\alpha\sim 1$, whereas the typically assumed value of $\alpha\sim 0.1$ holds far from the black hole.[abridged]

We consider a self-similar force-free wind flowing out of an infinitely thin disk located in the equatorial plane. On the disk plane, we assume that the magnetic stream function $P$ scales as $P\propto R^\nu$, where $R$ is the cylindrical radius. We also assume that the azimuthal velocity in the disk is constant: $v_\phi = Mc$, where $M<1$ is a constant. For each choice of the parameters $\nu$ and $M$, we find an infinite number of solutions that are physically well-behaved and have fluid velocity $\leq c$ throughout the domain of interest. Among these solutions, we show via physical arguments and time-dependent numerical simulations that the minimum-torque solution, i.e., the solution with the smallest amount of toroidal field, is the one picked by a real system. For $\nu \geq 1$, the Lorentz factor of the outflow increases along a field line as $\gamma \approx M(z/\Rfp)^{(2-\nu)/2} \approx R/R_{\rm A}$, where $\Rfp$ is the radius of the foot-point of the field line on the disk and $R_{\rm A}=\Rfp/M$ is the cylindrical radius at which the field line crosses the Alfven surface or the light cylinder. For $\nu < 1$, the Lorentz factor follows the same scaling for $z/\Rfp < M^{-1/(1-\nu)}$, but at larger distances it grows more slowly: $\gamma \approx (z/\Rfp)^{\nu/2}$. For either regime of $\nu$, the dependence of $\gamma$ on $M$ shows that the rotation of the disk plays a strong role in jet acceleration. On the other hand, the poloidal shape of a field line is given by $z/\Rfp \approx (R/\Rfp)^{2/(2-\nu)}$ and is independent of $M$. Thus rotation has neither a collimating nor a decollimating effect on field lines, suggesting that relativistic astrophysical jets are not collimated by the rotational winding up of the magnetic field.

We consider the distribution of local supermassive black hole Eddington ratios and accretion rates, accounting for the dependence of radiative efficiency and bolometric corrections on the accretion rate. We find that black hole mass growth, both of the integrated mass density and the masses of most individual objects, must be dominated by an earlier, radiatively efficient, high accretion rate stage, and not by the radiatively inefficient low accretion rate phase in which most local supermassive black holes are currently observed. This conclusion is particularly true of supermassive black holes in elliptical host galaxies, as expected if they have undergone merger activity in the past which would fuel quasar activity and rapid growth. We discuss models of the time evolution of accretion rates and show that they all predict significant mass growth in a prior radiatively efficient state. The only way to avoid this conclusion is through careful fine-tuning of the accretion/quasar timescale to a value that is inconsistent with observations. Our results agree with a wide range of observational inferences drawn from the quasar luminosity function and X-ray background synthesis models, but our approach has the virtue of being independent of the modeling of source populations. Models in which black holes spend the great majority of their time in low accretion rate phases are thus completely consistent both with observations implying mass gain in relatively short, high accretion rate phases and with the local distribution of accretion rates.

Based on a spectral analysis of the X-ray continuum that employs a fully relativistic accretion-disk model, we conclude that the compact primary of the binary X-ray source GRS 1915+105 is a rapidly-rotating Kerr black hole. We find a lower limit on the dimensionless spin parameter of a* greater than 0.98. Our result is robust in the sense that it is independent of the details of the data analysis and insensitive to the uncertainties in the mass and distance of the black hole. Furthermore, our accretion-disk model includes an advanced treatment of spectral hardening. Our data selection relies on a rigorous and quantitative definition of the thermal state of black hole binaries, which we used to screen all of the available RXTE and ASCA data for the thermal state of GRS 1915+105. In addition, we focus on those data for which the accretion disk luminosity is less than 30% of the Eddington luminosity. We argue that these low-luminosity data are most appropriate for the thin alpha-disk model that we employ. We assume that there is zero torque at the inner edge of the disk, as is likely when the disk is thin, although we show that the presence of a significant torque does not affect our results. Our model and the model of the relativistic jets observed for this source constrain the distance and black hole mass and could thus be tested by determining a VLBA parallax distance and improving the measurement of the mass function. Finally, we comment on the significance of our results for relativistic-jet and core-collapse models, and for the detection of gravitational waves.

We investigate the effect of the hot CNO cycle breakout reaction 15O(alpha,gamma)19Ne on the occurrence of type I X-ray bursts on accreting neutron stars. For f_rp <~ 0.1, where f_rp is a dimensionless factor by which we multiply the 15O(alpha,gamma)19Ne reaction rate of Caughlan & Fowler (1988), our model predicts that bursts should occur only for accretion rates below a critical value of approximately 0.3 times the Eddington limit. This agrees with observations. For larger values of f_rp, including the standard choice f_rp = 1, the model switches to a new regime in which bursts occur all the way up to roughly the Eddington limit. Since the latter regime disagrees with observations, we suggest that the true 15O(alpha,gamma)19Ne reaction rate is lower than usually assumed.

We construct a two-zone model to describe H and He burning on the surface of an accreting neutron star and use it to study the triggering of type I X-ray bursts. Although highly simplified, the model reproduces all of the bursting regimes seen in the more complete global linear stability analysis of Narayan & Heyl (2003), including the regime of delayed mixed bursts. The results are also consistent with observations of type I X-ray bursts. At accretion rates Mdot < 0.1 Mdot_Edd, thermonuclear He burning via the well-known thin-shell thermal instability triggers bursts. As Mdot increases, however, the trigger mechanism evolves from the fast thermal instability to a slowly growing overstability involving both H and He burning. The competition between nuclear heating via the beta-limited CNO cycle and the triple-alpha process on the one hand, and radiative cooling via photon diffusion and emission on the other hand, drives oscillations with a period approximately equal to the H-burning timescale. If these oscillations grow, the gradually rising temperature at the base of the helium layer eventually provokes a thin-shell thermal instability and hence a delayed mixed burst. For Mdot > 0.25 Mdot_Edd, there is no instability or overstability, and there are no bursts. Nearly all other theoretical models predict that bursts should occur for all Mdot < Mdot_Edd, in conflict with both our results and observations. We suggest that this discrepancy arises from the assumed strength of the hot CNO cycle breakout reaction 15O(alpha,gamma)19Ne in these other models. That observations agree much better with the results of Narayan & Heyl and our two-zone model, both of which neglect breakout reactions, may imply that the true 15O(alpha,gamma)19Ne cross section is much smaller than assumed in previous investigations.

We measure the rates of type I X-ray bursts, as a function of the bolometric luminosity, from a likely complete sample of 37 non-pulsing transients (1996-2004). Our goals are to test the burst model for neutron stars and to investigate whether black holes have event horizons. We find 135 type I bursts in 3.7 Ms of exposure for the neutron-star group, and the burst rate function is generally consistent with model predictions. However, for the black hole groups (18 sources), there are no confirmed type I bursts in 6.5 Ms of exposure, and the upper limits in the burst function are inconsistent with the model predictions for heavy compact objects with a solid surface. There are systematic spectral differences between the neutron-star and black-hole groups, supporting the presumption that physical differences underly the sample classifications. These results provide indirect evidence that black holes do have event horizons.

We consider a model in which Sgr A*, the 3.5x10^6 M_sun supermassive black hole candidate at the Galactic Center, is a compact object with a surface. Given the very low quiescent luminosity of Sgr A* in the near infrared, the existence of a hard surface, even in the limit in which the radius approaches the horizon, places severe constraints upon the steady mass accretion rate in the source, requiring dM/dt < 10^-12 M_sun/yr. This limit is well below the minimum accretion rate needed to power the observed submillimeter luminosity of Sgr A*. We thus argue that Sgr A* does not have a surface, i.e., it must have an event horizon. The argument could be made more restrictive by an order of magnitude with microarcsecond resolution imaging, e.g., with submillimeter VLBI.

Current X-ray missions are providing high-quality X-ray spectra from neutron stars (NSs) in quiescent low-mass X-ray binaries (qLMXBs). This has motivated us to calculate new hydrogen-atmosphere models, including opacity due to free-free absorption and Thomson scattering, thermal electron conduction, and self-irradiation by photons from the compact object. We have constructed a self-consistent grid of neutron star models covering a wide range of surface gravities as well as effective temperatures, which we make available to the scientific community. We present multi-epoch Chandra X-ray observations of the qLMXB X7 in the globular cluster 47 Tuc, which is remarkably nonvariable on timescales from minutes to years. Its high-quality X-ray spectrum is adequately fit by our hydrogen-atmosphere model without any hard power-law component or narrow spectral features. If a mass of 1.4 Msol is assumed, our spectral fits require that its radius be in the range R=14.5^{+1.8}_{-1.6} km (90% confidence), larger than expected from currently preferred models of NS interiors. If its radius is assumed to be 10 km, then a mass of M=2.20^{+0.03}_{-0.16} Msol is required. Using models with the appropriate surface gravity for each value of the mass and radius becomes important for interpretation of the highest quality data.

If the split, asymmetric molecular spectral line profiles that are seen in many starless cores are interpreted as indicative of global collapse or expansion of the core then one possible implication is that most starless cores have short lifetimes on the order of the collapse or sound crossing time scale. An alternative interpretation of the line profiles as indicative of perturbations on an underlying equilibrium structure leads to the opposite implication, that many cores have long lifetimes. While evidence suggests that some cores are collapsing on a free-fall time scale, we show that observations of some other starless cores can be reproduced by a model of non-radial oscillations about the equilibrium configuration of a pressure-bounded, thermally-supported sphere (Bonnor-Ebert sphere). We model the oscillations as linear perturbations following a standard analysis developed for stellar pulsations and compare the column densities and molecular spectral line profiles predicted from a particular model to observations of the Bok globule B68.

(abridged) We have investigated the physical conditions under which accreting neutron stars can both produce and preserve sufficient quantities of carbon fuel to trigger superbursts. Our models span the plausible ranges of neutron star thermal conductivities, core neutrino emission mechanisms, and areal radii, as well as the CNO abundances in the accreted material. We find that neutron stars that accrete hydrogen-rich material with CNO mass fractions <~ that of the Sun will not exhibit superbursts under any circumstances. Neutron stars that accrete material with CNO mass fractions >~ 4 times that of the Sun will exhibit superbursts at accretion rates in the observed range. On this basis, we suggest that the mass donors of superburst systems must have enhanced CNO abundances. The accreted CNO acts only as a catalyst for hydrogen burning via the hot CNO cycle, and therefore it is only the sum of the three elements' mass fractions that is important. Systems that exhibit superbursts are observed to differ from those that do not exhibit superbursts in the nature of their helium-triggered Type I X-ray bursts: the bursts have shorter durations and much greater alpha-values. Increasing the CNO abundance of the accreted material in our models reproduces both of these observations. Many compact binary systems have been observed in which the abundances of the accreting material are distinctly non-solar. Though abundance analyses of the systems that exhibit superbursts currently do not exist, Bowen fluorescence blend profiles of 4U 1636-536 and Ser X-1 suggest that the mass donor stars may indeed have non-solar CNO metallicities. More detailed abundance analyses of the accreting matter in systems that exhibit superbursts are needed to verify our assertion that the matter is rich in CNO elements.

We fit X-ray spectral data in the thermal dominant or high soft state of two dynamically confirmed black holes, GRO J1655-40 and 4U 1543-47, and estimate the dimensionless spin parameters a* = a/M of the two holes. For GRO J1655-40, using a spectral hardening factor computed for a non-LTE relativistic accretion disk, we estimate a* ~ 0.75 and ~ 0.65-0.75, respectively, from ASCA and RXTE data. For 4U 1543-47, we estimate a* ~ 0.75-0.85 from RXTE data. Thus, neither black hole has a spin approaching the theoretical maximum a* = 1.

We model thermal X-ray emission from the accreting supermassive black hole Sagittarius A* at the Galactic Center. For the region inside $1.^{\prime\prime}5$ of the center, we use a {generalized} radiatively inefficient accretion flow (RIAF) model, and for the region out to $10^{\prime\prime}$ we use published models of the ambient hot gas. We calculate the equivalent widths of Hydrogen-like and Helium-like emission lines of various elements, notably iron. We predict that a strong Helium-like iron line with an equivalent width $\sim1$ keV should be emitted by both the external medium and the RIAF. The equivalent width in the external medium is sensitive to the metallicity $Z$ of the gas as well as the mean temperature. For reasonable choices of these parameters, the calculated results agree with Chandra's detection of an iron line with an equivalent width of 1.3 keV within $10^{\prime\prime}$. The emission from within $1.^{\prime\prime}5$ is not sensitive to the external temperature, but is sensitive to the density and, especially, temperature profile inside the Bondi radius. For the range of profiles we consider, we calculate the equivalent width of the iron line to be $\sim0.6-1.5 (Z/Z_\odot)$ keV, where $Z_\odot$ is the solar metallicity. We present a new Chandra spectrum of the quiescent emission within $1.^{\prime\prime}5$ of Sgr A*. The measured equivalent width of the iron line is 0.7 keV. Although this measurement has a large uncertainty, it is consistent with our predictions, provided the metallicity of the gas is approximately solar.

We present a fitting function to describe the statistics of flux modulations caused by interstellar scintillation. The function models a very general quantity: the cross-correlation of the flux observed from a compact radio source of finite angular size observed at two frequencies and at two positions or times. The formula will be useful for fitting data from sources such as intra-day variables and gamma-ray burst afterglows. These sources are often observed at relatively high frequencies (several gigahertz) where interstellar scattering is neither very strong nor very weak, so that asymptotic formulae are inapplicable.

We study the possible origin of hydrodynamic turbulence in cold accretion disks such as those in star-forming systems and quiescent cataclysmic variables. As these systems are expected to have neutral gas, the turbulent viscosity is likely to be hydrodynamic in origin, not magnetohydrodynamic. Therefore MRI will be sluggish or even absent in such disks. Although there are no exponentially growing eigenmodes in a hydrodynamic disk, because of the non-normal nature of the eigenmodes, a large transient growth in the energy is still possible, which may enable the system to switch to a turbulent state. For a Keplerian disk, we estimate that the energy will grow by a factor of 1000 for a Reynolds number close to a million.

This article reviews the current status of black hole astrophysics, focusing on topics of interest to a physics audience. Astronomers have discovered dozens of compact objects with masses greater than 3 solar masses, the likely maximum mass of a neutron star. These objects are identified as black hole candidates. Some of the candidates have masses of 5 to 20 solar masses and are found in X-ray binaries, while the rest have masses from a million to a billion solar masses and are found in galactic nuclei. A variety of methods are being tried to estimate the spin parameters of the candidate black holes. There is strong circumstantial evidence that many of the objects have event horizons. Recent MHD simulations of magnetized plasma accreting on rotating black holes seem to hint that relativistic jets may be produced by a magnetic analog of the Penrose process.

Cold accretion disks such as those in star-forming systems, quiescent cataclysmic variables, and some active galactic nuclei, are expected to have neutral gas which does not couple well to magnetic fields. The turbulent viscosity in such disks must be hydrodynamic in origin, not magnetohydrodynamic. We investigate the growth of hydrodynamic perturbations in a linear shear flow sandwiched between two parallel walls. The unperturbed flow is similar to plane Couette flow but with a Coriolis force included. Although there are no exponentially growing eigenmodes in this system, nevertheless, because of the non-normal nature of the eigenmodes, it is possible to have a large transient growth in the energy of perturbations. For a constant angular momentum disk, we find that the perturbation with maximum growth has a wave-vector in the vertical direction. The energy grows by more than a factor of 100 for a Reynolds number R=300 and more than a factor of 1000 for R=1000. Turbulence can be easily excited in such a disk, as found in previous numerical simulations. For a Keplerian disk, on the other hand, similar vertical perturbations grow by no more than a factor of 4, explaining why the same simulations did not find turbulence in this system. However, certain other two-dimensional perturbations with no vertical structure do exhibit modest growth. For the optimum two-dimensional perturbation, the energy grows by a factor of ~100 for R~10^4.5 and by a factor of 1000 for R~10^6. It is conceivable that these two-dimensional disturbances might lead to self-sustained turbulence. The Reynolds numbers of cold astrophysical disks are much larger even than 10^6, therefore, hydrodynamic turbulence may be possible in disks.

Despite observational evidence for cold neutral astrophysical accretion disks, the viscous process which may drive the accretion in such systems is not yet understood. While molecular viscosity is too small to explain the observed accretion efficiencies by more than ten orders of magnitude, the absence of any linear instability in Keplerian accretion flows is often used to rule out the possibility of turbulent viscosity. Recently, the fact that some fine tuned disturbances of any inviscid shear flow can reach arbitrarily large transient growth has been proposed as an alternative route to turbulence in these systems. We present an analytic study of this process for 3D plane wave disturbances of a general rotating shear flow in Lagrangian coordinates, and demonstrate that large transient growth is the generic feature of non-axisymmetric disturbances with near radial leading wave vectors. The maximum energy growth is slower than quadratic, but faster than linear in time. The fastest growth occurs for two dimensional perturbations, and is only limited by viscosity, and ultimately by the disk vertical thickness. After including viscosity and vertical structure, we find that, as a function of the Reynolds number, R, the maximum energy growth is approximately 0.4 (R/log R)^{2/3}, and put forth a heuristic argument for why R > 10^4 is required to sustain turbulence in Keplerian disks. Therefore, assuming that there exists a non-linear feedback process to replenish the seeds for transient growth, astrophysical accretion disks must be well within the turbulent regime. However, large 3D numerical simulations running for many orbital times, and/or with fine tuned initial conditions, are required to confirm Keplerian hydrodynamic turbulence on the computer.

We suggest that part of the infalling material during the core-collapse of a massive star goes into orbit around the compact core to form a hot, dense, centrifugally-supported accretion disk whose evolution is strongly influenced by neutrino interactions. Under a wide range of conditions, this neutrino-dominated accretion flow will be advection-dominated and will develop a substantial outflowing wind. We estimate the energy carried out in the wind and find that it exceeds $10^{50}$ erg for a wide range of parameters and even exceeds 10^{51} erg for reasonable parameter choices. We propose that the wind energy will revive a stalled shock and will help produce a successful supernova explosion. We discuss the role of the disk wind in both prompt and delayed explosions. While both scenarios are feasible, we suggest that a delayed explosion is more likely, and perhaps even unavoidable. Finally, we suggest that the disk wind may be a natural site for r-process nucleosynthesis.

We carry out a general-relativistic global linear stability analysis of the amassed carbon fuel on the surface of an accreting neutron star to determine the conditions under which superbursts occur. We reproduce the general observational characteristics of superbursts, including burst fluences, recurrence times, and the absence of superbursts on stars with accretion rates below 10% of the Eddington limit. By comparing our results with observations, we are able to set constraints on neutron star parameters such as the stellar radius and neutrino cooling mechanism in the core. Specifically, we find that accreting neutron stars with ordered crusts and highly efficient neutrino emission in their cores (due to direct URCA or pionic reactions, for example) produce extremely energetic (> 10^44 ergs) superbursts which are inconsistent with observations, in agreement with previous investigations. Also, because of pycnonuclear burning of carbon, they do not have superbursts in the range of accretion rates at which superbursts are actually observed unless the crust is very impure. Stars with less efficient neutrino emission (due to modified URCA reactions, for example) produce bursts that agree better with observations. Stars with highly inefficient neutrino emission in their cores produce bursts that agree best with observations. All systems that accrete primarily hydrogen and in which superbursts are observed show evidence of H- and He-burning delayed mixed bursts. We speculate that delayed mixed bursts provide sufficient amounts of carbon fuel for superbursts and are thus a prerequisite for having superbursts. We compare our global stability analysis to approximate one-zone criteria used by other authors and identify a particular set of approximations that give accurate results for most choices of parameters. (abstract truncated)

Fallback disks are expected to form around new-born neutron stars following a supernova explosion. In almost all cases, the disk will pass through a propeller stage. If the neutron star is spinning rapidly (initial period $\sim 10$ ms) and has an ordinary magnetic moment ($\sim 10^{30}$ G cm$^3$), the rotational power transferred to the disk by the magnetic field of the neutron star will exceed the Eddington limit by many orders of magnitude, and the disk will be rapidly disrupted. Fallback disks can thus survive only around slow-born neutron stars and around black holes, assuming the latter do not torque their surrounding disks as strongly as do neutron stars. This might explain the apparent rarity of fallback disks around young compact objects.

Turbulent viscosity in cold accretion disks is likely to be hydrodynamic in origin. We investigate the growth of hydrodynamic perturbations in a small region of a disk, which we model as a linear shear flow with Coriolis force, between two parallel walls. Although there are no exponentially growing eigenmodes in this system, because of the non-normal nature of the modes, it is possible to have a large transient growth in the energy of certain perturbations. For a constant angular momentum disk, the energy grows by more than a factor of 1000 for a Reynolds number of only 1000, and so turbulence is easily excited. For a Keplerian disk, the growth is more modest, and energy growth by a factor of 1000 requires a Reynolds number of nearly a million. Accretion disks have even larger Reynolds numbers than this. Therefore, transient growth of perturbations could seed turbulence in such disks.

We investigate the dependence of observational properties of black hole X-ray binaries on the inclination angle i of their orbits. We find the following: (1) Transient black hole binaries show no trend in their quiescent X-ray luminosities as a function of i, suggesting that the radiation is not significantly beamed. This is consistent with emission from an accretion disk. If the X-rays are from a jet, then the Lorentz factor gamma of the jet is less than 1.24 at the 90% confidence level. (2) The X-ray binary 4U1543-47 with i of order 21 degrees has a surprisingly strong fluorescent iron line in the high soft state. Quantifying an earlier argument by Park et al. (2004), we conclude that if the continuum X-ray emission in this source is from a jet, then gamma < 1.04. (3) None of the known binaries has cos i < 0.25 or i > 75 degrees. This fact, plus the lack of eclipses among the 20 black hole binaries in our sample, strongly suggests at the 99.5% confidence level that systems with large inclination angles are hidden from view. The obscuration could be the result of disk flaring, as suggested by Milgrom (1978) for neutron star X-ray binaries. (4) Transient black hole binaries with i ~ 70-75 degrees have significantly more complex X-ray light curves than systems with i < 65 degrees. This may be the result of variable obscuration and/or variable height above the disk of the radiating gas.

We use a ray-tracing technique to compute the observed spectrum of a thin accretion disk around a Kerr black hole. We include all relativistic effects such as frame-dragging, Doppler boost, gravitational redshift, and bending of light by the gravity of the black hole. We also include self-irradiation of the disk as a result of light deflection. Assuming that the disk emission is locally blackbody, we show how the observed spectrum depends on the spin of the black hole, the inclination of the disk, and the torque at the inner edge of the disk. We find that the effect of a nonzero torque on the spectrum can, to a good approximation, be absorbed into a zero-torque model by adjusting the mass accretion rate and the normalization. We describe a computer model, called KERRBB, which we have developed for fitting the spectra of black hole X-ray binaries. Using KERRBB within the X-ray data reduction package XSPEC, and assuming a spectral hardening factor f_col = 1.7, we analyze the spectra of three black hole X-ray binaries: 4U1543-47, XTE J1550-564, and GRO J1655-40. We estimate the spin parameters of the black holes in 4U1543-47 and GRO J1655-40 to be a/M ~ 0.6 and ~ 0.6-0.7, respectively. If f_col ~ 1.5-1.6, as in a recent study, then we find a/M ~ 0.7-0.8 and ~ 0.8-0.9, respectively. These estimates are subject to additional uncertainties in the assumed black hole masses, distances and disk inclinations.

Multi-wavelength observations of the black hole X-ray binary XTE J1118+480 have offered abundant spectral and timing information about the source, and have thus provided serious challenges to theoretical models. We propose a coupled accretion-jet model to interpret the observations. We model the accretion flow as an outer standard thin accretion disk truncated at a transition radius by an inner hot accretion flow. The accretion flow accounts for the observed UV and X-ray emission, but it substantially under-predicts the radio and infrared fluxes, even after we allow for nonthermal electrons in the hot flow. We attribute the latter components to a jet. We model the jet emission by means of the internal shock scenario which is widely employed for gamma-ray bursts. In our accretion-jet model of XTE J1118+480, the jet dominates the radio and infrared emission, the thin disk dominates the UV emission, and the hot flow produces most of the X-ray emission. The optical emission has contributions from all three components: jet, thin disk, and hot flow. The model qualitatively accounts for timing features, such as the intriguing positive and negative time lags between the optical and X-ray emission, and the wavelength-dependent variability amplitude.

At luminosities below a few percent of Eddington, accreting black holes switch to a hard spectral state which is very different from the soft blackbody-like spectral state that is found at higher luminosities. The hard state is well-described by a two-temperature, optically thin, geometrically thick, advection-dominated accretion flow (ADAF) in which the ions are extremely hot (up to $10^{12}$ K near the black hole), the electrons are also hot ($\sim10^{9-10.5}$ K), and thermal Comptonization dominates the X-ray emission. The radiative efficiency of an ADAF decreases rapidly with decreasing mass accretion rate, becoming extremely low when a source reaches quiescence. ADAFs are expected to have strong outflows, which may explain why relativistic jets are often inferred from the radio emission of these sources. It has been suggested that most of the X-ray emission also comes from a jet, but this is less well established.

The spectra of many X-ray pulsars show, in addition to a power law, a low-energy component that has often been modeled as a blackbody with kT ~ 0.1 keV. However the physical origin of this soft excess has remained a mystery. We examine a sample of well-studied, bright X-ray pulsars, which have been observed using ROSAT, ASCA, Ginga, RXTE, BeppoSAX, Chandra, and XMM-Newton. In particular we consider the Magellanic Cloud pulsars SMC X-1, LMC X-4, XTE J0111.2-7317, and RX J0059.2-7138 and the Galactic sources Her X-1, 4U 1626-67, Cen X-3, and Vela X-1. We show that the soft excess is a very common if not ubiquitous feature intrinsic to X-ray pulsars. We evaluate several possible mechanisms for the soft emission, using theoretical arguments as well as observational clues such as spectral shapes, eclipses, pulsations of the soft component, and superorbital modulation of the source flux. We find that reprocessing of hard X-rays from the neutron star by the inner region of the accretion disk is the only process that can explain the soft excess in all the pulsars with Lx > 10^38 ergs/s. Other mechanisms, such as emission from diffuse gas in the system, are important in less luminous objects.

Detection of Jupiter mass companions to nearby solar type stars with precise radial velocity measurements is now routine, and Doppler surveys are moving towards lower velocity amplitudes. The detection of several Neptune-mass planets with orbital periods less than a week has been reported. The drive toward the search for close-in Earth-mass planets is on the agenda. Successful detection or meaningful upper limits will place important constraints on the process of planet formation. In this paper, we quantify the statistics of detection of low-mass planets in-close orbits, showing how the detection threshold depends on the number and timing of the observations. In particular, we consider the case of a low-mass planet close to but not on the 2:1 mean motion resonance with a hot jupiter. This scenario is a likely product of the core-accretion hypothesis for planet formation coupled with migration of jupiters in the protoplanetary disk. It is also advantageous for detection because the orbital period is well-constrained. Detection of few Earth mass rocky cores will require ~ 1 m/s velocity precision, and most importantly, a much better understanding of stellar radial velocity jitter.

The standard spectral model for analyzing the soft component of thermal emission from a thin accretion disk around a black hole is the multi-temperature blackbody (MTB) model. The widely used implementation of this model, which is known as ``diskbb,'' assumes nonzero torque at the inner edge of the accretion disk. This assumption is contrary to the classic and current literature on thin-disk accretion, which advocates the use of a zero-torque boundary condition. Consequently, we have written code for a zero-torque model, ``ezdiskbb,'' which we compare to the nonzero-torque model diskbb by fitting RXTE spectra of three well-known black hole binaries: 4U 1543-47, XTE J1550-564, and GRO J1655-40. The chief difference we find is that the zero-torque model gives a value for the inner disk radius that is about 2.2 times smaller than the value given by diskbb. This result has important implications, especially for the determination of black-hole angular momentum and mass accretion rate.

A soft component of thermal emission is very commonly observed from the surfaces of quiescent, accreting neutron stars. We searched with Chandra for such a surface component of emission from the dynamical black-hole candidate XTE J1118+480 (= J1118), which has a primary mass of about 8 solar masses. None was found, as one would expect if the compact X-ray source is a bona fide black hole that possesses an event horizon. The X-ray spectrum of J1118 is well-fitted by a simple power-law model. In our search for a thermal component, we fitted our Chandra data to a power-law model (with slope and column density fixed) plus a series of nine hydrogen-atmosphere models with radii ranging from 9/8 to 2.8 Schwarzschild radii. For the more compact models, we included the important effect of self-irradiation of the atmosphere. Because of the extremely low column density to J1118, we obtained very strong limits on a hypothetical thermal source: kT(infinity) < 0.011 keV and L(infinity) < 9.4E30 erg/s (99% confidence level). In analogy with neutron stars, we examine two possible sources of thermal radiation from a hypothetical surface of J1118: deep crustal heating and accretion in quiescence. We conclude that if J1118 possesses a surface composed of normal nuclear matter, we would have detected thermal radiation from this surface. The absence of thermal radiation indicates that J1118 is a true black hole with an event horizon.

Recent X-ray surveys by {\it Chandra} and {\it XMM-Newton} have revealed a population of X-ray bright, optically normal galaxies (XBONGs) at moderate redshifts. We propose that many XBONGs are powered by an inner radiatively inefficient accretion flow (RIAF) plus an outer radiatively efficient thin accretion disk. The absence of optical/UV activity in XBONGs is explained by the truncation of the thin disk near the black hole, while the relatively strong X-ray emission is explained as inverse Compton emission from the hot RIAF. As an example, we show that the spectra of two XBONGs can be fit fairly well with such a model. By comparing these two sources to other accreting black holes, we argue that XBONGs are intermediate in their characteristics between distant luminous active galactic nuclei and nearby low-luminosity nuclei.

We discuss hydrostatic models of galaxy clusters in which heat diffusion balances radiative cooling. We consider two different sources of diffusion, thermal conduction and turbulent mixing, parameterized by dimensionless coefficients, f and alpha_mix, respectively. Models with thermal conduction give reasonably good fits to the density and temperature profiles of several cooling flow clusters, but some clusters require unphysically large values of f>1. Models with turbulent mixing give good fits to all clusters, with reasonable values of alpha_mix ~ 0.01-0.03. Both types of models are found to be essentially stable to thermal perturbations. The mixing model reproduces the observed scalings of various cluster properties with temperature, and also explains the entropy floor seen in galaxy groups.

The existence of black holes remains open to doubt until other conceivable options are excluded. With this motivation, we consider a model of a compact star in which most of the mass consists of dark particles of some kind, and a small fraction of the mass is in the form of ordinary nucleonic gas. The gas does not interact with the dark matter other than via gravity, but collects at the center as a separate fermionic fluid component. Depending on whether the dark mass is made of fermions or bosons, the objects may be called fermion-fermion stars or boson-fermion stars, respectively. For appropriate choices of the mass of the dark matter particles, these objects are viable models of black hole candidates in X-ray binaries. We consider models with a dark mass of 10 solar masses and a range of gas mass from 10^{-6} to nearly one solar mass, and analyse the bursting properties of the models when they accrete gas. We show that all the models would experience thermonuclear Type I X-ray bursts at appropriate mass accretion rates. Since no Type I bursts have been reported from black hole candidates, the models are ruled out. The case for identifying black hole candidates in X-ray binaries as true black holes is thus strengthened.

Recent infrared (IR) observations of the center of our Galaxy indicate that the supermassive black hole source Sgr A* is strongly variable in the IR. The timescale for the variability, $\sim 30$ min, is comparable to that of the X-ray flares observed by {\em Chandra} and {\em XMM}, suggesting a common physical origin. In this paper, we investigate the nature of the IR emission in the context of models recently proposed to interpret the X-ray flares. We show that the IR emission in Sgr A* can be well explained by nonthermal synchrotron emission if a small fraction of the electrons in the innermost region of the accretion flow around the black hole are accelerated {into a broken power-law distribution}, perhaps due to transient events such as magnetic reconnection. The model predicts differences in the variability amplitudes of flares in the IR and X-rays, in general agreement with observations. It also predicts that the IR emission should be linearly polarized, as has indeed been observed during one epoch. IR and X-ray flares analogous to those observed in Sgr A* may be detectable from other accreting SMBHs, provided $L \lsim 10^{-8} L_{EDD}$; at higher luminosities the flaring emission is dominated by thermal synchrotron-self Compton emission, which is likely to be less variable.

Astronomers have discovered many candidate black holes in X-ray binaries and in the nuclei of galaxies. The candidate objects are too massive to be neutron stars, and for this reason they are considered to be black holes. While the evidence based on mass is certainly strong, there is no proof yet that any of the objects possesses the defining characteristic of a black hole, namely an event horizon. Type I X-ray bursts, which are the result of thermonuclear explosions when gas accretes onto the surface of a compact star, may provide important evidence in this regard. Type I bursts are commonly observed in accreting neutron stars, which have surfaces, but have never been seen in accreting black hole candidates. It is argued that the lack of bursts in black hole candidates is compelling evidence that these objects do not have surfaces. The objects must therefore possess event horizons.

We consider the interface between an accretion disk and a magnetosphere surrounding the accreting mass. We argue that such an interface can occur not only with a magnetized neutron star but also sometimes with an unmagnetized neutron star or a black hole. The gas at the magnetospheric interface is generally Rayleigh-Taylor unstable and may also be Kelvin-Helmholtz unstable. Because of these instabilities, modes with low azimuthal wavenumbers m are expected to grow to large amplitude. It is proposed that the resulting nonaxisymmetric structures contribute to the high frequency quasi-periodic oscillations that have been seen in neutron-star and black-hole X-ray binaries. The mode oscillation frequencies are calculated to be approximately equal to m \Omega_m, where \Omega_m is the angular velocity of the accreting gas at the magnetospheric radius. Thus, mode frequencies should often be in the approximate ratio 1:2:3, etc. If the pressure of the gas in the disk is not large, then the m = 1 mode will be stable. In this case, the mode frequencies should be in the approximate ratio 2:3, etc. There is some observational evidence for such simple frequency ratios.

We consider an accretion flow model originally proposed by Bisnovatyi-Kogan & Ruzmaikin (1974), which has been confirmed in recent 3D MHD simulations. In the model, the accreting gas drags in a strong poloidal magnetic field to the center such that the accumulated field disrupts the axisymmetric accretion flow at a relatively large radius. Inside the disruption radius, the gas accretes as discrete blobs or streams with a velocity much less than the free-fall velocity. Almost the entire rest mass energy of the gas is released as heat, radiation and mechanical/magnetic energy. Even for a non-rotating black hole, the efficiency of converting mass to energy is of order 50% or higher. The model is thus a practical analog of an idealized engine proposed by Geroch and Bekenstein.

We present a spherically-symmetric, steady-state model of galaxy clusters in which radiative cooling from the hot gas is balanced by heat transport through turbulent mixing. We assume that the gas is in hydrostatic equilibrium, and describe the turbulent heat diffusion by means of a mixing length prescription with a dimensionless parameter alpha_mix. Models with alpha_mix ~ 0.01-0.03 yield reasonably good fits to the observed density and temperature profiles of cooling core clusters. Making the strong simplification that alpha_mix is time-independent and that it is roughly the same in all clusters, the model reproduces remarkably well the observed scalings of X-ray luminosity, gas mass fraction and entropy with temperature. The break in the scaling relations at kT \~ 1-2 keV is explained by the break in the cooling function at around this temperature, and the entropy floor observed in galaxy groups is reproduced naturally.

We present a global linear stability analysis of nuclear fuel accumulating on the surface of an accreting neutron star and we identify the conditions under which thermonuclear bursts are triggered. The analysis reproduces all the recognized regimes of hydrogen and helium bursts, and in addition shows that at high accretion rates, near the limit of stable burning, there is a regime of ``delayed mixed bursts'' which is distinct from the more usual ``prompt mixed bursts.'' In delayed mixed bursts, a large fraction of the fuel is burned stably before the burst is triggered. Bursts thus have longer recurrence times, but at the same time have somewhat smaller fluences. Therefore, the parameter alpha, which measures the ratio of the energy released via accretion to that generated through nuclear reactions in the burst, is up to an order of magnitude larger than for prompt bursts. This increase in alpha near the threshold of stable burning has been seen in observations. We explore a wide range of mass accretion rates, neutron star radii and core temperatures, and calculate a variety of burst properties. From a preliminary comparison with data, we suggest that bursting neutron stars may have hot cores, with T_{core} >~ 10^{7.5} K, consistent with interior cooling via the modified URCA or similar low-efficiency process, rather than T_{core} ~ 10^7 K, as expected for the direct URCA process. There is also an indication that neutron star radii are somewhat small <~ 10 km. Both of these conclusions need to be confirmed by comparing more careful calculations with better data.

We investigate radiatively inefficient accretion flow models for Sgr A*, the supermassive black hole in our Galactic Center, in light of new observational constraints. Confirmation of linear polarization in the submm emission argues for accretion rates much less than the canonical Bondi rate. We consider models with low accretion rates and calculate the spectra produced by a hybrid electron population, consisting of both thermal and nonthermal particles. The thermal electrons produce the submm emission and can account for its linear polarization properties. As noted in previous work, the observed low-frequency radio spectrum can be explained if a small fraction ($\approx 1.5$%) of the electron thermal energy resides in a soft power-law tail. In the innermost region of the accretion flow, turbulence and/or magnetic reconnection events may occasionally accelerate a fraction of the electrons into a harder power-law tail. We show that the synchrotron emission from these electrons, or the Compton up-scattering of synchrotron photons by the same electrons, may account for the X-ray flares observed by {\it Chandra}.

We consider a model of galaxy clusters in which the hot gas is in hydrostatic equilibrium and maintains energy balance between radiative cooling and heating by thermal conduction. We analyze the thermal stability of the gas using a Lagrangian perturbation analysis. For thermal conductivity at the level of 20-40 % of Spitzer conductivity, consistent with previous estimates for cluster gas, we find that the growth rate of the most unstable global radial mode is ~ 6-9 times lower than the growth rate of local isobaric modes at the cluster center in the absence of conduction. The growth time in typical clusters is ~ 2-5 Gyr, which is comparable to the time since the last major merger episode, when the gas was presumably well mixed. Thus, we suggest that thermal instability is not dynamically significant in clusters, provided there is an adequate level of thermal conduction. On the other hand, if the heating of the gas is not the result of thermal conduction or any other diffusive process such as turbulent mixing, then the thermal instability has a growth time under a Gyr in the central regions of the cluster and is a serious threat to equilibrium. We also analyze local nonradial modes and show that the Lagrangian technique leads to the same dispersion relation as the Eulerian approach, provided that clusters are initially in strict thermal equilibrium. Because cluster gas is convectively stable, nonradial modes always have a smaller growth rate than equivalent radial modes.

A large number of neutron stars (NSs), ~10^9, populate the Galaxy, but only a tiny fraction of them is observable during the short radio pulsar lifetime. The majority of these isolated NSs, too cold to be detectable by their own thermal emission, should be visible in X-rays as a result of accretion from the interstellar medium. The ROSAT all sky survey has however shown that such accreting isolated NSs are very elusive: only a few tentative candidates have been identified, contrary to theoretical predictions that up to several thousands should be seen. We suggest that the fundamental reason for this discrepancy lies in the use of the standard Bondi formula to estimate the accretion rates. We compute the expected source counts using updated estimates of the pulsar velocity distribution, realistic hydrogen atmosphere spectra, and a modified expression for the Bondi accretion rate as suggested by recent MHD simulations, and supported by direct observations in the case of accretion around supermassive black holes in nearby galaxies and in our own. We find that, whereas the inclusion of atmospheric spectra partly compensates for the reduction in the counts due to the higher mean velocities of the new distribution, the modified Bondi formula dramatically suppresses the source counts. The new predictions are consistent with a null detection at the ROSAT sensitivity.

Medvedev & Narayan have described a hot accretion flow onto a spinning neutron star in which the gas viscously brakes the spin of the star. Their self-similar solution has the surprising property that the density, temperature and angular velocity of the gas at any radius are completely independent of the outer boundary conditions. Hence, the solution cannot be matched to a general external medium. We resolve this paradoxical situation by showing that there is a second self-similar solution which bridges the gap between the original solution and the external medium. This new solution has an extra degree of freedom which permits it to match general outer boundary conditions. We confirm the main features of the analytical results with a full numerical solution.

We present an X-ray/UV/optical spectrum of the black hole primary in the X-ray nova XTE J1118+480 in quiescence at L_x = 4 x 10^{-9} of the Eddington luminosity. The Chandra, HST and MMT spectroscopic observations were performed simultaneously on 2002 January 12 UT. Because this 4.1-hr binary is located at b = 62 deg, the transmission of the ISM is very high (e.g., 70% at 0.3 keV). We present many new results for the quiescent state, such as the first far-UV spectrum and evidence for an 0.35 mag orbital modulation in the near-UV flux. However, the centerpiece of our work is the multiwavelength spectrum of XTE J1118+480, which we argue represents the canonical spectrum of a stellar-mass black hole radiating at L_x = 4 x 10^{-8.5} of the Eddington luminosity. This spectrum is comprised of two apparently disjoint components: a hard X-ray spectrum with a photon index Gamma = 2.02 +/- 0.16, and an optical/UV continuum that resembles a 13,000 K disk blackbody spectrum punctuated by several strong emission lines. We present a model of the source in which the accretion flow has two components: (1) an X-ray-emitting interior region where the flow is advection-dominated, and (2) a thin, exterior accretion disk with a truncated inner edge (R_tr ~ 10^4 Schwarzschild radii) that is responsible for the optical/UV spectrum. For D = 1.8 kpc, the luminosity of the X-ray component is L_x = 3.5 x 10^{30} erg/s (0.3-7 keV); the bolometric luminosity of the optical/UV component is 20 times greater.

We present three-dimensional MHD simulations of rotating radiatively inefficient accretion flows onto black holes. In the simulations, we continuously inject magnetized matter into the computational domain near the outer boundary, and we run the calculations long enough for the resulting accretion flow to reach a quasi-steady state. We have studied two limiting cases for the geometry of the injected magnetic field: pure toroidal field and pure poloidal field. In the case of toroidal field injection, the accreting matter forms a nearly axisymmetric, geometrically-thick, turbulent accretion disk. The disk resembles in many respects the convection-dominated accretion flows found in previous numerical and analytical investigations of viscous hydrodynamic flows. Models with poloidal field injection evolve through two distinct phases. In an initial transient phase, the flow forms a relatively flattened, quasi-Keplerian disk with a hot corona and a bipolar outflow. However, when the flow later achieves steady state, it changes in character completely. The magnetized accreting gas becomes two-phase, with most of the volume being dominated by a strong dipolar magnetic field from which a thermal low-density wind flows out. Accretion occurs mainly via narrow slowly-rotating radial streams which `diffuse' through the magnetic field with the help of magnetic reconnection events.

We present a simple model of hot gas in galaxy clusters, assuming hydrostatic equilibrium and energy balance between radiative cooling and thermal conduction. For five clusters, A1795, A1835, A2199, A2390 and RXJ1347.5-1145, the model gives a good description of the observed radial profiles of electron density and temperature, provided we take the thermal conductivity $\kappa$ to be about 30% of the Spitzer conductivity. Since the required $\kappa$ is consistent with the recent theoretical estimate of Narayan & Medvedev (2001) for a turbulent magnetized plasma, we consider a conduction-based equilibrium model to be viable for these clusters. We further show that the hot gas is thermally stable because of the presence of conduction. For five other clusters, A2052, A2597, Hydra A, Ser 159-03 and 3C295, the model requires unphysically large values of $\kappa$ to fit the data. These clusters must have some additional source of heat, most likely an active galactic nucleus since all the clusters have strong radio galaxies at their centers. We suggest that thermal conduction, though not dominant in these clusters, may nevertheless play a significant role by preventing the gas from becoming thermally unstable.

It has been suggested that quasi-periodic oscillations of accreting X-ray sources may relate to the modes named in the title. We consider non-axisymmetric linear perturbations to an isentropic, isothermal, unmagnetized thin accretion disk. The radial wave equation, in which the number of vertical nodes (n) appears as a separation constant, admits a wave-action current that is conserved except, in some cases, at corotation. Waves without vertical nodes amplify when reflected by a barrier near corotation. Their action is conserved. As was previously known, this amplification allows the n=0 modes to be unstable under appropriate boundary conditions. In contrast, we find that waves with n >0 are strongly absorbed at corotation rather than amplified; their action is not conserved. Therefore, non-axisymmetric p-modes and g-modes with n>0 are damped and stable even in an inviscid disk. This eliminates a promising explanation for quasi-periodic oscillations in neutron-star and black-hole X-ray binaries.

In a recent paper (astro-ph/0207561, posted on July 25, 2002), Fromang & Balbus state that ``a uniformly rotating barotropic fluid in an external potential attains a true energy minimum if and only if the rotation profile is everywhere subsonic. If regions of supersonic rotation are present, fluid variations exist that could take the system to states of lower energy.'' We present a simple counter-example which demonstrates that the stablity criterion presented by Fromang & Balbus is incorrect.

Many models of Gamma Ray Bursts invoke a central engine consisting of a black hole of a few solar masses accreting matter from a disk at a rate of a fraction to a few solar masses per second. Popham et al. and Narayan et al. have shown that, for Mdot >~ 0.1 Msun/s, accretion proceeds via neutrino cooling and neutrinos can carry away a significant amount of energy from the inner regions of the disks. We improve on these calculations by including a simple prescription for neutrino transfer and neutrino opacities in such regions. We find that the flows become optically thick to neutrinos inside a radius R~6-40R_s for Mdot in the range of 0.1 -10 Msun/s, where R_s is the black hole Schwarzchild radius. Most of the neutrino emission comes from outside this region and, the neutrino luminosity stays roughly constant at a value L_{\nu} ~ 10^{53} erg/s. We show that, for Mdot > 1 Msun/s, neutrinos are sufficiently trapped that energy advection becomes the dominant cooling mechanism in the flow. These results imply that neutrino annihilation in hyperaccreting black holes is an inefficient mechanism for liberating large amounts of energy. Extraction of rotational energy by magnetic processes remains the most viable mechanism.

Type I X-ray bursts are very common in neutron star X-ray binaries, but no Type I burst has been seen in the dozen or so binaries in which the accreting compact star is too massive to be a neutron star and therefore is identified as a black hole candidate. We have carried out a global linear stability analysis of the accumulating fuel on the surface of a compact star to identify the conditions under which thermonuclear bursts are triggered. Our analysis, which improves on previous calculations, reproduces the gross observational trends of bursts in neutron star systems. It further shows that, if black hole candidates have surfaces, they would very likely exhibit instabilities similar to those that lead to Type I bursts on neutron stars. The lack of bursts in black hole candidates is thus significant, and indicates that these objects have event horizons. We discuss possible caveats to this conclusion.

X-ray lines have been recently detected in the afterglows of a few gamma-ray bursts. We derive constraints on the physical conditions in the line-emitting gas, using as an example the multiple K$_\alpha$ lines detected by Reeves et al. (2002) in GRB 011211. We argue that models previously discussed in the literature require either a very extreme geometry or too much mass in the line-emitting region. We propose a new model in which gamma-rays and radiation from the early x-ray afterglow are back-scattered by an electron-positron pair screen at a distance of about $10^{14}-10^{15}$ cm from the source and irradiate the expanding outer layers of the supernova ejecta, thereby producing x-ray lines. The model suffers from fewer problems compared to previous models. It also has the advantage of requiring only a single explosion to produce both the GRB and the supernova ejecta, in contrast to most other models for the lines which require the supernova to go off days or weeks prior to the GRB. The model, however, has difficulty explaining the $>10^{48}$ ergs of energy emitted in the x-ray lines, which requires somewhat extreme choices of model parameters. The difficulties associated with the various models are not particular to GRB 011211. They are likely to pose a problem for any GRB with x-ray lines.

Radiatively inefficient accretion flows onto black holes are unstable due to both an outwardly decreasing entropy (`convection') and an outwardly decreasing rotation rate (the `magnetorotational instability'; MRI). Using a linear magnetohydrodynamic stability analysis, we show that long-wavelength modes are primarily destabilized by the entropy gradient and that such `convective' modes transport angular momentum inwards. Moreover, the stability criteria for the convective modes are the standard Hoiland criteria of hydrodynamics. By contrast, shorter wavelength modes are primarily destabilized by magnetic tension and differential rotation. These `MRI' modes transport angular momentum outwards. The convection-dominated accretion flow (CDAF) model, which has been proposed for radiatively inefficient accretion onto a black hole, posits that inward angular momentum transport and outward energy transport by long-wavelength convective fluctuations are crucial for determining the structure of the accretion flow. Our analysis suggests that the CDAF model is applicable to a magnetohydrodynamic accretion flow provided the magnetic field saturates at a sufficiently sub-equipartition value (plasma beta >> 1), so that long-wavelength convective fluctuations can fit inside the accretion disk. Numerical magnetohydrodynamic simulations are required to determine whether such a sub-equipartition field is in fact obtained.

We model the 1990 giant X-ray flare of the quiescent galaxy NGC 5905 as the tidal disruption of a star by a supermassive black hole. From the observed rapid decline of the luminosity, over a timescale of a few years, we argue that the flare was powered by the fallback of debris rather than subsequent accretion via a thin disk. The fallback model allows constraints to be set on the black hole mass and the mass of debris. The latter must be very much less than a solar mass to explain the very low luminosity of the flare. The observations can be explained either as the partial stripping of the outer layers of a low-mass main sequence star or as the disruption of a brown dwarf or a giant planet. We find that the X-ray emission in the flare must have originated within a small patch rather than over the entire torus of circularized material surrounding the black hole. We suggest that the patch corresponds to the ``bright spot'' where the stream of returning debris impacts the torus. Interestingly, although the peak luminosity of the flare was highly sub-Eddington, the peak flux from the bright spot was close to the Eddington limit. We speculate on the implications of this result for observations of other flare events.

Roughly a dozen X-ray binaries are presently known in which the compact accreting primary stars are too massive to be neutron stars. These primaries are identified as black holes, though there is as yet no definite proof that any of the candidate black holes actually possesses an event horizon. We discuss how Type I X-ray bursts may be used to verify the presence of the event horizon in these objects. Type I bursts are caused by thermonuclear explosions when gas accretes onto a compact star. The bursts are commonly seen in many neutron star X-ray binaries, but they have never been seen in any black hole X-ray binary. Our model calculations indicate that black hole candidates ought to burst frequently if they have surfaces. Based on this, we argue that the lack of bursts constitutes strong evidence for the presence of event horizons in these objects.

Nearby galactic nuclei are observed to be very much dimmer than active galactic nuclei in distant galaxies. The Chandra X-ray Observatory has provided a definitive explanation for why this is so. With its excellent angular resolution, Chandra has imaged hot X-ray-emitting gas close to the gravitational capture radius of a handful of supermassive black holes, including Sgr A* in the nucleus of our own Galaxy. These observations provide direct and reliable estimates of the Bondi mass accretion rate in these nuclei. It is found that the Bondi accretion rate is significantly below the Eddington mass accretion rate, but this alone does not explain the dimness of the accretion flows. In all the systems observed so far, the accretion luminosity is orders of magnitude less than Mdot_{Bondi} times c^2, which means that the accretion must occur via a radiatively inefficient mode. This conclusion, which was strongly suspected for many years, is now inescapable. Furthermore, if the accretion in these nuclei occurs via either a Bondi flow or an advection-dominated accretion flow, the accreting plasma must be two-temperature at small radii, and the central mass must have an event horizon. Convection, winds and jets may play a role, but observations do not yet permit definite conclusions.

We estimate the thermal conductivity of a weakly collisional magnetized plasma with chaotic magnetic field fluctuations. When the fluctuation spectrum extends over two or more decades in wave-vector, we find that thermal conduction is very efficient; the conduction coefficient is only a factor ~5 below the classical Spitzer estimate. We suggest that conduction could play a significant role in cooling flows in clusters of galaxies.

We consider the radial structure of radiatively inefficient hydrodynamic accretion flows around black holes. We show that low-viscosity flows consist of two zones: an outer convection-dominated zone and an inner advection-dominated zone. The transition between these two zones occurs at ~50 Schwarzschild radii.

We discuss new observations of X-ray novae which provide strong evidence that black holes have event horizons. Optical observations of 13 X-ray novae indicate that these binary stars contain collapsed objects too heavy to be stable neutron stars. The objects have been identified as black hole candidates. X-ray observations of several of these X-ray novae in quiescence with the Chandra X-ray Observatory show that the systems are approximately 100 times fainter than nearly identical X-ray novae containing neutron stars. The advection-dominated accretion flow model provides a natural explanation for the difference. In this model, the accreting gas reaches the accretor at the center with a large amount of thermal energy. If the accretor is a black hole, the thermal energy will disappear through the event horizon, and the object will be very dim. If the accretor is a neutron star or any other object with a surface, the energy will be radiated from the surface, and the object will be bright. We discuss alternate interpretations of the data that eliminate the need for advection-dominated accretion. Most of these alternatives still require an event horizon to explain the unusually low X-ray luminosities of the black hole candidates. Some of the alternatives are also inconsistent with observations.

If the efficiency for producing $\gamma$-rays is the same in short duration ($\siml 2$ s) Gamma-Ray Bursts (GRBs) as in long duration GRBs, then the average kinetic energy of short GRBs must be $\sim 20$ times less than that of long GRBs. Assuming further that the relativistic shocks in short and long duration GRBs have similar parameters, we show that the afterglows of short GRBs will be on average 10--40 times dimmer than those of long GRBs. We find that the afterglow of a typical short GRB will be below the detection limit ($\siml 10 \microJy$) of searches at radio frequencies. The afterglow would be difficult to observe also in the optical, where we predict $R \simg 23$ a few hours after the burst. The radio and optical afterglow would be even fainter if short GRBs occur in a low-density medium, as expected in NS-NS and NS-BH merger models. The best prospects for detecting short-GRB afterglows are with early ($\siml 1$ day) observations in X-rays.

Recently, Medvedev & Narayan (2001) discovered a new type of accretion flow, a hot settling flow around a rapidly rotating neutron star. The flow is cooling-dominated and energetically similar to the Shapiro, Lightman, & Eardley (1976, SLE) solution. Since the SLE solution is known to be thermally unstable, one might suspect that the new solution would also be unstable. However, due to the very high temperature of the accreting gas, thermal conduction is very strong and could suppress the thermal instability. We analyze the role of thermal conduction in both the hot settling flow and the SLE solution. In the hot settling flow collisions are very rare. Therefore, thermal transport occurs via free streaming of electrons along tangled magnetic field lines. We find that conduction is strong enough to make the flow marginally stable. In contrast, in the cooler SLE solution, conduction is via collisional, Spitzer-type transport. In this case, conduction is weaker, and we find that the SLE solution is thermally unstable even in the presence of conduction.

We present three-dimensional numerical magnetohydrodynamic simulations of radiatively inefficient spherical accretion onto a black hole. The simulations are initialized with a Bondi flow, and with a weak, dynamically unimportant, large-scale magnetic field. The magnetic field is amplified as the gas flows in. When the magnetic pressure approaches equipartition with the gas pressure, the field begins to reconnect and the gas is heated up. The heated gas is buoyant and moves outward, causing line stretching of the frozen-in magnetic field. This leads to further reconnection, and more heating and buoyancy-induced motions, so that the flow makes a transition to a state of self-sustained convection. The radial structure of the flow changes dramatically from its initial Bondi profile, and the mass accretion rate onto the black hole decreases significantly. Motivated by the numerical results, we develop a simplified analytical model of a radiatively inefficient spherical flow in which convective transport of energy to large radii plays an important role. In this "convection-dominated Bondi flow" the accretion velocity is highly subsonic and the density varies with radius as ~R^{-1/2} rather than the standard Bondi scaling ~R^{-3/2}. We estimate that the mass accretion rate onto the black hole is significantly less than the Bondi accretion rate. Convection-dominated Bondi flows may be relevant for understanding many astrophysical phenomena, e.g. post-supernova fallback and radiatively inefficient accretion onto supermassive black holes, stellar-mass black holes and neutron stars.

We present analytical and numerical solutions which describe a hot, viscous, two-temperature accretion flow onto a rotating neutron star or any other rotating compact star with a surface. We assume Coulomb coupling between the protons and electrons, and free-free cooling from the electrons. Outside a thin boundary layer, where the accretion flow meets the star, we show that there is an extended settling region which is well-described by two self-similar solutions: (i) a two-temperature solution which is valid in an inner zone $r\le10^{2.5}$ ($r$ is in Schwarzchild units), and (ii) a one-temperature solution at larger radii. In both zones, $\rho\propto r^{-2}, \Omega\propto r^{-3/2}, v\propto r^0, T_p\propto r^{-1}$; in the two-temperature zone, $T_e\propto r^{-1/2}$. The luminosity of the settling zone arises from the rotational energy of the star as the star is braked by viscosity. Hence the luminosity and the flow parameters (density, temperature, angular velocity) are independent of $\dot M$. The settling solution described here is not advection-dominated, and is thus different from the self-similar ADAF found around black holes. When the spin of the star is small enough, however, the present solution transforms smoothly to a (settling) ADAF. We carried out a stability analysis of the settling flow. The flow is convectively and viscously stable and is unlikely to have strong winds or outflows. Unlike another cooling-dominated system --- the SLE disk, --- the settling flow is thermally stable provided that thermal conduction is taken into account. This strong saturated-like thermoconduction does not change the structure of the flow.

Many models of gamma-ray bursts (GRBs) involve accretion onto a compact object, usually a black hole, at a mass accretion rate of order a fraction of a solar mass per second. If the accretion disk is larger than a few tens or hundreds of Schwarzschild radii, the accretion will proceed via a convection-dominated accretion flow (CDAF) in which most of the matter escapes to infinity rather than falling onto the black hole. Models involving the mergers of black hole white dwarf binaries and black hole helium star binaries fall in this category. These models are unlikely to produce GRBs since very little mass reaches the black hole. If the accretion disk is smaller, then accretion will proceed via neutrino cooling in a neutrino-dominated accretion disk (NDAF) and most of the mass will reach the center. Models involving the mergers of double neutron star binaries and black hole neutron star binaries fall in this category and are capable of producing bright GRBs. If the viscosity parameter $\alpha$ in the NDAF has a standard value $\sim0.1$, these mergers can explain short GRBs with durations under a second, but they are unlikely to produce long GRBs with durations of tens or hundred of seconds. If the accretion disk is fed by fallback of material after a supernova explosion, as in the collapsar model, then the time scale of the burst is determined by fallback, not accretion. Such a model can produce long GRBs. Fallback models again require that the accretion should proceed via an NDAF rather than a CDAF in order for a significant amount of mass to reach the black hole. This condition imposes an upper limit on the radius of injection of the gas.

We estimate the coalescence rate of close binaries with two neutron stars (NS) and discuss the prospects for the detection of NS-NS inspiral events by ground-based gravitational-wave observatories, such as LIGO. We derive the Galactic coalescence rate using the observed sample of close NS-NS binaries (PSR B1913+16 and PSR B1534+12) and examine in detail each of the sources of uncertainty associated with the estimate. Specifically, we investigate (i) the dynamical evolution of NS-NS binaries in the Galactic potential and the vertical scale height of the population, (ii) the pulsar lifetimes, (iii) the effects of the faint end of the radio pulsar luminosity function and their dependence on the small number of observed objects, (iv) the beaming fraction, and (v) the extrapolation of the Galactic rate to extragalactic distances expected to be reachable by LIGO. We find that the dominant source of uncertainty is the correction factor (up to about 200) for faint (undetectable) pulsars. All other sources are much less important, each with uncertainty factors smaller than 2. Despite the relatively large uncertainty, the derived coalescence rate is approximately consistent with previously derived upper limits, and is more accurate than rates obtained from population studies. We obtain a most conservative lower limit for the LIGO II detection rate of 2 events per year. Our upper limit on the detection rate lies between 300 to more than 1000 events per year.

For small values of the dimensionless viscosity parameter, namely $\alpha\lesssim 0.1$, the dynamics of non-radiating accretion flows is dominated by convection; convection strongly suppresses the accretion of matter onto the central object and transports a luminosity $\sim 10^{-3}-10^{-2} \dot M c^2$ from small to large radii in the flow. A fraction of this convective luminosity is likely to be radiated at large radii via thermal bremsstrahlung emission. We show that this leads to a correlation between the frequency of maximal bremsstrahlung emission and the luminosity of the source, $\nu_{\rm peak} \propto L^{2/3}$. Accreting black holes with X-ray luminosities $10^{-4} L_{Edd}\gtrsim L_X(0.5-10{\rm keV}) \gtrsim 10^{-7}L_{Edd}$ are expected to have hard X-ray spectra, with photon indices $\Gamma\sim2$, and sources with $L_X\lesssim 10^{-9}L_{Edd}$ are expected to have soft spectra, with $\Gamma\sim3.5$. This is testable with {\it Chandra} and {\it XMM}.

Previously we claimed that Black Hole X-ray Novae (BHXN) in quiescence are much less luminous than equivalent Neutron Star X-ray Novae (NSXN). This claim was based on the quiescent detection of a single short period BHXN (A0620-00, P(orb)=7.8 hrs) and two longer period BHXN (GRO J1655-40, P(orb)=62.9 hrs; V404 Cyg, P(orb)=155.3 hrs), along with sensitive upper limits. We announce the detection of two more short period BHXN (GRO J0422+32, P(orb)=5.1 hrs; GS 2000+25, P(orb)=8.3 hrs), an upper limit for a third which is improved by two orders of magnitude (4U 1543-47, P(orb)=27.0 hrs) and a new, much lower quiescent measurement of GRO J1655-40. Taken together, these new Chandra measurements confirm that the quiescent X-ray luminosities of BHXN are significantly lower than those of NSXN. We argue that this provides strong evidence for the existence of event horizons in BHXN.

We consider hot, two-temperature, viscous accretion onto a rotating, unmagnetized neutron star. We assume Coulomb coupling betweenthe protons and electrons, and free-free cooling from the electrons. We show that the accretion flow has an extended settling region which can be described by means of two analytical self-similar solutions: a two-temperature solution which is valid in an inner zone, r<10^{2.5}, where r is the radius in Schwarzchild units; and a one-temperature solution which is valid in an outer zone, r>10^{2.5}. In both zones the density varies as \rho ~ r^{-2} and the angular velocity as \Omega ~ r^{-3/2}. We solve the flow equations numerically and confirm that the analytical solutions are accurate. The settling flow radiates the energy dissipated by viscosity; so it is not advection-dominated. Except for the radial velocity, all other gas properties - density, angular velocity, temperature, luminosity, angular momentum flux - are independent of the mass accretion rate; these quantities do depend sensitively on the spin of the neutron star. The angular momentum flux is outward under most conditions; therefore, the central star is nearly always spun-down. The luminosity of the settling zone arises from the rotational energy that is released as the star is braked by viscosity, and the contribution from gravity is small; hence the radiative efficiency can be arbitrarily large at low $\mdot$. For reasonable values of the gas adiabatic index, the Bernoulli parameter is negative; therefore, a strong outflow or wind is not expected. The flow is convectively stable, but may be thermally unstable. The described solution is not advection-dominated; however, when the spin of the star is small enough, it transforms smoothly to an advection-dominated branch of solution.

In some circumstances, the mass accretion rate M^dot onto a compact star may depend not only on external boundary conditions, but also on the radius R of the star. Writing the dependence as a power-law in which M^dot is proportional to R^p, we estimate the power-law index, p, for transient binary systems in a quiescent state. We use the observed luminosities $L$ of accreting neutron stars (R=10^6 cm) in soft X-ray transients and white dwarfs (R=10^9 cm) in similar cataclysmic variables, and estimate M^dot in each system through the relation L=(GM/R)*M^dot, where M is the mass of the star. From the available data we infer that p~0.9+-0.5. This radial dependence is consistent with radiatively inefficient accretion flows that either are convective or lose mass via a wind.

The spectacular first light observation by the Chandra X-Ray Observatory revealed an X-ray point source near the center of the 300 yr old Cas A supernova remnant. We present an analysis of the public X-ray spectral and timing data. No coherent pulsations were detected in the Chandra/HRC data. The 3-sigma upper limit on the pulsed fraction is <35% for P>20 ms. The Chandra/ACIS spectrum of the point source may be fit with an ideal blackbody (kT=0.5 keV), or with BB models modified by the presence of a NS atmosphere (kT=0.25-0.35 keV), but the temperature is higher and the inferred emitting area lower than expected for a 300 yr old NS according to standard cooling models. The spectrum may also be fit with a power law model (photon index 2.8-3.6). Both the spectral properties and the timing limits of the point source are inconsistent with a young Crab-like pulsar, but are quite similar to the properties of the anomalous X-ray pulsars. The spectral parameters are also very similar to those of the other radio-quiet X-ray point sources in the supernova remnants Pup A, RCW 103, and PKS 1209-52. Current limits on an optical counterpart for the Cas A point source rule out models that invoke fallback accretion onto a compact object if fallback disk properties are similar to those in quiescent low-mass X-ray binaries. However, the optical limits are marginally consistent with plausible alternative assumptions for a fallback disk. In this case, accreting NS models can explain the X-ray data, but an accreting BH model is not promising.

We compute X-ray emission lines from thermal plasma in hot accretion flows. We show that line profiles are strong probes of the gas dynamics, and we present line-ratio diagnostics which are sensitive to the distribution of mass with temperature in the flow. We show how these can be used to constrain the run of density with radius, and the size of the hot region. We also present diagnostics which are primarily sensitive to the importance of recombination versus collisional ionization, and which could help discriminate ADAFs from photoionization-dominated accretion disk coronae. We apply our results to the Galactic center source Sagittarius A* and to the nucleus of M87. We find that the brightest predicted lines are within the detection capability of current $X$-ray instruments.

Most models of anomalous X-ray pulsars (AXPs) account for the observed X-ray spectra and pulsations by means of radiation processes that occur on the surfaces of neutron stars. For any such model, general relativistic deflection of light severely suppresses the amplitude of the observed pulsations. We calculate the expected pulsation amplitudes of AXPs according to various models and compare the results with observations. We show that the high (<= 70%) pulse amplitudes observed in some AXPs can be accounted for only if the surface emission is localized (spot radius <40 degrees) and strongly beamed (cos^n[theta'] with n>2, where theta' is the angle to the normal). These constraints are incompatible with those cooling and magnetar models in which the observed X-rays originate as thermal emission from the neutron-star surface. Accretion models, on the other hand, are compatible with observations for a wide range of parameters. Finally, definitive conclusions cannot be reached on magnetospheric models, since their localization and beaming properties are not well understood.

We investigate the effect of a hybrid electron population, consisting of both thermal and non-thermal particles, on the synchrotron spectrum, image size, and image shape of a hot accretion flow onto a supermassive black hole. We find two universal features in the emitted synchrotron spectrum: (i) a prominent shoulder at low (< 10^11 Hz) frequencies that is weakly dependent on the shape of the electron energy distribution, and (ii) an extended tail of emission at high (> 10^13 Hz) frequencies whose spectral slope depends on the slope of the power-law energy distribution of the electrons. In the low-frequency shoulder, the luminosity can be up to two orders of magnitude greater than with a purely thermal plasma even if only a small fraction (< 1%) of the steady-state electron energy is in the non-thermal electrons. We apply the hybrid model to the Galactic center source, Sgr A*. The observed radio and IR spectra imply that at most 1% of the steady-state electron energy is present in a power-law tail in this source. This corresponds to no more than 10% of the electron energy injected into the non-thermal electrons and hence 90% into the thermal electrons. We show that such a hybrid distribution can be sustained in the flow because thermalization via Coulomb collisions and synchrotron self-absorption are both inefficient. The presence of non-thermal electrons enlarges the size of the radio image at low frequencies and alters the frequency dependence of the brightness temperature. A purely thermal electron distributions produces a sharp-edged image while a hybrid distribution causes strong limb brightening. These effects can be seen up to frequencies ~10^11 Hz and are accessible to radio interferometers.

Advection-dominated accretion flows (ADAFs) are known to be convectively unstable for low values of the viscosity parameter alpha. Two-dimensional axisymmetric hydrodynamic simulations of such flows reveal a radial density profile which is significantly flatter than the rho~r^{-3/2} expected for a canonical ADAF. The modified density profile is the result of inward transport of angular momentum by axisymmetric convective eddies. We present three-dimensional hydrodynamic simulations of convective ADAFs which are free of the assumption of axisymmetry. We find that the results are qualitatively and quantitatively similar to those obtained with two-dimensional simulations. In particular, the convective eddies are nearly axisymmetric and transport angular momentum inward.

We consider height-integrated equations of an advection-dominated accretion flow (ADAF), assuming that there is no mass outflow. We include convection through a mixing length formalism. We seek self-similar solutions in which the rotational velocity and sound speed scale as R^{-1/2}, where R is the radius, and consider two limiting prescriptions for the transport of angular momentum by convection. In one limit, the transport occurs down the angular velocity gradient, so convection moves angular momentum outward. In the other, the transport is down the specific angular momentum gradient, so convection moves angular momentum inward. We also consider general prescriptions which lie in between the two limits. When convection moves angular momentum outward, we recover the usual self-similar solution for ADAFs in which the mass density scales as rho ~ R^{-3/2}. When convection moves angular momentum inward, the result depends on the viscosity coefficient alpha. If alpha>alpha_{crit1} ~ 0.05, we once again find the standard ADAF solution. For alpha<alpha_{crit}, however, we find a non-accreting solution in which rho ~ R^{-1/2}. We refer to this as a "convective envelope" solution or a "convection-dominated accretion flow". Two-dimensional numerical simulations of ADAFs with values of alpha<0.03 have been reported by several authors. The simulated ADAFs exhibit convection. By virtue of their axisymmetry, convection in these simulations moves angular momentum inward, as we confirm by computing the Reynolds stress. The simulations give rho ~ R^{-1/2}, in good agreement with the convective envelope solution. The R^{-1/2} density profile is not a consequence of mass outflow.

We construct time-dependent models of accretion discs around black holes and neutron stars. We investigate the effect of evaporating the disc inner regions during quiescence on the predictions of the Disc Instability Model (DIM) for these systems. We do not include irradiation of the disc in the models. Removing the inner, most unstable parts of the accretion disc increases the predicted recurrence times. However, DIMs with values of the viscosity parameter alpha_hot ~ 0.1 and alpha_cold ~ 0.02 (typical of applications of the DIM to standard dwarf nova outbursts) fail to reproduce the long recurrence times of SXTs (unless we resort to fine-tuning of the parameters) independent of the evaporation strength. We show that models with evaporation and a smaller value of alpha_cold (~ 0.005) do reproduce the long recurrence times and the accretion rates at the level of the Eddington rate observed in outburst. The large difference between the values of alpha_hot and alpha_cold, if confirmed when disc irradiation is included, suggests that several viscosity mechanisms operate in these accretion discs. For some choices of parameters our models predict reflares during the decline from outburst. They are a physical property of the model and result from a heating front forming in the wake of an initial cooling front and subsequent, multiple front reflections. The reflares disappear in low-alpha models where front reflection can not occur.

The nature of the energy source powering anomalous X-ray pulsars is uncertain. Proposed scenarios involve either an ultramagnetized neutron star, or accretion onto a neutron star. We consider the accretion model proposed recently by Chatterjee, Hernquist & Narayan, in which a disk is fed by fallback material following a supernova. We compute the optical, infrared, and submillimeter emission expected from such a disk, including both viscous dissipation and reradiation of X-ray flux impinging on the disk from the pulsar. We find that it is possible with current instruments to put serious constraints on this and on other accretion models of AXPs. Fallback disks could also be found around isolated radio pulsars and we compute the corresponding spectra. We show that the excess emission in the R and I bands observed for the pulsar PSR 0656+14 is broadly consistent with emission from a disk.

We present a model for the anomalous X-ray pulsars (AXPs) in which the emission is powered by accretion from a fossil disk, established from matter falling back onto the neutron star following its birth. The time-dependent accretion drives the neutron star towards a ``tracking'' solution in which the rotation period of the star increases slowly, in tandem with the declining accretion rate. For appropriate choices of disk mass, neutron star magnetic field strength and initial spin period, we demonstrate that a rapidly rotating neutron star can be spun down to periods characteristic of AXPs on timescales comparable to the estimated ages of these sources. In other cases, accretion onto the neutron star switches off after a short time, and the star becomes an ordinary radio pulsar. Thus, in our picture, radio pulsars and AXPs are drawn from the same underlying population, in contrast to models involving neutron stars with ultrastrong magnetic fields, which require a new population of stars with very different properties.

We test whether the model proposed by Esin et al. to explain X-ray observations of the black hole soft X-ray transient GRS 1124-68 in outburst can also explain the optical lightcurves of a similar object A0620-00. We show that to reproduce the observed X-ray to optical flux ratio in A0620-00, we need to assume X-ray irradiation of the outer disk that is significantly in excess of what is expected for a standard planar disk. With enhanced irradiation, the Esin et al. model can reproduce the optical evolution of A0620-00 in outburst very well. Though we find that optical observations of GRS 1124-68 also imply enhanced X-ray heating of the outer disk, the irradiation appears to be a factor of roughly 3 weaker in this system than in A0620-00. This is surprising, since GRS 1124-68 has a larger disk. We speculate that enhanced irradiation may be due to disk warping, and that the degree of warping differs between the two binaries.

All axisymmetric self-similar equilibria of self-gravitating, rotating, isothermal systems are identified by solving the nonlinear Poisson equation analytically. There are two families of equilibria: (1) Cylindrically symmetric solutions in which the density varies with cylindrical radius as R^(-alpha), with 0 <= alpha <= 2. (2) Axially symmetric solutions in which the density varies as f(theta)/r^2, where `r' is the spherical radius and `theta' is the co-latitude. The singular isothermal sphere is a special case of the latter class with f(theta)=constant. The axially symmetric equilibrium configurations form a two-parameter family of solutions and include equilibria which are surprisingly asymmetric with respect to the equatorial plane. The asymmetric equilibria are, however, not force-free at the singular points r=0, infinity, and their relevance to real systems is unclear. For each hydrodynamic equilibrium, we determine the phase-space distribution of the collisionless analog.

M81 and NGC 4579 are two of the few low-luminosity active galactic nuclei which have an estimated mass for the central black hole, detected hard X-ray emission, and detected optical/UV emission. In contrast to the canonical ``big blue bump,'' both have optical/UV spectra which decrease with increasing frequency in a $\nu L_\nu$ plot. Barring significant reddening by dust and/or large errors in the black hole mass estimates, the optical/UV spectra of these systems require that the inner edge of a geometrically thin, optically thick, accretion disk lies at roughly 100 Schwarzschild radii. The observed X-ray radiation can be explained by an optically thin, two temperature, advection-dominated accretion flow at smaller radii.

Cooling flows in galaxy clusters and isolated elliptical galaxies are a source of mass for fueling accretion onto a central supermassive black hole. We calculate the dynamics of accreting matter in the combined gravitational potential of a host galaxy and a central black hole assuming a steady state, spherically symmetric flow (i.e., no angular momentum). The global dynamics depends primarily on the accretion rate. For large accretion rates, no simple, smooth transition between a cooling flow and an accretion flow is possible; the gas cools towards zero temperature just inside its sonic radius, which lies well outside the region where the gravitational influence of the central black hole is important. For accretion rates below a critical value, however, the accreting gas evolves smoothly from a radiatively driven cooling flow at large radii to a nearly adiabatic (Bondi) flow at small radii. We argue that this is the relevant parameter regime for most observed cooling flows. The transition from the cooling flow to the accretion flow should be observable in M87 with the {\it Chandra X-ray Observatory}.

The proposed Laser Interferometer Space Antenna is expected to detect gravitational waves from neutron stars and stellar-mass black holes spiraling into supermassive black holes in distant galactic nuclei. Analysis of the inspiral events will require careful comparison of the observed signals with theoretical waveform templates. The comparison could be seriously compromised if non-gravitational torques modify the orbit of the star. This paper estimates the torque exerted on an orbiting star as a result of hydrodynamic interactions with an accretion flow around the supermassive black hole. It is argued that the majority of inspiral events will take place in low luminosity galactic nuclei in which the mass accretion rate is low and the accretion occurs via an advection-dominated flow. The hydrodynamic torque is negligibly small in such systems and will have no effect on gravitational wave experiments.

The discovery of hard, power-law X-ray emission from a sample of six nearby elliptical galaxies, including the dominant galaxies of the Virgo,Fornax and Centaurus clusters (M87, NGC 1399 and NGC 4696, respectively), and NGC 4472, 4636 and 4649 in the Virgo cluster, has important implications for the study of quiescent supermassive black holes. We describe how the broad band spectral energy distributions of these galaxies, which accrete from their hot gaseous halos at rates comparable to their Bondi rates, can be explained by low-radiative efficiency accretion flows in which a significant fraction of the mass, angular momentum and energy is removed from the flows by winds. The observed suppression of the synchrotron component in the radio band and the systematically hard X-ray spectra, which are interpreted as thermal bremsstrahlung emission, support the conjecture that significant mass outflow is a natural consequence of systems accreting at low-radiative efficiencies. We briefly discuss an alternative model for the X-ray emission, namely that it is due to nonthermal synchrotron-self Compton processes in the accretion flow, or wind. This appears to require implausibly weak magnetic fields. Emission from a collimated jet viewed off axis should be distinguishable from the bremsstrahlung model by variability and thermal line emission studies. We argue that the difference in radiative efficiency between the nuclei of spiral and elliptical galaxies arises from the different manner in which interstellar gas is fed into the nuclei. In ellipticals, matter fed from the hot (slowly cooling) ISM is likely to be highly magnetized and with low specific angular momentum, both of which favor low-radiative efficiency accretion solutions and possibly the formation of the observed jets.

The radio source Sgr A* at the center of our Galaxy is believed to be a 2.6 x 10^6 solar mass black hole which accretes gas from the winds of nearby stars. We show that limits on the X-ray and infrared emission from the Galactic Center provide an upper limit of ~ 8 x 10^{-5} solar masses per year on the mass accretion rate in Sgr A*. The advection-dominated accretion flow (ADAF) model favors a rate < 10^{-5} solar masses per year. In comparison, the Bondi accretion rate onto Sgr A*, estimated using the observed spatial distribution of mass losing stars and assuming non-interacting stellar winds, is ~ 3 x 10^{-5} solar masses per year. There is thus rough agreement between the Bondi, the ADAF, and the X-ray inferred accretion rates for Sgr A*. We discuss uncertainties in these estimates, emphasizing the importance of upcoming observations by the Chandra X-ray observatory (CXO) for tightening the X-ray derived limits.

We consider the X-ray luminosity difference between neutron star and black hole soft X-ray transients (NS and BH SXTs) in quiescence. The current observational data suggest that BH SXTs are significantly fainter than NS SXTs. The luminosities of quiescent BH SXTs are consistent with the predictions of binary evolution models for the mass transfer rate if (1) accretion occurs via an ADAF in these systems and (2) the accreting compact objects have event horizons. The luminosities of quiescent NS SXTs are not consistent with the predictions of ADAF models when combined with binary evolution models, unless most of the mass accreted in the ADAF is prevented from reaching the neutron star surface. We consider the possibility that mass accretion is reduced in quiescent NS SXTs because of an efficient propeller and develop a model of the propeller effect that accounts for the observed luminosities. We argue that modest winds from ADAFs are consistent with the observations while strong winds are probably not.

Several models for the hard X-ray Background (XRB) suggest that it is due to the emission from heavily obscured AGN. Recent studies have revealed the presence of a new population of hard X-ray sources which must contribute significantly to the XRB. However, whether the majority of these sources are obscured AGN or some other type of object at present remains unclear. Here, we further examine the possibility that a significant fraction of the XRB comes from a population of galaxies undergoing advection-dominated accretion in the high-\mdot regime and thus produce intrinsically hard spectra. When the accretion rate is close to its critical value, above which an advection-dominated accretion flow (ADAF) no longer exists, the major contribution to X-ray emission is due to inverse Compton scattering of the soft seed photons produced by cyclo-synchrotron radiation. In this regime, the resulting ADAF spectra are relatively hard with a fairly constant X-ray spectral index \alpha \sim 0.2-0.4 and a spectral cut-off at \sim 200 keV. We show that the integrated emission from such sources can provide a good fit to the hard (> 2 keV) X-ray background, provided that the spectrum is dominated by the contribution from objects located at redshifts z \sim 2-3. The model requires most of the contribution to the XRB to be due to objects accreting at the critical rate.

We present model X-ray spectra of accreting black holes with advection-dominated accretion flows, paying attention to thermal emission lines from the hot plasma. We show that the Advanced X-ray Astrophysical Facility (AXAF) might be able to observe lines from X-ray binaries such as V404 Cyg in quiescence, the Galactic Center black hole Sagittarius A*, and the nuclei of nearby galaxies such as M87. Line intensities can provide new diagnostics to study the accreting plasma in these and related systems.

We study the thermal and viscous stability of accretion flows in Low Mass Black Hole Binaries (LMBHBs). We consider a model in which an inner advection-dominated accretion flow (ADAF) is surrounded by a geometrically thin accretion disk, the transition between the two zones occurring at a radius R_tr. In all the known LMBHBs, R_tr appears to be such that the outer disks could suffer from a global thermal-viscous instability. This instability is likely to cause the transient behavior of these systems. However, in most cases, if R_tr were slightly larger than the estimated values, the systems would be globally stable. This suggests that a population of faint persistent LMBHBs with globally stable outer disks could be present in the Galaxy. Such LMBHBs would be hard to detect because they would lack large amplitude outbursts, and because their ADAF zones would have very low radiative efficiencies, making the systems very dim. We present model spectra of such systems covering the optical and X-ray bands.

We calculate spectral models of advection-dominated accretion flows, taking into account the possibility that significant mass may be lost to a wind. We apply the models to the soft X-ray transient V404 Cyg in quiescence and the Galactic center source Sgr A*. We show that there are qualitative degeneracies between the mass loss rate in the wind and parameters characterizing the microphysics of the accretion flow; of particular importance is $\delta$, the fraction of the turbulent energy which heats the electrons. For small $\delta$, current observations of soft X-ray transients and Sgr A* suggest that at least $\sim 10 %$ of the mass originating at large radii must reach the central object. For large $\delta \sim 0.3$, however, models with significantly more mass loss are in agreement with the observations. We also discuss constraints on advection-dominated accretion flow models imposed by recent radio observations of NGC 4649 and other nearby elliptical galaxies. We conclude by highlighting future observations which may clarify the importance of mass loss in sub-Eddington accretion flows.

Using mean field MHD, we discuss the energetics of optically thin, two temperature, advection-dominated accretion flows (ADAFs). If the magnetic field is tangled and roughly isotropic, flux freezing is insufficient to maintain the field in equipartition with the gas. In this case, we expect a fraction of the energy generated by shear in the flow to be used to build up the magnetic field strength as the gas flows in; the remaining energy heats the particles. We argue that strictly equipartition magnetic fields are incompatible with a priori reasonable levels of particle heating; instead, the plasma $\beta$ in ADAFs (defined to be the gas pressure divided by magnetic/turbulent pressure) is likely to be $\gsim 5$; correspondingly, the viscosity parameter $\alpha$ is likely to be $\lsim 0.2$

We consider the implications of recent infrared and radio observations of the nucleus of NGC 4258. There is no direct evidence that the nucleus has been steadily accreting on the viscous timescale of the outer masing disk, which is >~ 10^9 \yr. Thus the mass accretion rate in the outer disk need not be the same as in the inner accretion flow where most of the gravitational binding energy is released. We show that an advection-dominated flow model with a transition radius of \sim (10-100) M (where M is the mass of the hole) and \mdot \approx 10^{-2}\msun \yr^{-1} is consistent with the observed spectrum from radio to X-rays. We also show that a thin (flat or warped) disk can fit the observed fluxes outside the X-ray band. The X-rays can be explained by means of a corona in such a model, but the absence of radio emission from the location of the putative central black hole provides a serious constraint on the properties of the corona. A wide range of accretion rates, 10^{-4} <~ \mdot <~ 10^{-2} \msun\yr^{-1}, can be made to fit the data, but the most ``natural'' models have 10^{-3} <~ \mdot <~ 10^{-2} \msun\yr^{-1}. Physical conditions in the observed VLBI jet features can also be related to conditions in the inner accretion flow. We conclude with a list of future observations that might help to constrain the accretion rate.

This article reviews the physics of advection-dominated accretion flows (ADAFs) and describes applications to several black hole X-ray binaries and galactic nuclei. The possibility of using ADAFs to explore the event horizons of black holes is highlighted.

It is often stated that gamma-ray bursts (GRBs) have typical energies of several hundreds $\keV$, where the typical energy may be characterized by the hardness H, the photon energy corresponding to the peak of $\nu F_{\nu}$. Among the 54 BATSE bursts analyzed by Band et al. (1993), and 136 analyzed by us, more then 60% have 50 keV < H < 300 keV. Is the narrow range of H a real feature of GRBs or is it due to an observational difficulty to detect harder bursts? We consider a population of standard candle bursts with a hardness distribution: rho(H) d log H \propto H^gamma d log H and no luminosity - hardness correlation. We model the detection algorithm of BATSE as a function of H, including cosmological effects, detector characteristics and triggering procedure, and we calculate the expected distribution of H in the observed sample for various values of gamma. Both samples shows a paucity of soft (X-ray) bursts, which may be real. However, we find that the observed samples are consistent with a distribution above H=120 keV with gamma \sim -0.5 (a slowly decreasing numbers of GRBs per decade of hardness). Thus, we suggest that a large population of unobserved hard gamma-ray bursts may exist.

Astronomers have discovered many potential black holes in X-ray binaries and galactic nuclei. These black holes are usually identified by the fact that they are too massive to be neutron stars. Until recently, however, there was no convincing evidence that the objects identified as black hole candidates actually have event horizons. This has changed with extensive applications of a class of accretion models for describing the flow of gas onto compact objects; for these solutions, called advection-dominated accretion flows (ADAFs), the black hole nature of the accreting star, specifically its event horizon, plays an important role. We review the evidence that, at low luminosities, accreting black holes in both X-ray binaries and galactic nuclei contain ADAFs rather than the standard thin accretion disk.

The recently discovered GRB afterglow is believed to be described reasonably well by synchrotron emission from a slowing down relativistic shell that collides with an external medium. To compare theoretical models with afterglow observations we calculate here the broad band spectrum and corresponding light curve of synchrotron radiation from a power-law distribution of electrons in an expanding relativistic shock. Both the spectrum and the light curve consist of several power-law segments. The light curve is constructed under two limiting models for the hydrodynamical evolution of the shock: fully adiabatic and fully radiative. We compare the results with observations of $\gamma$-ray burst afterglows.

We show that the model proposed by Esin, McClintock & Narayan (1997) for the low state, intermediate state and high state of the black hole soft X-ray transient, Nova Muscae 1991, is consistent with the spectral evolution of the black hole X-ray binary, Cyg X-1, during the hard-to-soft state transition observed in 1996. We also apply the model to the outbursts of two other black hole X-ray transients, GRO J0422+32 and GRO J1719-24.

These lectures give an introduction to Gravitational Lensing. We discuss lensing by point masses, lensing by galaxies, and lensing by clusters and larger-scale structures in the Universe. The relevant theory is developed and applications to astrophysical problems are discussed.

Sgr A* at the Galactic Center is a puzzling source. It has a mass M=(2.5+/-0.4) x 10^6 solar masses which makes it an excellent black hole candidate. Observations of stellar winds and other gas flows in its vicinity suggest a mass accretion rate approximately few x 10^{-6} solar masses per year. However, such an accretion rate would imply a luminosity > 10^{40} erg/s if the radiative efficiency is the usual 10 percent, whereas observations indicate a bolometric luminosity <10^{37} erg/s. The spectrum of Sgr A* is unusual, with emission extending over many decades of wavelength. We present a model of Sgr A* which is based on a two-temperature optically-thin advection-dominated accretion flow. The model is consistent with the estimated mass and accretion rate, and fits the observed fluxes in the cm/mm and X-ray bands as well as upper limits in the sub-mm and infrared bands; the fit is less good in the radio below 86 GHz and in gamma-rays above 100 MeV. The very low luminosity of Sgr A* is explained naturally in the model by means of advection. Most of the viscously dissipated energy is advected into the central mass by the accreting gas, and therefore the radiative efficiency is extremely low, approximately 5 x 10^{-6}. A critical element of the model is the presence of an event horizon at the center which swallows the advected energy. The success of the model could thus be viewed as confirmation that Sgr A* is a black hole.

Recently, Narayan, Garcia and McClintock showed that there is a difference in the quiescent luminosities of black hole and neutron star x-ray nova: when the quiescent, Eddington scaled luminosities are compared, every black hole x-ray nova (BHXN) is fainter than every the neutron star x-ray nova (NSXN). To this earlier work we add three more BHXN, GRO J1655-40, GRO J0422+32, and 4U1543-47. This increases the significance of the difference between the quiescent luminosities to the >99% level. Given the reasonable assumption that the (Eddington scaled) mass transfer swings in BHXN and NSXN are similar, this difference provides direct evidence for the existence of event horizons in the BHXN: they have larger luminosity swings because the quiescent mass transfer energy is being hidden behind the event horizon.

Combined ASCA and SIGMA data of 1E1740.7-2942 during the standard state (September 1993 and 1994) were fitted with two-phase models (ISMBB/Poutanen at al. and ADAF/Narayan et al.). In ISMBB's, the radius of the spherical hot (Te = 150 - 200 KeV) corona lies between 200 - 250 km where it joins the classical inner disc. The disc radiates 40 per cent of the total luminosity with 1.7 per cent of the Eddington accretion rate of 10 solar mass black hole. ADAF's need an extra component to reproduce the soft part of the spectrum. However, the origin of the soft excess remains somewhat uncertain, although special care was taken in the background elimination.

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

We present a self-consistent model of accretion flows which unifies four distinct spectral states observed in black hole X-ray binaries: quiescent, low, intermediate and high states. In the quiescent, low and intermediate states, the flow consists of an inner hot advection-dominated part extending from the black hole horizon to a transition radius and an outer thin disk. In the high state the thin disk is present at all radii. The model is essentially parameter-free and treats consistently the dynamics of the accretion flow, the thermal balance of the ions and electrons, and the radiation processes in the accreting gas. With increasing mass accretion rate, the model goes through a sequence of stages for which the computed spectra resemble very well observations of the four spectral states; in particular, the low-to-high state transition observed in black hole binaries is naturally explained as resulting from a decrease in the transition radius. We also make a tentative proposal for the very high state, but this aspect of the model is less secure. We test the model against observations of the soft X-ray transient Nova Muscae during its 1991 outburst. The model reproduces the observed lightcurves and spectra surprisingly well, and makes a number of predictions which can be tested with future observations.

High Li abundances have been reported in the late type secondaries of five soft X-ray transients (SXTs), V404 Cyg, A0620-00, GS2000+25, Nova Mus 1991, and Cen X-4. Since Li is likely to be depleted in stars of this type, the origin of the Li is puzzling. Li has not been seen in similar secondaries of cataclysmic variables, which suggests that the high Li abundance is not due to an anomalous suppression of Li depletion in close binaries. SXTs in the quiescent state have hot advection-dominated accretion flows (ADAFs) in which the ions are essentially at virial temperature. At such temperatures, Li production via $\alpha-\alpha$ spallation is possible. We show that quiescent SXTs can produce sufficient Li via spallation to explain the observations in V404 Cyg, A0620-00, GS2000+25, and Nova Mus 1991. Depending on the Li depletion time scale in the secondary, which may range between $10^7-10^9$ yr, the model requires $\sim 10^{-4}-10^{-6}$ of the accreted mass to be intercepted by the secondary after undergoing Li production and being ejected. In the case of Cen X-4, we can explain the observed Li only if the mass accretion rate is $\sim 10^{-3}$ times the Eddington rate and if there is enhanced ejection due to a propeller effect. We discuss possible observational tests of this proposal. Li production during outbursts could be quite important and may even dominate over the production during quiescence, but the estimate of the Li yield is uncertain. We calculate the expected luminosity in gamma-ray lines due to the production of excited Li and Be nuclei, but conclude that the line cannot be detected with current instruments.

We show that the time delay between the optical and X-ray outbursts of the black-hole soft X-ray transient source GRO J1655-40, observed in April 1996, requires that the accretion flow in this object must consist of two components: a cold outer accretion disk and an extremely hot inner advection-dominated accretion flow (ADAF). In quiescence, the model predicts a spectrum which is in good agreement with observations, with most of the observed radiation coming from the ADAF. By fitting the observed spectrum, we estimate the mass accretion rate of the quiescent system and the transition radius between the disk and the ADAF. We present a detailed numerical simulation of a dwarf-nova type instability in the outer disk. The resulting heat front reaches the ADAF cavity promptly; however, it must then propagate inward slowly on a viscous time scale, thereby delaying the onset of the X-ray flux. The model reproduces the observed optical and X-ray light curves of the April 1996 outburst, as well as the 6-day X-ray delay. Further, the model gives an independent estimate of the quiescent mass accretion rate which is in very good agreement with the rate estimated from fitting the quiescent spectrum. We show that a pure thin disk model without an ADAF zone requires significant tuning to explain the X-ray delay; moreover, such a model does not explain the quiescent X-ray emission of GRO J1655-40.

We present a model for the high-low state transition of the X-Ray Nova GS 1124-68 (Nova Muscae 1991) observed by Ginga. The model consists of an advection-dominated accretion flow (ADAF) near the central black hole surrounded by a thin accretion disk. During the rise phase of the outburst, as the mass accretion rate increases, the transition radius between the thin disk and the ADAF moves closer to the center, until the thin disk extends all the way in to the last stable orbit. The transition radius increases again during decline. We reproduce the basic features of the spectra taken during and after outburst, and the light curves in the soft and hard X-ray bands. We estimate that the accretion rate in Nova Muscae decreased exponentially with a time scale $\sim$ 95 days during decline.

The defining characteristic of a black hole is that it possesses an event horizon through which matter and energy can fall in but from which nothing escapes. Soft X-ray transients (SXTs), a class of X-ray binaries, appear to confirm this fundamental property of black holes. SXTs that are thought to contain accreting black holes display a large variation of luminosity between their bright and faint states, while SXTs with accreting neutron stars have a smaller variation. This difference is predicted if the former stars have horizons and the latter have normal surfaces.

We have analyzed archival ASCA data on the soft X-ray transient source V404 Cyg in quiescence. We find that in the energy range 0.7 to 8.5 keV the spectrum is a hard power-law with a photon spectral index between 1.8 and 2.6 (90% confidence limits). We present a model of V404 Cyg in which the accretion flow has two components: (1) an outer thin disk with a small annular extent, and (2) a large interior region where the flow is advection-dominated. Nearly all the radiation in the infrared, optical, UV and X-ray bands is from the advection-dominated zone; the thin disk radiates primarily in the infrared where it contributes about ten percent of the observed flux. The spectrum we calculate with this model is in excellent agreement with the ASCA X-ray data presented here, as well as with previous optical data. Moreover, the fit is very insensitive to the choice of parameters such as black hole mass, orbital inclination, viscosity coefficient $\alpha$, and magnetic field strength. We consider the success of the model to be strong support for the advection-dominated accretion paradigm, and further evidence of the black hole nature of V404 Cyg. We discuss strategies whereby systems with advection-dominated accretion could be used to prove the reality of event horizons in black holes.

Recent work on advection-dominated accretion flows (ADAFs) is reviewed. The article concentrates on an optically thin branch of ADAFs which is present at mass accretion rates below a critical value $\sim(10^{-2}-10^{-1})$ the Eddington rate. Models based on this branch have been quite successful at explaining a number of low-luminosity X-ray binaries and galactic nuclei, and some brighter systems. Some progress has also been made toward understanding the various spectral states of accreting black holes. It is argued that ADAFs may provide one of the best techniques for demonstrating the reality of event horizons in black holes.

We present global solutions that describe advection-dominated accretion flows around black holes. The solutions are obtained by numerically solving a set of coupled ordinary differential equations corresponding to a steady axisymmetric height-integrated flow. The solutions satisfy consistent boundary conditions at both ends. We obtain well-behaved transonic solutions for a wide range of values of the viscosity parameter $\alpha$, from 0.001 to 0.3. We do not find any need for shocks in our solutions, and disagree with previous claims that viscous accretion flows with low values of $\alpha$ must have shocks. We see a qualitative difference between solutions with low values of the viscosity parameter, $\alpha\ \sles\ 0.01$, and those with large values, $\alpha\sgreat\ 0.01$. The solutions with low $\alpha$ have their sonic transitions occurring close to the radius of the marginally bound orbit. These flows are characterized by regions of super-Keplerian rotation, and have pressure maxima outside the sonic point. The solutions are similar in many respects to the hydrostatic thick tori developed previously as models of active galactic nuclei. In contrast, the solutions with large $\alpha$ have sonic transitions farther out, close to or beyond the marginally stable orbit, and have no super-Keplerian rotation or pressure maxima. We believe these flows will be nearly quasi-spherical down to the sonic radius and will not have empty funnels along the rotation axis. The large $\alpha$ solutions are more likely to be representative of real systems.

We present solutions for the accretion disks and boundary layers in pre-main-sequence stars undergoing FU Orionis outbursts. These solutions differ from earlier disk solutions in that they include a self-consistent treatment of the boundary layer region. In a previous paper (Popham 1996), we showed that these stars should stop accreting angular momentum once they spin up to modest rotation rates. Here we show that for reasonable values of $\alpha$, these low angular momentum accretion rate solutions fit the spectra and line profiles observed in FU Orionis objects better than solutions with high rates of angular momentum accretion. We find solutions which fit the observations of FU Orionis and V1057 Cygni. These solutions have mass accretion rates of 2 and $1 \times 10^{-4} \msyr$, stellar masses of 0.7 and $0.5 \msun$, and stellar radii of 5.75 and $5.03 \rsun$, respectively. They also have modest stellar rotation rates $8-9 \xss$, comparable to the observed rotation rates of T Tauri stars, and angular momentum accretion rates of zero. This supports our earlier suggestion that FU Orionis outbursts may regulate the rotation rates of T Tauri stars.

We show that the Keplerian thin disk in quiescent Soft X-ray Transients cannot extend down to the last stable orbit around the central black hole. We analyse the properties of the Narayan, McClintock & Yi (1996) model of quiescent Soft X-ray Transients in which the cold Keplerian disk has its inner edge at a large transition radius and transforms to a hot, advection-dominated flow on the inside. We show that outbursts of transient sources could be triggered in this model either by a pure thermal accretion disk instability or by a disk instability generated by an enhanced mass transfer from the stellar companion. Both mechanisms operate in the outer thin disk and could be at work, either in different systems or in the same system at different epochs, depending on the mass transfer rate and the value of the viscosity parameter $\alpha_t$ of the thin disk. We show that the recurrence time between outbursts in SXTs can be explained with values of $\alpha_t$ similar to these required by the dwarf nova disk instability model instead of the unreasonably low values needed in the model in which the thin disks extends down to the last stable orbit. We extend the Narayan, McClintock & Yi (1996) model to the case when the outer disk is non-stationary. We show that such disk is too cold to account for the observed UV flux. This difficulty is common to all models in which the outer disk is assumed to be optically thick.

A leading mechanism for producing cosmological gamma-ray bursts (GRBs) is via ultra-relativistic particles in an expanding fireball. The kinetic energy of the particles is converted into thermal energy in a forward shock and a reverse shock, when the outward flowing particles encounter the interstellar medium. The thermal energy is radiated via synchrotron emission and Comptonization. We estimate the synchrotron cooling time scale of the shocked regions for electrons of various Lorentz factors, focusing in particular on those electrons whose radiation falls within the energy detection range of the BATSE detectors. We find that in order to produce the rapid variability observed in most bursts the magnetic energy density in the shocked material must be greater than about 1\% of the thermal energy density. Additionally, the electrons must be nearly in equipartition with the protons, otherwise the radiative efficiency is low. Inverse Compton scattering can increase the cooling rate of the relevant electrons but the Comptonized emission itself is never within the BATSE range. These arguments allow us to pinpoint the conditions within the radiating regions in GRBs and to determine the important radiation processes. In addition, the model predicts that the duty cycle of intensity variations in GRB light curves should be nearly independent of burst duration, and should scale inversely as the square root of the observed photon energy. Both correlations are in agreement with observations. The model also provides a plausible explanation for the bimodal distribution of burst durations. There is no explanation, however, for the characteristic break energy in GRB spectra.

The double QSO2345+007 comprises two optical components separated by 7.1\ arcseconds and is the most prominent `dark matter' gravitational lens candidate. Despite being known for more than a decade, optical spectroscopy and imaging have been unable to determine whether this double QSO is a binary QSO or a gravitational lens system. In this note we report a deep VLA observation of this system, yielding a map with a noise level of 8.5\ts $\mu$Jy per beam. We have a $4\sigma$ detection of a radio source within one arcsecond of the optical position of the brighter A-component of the QSO, but no significant detection of any radio counterpart of the B component. Given that the flux ratio in the optical waveband is $\sim$3--4, the gravitational lens hypothesis for this system would predict a radio flux of image B of $\lesssim 10$\ts $\mu$Jy. The non-detection of the B component is thus consistent with, but does not prove, the lens interpretation.

We describe hot, optically-thin solutions for one-temperature accretion disks around black holes. We include cooling by synchrotron, bremsstrahlung, and Comptonization. Our solutions are thermally and viscously stable, with gas temperatures on the order of $T \sim 10^9-10^{10.7}$K. The thermal stability is a direct result of the inclusion of synchrotron cooling. The new solution branch is related to the advection-dominated solution for a two-temperature gas described by Narayan \& Yi (1995b). It is present only for mass accretion rates less than some critical $\dot{M}_{crit}$ which depends on the radius $R$ and viscosity parameter $\alpha$. The solutions are advection-dominated for extremely low values of $\dot{M}$. However, for a range of intermediate accretion rates, the new solutions are both hot ($T \sim 10^{10}$K) and cooling-dominated. Because of this new feature, one-temperature solutions are significantly more luminous than the corresponding two temperature solutions. The radial profile of the new solutions is unusual. The inner parts of the flow are cooling-dominated and have a disk-like geometry, while the outer parts are fully advection-dominated and nearly quasi-spherical.

We present a complete solution to the cyclotron-synchrotron radiation due to an isotropic distribution of electrons moving in a magnetic field. We make no approximations in the calculations other than artificially broadening the harmonics by a small amount in order to facilitate the numerics. In contrast to previous calculations, we sum over all relevant harmonics and integrate over all particle and observer angles relative to the magnetic field. We present emission spectra for electron temperatures $T=5\times10^8$ K, $10^9$ K, $2\times10^9$ K to $3.2\times10^{10}$ K, and provide simple fitting formulae which give a fairly accurate representation of the detailed results. For $T\geq3.2\times10^{10}$ K, the spectrum is represented well by the asymptotic synchrotron formula, which is obtained by assuming that the radiating electrons have Lorentz factors large compared to unity. We give an improved fitting formula also for this asymptotic case.

We present a novel method to recontruct the mass distribution of galaxy clusters from their gravitational lens effect on background galaxies. The method is based on a least-chisquare fit of the two-dimensional gravitational cluster potential. The method combines information from shear and magnification by the cluster lens and is designed to easily incorporate possible additional information. We describe the technique and demonstrate its feasibility with simulated data. Both the cluster morphology and the total cluster mass are well reproduced.

We present a simple mass model for the lensing galaxy in the gravitationally lensed quasar 0957+561. The model is a generalization of the singular isothermal sphere and includes a core radius, $r_c$, and a power-law index, $\eta$, defined such that mass increases as $r^\eta$ at large radius. We approximate the galaxy cluster surrounding the lensing galaxy with a quadratic potential described by its convergence $\kappa$ and shear $\gamma$. We fit the model to a recent high resolution VLBI map of the two images of 0957+561. We obtain a tight constraint on the radial index, $1.07 < \eta <1.18$, which means that the lens galaxy is nearly isothermal with increasing mass-to-light ratio out to $15 h^{-1}$ kpc. We also obtain an upper limit on the core radius, $r_c < 330 h^{-1}$ pc. We use the model to calculate the Hubble constant $H_0$ as a function of the time delay $\Delta\tau_{BA}$ between the two images: $H_0 = \left({ 82.5^{+5.9} _{-3.0} }\right) (1 - \kappa) \left({ \Delta\tau_{BA}/1.1 \,{\rm yr} }\right)^{-1}$ km/s/Mpc. Once $\Delta \tau_{BA}$ is measured, this will provide an upper bound on $H_0$ since $\kappa$ cannot be negative. In addition, the model degeneracy due to $\kappa$ can be eliminated if the one-dimensional velocity dispersion $\sigma$ of the lensing galaxy is measured. We then have $H_0 = \left({ 82.5^{+8.7} _{-5.6} }\right) (\sigma / 322\,{\rm km/s})^2 \left({ \Delta \tau_{BA} /1.1 \,{\rm yr} }\right)^{-1}$ km/s/Mpc. We investigate the effects of ellipticity in the lensing galaxy and clumpiness in the lensing cluster and find that these cause little change in our results.

It is often stated that gamma-ray bursts (GRBs) have typical energies of several hundred keV. Is this a real feature of GRBs or is it due to an observational bias? We consider the possibility that bursts of a given bolometric luminosity occur with a hardness distribution $p(H)d \log H \propto H^\gamma d \log H$. We model the detection efficiency of BATSE as a function of $H$ and calculate the expected distribution of $H$ in the observed sample for various values of $\gamma$. We show that because the detection efficiency of BATSE falls steeply with increasing $H$, the paucity of hard bursts need not be real. We find that the observed sample is consistent with a distribution above $H = 100$ keV with $\gamma \approx 0$ or even $\gamma =0.5$. Thus, a large population of unobserved hard gamma-ray bursts may exist. It is important to extend the present analysis to a larger sample of BATSE bursts and to include the OSSE and COMPTEL limits. If the full sample is consistent with $\gamma\ \sgreat\ 0$, then it would be interesting to look for MeV bursts in the future.

General properties of advection-dominated accretion flows are discussed. Special emphasis is given to the optically thin branch of solutions, which has very high ion and electron temperatures and is thermally stable. This solution branch has been applied to a number of low-luminosity accreting black holes. The models have resolved some puzzles and have provided a straightforward explanation of the observed spectra. The success of the models confirms that the central objects in these low-luminosity sources are black holes. There is some indication that advection-dominated models may be relevant also for higher luminosity systems. The properties of the Low state of accreting black holes, the transition from the Low state to the High state, and the similarity of hard X-ray/gamma-ray spectra of black hole X-ray binaries and active galactic nuclei, are explained.

We present a new algorithm for inverting poorly resolved gravitational lens systems using the maximum entropy method (MEM). We test the method with simulations and then apply it to an 8 GHz VLA map of the radio ring lens MG1654+134. We model the lens as a singular isothermal sphere embedded in an external shear field and find the critical radius of the lens is $b=0\parcs9820$, the dimensionless shear is $\gamma=0.0771$, and the position angle of the shear is $\theta=100\pdeg8$. These results are consistent with the results obtained by Kochanek (1995) using a complementary inversion algorithm based on Clean.

The BATSE experiment has now observed more than 1100 gamma-ray bursts. The observed angular distribution is isotropic, while the brightness distribution of bursts shows a reduced number of faint events. These observations favor a cosmological burst origin. Alternatively, very extended Galactic Halo (EGH) models have been considered. In the latter scenario, the currently favored source of gamma-ray bursts involves high velocity pulsars ejected from the Galactic disk. To be compatible with the observed isotropy, most models invoke a sampling distance of 300 kpc, a turn-on delay of 30 Myrs, and a source life time of about 1 Gyr. We consider the global energy requirements of such models and show that the largest known resource. rotational kinetic energy, is insufficient by orders of magnitude to provide the observed burst rate. More exotic energy sources or differently tuned pulsar models may be able to get around the global energy constraint but at the cost of becoming contrived. Thus, while extended halo models are not ruled out, our arguments place a severe obstacle for such models and we encourage proponents of EGH models to clearly address the issue of global energetics.

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

It has been found that a class of optically-thin two-temperature advection-dominated accretion solutions explains many observations of low-luminosity accreting black holes. Here it is shown that these models give a satisfactory description also of higher luminosity systems, provided the viscosity parameter $\alpha$ is large. The models reproduce the spectra of black hole X-ray binaries in the Low State, and explain the transition from the Low State to the High State at a critical mass accretion rate. The models also show that the X-ray/$\gamma$-ray spectra of X-ray binaries and active galactic nuclei should be similar, as confirmed by observations.

Viscous rotating accretion flows around black holes become advection-dominated when the accretion rate $\dot M$ is sufficiently low. Most of the accretion energy in such flows is stored within the gas and advected radially inward. The temperature is therefore very high, and much of the radiation comes out in hard X-rays and $\gamma$-rays. We have constructed an advection-dominated accretion flow model for the Galactic Center source Sagittarius A$^*$. The model consists of a $7\times10^5M_\odot$ black hole accreting at $\dot M=1.2\times10^{-5}\alpha\, M_\odot{\rm yr^{-1}}$, where $\alpha$ is the usual viscosity parameter. The model spectrum fits the observations from radio to $\gamma$-rays quite well and explains the unusually low luminosity of the source. Since the model explicitly makes use of a horizon at the inner edge to swallow the advected energy, the success of the model strongly suggests that the central object in Sgr A$^*$ is a black hole. We further show that, if $\alpha$ is not much smaller than unity, then advection-dominated models can be applied even to higher luminosity black holes. The existence of Low and High States in black hole X-ray binaries, and the abrupt transition between the two states, find a natural explanation. The models also explain the close similarity in the hard X-ray/$\gamma$-ray spectra of black hole X-ray binaries and active galactic nuclei.

When the mass accretion rate onto a black hole (BH) falls below a critical rate, ${\dot M}_{crit}\sim \alpha^2 {\dot M}_{Edd}$, accretion can occur via a hot optically thin flow where most of the dissipated energy is advected inward. We present such an advection-dominated model for the soft X-ray transient (SXT) A0620-00. This source has a puzzlingly low X-ray luminosity in quiescence, $\sim 6\times 10^{30} ~{\rm erg\,s^{-1}}$, despite a relatively high mass accretion rate $\sim10^{-10} ~M_{\odot}{\rm yr^{-1}}$ deduced from its optical flux. The accreting gas makes a transition from a standard thin disc at large radii to an advection-dominated flow at small radii. The transition occurs when the effective temperature of the thin disc is $\sim 10^4$ K. Because of the very low accretion efficiency, $\sim 10^{-3}-10^{-4}$, in the inner flow, the model fits both the optical and X-ray data. We also present models for V404 Cyg and Nova Mus 1991 in quiescence. Quiescent neutron star (NS) transients are expected to appear very different from BH systems because the advected energy is re-radiated from the NS surface whereas a BH swallows the advected energy. We discuss models for NS SXTs.

Most accreting black holes emit a substantial fraction of their luminosity in hard X-rays/soft $\gamma$-rays of up to $\sim{\rm few}\times100$ keV. The observations suggest that the radiating gas is very hot, $T>10^9$ K, and optically thin. Theoretical attempts to explain the observations are reviewed and unresolved issues are highlighted.

We present models of the soft X-ray transients, A0620-00, V404 Cyg, and X-ray Nova Mus 1991, in quiescence. In each source, we postulate that there is an outer region, extending outward from about 3000 Schwarzschild radii, where the accretion flow is in the form of a standard thin disk. The outer disk produces most of the radiation we observe in the infrared, optical and UV bands. We propose that the disk undergoes an instability at its inner edge, perhaps by the mechanism described recently by Meyer \& Meyer-Hofmeister for cataclysmic variables. The accreting gas is thereby converted into a hot corona which flows into the black hole as a nearly virial two-temperature flow. We describe the hot inner flow by means of a recently discovered stable solution of optically thin advection-dominated accretion. In this flow, most of the thermal energy released by viscous dissipation is advected into the black hole and only a small fraction, $\sim10^{-4}-10^{-3}$, of the energy is radiated. The radiation is in the form of Comptonized synchrotron and bremsstrahlung emission, and has a broad spectrum extending from optical to soft gamma-rays. The models we present are consistent with all the available data in the three sources. In particular, the X-ray emission from the hot inner flow fits the observed flux and spectral index of A0620-00. We derive a mass accretion rate of $\sim10^{-11}\msyr$ in A0620-00 and Nova Mus, and $\sim{\rm few}\times10^{-10}\msyr$ in V404 Cyg. The best fit to the data is obtained for a viscosity parameter $\alpha\sim0.1-0.3$ in the hot flow. The models predict that all three sources must have substantial flux in hard X-rays and soft $\gamma$-rays. This prediction is testable in the case of V404 Cyg with current instruments. A necessary feature of our proposal is that most of the viscous

We describe new optically thin solutions for rotating accretion flows around black holes and neutron stars. These solutions are advection dominated, so that most of the viscously dissipated energy is advected radially with the flow. We model the accreting gas as a two temperature plasma and include cooling by bremsstrahlung, synchrotron, and Comptonization. We obtain electron temperatures $T_e\sim 10^{8.5}-10^{10}$K. The new solutions are present only for mass accretion rates $\dot M$ less than a critical rate $\dot M_{crit}$ which we calculate as a function of radius $R$ and viscosity parameter $\alpha$. For $\dot M<\dot M_{crit}$ we show that there are three equilibrium branches of solutions. One of the branches corresponds to a cool optically thick flow which is the well-known thin disk solution of Shakura \& Sunyaev (1973). Another branch corresponds to a hot optically thin flow, discovered originally by Shapiro, Lightman \& Eardley (1976, SLE). This solution is thermally unstable. The third branch corresponds to our new advection-dominated solution. This solution is hotter and more optically thin than the SLE solution, but is viscously and thermally stable. It is related to the ion torus model of Rees et al. (1982) and may potentially explain the hard X-ray and $\gamma$-ray emission from X-ray binaries and active galactic nuclei.

We propose a new technique, which we call the lens parallax method, to determine simultaneously the redshift distribution of the faint blue galaxies and the mass distributions of foreground clusters of galaxies. The method is based on gravitational lensing and makes use of the following: (1) the amplitude of lensing-induced distortions of background galaxies increases with redshift; (2) the surface brightnesses of galaxies decrease steeply with redshift. The distortions of galaxy images due to lensing are thus expected to be inversely correlated with surface brightness, allowing us to obtain relative distances to galaxies as a function of surface brightness. If the redshifts of the brightest galaxies are measured, then the relative distance scale can be converted to mean galaxy redshifts as a function of surface brightness. Further, by comparing the angular sizes of lensed galaxies with those of similar galaxies in empty control fields, it is possible to break the so-called mass sheet degeneracy inherent to cluster mass reconstruction techniques which are based purely on image ellipticities. This allows an unambiguous determination of the surface density of a lensing cluster. We describe an iterative algorithm based on these ideas and present numerical simulations which show that the proposed techniques are feasible with a sample of ~ 10 rich clusters at moderate redshifts ~ 0.3-0.4 and an equal number of control fields. The numerical tests show that the method can be used to determine the redshifts of galaxies with an accuracy of dz ~ 0.1-0.2 at z ~ 1-1.7, and to measure the masses of lensing clusters to about 5% accuracy.

We provide a unified description of thermal equilibria of black hole accretion disks, including the newly-discovered advection-dominated solutions. We classify the solutions on the basis of optical depth and importance of advection cooling. We demonstrate that only four physically distinct topological types of equilibria exist. Two of the types correspond to optically thin and optically thick equilibria, while the other two types are distinguished by whether advection is negligible or dominant. A stable Shakura-Sunyaev disk exists only for accretion rates $\dot M$ below a certain maximum. However, there is a critical viscosity parameter $\alpha_{\rm crit}$, which is a function of radius, such that for $\alpha>\alpha_{\rm crit}$ advection-dominated solutions exist for all $\dot M$. Even when $\alpha< \alpha_{\rm crit}$, the advection-dominated solutions are available for a wide range of $\dot M$ except for a gap around the Eddington rate. We therefore suggest that advection-dominated flows may be more common than standard thin disks in black hole systems. For certain ranges of radii and $\dot M$, no stable steady state solution is possible. In these cases, we suggest that limit cycle behavior may occur, leading to variability.

We generate random Gaussian turbulent velocity fields with a Kolmogorov spectrum and use these to obtain synthetic line-of-sight velocity profiles. The profiles are found to be similar to line profiles observed in molecular clouds. We suggest methods for analysing measured line profiles to test whether they might arise from Gaussian Kolmogorov turbulence.

The enigmatic radio source Sagittarius A* at the centre of our Galaxy appears to be a low-luminosity version of active galactic nuclei in other galaxies. By analogy with active galactic nuclei models, it has been proposed that Sgr A* may be a massive accreting black holes. The black hole hypothesis is, however, problematic because no model of Sgr A* has been able to explain the observed spectrum in any self-consistent way, and there is no consensus on either the mass of the black hole or the mass accretion rate. Sgr A* has been observed in the radio, sub-mm, infrared, and X-ray bands, and the various various detections and flux upper limits cover more than ten decades of frequency, from $\nu\sles\ 10^9$ Hz (meter wavelength radio bands) up to $\nu\ \sgreat\ 10^{19}$ Hz (150 keV X-rays). We present a robust model of Sgr A* in which a million solar mass black hole accretes at the rate of a few $\times10^{-6}$ solar masses per year. This model fits the entire spectrum self-consistently. The unique feature of the model is that the flow is advection-dominated, i.e. most of the energy viscously dissipated in the differentially rotating flow is carried along with the gas and lost through the horizon. The apparent success of this model in explaining the data may be considered ``proof'' that horizons are real and that a massive black hole does exist at the Galactic center.

We consider axisymmetric viscous accretion flows where a fraction f of the viscously dissipated energy is advected with the accreting gas as stored entropy and a fraction 1-f is radiated. When f is small (i.e. very little advection), our solutions resemble standard thin disks in many respects except that they have a hot tenuous corona above. In the opposite {\it advection-dominated} limit ($f\rightarrow1$), the solutions approach nearly spherical accretion. The gas is almost at virial temperature, rotates at much below the Keplerian rate, and the flow is much more akin to Bondi accretion than to disk accretion. We compare our exact self-similar solutions with approximate solutions previously obtained using a height-integrated system of equations. We conclude that the height- integration approximation is excellent for a wide range of conditions. We find that the Bernoulli parameter is positive in all our solutions, especially close to the rotation axis. This effect is produced by viscous transport of energy from small to large radii and from the equator to the poles. In addition, all the solutions are convectively unstable and the convection is especially important near the rotation axis. For both reasons we suggest that a bipolar outflow will develop along the axis of the flows, fed by material from the the surface layers of the equatorial inflow.

We propose a new technique, which we call the lens parallax method, to determine simultaneously the redshift distribution of the faint blue galaxies and the mass distribution of foreground clusters of galaxies. The method is based on gravitational lensing and uses the fact that the mean angular size of the faint blue galaxies is a well-determined function of surface brightness. The method requires first a calibration of the angular-size vs. surface brightness relation in unlensed blank fields and a determination of the mean redshift of the brightest of the faint blue galaxies. By combining this information with observations of the distorted images of background galaxies in the fields of foreground clusters, we show that it is possible to obtain the mean redshift of the galaxies as a function of their surface brightness. With a sample of about ten moderate redshift rich clusters and using ten bins in surface brightness, a bin-to-bin signal-to-noise ratio of $\sim3.5$ can be achieved. The method simultaneously allows a determination of the convergence and shear of the lensing cluster as a function of position, and through this the mass distribution of the lens can be obtained. The lens parallax method can be used in conjunction with, and will improve the accuracy of, other previously proposed techniques for mapping clusters.

We review and discuss information on the mass distribution of galaxy clusters obtained from gravitational lensing.

We consider viscous rotating accretion flows in which most of the viscously dissipated energy is stored as entropy rather than being radiated. Such advection-dominated flows may occur when the optical depth is either very small or very large. We obtain a family of self-similar solutions where the temperature of the accreting gas is nearly virial and the flow is quasi-spherical. The gas rotates at much less than the Keplerian angular velocity; therefore, the central stars in such flows will cease to spin up long before they reach the break-up limit. Further, the Bernoulli parameter is positive, implying that advection-dominated flows are susceptible to producing outflows. Convection is likely in many of these flows and, if present, will tend to enhance the above effects. We suggest that advection-dominated accretion may provide an explanation for the slow spin rates of accreting stars and the widespread occurrence of outflows and jets in accreting systems.

It is well known that the standard transport equations violate causality when gradients are large or when temporal variations are rapid. We derive a modified set of transport equations that satisfy causality. These equations are obtained from the underlying Boltzmann equation. We use a simple model for particle collisions which enables us to derive moment equations non-perturbatively, i.e. without making the usual assumption that the distribution function deviates only slightly from its equilibrium value. We apply the model to two problems: particle diffusion and viscous transport. In both cases we show that signals propagate at a finite speed and therefore that the formalism obeys causality. When the velocity gradient is large on the scale of a mean free path, the viscous shear stress is suppressed relative to the prediction of the standard diffusion approximation. The shear stress reaches a maximum at a finite value of the shear amplitude and then decreases as the velocity gradient increases. In the case of a steady Keplerian accretion disk with hydrodynamic turbulent viscosity, the stress-limit translates to an upper bound on the Shakura-Sunyaev $\alpha$-parameter, namely $\alpha<0.07$. The limit on $\alpha$ is much stronger in narrow boundary layers where the velocity shear is larger than Keplerian.

Following the discovery by Quashnock and Lamb (1993) of an apparent excess of $\gamma$-ray burst pairs with small angular separations, we reanalyze the angular distribution of the bursts in the BATSE catalogue. We find that in addition to an excess of close pairs, there is also a comparable excess of antipodal bursts, i.e pairs of bursts separated by about 180 degrees in the sky. Both excesses have only modest statistical significance. We reject the hypothesis put forward by Quashnock and Lamb that burst sources are repeaters, since it is obvious that this hypothesis does not predict an excess of antipodal coincidences. Lacking any physical model of bursts that can explain the antipodal pairs, we suggest that the two excesses seen in the data are either due to an unusual statistical fluctuation or caused by some unknown selection effect.

Kouveliotou et al. (1993) recently confirmed that gamma-ray bursts are bimodal in duration. In this paper we compute the statistical properties of the short ($\le 2$~s) and long ($>2$~s) bursts using a method of analysis that makes no assumption regarding the location of the bursts, whether in the Galaxy or at a cosmological distance. We find the 64 ms channel on BATSE to be more sensitive to short bursts and the 1024 ms channel is more sensitive to long bursts. We show that all the currently available data are consistent with the simple hypothesis that both short and long bursts have the same spatial distribution and that within each population the sources are standard candles. The rate of short bursts is $\sim 0.4$ of the rate of long bursts. Although the durations of short and long gamma-ray bursts span several orders of magnitude and the total energy of a typical short burst is smaller than that of a typical long burst by a factor of $\sim 20$, surprisingly the peak luminosities of the two kinds of bursts are equal to within a factor of $\sim 2$.

We examine analytically and numerically the evolution of a relativistic fireball. We show that, after an early rearrangement phase, most of the matter and energy in the fireball is concentrated within a narrow shell. The shell propagates at nearly the speed of light, with a frozen radial profile and according to a simple set of scaling laws. The spectrum of the escaping radiation is harder at early times and softer later on. The results are applicable to models of $\gamma$-ray bursts. Submitted to Monthly Notices of the RAS Figures available on request from the authors

It is proposed that gamma-ray bursts are created in the mergers of double neutron star binaries and black hole neutron star binaries at cosmological distances. Bursts with complex profiles and relatively long durations are the result of magnetic flares generated by the Parker instability in a post-merger differentially-rotating disk. Some bursts may also be produced through neutrino-antineutrino annihilation into electrons and positrons. In both cases, an optically thick fireball of size $\sles\ 100$ km is initially created, which expands ultrarelativistically to large radii before radiating. Several previous objections to the cosmological merger model are eliminated. It is predicted that $\gamma$-ray bursts will be accompanied by a burst of gravitational radiation from the spiraling-in binary which could be detected by LIGO.