Source author record

Mitchell C. Begelman

Mitchell C. Begelman appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.HE astro-ph.CO astro-ph physics.plasm-ph astro-ph.GA physics.space-ph

Catalog footprint

What is connected

64works

6topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

A radio, optical, UV and X-ray view of the enigmatic changing look Active Galactic Nucleus 1ES~1927+654 from its pre- to post-flare states

The nearby type-II AGN 1ES1927+654 went through a violent changing-look (CL) event beginning December 2017 during which the optical and UV fluxes increased by four magnitudes over a few months, and broad emission lines newly appeared in the optical/UV. By July 2018 the X-ray coronal emission had completely vanished, only to reappear a few months later. In this work we report the evolution of the radio, optical, UV and X-rays from the pre-flare state through mid-2021 with new and archival data from the Very Long Baseline Array (VLBA), the European VLBI Network, the Very Large Array (VLA), the Telescopio Nazionale Galileo (TNG), Gran Telescopio Canarias (GTC), The Neil Gehrels Swift observatory and XMM-Newton. The main results from our work are: (1) The source has returned to its pre-CL state in optical, UV, and X-ray; the disk-corona relation has been re-established as has been in the pre-CL state, with an $α_{\rm OX}\sim 1.02$. The optical spectra are dominated by narrow emission lines. (2) The UV light curve follows a shallower slope of $\propto t^{-0.91\pm 0.04}$ compared to that predicted by a tidal disruption event. We conjecture that a magnetic flux inversion event is the possible cause for this enigmatic event. (3) The compact radio emission which we tracked in the pre-CL (2014), during CL (2018) and post-CL(2021) at spatial scales $<1$ pc was at its lowest level during the changing look event in 2018, nearly contemporaneous with a low $2-10$ keV emission. The radio to X-ray ratio of the compact source $L_{\rm Radio}/L_{\rm X-ray}\sim 10^{-5.5}$, follows the Gudel-Benz relation, typically found in coronally active stars, and several AGN. (4) We do not detect any presence of nascent jets at the spatial scales of $\sim 5-10$ pc.

preprint2022arXiv

Kinetic Simulations of Instabilities and Particle Acceleration in Cylindrical Magnetized Relativistic Jets

Relativistic magnetized jets, such as those from AGN, GRBs and XRBs, are susceptible to current- and pressure-driven MHD instabilities that can lead to particle acceleration and non-thermal radiation. Here we investigate the development of these instabilities through 3D kinetic simulations of cylindrically symmetric equilibria involving toroidal magnetic fields with electron-positron pair plasma. Generalizing recent treatments by Alves et al. (2018) and Davelaar et al. (2020), we consider a range of initial structures in which the force due to toroidal magnetic field is balanced by a combination of forces due to axial magnetic field and gas pressure. We argue that the particle energy limit identified by Alves et al. (2018) is due to the finite duration of the fast magnetic dissipation phase. We find a rather minor role of electric fields parallel to the local magnetic fields in particle acceleration. In all investigated cases a kink mode arises in the central core region with a growth timescale consistent with the predictions of linearized MHD models. In the case of a gas-pressure-balanced (Z-pinch) profile, we identify a weak local pinch mode well outside the jet core. We argue that pressure-driven modes are important for relativistic jets, in regions where sufficient gas pressure is produced by other dissipation mechanisms.

preprint2022arXiv

Magnetization of Relativistic Current-Carrying Jets with Radial Velocity Shear

Astrophysical jets, launched from the immediate vicinity of accreting black holes, carry away large amounts of power in a form of bulk kinetic energy of jet particles and electromagnetic flux. Here we consider a simple analytical model for relativistic jets at larger distances from their launching sites, assuming a cylindrical axisymmetric geometry with a radial velocity shear, and purely toroidal magnetic field. We argue that, as long as the jet plasma is in magnetohydrostatic equilibrium, such outflows tend to be particle dominated, i.e. the ratio of the electromagnetic to particle energy flux, integrated over the jet cross-sectional area, is typically below unity, $σ< 1$. At the same time, for particular magnetic and radial velocity profiles, magnetic pressure may still dominate over particle pressure for certain ranges of the jet radius, i.e. the local jet plasma parameter $β_{pl} < 1$, and this may be relevant in the context of particle acceleration and production of high-energy emission in such systems. The jet magnetization parameter can be elevated up to the modest values $σ\lesssim \mathcal{O}(10)$ only in the case of extreme gradients or discontinuities in the gaseous pressure, and a significantly suppressed velocity shear. Such configurations, which consist of a narrow, unmagnetized jet spine surrounded by an extended, force-free layer, may require an additional poloidal field component to stabilize them against current-driven oscillations, but even this will not elevate substantially their $σ$ parameter.

preprint2022arXiv

Sgr A* X-ray flares from non-thermal particle acceleration in a magnetically arrested disc

