Source author record

Tsvi Piran

Tsvi Piran 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 gr-qc astro-ph astro-ph.CO astro-ph.GA astro-ph.SR hep-ph hep-th astro-ph.EP astro-ph.IM hep-ex nucl-th physics.plasm-ph quant-ph

Catalog footprint

What is connected

117works

14topics

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.

preprint2026arXiv

The Advanced X-ray Imaging Satellite (AXIS) Community Science Book

The AXIS Community Science Book represents the collective effort of 592 scientists worldwide to define the transformative science enabled by the Advanced X-ray Imaging Satellite (AXIS), a next-generation X-ray mission selected by NASA's Astrophysics Probe Program for Phase A study. AXIS will advance the legacy of high-angular-resolution X-ray astronomy with ~1.5'' imaging over a wide 24' field of view and an order of magnitude greater collecting area than Chandra in the 0.3-12 keV band. Combining sharp imaging, high throughput, and rapid response capabilities, AXIS will open new windows on virtually every aspect of modern astrophysics, exploring the birth and growth of supermassive black holes, the feedback processes that shape galaxies, the life cycles of stars and exoplanet environments, and the nature of compact stellar remnants, supernova remnants, and explosive transients. This book compiles 138 community-contributed science cases developed by five Science Working Groups focused on AGN and supermassive black holes, galaxy evolution and feedback, compact objects and supernova remnants, stellar physics and exoplanets, and time-domain and multi-messenger astrophysics. Together, these studies establish the scientific foundation for next-generation X-ray exploration in the 2030s and highlight strong synergies with facilities of the 2030s, such as JWST, Roman, Rubin/LSST, SKA, ALMA, ngVLA, and next-generation gravitational-wave and neutrino networks.

preprint2022arXiv

Analytic Modeling of Synchrotron-Self-Compton Spectra: Application to GRB 190114C

Observations of TeV emission from early gamma-ray burst (GRB) afterglows revealed the long sought for inverse Compton (IC) upscattering of the lower energy synchrotron. However, it turned out that the long hoped for ability to easily interpret the synchrotron-self-Compton (SSC) spectra didn't materialize. The TeV emission is in the Klein-Nishina (KN) regime and the simple Thomson regime SSC spectrum is modified, complicating the scene. We describe here a methodology to obtain an analytic approximation to an observed spectrum and infer the conditions at the emitting region. The methodology is general and can be used in any such source. As a test case, we apply it to the observations of GRB 190114C. We find that the procedure of fitting the model parameters using the analytic SSC spectrum suffers from some generic problems. However, at the same time, it conveniently gives a useful insight into the conditions that shape the spectrum. Once we introduce a correction to the standard KN approximation, the best fit solution is consistent with the one found in detailed numerical simulations. As in the numerical analysis, we find a family of solutions that provide a good approximation to the data and satisfy roughly $B\propto Γ^{-3}$ between the magnetic field and the bulk Lorentz factor, and we provide a tentative explanation why such a family arises.

preprint2022arXiv

What powers the radio emission in TDE AT2019dsg: a long-lived jet or the disruption itself?

The tidal disruption event AT2019dsg was observed from radio to X-rays and was possibly accompanied by a high-energy neutrino. Previous interpretations have focused on continued injection by a central engine as the source of energy for radio emission. We show that continuous energy injection is unnecessary; the radio data can be explained by a single ejection of plasma that supplies all the energy needed. To support this assertion, we analyze the synchrotron self-absorbed spectra in terms of the equipartition model. Similar to previous analyses, we find that the energy in the radio-emitting region increases approximately $\propto t^{0.7}$ and the lengthscale of this region grows $\propto t$ at a rate $\simeq0.06c$. This event resembles the earliest stage of a supernova remnant: because the ejected mass is much greater than the shocked external mass, its velocity remains unchanged, while the energy in shocked gas grows with time. The radio-emitting material gains energy from the outflow, not continuing energy injection by the central object. Although energy injection from an accreting BH cannot be completely excluded, the energy injection rate is very different from the fallback luminosity, and maintaining constant outflow velocity requires fine-tuning demanding further physical explanation. If the neutrino association is real, the energy injection needed is much greater than for the radio emission, suggesting that the detected neutrino did not arise from the radio-emitting region.

preprint2021arXiv

High-energy emission from tidal disruption events in active galactic nuclei

Tidal disruption events (TDEs) taking place in active galactic nuclei (AGNs) are different from ordinary TDEs. In these events, the returning tidal debris stream drills through the pre-existing AGN accretion disk near the stream pericenter, destroying the inner disk in the process, and then intersects with the disk a second time at radii ranging from a few to hundreds of times the pericenter distance. The debris dynamics of such TDEs, and hence their appearance, are distinct from ordinary TDEs. Here we explore the observational signatures of this "second impact" of the stream with the disk. Strong shocks form as the dilute stream is stopped by the denser disk. Compton cooling of the shocked material produces hard X-rays, even soft gamma-rays, with most of the energy emitted between ~10 keV and 1 MeV. The luminosity follows the mass-return rate, peaking between ~$10^{42}$ and $10^{44}$ erg/s. The X-ray hardness and the smoothness of the light curve provide possible means for distinguishing the second impact from ordinary AGN flares, which exhibit softer spectra and more irregular light curves.

preprint2021arXiv

Limits on mass outflow from optical tidal disruption events

The discovery of optical/UV tidal disruption events (TDEs) was surprising. The expectation was that, upon returning to the pericenter, the stellar-debris stream will form a compact disk that will emit soft X-rays. Indeed the first TDEs were discovered in this energy band. A common explanation for the optical/UV events is that surrounding optically-thick matter reprocesses the disk's X-ray emission and emits it from a large photosphere. If accretion follows the super-Eddington mass infall rate it would inevitably result in an energetic outflow, providing naturally the reprocessing matter. We describe here a new method to estimate, using the observed luminosity and temperature, the mass and energy of outflows from optical transients. When applying this method to a sample of supernovae our estimates are consistent with a more detailed hydrodynamic modeling. For the current sample of a few dozen optical TDEs the observed luminosity and temperature imply outflows that are significantly more massive than typical stellar masses, posing a problem to this common reprocessing picture.

preprint2020arXiv

Generalized compactness limit from an arbitrary viewing angle

A $γ$-ray source must have a limited optical depth to pair production. This simple condition, called compactness, implies that gamma-ray bursts (GRBs) must involve a highly relativistic motion ($Γ\gtrsim 100$) giving the first and most important clue on their nature. So far, this condition has been discussed under the assumption that the $γ$-ray sources are viewed on-axis, that is, by an observer within the beaming cone of the relativistic source. Recently, following the detection of the weak short GRB 170817A, an extensive interest arose in the possibility that some $γ$-ray sources are viewed off-axis. We generalize here the compactness formalism for an arbitrary viewing angle taking several possible opacity processes and $γ$-ray spectra into account. We find that for a given observables (peak luminosity, temporal variability, and spectra) the minimal Lorentz factor, $Γ_{\rm min}$, is obtained, as expected, for an on-axis observer. More remarkably we find that compactness dictates also a maximal viewing angle, $θ_{\rm max} \simeq 1/2Γ_{\rm min}$. Our limit implies for regular GRBs a very small allowed viewing angle ($\lesssim10^{-2}\,\rm rad$), making it extremely unlikely that they are viewed off-axis. For GRB 170817A we confirm earlier results that rule out the possibility that the observed $γ$-rays were seen by an on-axis observer as a regular short GRB. The short GRB 150101B was also suggested to be an off-axis event. We show that its maximal viewing angle $\lesssim0.05\,\rm rad$, which is inconsistent with the off-axis model. Finally we show that for low luminosity GRBs, compactness does not exclude by itself an off-axis model, but when combined with other consideration this option is strongly disfavored.

preprint2020arXiv

On short GRBs similar to GRB 170817A detected by Fermi GBM

von Kienlin et al. (2019) selected 11 short gamma-ray bursts (sGRBs) whose characteristics are similar to GRB 170817A. These bursts, like GRB 170817A, have a hard spike followed by a soft thermal tail. However, as their redshifts are unknown it is not clear if their luminosities are as low as that of GRB 170817A. Comparing the positions in the $ε_{\rm p}$-$E_{\rmγ,iso}$ (spectral peak energy - isotropic-equivalent energy) plane and using compactness arguments to estimate the minimal Lorentz factor, $Γ$, we find that all the bursts in this sample are consistent with being regular sGRBs if they are located at $z\simeq0.3-3$. They are also consistent with being similar to GRB 170817A if they are located at $z\lesssim0.1$. Even in the latter case, the events must involve at least mildly relativistic ($Γ\gtrsim 2$) motion within the sources. We, further, find that at most one or two bursts in the sample are consistent with the cocoon shock-breakout model. Finally, we calculate the event rate of either off-axis emission from a jet core or from a jet-wing (surrounding the core). We find that the off-axis emission model as an origin of the sample is rejected as it predicts too small event rate. The wing model can be consistent with the observed rate but the model parameters cannot be constrained by the current observations.

preprint2020arXiv

Radiative Emission Mechanisms of Tidal Disruption Events

We describe how the various outcomes of stellar tidal disruption give rise to observable radiation. We separately consider the cases where gas circularizes rapidly into an accretion disc, as well as the case when shocked debris streams provide the observable emission without having fully circularized. For the rapid circularization case, we describe how outflows, absorption by reprocessing layers, and Comptonization can cause the observed radiation to depart from that of a bare disc, possibly giving rise to the observed optical/UV emission along with soft X-rays from the disc. If, instead, most of the debris follows highly eccentric orbits for a significant time, many properties of the observed optical/UV emission can be explained by the scale of those eccentric orbits and the shocks embedded in the debris flow near orbital apocenter. In this picture, soft X-ray emission at early times results from the smaller amount of debris mass deflected into a compact accretion disc by weak shocks near the stellar pericenter. A general proposal for the near-constancy of the ultraviolet/optical color temperatures is provided, by linking it to incomplete thermalization of radiation in the atmosphere of the emitting region. We also briefly discuss the radio signals from the interaction of unbound debris and jets with the black hole environment.

preprint2020arXiv

Shock within a shock: revisiting the radio flares of NS merger ejecta and GRB-supernovae

Fast ejecta expelled in binary neutron star (NS) mergers or energetic supernovae (SNe) should produce late-time synchrotron radio emission as the ejecta shocks into the surrounding ambient medium. Models for such radio flares typically assume the ejecta expands into an unperturbed interstellar medium (ISM). However, it is also well-known that binary NS mergers and broad-lined Ic SNe can harbor relativistic jetted outflows. In this work, we show that such jets shock the ambient ISM ahead of the ejecta, thus evacuating the medium into which the ejecta subsequently collides. Using an idealized spherically-symmetric model, we illustrate that this inhibits the ejecta radio flare at early times $t < t_{\rm col} \approx 12 \, {\rm yr} \, (E_{\rm j}/10^{49} \, {\rm erg})^{1/3} (n/1 \, {\rm cm}^{-3})^{-1/3} (v_{\rm ej}/0.1c)^{-5/3}$ where $E_{\rm j}$ is the jet energy, $n$ the ISM density, and $v_{\rm ej}$ the ejecta velocity. We also show that this can produce a sharply peaked enhancement in the light-curve at $t = t_{\rm col}$. This has implications for radio observations of GW170817 and future binary NS mergers, gamma-ray burst (GRB) SNe, decade-long radio transients such as FIRST J1419, and possibly other events where a relativistic outflow precedes a slower-moving ejecta. Future numerical work will extend these analytic estimates and treat the multi-dimensional nature of the problem.

preprint2020arXiv

The Origin of Binary Black Holes Mergers

Recently Venumadhav et al. (2019) proposed a new pipeline to analyze LIGO-Virgo's O1-O2 data and discovered eight new binary black hole (BBH) mergers, including a high effective spin, chi_{eff}, one. This discovery sheds new light on the origin of the observed BBHs and the dynamical capture vs. field binaries debate. Using a new statistical approach, we show that, while isotropic models are not ruled out, the observed chi_{eff} distribution favors field binaries, whose chi_{eff} is determined by tidal forces and wind losses, over capture. However, given the current limited sample size, capture scenarios (isotropic models) cannot be ruled out. Observations of roughly a hundred merges will enable us to distinguish between the different formation scenarios. However, if as expected, both formation channels operate it may be difficult to resolve their exact fraction.

preprint2020arXiv

The origin of hotspots around Sgr A*: Orbital or pattern motion?

The Gravity Collaboration detected a near-infrared hotspot moving around Sgr A* during the 2018 July 22 flare. They fitted the partial loop the hotspot made on the sky with a circular Keplerian orbit of radius $\simeq7.5\,r_{\rm g}$ around the supermassive black hole (BH), where $r_{\rm g}$ is the gravitational radius. However, because the hotspot traversed the loop in a short time, models in which the hotspot tracks the motion of some fluid element tend to produce a best-fit trajectory smaller than the observed loop. This is true for a circular Keplerian orbit, even when BH spin is accounted for, and for motion along a RIAF streamline. A marginally bound geodesic suffers from the same problem; in addition, it is not clear what the origin of an object following the geodesic would be. The observed hotspot motion is more likely a pattern motion. Circular motion with $r\simeq12.5\,r_{\rm g}$ and a super-Keplerian speed $\simeq0.8\,c$ is a good fit. Such motion must be pattern motion because it cannot be explained by physical forces. The pattern speed is compatible with magnetohydrodynamic perturbations, provided that the magnetic field is sufficiently strong. Circular pattern motion of radius $\sim20\, r_{\rm g}$ on a plane above the BH is an equally good alternative; in this case, the hotspot may be caused by a precessing outflow interacting with a surrounding disk. As all our fits have relatively large radii, we cannot constrain the BH spin using these observations.

preprint2020arXiv

Tidal disruption events in active galactic nuclei

A fraction of tidal disruption events (TDEs) occur in active galactic nuclei (AGNs) whose black holes possess accretion disks; these TDEs can be confused with common AGN flares. The disruption itself is unaffected by the disk, but the evolution of the bound debris stream is modified by its collision with the disk when it returns to pericenter. The outcome of the collision is largely determined by the ratio of the stream mass current to the azimuthal mass current of the disk rotating underneath the stream footprint, which in turn depends on the mass and luminosity of the AGN. To characterize TDEs in AGNs, we simulated a suite of stream--disk collisions with various mass current ratios. The collision excites shocks in the disk, leading to inflow and energy dissipation orders of magnitude above Eddington; however, much of the radiation is trapped in the inflow and advected into the black hole, so the actual bolometric luminosity may be closer to Eddington. The emergent spectrum may not be thermal, TDE-like, or AGN-like. The rapid inflow causes the disk interior to the impact point to be depleted within a fraction of the mass return time. If the stream is heavy enough to penetrate the disk, part of the outgoing material eventually hits the disk again, dissipating its kinetic energy in the second collision; another part becomes unbound, emitting synchrotron radiation as it shocks with surrounding gas.

preprint2020arXiv

Tidal radii of main sequence stars -- I. Physical tidal radius, semi-analytic model and their implications

A star is tidally disrupted by a supermassive black hole when their separation is shorter than the "tidal radius". This quantity is often estimated on an order-of-magnitude basis without reference to the star's internal structure. Using MESA models for main sequence stars and fully general relativistic dynamics, we find the physical tidal radius for complete disruption $\cal{R}_t$ for a $10^6M_\odot$ black hole (BH). We find that across a factor $\sim20$ in stellar mass $M_*$, i.e., $0.15M_{\odot}\leq M_*\leq3M_\odot$, $\cal{R}_t\sim27\times$(BH's gravitational radius). When comparing $\cal{R}_t$ with the commonly used order-of-magnitude estimate $r_t$, we find that $\cal{R}_t\sim1.05-1.45r_t$ for $0.15M_\odot\leq M_*\leq0.5M_\odot$, but between $0.5 M_\odot$ and $1 M_\odot$, $\cal{R}_t$ drops to $\sim 0.45r_t$, and it remains at this value up to $10 M_\odot$. The near-constancy of $\cal{R}_t$ implies a weaker dependence of the full disruption rate on $M_*$ than when predicted with $r_t$. The characteristic energy width of the debris $ΔE$ ranges from $\sim1.2Δ\cal{E}$ for low-mass stars to $\sim 0.35Δ\cal{E}$ for higher-mass stars, where $Δ\cal{E}=GM_{\rm BH}R_*/\cal{R}_t^{2}$. We present analytic fits for the $M_*$ dependence of $\cal{R}_t$ and $ΔE$; these fits lead to analytic expressions for the time of peak mass fallback rate and the maximal mass fallback rate. Our results also bear on the fraction of events leading to fast or slow circularization, as well as on the character of the tidal event occurring when the remnant of a partial disruption returns to the black hole. Using a semi-analytic model, we show that $\cal{R}_t$ is primarily determined by the star's central density rather than its mean density. For high-mass stars, the full disruption rate is roughly 1/4 the partial disruption rate, while this ratio is close to unity for low-mass stars.

