Source author record

Dimitrios Giannios

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

45works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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.

Published work

45 published item(s)

preprint2026arXiv

The effect of inverse Compton losses on particle acceleration in three-dimensional relativistic reconnection

Relativistic magnetic reconnection is a key mechanism for dissipating magnetic energy and accelerating particles in astrophysics. In the absence of radiative cooling, recent particle-in-cell (PIC) simulations have shown that high-energy particles gain most of their energy in the upstream region, during a short-lived "free phase" where they meander between the two sides of the layer; when they get captured/trapped by the downstream flux ropes, they undergo a "trapped phase", where no significant energization occurs. Here, we perform a suite of 3D PIC simulations of relativistic reconnection including inverse Compton (IC) losses in the weakly cooled regime in which the radiation-reaction-limited Lorentz factor $γ_{\rm rad}$ exceeds the magnetization $σ$. We show that electron cooling losses do not appreciably alter the reconnection rate, the structure of the layer, and the physics of particle acceleration in the free phase, so the spectrum of free electrons is $dN_{\rm free}/dγ\propto γ^{-1}$, as in the uncooled case. The spectrum of trapped electrons above the cooling break $γ_{\rm cool}$ (in the range $γ_{\rm cool}<γ<γ_{\rm rad}$) is $dN/dγ\propto γ^{-3}$, steeper than the scaling $dN/dγ\propto γ^{-2}$ of uncooled simulations. This confirms that no significant particle energization occurs during the trapped phase. Our results validate the model by arXiv:2302.12269 for particle acceleration in 3D relativistic reconnection, and imply that radiative emission models of reconnection-powered astrophysical sources should employ a two-zone structure, that differentiates between free, rapidly accelerating particles and trapped, passively cooling particles.

preprint2022arXiv

Balancing Turbulent Heating with Radiative Cooling in Blazars

Recently, particle in cell (PIC) simulations have shown that relativistic turbulence in collisionless plasmas can result in an equilibrium particle distribution function where turbulent heating is balanced by radiative cooling of electrons. Strongly magnetized plasmas are characterized by higher energy peaks and broader particle distributions. In relativistically moving astrophysical jets, it is believed that the flow is launched Poynting flux dominated and that the resulting magnetic instabilities may create a turbulent environment inside the jet, i.e., the regime of relativistic turbulence. In this paper, we extend previous PIC simulation results to larger values of plasma magnetization by linearly extrapolating the diffusion and advection coefficients relevant for the turbulent plasmas under consideration. We use these results to build a single zone turbulent jet model that is based on the global parameters of blazar emission region, and consistently calculate the particle distribution and resulting synchrotron and inverse Compton emission spectra. We then test our model by comparing its predictions with the broad-band quiescent emission spectra of a dozen blazars. Our results show good agreement with observations of low-synchrotron peaked (LSP) sources and find that LSPs are moderately Poynting flux dominated with magnetization $1\lesssim σ\lesssim 5$, have bulk Lorentz factor $Γ\sim 10-30$, and that the turbulent region is located at the edge, or just beyond, the broad line region (BLR). The turbulence is found to be driven at an area comparable to the jet cross section.

preprint2022arXiv

Black hole to photosphere: 3D GRMHD simulations of collapsars reveal wobbling and hybrid composition jets

Long-duration $γ$-ray bursts (GRBs) accompany the collapse of massive stars and carry information about the central engine. However, no 3D models have been able to follow these jets from their birth by a black-hole (BH) to the photosphere. We present the first such 3D general-relativity magnetohydrodynamic simulations, which span over 6 orders of magnitude in space and time. The collapsing stellar envelope forms an accretion disk, which drags inwardly the magnetic flux that accumulates around the BH, becomes dynamically important and launches bipolar jets. The jets reach the photosphere at $\sim10^{12}$ cm with an opening angle $θ_j\sim6^\circ$ and a Lorentz factor $Γ_j\lesssim 30$, unbinding $\gtrsim90\%$ of the star. We find that (i) the disk-jet system spontaneously develops misalignment relative to the BH rotational axis. As a result, the jet wobbles with an angle $θ_t\sim12^\circ$, which can naturally explain quiescent times in GRB lightcurves. The effective opening angle for detection $θ_j+θ_t$ suggests that the intrinsic GRB rate is lower by an order of magnitude than standard estimates. This suggests that successful GRBs may be rarer than currently thought and emerge in only $\sim 0.1\%$ of supernovae Ib/c, implying that jets are either not launched or choked inside most supernova Ib/c progenitors. (ii) The magnetic energy in the jet decreases due to mixing with the star, resulting in jets with a hybrid composition of magnetic and thermal components at the photosphere, where $\sim 10\%$ of the gas maintains magnetization $σ\gtrsim 0.1$. This indicates that both a photospheric component and reconnection may play a role in the prompt emission.

preprint2022arXiv

FRB-periodicity: mild pulsars in tight O/B-star binaries

Periodicities observed in two Fast Radio Burst (FRB) sources (16 days in FRB 180916.J0158+65 and 160 days in FRB 121102) are consistent with that of tight, stellar mass binary systems. In the case of FRB 180916.J0158+65 the primary is an early OB-type star with mass loss rate $\dot{M} \sim 10^{-8}- 10^{-7} M_\odot$ yr$^{-1}$, and the secondary a neutron star. The observed periodicity is not intrinsic to the FRB's source, but is due to the orbital phase-dependent modulation of the absorption conditions in the massive star's wind. The observed relatively narrow FRB activity window implies that the primary's wind dynamically dominates that of the pulsar, $η= L_{sd}/(\dot{M} v_w c) \leq 1$, where $L_{sd} $ is pulsar spin-down, $\dot{M}$ is the primary's wind mass loss rate and $v_w$ is its velocity. The condition $η\leq 1$ requires mildly powerful pulsar with $L_{sd} \lesssim 10^{37}$ erg $s^{-1}$. The observations are consistent with magnetically-powered radio emission originating in the magnetospheres of strongly magnetized neutron stars, the classical magnetars.

preprint2021arXiv

Blazar jets launched with similar energy per baryon, independently of their power

The most extreme active galactic nuclei (AGN) are the radio active ones whose relativistic jet propagates close to our line of sight. These objects were first classified according to their emission line features into flat-spectrum radio quasars (FSRQs) and BL Lacertae objects (BL Lacs). More recently, observations revealed a trend between these objects known as the \emph{blazar sequence}, along with an anti-correlation between the observed power and the frequency of the synchrotron peak. In the present work, we propose a fairly simple idea that could account for the whole blazar population: all jets are launched with similar energy per baryon, independently of their power. In the case of FSRQs, the most powerful jets, manage to accelerate to high bulk Lorentz factors, as observed in the radio. As a result, they have a rather modest magnetization in the emission region, resulting in magnetic reconnection injecting a steep particle energy distribution and, consequently, steep emission spectra in the $γ$-rays. For the weaker jets, namely BL Lacs, the opposite holds true; i.e., the jet does not achieve a very high bulk Lorentz factor, leading to more magnetic energy available for non-thermal particle acceleration, and harder emission spectra at frequencies $\gtrsim$ GeV. In this scenario, we recover all observable properties of blazars with our simulations, including the \emph{blazar sequence} for models with mild baryon loading ($50 \lesssim μ\lesssim 80$). This interpretation of the blazar population, therefore, tightly constrains the energy per baryon of blazar jets regardless of their accretion rate.

