Person:

Narayan, Ramesh

A popular paradigm to explain the rapid temporal variability observed in gamma-ray burst (GRB) light curves 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 light curve are produced at roughly the same distance (R) from the centre 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 light curve is (\sim R/c \Gamma^2{\gamma^2}_t). The model naturally produces highly variable light curves with (\sim {\gamma^2}_t) individual pulses. Even though the model assumes highly inhomogeneous conditions, nevertheless the efficiency for converting jet energy to radiation is high.

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 (M_B) 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 outwards 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 (M_B). These solutions assume that the atmosphere beyond the Bondi radius rotates with a sub-Keplerian velocity and that the viscosity parameter is large, α≥ 0.1, both of which are reasonable for the problem at hand. 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 (M_{B}c^{2}) are expected if either a fraction of the accretion power is channelled 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. We discuss the Bernoulli parameter of the flow, the role of convection and the possibility that these as well as magnetohydrodynamic effects may invalidate the model. If the latter comes to pass, it would imply that the rough agreement between observed jet powers and the Bondi accretion rate is a coincidence and jet power is determined by factors other than the mass accretion rate.

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 * > 0.95 (3σ). For a less probable (synchronous) dynamical model, we find a * > 0.92 (3σ). 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 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.

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 γ-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.

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 pre-supernova 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 γ ∼ 400 and a half-opening angle θj ∼ 2◦, consistent with observations of many longduration GRBs. Other models with slightly different parameters give γ in the range 100–5000 and θj from 0. ◦1 to 10◦, 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 θ ∼ θj with the jet power concentrated in a hollow cone, while the maximum in the Lorentz factor occurs at an angle θ substantially smaller than θj also in a hollow cone. We derive approximate analytical formulae for the radial and angular distribution of γ and the radial dependence of θj . These formulae reproduce the simulation results and allow us to predict the outcome of models beyond those simulated. We also briefly discuss applications to active galactic nuclei, X-ray binaries and short-duration GRBs.

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.3 mm. 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 (. 1 hr) 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.

With an inferred bolometric luminosity exceeding (10^{42};erg;s^{–1}), 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;M_{☉} \lesssim M \lesssim 3 × 10^{5} M_{☉}), 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 modeling of the XMM-Newton data with more luminous Swift and Chandra observations increases the lower bound to (6000;M_☉), 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 \lesssim 10^{5} M_{☉}) if we limit (i \lesssim 75°).

Slim-disk models describe advective 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. The results also illustrate that advection alone does not solve the problem of the origin of the surprisingly soft high-luminosity spectra in LMC X-3. We show that X-ray continuum-fitting in the high accretion rate regime can powerfully test black-hole accretion disk models.

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.