preprint2020arXiv

Tidal radii of main sequence stars -- II. Simulation methodology and the character of full tidal disruptions

This is the second in a series of papers presenting the results of fully general relativistic simulations of stellar tidal disruptions in which the stars' initial states are realistic main-sequence models. We consider eight different stellar masses, from $0.15~{\rm M}_\odot$ to $10~{\rm M}_\odot$. In the first paper (Ryu et al. 2019a), we gave an overview of this program and discussed the principal observational implications of our work. Here we describe our calculational method and provide details about the outcomes of full disruptions. We find that, relative to the traditional order-of-magnitude estimate $r_{\rm t}$, the physical tidal radius of low-mass stars is larger by tens of percent, while for high-mass stars ($M_{\star} \gtrsim1~ {\rm M}_\odot$) it is smaller by a factor $2-2.5$. The traditional estimate of the range of energies found in the debris is approximately accurate for low-mass stars, but is a factor $\sim 2$ too small for high-mass stars; in addition, the energy distribution for high-mass stars has significant wings. For all stars undergoing tidal encounters, we find that mass-loss continues for a long time because the ${\it instantaneous}$ tidal radius, the distance out to which the black hole's tidal gravity competes with the instantaneous stellar gravity at the star's surface, stays comparable to the distance to the black hole until the star has reached $O(10)~r_{\rm t}$. These findings indicate significant failings in the popular "frozen-in" approximation.

preprint2020arXiv

Tidal radii of main sequence stars -- III. Partial disruptions

In this paper, the third in this series, we continue our study of tidal disruption events of main-sequence stars by a non-spinning $10^{6}~\rm{M}_\odot$ supermassive black hole. Here we focus on the outcomes of partial disruptions. As the encounter becomes weaker, the debris mass is increasingly concentrated near the outer edges of the energy distribution. As a result, the mass fallback rate can deviate substantially from a $t^{-5/3}$ power-law, becoming more like a single peak with a tail declining as $t^{-p}$ with $p\simeq2-5$. Surviving remnants are spun-up in the prograde direction and are hotter than MS stars of the same mass. Their specific orbital energy is $\simeq10^{-3}\times$ that of the debris (but of either sign with respect to the black hole potential) while their specific angular momentum is close to that of the original star. Even for strong encounters, remnants have speeds at infinity relative to the black hole potential $\lesssim 300$ km s$^{-1}$, so they are unable to travel far out into the galactic bulge. Remnants bound to the black hole can possibly go through a second tidal disruption event.

preprint2018arXiv

Magnetorotational instability in eccentric disks

Eccentric disks arise in such astrophysical contexts as tidal disruption events, but it is unknown whether the magnetorotational instability (MRI), which powers accretion in circular disks, operates in eccentric disks as well. We examine the linear evolution of unstratified, incompressible MRI in an eccentric disk orbiting a point mass. We consider vertical modes of wavenumber $k$ on a background flow with uniform eccentricity $e$ and vertical Alfvén speed $v_\mathrm A$ along an orbit with mean motion $n$. We find two mode families, one with dominant magnetic components, the other with dominant velocity components; the former is unstable at $(1-e)^3f^2\lesssim3$, where $f\equiv kv_\mathrm A/n$, the latter at $e\gtrsim0.8$. For $f^2\lesssim3$, MRI behaves much like in circular disks, but the growth per orbit declines slowly with increasing $e$; for $f^2\gtrsim3$, modes grow by parametric amplification, which is resonant for $0<e\ll1$. MRI growth and the attendant angular momentum and energy transport happen chiefly near pericenter, where orbital shear dominates magnetic tension.

preprint2016arXiv

A Revised Analysis of Gamma Ray Bursts' prompt efficiencies

The prompt Gamma-Ray Bursts' (GRBs) efficiency is an important clue on the emission mechanism producing the $γ$-rays. Previous estimates of the kinetic energy of the blast waves, based on the X-ray afterglow luminosity $L_X$, suggested that this efficiency is large, with values above 90\% in some cases. This poses a problem to emission mechanisms and in particular to the internal shocks model. These estimates are based, however, on the assumption that the X-ray emitting electrons are fast cooling and that their Inverse Compton (IC) losses are negligible. The observed correlations between $L_X$ (and hence the blast wave energy) and $E_{γ\rm ,iso}$, the isotropic equivalent energy in the prompt emission, has been considered as observational evidence supporting this analysis. It is reasonable that the prompt gamma-ray energy and the blast wave kinetic energy are correlated and the observed correlation corroborates, therefore, the notion $L_X$ is indeed a valid proxy for the latter. Recent findings suggest that the magnetic field in the afterglow shocks is significantly weaker than was earlier thought and its equipartition fraction, $ε_B$, could be as low as $10^{-4}$ or even lower. Motivated by these findings we reconsider the problem, taking now IC cooling into account. We find that the observed $L_X-E_{γ\rm ,iso}$ correlation is recovered also when IC losses are significant. For small $ε_B$ values the blast wave must be more energetic and we find that the corresponding prompt efficiency is significantly smaller than previously thought. For example, for $ε_B\sim10^{-4}$ we infer a typical prompt efficiency of $\sim15\%$.

preprint2016arXiv

Are Ultra-Long Gamma-Ray Bursts Caused by Blue Supergiant Collapsars, Newborn Magnetars, or White Dwarf Tidal Disruption Events?

Ultra-long gamma-ray bursts (ulGRBs) are a new population of GRBs with extreme durations of $\sim 10^{4}$ s. Leading candidates for their origin are blue supergiant collapsars, magnetars, and white dwarf tidal disruption events (WD-TDEs) caused by massive black holes (BHs). Recent observations of supernova-like (SN-like) bumps associated with ulGRBs challenged both the WD-TDE and the blue supergiant models because of the detection of SNe and the absence of hydrogen lines, respectively. We propose that WD-TDEs can accommodate the observed SN-like bumps if the fallback WD matter releases energy into the unbound WD ejecta. The observed ejecta energy, luminosity, and velocity are explained by the gravitational energy, Eddington luminosity, and escape velocity of the formed accretion disk, respectively. We also show that the observed X-rays can ionize the ejecta, eliminating lines. The SN-like light curves (SN 2011kl) for the ulGRB 111209A are consistent with all three models, although a magnetar model is unnatural because the spin-down time required to power the SN-like bump is a hundred times longer than the GRB. Our results imply that TDEs are a possible energy source for SN-like events in general and for ulGRBs in particular.

preprint2016arXiv

ASASSN-14li: A Model Tidal Disruption Event

ASASSN-14li is a recently-discovered tidal disruption event with an exceptionally rich data-set: spectra and lightcurves in soft X-rays, UV, optical, and radio. To understand its emission properties in all these bands, we have extended our model for post-tidal disruption accretion and photon production to estimate both soft X-ray radiation produced by the prompt accretion phase and synchrotron emission associated with the bow shock driven through an external medium by the unbound tidal debris, as well as optical and UV light. We find that fiducial values of the stellar mass ($1 M_\odot$) and black hole mass ($10^{6.5} M_{\odot}$) yield: quantitative agreement with the optical/UV luminosity, lightcurve, and color temperature; approximate agreement with the somewhat uncertain soft X-ray spectrum and lightcurve; and quantitative agreement with the radio luminosity, spectrum and lightcurve. Equipartition analysis of the radio data implies that the radio-emitting region expands with a constant speed, and its magnitude is comparable to the speed expected for the unbound stellar ejecta. Both facts provide strong support to our model. We find that the disruption event took place in mid-September 2014. Two independent parameters, the magnitude and logarithmic radial gradient of the ambient gas density near the black hole, must be fit to the data to explain the radio emission; their inferred values are comparable to those found near both Sgr A* and the TDE candidate Swift J1644.

preprint2016arXiv

Beaming of particles and synchrotron radiation in relativistic magnetic reconnection

Relativistic reconnection has been invoked as a mechanism for particle acceleration in numerous astrophysical systems. According to idealised analytical models reconnection produces a bulk relativistic outflow emerging from the reconnection sites (X-points). The resulting radiation is therefore highly beamed. Using two-dimensional particle-in-cell (PIC) simulations, we investigate particle and radiation beaming, finding a very different picture. Instead of having a relativistic average bulk motion with isotropic electron velocity distribution in its rest frame, we find that the bulk motion of particles in X-points is similar to their Lorentz factor gamma, and the particles are beamed within about 5/gamma. On the way from the X-point to the magnetic islands, particles turn in the magnetic field, forming a fan confined to the current sheet. Once they reach the islands they isotropise after completing a full Larmor gyration and their radiation is not strongly beamed anymore. The radiation pattern at a given frequency depends on where the corresponding emitting electrons radiate their energy. Lower energy particles that cool slowly spend most of their time in the islands, and their radiation is not highly beamed. Only particles that quickly cool at the edge of the X-points generate a highly beamed fan-like radiation pattern. The radiation emerging from these fast cooling particles is above the burn-off limit (about 100 MeV in the overall rest frame of the reconnecting plasma.) This has significant implications for models of GRBs and AGNs that invoke beaming in that frame at much lower energies.

preprint2016arXiv

Cosmic explosions, life in the Universe and the Cosmological Constant

Galactic Gamma-Ray Bursts (GRBs) are copious sources of gamma-rays that can pose a threat to complex life. Using recent determinations of their rate and the probability of GRBs causing massive extinction, we explore what type of universes are most likely to harbour advanced forms of life. For this purpose we use cosmological N-body simulations to determine at what time and for what value of the cosmological constant ($Λ$) the chances of life being unaffected by cosmic explosions are maximised. We find that $Λ-$dominated universes favour the survival of life against GRBs. Within a $Λ$CDM cosmology, the parameters that govern the likelihood of life survival to GRBs are dictated by the value of $Λ$ and the age of the Universe. We find that we seem to live in a favorable point in this parameter phase space which minimises the exposure to cosmic explosions, yet maximises the number of main sequence (hydrogen-burning) stars around which advanced life forms can exist.

preprint2016arXiv

Formation of Double Neutron Star systems as implied by observations

Double Neutron Stars (DNS) have to survive two supernovae and still remain bound. This sets strong limits on the nature of the second collapse in these systems. We consider the masses and orbital parameters of the DNS population and constrain the two distributions of mass ejection and kick velocities directly from observations with no a-priori assumptions regarding evolutionary models and/or the types of the supernovae involved. We show that there is strong evidence for two distinct types of supernovae in these systems, where the second collapse in the majority of the observed systems involved small mass ejection ($ΔM\lesssim 0.5M_{\odot}$) and a corresponding low-kick velocity ($v_{k}\lesssim 30 $km\,s$^{-1}$). This formation scenario is compatible, for example, with an electron capture supernova. Only a minority of the systems have formed via the standard SN scenario involving larger mass ejection of $\sim 2.2 M_{\odot}$ and kick velocities of up to $400$km\,s$^{-1}$. Due to the typically small kicks in most DNS (which are reflected by rather low proper motion), we predict that most of these systems reside close to the galactic disc. In particular, this implies that more NS-NS mergers occur close to the galactic plane. This may have non-trivial implications to the estimated merger rates of DNS and to the rate of LIGO / VIRGO detections.

preprint2016arXiv

Natal Kicks and Time Delays in Merging Neutron Star Binaries - Implications for r-process nucleosynthesis in Ultra Faint Dwarfs and in the Milky Way

Merging neutron star binaries are prime candidate sources for heavy r-process nucleosynthesis. The amount of heavy r-process material is consistent with the mass ejection and rates of mergers, and abundances of relic radioactive materials suggest that heavy r-process material is produced in rare events. Observations of possible macronovae provide further support for this model. Still, some concerns remain. One is the observation of heavy r-process elements in Ultra Faint Dwarf (UFD) galaxies. The escape velocities from UFDs are so small that the natal kicks, taking place at neutron stars birth, might eject such binaries from UFDs. Furthermore the old stellar populations of UFDs requires that r-process nucleosynthesis must have taken place very early on, while it may take several Gyr for compact binaries to merge. This last problem arises also within the Milky Way where heavy r-process materials has been observed in some low metallicity stars. We show here that since a significant fraction of neutron star binaries form with a very small proper motion those won't be ejected even from a UFD. Furthermore, a significant fraction of these binaries merge within a few hundred Myr and in some cases the merger time is even shorter than a hundred Myr. This population of "rapid mergers" explains the appearance of heavy r-process material in both UFDs and in the Milky Way.

preprint2016arXiv

On the maximal efficiency of the collisional Penrose process

The center of mass (CM) energy in a collisional Penrose process - a collision taking place within the ergosphere of a Kerr black hole - can diverge under suitable extreme conditions (maximal Kerr, near horizon collision and suitable impact parameters). We present an analytic expression for the CM energy, refining expressions given in the literature. Even though the CM energy diverges, we show that the maximal energy attained by a particle that escapes the black hole's gravitational pull and reaches infinity is modest. We obtain an analytic expression for the energy of an escaping particle resulting from a collisional Penrose process, and apply it to derive the maximal energy and the maximal efficiency for several physical scenarios: pair annihilation, Compton scattering, and the elastic scattering of two massive particles. In all physically reasonable cases (in which the incident particles initially fall from infinity towards the black hole) the maximal energy (and the corresponding efficiency) are only one order of magnitude larger than the rest mass energy of the incident particles. The maximal efficiency found is ~ 13.92 and it is obtained for the scattering of an outgoing massless particle by a massive particle.

preprint2016arXiv

Particle acceleration, magnetization and radiation in relativistic shocks

What are the mechanisms of particle acceleration and radiation, as well as magnetic field build up and decay in relativistic shocks are open questions with important implications to various phenomena in high energy astrophysics. While the Weibel instability is possibly responsible for magnetic field build up and diffusive shock acceleration is a model for acceleration, both have problems and current PIC simulation show that particles are accelerated only under special conditions and the magnetic field decays on a short length scale. We present here a novel model for the structure and the emission of highly relativistic collisionless shocks. The model takes into account (and is based on) non-local energy and momentum transport across the shock front via emission and absorption of high-energy photons. This leads to a pre-acceleration of the fluid and pre-amplificaiton of the magnetic fields in the upstream region. Both have drastic implications on the shock structure. The model explains the persistence of the shock generated magnetic field at large distances from the shock front. The dissipation of this magnetic field results in a continuous particle acceleration within the downstream region. The model suggests two non-uniform emission zones (the downstream and the upstream), that give rise to three emission components with different spectral and temporal properties. A unique feature of the model is the existence of an "attractor", toward which any shock will evolve. This enables us to estimate from first principles the synchrotron and inverse Compton spectrum of the downstream emission. The model is applicable to any relativistic shock, but its distinctive features show up only for large compactness. We demonstrate that prompt and afterglow Gamma-Ray Bursts' shocks satisfy the relevant conditions and compare their observations with the predictions.

preprint2016arXiv

r-process Production Sites as inferred from Eu Abundances in Dwarf Galaxies

Recent observations of $r$-process material in ultra-faint dwarf galaxies (UFDs) shed light on the sources of these elements. Strong upper limits on the Eu mass in some UFDs combined with detections of much larger masses in a UFD, Reticulum II, and other dwarf galaxies imply that Eu production is dominated by rare events, and that the minimal Eu mass observed in any UFD is approximately the amount of Eu mass produced per event. This is consistent with other independent observations in the Galaxy. We estimate, using a model independent likelihood analysis, the rate and Eu (Fe) mass produced per $r$-process (Fe production) event in dwarf galaxies including classical dwarfs and UFDs. The mass and rate of the Fe production events are consistent with the normal core-collapse supernova~(ccSN) scenario. The Eu mass per event is $3\times 10^{-5}M_{\odot}<\tilde{m}_{\rm Eu}<2\times 10^{-4}M_{\odot}$, corresponding to a total $r$-process mass per event of $6\times 10^{-3}M_{\odot}<\tilde{m}_{r-process}<4\times 10^{-2}M_{\odot}$. The rate of $r$-process events is $2.5\times 10^{-4}<R_{rp/SN}<1.4\times 10^{-3}$ as compared with the ccSNe rate. These values are consistent with the total Eu mass observed in our own Galaxy, suggesting that the same mechanism is behind the production of $r$-process events in both dwarf galaxies and the Milky Way, and that it may be the dominant mechanism for production of $r$-process elements in the Universe. The results are consistent with neutron star mergers estimates, but cannot rule out other rare core collapse scenarios, provided that they produce a significant amounts of $r$-process material per event.

preprint2016arXiv

Radio Counterparts of Compact Binary Mergers detectable in Gravitational Waves: A Simulation for an Optimized Survey