preprint2021arXiv

Multimessenger parameter estimation of GW170817: from jet structure to the Hubble constant

The electromagnetic radiation that followed the neutron star merger event GW170817 revealed that gamma-ray burst afterglows from jets misaligned with our line of sight exhibit a light curve with slowly rising flux. The slope of the rising light curve depends sensitively on the angle of the observer with respect to the jet axis, which is likely to be perpendicular to the merger plane of the neutron star binary. Therefore, the afterglow emission can be used to constrain the inclination of the merging system. Here, we calculate the gamma-ray burst afterglow emission based on the realistic jet structure derived from general-relativistic magnetohydrodynamical simulations of a black hole torus system for the central engine of the gamma-ray burst. Combined with gravitational wave parameter estimation, we fit the multi-epoch afterglow emission of GW170817. We show that with such a jet model, the observing angle can be tightly constrained by multi messenger observations. The best fit observing angle of GW170817 is $θ_{\rm v} = 0.38\pm 0.02$ rad. With such a constraint, we can break the degeneracy between inclination angle and luminosity distance in gravitational wave parameter estimation, and substantially increase the precision with which the Hubble constant is constrained by the standard siren method. Our estimation of the distance is $D_{\rm L}=43.4\pm 1\ \rm Mpc$ and the Hubble constant constraint is $69.5\pm 4\ \mathrm{km\ s^{-1}\ Mpc^{-1}}$. As a result, multimessenger observations of short-duration gamma-ray bursts, combined with a good theoretical understanding of the jet structure, can be powerful probes of cosmological parameters.

preprint2021arXiv

Radiation and Polarization Signatures from Magnetic Reconnection in Relativistic Jets--II. Connection with $γ$-rays

It is commonly believed that blazar jets are relativistic magnetized plasma outflows from supermassive black holes. One key question is how the jets dissipate magnetic energy to accelerate particles and drive powerful multi-wavelength flares. Relativistic magnetic reconnection has been proposed as the primary plasma physical process in the blazar emission region. Recent numerical simulations have shown strong acceleration of nonthermal particles that may lead to multi-wavelength flares. Nevertheless, previous works have not directly evaluated $γ$-ray signatures from first-principle simulations. In this paper, we employ combined particle-in-cell and polarized radiation transfer simulations to study multi-wavelength radiation and optical polarization signatures under the leptonic scenario from relativistic magnetic reconnection. We find harder-when-brighter trends in optical and {\it Fermi-LAT} $γ$-ray bands as well as closely correlated optical and $γ$-ray flares. The optical polarization angle swings are also accompanied by $γ$-ray flares with trivial time delays. Intriguingly, we find highly variable synchrotron self Compton signatures due to inhomogeneous particle distributions during plasmoid mergers. This feature may result in fast $γ$-ray flares or orphan $γ$-ray flares under the leptonic scenario, complementary to the frequently considered mini-jet scenario. It may also infer neutrino emission with low secondary synchrotron flux under the hadronic scenario, if plasmoid mergers can accelerate protons to very high energy.

preprint2021arXiv

Shock-powered radio precursors of neutron star mergers from accelerating relativistic binary winds

During the final stages of a compact object merger, if at least one of the binary components is a magnetized neutron star (NS), then its orbital motion substantially expands the NS's open magnetic flux -- and hence increases its wind luminosity -- relative to that of an isolated pulsar. As the binary orbit shrinks due to gravitational radiation, the power and speed of this binary-induced inspiral wind may (depending on pair loading) secularly increase, leading to self-interaction and internal shocks in the outflow beyond the binary orbit. The magnetized forward shock can generate coherent radio emission via the synchrotron maser process, resulting in an observable radio precursor to binary NS merger. We perform 1D relativistic hydrodynamical simulations of shock interaction in the accelerating binary NS wind, assuming that the inspiral wind efficiently converts its Poynting flux into bulk kinetic energy prior to the shock radius. This is combined with the shock maser spectrum from particle-in-cell simulations, to generate synthetic radio light curves. The precursor burst with a fluence of $\sim1$ Jy$\cdot$ms at $\sim$GHz frequencies lasts $\sim 1-500$ ms following the merger for a source at $\sim3$ Gpc ($B_{\rm d}/10^{12}$ G)$^{8/9}$, where $B_{\rm d}$ is the dipole field strength of the more strongly-magnetized star. Given an outflow geometry concentrated along the binary equatorial, the signal may be preferentially observable for high-inclination systems, i.e. those least likely to produce a detectable gamma-ray burst.

preprint2020arXiv

Inverse Compton Signatures of Gamma-Ray Burst Afterglows

The afterglow emission from gamma-ray bursts (GRBs) is believed to originate from a relativistic blast wave driven into the circumburst medium. Although the afterglow emission from radio up to X-ray frequencies is thought to originate from synchrotron radiation emitted by relativistic, non-thermal electrons accelerated by the blast wave, the origin of the emission at high energies (HE; $\gtrsim$~GeV) remains uncertain. The recent detection of sub-TeV emission from GRB~190114C by MAGIC raises further debate on what powers the very high-energy (VHE; $\gtrsim 300$GeV) emission. Here, we explore the inverse Compton scenario as a candidate for the HE and VHE emissions, considering two sources of seed photons for scattering: synchrotron photons from the blast wave (synchrotron self-Compton or SSC) and isotropic photon fields external to the blast wave (external Compton). For each case, we compute the multi-wavelength afterglow spectra and light curves. We find that SSC will dominate particle cooling and the GeV emission, unless a dense ambient infrared photon field, typical of star-forming regions, is present. Additionally, considering the extragalactic background light attenuation, we discuss the detectability of VHE afterglows by existing and future gamma-ray instruments for a wide range of model parameters. Studying GRB~190114C, we find that its afterglow emission in the \fermi-LAT band is synchrotron-dominated.The late-time \fermi-LAT measurement (i.e., $t\sim 10^4$~s), and the MAGIC observation also set an upper limit on the energy density of a putative external infrared photon field (i.e. $\lesssim 3\times 10^{-9}\,{\rm erg\,cm^{-3}}$), making the inverse Compton dominant in the sub-TeV energies.

preprint2020arXiv

Kink Instabilities In Relativistic Jets Can Drive Quasi-Periodic Radiation Signatures

