BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Ehud Nakar

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

ResearcherUnclaimed source record
astro-ph.HEastro-ph.SRastro-phastro-ph.COgr-qcastro-ph.EPastro-ph.IMphysics.plasm-phastro-ph.GAphysics.flu-dyn

Catalog footprint

What is connected

61works
10topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

61 published item(s)

preprint2022arXiv

Observational signatures of stellar explosions driven by relativistic jets

The role of relativistic jets in unbinding the stellar envelope during a supernova (SN) associated with a gamma-ray burst (GRB) is unclear. To study that, we explore observational signatures of stellar explosions that are driven by jets. We focus on the final velocity distribution of the outflow in such explosions and compare its observational imprints to SN/GRB data. We find that jet driven explosions produce an outflow with a flat distribution of energy per logarithmic scale of proper velocity. The flat distribution seems to be universal as it is independent of the jet and the progenitor properties that we explored. The velocity range of the flat distribution for typical GRB parameters is $γβ\approx 0.03-3$, where $γ$ is the outflow Lorentz factor and $β$ is its dimensionless velocity. A flat distribution is seen also for collimated choked jets where the highest outflow velocity decreases with the depth at which the jet is choked. Comparison to observations of SN/GRBs rules out jets as the sole explosion source in these events. Instead, in SN/GRB the collapsing star must deposit its energy into two channels - a quasi-spherical (or wide angle) channel and a narrowly collimated one. The former carries most of the energy and is responsible for the SN sub-relativistic ejecta while the latter carries 0.01-0.1 of the total outflow energy and is the source of the GRB. Intriguingly, the same two channels, with a similar energy ratio, were seen in the binary neutron star merger GW170817, suggesting that similar engines are at work in both phenomena.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The role of plasma instabilities in relativistic radiation mediated shocks: stability analysis and particle-in-cell simulations

Relativistic radiation mediated shocks (RRMS) likely form in prodigious cosmic explosions. The structure and emission of such shocks is regulated by copious production of electron-positron pairs inside the shock transition layer. It has been pointed out recently that substantial abundance of positrons inside the shock leads to a velocity separation of the different plasma constituents, which is expected to induce a rapid growth of plasma instabilities. In this paper, we study the hierarchy of plasma microinstabilities growing in an electron-ion plasma loaded with pairs and subject to a radiation force. Linear stability analysis indicates that such a system is unstable to the growth of various plasma modes which ultimately become dominated by a current filamentation instability driven by the relative drift between the ions and the pairs. These results are validated by particle-in-cell simulations that further probe the nonlinear regime of the instabilities, and the pair-ion coupling in the microturbulent electromagnetic field. Based on this analysis, we derive a reduced transport equation for the particles via pitch angle scattering in the microturbulence and demonstrate that it can couple the different species and lead to nonadiabatic compression via a Joule-like heating. The heating of the pairs and, conceivably, the formation of nonthermal distributions, arising from the microturbulence, can affect the observed shock breakout signal in ways unaccounted for by current single-fluid models.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A transient radio source consistent with a merger-triggered core collapse supernova

A core-collapse supernova occurs when exothermic fusion ceases in the core of a massive star, typically due to exhaustion of nuclear fuel. Theory predicts that fusion could be interrupted earlier, by merging of the star with a compact binary companion. We report a luminous radio transient, VT J121001+495647, found in the Very Large Array Sky Survey. The radio emission is consistent with supernova ejecta colliding with a dense shell of material, potentially ejected by binary interaction in the centuries prior to explosion. We associate the supernova with an archival X-ray transient, which implies a relativistic jet was launched during the explosion. The combination of an early relativistic jet and late-time dense interaction is consistent with expectations for a merger-driven explosion.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Physics of radiation mediated shocks and its applications to GRBs, supernovae, and neutron star mergers

The first electromagnetic signal observed in different types of cosmic explosions is released upon emergence of a shock created in the explosion from the opaque envelope enshrouding the central source. Notable examples are the early emission from various types of supernovae and low luminosity GRBs, the prompt photospheric emission in long GRBs, and the gamma-ray emission that accompanied the gravitational wave signal in neutron star mergers. In all of these examples, the shock driven by the explosion is mediated by the radiation trapped inside it, and its velocity and structure, that depend on environmental conditions, dictate the characteristics of the observed electromagnetic emission at early times, and potentially also their neutrino emission. Much efforts have been devoted in recent years to develop a detailed theory of radiation mediated shocks in an attempt to predict the properties of the early emission in the aforementioned systems. These efforts are timely in view of the anticipated detection rate of shock breakout candidates by upcoming transient factories, and the potential detection of a gamma-ray flash from shock breakout in neutron star mergers like GW170817. This review aims at providing a comprehensive overview of the theory and applications of radiation mediated shocks, starting from basic principles. The classification of shock solutions, which are governed by the conditions prevailing in each class of objects, and the methods used to solve the shock equations in different regimes will be described, with particular emphasis on the observational diagnostics. The applications to supernovae, low-luminosity GRBs, long GRBs, neutron star mergers, and neutrino emission will be highlighted.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Shock breakouts from red supergiants: analytical and numerical predictions

Shock breakout (SBO) signal is the first signature of the supernova explosion apart from gravitational waves and neutrinos. Observational properties of SBO, such as bolometric luminosity and colour temperature, connect to the supernova progenitor and explosion parameters. Detecting SBO or SBO-cooling will constrain the progenitor and explosion models of collapsing stars. In the light of recently launched eROSITA telescope, the rate for detection of SBO is a few events during a year. In the current study, we examine the analytic formulae derived by Shussman et al. (2016). We use four red supergiant models from their study, while running explosions with the radiation hydrodynamics code STELLA. We conclude that there is a good agreement between analytic and numerical approaches for bolometric luminosity and colour temperature during SBO. The analytic formulae for the SBO signal based on the global supernova parameters can be used instead of running time-consuming numerical simulations. We define spectral range where analytic formulae for the SBO spectra are valid. We provide improved analytical expression for the SBO spectral energy distribution. We confirm dependence of colour temperature on radius derived by analytical studies and suggest to use early time observations to confine the progenitor radius. Additionally we show the prediction for the SBO signal from red supergiants as seen by eROSITA instrument.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Spherical shocks in a steep density gradient of expanding media