Sgr A* exhibits flares in the near-infrared and X-ray bands, with the luminosity in these bands increasing by factors of $10-100$ for $\approx 60$ minutes. One of the models proposed to explain these flares is synchrotron emission of non-thermal particles accelerated by magnetic reconnection events in the accretion flow. We use the results from PIC simulations of magnetic reconnection to post-process 3D two-temperature GRMHD simulations of a magnetically arrested disc (MAD). We identify current sheets, retrieve their properties, estimate their potential to accelerate non-thermal particles and compute the expected non-thermal synchrotron emission. We find that the flux eruptions of MADs can provide suitable conditions for accelerating non-thermal particles to energies $γ_e \lesssim 10^6$ and producing simultaneous X-ray and near-infrared flares. For a suitable choice of current-sheet parameters and a simpified synchrotron cooling prescription, the model can simultaneously reproduce the quiescent and flaring X-ray luminosities as well as the X-ray spectral shape. While the near-infrared flares are mainly due to an increase in the temperature near the black hole during the MAD flux eruptions, the X-ray emission comes from narrow current sheets bordering highly magnetized, low-density regions near the black hole and equatorial current sheets where the flux on the black hole reconnects. As a result, not all infrared flares are accompanied by X-ray ones. The non-thermal flaring emission can extend to very hard ($\lesssim 100$ keV) X-ray energies.

preprint2022arXiv

What really makes an accretion disc MAD

Magnetically arrested accretion discs (MADs) around black holes (BH) have the potential to stimulate the production of powerful jets and account for recent ultra-high-resolution observations of BH environments. Their main properties are usually attributed to the accumulation of dynamically significant net magnetic (vertical) flux throughout the arrested region, which is then regulated by interchange instabilities. Here we propose instead that it is mainly a dynamically important {\it toroidal} field -- the result of dynamo action triggered by the significant but still relatively weak vertical field -- that defines and regulates the properties of MADs. We suggest that rapid convection-like instabilities, involving interchange of toroidal flux tubes and operating concurrently with the magnetorotatonal instability (MRI), can regulate the structure of the disc and the escape of net flux. We generalize the convective stability criteria and disc structure equations to include the effects of a strong toroidal field and show that convective flows could be driven towards two distinct marginally stable states, one of which we associate with MADs. We confirm the plausibility of our theoretical model by comparing its quantitative predictions to simulations of both MAD and SANE (strongly magnetized but not "arrested") discs, and suggest a set of criteria that could help to distinguish MADs from other accretion states. Contrary to previous claims in the literature, we argue that MRI is not suppressed in MADs and is probably responsible for the existence of the strong toroidal field.

preprint2020arXiv

First-principles demonstration of diffusive-advective particle acceleration in kinetic simulations of relativistic plasma turbulence

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

preprint2020arXiv

Kinetic turbulence in shining pair plasma: intermittent beaming and thermalization by radiative cooling

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

preprint2020arXiv

Magnetic flux inversion in a peculiar changing look AGN

We argue that the changing-look event in the active galactic nucleus 1ES 1927+654, followed by a dip of 3 orders of magnitude in the X-ray luminosity, is controlled by a change in the accretion rate and an inversion of magnetic flux in a magnetically arrested disk (MAD). Before the changing-look event, strong magnetic flux on the black hole powers X-ray emission via the Blandford-Znajek process while the UV emission is produced by a radiatively inefficient magnetized disk. An advection event, bringing flux of the opposite polarity, propagates inward leading, first, to a rise in the UV/optical luminosity and, then, to a dip in the X-ray luminosity when it reaches the black hole. We estimate the timescale for magnetic flux advection and find that the observed timescale between the beginning of the changing-look event and the minimum in the X-ray luminosity, $\approx200$ days, is in agreement with the time needed to cancel the magnetic flux in a MAD extending to $\approx180\:r_g$. Although flux inversion events might be rare due to the large ratio of flux-to-mass that is needed, we argue that AGN showing an unusually high ratio of X-ray to UV luminosity are prime candidates for such events. We also suggest that similar events may lead to jet interruptions in radio-loud objects.

preprint2020arXiv

Strongly magnetized accretion disks: structure and accretion from global magnetohydrodynamic simulations

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

preprint2016arXiv

A spine-sheath model for strong-line blazars

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

preprint2016arXiv

Accretion disc dynamo activity in local simulations spanning weak-to-strong net vertical magnetic flux regimes

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

preprint2016arXiv

Far-UV Emission Properties of FR1 Radio Galaxies

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

preprint2016arXiv

Magnetically elevated accretion disks in active galactic nuclei: broad emission line regions and associated star formation

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

preprint2016arXiv

On the structure of tidally-disrupted stellar debris streams

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

preprint2016arXiv