Relativistic jets are highly collimated plasma outflows emerging from accreting black holes. They are launched with a significant amount of magnetic energy, which can be dissipated to accelerate nonthermal particles and give rise to electromagnetic radiation at larger scales. Kink instabilities can be an efficient mechanism to trigger dissipation of jet magnetic energy. While previous works have studied the conditions required for the growth of kink instabilities in relativistic jets, the radiation signatures of these instabilities have not been investigated in detail. In this paper, we aim to self-consistently study radiation and polarization signatures from kink instabilities in relativistic jets. We combine large-scale relativistic magnetohydrodynamic (RMHD) simulations with polarized radiation transfer of a magnetized jet, which emerges from the central engine and propagates through the surrounding medium. We observe that a localized region at the central spine of the jet exhibits the strongest kink instabilities, which we identify as the jet emission region. Very interestingly, we find quasi-periodic oscillation (QPO) signatures in the light curve from the emission region. Additionally, the polarization degree appears to be anti-correlated to flares in the light curves. Our analyses show that these QPO signatures are intrinsically driven by kink instabilities, where the period of the QPOs is associated to the kink growth time scale. The latter corresponds to weeks to months QPOs in blazars. The polarization signatures offer unique diagnostics for QPOs driven by kink instabilities.

preprint2020arXiv

Radiation Signatures From Striped Blazar Jet

Relativistic jets from supermassive black holes are among the most powerful and luminous astrophysical systems in Universe. We propose that the open magnetic field lines through the black hole, which drive a strongly magnetized jet, may have their polarity reversing over time scales related to the growth of the magneto-rotational dynamo in the disc, resulting in dissipative structures in the jet characterized by reversing toroidal field polarities, referred to as ``stripes''. Magnetic reconnection between the stripes dissipates the magnetic energy and powers jet acceleration. The striped jet model can explain the jet acceleration, large-scale jet emission, and blazar emission signatures consistently in a unified physical picture. Specifically, we find that the jet accelerates to the bulk Lorentz factor $Γ\gtrsim 10$ within one parsec distance from the central engine. The acceleration slows down but continues at larger distances, with intrinsic acceleration rate $\dotΓ/Γ$ between $0.0005~\rm{yr^{-1}}$ and $0.005~\rm{yr^{-1}}$ at tens of parsecs, which is in very good agreement with recent radio observations. Magnetic reconnection continuously accelerates nonthermal particles over large distances from the central engine, resulting in the core-shift effect and overall flat-to-inverted synchrotron spectrum. The large-scale spectral luminosity peak $ν_{\rm peak}$ is anti-proportional to the location of the peak of the dissipation, which is set by the minimal stripe width $l_{\rm min}$. The blazar zone is approximately at the same location. At this distance, the jet is moderately magnetized, with the comoving magnetic field strength and dissipation power consistent with typical leptonic blazar model parameters.

preprint2020arXiv

Ready, set, launch: time interval between BNS merger and short GRB jet formation

The joint detection of GW~170817/GRB 170817 confirmed the long-standing theory that binary neutron star mergers produce short gamma-ray burst (sGRB) jets that can successfully break out of the surrounding ejecta. At the same time, the association with a kilonova provided unprecedented information regarding the physical properties (such as masses and velocities) of the different ejecta constituents. Combining this knowledge with the observed luminosities and durations of cosmological sGRBs detected by the Burst Alert Telescope (BAT) onboard the Neil Gehrels Swift Observatory, we revisit the breakout conditions of sGRB jets. Assuming self-collimation of sGRB jets does not play a critical role, we find that the time interval between the binary merger and the launching of a typical sGRB jet is $\lesssim0.1$~s. We also show that for a fraction of at least $\sim 30\%$ of sGRBs, the usually adopted assumption of static ejecta is inconsistent with observations, even if the polar ejecta mass is an order of magnitude smaller than the one in GRB 170817. Our results disfavour magnetar central engines for powering cosmological sGRBs, limit the amount of energy deposited in the cocoon prior to breakout, and suggest that the observed delay of $\sim 1.$7~s in GW 170817 /GRB 170817 between the gravitational wave and $γ$-ray signals is likely dominated by the propagation time of the jet to the $γ$-ray production site.

preprint2020arXiv

The blazar sequence revised

We propose and test a fairly simple idea that could account for the blazar sequence: all jets are launched with similar energy per baryon, independently of their power. For instance, flat-spectrum radio quasars (FSRQs), the most powerful jets, manage to accelerate to high bulk Lorentz factor, as observed in the radio. As a result, the emission region will have a rather modest magnetization which will induce a steep particle spectra therein and a rather steep emission spectra in the gamma-rays; particularly in the \textit{Fermi}-LAT band. For the weaker jets, namely BL Lacertae objects (BL Lacs), the opposite holds true; i.e., the jet does not achieve a very high bulk Lorentz factor, leading to more magnetic energy available for non-thermal particle acceleration and harder emission spectra. Moreover, this model requires but a handful of parameters. By means of numerical simulations we have accomplished to reproduce the spectral energy distributions and light-curves from fiducial sources following the aforementioned model. With the a complete evolution of the broadband spectra we were able to study in detail the spectral features at any particular frequency band at any given stage. Finally numerical results are compared and contrasted with observations.

preprint2016arXiv

Blazar flares powered by plasmoids in relativistic reconnection

Powerful flares from blazars with short ($\sim$ min) variability timescales are challenging for current models of blazar emission. Here, we present a physically motivated ab initio model for blazar flares based on the results of recent particle-in-cell (PIC) simulations of relativistic magnetic reconnection. PIC simulations demonstrate that quasi-spherical plasmoids filled with high-energy particles and magnetic fields are a self-consistent by-product of the reconnection process. By coupling our PIC-based results (i.e., plasmoid growth, acceleration profile, particle and magnetic content) with a kinetic equation for the evolution of the electron distribution function we demonstrate that relativistic reconnection in blazar jets can produce powerful flares whose temporal and spectral properties are consistent with the observations. In particular, our model predicts correlated synchrotron and synchrotron self-Compton flares of duration of several hours--days powered by the largest and slowest moving plasmoids that form in the reconnection layer. Smaller and faster plasmoids produce flares of sub-hour duration with higher peak luminosities than those powered by the largest plasmoids. Yet, the observed fluence in both types of flares is similar. Multiple flares with a range of flux-doubling timescales (minutes to several hours) observed over a longer period of flaring activity (days or longer) may be used as a probe of the reconnection layer's orientation and the jet's magnetization. Our model shows that blazar flares are naturally expected as a result of magnetic reconnection in a magnetically-dominated jet.

preprint2016arXiv

Flares from Galactic centre pulsars: a new class of X-ray transients?