We study the propagation of a Newtonian shock in a spherically symmetric, homologously expanding ejecta. We focus on media with a steep power-law density profile of the form $ρ\propto t^{-3}v^{-α}$, with $α>5$, where $v$ is the velocity of the expanding medium and $t$ is time. Such profiles are expected in the leading edge of supernovae ejecta and sub-relativistic outflows from binary neutron star mergers. We find that such shocks always accelerate in the observer frame and lose causal contact with the bulk of the driver gas, owing to the steep density profile. However, the prolonged shock evolution exhibits two distinct pathways: In one, the shock strength diminishes with time until the shock eventually dies out. In the other, the shock strength steadily increases, and the solution approaches the self-similar solution of a shock is a static medium. By mapping the parameter space of shock solutions, we find that the evolutionary pathways are dictated by $α$ and by the initial ratio between the shock velocity and the local upstream velocity. We find that for $α<ω_c$ ($ω_c \approx 8$), the shock always decays, and that for $α>ω_c$ the shock may decay or grow stronger depending on the initial value of the velocity ratio. These two branches bifurcate from a self-similar solution derived analytically for a constant velocity ratio. We analyze properties of the solutions that may have an impact on the observational signatures of such systems, and assess the conditions required for decaying shocks to break out from a finite medium.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The structure of weakly-magnetized $ γ$-ray burst jets

The interaction of gamma-ray burst (GRB) jets with the dense media into which they are launched promote the growth of local hydrodynamic instabilities along the jet boundary. In a companion paper we study the evolution of hydrodynamic (unmagnetized) jets, finding that mixing of jet-cocoon material gives rise to an interface layer, termed jet-cocoon interface (JCI), which contains a significant fraction of the system energy. We find that the angular structure of the jet + JCI, when they reach the homologous phase, can be approximated by a flat core (the jet) + a power-law function (the JCI) with indices that depend on the degree of mixing. In this paper we examine the effect of subdominant toroidal magnetic fields on the jet evolution and morphology. We find that weak fields can stabilize the jet against local instabilities. The suppression of the mixing diminishes the JCI and thus reshapes the jet's post-breakout structure. Nevertheless, the overall shape of the outflow can still be approximated by a flat core + a power-law function, although the JCI power-law decay is steeper. The effect of weak fields is more prominent in long GRB jets, where the mixing in hydrodynamic jets is stronger. In short GRB jets there is small mixing in both weakly magnetized and unmagnetized jets. This result influences the expected jet emission which is governed by the jet's morphology. Therefore, prompt and afterglow observations in long GRBs may be used as probes for the magnetic nature at the base of the jets.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

GROWTH on S190814bv: Deep Synoptic Limits on the Optical/Near-Infrared Counterpart to a Neutron Star-Black Hole Merger

On 2019 August 14, the Advanced LIGO and Virgo interferometers detected the high-significance gravitational wave (GW) signal S190814bv. The GW data indicated that the event resulted from a neutron star--black hole (NSBH) merger, or potentially a low-mass binary black hole merger. Due to the low false alarm rate and the precise localization (23 deg$^2$ at 90\%), S190814bv presented the community with the best opportunity yet to directly observe an optical/near-infrared counterpart to a NSBH merger. To search for potential counterparts, the GROWTH collaboration performed real-time image subtraction on 6 nights of public Dark Energy Camera (DECam) images acquired in the three weeks following the merger, covering $>$98\% of the localization probability. Using a worldwide network of follow-up facilities, we systematically undertook spectroscopy and imaging of optical counterpart candidates. Combining these data with a photometric redshift catalog, we ruled out each candidate as the counterpart to S190814bv and we placed deep, uniform limits on the optical emission associated with S190814bv. For the nearest consistent GW distance, radiative transfer simulations of NSBH mergers constrain the ejecta mass of S190814bv to be $M_\mathrm{ej} < 0.04$~$M_{\odot}$ at polar viewing angles, or $M_\mathrm{ej} < 0.03$~$M_{\odot}$ if the opacity is $κ< 2$~cm$^2$g$^{-1}$. Assuming a tidal deformability for the neutron star at the high end of the range compatible with GW170817 results, our limits would constrain the BH spin component aligned with the orbital momentum to be $ χ< 0.7$ for mass ratios $Q < 6$, with weaker constraints for more compact neutron stars. We publicly release the photometry from this campaign at http://www.astro.caltech.edu/~danny/static/s190814bv.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Radioactive heating rate of r-process elements and macronova light curve

We study the heating rate of r-process nuclei and thermalization of decay products in neutron star merger ejecta and macronova (kilonova) light curves. Thermalization of charged decay products, i.e., electrons, $α$-particles, and fission fragments is calculated according to their injection energy. The $γ$-ray thermalization processes are also properly calculated by taking the $γ$-ray spectrum of each decay into account. We show that the $β$-decay heating rate at later times approaches a power-law decline as $\propto t^{-2.8}$, which agrees with the result of Waxman et al. (2019). We present a new analytic model to calculate macronova light curves, in which the density structure of the ejecta is accounted for. We demonstrate that the observed bolometric light curve and temperature evolution of the macronova associated with GW170817 are reproduced well by the $β$-decay heating rate with the solar r-process abundance pattern. We interpret the break in the observed bolometric light curve around a week as a result of the diffusion wave crossing a significant part of the ejecta rather than a thermalization break. We also show that the time-weighted integral of the bolometric light curve (Katz integral) is useful to provide an estimate of the total r-process mass from the observed data, which is independent of the highly uncertain radiative transfer. For the macronova in GW170817, the ejecta mass is robustly estimated as $\approx 0.05M_{\odot}$ for $A_{\rm min}\leq 72$ and $85\leq A_{\rm min}\leq 130$ with the solar r-process abundance pattern. The code for computation of the heating rate and light curve for given initial nuclear abundances is publicly available.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The importance of 56Ni in shaping the light curves of type II supernovae