Mergers of binary neutron stars and black hole-neutron star binaries produce gravitational-wave (GW) emission and outflows with significant kinetic energies. These outflows result in radio emissions through synchrotron radiation. We explore the detectability of these synchrotron generated radio signals by follow-up observations of GW merger events lacking a detection of electromagnetic counterparts in other wavelengths. We model radio light curves arising from (i) sub-relativistic merger ejecta and (ii) ultra-relativistic jets. The former produces radio remnants on timescales of a few years and the latter produces $γ$-ray bursts in the direction of the jet and orphan-radio afterglows extending over wider angles on timescales of weeks. Based on the derived light curves, we suggest an optimized survey at $1.4$ GHz with five epochs separated by a logarithmic time interval. We estimate the detectability of the radio counterparts of simulated GW-merger events to be detected by advanced LIGO and Virgo by current and future radio facilities. The detectable distances for these GW merger events could be as high as 1 Gpc. $20$--$60\%$ of the long-lasting radio remnants will be detectable in the case of the moderate kinetic energy of $3\cdot 10^{50}$ erg and a circum-merger density of $0.1 {\rm cm^{-3}}$ or larger, while $5$--$20\%$ of the orphan radio afterglows with kinetic energy of $10^{48}$ erg will be detectable. The detection likelihood increases if one focuses on the well-localizable GW events. We discuss the background noise due to radio fluxes of host galaxies and false positives arising from extragalactic radio transients and variable Active Galactic Nuclei and we show that the quiet radio transient sky is of great advantage when searching for the radio counterparts.

preprint2016arXiv

Testing the Magnetar Model via Late Time Radio Observations of Two Macronova Candidates

Compact binary mergers may have already been observed as they are the leading model for short gamma-ray bursts (sGRBs). Radioactive decay within the ejecta from these mergers is expected to produce an infra-red flare, dubbed macronova (or kilonova), on a time scale of a week. Recently two such macronova candidates were identified in followup observations of sGRBs, strengthening the possibility that those indeed arise from mergers. The same ejecta will also produce a long term (months to years) radio emission due to its interaction with the surrounding ISM. In search for this emission, we observed the two macronova candidates, GRB 130603B and GRB 060614 with the Jansky very large array (VLA) and the Australia Telescope Compact Array (ATCA). Our observations resulted in null-detections, putting strong upper limits on the kinetic energy and mass of the ejecta. A possible outcome of a merger is a highly magnetized neutron star (a magnetar), which has been suggested as the central engine for GRBs. Such a magnetar will deposit a significant fraction of its energy into the ejecta leading to a brighter radio flare. Our results, therefore, rule out magnetars in these two events.

preprint2016arXiv

The B/C and sub-Iron/Iron Cosmic ray ratios - further evidence in favor of the spiral arm diffusion model

The Boron to Carbon (B/C) and sub-Fe/Fe ratios provides an important clue on Cosmic Ray (CR) propagation within the Galaxy. These ratios estimate the grammage that the CR traverse as they propagate from their sources to Earth. Attempts to explain these ratios within the standard CR propagation models require ad hoc modifications and even with those these models necessitate inconsistent grammages to explain both ratios. As an alternative, physically motivated model, we have proposed that CR originate preferably within the galactic spiral arms. CR propagation from dynamic spiral arms has important imprints on various secondary to primary ratios, such as the B/C ratio and the positron fraction. We use our spiral arm diffusion model with the spallation network extended up to Nickel to calculate the sub-Fe/Fe ratio. We show that without any additional parameters the spiral arm model consistently explains both ratios with the same grammage, providing further evidence in favor of this model.

preprint2016arXiv

The dynamics of radiation driven, optically thick winds

Recent observation of some luminous transient sources with low color temperatures suggests that the emission is dominated by optically thick winds driven by super-Eddington accretion. We present a general analytical theory of the dynamics of radiation pressure-driven, optically thick winds. Unlike the classical adiabatic stellar wind solution whose dynamics are solely determined by the sonic radius, here the loss of the radiation pressure due to photon diffusion also plays an important role. We identify two high mass loss rate regimes ($\dot{M} > L_{\rm Edd\,}/c^2$). In the large total luminosity regime the solution resembles an adiabatic wind solution. Both the radiative luminosity, $L$, and the kinetic luminosity, $L_k$, are super-Eddington with $L < L_k$ and $L \propto L_k^{1/3}$. In the lower total luminosity regime most of the energy is carried out by the radiation with $L_k < L \approx L_{\rm Edd\,}$. In a third, low mass loss regime ($\dot{M} < L_{\rm Edd\,}/c^2$), the wind becomes optically thin early on and, unless gas pressure is important at this stage, the solution is very different from the adiabatic one. The results are independent from the energy generation mechanism at the foot of the wind, therefore they are applicable to a wide range of mass ejection systems, from black hole accretion, to planetary nebulae, and to classical novae.

preprint2016arXiv

The effect of cooling on particle trajectories and acceleration in relativistic magnetic reconnection

The maximum synchrotron burnoff limit of 160 MeV represents a fundamental limit to radiation resulting from electromagnetic particle acceleration in one-zone ideal plasmas. In magnetic reconnection, however, particle acceleration and radiation are decoupled because the electric field is larger than the magnetic field in the diffusion region. We carry out two-dimensional particle-in-cell simulations to determine the extent to which magnetic reconnection can produce synchrotron radiation above the burnoff limit. We use the test particle comparison (TPC) method to isolate the effects of cooling by comparing the trajectories and acceleration efficiencies of test particles incident on such a reconnection region with and without cooling them. We find that the cooled and uncooled particle trajectories are typically similar during acceleration in the reconnection region, and derive an effective limit on particle acceleration that is inversely proportional to the average magnetic field experienced by the particle during acceleration. Using the calculated distribution of this average magnetic field as a function of uncooled final particle energy, we find analytically that cooling does not affect power-law particle energy spectra except at energies far above the synchrotron burnoff limit. Finally, we compare fully cooled and uncooled simulations of reconnection, confirming that the synchrotron burnoff limit does not produce a cutoff in the particle energy spectrum. Our results indicate that the TPC method accurately predicts the effects of cooling on particle acceleration in relativistic reconnection, and that even far above the burnoff limit, the synchrotron energy of radiation produced in reconnection is not limited by cooling.

preprint2016arXiv

The Macronova in GRB 050709 and the GRB/macronova connection

GRB 050709 was the first short Gamma-ray Burst (sGRB) with an identified optical afterglow. In this work we report a re-analysis of the publicly available data of this event and the discovery of a Li-Paczynski macronova/kilonova that dominates the optical/IR signal at t $>$ 2.5 days. Such a signal would arise from ~0.05 M$_\odot$ r-process material launched by a compact binary merger. The implied mass ejection supports the suggestion that compact binary mergers are significant and possibly main sites of heavy r-process nucleosynthesis. We have re-analyzed all afterglow data from nearby short and hybrid GRBs. A statistical study of sGRB/macronova connection reveals that macronova may have taken place in all these GRBs though the fraction as low as 0.18 cannot be ruled out. The identification of two of the three macronova candidates in the I-band implies a more promising detection prospect for ground-based surveys.

preprint2016arXiv

The Observable Signatures of GRB Cocoons

As a long GRB jet propagates within the surrounding stellar atmosphere it creates a cocoon composed of an outer Newtonian shocked stellar material and an inner (possibly relativistic) shocked jet material. The jet deposits $10^{51}-10^{52}$ erg into this cocoon. This energy is comparable to the GRB's energy and to the energy of the accompanying supernova, yet its signature has been largely neglected so far. A fraction of the cocoon energy is released during its expansion following the breakout from the star and later as it interacts with the surrounding matter. We explore here the possible signatures of the cocoon emission and outline a framework to calculate them from the conditions of the cocoon at the time of the jet breakout. We show that the cocoon signature depends strongly on the level of mixing between the shocked jet and shocked stellar material that fills it, which is currently unknown. We find that if there is no mixing at all then the $γ$-ray emission from the cocoon is so bright that it should have been already detected, and the lack of such detections indicates that mixing at some level must take place. We calculate also the expected signal for partial and full mixing. While the typical signals are weaker than GRBs' afterglows, the latter are highly beamed while the former have wide angles. We predict that future optical, UV and X-ray transient searches, like LSST, ZTF, ULTRASAT, ISS-Lobster and others will most likely detect such signals, providing a wealth of information on the progenitors and jets of GRBs. While we focus on long GRBs, we note that analogous (but weaker) cocoons may arise in short GRBs as well. Their signatures might be the most promising electromagnetic counterparts for gravitational waves merger's signals.

preprint2015arXiv

"Circularization" vs. Accretion -- What Powers Tidal Disruption Events?

A tidal disruption event (TDE) takes place when a star passes near enough to a massive black hole to be disrupted. About half the star's matter is given elliptical trajectories with large apocenter distances, the other half is unbound. To "circularize", i.e., to form an accretion flow, the bound matter must lose a significant amount of energy, with the actual amount depending on the characteristic scale of the flow measured in units of the black hole's gravitational radius ($\sim 10^{51} (R/1000R_g)^{-1}$~erg). Recent numerical simulations \citep{Shiokawa+2015} have revealed that the circularization scale is close to the scale of the most-bound initial orbits, $\sim 10^3 M_{BH,6.5}^{-2/3} R_g \sim 10^{15} M_{BH,6.5}^{1/3}$~cm from the black hole, and the corresponding circularization energy dissipation rate is $\sim 10^{44} M_{BH,6.5}^{-1/6}$~erg/s. We suggest that the energy liberated during circularization, rather then energy liberated by accretion onto the black hole, powers the observed optical TDE candidates. The observed rise times, luminosities, temperatures, emission radii, and line widths seen in these TDEs \citep[e.g.][]{Arcavi+2014} are all more readily explained in terms of heating associated with circularization than in terms of accretion.

preprint2015arXiv