Despite intensive searches, the only pulsar within 0.1 pc of the central black hole in our Galaxy, Sgr A*, is a radio-loud magnetar. Since magnetars are rare among the Galactic neutron star population, and a large number of massive stars are already known in this region, the Galactic centre (GC) should harbor a large number of neutron stars. Population syntheses suggest several thousand neutron stars may be present in the GC. Many of these could be highly energetic millisecond pulsars which are also proposed to be responsible for the GC gamma-ray excess. We propose that the presence of a neutron star within 0.03~pc from Sgr~A* can be revealed by the shock interactions with the disk around the central black hole. As we demonstrate, these interactions result in observable transient non-thermal X-ray and gamma-ray emission over timescales of months, provided that the spin down luminosity of the neutron star is L_{sd}~10^{35} erg/s. Current limits on the population of normal and millisecond pulsars in the GC region suggest that a number of such pulsars are present with such luminosities.

preprint2016arXiv

GRB off-axis afterglows and the emission from the accompanying supernovae

Gamma-ray burst (GRB) afterglows are likely produced in the shock that is driven as the GRB jet interacts with the external medium. Long-duration GRBs are also associated with powerful supernovae (SNe). We consider the optical and radio afterglows of long GRBs for both blasts viewed along the jet axis ("on-axis" afterglows) and misaligned observes ("off-axis" afterglows). Comparing the optical emission from the afterglow with that of the accompanying SN, using SN 1998bw as an archetype, we find that only a few percent of afterglows viewed off-axis are brighter than the SN. For observable optical off-axis afterglows, the viewing angle is at most twice the half-opening angle of the GRB jet. Radio off-axis afterglows should be detected with upcoming radio surveys within a few hundred Mpc. We propose that these surveys will act as "radio triggers", and that dedicated radio facilities should follow-up these sources. Follow-ups can unveil the presence of the radio SN remnant, if present. In addition, they can probe the presence of a mildly relativistic component, either associated with the GRB jet or the SN ejecta, expected in these sources.

preprint2016arXiv

Radio SNRs in the Magellanic Clouds as probes of shock microphysics

A large number of supernova remnants (SNRs) in our Galaxy and galaxies nearby have been resolved in various radio bands. This radio emission is thought to be produced via synchrotron emission from electrons accelerated by the shock that the supernova ejecta drives into the external medium. Here we consider the sample of radio SNRs in the Magellanic Clouds. Given the size and radio flux of a SNR, we seek to constrain the fraction of shocked fluid energy in non-thermal electrons ($ε_e$) and magnetic field ($ε_B$), and find $ε_e ε_B \sim 10^{-3}$. These estimates do not depend on the largely uncertain values of the external density and the age of the SNR. We develop a Monte Carlo scheme that reproduces the observed distribution of radio fluxes and sizes of the population of radio SNRs in the Magellanic Clouds. This simple model provides a framework that could potentially be applied to other galaxies with complete radio SNRs samples.

preprint2016arXiv

The TeV emission of Ap Librae: a hadronic interpretation and prospects for CTA

Ap Librae is one out of a handful of low-frequency peaked blazars to be detected at TeV $γ$-rays and the only one with an identified X-ray jet. Combined observations of Fermi-LAT at high energies (HE) and of H.E.S.S at very high energies (VHE) revealed a striking spectral property of Ap Librae; the presence of a broad high-energy component that extends more than nine orders of magnitude in energy and is, therefore, hard to be explained by the usual single-zone synchrotron self-Compton model. We show that the superposition of different emission components related to photohadronic interactions can explain the $γ$-ray emission of Ap Librae without invoking external radiation fields. We present two indicative model fits to the spectral energy distribution of Ap Librae where the VHE emission is assumed to originate from a compact, sub-pc scale region of the jet. A robust prediction of our model is VHE flux variability on timescales similar to those observed at X-rays and HE $γ$-rays, which can be further used to distinguish between a sub-pc or kpc scale origin of the TeV emission. We thus calculate the expected variability signatures at X-rays, HE and VHE $γ$-rays and show that quasi-simultaneous flares are expected, with larger amplitude flares appearing at $γ$-rays. We assess the detectability of VHE variability from Ap Librae with CTA, next generation of IACTs. We show that $\sim$hr timescale variability at $E_γ>0.1$ TeV could be detectable at high significance with shorter exposure times than current Cherenkov telescopes.

preprint2015arXiv

An anisotropic minijets model for the GRB prompt emission

In order to explain rapid light curve variability without invoking a variable source, several authors have proposed "minijets" that move relativistically relative to the main flow of the jet. Here we consider the possibility that these minijets, instead of being isotropically distributed in the comoving frame of the jet, form primarily perpendicular to the direction of the flow, as the jet dissipates its energy at a large emission radius. This yields two robust features. First, the emission is significantly delayed compared with the isotropic case. This delay allows for the peak of the afterglow emission to appear while the source is still active, in contrast to the simplest isotropic model. Secondly, the flux decline after the source turns off is steeper than the isotropic case. We find that these two features are realized in gamma-ray bursts (GRBs): 1. The peak of most GeV light curves (ascribed to the external shock) appears during the prompt emission phase. 2. Many X-ray light curves exhibit a period of steep decay, which is faster than that predicted by the standard isotropic case. The gamma-ray generation mechanism in GRBs, and possibly in other relativistic flows, may therefore be anisotropic.

preprint2015arXiv

Radio rebrightening of the GRB afterglow by the accompanying supernova

The gamma-ray burst (GRB) jet powers the afterglow emission by shocking the surrounding medium, and radio afterglow can now be routinely observed to almost a year after the explosion. Long-duration GRBs are accompanied by supernovae (SNe) that typically contain much more energy than the GRB jet. Here we consider the fact that the SN blast wave will also produce its own afterglow (supernova remnant emission), which will peak at much later time (since it is non-relativistic), when the SN blast wave transitions from a coasting phase to a decelerating Sedov-Taylor phase. We predict that this component will peak generally a few tens of years after the explosion and it will outshine the GRB powered afterglow well-before its peak emission. In the case of GRB 030329, where the external density is constrained by the $\sim 10$-year coverage of the radio GRB afterglow, the radio emission is predicted to start rising over the next decade and to continue to increase for the following decades up to a level of $\sim$ mJy. Detection of the SN-powered radio emission will greatly advance our knowledge of particle acceleration in $ \sim 0.1$c shocks.

preprint2015arXiv

Relativistic Jets Shine through Shocks or Magnetic Reconnection?

Observations of gamma-ray-bursts and jets from active galactic nuclei reveal that the jet flow is characterized by a high radiative efficiency and that the dissipative mechanism must be a powerful accelerator of non-thermal particles. Shocks and magnetic reconnection have long been considered as possible candidates for powering the jet emission. Recent progress via fully-kinetic particle-in-cell simulations allows us to revisit this issue on firm physical grounds. We show that shock models are unlikely to account for the jet emission. In fact, when shocks are efficient at dissipating energy, they typically do not accelerate particles far beyond the thermal energy, and vice versa. In contrast, we show that magnetic reconnection can deposit more than 50% of the dissipated energy into non-thermal leptons as long as the energy density of the magnetic field in the bulk flow is larger than the rest mass energy density. The emitting region, i.e., the reconnection downstream, is characterized by a rough energy equipartition between magnetic fields and radiating particles, which naturally accounts for a commonly observed property of blazar jets.