What intrinsic properties shape the light curves of Type II supernovae (SNe)? To address this question we derive observational measures that are robust (i.e., insensitive to detailed radiative transfer) and constrain the contribution from $^{56}$Ni, as well as a combination of the envelope mass, progenitor radius, and explosion energy. By applying our methods to a sample of type II SNe from the literature we find that $^{56}$Ni contribution is often significant. In our sample its contribution to the time weighted integrated luminosity during the photospheric phase ranges between 8% and 72% with a typical value of 30%. We find that the $^{56}$Ni relative contribution is anti-correlated with the luminosity decline rate. When added to other clues, this in turn suggests that the flat plateaus often observed in type II SNe are not a generic feature of the cooling envelope emission, and that without $^{56}$Ni many of the SNe that are classified as II-P would have shown a decline rate that is steeper by up to 1 mag/100 d. Nevertheless, we find that the cooling envelope emission, and not $^{56}$Ni contribution, is the main driver behind the observed range of decline rates. Furthermore, contrary to previous suggestions, our findings indicate that fast decline rates are not driven by lower envelope masses. We therefore suggest that the difference in observed decline rates is mainly a result of different density profiles of the progenitors.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Type II supernovae Early Light Curves

Observations of type II supernova early light, from breakout until recombination, can be used to constrain the explosion energy and progenitor properties. Currently available for this purpose are purely analytic models, which are accurate only to within an order of magnitude, and detailed numerical simulations, which are more accurate but are applied to any event separately. In this paper we derive an analytic model that is calibrated by numerical simulations. This model is much more accurate than previous analytic models, yet it is as simple to use. To derive the model we analyze simulated light curves from numerical explosion of $124$ red supergiant progenitors, calculated using the stellar evolution code MESA. We find that although the structure of the progenitors we consider varies, the resulting light curves can be described rather well based only on the explosion energy, ejecta mass and progenitor radius. Our calibrated analytic model, which is based on these three parameters, reproduces the bolometric luminosity within $25\%-35\%$ accuracy and the observed temperature within $15\%$ accuracy (compared to previous analytic models which are indeed found to be accurate only to within an order of magnitude). We also consider deviations of the early time spectrum from blackbody, and find that the Rayleigh-Jeans regime is slightly shallower (roughly $L_ν\propto ν^{1.4}$). This modified spectrum affects the optical/near-UV light curve mostly during the first day when the typical observed temperature is $\gg 10^4 ~^\circ$K. We use our results to study the optical and near-UV early light curves from first light until recombination and briefly discuss what can be learned from current and future observations. Light curves generated using our calibrated model can be downloaded at http://www.astro.tau.ac.il/~tomersh/.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Type II supernovae progenitor and ejecta properties from the total emitted light, ET

It was recently shown that the bolometric light curves of type II supernovae (SNe) allow an accurate and robust measurement of the product of the radiation energy in the ejecta, $E_r$, and the time since the explosion, $t$, at early phases ($t\lesssim 10d$) of the homologous expansion. This observable, denoted here $ET \equiv E_rt$ is constant during that time and depends only on the progenitor structure and explosion energy. We use a 1D hydrodynamic code to find $ET$ of simulated explosions of 145 red supergiant progenitors obtained using the stellar evolution code MESA, and relate this observable to the properties of the progenitor and the explosion energy. We show that $ET$ probes only the properties of the envelope (velocity, mass and initial structure), similarly to other observables that rely on the photospheric phase emission. Nevertheless, for explosions where the envelope dominates the ejected mass, $M_{env}/M_{ej} \gtrsim 0.6$, $ET$ is directly related to the explosion energy $E_{exp}$ and ejected mass $M_{ej}$ through the relation $ET \approx 0.15 E_{exp}^{1/2} R_* M_{ej}^{1/2}$, where $R_*$ is the progenitor radius, to an accuracy better than $30\%$. We also provide relations between $ET$ and the envelope properties that are accurate (to within 20\%) for all the progenitors in our sample, including those that lost most of their envelope. We show that when the envelope velocity can be reasonably measured by line shifts in observed spectra, the envelope is directly constrained from the bolometric light curve (independent of $E_{exp}$). We use that to compare observations of 11 SNe with measured $ET$ and envelope velocity to our sample of numerical progenitors. This comparison suggests that many SNe progenitors have radii that are $\lesssim 500~R_\odot$. In the framework of our simulations this indicates, most likely, a rather high value of the mixing length parameter.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A unified picture for low-luminosity and long gamma-ray bursts based on the extended progenitor of llgrb 060218/SN 2006aj

The relation between long gamma-ray bursts (LGRBs) and low-luminosity GRBs (llgrbs) is a long standing puzzle -- on the one hand their high energy emission properties are fundamentally different, implying a different gamma-ray source, yet both are associated with similar supernovae of the same peculiar type (broad-line Ic), pointing at a similar progenitor and a similar explosion mechanism. Here we analyze the multi-wavelength data of the particularly well-observed SN 2006aj, associated with llgrb 060218, finding that its progenitor star is sheathed in an extended ($>100R_\odot$), low-mass ($\sim 0.01M_\odot$) envelope. This progenitor structure implies that the gamma-ray emission in this llgrb is generated by a mildly relativistic shock breakout. It also suggests a unified picture for llgrbs and LGRBs, where the key difference is the existence of an extended low-mass envelope in llgrbs and its absence in LGRBs. The same engine, which launches a relativistic jet, can drive the two explosions, but, while in LGRBs the ultra-relativistic jet emerges from the bare progenitor star and produces the observed gamma-rays, in llgrbs the extended envelope smothers the jet and prevents the generation of a large gamma-ray luminosity. Instead, the jet deposits all its energy in the envelope, driving a mildly relativistic shock that upon breakout produces a llgrb. In addition for giving a unified view of the two phenomena, this model provides a natural explanation to many observed properties of llgrbs. It also implies that llgrbs are a viable source of the observed extra-galactic diffuse neutrino flux and that they are promising sources for future gravitational wave detectors.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

X-ray-powered macronovae