A possible Macronova in the late afterglow of the `long-short' burst GRB 060614

Long-duration ($>2$ s) $γ$-ray bursts that are believed to originate from the death of massive stars are expected to be accompanied by supernovae. GRB 060614, that lasted 102 s, lacks a supernova-like emission down to very stringent limits and its physical origin is still debated. Here we report the discovery of near-infrared bump that is significantly above the regular decaying afterglow. This red bump is inconsistent with even the weakest known supernova. However, it can arise from a Li-Paczyński macronova $-$ the radioactive decay of debris following a compact binary merger. If this interpretation is correct GRB 060614 arose from a compact binary merger rather than from the death of a massive star and it was a site of a significant production of heavy r-process elements. The significant ejected mass favors a black hole-neutron star merger but a double neutron star merger cannot be ruled out.

preprint2015arXiv

Energies of GRB blast waves and prompt efficiencies as implied by modeling of X-ray and GeV afterglows

We consider a sample of ten GRBs with long lasting ($\gtrsim10^2\rm\,sec$) emission detected by Fermi/LAT and for which X-ray data around $1\,$day are also available. We assume that both the X-rays and the GeV emission are produced by electrons accelerated at the external forward shock, and show that the X-ray and the GeV fluxes lead to very different estimates of the initial kinetic energy of the blast wave. The energy estimated from GeV is on average $\sim50$ times larger than the one estimated from X-rays. We model the data (accounting also for optical detections around $1\,$day, if available) to unveil the reason for this discrepancy and find that good modelling within the forward shock model is always possible and leads to two possibilities: either the X-ray emitting electrons (unlike the GeV emitting electrons) are in the slow cooling regime or ii) the X-ray synchrotron flux is strongly suppressed by Compton cooling, whereas, due to the Klein-Nishina suppression, this effect is much smaller at GeV energies. In both cases the X-ray flux is no longer a robust proxy for the blast wave kinetic energy. On average, both cases require weak magnetic fields ($10^{-6}\lesssim ε_B \lesssim 10^{-3}$) and relatively large isotropic kinetic blast wave energies $10^{53}\rm\,erg<E_{0,kin}<10^{55}\rm\,erg$ corresponding to large lower limits on the collimated energies, in the range $10^{52}\rm\,erg<E_{θ,kin}<5\times10^{52}\rm\,erg$ for an ISM environment with $n\sim 1\mbox{cm}^{-3}$ and $10^{52}\rm\,erg<E_{θ,kin}<10^{53}\rm\,erg$ for a wind environment with $A_* \sim 1$. These energies are larger than those estimated from the X-ray flux alone, and imply smaller inferred values of the prompt efficiency mechanism, reducing the efficiency requirements on the still uncertain mechanism responsible for prompt emission.

preprint2015arXiv

Gamma-Ray Bursts as Sources of Strong Magnetic Fields

Gamma-Ray Bursts (GRBs) are the strongest explosions in the Universe, which due to their extreme character likely involve some of the strongest magnetic fields in nature. This review discusses the possible roles of magnetic fields in GRBs, from their central engines, through the launching, acceleration and collimation of their ultra-relativistic jets, to the dissipation and particle acceleration that power their $γ$-ray emission, and the powerful blast wave they drive into the surrounding medium that generates their long-lived afterglow emission. An emphasis is put on particular areas in which there have been interesting developments in recent years.

preprint2015arXiv

General Relativistic Hydrodynamic Simulation of Accretion Flow from a Stellar Tidal Disruption

We study how the matter dispersed when a supermassive black hole tidally disrupts a star joins an accretion flow. Combining a relativistic hydrodynamic simulation of the stellar disruption with a relativistic hydrodynamics simulation of the tidal debris motion, we track such a system until ~80% of the stellar mass bound to the black hole has settled into an accretion flow. Shocks near the stellar pericenter and also near the apocenter of the most tightly-bound debris dissipate orbital energy, but only enough to make the characteristic radius comparable to the semi-major axis of the most-bound material, not the tidal radius as previously thought. The outer shocks are caused by post-Newtonian effects, both on the stellar orbit during its disruption and on the tidal forces. Accumulation of mass into the accretion flow is non-monotonic and slow, requiring ~3--10x the orbital period of the most tightly-bound tidal streams, while the inflow time for most of the mass may be comparable to or longer than the mass accumulation time. Deflection by shocks does, however, remove enough angular momentum and energy from some mass for it to move inward even before most of the mass is accumulated into the accretion flow. Although the accretion rate rises sharply and then decays roughly as a power-law, its maximum is ~0.1x the previous expectation, and the duration of the peak is ~5x longer than previously predicted. The geometric mean of the black hole mass and stellar mass inferred from a measured event timescale is therefore ~0.2x the value given by classical theory.

preprint2015arXiv

Jet and disk luminosities in tidal disruption events

Tidal disruption events (TDE) in which a star is devoured by a massive black hole at a galac- tic center pose a challenge to our understanding of accretion processes. Within a month the accretion rate reaches super-Eddington levels. It then drops gradually over a time scale of a year to sub-Eddington regimes. The initially geometrically thick disk becomes a thin one and eventually an ADAF at very low accretion rates. As such, TDEs explore the whole range of accretion rates and configurations. A challenging question is what the corresponding light curves of these events are. We explore numerically the disk luminosity and the conditions within the inner region of the disk using a fully general relativistic slim disk model. Those conditions determine the magnitude of the magnetic field that engulfs the black hole and this, in turn, determines the Blandford-Znajek jet power. We estimate this power in two different ways and show that they are self-consistent. We find, as expected earlier from analytic argu- ments (Krolik & Piran 2012), that neither the disk luminosity nor the jet power follows the accretion rate throughout the disruption event. The disk luminosity varies only logarithmi- cally with the accretion rate at super-Eddington luminosities. The jet power follows initially the accretion rate but remains a constant after the transition from super- to sub- Eddington. At lower accretion rates at the end of the MAD phase the disk becomes thin and the jet may stop altogether. These new estimates of the jet power and disk luminosity that do not simply follow the mass fallback rate should be taken into account when searching for TDEs and analysing light curves of TDE candidates. Identification of some of the above mentioned transitions may enable us to estimate better TDE parameters.

preprint2015arXiv

Linear and circular polarization in GRB afterglows

A certain degree of linear polarization has been measured in several GRB afterglows. More surprisingly, circular polarization has been recently measured in GRB121024A. For synchrotron emission, the polarization level depends on: (i) the local magnetic field orientation (ii) the geometry of the emitting region with respect to the line of sight and (iii) the electron pitch-angle distribution. For this reason, polarization measurements are a valuable tool to probe afterglow micro-physics. We present numerical estimates of linear and circular polarization for different configurations (i.e., magnetic fields, geometries and pitch-angle distributions). For each different scenario, we study the conditions for reaching the maximum and minimum linear and circular polarization and provide their values. We discuss the implication of our results to the micro-physics of GRB afterglows in view of recent polarization measurements.

preprint2015arXiv

Linear and Circular polarization in ultra-relativistic synchrotron sources - implications to GRB afterglows

Polarization measurements from relativistic outflows are a valuable tool to probe the geometry of the emission region and the microphysics of the particle distribution. Indeed, the polarization level depends on: (i) the local magnetic field orientation, (ii) the geometry of the emitting region with respect to the line of sight, and (iii) the electron pitch-angle distribution. Here we consider optically thin synchrotron emission and we extend the theory of circular polarization from a point source to an extended radially expanding relativistic jet. We present numerical estimates for both linear and circular polarization in such systems. We consider different configurations of the magnetic field, spherical and jetted outflows, isotropic and anisotropic pitch-angle distributions, and outline the difficulty in obtaining the reported high level of circular polarization observed in the afterglow of GRB 121024A. We conclude that the origin of the observed polarization cannot be intrinsic to an optically thin synchrotron process, even when the electron pitch-angle distribution is extremely anisotropic.

preprint2015arXiv

Mass ejection from neutron star mergers: different components and expected radio signals

In addition to producing a strong gravitational signal, a short gamma-ray burst (GRB), and a compact remnant, neutron star mergers eject significant masses at significant kinetic energies. This mass ejection takes place via dynamical mass ejection and a GRB jet but other processes have also been suggested: a shock-breakout material, a cocoon resulting from the interaction of the jet with other ejecta, and viscous and neutrino driven winds from the central remnant or the accretion disk. The different components of the ejected masses include up to a few percent of a solar mass, some of which is ejected at relativistic velocities. The interaction of these ejecta with the surrounding interstellar medium will produce a long lasting radio flare, in a similar way to GRB afterglows or to radio supernovae. The relative strength of the different signals depends strongly on the viewing angle. An observer along the jet axis or close to it will detect a strong signal at a few dozen days from the radio afterglow (or the orphan radio afterglow) produced by the highly relativistic GRB jet. For a generic observer at larger viewing angles, the dynamical ejecta, whose contribution peaks a year or so after the event, will generally dominate. Depending on the observed frequency and the external density, other components may also give rise to a significant contribution. We also compare these estimates with the radio signature of the short GRB 130603B. The radio flare from the dynamical ejecta might be detectable with the EVLA and the LOFAR for the higher range of external densities $n\gtrsim 0.5{\rm cm^{-3}}$.

preprint2015arXiv

Radio Flares of Compact Binary Mergers: the Effect of Non-Trivial Outflow Geometry

The next generation gravitational waves (GW) detectors are most sensitive to GW emitted by compact (neutron star/black hole) binary mergers. If one of those is a neutron star the merger will also emit electromagnetic radiation via three possible channels: Gamma-ray bursts and their (possibly orphan) afterglows (Eichler et al. 1989), Li-Paczynski Macronovae (Li & Paczynski 1998) and radio flares (Nakar & Piran 2011). This accompanying electromagnetic radiation is vitally important in confirming the GW detections (Kochanek & Piran 1993). It could also reveal a wealth of information regarding the merger and will open a window towards multi-messenger astronomy. Identifying and characterizing these counterparts is therefore of utmost importance. In this work we explore late time radio flares emitted by the dynamically ejected outflows. We build upon previous work and consider the effect of the outflow's non-trivial geometry. Using an approximate method we estimate the radio light-curves for several ejected matter distributions obtained in numerical simulations. Our method provides an upper limit to the effect of non-sphericity. Together with the spherical estimates the resulting light curves bound the actual signal. We find that while non-spherical geometries can in principle lead to an enhanced emission, in most cases they result in an increase in the timescale compared with a corresponding spherical configuration. This would weaken somewhat these signals and might decrease the detection prospects.

preprint2015arXiv

Radioactive decay products in neutron star merger ejecta: heating efficiency and $γ$-ray emission

The radioactive decay of the freshly synthesized $r$-process nuclei ejected in compact binary mergers power optical/infrared macronovae (kilonovae) that follow these events. The light curves depend critically on the energy partition among the different products of the radioactive decay and this plays an important role in estimates of the amount of ejected $r$-process elements from a given observed signal. We study the energy partition and $γ$-ray emission of the radioactive decay. We show that $20$-$50\%$ of the total radioactive energy is released in $γ$-rays on timescales from hours to a month. The number of emitted $γ$-rays per unit energy interval has roughly a flat spectrum between a few dozen keV and $1$ MeV so that most of this energy is carried by $\sim 1$ MeV $γ$-rays. However at the peak of macronova emission the optical depth of the $γ$-rays is $\sim 0.02$ and most of the $γ$-rays escape. The loss of these $γ$-rays reduces the heat deposition into the ejecta and hence reduces the expected macronova signals if those are lanthanides dominated. This implies that the ejected mass is larger by a factor of $2$-$3$ than what was previously estimated. Spontaneous fission heats up the ejecta and the heating rate can increase if a sufficient amount of transuranic nuclei are synthesized. Direct measurements of these escaping $γ$-rays may provide the ultimate proof for the macronova mechanisms and an identification of the $r$-process nucleosynthesis sites. However, the chances to detect these signals are slim with current X-ray and $γ$-ray missions. New detectors, more sensitive by at least a factor of ten, are needed for a realistic detection rate.

preprint2015arXiv

Shedding light on the prompt high efficiency paradox - self consistent modeling of GRB afterglows

We examine GRBs with both Fermi-LAT and X-ray afterglow data. Assuming that the 100MeV (LAT) emission is radiation from cooled electrons accelerated by external shocks, we show that the kinetic energy of the blast wave estimated from the 100MeV flux is 50 times larger than the one estimated from the X-ray flux. This can be explained if either: i) electrons radiating at X-rays are significantly cooled by SSC (suppressing the synchrotron flux above the cooling frequency) or ii) if the X-ray emitting electrons, unlike those emitting at 100MeV energies, are in the slow cooling regime. In both cases the X-ray flux is no longer an immediate proxy of the blast wave kinetic energy. We model the LAT, X-ray and optical data and show that in general these possibilities are consistent with the data, and explain the apparent disagreement between X-ray and LAT observations. All possible solutions require weak magnetic fields: $10^{-6}< ε_B < 10^{-3}$ (where $ε_B$ is the fraction of shocked plasma energy in magnetic fields). Using the LAT emission as a proxy for the blast wave kinetic energy we find that the derived prompt efficiencies are of order 15%. This is considerably lower compared with previous estimates (87% and higher for the same bursts). This provides at least a partial solution to the "prompt high efficiency paradox".

preprint2015arXiv

Spectral signatures of compact sources in the inverse Compton catastrophe limit

The inverse Compton catastrophe is defined as a dramatic rise in the luminosity of inverse Compton scattered photons. It is described by a non-linear loop of radiative processes that sets in for high values of the electron compactness and is responsible for the efficient transfer of energy from electrons to photons, predominantly through inverse Compton scatterings. We search for the conditions that drive a magnetized non-thermal source to the inverse Compton catastrophe regime and study its multi-wavelength (MW) photon spectrum. We develop a generic analytical framework and use numerical calculations as a backup to the analytical predictions. We find that the escaping radiation from a source in the Compton catastrophe regime bears some unique features. The MW photon spectrum is a broken power law with a break at $\sim m_e c^2$ due to the onset of the Klein-Nishina suppression. The spectral index below the break energy depends on the electron and magnetic compactnesses logarithmically, while it is independent of the electron power-law index ($s$). The maximum radiating power emerges typically in the $γ$-ray regime, at energies $\sim m_e c^2$ ($\sim γ_{\max} m_e c^2$ ) for $s>2$ ($s\lesssim 2$), where $γ_{\max}$ is the maximum Lorentz factor of the injected electron distribution. We apply the principles of the inverse Compton catastrophe to blazars and $γ$-ray bursts using the analytical framework we developed, and show how these can be used to impose robust constraints on the source parameters.

preprint2015arXiv

Super-Penrose collisions are inefficient - a Comment on: Black hole fireworks: ultra-high-energy debris from super-Penrose collisions

In a paper posted on the arXiv a few weeks ago Berti, Brito and Cardoso \cite{Berti+14} suggest that ultra-high-energy particles can emerge from collisions in a black hole's ergosphere. This can happen if the process involves a particle on an outgoing trajectory very close to the black hole. Clearly such a particle cannot emerge from the black hole. It is argued \cite{Berti+14} that this particle can arise in another collision. Thus the process involves two collisions: one in which an outgoing particle is produced extremely close to the horizon, and a second one in which energy is gained. The real efficiency of this process should take into account, therefore, the energy needed to produce the first particle. We show here that while this process is kinematically possible, it requires a deposition of energy that is divergently large compared with the energy of the escaping particle. Thus, in contradiction to claims of infinitely high efficiencies, the efficiency of the combined process is in fact extremely small, approaching zero for very high output energies. Even under more general conditions than those considered in \cite{Berti+14} the total energy gain never diverges, and is larger only by a factor of a few than the energy gain of the original collisional Penrose process that takes place between two infalling particles \cite{Piran+75,PiranShaham77,Bejger+12}

preprint2015arXiv

The afterglow of a relativistic shock breakout and low-luminosity GRBs

The prompt emission of low-luminosity gamma-ray bursts (llGRBs) indicates that these events originate from a relativistic shock breakout. In this case we can estimate, based on the properties of the prompt emission, the energy distribution of the ejecta. We develop a general formalism to estimate the afterglow produced by synchrotron emission from the forward shock resulting from the interaction of this ejecta with the circum-burst matter. We assess whether this emission can produce the observed radio and X-ray afterglows of the available sample of 4 llGRBs. All 4 radio afterglows can be explained within this model, providing further support for shock breakouts being the origin of llGRBs. We find that in one of the llGRBs (GRB 031203) the predicted X-ray emission, using the same parameters that fit the radio, can explain the observed one. In another one (GRB 980425) the observed X-rays can be explained if we allow for a slight modification of the simplest model. For the last two cases (GRBs 060218 and 100316D), we find that, as is the case for previous attempts to model these afterglows, the simplest model that fits the radio emission underpredicts the observed X-ray afterglows. Using general arguments, we show that the most natural location of the X-ray source is, like the radio source, within the ejecta-external medium interaction layer but that emission is due to a different population of electrons or to a different emission process.

preprint2015arXiv

The fate of supernova remnants near quiescent supermassive black holes

There is mounting observational evidence that most galactic nuclei host both supermassive black holes (SMBHs) and young populations of stars. With an abundance of massive stars, core-collapse supernovae are expected in SMBH spheres of influence. We develop a novel numerical method, based on the Kompaneets approximation, to trace supernova remnant (SNR) evolution in these hostile environments, where radial gas gradients and SMBH tides are present. We trace the adiabatic evolution of the SNR shock until 50% of the remnant is either in the radiative phase or is slowed down below the SMBH Keplerian velocity and is sheared apart. In this way, we obtain shapes and lifetimes of SNRs as a function of the explosion distance from the SMBH, the gas density profile and the SMBH mass. As an application, we focus here exclusively on quiescent SMBHs, because their light may not hamper detections of SNRs and because we can take advantage of the unsurpassed detailed observations of our Galactic Centre. Assuming that properties such as gas and stellar content scale appropriately with the SMBH mass, we study SNR evolution around other quiescent SMBHs. We find that, for SMBH masses over ~10^7 Msun, tidal disruption of SNRs can occur at less than 10^4 yr, leading to a shortened X-ray emitting adiabatic phase, and to no radiative phase. On the other hand, only modest disruption is expected in our Galactic Centre for SNRs in their X-ray stage. This is in accordance with estimates of the lifetime of the Sgr A East SNR, which leads us to expect one supernova per 10^4 yr in the sphere of influence of Sgr A*.

preprint2015arXiv

The rate, luminosity function and time delay of non-Collapsar short GRBs

We estimate the rate and the luminosity function of short (hard) Gamma-Ray Bursts (sGRBs) that are non-Collapsars, using the peak fluxes and redshifts of BATSE, Swift and Fermi GRBs. Following Bromberg2013 we select a sub-sample of Swift bursts which are most likely non-Collapsars. We find that these sGRBs are delayed relative to the global star formation rate (SFR) with a typical delay time of a 3-4 Gyr (depending on the SFR model). However, if two or three sGRB at high redshifts have been missed because of selection effects, a distribution of delay times of ~1/t would be also compatible. The current event rate of these non-Collapsar sGRBs with L_iso > 5*10^49 erg/s is 4.1(-1.9,+2.3)Gpc^-3 yr^-1. The rate was significantly larger around z ~ 1 and it declines since that time. The luminosity function we find is a broken power law with a break at 2.0(-0.4,+1.4) * 10^52~erg/s and power-law indices 0.95(-0.12,+0.12) and 2.0(-0.8,+1.0). When considering the whole Swift sGRB sample we find that it is composed of two populations: One group (~ 60%-80% of Swift sGRBs) with the above rate and time delay and a second group (~ 20%-40% of Swift sGRBs) of potential "impostors" that follow the SFR with no delay. These two populations are in very good agreement with the division of sGRBs to non-Collapsars and Collapsars suggested recently by Bromberg2013. If non-Collapsar sGRBs arise from neutron star merger this rate suggest a detection rate of 3-100 yr^-1 by a future gravitational wave detectors (e.g. Advanced Ligo/Virgo with detection horizon on 300 Mpc), and a co-detection with Fermi (Swift) rate of 0.1-1 yr^-1 (0.02-0.14 yr^-1). We estimate that about 4 * 10^5 (f_b^-1 / 30) mergers took place in the Milky Way. If $0.025 m_\odot$ were ejected in each event this would have been sufficient to produce all the heavy r-process material in the Galaxy.

preprint2014arXiv

Implications of GRB 130603B and its macronova for r-process nucleosynthesis

The tentative identification of a Li-Paczynski macronova following the short GRB 130603B indicated that a few hundredths of a solar mass of neutron star matter were ejected and that this ejected mass has radioactively decayed into heavy r-process elements. If correct, this confirms long standing predictions (Eichler, Livio, Piran and Schramm 1989) that on the one hand, sGRBs are produced in compact binary mergers (CBMs) and on the other hand that these events are significant and possibly dominant sources of the heavy (A>130) r-process nuclei. Assuming that this interpretation is correct we obtain a lower limit of 0.02 m_sun on the ejected mass. Using the current estimates of the rate of sGRBs and with a beaming factor of 50, mergers associated with sGRBs can produce all the observed heavy r-process material in the Universe. We confront this conclusion with cosmochemistry and show that even though such events are rare, mixing is sufficient to account for the current homogeneous distribution of r-process material in the Galaxy. However, the appearance of significant amounts of Eu in some very low metallicity stars requires that some mergers took place very early on, namely with a very short time delay after the earliest star formation episodes. Alternatively, an additional early r-processsource may have contributed at that early stage. Finally, we note that evidence for short lived {244}^Pu in the very early solar system suggests that a merger of this kind took place within the vicinity of the solar system shortly (a few hundred million years) before its formation.

preprint2014arXiv

On the Composition of GRBs' Collapsar Jets

The duration distribution of long Gamma Ray Bursts reveals a plateau at durations shorter than ~20 s (in the observer frame) and a power-law decline at longer durations (Bromberg et al., 2012). Such a plateau arises naturally in the Collapsar model. In this model the engine has to operate long enough to push the jet out of the stellar envelope and the observed duration of the burst is the difference between the engine's operation time and the jet breakout time. We compare the jet breakout time inferred from the duration distribution (~10 s in the burst's frame) to the breakout time of a hydrodynamic jet (~10 s for typical parameters) and of a Poynting flux dominated jet with the same overall energy (<~1 s). As only the former is compatible with the duration of the plateau in the GRB duration distribution, we conclude that the jet is hydrodynamic during most of the time that its head is within the envelope of the progenitor star and around the time when it emerges from the star. This would naturally arise if the jet forms as a hydrodynamic jet in the first place or if it forms Poynting flux dominated but dissipates most of its magnetic energy early on within the progenitor star and emerges as a hydrodynamic jet.

preprint2014arXiv

On the role of GRBs on life extinction in the Universe

As a copious source of gamma-rays, a nearby Galactic Gamma-Ray Burst (GRB) can be a threat to life. Using recent determinations of the rate of GRBs, their luminosity function and properties of their host galaxies, we estimate the probability that a life-threatening (lethal) GRB would take place. Amongst the different kinds of GRBs, long ones are most dangerous. There is a very good chance (but no certainty) that at least one lethal GRB took place during the past 5 Gyr close enough to Earth as to significantly damage life. There is a 50% chance that such a lethal GRB took place during the last 500 Myr causing one of the major mass extinction events. Assuming that a similar level of radiation would be lethal to life on other exoplanets hosting life, we explore the potential effects of GRBs to life elsewhere in the Galaxy and the Universe. We find that the probability of a lethal GRB is much larger in the inner Milky Way (95% within a radius of 4 kpc from the galactic center), making it inhospitable to life. Only at the outskirts of the Milky Way, at more than 10 kpc from the galactic center, this probability drops below 50%. When considering the Universe as a whole, the safest environments for life (similar to the one on Earth) are the lowest density regions in the outskirts of large galaxies and life can exist in only ~ 10% of galaxies. Remarkably, a cosmological constant is essential for such systems to exist. Furthermore, because of both the higher GRB rate and galaxies being smaller, life as it exists on Earth could not take place at $z > 0.5$. Early life forms must have been much more resilient to radiation.

preprint2014arXiv

The dynamics of a highly magnetized jet propagating inside a star

The collapsar model explains the association of long duration gamma-Ray Bursts (GRBs) with stellar collapse. It involves a relativistic jet that forms at the core of a collapsing massive star. The jet penetrates the stellar envelope and the prompt GRB emission is produced once the jet is well outside the star. Most current models for generation of relativistic jets involve Poynting flux dominated outflows. We explore here the propagation of such a jet through a stellar envelope. The jet forms a bow shock around it. Energy dissipation at the head of this shock supplies energy to a cocoon that surrounds the jet. This cocoon exerts pressure on the jet and collimates it. While this description resembles the propagation of a hydrodynamic jets there are significant qualitative differences. Two Strong shocks, the reverse shock that slows down the hydrodynamic jet and the collimation shock that collimates it, cannot form within the Poynting flux dominated jet. As a result this jet moves much faster and dissipates much less energy while it crosses the stellar envelope. We construct here a simple analytic model that explores, self consistently, the jet-cocoon interaction and dynamics. Using this model we determine the properties of the jet, including its velocity, propagation time and shape.

preprint2014arXiv

The emission mechanism in magnetically dominated GRB outflows

We consider the conditions within a Gamma-Ray Burst (GRB) emission region that is Poynting flux dominated. Due to the enormous magnetic energy density, relativistic electrons will cool in such a region extremely rapidly via synchrotron. As there is no known mechanism that can compete with synchrotron it must be the source of the prompt sub-MeV emission. This sets strong limits on the size and Lorentz factor of the outflow. Furthermore, synchrotron cooling is too efficient. It overproduces optical and X-ray as compared with the observations. This overproduction of low energy emission can be avoided if the electrons are re-accelerated many times ($\gtrsim 5\times 10^4$)during each pulse (or are continuously heated) or if they escape the emitting region before cooling down. We explore the limitations of both models practically ruling out the later and demonstrating that the former requires two different acceleration mechanisms as well as an extremely large magnetic energy to Baryonic energy ratio. To be viable, any GRB model based on an emission region that is Poynting flux dominated must demonstrate how these conditions are met. We conclude that if GRB jets are launched magnetically dominated they must dissipate somehow most of their magnetic energy before they reach the emission region.

preprint2014arXiv

The long-term evolution of neutron star merger remnants - II. Radioactively powered transients

We use 3D hydrodynamic simulations of the long-term evolution of neutron star merger ejecta to predict the light curves of electromagnetic transients that are powered by the decay of freshly produced r-process nuclei. For the dynamic ejecta that are launched by tidal and hydrodynamic interaction, we adopt grey opacities of 10 cm$^2$/g, as suggested by recent studies. For our reference case of a 1.3-1.4 $M_\odot$ merger, we find a broad IR peak 2-4 d after the merger. The peak luminosity is $\approx 2\times 10^{40}$ erg/s for an average orientation, but increased by up to a factor of 4 for more favourable binary parameters and viewing angles. These signals are rather weak and hardly detectable within the large error box (~100 deg$^2$) of a gravitational wave trigger. A second electromagnetic transient results from neutrino-driven winds. These winds produce `weak' r-process material with $50 < A < 130$ and abundance patterns that vary substantially between different merger cases. For an adopted opacity of 1 cm$^2$/g, the resulting transients peak in the UV/optical about 6 h after the merger with a luminosity of $\approx 10^{41}$ erg/s (for a wind of 0.01 $M_\odot$) These signals are marginally detectable in deep follow-up searches (e.g. using Hypersuprime camera on Subaru). A subsequent detection of the weaker but longer lasting IR signal would allow an identification of the merger event. We briefly discuss the implications of our results to the recent detection of an nIR transient accompanying GRB 130603B.