preprint2014arXiv

Black-hole jets without large-scale net magnetic flux

We propose a scenario for launching relativistic jets from rotating black holes, in which small-scale magnetic flux loops, sustained by disc turbulence, are forced to inflate and open by differential rotation between the black hole and the accretion flow. This mechanism does not require a large-scale net magnetic flux in the accreting plasma. Estimates suggest that the process could operate effectively in many systems, and particularly naturally and efficiently when the accretion flow is retrograde. We present the results of general-relativistic force-free electrodynamic simulations demonstrating the time evolution of the black hole's magnetosphere, the cyclic formation of jets, and the effect of magnetic reconnection. The jets are highly variable on timescales ~ 10-10^3 r_ g/c, where r_g is the black hole's gravitational radius. The reconnecting current sheets observed in the simulations may be responsible for the hard X-ray emission from accreting black holes.

preprint2014arXiv

Hadronic supercriticality as a trigger for GRB emission

We explore a one-zone hadronic model that may be able to reproduce $γ$-ray burst (GRB) prompt emission with a minimum of free parameters. Assuming only that GRBs are efficient high-energy proton accelerators and without the presence of an ab initio photon field, we investigate the conditions under which the system becomes supercritical, i.e. there is a fast, non-linear transfer of energy from protons to secondary particles initiated by the spontaneous quenching of proton-produced $γ$-rays. We first show analytically that the transition to supercriticality occurs whenever the proton injection compactness exceeds a critical value, which favours high proton injection luminosities and a wide range of bulk Lorentz factors. The properties of supercriticality are then studied with a time-dependent numerical code that solves concurrently the coupled equations of proton, photon, electron, neutron and neutrino distributions. For conditions that drive the system deep into the supercriticality we find that the photon spectra obtain a Band-like shape due to Comptonization by cooled pairs and that the energy transfer efficiency from protons to $γ$-rays and neutrinos is high reaching $\sim 0.3$. Although some questions concerning its full adaptability to the GRB prompt emission remain open, supercriticality is found to be a promising process in that regard.

preprint2014arXiv

Implications of a PeV neutrino spectral cutoff in GRB models

The recent discovery of extragalactic PeV neutrinos opens a new window to the exploration of cosmic-ray accelerators. The observed PeV neutrino flux is close to the Waxman-Bahcall upper bound implying that gamma-ray bursts (GRBs) may be the source of ultra-high energy cosmic rays (UHECRs). Starting with the assumption of the GRB-UHECR connection, we show using both analytical estimates and numerical simulations that the observed neutrinos can originate at the jet as a result of photopion interactions with the following implications: the neutrino spectra are predicted to have a cutoff at energy $\le 10$ PeV; the dissipation responsible for the GRB emission and cosmic-ray acceleration takes place at distances $r_{\rm diss} \approx 3 \times 10^{11}-3 \times 10^{13} \ {\rm cm}$ from the central engine; the Thomson optical depth at the dissipation region is $τ_{\rm T} \sim 1$; the jet carries a substantial fraction of its energy in the form of Poynting flux at the dissipation region, and has a Lorentz factor $Γ\simeq 100-500$. The non-detection of $\sim$PeV neutrinos coincident with GRBs will indicate that GRBs are either poor cosmic accelerators or the dissipation takes place at small optical depths in the jet.

preprint2014arXiv

Magnetic Flux of Progenitor Stars Sets Gamma-ray Burst Luminosity and Variability

Long-duration gamma-ray bursts (GRBs) are thought to come from the core-collapse of Wolf-Rayet stars. Whereas their stellar masses $M_*$ have a rather narrow distribution, the population of GRBs is very diverse, with gamma-ray luminosities $L_γ$ spanning several orders of magnitude. This suggests the existence of a "hidden" stellar variable whose burst-to-burst variation leads to a spread in $L_γ$. Whatever this hidden variable is, its variation should not noticeably affect the shape of GRB lightcurves, which display a constant luminosity (in a time-average sense) followed by a sharp drop at the end of the burst seen with Swift/XRT. We argue that such a hidden variable is progenitor star's large-scale magnetic flux. Shortly after the core collapse, most of stellar magnetic flux accumulates near the black hole (BH) and remains there. The flux extracts BH rotational energy and powers jets of roughly a constant luminosity, $L_j$. However, once BH mass accretion rate $\dot M$ falls below $\sim L_j/c^2$, the flux becomes dynamically important and diffuses outwards, with the jet luminosity set by the rapidly declining mass accretion rate, $L_j\sim \dot M c^2$. This provides a potential explanation for the sharp end of GRBs and the universal shape of their lightcurves. During the GRB, gas infall translates spatial variation of stellar magnetic flux into temporal variation of $L_j$. We make use of the deviations from constancy in $L_j$ to perform stellar magnetic flux "tomography". Using this method, we infer the presence of magnetised tori in the outer layers of progenitor stars for GRB 920513 and GRB 940210.

preprint2014arXiv

Relativistic magnetic reconnection in pair plasmas and its astrophysical applications

This review discusses the physics of magnetic reconnection, a process in which the magnetic field topology changes and magnetic energy is converted to kinetic energy, in pair plasmas in the relativistic regime. We focus on recent progress in the field driven by theory advances and the maturity of particle-in-cell codes. This work shows that fragmentation instabilities at the current sheet can play a critical role in setting the reconnection speed and affect the resulting particle acceleration, anisotropy, bulk flows, and radiation. Then, we discuss how this novel understanding of relativistic reconnection can be applied to high-energy astrophysical phenomena, with an emphasis on pulsars, pulsar wind nebulae, and active galactic nucleus jets.

preprint2014arXiv

Testing the neutrino annihilation model for launching GRB jets

The mechanism behind the launching of gamma-ray-burst (GRB) jets remains debated resulting in large uncertainty over the jet composition. Both magnetohydrodynamical and neutrino annihilation models have been proposed for the energy extraction in a black hole/accretion-disc central engine. In particular, for the extreme accretion rates $\dot M\sim 0.1-1$~ M$_\odot$s$^{-1}$ expected for bursts of duration $T\lesssim 100$~s, the disc can be an efficient neutrino emitter. Neutrino-antineutrino annihilation results in an energy deposition rate at the jet that can, in principle, account for the burst's energetics. Recent discoveries of X-ray flares hours after the burst and of ultra-long GRBs suggest that GRB activity can last for $\sim 10^4$~s or longer. These long-lived events have fluence similar to that of classical GRBs. In view of these findings, we re-evaluate the neutrino annihilation model. We derive the maximum possible energy of a neutrino-powered jet as a function of the burst duration and show that the available energy drops fast for longer bursts. For a standard choice of the parameters, the model falls short by three to four orders of magnitude in explaining the observed energetics of events that last longer than $\sim 10^3$~s.

preprint2013arXiv

A Late-Time Flattening of Light Curves in Gamma-Ray Burst Afterglows