A macronova (or kilonova) was observed as an infrared excess several days after short gamma-ray burst, GRB 130603B. Although the $r$-process radioactivity is widely discussed as an energy source, it requires huge mass of ejecta from a neutron star (NS) binary merger. We propose a new model that the X-ray excess gives rise to the simultaneously observed infrared excess via thermal re-emission and explore what constraints this would place on the mass and velocity of the ejecta. This X-ray-powered model explains both the X-ray and infrared excesses with a single energy source by the central engine like a black hole, and allows for broader parameter region, in particular smaller ejecta mass $\sim10^{-3}-10^{-2}M_{\odot}$ with iron mixed as suggested by general relativistic simulations for typical NS-NS mergers, than the previous models. We also discuss the other macronova candidates in GRB 060614 and GRB 080503, and implications for the search of electromagnetic counterparts to gravitational waves.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Constraints on Shallow 56Ni from the Early Lightcurves of Type Ia Supernovae

Ongoing transient surveys are presenting an unprecedented account of the rising lightcurves of Type Ia supernovae (SNe Ia). This early emission probes the shallowest layers of the exploding white dwarf, which can provide constraints on the progenitor star and the properties of the explosive burning. We use semi-analytic models of radioactively-powered rising lightcurves to analyze these observations. As we have summarized in previous work, the main limiting factor in determining the surface distribution of 56Ni is the lack of an unambiguously identified time of explosion, as would be provided by detection of shock breakout or shock-heated cooling. Without this the SN may in principle exhibit a "dark phase" for a few hours to days, where the only emission is from shock-heated cooling that is too dim to be detected. We show that by assuming a theoretically-motivated time-dependent velocity evolution, the explosion time can be better constrained, albeit with potential systematic uncertainties. This technique is used to infer the surface 56Ni distributions of three recent SNe Ia that were caught especially early in their rise. In all three we find fairly similar 56Ni distributions. Observations of SN 2011fe and SN 2012cg probe shallower depths than SN 2009ig, and in these two cases 56Ni is present merely ~0.01Msun from the WDs' surfaces. The uncertainty in this result is up to an order of magnitude given the difficulty of precisely constraining the explosion time. We also use our conclusions about the explosion times to reassess radius constraints for the progenitor of SN 2011fe, as well as discuss the roughly t^2 power law that is inferred for many observed rising lightcurves.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Photometric and Spectroscopic Properties of Type II-P Supernovae

We study a sample of 23 Type II Plateau supernovae (SNe II-P), all observed with the same set of instruments. Analysis of their photometric evolution confirms that their typical plateau duration is 100 days with little scatter, showing a tendency to get shorter for more energetic SNe. The rise time from explosion to plateau does not seem to correlate with luminosity. We analyze their spectra, measuring typical ejecta velocities, and confirm that they follow a well behaved power-law decline. We find indications of high-velocity material in the spectra of six of our SNe. We test different dust extinction correction methods by asking the following -- does the uniformity of the sample increase after the application of a given method? A reasonably behaved underlying distribution should become tighter after correction. No method we tested made a significant improvement.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Recombination Effects on Supernovae Light-Curves

Supernovae of type IIP are marked by the long plateau seen in their optical light curves. The plateau is believed to be the result of a recombination wave that propagates through the outflowing massive hydrogen envelope. Here, we analytically investigate the transition from a fully ionized envelope to a partially recombined one and its effects on the SN light curve. The motivation is to establish the underlying processes which dominate the evolution at late times when recombination takes place in the envelope, yet early enough so that $^{56}$Ni decay is a negligible source of energy. We assume a simple, yet adequate, hydrodynamic profile of the envelope and study the mechanisms which dominate the energy emission and the observed temperature. We consider the diffusion of photons through the envelope while analyzing the ionization fraction and the coupling between radiation and gas. We find that once recombination starts, the observed temperature decreases slowly in time. However, in a typical red supergiant (RSG) explosion, the recombination wave does not affect the bolometric luminosity immediately. Only at later times, the cooling wave may reach layers that are deep enough to affect the luminosity. We find that the plateau is not a generic result of a recombination process in expanding gas. Instead it depends on the density profile of the parts of the envelope which undergo recombination. Our results are useful to investigate the light curves of RSG explosions. We show the resulting light curves of two examples of RSG explosions according to our model and discuss their compatibility with observations. In addition, we improve the analytical relations between the plateau luminosity and plateau duration to the properties of the pre-explosion progenitor (Arnett 1980; Popov 1993).

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

SN 2008D: A Wolf-Rayet explosion through a thick wind

Supernova (SN) 2008D/XRT 080109 is considered to be the only direct detection of a shock breakout from a regular SN to date. While a breakout interpretation was favored by several papers, inconsistencies remain between the observations and current SN shock breakout theory. Most notably, the duration of the luminous X-ray pulse is considerably longer than expected for a spherical breakout through the surface of a type Ibc SN progenitor, and the X-ray radiation features, mainly its flat spectrum and its luminosity evolution, are enigmatic. We apply a recently developed theoretical model for the observed radiation from a Wolf-Rayet SN exploding through a thick wind and show that it naturally explains all the observed features of SN 2008D X-ray emission, including the energetics, the spectrum and the detailed luminosity evolution. We find that the inferred progenitor and SN parameters are typical for an exploding Wolf-Rayet. A comparison of the wind density found at the breakout radius to the density at much larger radii, as inferred by late radio observations, suggests an enhanced mass loss rate taking effect about ten days or less prior to the SN explosion. This finding joins accumulating evidence for a possible late phase in the stellar evolution of massive stars, involving vigorous mass loss a short time before the SN explosion.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Spectrum and light curve of a supernova shock breakout through a thick Wolf-Rayet wind