preprint2014arXiv

The Nature of ULX Source M101 X-1: Optically Thick Outflow from A Stellar Mass Black Hole

The nature of ultra-luminous X-ray sources (ULXs) has long been plagued by an ambiguity about whether the central compact objects are intermediate-mass (IMBH, >~ 10^3 M_sun) or stellar-mass (a few tens M_sun) black holes (BHs). The high luminosity (~ 10^39 erg/s) and super-soft spectrum (T ~ 0.1 keV) during the high state of the ULX source X-1 in the galaxy M101 suggest a large emission radius (>~ 10^9 cm), consistent with being an IMBH accreting at a sub-Eddington rate. However, recent kinematic measurement of the binary orbit of this source and identification of the secondary as a Wolf-Rayet star suggest a stellar-mass BH primary with a super-Eddington accretion. If that is the case, a hot, optically thick outflow from the BH can account for the large emission radius and the soft spectrum. By considering the interplay of photons' absorption and scattering opacities, we determine the radius and mass density of the emission region of the outflow and constrain the outflow mass loss rate. The analysis presented here can be potentially applied to other ULXs with thermally dominated spectra, and to other super-Eddington accreting sources.

preprint2014arXiv

Tidal disruption jets as the source of Ultra-High Energy Cosmic Rays

Observations of the spectacular, blazar-like tidal disruption event (TDE) candidates Swift J1644+57 and J2058+05 show that the conditions required for accelerating protons to 10^{20} eV appear to be realized in the outer jet, and possibly in the inner jet as well. Direct and indirect estimates of the rate of jetted-TDEs, and of the energy they inject, are compatible with the observed flux of ultra-high energy cosmic rays (UHECRs) and the abundance of presently contributing sources. Thus TDE-jets can be a major source of UHECRs, even compabile with a pure proton composition.

preprint2013arXiv

Binary Pulsar J0737-3039: Evidence for a new core collapse and neutron star formation mechanism

The binary pulsar J0737-3039 is the only known system having two observable pulsars, thus offering a unique laboratory to test general relativity and explore pulsar physics. Based on the low eccentricity and the position within the galactic plane, Piran & Shaviv (2004, 2005) argued that pulsar B had a non-standard formation scenario with little or no mass ejection. They have also predicted that the system would have a very slow proper motion. Pulsar timing measurements (Kramer et al. 2006; Deller et al. 2009) confirmed this prediction. The recent observations of the alignment between the spin of pulsar A and the binary orbit is also in agreement with this scenario. Detailed simulations of the formation process of pulsar B enable us to show that its progenitor, just before the collapse, was a massive O-Ne-Mg white dwarf surrounded by a tenuous, 0.1-0.16 M_sun, envelope. This envelope was ejected when the white dwarf collapsed to form a neutron star. Pulsar B was born as a slow rotator (spin period ~ 1 s) and a kick received when the pulsar formed changed its spin direction to the current one. This realization sheds light on the angular momentum evolution of the progenitor star, a process which is strongly affected by interaction with the binary companion. The slow proper motion of the system also implies that the system must have undergone a phase of mass transfer in which Star A shed a significant fraction of its mass onto B.

preprint2013arXiv

Constraints on the Synchrotron Emission Mechanism in GRBs

We reexamine the general synchrotron model for GRBs' prompt emission and determine the regime in the parameter phase space in which it is viable. We characterize a typical GRB pulse in terms of its peak energy, peak flux and duration and use the latest Fermi observations to constrain the high energy part of the spectrum. We solve for the intrinsic parameters at the emission region and find the possible parameter phase space for synchrotron emission. Our approach is general and it does not depend on a specific energy dissipation mechanism. Reasonable synchrotron solutions are found with energy ratios of $10^{-4}<ε_B/ε_e<10$, bulk Lorentz factor values of $300<Γ<3000$, typical electrons' Lorentz factor values of $3\times 10^3<γ_e<10^5$ and emission radii of the order $10^{15}cm<R<10^{17}$cm. Most remarkable among those are the rather large values of the emission radius and the electron's Lorentz factor. We find that soft (with peak energy less than 100KeV) but luminous (isotropic luminosity of $~1.5 \times 10^{53}$) pulses are inefficient. This may explain the lack of strong soft bursts. In cases when most of the energy is carried out by the kinetic energy of the flow, such as in the internal shocks, the synchrotron solution requires that only a small fraction of the electrons are accelerated to relativistic velocities by the shocks. We show that future observations of very high energy photons from GRBs by CTA, could possibly determine all parameters of the synchrotron model or rule it out altogether.

preprint2013arXiv

Gravitational Wave Memory from Gamma Ray Bursts' Jets

While the possible roles of GRBs' progenitors as Gravitational Waves (GW) sources have been extensively studied, little attention has been given to the jet itself as a GW source. We expect the jet's acceleration to produce a GW Memory signal. While all relativistic jet models display anti-beaming of GW radiation away from the jet axis, thus radiating away from directions of GRBs' gamma radiation, this effect is not overwhelming. The decrease of the signal amplitude towards the cone of gamma-ray detectability is weak, and for some models and parameters the GW signal reaches a significant amplitude for much of the gamma-ray cone. Thus both signals may be jointly detected. We find different waveforms and fourier signatures for uniform jets and structured jet models - thus offering a method of using GW signatures to probe the internal structure and acceleration of GRB jets. The GW signal peaks just outside the jet (core) of a uniform (structuted) jet. Within the jet (core) the GW signal displays wiggles, due to a polarization effect; thus for a uniform jet, the peak amplitude accompanies a smoother signal than the peak of a structured jet. For the most probable detection angle and for typical GRB parameters, we expect frequencies < ~600Hz and amplitudes h~10^-25. Our estimates of the expected signals suggest that the signals are not strong enough for a single cluster of DECIGO nor for aLIGO's sensitivities. However, sensitivies of ~10^-25/sqrt(Hz) in the DECIGO band should suffice to detect typical long GRBs at 2Gpc and short GRBs at 200Mpc, implying a monthly event of a long GRB and a detection of a short GRB every decade. In addition, we expect much more frequent detection of GW from GRBs directed away from us, including orphan afterglows. The ultimate DECIGO sensitivy should increase the range and enable detecting these signals in all models even to high cosmological z.

preprint2013arXiv

On the origin of the radio emission of Sw 1644+57

We apply relativistic equipartition synchrotron arguments to the radio data of the tidal disruption event candidate Sw 1644+57. We find that, regardless of the details of the equipartition scenario considered, the energy required to produce the observed radio (i.e., energy in magnetic field and radio emitting electrons) must increase by a factor of ~20 during the first 200 days. It then saturates. This energy increase cannot be alleviated by a varying geometry of the system. The radio data can be explained by: (i) An afterglow like emission of the X-ray emitting narrow relativistic jet. The additional energy can arise here from a slower moving material ejected in the first few days that gradually catches up with the slowing down blast wave (Berger et al. 2012). However, this requires at least ~4x10^{53} erg in the slower moving outflow. This is much more than the energy of the fast moving outflow that produced the early X-rays and it severely constrains the overall energy budget. (ii) Alternatively, the radio may arise from a mildly relativistic quasi-spherical outflow. Here, the energy for the radio emission increases with time to at least ~10^{51} erg after 200 days. This scenario requires, however, a second X-ray emitting collimated relativistic component. Given these results, it is worthwhile to consider models in which the energy of the magnetic field and/or of the radio emitting electrons increases with time without a continuous energy supply to the blast wave. This can happen, for example, if the energy is injected initially mostly in one form (Poynting flux or baryonic) and it is gradually converted to the other form, leading to a deviation from equipartition. Another intriguing possibility is that a gradually decreasing Inverse Compton cooling modifies the synchrotron emission and leads to an increase of the available energy in the radio emitting electrons (Kumar et al. 2013).

preprint2013arXiv

Possible relevance of quantum spacetime for neutrino-telescope data analyses

One of the primary goals of neutrino telescopes, such as IceCube, is the discovery of neutrinos emitted by gamma-ray bursts (GRBs). Another source of interest in the results obtained by these telescopes is their possible use for tests of the applicability of Einstein's Special Relativity to neutrinos, particularly with respect to modifications that lead to Lorentz invariance violation that have been conjectured by some models of quantum space-time. We examine here the fascinating scenario in which these two aspects of neutrino-telescope physics require a combined analysis. We discuss how neutrinos that one would not associate to a GRB, when assuming a classical spacetime picture, may well be GRB neutrinos if the possibility that Lorentz invariance is broken at very high energies is taken into account. As an illustrative example we examine three IceCube high energy neutrinos that arrived hours before GRBs (but from the same direction) and we find that the available, IceCube data, while inconclusive, is compatible with a scenario in which one or two of these neutrinos were GRB neutrinos and their earlier arrival reflects Lorentz invariance violation. We outline how future analyses of neutrino data should be done in order to systematically test this possibility.

preprint2013arXiv

Radius constraints and minimal equipartition energy of relativistically moving synchrotron sources

A measurement of the synchrotron self-absorption flux and frequency provides tight constraints on the physical size of the source and a robust lower limit on its energy. This lower limit is also a good estimate of the magnetic field and electrons' energy, if the two components are at equipartition. This well-known method was used for decades to study numerous astrophysical sources moving at non-relativistic (Newtonian) speeds. Here we generalize the Newtonian equipartition theory to sources moving at relativistic speeds including the effect of deviation from spherical symmetry expected in such sources. Like in the Newtonian case, minimization of the energy provides an excellent estimate of the emission radius and yields a useful lower limit on the energy. We find that the application of the Newtonian formalism to a relativistic source would yield a smaller emission radius, and would generally yield a larger lower limit on the energy (within the observed region). For sources where the Synchrotron-self-Compton component can be identified, the minimization of the total energy is not necessary and we present an unambiguous solution for the parameters of the system.

preprint2013arXiv

Reconciling the GRB rate and star formation histories

While there are numerous indications that GRBs arise from the death of massive stars, the GRB rate does not follow the global cosmic star formation rate and, within their hosts, GRBs are more concentrated in regions of very high star formation. We explain both puzzles here. Using the publicly available VESPA database of SDSS Data Release 7 spectra, we explore a multi-parameter space in galaxy properties, like stellar mass, metallicity, dust etc. to find the sub-set of galaxies that reproduce the recently obtained GRB rate by Wanderman & Piran (2010). We find that only galaxies with present stellar masses below < 10^{10} Msun and low metallicity reproduce the observed GRB rate. This is consistent with direct observations of GRB hosts and provides an independent confirmation of the nature of GRB hosts. Because of the significantly larger sample of SDSS galaxies, we compute their correlation function and show that they are anti-biased with respect to the dark matter: they are in filaments and voids. Using recent observations of massive stars in local dwarfs we show how the fact that GRB hosts galaxies are dwarfs can explain the observation that GRBs are more concentrated in regions of high star formation than SNe. Finally we explain these results using new theoretical advances in the field of star formation.

preprint2013arXiv

The Discovery of a New Instability in a Hyperaccretion Flow and its Implication for Gamma-ray Bursts

A hyperaccretion flow around a stellar mass black hole is thought to be the most plausible engine that powers gamma-ray bursts (GRBs). The flow efficiently cools via neutrino emission at >~ 0.003-0.01 M_sun s^{-1} (corresponding to a luminosity of ~10^50 erg s^{-1}), while neither neutrino nor photon emission is efficient below this rate, so the flow should be advection dominated. We carefully solve how a transition occurs from the advection-dominated to the neutrino-dominated branches, and find that the slope of the thermal equilibrium curve is negative in the surface density - accretion rate (Sigma-Mdot) plane, a condition for viscous instability, at radii smaller than ~12R_g (with R_g being the gravitational radius). We also confirm that the flow is thermally stable. The consequence of this instability is the formation of a clumpy structure in the flow. This is because the larger (respectively smaller) surface density is, the smaller (respectively larger) the mass accretion rate from the region in question becomes, leading to growth of the density contrast. The timescale for clump formation is estimated to be shorter than 0.1s. The observational implication is discussed in the context of GRBs. We suggest that this might explain the origin of the large variability observed in the prompt emission of GRBs.

preprint2012arXiv

Choptuik's Critical Phenomenon in Einstein-Gauss-Bonnet Gravity

We investigate the effects of higher order curvature corrections to Einstein's Gravity on the critical phenomenon near the black hole threshold, namely the Choptuik phenomenon. We simulate numerically a five dimensional spherically symmetric gravitational collapse of massless scalar field in Einstein-Gauss-Bonnet gravity towards a black hole formation threshold. When the curvature is sufficiently large the additional higher order terms, affect the evolution of the whole system. Since high curvature characterizes the region when the critical behavior takes place this critical behavior is destroyed. Both the self similarity and the mass scaling relation disappear. Instead we find a different behavior near the black hole threshold, which depends on the coupling constant of the higher order terms. The new features include a change of the sign of the Ricci scalar on the origin which indicates changes in the local geometry of space-time, and never occurs in classical general relativity collapse, and oscillations with a constant rather than with a diminishing length scale.

preprint2012arXiv

Collisional Penrose process near the horizon of extreme Kerr black holes

Collisions of particles in black holes' ergospheres may result in an arbitrarily large center of mass energy. This led recently to the suggestion (Banados et al., 2009) that black holes can act as ultimate particle accelerators. If the energy of an outgoing particle is larger than the total energy of the infalling particles the energy excess must come from the rotational energy of the black hole and hence this must involve a Penrose process. However, while the center of mass energy diverges the position of the collision makes it impossible for energetic particles to escape to infinity. Following an earlier work on collisional Penrose processes (Piran & Shaham 1977) we show that even under the most favorable idealized conditions the maximal energy of an escaping particle is only a modest factor above the total initial energy of the colliding particles. This implies that one shouldn't expect collisions around a black hole to act as spectacular cosmic accelerators.

preprint2012arXiv