The afterglow emission from Gamma-Ray Bursts (GRBs) is usually interpreted as synchrotron radiation from relativistic electrons accelerated at the GRB external shock, that decelerates from ultra-relativistic to non-relativistic speeds as it sweeps up the surrounding medium. We investigate the temporal decay of the emission from GRB afterglows at late times, when the bulk of the shock-accelerated electrons are non-relativistic. For a uniform circumburst medium, we show that such "deep Newtonian phase" begins at t_{DN} ~ 3 epsilon_{e,-1}^{5/6} t_{ST}, where t_{ST} marks the transition of the blast wave to the non-relativistic spherically-symmetric Sedov-Taylor solution, and epsilon_e = 0.1 epsilon_{e,-1} quantifies the amount of shock energy transferred to the post-shock electrons. For typical parameters, the deep Newtonian stage starts ~ 0.5-several years after the GRB. The radio flux in this phase decays as F_nu ~ t^{-3(p+1)/10}, with decay slope between -0.9 and -1.2, for a power-law distribution of shock-accelerated electrons with index 2 < p < 3. This is shallower than the commonly assumed scaling F_nu ~ t^{-3(5p-7)/10} (with decay slope between -0.9 and -2.4) derived by Frail et al. (2000), which only applies if the GRB shock is non-relativistic, but the electron distribution still peaks at ultra-relativistic energies (a regime that is relevant for a narrow time interval, and only if t_{DN} > t_{ST}, or epsilon_e > 0.03). We discuss how the deep Newtonian phase can be reliably used for GRB calorimetry, and we comment on the good detection prospects of trans-relativistic blast waves at 0.1-10 GHz with EVLA and LOFAR.

preprint2013arXiv

Reconnection-driven plasmoids in blazars: fast flares on a slow envelope

Minute-timescale TeV flares have been observed in several blazars. The fast flaring requires compact regions in the jet that boost their emission towards the observer at an extreme Doppler factor of delta>50. For TeV photons to avoid annihilation in the broad line region of PKS 1222+216, the flares must come from large (pc) scales challenging most models proposed to explain them. Here we elaborate on the magnetic reconnection minijet model for the blazar flaring, focusing on the inherently time-dependent aspects of the process of magnetic reconnection. I argue that, for the physical conditions prevailing in blazar jets, the reconnection layer fragments leading to the formation a large number of plasmoids. Occasionally a plasmoid grows to become a large, "monster" plasmoid. I show that radiation emitted from the reconnection event can account for the observed "envelope" of ~ day-long blazar activity while radiation from monster plasmoids can power the fastest TeV flares. The model is applied to several blazars with observed fast flaring. The inferred distance of the dissipation zone from the black hole and the typical size of the reconnection regions are R_diss~0.3-1 pc and l' \simless 10^{16} cm, respectively. The required magnetization of the jet at this distance is modest: sigma ~a few. Such dissipation distance R_diss and reconnection size l' are expected if the jet contains field structures with size of the order of the black-hole horizon.

preprint2013arXiv

Relativistic Pair Beams from TeV Blazars: A Source of Reprocessed GeV Emission rather than IGM Heating

The interaction of TeV photons from blazars with the extragalactic background light produces a relativistic beam of electron-positron pairs streaming through the intergalactic medium (IGM). The fate of the beam energy is uncertain. By means of two- and three-dimensional particle-in-cell simulations, we study the non-linear evolution of dilute ultra-relativistic pair beams propagating through the IGM. We explore a wide range of beam Lorentz factors gamma_b>>1 and beam-to-plasma density ratios alpha<<1, so that our results can be extrapolated to the extreme parameters of blazar-induced beams (gamma_b~10^6 and alpha~10^(-15), for the most powerful blazars). For cold beams, we show that the oblique instability governs the early stages of evolution, but its exponential growth terminates - due to self-heating of the beam in the transverse direction - when only a negligible fraction ~(alpha/gamma_b)^(1/3)~10^(-7) of the beam energy has been transferred to the IGM plasma. Further relaxation of the beam proceeds through quasi-longitudinal modes, until the momentum dispersion in the direction of propagation saturates at DeltaP_{b,parallel}/gamma_b m_e c ~ 0.2. This corresponds to a fraction ~10% of the beam energy being ultimately transferred to the IGM plasma, irrespective of gamma_b or alpha. If the initial dispersion in beam momentum satisfies DeltaP_{b,parallel}/gamma_b m_e c > 0.2 (as typically expected for blazar-induced beams), the fraction of beam energy deposited into the IGM is much smaller than ~10%. It follows that at least ~90% of the beam energy is still available to power the GeV emission produced by inverse Compton up-scattering of the Cosmic Microwave Background by the beam pairs.

preprint2013arXiv

Swift J1644+57 gone MAD: the case for dynamically-important magnetic flux threading the black hole in a jetted tidal disruption event

The unusual transient Swift J1644+57 likely resulted from a collimated relativistic jet powered by accretion onto a massive black hole (BH) following the tidal disruption (TD) of a star. Several mysteries cloud the interpretation of this event: (1) extreme flaring and `plateau' shape of the X-ray/gamma-ray light curve during the first 10 days after the gamma-ray trigger; (2) unexpected rebrightening of the forward shock radio emission months after trigger; (3) no obvious evidence for jet precession, despite misalignment typically expected between the angular momentum of the accretion disk and BH; (4) recent abrupt shut-off in jet X-ray emission after 1.5 years. Here we show that all of these seemingly disparate mysteries are naturally resolved by one assumption: the presence of strong magnetic flux Phi threading the BH. Initially, Phi is weak relative to high fall-back mass accretion rate, Mdot, and the disk and jets precess about the BH axis = our line of sight. As Mdot drops, Phi becomes dynamically important and leads to a magnetically-arrested disk (MAD). MAD naturally aligns disk and jet axis along the BH spin axis, but only after a violent rearrangement phase (jet wobbling). This explains the erratic light curve at early times and the lack of precession at later times. We use our model for Swift J1644+57 to constrain BH and disrupted star properties, finding that a solar-mass main sequence star disrupted by a relatively low mass, M~10^5-10^6 Msun, BH is consistent with the data, while a WD disruption (though still possible) is disfavored. The magnetic flux required to power Swift J1644+57 is too large to be supplied by the star itself, but it could be collected from a quiescent `fossil' accretion disk present in the galactic nucleus prior to the TD. The presence (lack of) of such a fossil disk could be a deciding factor in what TD events are accompanied by powerful jets.[abridged]

preprint2013arXiv

The S2 star as a probe of the accretion disk of Sgr A*

How accretion proceeds around the massive black hole in the Galactic center and other highly sub-Eddington accretors remains poorly understood. The orbit of the S2 star in the Galactic center passes through the accretion disk of the massive black hole and any observational signature from such interaction may be used as an accretion probe. Because of its early stellar type, S2 is expected to possess a fairly powerful wind. We show here that the ram pressure of the accretion disk shocks the stellar wind fairly close to the star. The shocked fluid reaches a temperature of ~ 1 keV and cools efficiently through optically thin, thermal bremsstrahlung emission. The radiation from the shocked wind peaks around the epoch of the pericenter passage of the star at a luminosity potentially comparable to the quiescent emission detected from Sgr A*. Detection of shocked wind radiation can constrain the density of the accretion disk at a distance of several thousands of gravitational radii from the black hole.