Wolf-Rayet stars are known to eject winds. Thus, when a Wolf-Rayet star explodes as a supernova, a fast, $>30,000$ km/s, shock is expected to be driven through a wind. We study the signal expected from a fast supernova shock propagating through an optically thick wind, and find that the electrons behind the shock driven into the wind are cooled efficiently, by inverse Compton over soft photons that were deposited by the radiation mediated shock that crossed the star. Therefore, the bolometric luminosity is comparable to the kinetic energy flux through the shock, and the spectrum is found to be a power-law, which slope and frequency range depend on the number flux of soft photons available for cooling. Wolf-Rayet supernovae that explode through a thick wind have a high flux of soft photons, producing a flat spectrum, $νF_ν=Const$, in the X-ray range $0.1\lesssim T\lesssim 50$ keV. As the shock expands into an optically thin wind, the soft photons are no longer able to cool the shock that plows through the wind, and the bulk of the emission takes the form of a standard core-collapse supernova (without a wind). However, a small fraction of the soft photons is upscattered by the shocked wind and produces a transient unique X-ray signature.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Supernovae with Two Peaks in the Optical Light Curve and the Signature of Progenitors with Low-mass Extended Envelopes

Early observations of supernova light curves are powerful tools for shedding light on the pre-explosion structures of their progenitors and their mass-loss histories just prior to explosion. Some core-collapse supernovae that are detected during the first days after the explosion prominently show two peaks in the optical bands, including the $R$ and $I$ bands, where the first peak appears to be powered by the cooling of shocked surface material and the second peak is clearly powered by radioactive decay. Such light curves have been explored in detail theoretically for SN 1993J and 2011dh, where it was found that they may be explained by progenitors with extended, low-mass envelopes. Here we generalize these results. We explore, first, whether any double-peaked light curve of this type can be generated by a progenitor with a "standard" density profile, such as a red supergiant or a Wolf-Rayet star. We show that a standard progenitor (1) cannot produce a double-peaked light curve in the $R$ and $I$ bands, and (2) cannot exhibit a fast drop in the bolometric luminosity as is seen after the first peak. We then explore the signature of a progenitor with a compact core surrounded by extended, low-mass material. This may be a hydrostatic low-mass envelope or material ejected just prior to the explosion. We show that it naturally produces both of these features. We use this result to provide simple formulae to estimate (1) the mass of the extended material from the time of the first peak, (2) the extended material radius from the luminosity of the first peak, and (3) an upper limit on the core radius from the luminosity minimum between the two peaks.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A Wolf in Sheepskin: Extraordinary Supernova 2012au Veiled Behind Ordinary Radio Emission

We present extensive radio and X-ray observations of SN\, 2012au, the energetic radio luminous supernova of type Ib that may be a link between subsets of hydrogen-poor superluminous and normal core-collapse supernovae. The observations closely follow models of synchrotron emission from shock heated circum-burst medium that has a wind density profile ($ρ\propto r^{-2}$). We infer a sub-relativistic velocity for the shock wave $v \approx 0.2\,c$ and a radius of $r \approx 1.4 \times 10^{16} \rm cm$ at 25 days after the estimated date of explosion. For a constant wind velocity of 1000 km/s we determine the constant mass loss rate of the progenitor to be $\dot{M} = 3.6 \times 10^{-6} \rm M_{\odot} yr^{-1}$, consistent with the estimates from X-ray observations. We estimate the total internal energy of the radio emitting material to be $E \approx 10^{47} \rm erg$, which is intermediate to SN\,1998bw and SN\,2002ap. Evolution of the radio light curves of SN\,2012au is consistent with interaction with a smoothly distributed circum-burst medium and absence of stellar shells ejected from previous outbursts out to $r \approx 10^{17} \rm cm$ from the supernova site. Based on this we conclude that the evolution of the SN\,2012au progenitor star was relatively quiet during the final years preceding explosion. We find that the bright radio emission from SN2012au was not dissimilar from other core collapse supernovae despite it's extraordinary optical properties. We speculate that it was the nature of the explosion that led to the unusual demise of the SN2012au progenitor star.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On particle acceleration rate in GRB afterglows

It is well known that collisionless shocks are major sites of particle acceleration in the Universe, but the details of the acceleration process are still not well understood. The particle acceleration rate, which can shed light on the acceleration process, is rarely measured in astrophysical environments. Here we use observations of gamma-ray burst afterglows, which are weakly magnetized relativistic collisionless shocks in ion-electron plasma, to constrain the rate of particle acceleration in such shocks. We find, based on X-ray and GeV afterglows, an acceleration rate that is most likely very fast, approaching the Bohm limit, when the shock Lorentz factor is in the range of 10-100. In that case X-ray observations may be consistent with no amplification of the magnetic field in the shock upstream region. We examine the X-ray afterglow of GRB 060729, which is observed for 642 days showing a sharp decay in the flux starting about 400 days after the burst, when the shock Lorentz factor is about 5. We find that inability to accelerate X-ray emitting electrons at late time provides a natural explanation for the sharp decay, and that also in that case acceleration must be rather fast, and cannot be more than a 100 times slower than the Bohm limit. We conclude that particle acceleration is most likely fast in GRB afterglows, at least as long as the blast wave is ultra-relativistic.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Optical to X-rays supernovae light curves following shock breakout through a thick wind

Recent supernovae (SNe) observations have motivated renewed interest in SN shock breakouts from stars surrounded by thick winds. In such events the interaction with the wind powers the observed luminosity, and predictions include observable hard X-rays. Wind breakouts on timescales of a day or longer are currently the most probable for detection. Here we study the signal that follows such events. We start from the breakout of the radiation mediated shock, finding that the breakout temperature can vary significantly from one event to another ($10^4 - 5 \times 10^6$ K) due to possible deviation from thermal equilibrium. In general, events with longer breakout pulse duration, $t_{bo}$, are softer. We follow the observed radiation through the evolution of the collisionless shock which forms after the breakout of the radiation mediated shock. We restrict the study of the collisionless shock evolution to cases where the breakout itself is in thermal equilibrium, peaking in optical/UV. In these cases the post-breakout emission contains two spectral components - soft (optical/UV) and hard (X-rays and possibly soft $γ$-rays). Right after the breakout pulse X-rays are strongly suppressed, and they carry only a small fraction of the total luminosity. The hard component becomes harder and its luminosity rises quickly afterwards, gaining dominance at $\sim 10-50 \,t_{bo}$. The ratio of the peak optical/UV to the peak X-ray luminosity depends mostly on the breakout time. In terms of prospects for X-ray and soft gamma-ray detections, it is best to observe 100-500 days after explosions with breakout timescales between a week and a month.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