Implications of the UHECRs penetration depth measurements

The simple interpretation of PAO's UHECRs' penetration depth measurements suggests a transition at the energy range $1.1 - 35 \cdot 10^{18} $ eV from protons to heavier nuclei. A detailed comparison of this data with air shower simulations reveals strong restrictions on the amount of light nuclei (protons and He) in the observed flux. We find a robust upper bound on the observed proton fraction of the UHECRs flux and we rule out a composition dominated by protons and He. Acceleration and propagation effects lead to an observed composition that is different from the one at the source. Using a simple toy model that take into account these effects, we show that the observations requires an extreme metallicity at the sources with metals to protons mass ratio of 1:1, a ratio that is larger by a factor of a hundred than the solar abundance. This composition imposes an almost impossible constraint on all current astrophysical models for UHECRs accelerators. This may provide a first hint towards new physics that emerges at $\sim 100$ TeV and leads to a larger proton cross section at these energies.

preprint2012arXiv

Jet Luminosity From Neutrino-Dominated Accretion Flows in Gamma-Ray Bursts

A hyperaccretion disk formed around a stellar mass black hole is a plausible model for the central engine that powers gamma-ray bursts (GRBs). If the central black hole rotates and a poloidal magnetic field threads its horizon, a powerful relativistic jet may be driven by a process resembling the Blandford-Znajek mechanism. We estimate the luminosity of such a jet assuming that the poloidal magnetic field strength is comparable to the inner accretion disk pressure. We show that the jet efficiency attains its maximal value when the accretion flow is cooled via optically-thin neutrino emission. The jet luminosity is much larger than the energy deposition through neutrino-antineutrino annihilation provided that the black hole is spinning rapidly enough. When the accretion rate onto a rapidly spinning black hole is large enough (> 0.003-0.01M_sun/sec), the predicted jet luminosity is sufficient to drive a GRB.

preprint2012arXiv

Jets from Tidal Disruptions of Stars by Black Holes

Tidal disruption of main sequence stars by black holes has generally been thought to lead to a signal dominated by UV emission. If, however, the black hole spins rapidly and the poloidal magnetic field intensity on the black hole horizon is comparable to the inner accretion disk pressure, a powerful jet may form whose luminosity can easily exceed the thermal UV luminosity. When the jet beam points at Earth, its non-thermal luminosity can dominate the emitted spectrum. The thermal and non-thermal components decay differently with time. In particular, the thermal emission should remain roughly constant for a significant time after the period of maximum accretion, beginning to diminish only after a delay, whereas after the peak accretion rate, the non-thermal jet emission decays, but then reaches a plateau. Both transitions are tied to a characteristic timescale $t_{\rm Edd}$ at which the accretion rate falls below Eddington. Making use of this timescale in a new parameter-inference formalism for tidal disruption events with significant emission from a jet, we analyze the recent flare source Swift J2058. It is consistent with an event in which a main sequence solar-type staris disrupted by a black hole of mass $\sim 4 \times 10^7 M_{\odot}$. The beginning of the flat phase in the non-thermal emission from this source can possibly be seen in the late-time lightcurve. Optical photometry over the first $\simeq 40$ d of this flare is also consistent with this picture, but is only weakly constraining because the bolometric correction is very uncertain. We suggest that future searches for main sequence tidal disruptions use methods sensitive to jet radiation as well as to thermal UV radiation.

preprint2012arXiv

On the lateral expansion of GRB jets

The dynamics of GRB jets during the afterglow phase have an important effect on the interpretation of their observations and for inferring key physical parameters such as their true energy and event rate. Semi-analytic models generally predict a fast lateral expansion, where the jet opening angle asymptotically grows exponentially with its radius. Numerical simulations, however, show a much more modest lateral expansion, where the jet retains memory of its initial opening angle for a very long time, and the flow remains non-spherical until it becomes sub-relativistic, and only then gradually approaches spherical symmetry. Here we suggest a new analytic model based on a new physically derived recipe for the lateral expansion. We also generalize the model by relaxing the common approximations of ultra-relativistic motion and a narrow jet opening angle. We find that the new analytic model fits much better the results of numerical simulations, mainly because it remains valid also in the mildly relativistic, quasi spherical regime. This model shows that for modest initial jet half-opening angles, θ_0, the outflow is not sufficiently ultra-relativistic when its Lorentz factor reaches Γ= 1/θ_0 and therefore the sideways expansion is rather slow, showing no rapid, exponential phase. On the other hand, we find that jets with an extremely narrow initial half-opening angle, of about θ_0 << 10^{-1.5} or so, which are still sufficiently ultra-relativistic at Γ= 1/θ_0, do show a phase of rapid, exponential lateral expansion. However, even such jets that expand sideways exponentially are still not spherical when they become sub-relativistic.

preprint2012arXiv

On thermalization in gamma-ray burst jets and the peak energies of photospheric spectra

The low energy spectral slopes of the prompt emission of most gamma-ray bursts (GRBs) are difficult to reconcile with radiatively efficient optically thin emission models irrespective of the radiation mechanism. An alternative is to ascribe the radiation around the spectral peak to a thermalization process occurring well inside the Thomson photosphere. This quasi-thermal spectrum can evolve into the observed non-thermal shape by additional energy release at moderate to small Thomson optical depths, which can readily give rise to the hard spectral tail. The position of the spectral peak is determined by the temperature and Lorentz factor of the flow in the termalization zone, where the total number of photons carried by the jet is established. To reach thermalization, dissipation alone is not sufficient and photon generation requires an efficient emission/absorption process in addition to scattering. We perform a systematic study of all relevant photon production mechanisms searching for possible conditions in which thermalization can take place. We find that a significant fraction of the available energy should be dissipated at intermediate radii, $\sim 10^{10}$-- a few$\times 10^{11}$ cm and the flow there should be relatively slow: the bulk Lorentz factor could not exceed a few tens for all but the most luminous bursts with the highest $\Epk$-s. The least restrictive constraint for successful thermalization, $Γ\lesssim 20$, is obtained if synchrotron emission acts as the photon source. This requires, however, a non-thermal acceleration deep below the Thomson photosphere transferring a significant fraction of the flow energy to relativistic electrons with Lorentz factors between 10 and 100. Other processes require bulk flow Lorentz factors of order of a few for typical bursts. We examine the implications of these results to different GRB photospheric emission models.

preprint2012arXiv

Short vs Long and Collapsars vs. non-Collapsar: a quantitative classification of GRBs

Gamma-Ray Bursts (GRBs) are traditionally divided to long and short according to their durations (>/< 2 sec). It was generally believed that this reflects a different physical origin: Collapsars (long) and non-Collapsars (short). We have recently shown that the duration distribution of Collapsars is flat, namely independent of the duration, at short durations. Using this model for the distribution of Collapsars we determine the duration distribution of non-Collapsars and estimate the probability that a burst with a given duration (and hardness) is a Collapsar or not. We find that this probability depends strongly on the spectral window of the observing detector. While the commonly used limit of 2 sec is conservative and suitable for BATSE bursts, 40% of swift's bursts shorter than 2 sec are Collapsars and division >/<0.8 sec is more suitable for swift. We find that the duration overlap of the two populations is very large. On the one hand there is a non-negligible fraction of non-Collapsars longer than 10 sec, while on the other hand even bursts shorter than 0.5 sec in the swift sample have a non-negligible probability to be Collapsars. Our results enable the construction of non-Collapsar samples while controlling the Collapsar contamination. They also highlight that no firm conclusions can be drawn based on a single burst and they have numerous implications concerning previous studies of non-Collapsar properties that were based on the current significantly contaminated swift samples of localized short GRBs. Specifically: (i) all known short bursts with z>1 are most likely Collapsars, (ii) the only short burst with a clear jet break is most likely a Collapsar, indicating our lack of knowledge concerning non-Collapsar beaming (iii) the existence of non-Collapsars with durations up to 10 sec impose new challenges to non-Collapsar models.

preprint2012arXiv

The Long, the Short and the Weak - the origin of GRBs

The origin of Gamma-Ray Bursts is one of the most interesting puzzles in recent astronomy. During the last decade a consensus formed that long GRBs (LGRBs) arise from the collapse of massive stars and that short GRBs (SGRBs) have a different origin, most likely neutron star mergers. A key ingredient of the Collapsar model that explains how the collapse of massive stars produces a GRB is the emergence of a relativistic jet that penetrates the stellar envelope. The condition that the emerging jet penetrates the envelope poses strong constraints on the system. Using these constraints we show that: (i) Low luminosity GRBs (llGRBs), a sub population of GRBs with a very low luminosities (and other peculiar properties: single peaked, smooth and soft) cannot be formed by Collapsars. llGRBs must have a different origin (most likely a shock breakout). (ii) On the other hand regular LGRBs must be formed by Collapsars. (iii) While for BATSE the dividing duration between Collapsars and non-Collapsar is indeed at $\sim 2$ sec, the dividing duration is different for other GRBs detectors. In particular most Swift bursts longer than 0.8 sec are of a Collapsar origin. This last results requires a revision of many conclusions concerning the origin of Swift SGRBs which were based on the commonly used 2 sec limit.

preprint2012arXiv

The Swift short gamma-ray burst rate density: implications for binary neutron star merger rates

Short gamma-ray bursts (SGRBs) observed by {\it Swift} are potentially revealing the first insight into cataclysmic compact object mergers. To ultimately acquire a fundamental understanding of these events requires pan-spectral observations and knowledge of their spatial distribution to differentiate between proposed progenitor populations. Up to April 2012 there are only some 30% of SGRBs with reasonably firm redshifts, and this sample is highly biased by the limited sensitivity of {\it Swift} to detect SGRBs. We account for the dominant biases to calculate a realistic SGRB rate density out to $z\approx0.5$ using the {\it Swift} sample of peak fluxes, redshifts, and those SGRBs with a beaming angle constraint from X-ray/optical observations. We find an SGRB lower rate density of $8^{+5}_{-3}$ $\mathrm{Gpc}^{-3}\mathrm{yr}^{-1}$ (assuming isotropic emission), and a beaming corrected upper limit of $1100^{+700}_{-470} $ $\mathrm{Gpc}^{-3}\mathrm{yr}^{-1}$. Assuming a significant fraction of binary neutron star mergers produce SGRBs, we calculate lower and upper detection rate limits of $(1-180)$ yr$^{-1}$ by an aLIGO and Virgo coincidence search. Our detection rate is similar to the lower and realistic rates inferred from extrapolations using Galactic pulsar observations and population synthesis.

preprint2012arXiv

The Swift short gamma-ray burst rate density: prospects for detecting binary neutron star mergers by aLIGO

Presently only 30% of short gamma ray bursts (SGRBs) have accurate redshifts, and this sample is highly biased by the limited sensitivity of {\it Swift} to detect SGRBs. We account for the dominant biases to calculate a realistic SGRB rate density out to $z = 0.5$ using the {\it Swift} sample of peak fluxes, redshifts, and those SGRBs with a beaming angle constraint from X-ray/optical observations. Assuming a significant fraction of binary neutron star mergers produce SGRBs, we calculate lower and upper detection rate limits of (1-180) per Yr by an advanced LIGO and Virgo coincidence search. Our detection rate is compatible with extrapolations using Galactic pulsar observations and population synthesis.

preprint2012arXiv

Time Scales in Tidal Disruption Events

We explore the temporal structure of tidal disruption events pointing out the corresponding transitions in the lightcurves of the thermal accretion disk and of the jet emerging from such events. The hydrodynamic time scale of the disrupted star is the minimal time scale of building up the accretion disk and the jet and it sets a limit on the rise time. This suggest that Swift J1644+57, that shows several flares with a rise time as short as a few hundred seconds could not have arisen from a tidal disruption of a main sequence star whose hydrodynamic time is a few hours. The disrupted object must have been a white dwarf. A second important time scale is the Eddington time in which the accretion rate changes form super to sub Eddington. It is possible that such a transition was observed in the light curve of Swift J2058+05. If correct this provides intersting constraints on the parameters of the system.

preprint2011arXiv

An observational imprint of the Collapsar model of long Gamma Ray Bursts

The Collapsar model provides a theoretical framework for the well known association between long gamma-ray bursts (GRBs) and collapsing massive stars. A bipolar relativistic jet, launched at the core of a collapsing star, drills its way through the stellar envelope and breaks out of the surface before producing the observed gamma-rays. While a wealth of observations associate GRBs with the death of massive stars, as yet there is no direct evidence for the Collapsar model itself. Here we show that a distinct signature of the Collapsar model is the appearance of a plateau in the duration distribution of the prompt GRB emission at times much shorter than the typical breakout time of the jet. This plateau is evident in the data of all three major satellites. These findings provide an evidence that directly supports the Collapsar model. Additionally, it suggests the existence of a large population of choked (failed) GRBs and that the 2 s duration commonly used to separate Collapsars and non-Collapasars is inconsistent with the duration distributions of {\it Swift} and Fermi GRBs and only holds for BATSE GRBs.

preprint2011arXiv

Are low luminosity GRBs generated by relativistic jets?

Low luminosity gamma-ray bursts (ll-GRBs) constitute a sub-class of gamma-ray bursts (GRBs) that plays a central role in the GRB-supernova connection. While ll-GRBs differ from typical long GRBs (LGRBs) in many aspects, they also share some common features. Therefore, the question whether the gamma-ray emission of ll-GRBs and LGRBs has a common origin is of great interest. Here we address this question by testing whether ll-GRBs, like LGRBs according to the Collapsar model, can be generated by relativistic jets that punch holes in the envelopes of their progenitor stars. The collapsar model predicts that the durations of most observed bursts will be comparable to, or longer than, the time it takes the jets to breakout of the star. We calculate the jet breakout times of ll-GRBs and compare them to the observed durations. We find that there is a significant access of ll-GRBs with durations that are much shorter than the jet breakout time and that these are inconsistent with the Collapsar model. We conclude that the processes that dominate the gamma-ray emission of ll-GRBs and of LGRBs are most likely fundamentally different.

preprint2011arXiv

Effects of a low electron distribution cutoff on multiwavelength spectra and light curves of GRB afterglows

Aims: We investigate the behavior of the frequency-centered light curves expected within the standard model of Gamma Ray Bursts allowing the maximum electron energy to be a free parameter permitted to take low values. Methods: We solve the spatially averaged kinetic equations which describe the simultaneous evolution of particles and photons, obtaining the multi-wavelength spectra as a function of time. From these we construct the frequency-centered light curves giving emphasis in the X-ray and optical bands. Results: We show that in cases where the maximum electron energy takes low values, the produced X-ray light curves show a plateau as the synchrotron component gives its place to the Synhro Self-Compton one in the X-ray band.

preprint2011arXiv

New insights on the duration distribution of long GRBs from Collapsars

According the Collapsar model long gamma-ray bursts (LGRBs) involve relativistic jets that puncture the envelope of a collapsing star, and produced the γ-rays after they break out. This model provides a theoretical framework for the well known association between LGRBs and massive stars. However although this association is supported by a wealth of observations, to this date there is no direct observational evidence for the emergence of the jet from the star. In other words there is no direct evidence for the Collapsar model. Here we show that a distinct signature of the Collapsar model is the appearance of a plateau in the duration distribution of the prompt GRB emission at times much shorter than the typical breakout time of the jet. This plateau is evident in the data of all major GRB satellites, and provides a direct evidence supporting the Collapsar model. It also enables us to place limits on the sizes and masses of LGRB progenitors; suggests the existence of a large population of choked (failed) GRBs; and indicates that the 2 s duration commonly used to separate Collapsars and non-Collapsars holds for BATSE and possibly Fermi GBM GRBs, but it is inconsistent with the duration distributions of Swift GRBs.

preprint2011arXiv

Radio Remnants of Compact Binary Mergers - the Electromagnetic Signal that will follow the Gravitational Waves

The question "what is the observable electromagnetic (EM) signature of a compact binary merger?" is an intriguing one with crucial consequences to the quest for gravitational waves (GW). Compact binary mergers are prime sources of GW, targeted by current and next generation detectors. Numerical simulations have demonstrated that these mergers eject energetic sub-relativistic (or even relativistic) outflows. This is certainly the case if the mergers produce short GRBs, but even if not, significant outflows are expected. The interaction of such outflows with the surround matter inevitably leads to a long lasting radio signal. We calculate the expected signal from these outflows (our calculations are also applicable to short GRB orphan afterglows) and we discuss their detectability. We show that the optimal search for such signal should, conveniently, take place around 1.4 GHz. Realistic estimates of the outflow parameters yield signals of a few hundred $μ$Jy, lasting a few weeks, from sources at the detection horizon of advanced GW detectors. Followup radio observations, triggered by GW detection, could reveal the radio remnant even under unfavorable conditions. Upcoming all sky surveys can detect a few dozen, and possibly even thousands, merger remnants at any give time, thereby providing robust merger rate estimates even before the advanced GW detectors become operational. In fact, the radio transient RT 19870422 fits well the overall properties predicted by our model and we suggest that its most probable origin is a compact binary merger radio remnant.