preprint2012arXiv

Afterglow Model for the Radio Emission from the Jetted Tidal Disruption Candidate Swift J1644+57

The recent transient event Swift J1644+57 has been interpreted as emission from a collimated relativistic jet, powered by the sudden onset of accretion onto a supermassive black hole following the tidal disruption of a star. Here we model the radio-microwave emission as synchrotron radiation produced by the shock interaction between the jet and the gaseous circumnuclear medium (CNM). At early times after the onset of the jet (t < 5-10 days) a reverse shock propagates through and decelerates the ejecta, while at later times the outflow approaches the Blandford-McKee self-similar evolution (possibly modified by additional late energy injection). The achromatic break in the radio light curve of J1644+57 is naturally explained as the transition between these phases. We show that the temporal indices of the pre- and post-break light curve are consistent with those predicted if the CNM has a wind-type radial density profile n ~ 1/r^2. The observed synchrotron frequencies and self-absorbed flux constrain the fraction of the post-shock thermal energy in relativistic electrons epsilon_e ~ 0.03-0.1; the CNM density at 1e18 cm ~ 1-10 1/cm^3; and the initial Lorentz factor Gamma_j ~ 10-20 and opening angle theta_j ~ (0.1-1)Gamma_j^(-1) ~ 0.01-0.1 of the jet. Radio modeling thus provides robust independent evidence for a narrowly collimated outflow. Extending our model to the future evolution of J1644+57, we predict that the radio flux at low frequencies (<~ few GHz) will begin to brighten more rapidly once the characteristic frequency crosses below the self-absorption frequency on a timescale of months (indeed, such a transition may already have begun). Our results demonstrate that relativistic outflows from tidal disruption events provide a unique probe of the conditions in distant, previously inactive galactic nuclei, complementing studies of normal AGN.

preprint2012arXiv

Optical and X-ray Transients from Planet-Star Mergers