What Can We Learn from the Rising Lightcurves of Radioactively-Powered Supernovae?

The lightcurve of the explosion of a star with a radius <10-100Rsun is powered mostly by radioactive decay. Observationally such events are dominated by hydrogen deficient progenitors and classified as Type I supernovae: white dwarf thermonuclear explosions (Type Ia) and core collapses of hydrogen-stripped massive stars (Type Ibc). Transient surveys are finding SNe I in increasing numbers and at earlier times, allowing their early emission to be studied in unprecedented detail. Motivated by these developments, we summarize the physics that produces their rising lightcurves and discuss how observations can be utilized to study these exploding stars. The early radioactive-powered lightcurves probe the shallowest 56Ni deposits. If the amount of 56Ni mixing can be deduced, then it places constraints on the progenitor and properties of the explosive burning. In practice we find it is difficult to disentangle whether the explosion occurred recently and one is seeing 56Ni heating near the surface or whether the explosion began in the past and 56Ni heating is deeper. In the latter case there is a "dark phase" between the moment of explosion and the first light observed from the shallowest layers of 56Ni. Because of this, simply extrapolating a lightcurve from 56Ni back in time is not a reliable method for estimating the explosion time. The best solution is to directly identify the moment of explosion, by either observing shock breakout or shock-heated surface cooling, so the depth being probed by the rising lightcurve is known. Since this is typically not available, we identify other diagnostics that are helpful for deciphering how recently an explosion occurred. As an example we apply these arguments to the SN Ic PTF 10vgv. We demonstrate that just a single measurement of the photospheric velocity and temperature during the rise places constraints on its explosion time, radius, and 56Ni mixing.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Relativistic shock breakouts - a variety of gamma-ray flares: from low luminosity gamma-ray bursts to type Ia supernovae

A large variety of explosions result in mildly or ultra relativistic shock breakouts. Here we calculate the luminosity and spectrum that these breakouts produce. In order to do so we improve an analytic description of relativistic radiation mediated shocks and follow the system from the breakout itself, through the planar phase and into the spherical phase. We limit our calculation to cases where the post breakout acceleration of the gas ends during the planar phase. (i.e., the final gas Lorentz factor <~30). We find that spherical relativistic breakouts produce a flash of gamma-rays with energy, temperature, and duration that always satisfy a specific relation between them as well as providing the breakout radius and Lorentz factor. The breakout flare is typically followed, on longer time scales, by X-rays that carry a comparable energy. We apply our model to a variety of explosions, including Ia and .Ia supernovae (SNe), AIC, energetic SNe and gamma-ray bursts (GRBs). We find that all these events produce detectable gamma-ray signals, some of which may have already been seen. Some particular examples are: (i) Relativistic shock breakouts provide a natural explanation to the energy, temperature and time scales, as well as many other observed features, of all low luminosity GRBs. (ii) Nearby broad-line Ib/c (like SN 2002ap) may produce a detectable gamma-ray signal. (iii) Galactic Ia SNe produce detectable gamma-ray flares, if their progenitors are single degenerate. We conclude that relativistic shock breakout is a generic process for the production of gamma-ray flares, which opens a new window for the study and detection of a variety of stellar explosions.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

SGR 1806-20 distance and dust properties in molecular clouds by analysis of a flare x-ray echoes

The soft gamma repeater SGR 1806-20 is most famous for its giant flare from 2004, which yielded the highest gamma-ray flux ever observed on Earth. The flare emphasized the importance of determining the distance to the SGR, thus revealing the flare's energy output, with implications on SGRs energy budget and giant flare rates. We analyze x-ray scattering echoes observed by Swift/XRT following the 2006 August 6 intermediate burst of SGR 1806-20. Assuming positions and opacities of the molecular clouds along the line-of-sight from previous works, we derive direct constrains on the distance to SGR 1806-20, setting a lower limit of 9.4 kpc and an upper limit of 18.6 kpc (90% confidence), compared with a 6-15 kpc distance range by previous works. This distance range matches an energy output of ~10^46 erg/s for the 2004 giant flare. We further use, for the first time, the x-ray echoes in order to study the dust properties in molecular clouds. Analyzing the temporal evolution of the observed flux using a dust scattering model, which assumes a power-law size distribution of the dust grains, we find a power-law index of -3.3_{-0.7}^{+0.6} (1 sigma) and a lower limit of 0.1 micron (2 sigma) on the dust maximal grain size, both conforming to measured dust properties in the diffused interstellar medium (ISM). We advocate future burst follow-up observations with Swift, Chandra and the planned NuSTAR telescopes, as means of obtaining much superior results from such an analysis.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Two-stream-like instability in dilute hot relativistic beams and astrophysical relativistic shocks

Relativistic collisionless shocks are believed to be efficient particle accelerators. Nonlinear outcome of the interaction of accelerated particles that run ahead of the shock, the so-called "precursor", with the unperturbed plasma of the shock upstream, is thought to facilitate additional acceleration of these particles and to possibly modify the hydrodynamic structure of the shock. We explore here the linear growth of kinetic modes appearing in the precursor-upstream interaction in relativistic shocks propagating in non and weakly magnetized plasmas: electrostatic two-stream parallel mode and electrostatic oblique modes. These modes are of particular interest because they are the fastest growing modes known in this type of system. Using a simplified distribution function for a dilute ultra-relativistic beam that is relativistically hot in its own rest frame, yet has momenta that are narrowly collimated in the frame of the cold upstream plasma into which it propagates, we identify the fastest growing mode in the full $k$-space and calculate its growth rate. We consider all types of plasma (pairs and ions-electrons) and beam (charged and charge-neutral). We find that unstable electrostatic modes are present in any type of plasma and for any shock parameters. We further find that two modes, one parallel ($k_\perp=0$) and the other one oblique ($k_\perp \sim k_\|$), are competing for dominance and that either one may dominate the growth rate in different regions of the phase space. The dominant mode is determined mostly by the perpendicular spread of the accelerated particle momenta in the upstream frame, which reflects the shock Lorentz factor. The parallel mode becomes more dominant in shocks with lower Lorentz factors (i.e., with larger momentum spreads). We briefly discuss possible implications of our results for external shocks in gamma-ray burst sources.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Early supernovae light-curves following the shock-breakout