Strongly magnetized accretion discs require poloidal flux

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

preprint2016arXiv

The radiation hydrodynamics of relativistic shear flows

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

preprint2015arXiv

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

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

preprint2015arXiv

Post-periapsis pancakes: sustenance for self-gravity in tidal disruption events

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

preprint2015arXiv

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

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

preprint2015arXiv

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

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

preprint2014arXiv

A mechanism for hysteresis in black hole binary state transitions

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

preprint2014arXiv

Accreting Black Holes

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

preprint2014arXiv

Constraining the Location of Gamma-Ray Flares in Luminous Blazars

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

preprint2014arXiv

Entropy Production in Relativistic Jet Boundary Layers

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

preprint2014arXiv

Gamma-ray flares in the Crab Nebula: A case of relativistic reconnection?

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

preprint2014arXiv

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

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

preprint2014arXiv

The general relativistic equations of radiation hydrodynamics in the viscous limit

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

preprint2014arXiv

Thermal stability of thin disk with magnetically driven winds

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

preprint2014arXiv

Three-dimensional relativistic pair plasma reconnection with radiative feedback in the Crab Nebula

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

preprint2013arXiv

Collimation and confinement of magnetic jets by external media

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

preprint2013arXiv

Hyperaccretion during tidal disruption events: weakly bound debris envelopes and jets

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

preprint2013arXiv

Jet Emission in Young Radio Sources: a Fermi-LAT Gamma-ray View

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

preprint2013arXiv

Magnetic flux paradigm for radio-loudness of AGN

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

preprint2013arXiv

Polarimetry and the High-Energy Emission Mechanisms in Quasar Jets. The Case of PKS 1136-135

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

preprint2013arXiv

Quantifying energetics and dissipation in magnetohydrodynamic turbulence

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

preprint2013arXiv

Simulations of particle acceleration beyond the classical synchrotron burnoff limit in magnetic reconnection: An explanation of the Crab flares

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

preprint2013arXiv

Spectral hardening as a viable alternative to disc truncation in black hole state transitions

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

preprint2013arXiv

Supermassive Black Hole Formation at High Redshifts via Direct Collapse: Physical Processes in the Early Stage

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

preprint2012arXiv

Beaming and rapid variability of high-energy radiation from relativistic pair plasma reconnection

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

preprint2012arXiv

Energetic Constraints on a Rapid Gamma-Ray Flare in PKS 1222+216

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

preprint2012arXiv

Extreme particle acceleration in magnetic reconnection layers. Application to the gamma-ray flares in the Crab Nebula

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

preprint2012arXiv

Force-feeding Black Holes

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

preprint2012arXiv

Local Simulations of Instabilities in Relativistic Jets I: Morphology and Energetics of the Current-Driven Instability

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

preprint2012arXiv

Magnetic Domination of Recollimation Boundary Layers in Relativistic Jets

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

preprint2012arXiv

Recollimation Boundary Layers in Relativistic Jets

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

preprint2012arXiv

The effect of poloidal velocity shear on the local development of current-driven instabilities

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

preprint2011arXiv

High-efficiency photospheric emission of long-duration gamma-ray burst jets: the effect of the viewing angle

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

preprint2011arXiv

Radiatively Inefficient Accretion: Breezes, Winds and Hyperaccretion

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

preprint2011arXiv

Reconnection-Powered Linear Accelerator and Gamma-Ray Flares in the Crab Nebula

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

preprint2010arXiv

Non-Thermal emission from the photospheres of Gamma-Ray Burst outflows. I: High frequency tails

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

preprint2010arXiv

Quasistars and the cosmic evolution of massive black holes

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

preprint2010arXiv

Radiative Properties of Reconnection-Powered Minijets in Blazars

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

preprint2010arXiv

Short-duration gamma-ray bursts from off-axis collapsars

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

preprint2010arXiv

The origin and propagation of variability in the outflows of long duration gamma-ray bursts

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

preprint2009arXiv

Evolution of supermassive stars as a pathway to black hole formation

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

preprint2009arXiv

Fast TeV variability from misaligned minijets in the jet of M87

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

preprint2008arXiv

Hot self-similar relativistic MHD flows

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

preprint2007arXiv

Temporal and Angular Properties of GRB Jets Emerging from Massive Stars

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

preprint2006arXiv

The physics of GRB jets and their interaction with the progenitor star

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

preprint2002arXiv

Heating, conduction and minimum temperatures in cooling flows

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

preprint1998arXiv

RXTE Observation of Cygnus X--1: III. Implications for Compton Corona and ADAF Models

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

preprint1998arXiv

RXTE Observation of Cygnus X-1: II. Timing Analysis

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

preprint1995arXiv

ON THE GEOMETRY OF THE X-RAY EMITTING REGION IN SEYFERT GALAXIES

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

preprint1993arXiv