preprint2011arXiv

Swift J1644+57: A White Dwarf Tidally Disrupted by a 10^4 M_{odot} Black Hole?

We propose that the remarkable object Swift J1644+57, in which multiple recurring hard X-ray flares were seen over a span of several days, is a system in which a white dwarf was tidally disrupted by an intermediate mass black hole. Disruption of a white dwarf rather than a main sequence star offers a number of advantages in understanding the multiple, and short, timescales seen in the light curve of this system. In particular, the short internal dynamical timescale of a white dwarf offers a more natural way of understanding the short rise times ($\sim 100$ s) observed. The relatively long intervals between flares ($\sim 5 \times 10^4$ s) may also be readily understood as the period between successive pericenter passages of the remnant white dwarf. In addition, the expected jet power is larger when a white dwarf is disrupted. If this model is correct, the black hole responsible must have mass $\lesssim 10^5 M_{\odot}$.

preprint2011arXiv

The propagation of relativistic jets in external media

Relativistic jets are ubiquitous in astrophysical systems that contain compact objects. They transport large amounts of energy to large distances from the source, and their interaction with the ambient medium has a crucial effect on the evolution of the system. The propagation of the jet is characterized by the formation of a shocked "head" at the front of the jet which dissipates the jet's energy and a cocoon that surrounds the jet and potentially collimates it. We present here a self consistent, analytic model that follows the evolution of the jet and its cocoon, and describes their interaction. We show that the critical parameter that determines the properties of the jet-cocoon system is the dimensionless ratio between the jet's energy density and the rest-mass energy density of the ambient medium. This parameter, together with the jet's injection angle, also determines whether the jet is collimated by the cocoon or not. The model is applicable to relativistic, unmagnetized, jets on all scales and may be used to determine the conditions in AGNs jets as well as in GRBs or microquasars. It shows that AGN and microquasar jets are hydrodynamically collimated due to the interaction with the ambient medium, while GRB jets can be collimated only inside a star and become uncollimated once they breakout.

preprint2011arXiv

Variability in Blazars: Clues from PKS 2155-304

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

preprint2010arXiv

A New Limit on the Distances of Nuclei UHECRs Sources

Recent evidence from the Pierre Auger Observatory suggests a transition, at 5 EeV-10EeV in the composition of Ultra High Energy Cosmic Rays (UHECRs), from protons to heavier nuclei such as iron. I consider here the implications of the heavier composition on the sources of UHECRs. The smaller magnetic rigidity implies that nuclei UHECRs are: (i) More easily accelerated (ii) Local, as they can diffuse only a few Mpc from their sources before disintegrating (iii) Isotropic, because large deflections in the extra galactic and the galactic magnetic fields erase the directional information. Uncertainties in the strength and structure of the the extra galactic magnetic field (EGMF) makes it difficult to estimate the overall effects. However, with typical reasonable parameters of a few nG and a coherence distance of a Mpc the distance a nuclei UHECR above the GZK energy traverses before photodisintegrating is only a few Mpc. In spite of the significantly weaker limits on the luminosity, Cen A is the only currently active potential source of nuclei UHECRs within this distance. The large deflections erases the directional anisotropy expected from a single source. If indeed the composition of above-GZK-UHECRs is iron and if the EGMF is not too small then Cen A is the dominant source of observed nuclei UHECRs above the GZK limit.

preprint2010arXiv

A Revised Limit of the Lorentz Factors of GRBs with Two Emitting Regions

Fermi observations of GeV emission from GRBs have suggested that the Lorentz factor of some GRBs is around a thousand or even higher. At the same time the same Fermi observations have shown an extended GeV emission indicating that this higher energy emission might be a part of the afterglow and it does not come from the same region as the lower energy prompt emission. If this interpretation is correct than we should reconsider the opacity limits on the Loretnz factor of the emitting regions which are based on a one-zone model. We describe here a two-zone model in which the GeV photons are emitted in a larger radius than the MeV photons and we calculate the optical depth for pair creation of a GeV photon passing the lower energy photons shell. We find that, as expected, the new two-zone limits on the Lorentz factor are significantly lower. The corresponding limits for the Fermi bursts are lower by a factor of five compared to the one-zone model and it is possible that both the MeV and GeV regions have relatively modest Lorentz factors (~200 - 400).

preprint2010arXiv

Modifications to Lorentz invariant dispersion in relatively boosted frames

We investigate the implications of energy-dependence of the speed of photons, one of the candidate effects of quantum-gravity theories that has been most studied recently, from the perspective of observations in different reference frames. We examine how a simultaneous burst of photons would be measured by two observers with a relative velocity, establishing some associated conditions for the consistency of theories. For scenarios where the Lorentz transformations remain valid these consistency conditions allow us to characterize the violations of Lorentz symmetry through an explicit description of the modification of the quantum-gravity scale in boosted frames with respect to its definition in a preferred frame. When applied to relativistic scenarios with a deformation of Lorentz invariance that preserves the equivalence of inertial observers, we find an insightful characterization of the necessity to adopt in such frameworks non-classical features of spacetime geometry, e.g. events that are at the same spacetime point for one observer cannot be considered at the same spacetime point for other observers. Our findings also suggest that, at least in principle (and perhaps one day even in practice), measurements of the dispersion of photons in relatively boosted frames can be particularly valuable for the purpose of testing these scenarios.

preprint2010arXiv

On The External Shock Synchrotron Model for GRBs' GeV Emission

The dominant component of the (100 MeV - 50 GeV) GRB emission detected by LAT starts with a delay relative to the prompt soft (sub-MeV) gamma-rays and lasts long after the soft component fades. This has lead to the intriguing suggestion that this high energy emission is generated via synchrotron emission of relativistic electrons accelerated by the external shock. Moreover, the limits on the MeV afterglow emission lead to the suggestion that, at least in bright GeV bursts the field is not amplified beyond compression in the shock. We show here that considerations of confinement (within the decelerating shock), efficiency and cooling of the emitting electrons constrain, within this model, the magnetic fields that arise in both the upstream (circum burst) and downstream (ejecta) regions, allowing us to obtain a direct handle on their values. The well known limit on the maximal synchrotron emission, when combined with the blast wave evolution, implies that late photons (arriving more than ~100 s after the burst) with energies higher than ~ 10GeV do not arise naturally from external shock synchrotron and almost certainly have a different origin. Finally, even a modest seed flux (a few mJy) at IR-optical would quench, via Inverse Compton cooling, the GeV emission unless the magnetic field is significantly amplified behind the shock. An observation of a burst with simultaneous IR-optical and GeV emission will rule out this model.

preprint2010arXiv

The late jet in gamma-ray bursts and its interactions with a supernova ejecta and a cocoon

Late X-ray flares observed in X-ray afterglows of gamma-ray bursts (GRBs) suggest late central engine activities at a few minuets to hours after the burst. A few unambiguously confirmed cases of supernova associations with nearby long GRBs imply that an accompanying supernova-like component might be a common feature in all long GRB events. These motivate us to study the interactions of a late jet, responsible for a x-ray flare, with various components in a stellar explosion, responsible for a GRB. These components include a supernova shell-like ejecta, and a cocoon that was produced when the main jet producing the GRB itself was propagating through the progenitor star. We find that the interaction between the late jet and the supernova ejecta may produce a luminous (up to 10^49 erg s^-1) thermal X-ray transient lasting for ~ 10 s. The interaction between the late jet and the cocoon produces synchrotron self-absorbed non-thermal emission, with the observed peak X-ray flux density from 0.001 micro Jy to 1 mJy at 1 keV and a peak optical flux density from 0.01 micro Jy to 0.1 Jy (for a redshift z= 2). The light curve due to the late jet - cocoon interaction has very small pulse-width-to-time ratio, Δt / t \approx 0.01 - 0.5, where t is the pulse peak time since the burst trigger. Identifying these features in current and future observations would open a new frontier in the study of GRB progenitor stars.

preprint2010arXiv

The luminosity function and the rate of Swift's Gamma Ray Bursts

We invert directly the redshift - luminosity distribution of observed long Swift GRBs to obtain their rate and luminosity function. Our best fit rate is described by a broken power law that rises like (1+z)^2.1{+0.5-0.6} for 0<z<3 and decrease like (1+z)^-1.4{+2.4-1.0} for z>3. The local rate is 1.3^{+0.6-0.7} [Gpc^-3 yr^-1]. The luminosity function is well described by a broken power law with a break at L* = 10^52.5{+-0.2}[erg/sec] and with indices alpha = 0.2^{+0.2-0.1} and beta = 1.4^{+0.3-0.6}. The recently detected GRB 090423, with redshift ~8, fits nicely into the model's prediction, verifying that we are allowed to extend our results to high redshifts. While there is a possible agreement with the star formation rate (SFR) for z<3, the high redshift slope is shallower than the steep decline in the SFR for 4<z. However we cannot rule out a GRB rate that follows one of the recent SFR models.

preprint2009arXiv

Astro2010 Decadal Survey Whitepaper: Coordinated Science in the Gravitational and Electromagnetic Skies

It is widely expected that the coming decade will witness the first direct detection of gravitational waves (GWs). The ground-based LIGO and Virgo GW observatories are being upgraded to advanced sensitivity, and are expected to observe a significant binary merger rate. The launch of The Laser Interferometer Space Antenna (LISA) would extend the GW window to low frequencies, opening new vistas on dynamical processes involving massive (M >~ 10^5 M_Sun) black holes. GW events are likely to be accompanied by electromagnetic (EM) counterparts and, since information carried electromagnetically is complementary to that carried gravitationally, a great deal can be learned about an event and its environment if it becomes possible to measure both forms of radiation in concert. Measurements of this kind will mark the dawn of trans-spectral astrophysics, bridging two distinct spectral bands of information. The aim of this whitepaper is to articulate future directions in both theory and observation that are likely to impact broad astrophysical inquiries of general interest. What will EM observations reflect on the nature and diversity of GW sources? Can GW sources be exploited as complementary probes of cosmology? What cross-facility coordination will expand the science returns of gravitational and electromagnetic observations?

preprint2009arXiv

GRB Light Curves in the Relativistic Turbulence Model

Randomly oriented relativistic emitters in a relativistically expanding shell provides an alternative to internal shocks as a mechanism for producing GRBs' variable light curves with efficient conversion of energy to radiation. In this model the relativistic outflow is broken into small emitters moving relativistically in the outflow's rest frame. Variability arises because an observer sees an emitter only when its velocity points towards him so that only a small fraction of the emitters are seen by a given observer. Models with significant relativistic random motions require converting and maintaining a large fraction of the overall energy into these motions. While it is not clear how this is achieved, we explore here, using two toy models, the constraints on parameters required to produce light curves comparable to the observations. We find that a tight relation between the size of the emitters and the bulk and random Lorentz factors is needed and that the random Lorentz factor determines the variability. While both models successfully produce the observed variability there are several inconsistencies with other properties of the light curves. Most of which, but not all, might be resolved if the central engine is active for a long time producing a number of shells, resembling to some extent the internal shocks model.

preprint2009arXiv

Lorentz Factor Constraint from the very early external shock of the gamma-ray burst ejecta

While it is generally agreed that the emitting regions in Gamma-Ray Bursts (GRBs) move ultra relativistically towards the observer, different estimates of the initial Lorentz factors, $Γ_0$, lead to different, at times conflicting estimates. We show here that the quiet periods in which the signals goes down below the instrumental thresholds, put strong upper limits on the values of $Γ_0$. According to the standard internal-external shocks model an external shock should develop during the prompt stage. This external shock radiates in the hard X-rays to soft gamma-rays bands and this emission should be seen as a smooth background signal. The observed deep minima indicate that this contribution is negligible. This limits, in turn, $Γ_0$. We obtain upper limits on $Γ_0$ for several bursts with typical values around hundreds. We compare these values with those obtained by the other methods, which typically yield lower limits. The results are marginally consistent leaving only a narrow range of allowed values for $Γ_0$.

preprint2006arXiv

GRB 051221A and Tests of Lorentz Symmetry

Various approaches to quantum gravity suggest the possibility of violation of Lorentz symmetry at very high energies. In these cases we expect a modification at low energies of the dispersion relation of photons that contains extra powers of the momentum suppressed by a high energy scale. These terms break boost invariance and can be tested even at relatively low energies. We use the light curves of the very bright short Gamma-Ray Burst GRB 051221A and compare the arrival times of photons at different energies with the expected time delay due to a modified dispersion relation. As no time delay was observed, we set a lower bound of 0.0066 E_{pl} \sim 0.66 10^{17} GeV on the scale of Lorentz invariance violation.

preprint2005arXiv

Early GRB afterglow from a reverse shock as a tracer of the prompt gamma-ray light curve

We discuss the optical and radio early afterglow emission of the reverse shock that crosses a baryonic ejecta as it interacts with the external interstellar medium (ISM). We show that the peak of the optical flash divides the light curve of the reverse shock into two distinctive phases. The emission after the peak depends weakly on the initial conditions of the ejecta and therefore it can be used as an identifiable signature of a reverse shock emission. On the other hand, the emission before the optical peak is highly sensitive to the initial conditions and therefore can be used to investigate the initial hydrodynamic profile of the ejecta. In particular, if the prompt $γ$-ray emission results from internal shocks, the early reverse shock emission should resemble a smoothed version of the prompt $γ$-ray light curve.

preprint2005arXiv

The Beaming Factor and Other Open Questions In GRB Jets

I review several central open questions concerning GRB Jets. 1) I discuss a new estimates of the beaming correction for the rate of GRBs ~ 75 +- 20. 2) I discuss the universal structured jet (USJ) model and conclude that while jets might be structured they are less likely to be universal. 3) I discuss recent observations of a sideways expansion of a GRB afterglow and compare these with current hydrodynamics simulations of jet evolution. 4) I discuss the implications of resent outliers to the energy-angle relation.

preprint2005arXiv

The distances of short-hard GRBs and the SGR connection

We present a search for nearby (D<100 Mpc) galaxies in the error boxes of six well-localized short-hard gamma-ray bursts (GRBs). None of the six error boxes reveals the presence of a plausible nearby host galaxy. This allows us to set lower limits on the distances and, hence, the isotropic-equivalent energy of these GRBs. Our lower limits are around $1 \times 10^{49}$ erg (at $2σ$ confidence level); as a consequence, some of the short-hard GRBs we examine would have been detected by BATSE out to distances greater than 1 Gpc and therefore constitute a bona fide cosmological population. Our search is partially motivated by the December 27, 2004 hypergiant flare from SGR 1806-20, and the intriguing possibility that short-hard GRBs are extragalactic events of a similar nature. Such events would be detectable with BATSE to a distance of \~50 Mpc, and their detection rate should be comparable to the actual BATSE detection rate of short-hard GRBs. The failure of our search, by contrast, suggests that such flares constitute less than 15% of the short-hard GRBs (<40% at 95% confidence). We discuss possible resolutions of this discrepancy.

preprint2001arXiv

Formation and Evaporation of Charged Black Holes

We investigate the dynamical formation and evaporation of a spherically symmetric charged black hole. We study the self-consistent one loop order semiclassical back-reaction problem. To this end the mass-evaporation is modeled by an expectation value of the stress-energy tensor of a neutral massless scalar field, while the charge is not radiated away. We observe the formation of an initially non extremal black hole which tends toward the extremal black hole $M=Q$, emitting Hawking radiation. If also the discharge due to the instability of vacuum to pair creation in strong electric fields occurs, then the black hole discharges and evaporates simultaneously and decays regularly until the scale where the semiclassical approximation breaks down. We calculate the rates of the mass and the charge loss and estimate the life-time of the decaying black holes.

preprint2000arXiv

Cosmic Censorship: The Role of Quantum Gravity

The cosmic censorship hypothesis introduced by Penrose thirty years ago is still one of the most important open questions in {\it classical} general relativity. In this essay we put forward the idea that cosmic censorship is intrinsically a {\it quantum gravity} phenomena. To that end we construct a gedanken experiment in which cosmic censorship is violated within the purely {\it classical} framework of general relativity. We prove, however, that {\it quantum} effects restore the validity of the conjecture. This suggests that classical general relativity is inconsistent and that cosmic censorship might be enforced only by a quantum theory of gravity.

preprint2000arXiv

Cosmic rays and TeV photons as probes of quantum properties of space-time

It has been recently observed that small violations of Lorentz invariance, of a type which may arise in quantum gravity, could explain both the observations of cosmic rays above the GZK cutoff and the observations of 20-TeV gamma rays from Markarian 501. We show here that different pictures of the short-distance structure of space-time would lead to different manifestations of Lorentz-invariance violation. Specifically, the deformation of Lorentz invariance needed to resolve these observational paradoxes can only arise within commutative short-distance pictures of space-time. In noncommutative space-times there is no anomalous effect, at least at leading order. Also exploiting the fact that arrival-time delays between high energy photons with different energies would arise in both the commutative and the noncommutative Lorentz-violation pictures, we describe an experimental programme, based on time-of-arrival analysis of high energy photons and searches of violations of GZK and TeV-photon limits, which could discriminate between alternative scenarios of Lorentz-invariance breakdown and could provide and unexpected window on the (quantum) nature of space-time at very short distances.