We evaluate the prompt observational signatures of the merger between a massive close-in planet (a `hot Jupiter') and its host star, events with an estimated Galactic rate of ~0.1-1/yr. Depending on the ratio of the mean density of the planet rho_p to that of the star rho_star, a merger results in three possible outcomes. If rho_p/rho_star > 5, then the planet directly plunges below the stellar atmosphere before being disrupted by tidal forces. Dissipation of orbital energy creates a hot wake behind the planet, producing a EUV/soft X-ray transient as the planet sinks below the stellar surface. The peak luminosity L_X ~ 1e36 erg/s is achieved weeks-months prior to merger, after which the stellar surface is enshrouded by an outflow. The final inspiral is accompanied by an optical transient powered by the recombination of hydrogen in the outflow, which peaks at L~1e37-38 erg/s on a timescale ~days. If instead rho_planet/rho_star < 5, then Roche Lobe overflow occurs above the stellar surface. For rho_p/rho_star < 1, mass transfer is stable, resulting the planet being accreted on a relatively slow timescale. However, for 1 < rho_p/rho_star < 5, mass transfer may instead be unstable, resulting in the planet being dynamically disrupted into an accretion disk around the star. Super-Eddington outflows from the disk power an optical transient with L~1e37-38 erg/s and characteristic duration ~week-months. The disk itself becomes visible once the accretion rate become sub-Eddington, resulting in a bolometric brightening and spectral shift to the UV. Optical transients from planet merger events may resemble classical novae, but are distinguished by lower ejecta mass and velocity ~100s km/s, and by hard pre- and post-cursor emission, respectively. Promising search strategies include combined optical, UV, and X-ray surveys of nearby massive galaxies with cadences from days to months.

preprint2012arXiv

The peak energy of dissipative GRB photospheres

The radiation released at the transparency radius of an ultrarelativistic flow can account for the observed properties of gamma-ray bursts (GRBs) provided that sufficient energy is dissipated in the sub-photospheric region. Here, I investigate how the peak energy of the E*f(E) spectrum and its overall shape depend on the properties of the jet for various "dissipative photospheres". I find that continuous energy release which results in electron heating over a wide range of distances may be the key to explain the GRB emission. In this picture, the peak of the spectrum forms at a Thomson optical depth of several tens. The peak depends mainly on the bulk Lorentz factor Gamma of the flow and can, therefore, be used to determine it. The Gamma is predicted to range from ~10 to 1000 from X-ray flashes to the brightest observed GRBs in agreement with recent observational inferences. The Amati relation can be understood if the brightest bursts are the least baryon loaded ones. Implications from this interpretation of the GRB emission for the central engine are discussed.

preprint2011arXiv

A relativistic jetted outburst from a massive black hole fed by a tidally disrupted star

While gas accretion onto some massive black holes (MBHs) at the centers of galaxies actively powers luminous emission, the vast majority of MBHs are considered dormant. Occasionally, a star passing too near a MBH is torn apart by gravitational forces, leading to a bright panchromatic tidal disruption flare (TDF). While the high-energy transient Swift J164449.3+573451 ("Sw 1644+57") initially displayed none of the theoretically anticipated (nor previously observed) TDF characteristics, we show that the observations (Levan et al. 2011) suggest a sudden accretion event onto a central MBH of mass ~10^6-10^7 solar masses. We find evidence for a mildly relativistic outflow, jet collimation, and a spectrum characterized by synchrotron and inverse Compton processes; this leads to a natural analogy of Sw 1644+57 with a smaller-scale blazar. The phenomenologically novel Sw 1644+57 thus connects the study of TDFs and active galaxies, opening a new vista on disk-jet interactions in BHs and magnetic field generation and transport in accretion systems.

preprint2011arXiv

GRB afterglow light curves from realistic density profiles

The afterglow emission that follows gamma-ray bursts (GRBs) contains valuable information about the circumburst medium and, therefore, about the GRB progenitor. Theoretical studies of GRB blast waves, however, are often limited to simple density profiles for the external medium (mostly constant density and power-law R^{-k} ones). We argue that a large fraction of long-duration GRBs should take place in massive stellar clusters where the circumburst medium is much more complicated. As a case study, we simulate the propagation of a GRB blast wave in a medium shaped by the collision of the winds of O and Wolf-Rayet stars, the typical distance of which is d /sim 0.1 - 1 pc. Assuming a spherical blast wave, the afterglow light curve shows a flattening followed by a shallow break on a timescale from hours up to a week after the burst, which is a result of the propagation of the blast wave through the shocked wind region. If the blast wave is collimated, the jet break may, in some cases, become very pronounced with the post break decline of the light curve as steep as t-5. Inverse Compton scattering of ultra-violet photons from the nearby star off energetic electrons in the blast wave leads to a bright \simGeV afterglow flare that may be detectable by Fermi.

preprint2011arXiv

Heavy Nuclei Synthesized in Gamma-Ray Burst Outflows as the Source of UHECRs

Recent measurements by the Pierre Auger Observatory suggest that the composition of ultra-high energy cosmic rays (UHECRs) becomes dominated by heavy nuclei at high energies. However, until now there has been no astrophysical motivation for considering a source highly enriched in heavy elements. Here we demonstrate that the outflows from Gamma-Ray Bursts (GRBs) may indeed be composed primarily of nuclei with masses A ~ 40-200, which are synthesized as hot material expands away from the central engine. In particular, if the jet is magnetically-dominated (rather than a thermally-driven fireball) its low entropy enables heavy elements to form efficiently. Adopting the millisecond proto-magnetar model for the GRB central engine, we show that heavy nuclei are both synthesized in proto-magnetar winds and can in principle be accelerated to energies >1e20 eV in the shocks or regions of magnetic reconnection that are responsible for powering the GRB. Similar results may apply to accretion-powered GRB models if the jet originates from a magnetized disk wind. Depending on the precise distribution of nuclei synthesized, we predict that the average primary mass may continue to increase beyond Fe group elements at the highest energies, possibly reaching the A ~ 90 (Zirconium), A ~ 130 (Tellurium), or even A ~ 195 (Platinum) peaks. Future measurements of the UHECR composition at energies >~ 1e20 eV can thus confirm or constrain our model and, potentially, probe the nature of GRB outflows. The longer attenuation length of ultra-heavy nuclei through the extragalactic background light greatly expands the volume of accesible sources and alleviates the energetic constraints on GRBs as the source of UHECRs.

preprint2011arXiv

Radio transients from stellar tidal disruption by massive black holes

The tidal disruption of a star by a supermassive black hole provides us with a rare glimpse of these otherwise dormant beasts. It has long been predicted that the disruption will be accompanied by a thermal `flare', powered by the accretion of bound stellar debris. Several candidate disruptions have been discovered in this manner at optical, UV and X-ray wavelengths. Here we explore the observational consequences if a modest fraction of the accretion power is channeled into an ultra-relativistic outflow. We show that a relativistic jet decelerates due to its interaction with the interstellar medium at sub-parsec distances from the black hole. Synchrotron radiation from electrons accelerated by the reverse shock powers a bright radio-infrared transient that peaks on a timescale ~1 yr after disruption. Emission from the forward shock may be detectable for several years after the peak. Deep radio follow-up observations of tidal disruption candidates at late times can test for the presence of relativistic ejecta. Upcoming radio transient surveys may independently discover tens to hundreds of tidal disruptions per year, complimenting searches at other wavelengths. Non-thermal emission from tidal disruption probes the physics of jet formation under relatively clean conditions, in which the flow parameters are independently constrained.

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

Superflares from magnetars revealing the GRB central engine

Long-duration gamma-ray bursts (GRBs) may be powered by the rotational energy of a millisecond magnetar. I argue that the GRB-driving magnetars lie at the high end of the distribution of magnetic field strengths of magnetars. The field of GRB magnetars decays on timescale of hundreds of years and can power SGR-like flares up to ~100 times more powerful than the 2004 event of SGR 1806--20. A few of these flares per year may have been observed by {\it BATSE} and classified as short-duration GRBs. Association of one of these superflares with a nearby d_L\simless 250 Mpc galaxy and the discovery of a, coincident in space, 100-year-old GRB afterglow (observed in the radio) will be the characteristic signature of the magnetar model for GRBs.

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.

preprint2009arXiv

Signatures of a Maxwellian Component in Shock-Accelerated Electrons in GRBs

Recent particle-in-cell simulations suggest that a large fraction of the energy dissipated in a relativistic shock is deposited into a Maxwellian distribution of electrons that is connected to the high-energy power-law tail. Here, we explore the observational implications of such a mixed thermal-nonthermal particle distribution for the afterglow and prompt emission of gamma-ray bursts. When the Maxwellian component dominates the energy budget, the afterglow lightcurves show a very steep decline phase followed by a more shallow decay when the characteristic synchrotron frequency crosses the observed band. The steep decay appears in the X-rays at ~100 sec after the burst and is accompanied by a characteristic hard-soft-hard spectral evolution that has been observed in a large number of early afterglows. If internal shocks produce a similar mixed electron distribution, a bump is expected at the synchrotron peak of the nu*f_nu spectrum.

preprint2008arXiv

Structure of Neutron Stars in Tensor-Vector-Scalar Theory

Bekenstein's Tensor-Vector-Scalar (TeVeS) theory has had considerable success in explaining various phenomena without the need for dark matter. However, it is difficult to observationally discern the differences between TeVeS and predictions made within the Lambda-cold dark matter concordance model. This implies that alternative tests are required that independently verify which theory is correct. For this we turn to the strong-field regime of TeVeS. In particular, we solve the spherically symmetric equations of hydrostatic equilibrium for a perfect fluid with a realistic equation of state to build models of neutron stars in TeVeS. We show that causality within the neutron star is only maintained for certain cosmological values of the scalar field, which allows us to put constraints on this value independently of cosmological observations. We also discuss in detail the internal structure of neutron stars and how each of the free parameters in the theory effects the overall size and mass of the neutron stars. In particular, the radii of neutron stars in TeVeS can significantly differ from those in General Relativity for certain values of the vector field coupling, which allows us to also place extra constraints on this parameter. Finally, we discuss future observations of neutron stars using both the electromagnetic and gravitational wave spectrums that will allow for tests of the appropriate theory of gravity.

preprint2005arXiv

Spherically symmetric, static spacetimes in a tensor-vector-scalar theory

Recently, a relativistic gravitation theory has been proposed [J. D. Bekenstein, Phys. Rev. D {\bf 70}, 083509 (2004)] that gives the Modified Newtonian Dynamics (or MOND) in the weak acceleration regime. The theory is based on three dynamic gravitational fields and succeeds in explaining a large part of extragalactic and gravitational lensing phenomenology without invoking dark matter. In this work we consider the strong gravity regime of TeVeS. We study spherically symmetric, static and vacuum spacetimes relevant for a non-rotating black hole or the exterior of a star. Two branches of solutions are identified: in the first the vector field is aligned with the time direction while in the second the vector field has a non-vanishing radial component. We show that in the first branch of solutions the βand γPPN coefficients in TeVeS are identical to these of general relativity (GR) while in the second the βPPN coefficient differs from unity violating observational determinations of it (for the choice of the free function $F$ of the theory made in Bekenstein's paper). For the first branch of solutions, we derive analytic expressions for the physical metric and discuss their implications. Applying these solutions to the case of black holes, it is shown that they violate causality (since they allow for superluminal propagation of metric, vector and scalar waves) in the vicinity of the event horizon and/or that they are characterized by negative energy density carried by the fields.