The first light from a supernova (SN) emerges once the SN shock breaks out of the stellar surface. The first light, typically a UV or X-ray flash, is followed by a broken power-law decay of the luminosity generated by radiation that leaks out of the expanding gas sphere. Motivated by recent detection of emission from very early stages of several SNe, we revisit the theory of shock breakout and the following emission. We derive analytic light curves, paying special attention to the photon-gas coupling and deviations from thermal equilibrium. We then consider the breakout from several SNe progenitors. We find that for more compact progenitors, white dwarfs, Wolf-Rayet stars (WRs) and possibly more energetic blue-supergiant explosions, the observed radiation is out of thermal equilibrium at the breakout, during the planar phase (i.e., before the expanding gas doubles its radius), and during the early spherical phase. Therefore, during these phases we predict significantly higher temperatures than previous analysis that assumed equilibrium. When thermal equilibrium prevails, we find the location of the thermalization depth and its temporal evolution. Our results are useful for interpretation of early SN light curves. Some examples are: (i) Red supergiant SNe have an early bright peak in optical and UV flux, less than an hour after breakout. It is followed by a minimum at the end of the planar phase (about 10 hr), before it peaks again once the temperature drops to the observed frequency range. In contrast WRs show only the latter peak in optical and UV. (ii) Bright X-ray flares are expected from all core-collapse SNe types. (iii) The light curve and spectrum of the initial breakout pulse holds information on the explosion geometry and progenitor wind opacity. Its spectrum in compact progenitors shows a (non-thermal) power-law.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Radio and X-ray Observations of the Type Ic SN 2007gr Reveal an Ordinary, Non-relativistic Explosion

We present extensive radio and X-ray observations of the nearby Type Ic SN 2007gr in NGC 1058 obtained with the Very Large Array and the Chandra X-ray Observatory and spanning 5 to 150 days after explosion. Through our detailed modeling of these data, we estimate the properties of the blastwave and the circumstellar environment. We find evidence for a freely-expanding and non-relativistic explosion with an average blastwave velocity, v~0.2c, and a total internal energy for the radio emitting material of E ~ 2 x 10^46 erg assuming equipartition of energy between electrons and magnetic fields (epsilon_e=epsilon_B=0.1). The temporal and spectral evolution of the radio emission points to a stellar wind-blown environment shaped by a steady progenitor mass loss rate of Mdot ~ 6 x 10^-7 solar masses per year (wind velocity, v_w=10^3 km/s). These parameters are fully consistent with those inferred for other SNe Ibc and are in line with the expectations for an ordinary, homologous SN explosion. Our results are at odds with those of Paragi et al. (2010) who recently reported evidence for a relativistic blastwave in SN 2007gr based on their claim that the radio emission was resolved away in a low signal-to-noise Very Long Baseline Interferometry (VLBI) observation. Here we show that the exotic physical scenarios required to explain the claimed relativistic velocity -- extreme departures from equipartition and/or a highly collimated outflow -- are excluded by our detailed Very Large Array radio observations. Moreover, we present an independent analysis of the VLBI data and propose that a modest loss of phase coherence provides a more natural explanation for the apparent flux density loss which is evident on both short and long baselines. We conclude that SN 2007gr is an ordinary Type Ibc supernova.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

What do we know about gamma-ray bursts?

Decades of improving data and extensive theoretical research have led to a popular model of gamma-ray bursts. According to this model, a catastrophic event in a stellar system results in the formation of a compact central engine, which releases a fraction of a solar rest-mass energy within seconds in the form of ultra-relativistic jets. Dissipation of the jets energy leads first to prompt gamma-ray emission and later to a long lasting afterglow. Here I summarize the introduction that I gave to the debate "where do we stand?" in the conference "The Shocking Universe" held in Venice. This is a very brief summary of my view of the facts that we are (almost) certain about, models that are popular but may need rethinking, and main open questions.

0 citations0 reviewsNo code linkedOpen access
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?

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Detectability of Oort cloud objects using Kepler

The size distribution and total mass of objects in the Oort Cloud have important implications to the theory of planets formation, including the properties of, and the processes taking place in the early solar system. We discuss the potential of space missions like Kepler and CoRoT, designed to discover transiting exo-planets, to detect Oort Cloud, Kuiper Belt and main belt objects by occultations of background stars. Relying on published dynamical estimates of the content of the Oort Cloud, we find that Kepler's main program is expected to detect between 0 and ~100 occultation events by deca-kilometer-sized Oort Cloud objects. The occultations rate depends on the mass of the Oort cloud, the distance to its "inner edge", and the size distribution of its objects. In contrast, Kepler is unlikely to find occultations by Kuiper Belt or main belt asteroids, mainly due to the fact that it is observing a high ecliptic latitude field. Occultations by Solar System objects will appear as a photometric deviation in a single measurement, implying that the information regarding the time scale and light-curve shape of each event is lost. We present statistical methods that have the potential to verify the authenticity of occultation events by Solar System objects, to estimate the distance to the occulting population, and to constrain their size distribution. Our results are useful for planning of future space-based exo-planet searches in a way that will maximize the probability of detecting solar system objects, without hampering the main science goals.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

LSST Science Book, Version 2.0

A survey that can cover the sky in optical bands over wide fields to faint magnitudes with a fast cadence will enable many of the exciting science opportunities of the next decade. The Large Synoptic Survey Telescope (LSST) will have an effective aperture of 6.7 meters and an imaging camera with field of view of 9.6 deg^2, and will be devoted to a ten-year imaging survey over 20,000 deg^2 south of +15 deg. Each pointing will be imaged 2000 times with fifteen second exposures in six broad bands from 0.35 to 1.1 microns, to a total point-source depth of r~27.5. The LSST Science Book describes the basic parameters of the LSST hardware, software, and observing plans. The book discusses educational and outreach opportunities, then goes on to describe a broad range of science that LSST will revolutionize: mapping the inner and outer Solar System, stellar populations in the Milky Way and nearby galaxies, the structure of the Milky Way disk and halo and other objects in the Local Volume, transient and variable objects both at low and high redshift, and the properties of normal and active galaxies at low and high redshift. It then turns to far-field cosmological topics, exploring properties of supernovae to z~1, strong and weak lensing, the large-scale distribution of galaxies and baryon oscillations, and how these different probes may be combined to constrain cosmological models and the physics of dark energy.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Cluster Merger Shock Constraints on Particle Acceleration and Nonthermal Pressure in the Intracluster Medium