preprint2000arXiv

Effects of Pair Creation on Charged Gravitational Collapse

We investigate the effects of pair creation on the internal geometry of a black hole, which forms during the gravitational collapse of a charged massless scalar field. Classically, strong central Schwarzschild-like singularity forms, and a null, weak, mass-inflation singularity arises along the Cauchy horizon, in such a collapse. We consider here the discharge, due to pair creation, below the event horizon and its influence on the {\it dynamical formation} of the Cauchy horizon. Within the framework of a simple model we are able to trace numerically the collapse. We find that a part of the Cauchy horizon is replaced by the strong space-like central singularity. This fraction depends on the value of the critical electric field, $E_{\rm cr}$, for the pair creation.

preprint2000arXiv

Planck-scale deformation of Lorentz symmetry as a solution to the UHECR and the TeV-$γ$ paradoxes

One of the most puzzling current experimental physics paradoxes is the arrival on Earth of Ultra High Energy Cosmic Rays with energies above the GZK threshold. The recent observation of 20TeV photons from Mk 501 is another somewhat similar paradox. Several models have been proposed for the UHECR paradox. No solution has yet been proposed for the TeV-$γ$ paradox. Remarkably, the drastic assumption of a violation of ordinary Lorentz invariance would resolve both paradoxes. We present a formalism for the description of the type of Lorentz-invariance deformation (LID) that could be induced by non-trivial short-distance structure of space-time, and we show that this formalism is well suited for comparison of experimental data with LID predictions. We use the UHECR and TeV-$γ$ data, as well as bounds on time-of-flight differences between photons of different energies, to constrain the LID parameter space. A model with only two parameters, an energy scale and a dimensionless parameter characterizing the functional dependence on the energy scale, is shown to be sufficient to solve both the UHECR and the TeV-$γ$ threshold anomalies while satisfying the time-of-flight bounds. The allowed region of the two-parameter space is relatively small, but, remarkably, it fits perfectly the expectations of the quantum-gravity-motivated space-time models known to support such deformations of Lorentz invariance: integer value of the dimensionless parameter and characteristic energy scale constrained to a narrow interval in the neighborhood of the Planck scale.

preprint1999arXiv

The Inner Structure of Black Holes

We study the gravitational collapse of a self-gravitating charged scalar-field. Starting with a regular spacetime, we follow the evolution through the formation of an apparent horizon, a Cauchy horizon and a final central singularity. We find a null, weak, mass-inflation singularity along the Cauchy horizon, which is a precursor of a strong, spacelike singularity along the $r=0$ hypersurface. The inner black hole region is bounded (in the future) by singularities. This resembles the classical inner structure of a Schwarzschild black hole and it is remarkably different from the inner structure of a charged static Reissner-Nordström or a stationary rotating Kerr black holes.

preprint1998arXiv

Gamma-Ray Bursts and Related Phenomena

Gamma-ray bursts (GRBs) have puzzled astronomers since their accidental discovery in the sixties. The BATSE detector on the COMPTON-GRO satellite has been detecting one burst per day for the last six years. Its findings have revolutionized our ideas about the nature of these objects. They have shown that GRBs are at cosmological distances. This idea was accepted with difficulties at first. However, the recent discovery of an x-ray afterglow by the Italian/Dutch satellite BeppoSAX led to a detection of high red-shift absorption lines in the optical afterglow of GRB970508 and to a confirmation of its cosmological origin. The simplest and practically inevitable interpretation of these observations is that GRBs result from the conversion of the kinetic energy of ultra-relativistic particles flux to radiation in an optically thin region. The "inner engine" that accelerates the particles or generates the Poynting flux is hidden from direct observations. Recent studies suggest the ``internal-external'' model: internal shocks that take place within the relativistic flow produce the GRB while the subsequent interaction of the flow with the external medium produce the afterglow. The ``inner engine'' that produces the flow is, however, hidden from direct observations. We review this model with a specific emphasis on its implications to underground physics.

preprint1998arXiv

Images and Spectra From the Interior of a Relativistic Fireball

The detection of an afterglow, following a gamma-ray burst (GRB), can be described reasonably well by synchrotron emission from a relativistic spherical expanding blast wave, driven by an expanding fireball. We perform detailed calculations considering the emission from the whole region behind the shock front. We use the Blandford McKee self similar adiabatic solution to describe the fluid behind the shock. Using this detailed model, we derive expressions for the peak flux, and the peak frequency at a given observed time. These expressions provide important numerical corrections to previous, more simplified models. We calculate the observed light curve and spectra for several magnetic field models. We show that both the light curve and the spectra are flat near the peak. This rules out the interpretation of the optical peak of GRB970508 as the peak of the light curve, predicted by the existing fireball models. We calculate the observed image of a GRB afterglow. The observed image is bright near the edge and dimmer at the center, thus creating a ring. The contrast between the edge and the center is larger at high frequencies and the width of the ring is smaller.

preprint1997arXiv

A catalogue of the voids in the IRAS 1.2-Jy survey

Using the VOID FINDER algorithm we have compiled a catalogue of voids in the IRAS 1.2-Jy sample. The positions of the voids correspond well to underdense regions seen in the IRAS smoothed density map. However, since in our analysis no smoothing is used, all structures appear much sharper: walls are not smeared and the voids are not artificially reduced by them. Therefore the current method based on the galaxy point distribution is better suited to determine the diameter of voids in the galaxy distribution. We have identified 24 voids, covering more than 30% of the volume considered. By comparing the results with equivalent random catalogues we have determined that 12 voids are significant at a 0.95 confidence level, having an average diameter of 40+-6 h^{-1} Mpc. Our results serve not only for charting the cosmography of the nearby Universe, but also to give support to the results recently obtained with the SSRS2 sample, suggesting a void-filled Universe. Moreover, our results indicate that the voids detected have a similar scale, demonstrating that both optically and IRAS-selected galaxies delineate the same large-scale structures.

preprint1997arXiv

Can internal shocks produce the variability in GRBs?

We discuss the possibility that gamma-ray bursts result from internal shocks in an ultra-relativistic matter. Using a simple model we calculate the temporal structure and we estimate the efficiency of this process. In this model the ultra-relativistic matter flow is represented by a succession of shells with random values of the Lorentz factor. We calculate the shocks that take place between those shells and we estimate the resulting emission. Internal shocks can produce the highly variable temporal structure observed in most of the bursts provided that the source emitting the relativistic flow is highly variable. The observed peaks are in almost one to one correlation to the activity of the emitting source. A large fraction of the kinetic energy is converted to radiation. The most efficient case is when an inner engine produces shells with comparable energy but very different Lorentz factors. It also gives the most preferable temporal structure.

preprint1997arXiv

On Gravitational Repulsion

The concepts of negative gravitational mass and gravitational repulsion are alien to general relativity. Still, we show here that small negative fluctuations - small dimples in the primordial density field - that act as if they have an effective negative gravitational mass, play a dominant role in shaping our Universe. These initially tiny perturbations repel matter surrounding them, expand and grow to become voids in the galaxy distribution. These voids - regions with a diameter of $40 \h$ Mpc which are almost devoid of galaxies - are the largest object in the Universe.

preprint1997arXiv

Variability in GRBs - A Clue

We show that external shocks cannot produce a variable GRB, unless they are produced by an extremely narrow jets (angular opening of < ~10^{-4}) or if only a small fraction of the shell emits the radiation and the process is very inefficient. Internal shocks can produce the observed complex temporal structure provided that the source itself is variable. In this case, the observed temporal structure reflects the activity of the ``inner engine'' that drives the bursts. This sets direct constraints on it.

preprint1996arXiv

Cosmological GRBs: Internal vs. External Shocks

Internal (IS) or external (ES) shocks in relativistic expanding shells are currently the best known mechanism for producing GRBs. We calculate the hydrodynamic conditions and the cooling processes in internal shocks, that occur when one layer overtakes another. We compare these conditions with the conditions in external shocks that occur when the shell encounters the ISM. We find that the complicated temporal structure observed in GRBs could be easily explained in IS scenario. Unlike ES, the observed temporal structure simply reflects the temporal behavior of the internal engine that produces the shell. We also find that opacity and in particular the opacity for pair production poses strong constraints on the parameter space and consequently bursts with very narrow peaks are expected to be optically thick and would not contain GeV photons. Finally we find that as in ES synchrotron is the most likely radiation process.

preprint1996arXiv

Critical Behaviour and Universality in Gravitational Collapse of a Charged Scalar Field

We summarize results from a study of spherically symmetric collapse of a {\it charged} (complex) massless scalar-field \cite{Hod}. We present an analytic argument which conjecture the generalization of the mass-scaling relation and echoing phenomena, originally discovered by Choptuik, for the {\it charged} case. Furthermore, we study the behaviour of the self-similar critical solution under {\it external} perturbations -- addition of a cosmological constant $Λ$ and a charge-conjugation $e$. Finally, we study the scaling-relation of the black-hole charge. Using an analytic argument we conjecture that black-holes of infinitesimal mass are neutral or obey the relation $Q_{BH} \ll M_{BH}$. We verify our predictions with numerical results.

preprint1996arXiv

The X-Ray Background as a Probe of Density Fluctuations at High Redshift

The X-Ray Background (XRB) probes structure on scales intermediate between those explored by local galaxy redshift surveys and by the COBE Microwave Background measurements. We predict the large scale angular fluctuations in the XRB, expressed in terms of spherical harmonics for a range of assumed power-spectra and evolution scenarios. The dipole is due to large scale structure as well as to the observer's motion (the Compton-Getting effect). For a typical observer the two effects turn out to be comparable in amplitude. The coupling of the two effects makes it difficult to use the XRB for independent confirmation of the CMB dipole being due to the observer's motion. The large scale structure dipole (rms per component) relative to the monopole is in the range $a_{1m}/a_{00} \sim (0.5-9.0) \times 10^{-3} $. The spread is mainly due to the assumed redshift evolution scenarios of the X-ray volume emissivity $ρ_x(z)$. The dipole's prediction is consistent with a measured dipole in the HEAO1 XRB map. Typically, the harmonic spectrum drops with $l$ like $a_{lm} \sim l^{-0.4}$. This behaviour allows us to discriminate a true clustering signal against the flux shot noise, which is constant with $l$, and may dominate the signal unless bright resolved sources are removed from the XRB map. We also show that Sachs-Wolfe and Doppler (due to the motion of the sources) effects in the XRB are negligible. Although our analysis focuses on the XRB, the formalism is general and can be easily applied to other cosmological backgrounds.

preprint1994arXiv

The Internal Geometry of an Evaporating Black Hole

We present a semi-classical model for the formation and evaporation of a four dimensional black hole. We solve the equations numerically and obtain solutions describing the entire the space-time geometry from the collapse to the end of the evaporation. The solutions satisfy the evaporation law: $\dot M \propto -M^{-2}$ which confirms dynamically that black holes do evaporate thermally. We find that the evaporation process is in fact the shrinking of a throat that connects a macroscopic interior ``universe" to the asymptotically flat exterior. It ends either by pinching off the throat leaving a closed universe and a Minkowskian exterior or by freezing up when the throat's radius approaches a Planck size. In either case the macroscopic inner universe is the region where the information lost during the evaporation process is hidden.

preprint1993arXiv

Do Gamma-Ray Burst Sources Repeat

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

preprint1993arXiv

Hydrodynamics of Relativistic Fireballs

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

Tsvi Piran

What is connected

Connect this record

See the researcher in context

Building this map preview

117 published item(s)

The Advanced X-ray Imaging Satellite (AXIS) Community Science Book

Analytic Modeling of Synchrotron-Self-Compton Spectra: Application to GRB 190114C

What powers the radio emission in TDE AT2019dsg: a long-lived jet or the disruption itself?

High-energy emission from tidal disruption events in active galactic nuclei

Limits on mass outflow from optical tidal disruption events

Generalized compactness limit from an arbitrary viewing angle

On short GRBs similar to GRB 170817A detected by Fermi GBM

Radiative Emission Mechanisms of Tidal Disruption Events

Shock within a shock: revisiting the radio flares of NS merger ejecta and GRB-supernovae

The Origin of Binary Black Holes Mergers

The origin of hotspots around Sgr A*: Orbital or pattern motion?

Tidal disruption events in active galactic nuclei

Tidal radii of main sequence stars -- I. Physical tidal radius, semi-analytic model and their implications

Tidal radii of main sequence stars -- II. Simulation methodology and the character of full tidal disruptions

Tidal radii of main sequence stars -- III. Partial disruptions

Magnetorotational instability in eccentric disks

A Revised Analysis of Gamma Ray Bursts' prompt efficiencies

Are Ultra-Long Gamma-Ray Bursts Caused by Blue Supergiant Collapsars, Newborn Magnetars, or White Dwarf Tidal Disruption Events?

ASASSN-14li: A Model Tidal Disruption Event

Beaming of particles and synchrotron radiation in relativistic magnetic reconnection

Cosmic explosions, life in the Universe and the Cosmological Constant

Formation of Double Neutron Star systems as implied by observations

Natal Kicks and Time Delays in Merging Neutron Star Binaries - Implications for r-process nucleosynthesis in Ultra Faint Dwarfs and in the Milky Way

On the maximal efficiency of the collisional Penrose process

Particle acceleration, magnetization and radiation in relativistic shocks

r-process Production Sites as inferred from Eu Abundances in Dwarf Galaxies

Radio Counterparts of Compact Binary Mergers detectable in Gravitational Waves: A Simulation for an Optimized Survey

Testing the Magnetar Model via Late Time Radio Observations of Two Macronova Candidates

The B/C and sub-Iron/Iron Cosmic ray ratios - further evidence in favor of the spiral arm diffusion model

The dynamics of radiation driven, optically thick winds

The effect of cooling on particle trajectories and acceleration in relativistic magnetic reconnection

The Macronova in GRB 050709 and the GRB/macronova connection

The Observable Signatures of GRB Cocoons

"Circularization" vs. Accretion -- What Powers Tidal Disruption Events?

A possible Macronova in the late afterglow of the `long-short' burst GRB 060614

Energies of GRB blast waves and prompt efficiencies as implied by modeling of X-ray and GeV afterglows

Gamma-Ray Bursts as Sources of Strong Magnetic Fields

General Relativistic Hydrodynamic Simulation of Accretion Flow from a Stellar Tidal Disruption

Jet and disk luminosities in tidal disruption events

Linear and circular polarization in GRB afterglows

Linear and Circular polarization in ultra-relativistic synchrotron sources - implications to GRB afterglows

Mass ejection from neutron star mergers: different components and expected radio signals

Radio Flares of Compact Binary Mergers: the Effect of Non-Trivial Outflow Geometry

Radioactive decay products in neutron star merger ejecta: heating efficiency and $γ$-ray emission

Shedding light on the prompt high efficiency paradox - self consistent modeling of GRB afterglows

Spectral signatures of compact sources in the inverse Compton catastrophe limit

Super-Penrose collisions are inefficient - a Comment on: Black hole fireworks: ultra-high-energy debris from super-Penrose collisions

The afterglow of a relativistic shock breakout and low-luminosity GRBs

The fate of supernova remnants near quiescent supermassive black holes

The rate, luminosity function and time delay of non-Collapsar short GRBs

Implications of GRB 130603B and its macronova for r-process nucleosynthesis

On the Composition of GRBs' Collapsar Jets

On the role of GRBs on life extinction in the Universe

The dynamics of a highly magnetized jet propagating inside a star

The emission mechanism in magnetically dominated GRB outflows

The long-term evolution of neutron star merger remnants - II. Radioactively powered transients

The Nature of ULX Source M101 X-1: Optically Thick Outflow from A Stellar Mass Black Hole

Tidal disruption jets as the source of Ultra-High Energy Cosmic Rays

Binary Pulsar J0737-3039: Evidence for a new core collapse and neutron star formation mechanism

Constraints on the Synchrotron Emission Mechanism in GRBs

Gravitational Wave Memory from Gamma Ray Bursts' Jets

On the origin of the radio emission of Sw 1644+57

Possible relevance of quantum spacetime for neutrino-telescope data analyses

Radius constraints and minimal equipartition energy of relativistically moving synchrotron sources

Reconciling the GRB rate and star formation histories

The Discovery of a New Instability in a Hyperaccretion Flow and its Implication for Gamma-ray Bursts

Choptuik's Critical Phenomenon in Einstein-Gauss-Bonnet Gravity

Collisional Penrose process near the horizon of extreme Kerr black holes

Implications of the UHECRs penetration depth measurements

Jet Luminosity From Neutrino-Dominated Accretion Flows in Gamma-Ray Bursts

Jets from Tidal Disruptions of Stars by Black Holes

On the lateral expansion of GRB jets

On thermalization in gamma-ray burst jets and the peak energies of photospheric spectra

Short vs Long and Collapsars vs. non-Collapsar: a quantitative classification of GRBs