Gamma-Ray Bursts from Blast Waves Around Galactic Nuetron Stars

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

Mitchell C. Begelman

What is connected

Connect this record

See the researcher in context

Building this map preview

64 published item(s)

A radio, optical, UV and X-ray view of the enigmatic changing look Active Galactic Nucleus 1ES~1927+654 from its pre- to post-flare states

Kinetic Simulations of Instabilities and Particle Acceleration in Cylindrical Magnetized Relativistic Jets

Magnetization of Relativistic Current-Carrying Jets with Radial Velocity Shear

Sgr A* X-ray flares from non-thermal particle acceleration in a magnetically arrested disc

What really makes an accretion disc MAD

First-principles demonstration of diffusive-advective particle acceleration in kinetic simulations of relativistic plasma turbulence

Kinetic turbulence in shining pair plasma: intermittent beaming and thermalization by radiative cooling

Magnetic flux inversion in a peculiar changing look AGN

Strongly magnetized accretion disks: structure and accretion from global magnetohydrodynamic simulations

A spine-sheath model for strong-line blazars

Accretion disc dynamo activity in local simulations spanning weak-to-strong net vertical magnetic flux regimes

Far-UV Emission Properties of FR1 Radio Galaxies

Magnetically elevated accretion disks in active galactic nuclei: broad emission line regions and associated star formation

On the structure of tidally-disrupted stellar debris streams

Strongly magnetized accretion discs require poloidal flux

The radiation hydrodynamics of relativistic shear flows

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

Post-periapsis pancakes: sustenance for self-gravity in tidal disruption events

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

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

A mechanism for hysteresis in black hole binary state transitions

Accreting Black Holes

Constraining the Location of Gamma-Ray Flares in Luminous Blazars

Entropy Production in Relativistic Jet Boundary Layers

Gamma-ray flares in the Crab Nebula: A case of relativistic reconnection?

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

The general relativistic equations of radiation hydrodynamics in the viscous limit

Thermal stability of thin disk with magnetically driven winds

Three-dimensional relativistic pair plasma reconnection with radiative feedback in the Crab Nebula

Collimation and confinement of magnetic jets by external media

Hyperaccretion during tidal disruption events: weakly bound debris envelopes and jets

Jet Emission in Young Radio Sources: a Fermi-LAT Gamma-ray View

Magnetic flux paradigm for radio-loudness of AGN

Polarimetry and the High-Energy Emission Mechanisms in Quasar Jets. The Case of PKS 1136-135

Quantifying energetics and dissipation in magnetohydrodynamic turbulence

Simulations of particle acceleration beyond the classical synchrotron burnoff limit in magnetic reconnection: An explanation of the Crab flares

Spectral hardening as a viable alternative to disc truncation in black hole state transitions

Supermassive Black Hole Formation at High Redshifts via Direct Collapse: Physical Processes in the Early Stage

Beaming and rapid variability of high-energy radiation from relativistic pair plasma reconnection

Energetic Constraints on a Rapid Gamma-Ray Flare in PKS 1222+216

Extreme particle acceleration in magnetic reconnection layers. Application to the gamma-ray flares in the Crab Nebula

Force-feeding Black Holes

Local Simulations of Instabilities in Relativistic Jets I: Morphology and Energetics of the Current-Driven Instability

Magnetic Domination of Recollimation Boundary Layers in Relativistic Jets

Recollimation Boundary Layers in Relativistic Jets

The effect of poloidal velocity shear on the local development of current-driven instabilities

High-efficiency photospheric emission of long-duration gamma-ray burst jets: the effect of the viewing angle

Radiatively Inefficient Accretion: Breezes, Winds and Hyperaccretion

Reconnection-Powered Linear Accelerator and Gamma-Ray Flares in the Crab Nebula

Non-Thermal emission from the photospheres of Gamma-Ray Burst outflows. I: High frequency tails

Quasistars and the cosmic evolution of massive black holes

Radiative Properties of Reconnection-Powered Minijets in Blazars

Short-duration gamma-ray bursts from off-axis collapsars

The origin and propagation of variability in the outflows of long duration gamma-ray bursts

Evolution of supermassive stars as a pathway to black hole formation

Fast TeV variability from misaligned minijets in the jet of M87

Hot self-similar relativistic MHD flows

Temporal and Angular Properties of GRB Jets Emerging from Massive Stars

The physics of GRB jets and their interaction with the progenitor star

Heating, conduction and minimum temperatures in cooling flows

RXTE Observation of Cygnus X--1: III. Implications for Compton Corona and ADAF Models

RXTE Observation of Cygnus X-1: II. Timing Analysis

ON THE GEOMETRY OF THE X-RAY EMITTING REGION IN SEYFERT GALAXIES

Gamma-Ray Bursts from Blast Waves Around Galactic Nuetron Stars