X-ray observations of galaxy cluster merger shocks can be used to constrain nonthermal processes in the intracluster medium (ICM). The presence of nonthermal pressure components in the ICM, as well as the shock acceleration of particles and their escape, all affect shock jump conditions in distinct ways. Therefore, these processes can be constrained using X-ray surface brightness and temperature maps of merger shock fronts. Here we use these observations to place constraints on particle acceleration efficiency in intermediate Mach number (M ~ 2-3) shocks and explore the potential to constrain the contribution of nonthermal components (e.g., cosmic rays, magnetic field, and turbulence) to ICM pressure in cluster outskirts. We model the hydrodynamic jump conditions in merger shocks discovered in the galaxy clusters A520 (M ~ 2) and 1E 0657-56 (M ~ 3) using a multifluid model comprised of a thermal plasma, a nonthermal plasma, and a magnetic field. Based on the published X-ray spectroscopic data alone, we find that the fractional contribution of cosmic rays accelerated in these shocks is lower than about 10% of the shock downstream pressure. Current observations do not constrain the fractional contribution of nonthermal components to the pressure of the undisturbed shock upstream. Future X-ray observations, however, have the potential to either detect particle acceleration in these shocks through its effect on the shock dynamics, or to place a lower limit on the nonthermal pressure contributions in the undisturbed ICM. We briefly discuss implications for models of particle acceleration in collisionless shocks and the estimates of galaxy cluster masses derived from X-ray and Sunyaev-Zel'dovich effect observations.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Short-Hard Gamma-Ray Bursts

Two types of Gamma-ray bursts (GRBs) are observed: short duration hard spectrum GRBs and long duration soft spectrum GRBs. For many years long GRBs were the focus of intense research while the lack of observational data limited the study of short-hard GRBs (SHBs). In 2005 a breakthrough occurred following the first detections of SHB afterglows, longer wavelength emission that follows the burst of gamma-rays. Similarly to long GRBs, afterglow detections led to the identification of SHB host galaxies and measurement of their redshifts. These observations established that SHBs are cosmological relativistic sources that, unlike long GRBs, do not originate from the collapse of massive stars, and therefore constitute a distinct physical phenomenon. One viable model for SHB origin is the coalescence of compact binary systems (double neutron stars or a neutron star and a black hole), in which case SHBs are the electromagnetic counterparts of strong gravitational-wave sources. The theoretical and observational study of SHBs following the recent pivotal discoveries is reviewed, along with new theoretical results that are presented here for the first time.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

A Complete Survey of the Transient Radio Sky and Implications for Gamma-Ray Bursts, Supernovae, and other Relativistic Explosions

We had previously reported on a survey for radio transients, used to set an upper limit on the number of orphan gamma-ray burst (GRB) radio afterglows, and thus a lower limit on the typical GRB beaming factor. Here we report radio and optical follow-up observations of these possible transients, achieving the first full characterization of the transient radio sky. We find that only two source are likely to be real radio transients, an optically obscured radio supernova (SN) in the nearby galaxy NGC 4216, and a source not associated with a bright host galaxy, which is too radio luminous to be a GRB afterglow. We speculate that this may be a flare from a peculiar active galactic nucleus, or a burst from an unusual Galactic compact object. We place an upper limit of 65 radio transients above 6 mJy over the entire sky at the 95% confidence level. The implications are as follows. First, we derive a limit on the typical beaming of GRBs; we find f_b^{-1} >~ 60, ~5 times higher than our earlier results. Second, we impose an upper limit on the rate of events that eject >~ 10^{51} erg in unconfined relativistic ejecta, whether or not accompanied by detectable emission in wavebands other than the radio. Our estimated rate, <=1000/y/Gpc, is about two orders of magnitude smaller than the rate of core-collapse SNe (and type Ib/c events in particular), indicating that only a minority of such events eject significant amounts of relativistic material, which are required by fireball models of long-soft GRBs. Finally, we show that future wider and/or deeper radio variability surveys are expected to detect numerous orphan radio GRB afterglows. Our survey also illustrates the great potential of sensitive surveys with new instruments to revolutionize the study of nearby SNe (abridged).

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

The Progenitors of Short-Hard Gamma-Ray Bursts from an Extended Sample of Events

The detection of the afterglow emission and host galaxies of short-hard gamma-ray bursts (SHBs) is one of the most exciting recent astronomical discoveries. Indications that SHB progenitors belong to old stellar populations, in contrast to those of the long-soft GRBs, provide a strong clue about their physical nature. Definitive conclusions however are limited by the small number of SHBs with known hosts. Here, we present our investigation of SHBs previously localized by the interplanetary network (IPN) using new and archival optical and X-ray observations. We show that we can likely identify the host galaxies/clusters for additional two bursts, significantly increasing the sample of SHBs with known hosts and/or distances. In particular, we determine that the bright SHB 790613 occurred within the rich galaxy cluster Abell 1892, making it probably the nearest SHB currently known. We show that the brightest galaxy within the error box of SHB 000607, at z=0.14, is most likely the host galaxy of this event. Additionally, we rule out the existence of galaxy overdensities (down to ~21mag) near the locations of two other SHBs, and set a lower limit on their probable redshift. We combine our SHB sample with events discovered recently by the Swift and HETE-2 missions, and investigate the properties of the extended sample. Comparison to SNe Ia shows that the progenitors of SHBs are typically older, implying a typical life time of several Gy. We also show that it is unlikely that there is a significant population of progenitors with life time shorter than 1Gy. This result disfavors the popular model of NS-NS mergers as the progenitors of SHBs.

0 citations0 reviewsNo code linkedOpen access