One of the key properties of any binary is its viewing angle (i.e., inclination), $\theta_{\rm obs}$. In binary neutron star (BNS) mergers it is of special importance due to the role that it plays in the measurement of the Hubble constant, $H_0$. The opening angle of the jet that these mergers launch, $\theta_j$, is also of special interest. Following the detection of the first BNS merger, GW170817, there were numerous attempts to estimate these angles using the afterglow light curve, finding a wide range of values. Here we provide a simple formula for the ratio $\theta_{\rm obs}/\theta_j$ based on the afterglow light curve and show that this is the only quantity that can be determined from the light curve alone. Namely, it is impossible to determine each of the angles separately without additional information. Our result explains the inconsistency of the values found by the various studies of GW170817 that were largely driven by the different priors taken in each study. Among the additional information that can be used to estimate $\theta_{\rm obs}$ and $\theta_j$, the most useful is a VLBI measurement of the afterglow image superluminal motion. An alternative is an identification of the afterglow transition to the sub-relativistic phase. These observations are possible only for mergers observed at small viewing angles, whose afterglow is significantly brighter than the detector's threshold. We discuss the implications of these results to measurements of $H_0$ using GW observations. We show that while the viewing angle will be measured only in a small fraction of future BNS mergers, it can significantly reduce the uncertainty in $H_0$ in each one of these events, possibly to a level of 4-5\%. A minority of the mergers with high precision measurements of this kind may dominate in the future the precision in which $H_0$ will be measured using this method.

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

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

The black hole of an active galactic nucleus is encircled by an accretion disk. The surface density of the disk is always too low to affect the tidal disruption of a star, but it can be high enough that a vigorous interaction results when the debris stream returns to pericenter and punches through the disk. Shocks excited in the disk dissipate the kinetic energy of the disk interior to the impact point and expedite inflow toward the black hole. Radiatively efficient disks with luminosity $\gtrsim10^{-3}$ Eddington have high enough surface density that the initial stream-disk interaction leads to energy dissipation at a super-Eddington rate. Because of the rapid inflow, only part of this dissipated energy emerges as radiation, while the rest is advected into the black hole. Dissipation, inflow, and cooling balance to keep the bolometric luminosity at an Eddington-level plateau that lasts tens of days. After the plateau, the luminosity decreases in proportion to the disk surface density, with a power-law index between $-3$ and $-2$ at earlier times, and possibly a steeper index at later times.

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

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

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

Tidal disruption events (TDEs), events in which a star passes very close to a supermassive black hole, are generally imagined as leading either to the star's complete disruption or to its passage directly into the black hole. In the former case it is widely believed that in all cases the bound portion of the debris quickly "circularizes" due to relativistic apsidal precession, i.e., forms a compact accretion disk, and emits a flare of standardized lightcurve and spectrum. We show here that TDEs are more diverse and can be grouped into several distinct categories on the basis of stellar pericenter distance $r_p$; we calculate the relative frequency of these categories. In particular, because rapid circularization requires $r_p \lesssim 10r_g$ ($r_g \equiv GM_{\rm BH}/c^2$), it can happen in only a minority of total disruptions, $\lesssim 1/4$ when the black hole has mass $M_{\rm BH} = 10^6 M_\odot$. For larger pericenter distances, $10 < r_p/r_g < 27$ (for $M_{\rm BH}=10^6M_\odot$), main sequence stars are completely disrupted, but the bound debris orbits are highly eccentric and possess semimajor axes $\sim 100\times$ the scale of the expected compact disk. Partial disruptions with fractional mass-loss $\gtrsim 10\%$ should occur with a rate similar to that of total disruptions; for fractional mass-loss $\gtrsim 50\%$, the rate is $\approx 1/3$ as large. Partial disruptions---which must precede total disruptions when the stars' angular momenta evolve in the "empty loss-cone" regime---change the orbital energy by factors $\gtrsim O(1)$. Remnants of partial disruptions are in general far from thermal equilibrium. Depending on the orbital energy of the remnant and conditions within the stellar cluster surrounding the SMBH, it may return after hundreds or thousands of years and be fully disrupted, or it may rejoin the stellar cluster.

This paper introduces a series of papers presenting a quantitative theory for the tidal disruption of main sequence stars by supermassive black holes. Using fully general relativistic hydrodynamics simulations and MESA-model initial conditions, we explore the pericenter-dependence of tidal disruption properties for eight stellar masses ($0.15 \leq M_*/M_\odot \leq 10$) and six black hole masses ($10^5 \leq M_{BH}/M_\odot \leq 5 \times 10^7$). We present here the results most relevant to observations. The effects of internal stellar structure and relativity decouple for both the disruption cross section and the characteristic energy width of the debris. Moreover, the full disruption cross section is almost independent of $M_*$ for $M_*/M_\odot \lesssim 3$. Independent of $M_*$, relativistic effects increase the critical pericenter distance for full disruptions by up to a factor $\sim 3$ relative to the Newtonian prediction. The probability of a direct capture is also independent of $M_*$; at $M_{BH}/M_\odot \simeq 5 \times 10^6$ this probability is equal to that of a complete disruption. The width of the debris energy distribution $\Delta E$ can differ from the standard estimate by factors from 0.35 to 2, depending on $M_*$ and $M_{BH}$, implying a corresponding change in the characteristic mass-return timescale. The "frozen-in approximation" is inconsistent with $\Delta E$, and mass-loss continues over a long span of time. We provide analytic forms, suitable for use in both event rate estimates and parameter inference, to describe all these trends. For partial disruptions, we find a nearly-universal relation between the star's angular momentum and the fraction of $M_*$ remaining. Within the "empty loss-cone" regime, partial disruptions must precede full disruptions. These partial disruptions can drastically affect the rate and appearance of subsequent total disruptions.

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

Using a suite of fully relativistic hydrodynamic simulations applied to main-sequence stars with realistic internal density profiles, we examine full and partial tidal disruptions across a wide range of black hole mass ($10^{5}\leq M_{\rm BH}/\mathrm{M}_{\odot}\leq 5\times 10^{7}$) and stellar mass ($0.3 \leq M_{\star} /\mathrm{M}_{\odot}\leq 3$) as larger $M_{\rm BH}$ leads to stronger relativistic effects. For fixed $M_{\star}$, as $M_{\rm BH}$ increases, the ratio of the maximum pericenter distance yielding full disruptions ($\mathcal{R}_{\rm t}$) to its Newtonian prediction rises rapidly, becoming triple the Newtonian value for $M_{\rm BH} = 5\times10^{7}~{\rm M}_\odot$, while the ratio of the energy width of the stellar debris for full disruptions to the Newtonian prediction decreases steeply, resulting in a factor of two correction at $M_{\rm BH} = 5 \times 10^7~{\rm M}_\odot$. We find that for partial disruptions, the fractional remnant mass for a given ratio of the pericenter to $\mathcal{R}_{\rm t}$ is higher for larger $M_{\rm BH}$. These results have several implications. As $M_{\rm BH}$ increases above $\sim 10^7~{\rm M}_\odot$, the cross section for complete disruptions is suppressed by competition with direct capture. However, the cross section ratio for partial to complete disruptions depends only weakly on $M_{\rm BH}$. The relativistic correction to the debris energy width delays the time of peak mass-return rate and diminishes the magnitude of the peak return rate. For $M_{\rm BH} \gtrsim 10^7~{\rm M}_\odot$, the $M_{\rm BH}$-dependence of the full disruption cross section and the peak mass-return rate and time is influenced more by relativistic effects than by Newtonian dynamics.

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

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

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

A $\gamma$-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 ($\Gamma \gtrsim 100$) giving the first and most important clue on their nature. So far, this condition has been discussed under the assumption that the $\gamma$-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 $\gamma$-ray sources are viewed off-axis. We generalize here the compactness formalism for an arbitrary viewing angle taking several possible opacity processes and $\gamma$-ray spectra into account. We find that for a given observables (peak luminosity, temporal variability, and spectra) the minimal Lorentz factor, $\Gamma_{\rm min}$, is obtained, as expected, for an on-axis observer. More remarkably we find that compactness dictates also a maximal viewing angle, $\theta_{\rm max} \simeq 1/2\Gamma_{\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 $\gamma$-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.

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

GrailQuest (Gamma-ray Astronomy International Laboratory for Quantum Exploration of Space-Time) is an ambitious astrophysical mission concept that uses a fleet of small satellites, whose scientific objectives are discussed below. Within Quantum Gravity theories, different models for space-time quantisation predict an energy dependent speed for photons. Although the predicted discrepancies are minuscule, Gamma-Ray Bursts, occurring at cosmological distances, could be used to detect this signature of space-time granularity with a new concept of modular observatory of huge overall collecting area consisting in a fleet of small satellites in low orbits, with sub-microsecond time resolution and wide energy band (keV-MeV). The enormous number of collected photons will allow to effectively search these energy dependent delays. Moreover, GrailQuest will allow to perform temporal triangulation of high signal-to-noise impulsive events with arc-second positional accuracies: an extraordinary sensitive X-ray/Gamma all-sky monitor crucial for hunting the elusive electromagnetic counterparts of Gravitational Waves. A pathfinder of GrailQuest is already under development through the HERMES (High Energy Rapid Modular Ensemble of Satellites) project: a fleet of six 3U cube-sats to be launched by the end of 2021.

Binary neutron stars (BNS) mergers are prime sites for $r$-process nucleosynthesis. Their rate determines the chemical evolution of heavy elements in the Milky Way. The merger rate of BNS is a convolution of their birth rate and the gravitational radiation spiral-in delay time. Using the observed population of Galactic BNS we show here that the lifetimes of pulsars in observed BNSs are sufficiently short that the ages of BNSs have little to no effect on the observed merger time distribution. We find that at late times ($t\gtrsim 1$ Gyr) the gravitational wave delay time distribution (DTD) follows the expected $ t^{-1}$. However, a significant excess of rapidly merging systems (between $40-60\%$ of the entire population) is apparent at shorter times. Although the exact shape of the DTD cannot be determined with the existing data, in all models that adequately describe the data we find at least $40\%$ of BNSs with merger times less than 1Gyr. This population of rapid mergers implies a declining deposition rate of $r$-process materials that is consistent with several independent observations of heavy element abundances in the Milky Way. At the same time this population that requires initial binary separations of roughly one solar radius clearly indicates that these binaries had common envelope progenitors. Our results suggest that a significant fraction of future LIGO/Virgo BNS mergers would reside in star forming galaxies.

The MAGIC observations of late sub-TeV photons from GRB190114C enable us, for the first time, to determine the details of the emission process in a GRB afterglow and to pin down the physical parameters, such as the bulk Lorentz factor and the Lorentz factor of the emitting electrons as well as some of the microphysical parameters. We find that the sub-TeV emission is an Inverse Compton radiation produced at the early afterglow stage. Combining the sub-TeV and X-ray observations we narrow uncertainties in the conditions inside the emitting zone, almost eliminating them for some parameters. Seventy seconds after the trigger the forward shock had Lorentz factor $\simeq 100$ and the electrons producing the observed sub-TeV radiation had a Lorentz factor $\simeq 10^4$, so that the sub-TeV radiation originates from Comptonization of X-ray photons at the border between Thomson and Klein-Nishina regimes. The inferred conditions within the emitting zone are at odds with theoretical expectations unless one assumes moderate (with $\tau \simeq 2$) absorption of sub-TeV photons inside the source. With this correction the conditions are in good agreement with predictions of the pair-balance model, but are also acceptable for generic afterglow model as one of many possibilities. However, the temporal evolution of the Inverse Compton peak energyis very different in these two models. This opens a way to discriminate between the models once late-time detection in the TeV range become available.

Observations of several Fanaroff-Riley (FR) type I sources reveal outflowing bipolar bubbles of hot gas surrounded by a weak forward shock. We consider the possibility that these bubbles were driven by choked relativistic jets which failed to penetrate the ambient intracluster medium (ICM). Using new results on choked jets linking the geometry of the forward shock to the jet properties, we infer robust limits on the radius $R_{ch}$ at which the jet was quenched in 5 well-studied FRI sources, finding typically $R_{ch}\sim 10$ kpc. We further show that, in order to reach this radius in less than the current age of the system, the jet must have been tightly collimated, with the jet head subtending an angle of $\theta_h < 2^\circ$. The ambient pressure is not high enough to explain this collimation, suggesting that the jet was collimated by interaction with its own cocoon. Although the choking radius is well-constrained, we find a degeneracy between the initial jet opening angle before collimation, $\theta_0$, and the duration of jet activity, $t_b$, with $(t_b/1\rm{Myr})(\theta_0/5^\circ)^{-2} \sim 0.1$. We speculate that the working time and/or opening angle of the jet may be important factors contributing to the FR type I/type II morphology in galaxy clusters, with short-lived or wide jets being choked to form bipolar bubbles filled with diffuse radio emission, and longer-lived or narrow jets successfully escaping the cluster core to produce cocoons with radio hotspots.

Observations of both gamma-ray bursts (GRBs) and active galactic nuclei (AGNs) point to the idea that some relativistic jets are suffocated by their environment before we observe them. In these "choked" jets, all the jet's kinetic energy is transferred into a hot and narrow cocoon of near-uniform pressure. We consider the evolution of an elongated, axisymmetric cocoon formed by a choked jet as it expands into a cold power-law ambient medium $\rho \propto R^{-\alpha}$, in the case where the shock is decelerating ($\alpha<3$). The evolution proceeds in three stages, with two breaks in behaviour: the first occurs once the outflow has doubled its initial width, and the second once it has doubled its initial height. Using the Kompaneets approximation, we derive analytical formulae for the shape of the cocoon shock, and obtain approximate expressions for the height and width of the outflow versus time in each of the three dynamical regimes. The asymptotic behaviour is different for flat ($\alpha \le 2$) and steep ($2 < \alpha < 3$) density profiles. Comparing the analytical model to numerical simulations, we find agreement to within $\sim15$ per cent out to 45 degrees from the axis, but discrepancies of a factor of 2-3 near the equator. The shape of the cocoon shock can be measured directly in AGNs, and is also expected to affect the early light from failed GRB jets. Observational constraints on the shock geometry provide a useful diagnostic of the jet properties, even long after jet activity ceases.

VLBI and JVLA observations revealed that GW170817 involved a narrow jet ($ \theta_j \approx 4^\circ $) that dominated the afterglow peak at our viewing angle, $ \theta_{\rm obs} \approx 20^\circ $. This implies that at the time of the afterglow peak, the observed signal behaved like an afterglow of a top-hat jet seen at $ \theta_{\rm obs} \gg \theta_j $, and it can be modeled by analytic expressions that describe such jets. We use a set of numerical simulations to calibrate these analytic relations and obtain generic equations for the peak time and flux of such an afterglow as seen from various observing angles. Using the calibrated equations and the estimated parameters of GW170817, we estimate the detectability of afterglows from future double neutron star mergers during the Advanced LIGO/Virgo observation run O3. GW170817 took place at a relatively low-density environment. Afterglows of similar events will be detectable only at small viewing angles, $ \theta_{\rm obs} \lesssim 20^\circ $, and only $\sim 20\% $ of the GW detections of these events will be accompanied by a detectable afterglow. At higher densities, more typical to sGRB sites, up to $ 70\% $ of the GW detections are expected to be followed by a detectable afterglow, typically at $ \theta_{\rm obs} \sim 30^\circ $. We also provide the latest time one should expect an afterglow detection. We find that for typical parameters, if the jet emission had not been detected within about a year after the merger, it is unlikely to be ever detected.

When a star gets too close to a supermassive black hole, it is torn apart by the tidal forces. Roughly half of the stellar mass becomes unbound and flies away at tremendous velocities - around $10^4$ km/s. In this work we explore the idea that the shock produced by the interaction of the unbound debris with the ambient medium gives rise to the synchrotron radio emission observed in several TDEs. We use a moving mesh numerical simulation to study the evolution of the unbound debris and the bow shock around it. We find that as the periapse distance of the star decreases, the outflow becomes faster and wider. A tidal disruption event whose periapse distance is a factor of 7 smaller than the tidal radius can account for the radio emission observed in ASASSN-14li. This model also allows us to obtain a more accurate estimate for the gas density around the centre of the host galaxy of ASASSN-14li.

Motivated by the ~7 % dipole anisotropy in the distribution of ultra-high energy cosmic-rays (UHECRs) above 8 EeV, we explore the anisotropy induced by the large scale structure, using constrained simulations of the local Universe and taking into account the effect of magnetic fields. The value of the intergalactic magnetic field (IGMF) is critical as it determines the UHECR cosmic horizon. We calculate the UHECR sky maps for different values of the IGMF variance and show the effect of the UHECR horizon on the observed anisotropy. The footprint of the local ($\lesssim350$ Mpc) Universe on the UHECR background, a small angular scale enhancement in the Northern Hemisphere, is seen. At 11.5 EeV (the median value of the energy bin at which the dipole has been reported), the LSS-induced dipole amplitude is $A_1\sim10\%$, for IGMF in the range [0.3-3] nG for protons, helium and nitrogen, compatible with the rms value derived from the cosmic power spectrum. However at these energies the UHECRs are also influenced by the Galactic Magnetic Field (GMF) and we discuss its effect on the LSS-induced anisotropy.

Nearby pulsars have been suggested as sources of $\sim$TeV $e^+/e^-$ Cosmic Ray (CR) excess on Earth. The High-Altitude Water Cherenkov Observatory (HAWC) detected extended TeV emission regions in the direction of two nearby middle-aged pulsars, Geminga and PSR B0656+14. By modeling the TeV emission as inverse Compton emission from electron-positron pairs diffusing in the interstellar medium (ISM), the HAWC collaboration derives a diffusion coefficient much smaller than the standard value in the vicinity of the two pulsars, which make them unlikely the origin of the positron excess. We propose that the observed $\gamma$-ray emission originate from the relic pulsar wind nebula. A two zone diffusion model with a slow diffusion in the nebula and a fast diffusion in the ISM can explain the HAWC surface brightness profile and the positron excess simultaneously. Inefficient diffusion in the $\gamma$-ray emission region surrounding a middle-aged pulsar maybe a common phenomenon that can be tested by future observation. The implied diffusion coefficient in the ISM is smaller than the one suggested by the standard CR propagation model, but it is fully consistent with the predictions of the spiral arm model.

The gravitational waves from the neutron star merger event GW170817 were accompanied by an unusually weak short GRB 170817A, by an optical/IR macronova/kilonova and by a long lasting radio to X-rays counterpart. While association of short GRBs with mergers was predicted a long time ago, the luminosity of this prompt {\gamma}-ray emission was weaker by a few orders of magnitude than all known previous sGRBs and it was softer than typical sGRBs. This raise the question whether the {\gamma}-rays that we have seen were a regular sGRB viewed off-axis. We revisit this question following recent refined analyses of the {\gamma}-ray signal and the VLBI observations that revealed the angular structure of the relativistic outflow: observing angle of $\sim\,20^\circ$, a narrow jet with core $\lesssim5^\circ$ and $E_{\rm iso}>10^{52}\,{\rm ergs}$. We show here that: (i) The region emitting the observed $\gamma$-rays must have been moving with a Lorentz factor $\Gamma\gtrsim5$; (ii) The observed {\gamma}-rays were not "off-axis" emission (viewing angle $>1/\Gamma$) emerging from the core of the jet, where a regular sGRB was most likely produced; (iii) The $\gamma$-ray emission region was either "on-axis" (at an angle $<1/\Gamma$) or if it was "off-axis" then the observing angle must have been small ($<5^\circ$) and the on-axis emission from this region was too faint and too hard to resemble a regular sGRB.

LIGO's detection of gravitational waves from binary black hole mergers was an unexpected surprise that immediately raised the question - what is the origin of these black hole binaries? The "simplest" scenario is evolution of field massive stellar binaries. However, other possibilities involving capture have been proposed. We explore here one of the more interesting clues on this puzzle: the relatively modest spins of the resulting black holes that imply that the progenitor black holes were not spinning rapidly. More specifically we consider the implication of observed distribution of, $\chi_{\rm eff}$, the mass weighted projected (along the orbital axis) spins on the field evolution scenario. In all cases $\chi_{\rm eff}$ is small and in two of the cases the best fit value is negative. Only in one event the spin is positive at 90\% credible. These observations are puzzling within the field binary scenario in which positive higher spins ($\chi_{\rm eff} \ge 0.5$) are expected. At first sight one may expect that this rules out the field evolutionary scenario. Indeed we show that with typical parameters a significant fraction ($\ge 25\%$) of the mergers should have high effective spin values. However, uncertainties in the outcome of the common envelope phase (the typical separation and whether the stars are rotating or not) and in the late stages of massive star evolution (the strength of the winds) make it impossible to rule out, at present, these scenarios. While observations of mergers with high effective spin will support this scenario, future observations of negative spin mergers would rule it out.

We study the prospect of the Gaia satellite to identify black hole binary systems by detecting the orbital motion of the companion stars. Taking into account the initial mass function, mass transfer, common envelope phase, interstellar absorption and identifiability of black holes, we estimate the number of black hole binaries detected by Gaia and their distributions with respect to the black hole mass for several models with different parameters. We find that $\sim 200-1000$ black hole binaries will be detected by Gaia during its 5 years operation. The shape of distribution function of the black hole mass is affected most severely by the relation between the zero-age main sequence (ZAMS) stellar mass and the black hole mass in parameters adopted in this paper, which implies that black holes detected with Gaia enables us to constrain the mass relation.

The short Gamma-Ray Burst, GRB170817A, that followed the binary neutron star merger gravitational waves signal, GW170817, is not a usual sGRB. It is weaker by three orders of magnitude than the weakest sGRB seen before and its spectra, showing a hard early signal followed by a softer thermal spectrum, is unique. We show, first, that the $\gamma$-rays must have emerged from at least mildly relativistic outflow, implying that a relativistic jet was launched following the merger. We then show that the observations are consistent with the predictions of a mildly relativistic shock breakout: a minute $\gamma$-ray energy as compared with the total energy and a rather smooth light curve with a hard to soft evolution. We present here a novel analytic study and detailed numerical 2D and 3D relativistic hydrodynamic and radiation simulations that support the picture in which the observed $\gamma$-rays arose from a shock breakout of a cocoon from the merger's ejecta (Kasliwal 2017). The cocoon can be formed by either a choked jet which does not generate a sGRB (in any direction) or by a successful jet which generates an undetected regular sGRB along the system's axis pointing away from us. Remarkably, for the choked jet model, the macronova signal produced by the ejecta (which is partially boosted to high velocities by the cocoon's shock) and the radio that is produced by the interaction of the shocked cocoon material with the surrounding matter, agree with the observed UV/optical/IR emission and with current radio observations. Finally, we discuss the possibility that the jet propagation within the ejecta may photodissociate some of of the heavy elements and may affect the composition of a fraction of ejecta and, in turn, the opacity and the early macronova light.

The radio and X-rays that followed GW170817 increased gradually over $\sim~100$ days, resembling the radio flare predicted to arise from the interaction of a binary neutron star merger outflow with the ISM \citep{nakar2011}. Considering a blast wave moving with a Lorentz factor $\Gamma$, we show that an off-axis observer, namely an observer at $\theta_{obs}>1/\Gamma$, sees a light curve rising faster than $F_\nu\propto t^3$. Therefore, the observed rise, $F_\nu\propto~t^{0.78}$, implies that at all times we have seen an on-axis emission. Namely, the emitting matter was within $\theta_{obs}<1/\Gamma$ at the time of observations (even if it was off-axis beforehand). The observations tightly constrain the blast wave Lorentz factor: $\Gamma\sim(1.5-7)(t/10~{\rm day})^{-0.21}$.The isotropic equivalent energy in the observed region increases with time, $E_{iso}\sim~10^{50}{\rm~erg} ~(t/150~{\rm day})^{1.3}$, implying that the outflow was structured. The energy increase can arise from a slower material moving behind the blast wave or from a matter moving at larger angles that has slowed down. Such a structure could have different origins, however, the only physically motivated one, proposed so far, is the interaction of a relativistic jet with the ejecta and the resulting cocoon. The jet could have been choked or successful. In the latter case, it has produced a short GRB pointing elsewhere (this successful jet-cocoon system is sometimes called a "structured jet"). Although circumstantial evidence disfavors a successful jet, the jet fate cannot be decisively determined from current observations. The light curve alone may not be sufficient to resolve this question, since both chocked and successful jets can lead to a gradual rise to a peak, followed by a decay. Therefore, the recent turnover of the light curve does not necessitate a successful jet.

THESEUS is a space mission concept aimed at exploiting Gamma-Ray Bursts for investigating the early Universe and at providing a substantial advancement of multi-messenger and time-domain astrophysics. These goals will be achieved through a unique combination of instruments allowing GRB and X-ray transient detection over a broad field of view (more than 1sr) with 0.5-1 arcmin localization, an energy band extending from several MeV down to 0.3 keV and high sensitivity to transient sources in the soft X-ray domain, as well as on-board prompt (few minutes) follow-up with a 0.7 m class IR telescope with both imaging and spectroscopic capabilities. THESEUS will be perfectly suited for addressing the main open issues in cosmology such as, e.g., star formation rate and metallicity evolution of the inter-stellar and intra-galactic medium up to redshift $\sim$10, signatures of Pop III stars, sources and physics of re-ionization, and the faint end of the galaxy luminosity function. In addition, it will provide unprecedented capability to monitor the X-ray variable sky, thus detecting, localizing, and identifying the electromagnetic counterparts to sources of gravitational radiation, which may be routinely detected in the late '20s / early '30s by next generation facilities like aLIGO/ aVirgo, eLISA, KAGRA, and Einstein Telescope. THESEUS will also provide powerful synergies with the next generation of multi-wavelength observatories (e.g., LSST, ELT, SKA, CTA, ATHENA).

The gravitational waves from the first binary neutron star merger, GW170817, were accompanied by a multi-wavelength electromagnetic counterpart, from $\gamma$-rays to radio. The accompanying gamma-rays, seems at first to confirm the association of mergers with short gamma-ray bursts (sGRBs). The common interpretation was that we see an emission from an sGRB jet seen off-axis. However, a closer examination of the sub-luminous $\gamma$-rays and the peculiar radio afterglow were inconsistent with this simple interpretation. Here we present results of 3D and 2D numerical simulations that follow the hydrodynamics and emission of the outflow from a neutron star merger form its ejection and up to its deceleration by the circum-merger medium. Our results show that the entire set of $\gamma$-rays, X-rays and radio observations can be explained by the emission from a mildly relativistic cocoon material (Lorentz factor $\sim$2-5) that was formed while a jet propagated through the material ejected during the merger. The $\gamma$-rays are generated when the cocoon breaks out from the engulfing ejecta while the afterglow is produced by interaction of the cocoon matter with the interstellar medium. The strong early uv/optical signal may be a Lorentz boosted macronova/kilonova. The fate of the jet itself is currently unknown, but our full-EM models define a path to resolving between successful and choked jet scenarios, outputting coupled predictions for the image size, morphology, observed time-dependent polarization and light curve behavior from radio to X-ray. The predictive power of these models will prove key in interpreting the on-going multi-faceted observations of this unprecedented event.

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

We find, using high resolution numerical relativistic simulations, that the tail of the dynamical ejecta of neutron star mergers extends to mildly relativistic velocities faster than $0.7c$. The kinetic energy of this fast tail is $\sim 10^{47}$--$10^{49}$ erg, depending on the neutron star equation of state and on the binary masses. The synchrotron flare arising from the interaction of this fast tail with the surrounding ISM can power the observed non-thermal emission that followed GW170817, provided that the ISM density is $\sim 10^{-2}\,{\rm cm^{-3}}$, the two neutron stars had roughly equal masses and the neutron star equation of state is soft (small neutron star radii). One of the generic predictions of this scenario is that the cooling frequency crosses the X-ray band on a time scale of a few months to a year, leading to a cooling break in the X-ray light curve. If this dynamical ejecta scenario is correct, we expect that the synchrotron radio flare from the ejecta that have produced the macronova/kilonova emission will be observable on time scales of $10^3$ to $10^5$ days. Further multi-frequency observations will confirm or rule out this dynamical ejecta scenario.

Neutron star mergers have been long considered as promising sites of heavy $r$-process nucleosynthesis. We overview observational evidence supporting this scenario including: the total amount of $r$-process elements in the Galaxy, extreme metal poor stars, geological radioactive elemental abundances, dwarf galaxies, and short gamma-ray bursts (sGRBs). Recently, the advanced LIGO and Virgo observatories discovered a gravitational-wave signal of a neutron star merger, GW170817, as well as accompanying multi-wavelength electromagnetic (EM) counterparts. The ultra-violet, optical, and near infrared observations point to $r$-process elements that have been synthesized in the merger ejecta. The rate and ejected mass inferred from GW170817 and the EM counterparts are consistent with other observations. We find however that, within simple one zone chemical evolution models (based on merger rates with reasonable delay time distributions as expected from evolutionary models, or from observations of sGRBs), it is difficult to reconcile the current observations of the europium abundance history of Galactic stars for [Fe/H] $\gtrsim -1$. This implies that to account for the role of mergers in the Galactic chemical evolution, we need a Galactic model with multiple populations that have different spatial distributions and/or varying formation rates.

GW170817 is the first gravitational wave detection of a binary neutron star merger. It was accompanied by radiation across the electromagnetic spectrum and localized to the galaxy NGC 4993 at a distance of 40 Mpc. It has been proposed that the observed gamma-ray, X-ray and radio emission is due to an ultra-relativistic jet launched during the merger, directed away from our line of sight. The presence of such a jet is predicted from models positing neutron star mergers as the central engines driving short-hard gamma-ray bursts (SGRBs). Here we show that the radio light curve of GW170817 has no direct signature of an off-axis jet afterglow. While we cannot rule out the existence of a jet pointing elsewhere, the observed gamma-rays could not have originated from such a jet. Instead, the radio data requires a mildly relativistic wide-angle outflow moving towards us. This outflow could be the high velocity tail of the neutron-rich material dynamically ejected during the merger or a cocoon of material that breaks out when a jet transfers its energy to the dynamical ejecta. The cocoon scenario can explain the radio light curve of GW170817 as well as the gamma-rays and X-rays (possibly also ultraviolet and optical emission), and hence is the model most consistent with the observational data. Cocoons may be a ubiquitous phenomenon produced in neutron star mergers, giving rise to a heretofore unidentified population of radio, ultraviolet, X-ray and gamma-ray transients in the local universe.

We investigate the saturation of particle acceleration in relativistic reconnection using two-dimensional particle-in-cell simulations at various magnetizations \sigma. We find that the particle energy spectrum produced in reconnection quickly saturates as a hard power law that cuts off at \gamma~4\sigma, confirming previous work. Using particle tracing, we find that particle acceleration by the reconnection electric field in X-points determines the shape of the particle energy spectrum. By analyzing the current sheet structure, we show that physical cause of saturation is the spontaneous formation of secondary magnetic islands that can disrupt particle acceleration. By comparing the size of acceleration regions to the typical distance between disruptive islands, we show that the maximum Lorentz factor produced in reconnection is \gamma ~ 5 \sigma, which is very close to what we find in our particle energy spectra. We also show that the dynamic range in Lorentz factor of the power law spectrum in reconnection is < 40. The hardness of the power law combined with its narrow dynamic range implies that relativistic reconnection is capable of producing the hard narrowband flares observed in the Crab Nebula but has difficulty producing the softer broadband prompt GRB emission.

We explore the possibility that the recently detected dipole anisotropy in the arrival directions of~$>8$~EeV ultra-high energy cosmic-rays (UHECRs) arises due to the large-scale structure (LSS). We assume that the cosmic ray sources follow the matter distribution and calculate the flux-weighted UHECRs' RMS dipole amplitude taking into account the diffusive transport in the intergalactic magnetic field (IGMF). We find that the flux-weighted RMS dipole amplitude is $\sim8$% before entering the Galaxy. The amplitude in the [4-8] EeV is only slightly lower $\sim 5$%. The required IGMF is of the order of {5-30 nG}, and the UHECR sources must be relatively nearby, within $\sim$300 Mpc. The absence of statistically significant signal in the lower energy bin can be explained if the same nuclei specie dominates the composition in both energy bins and diffusion in the Galactic magnetic field (GMF) reduces the dipole of these lower rigidity particles. Photodisintegration of higher energy UHECRs could also reduce somewhat the lower energy dipole.

We develop a fully three dimensional numerical code describing the diffusion of cosmic rays in the Milky Way. It includes the nuclear spallation chain up to Oxygen, and allows the study of various cosmic ray properties, such as the CR age, grammage traversed, and the ratio between secondary and primary particles. This code enables us to explore a model in which a large fraction of the cosmic ray acceleration takes place in the vicinity of galactic spiral arms and that these spiral arms are dynamic. We show that the effect of having dynamic spiral arms is to limit the age of cosmic rays at low energies. This is because at low energies the time since the last spiral arm passage governs the Cosmic Ray (CR) age, and not diffusion. Using the model, the observed spectral dependence of the secondary to primary ratio is recovered without requiring any further assumptions such as a galactic wind, re-acceleration or various assumptions on the diffusivity. In particular, we obtain a secondary to primary ratio which increases with energy below about 1 GeV.

Electron capture (EC) isotopes are known to provide constraints on the low energy behavior of cosmic rays (CRs), such as re-acceleration. Here we study the EC isotopes within the framework of the dynamic spiral-arms CR propagation model in which most of the CR sources reside in the galactic spiral arms. The model was previously used to explain the B/C and sub-Fe/Fe ratios \citep{BoverC,Iron}. We show that the known inconsistency between the $^{49}$Ti/$^{49}$V and $^{51}$V/$^{51}$Cr ratios remains also in the spiral-arms model. On the other hand, unlike the general wisdom in which the isotope ratios depend primarily on reacceleration, we find here that the ratio also depends on the halo size ($z_\mathrm{h}$) and in spiral-arms models also on the time since the last spiral arm passage ($\tau_\mathrm{arm}$). Namely, EC isotopes can in principle provide interesting constraints on the diffusion geometry. However, with the present uncertainties in the lab measurements of both the electron attachment rate and the fragmentation cross-sections, no meaningful constraint can be placed.

Short Gamma-Ray Bursts (SGRBs) are believed to arise from compact binary mergers (either neutron star-neutron star or black hole-neutron star). If so their jets must penetrate outflows that are ejected during the merger. As a jet crosses the ejecta it dissipates its energy, producing a hot cocoon which surrounds it. We present here 3D numerical simulations of jet propagation in mergers' outflows and we calculate the resulting emission. This emission consists of two components: the cooling emission, the leakage of the thermal energy of the hot cocoon, and the cocoon macronova that arises from the radioactive decay of the cocoon's material. This emission gives a brief (~ one hour) blue, wide angle signal. While the parameters of the outflow and jet are uncertain, for the configurations we have considered the signal is bright (~ -14 $-$ -15 absolute magnitude) and outshines all other predicted UV-optical signals. The signal is brighter when the jet breakout time is longer and its peak brightness does not depend strongly on the highly uncertain opacity. A rapid search for such a signal is a promising strategy to detect an electromagnetic merger counterpart. A detected candidate could be then followed by deep IR searches for the longer but weaker macronova arising from the rest of the ejecta.

Merging neutron stars offer an exquisite laboratory for simultaneously studying strong-field gravity and matter in extreme environments. We establish the physical association of an electromagnetic counterpart EM170817 to gravitational waves (GW170817) detected from merging neutron stars. By synthesizing a panchromatic dataset, we demonstrate that merging neutron stars are a long-sought production site forging heavy elements by r-process nucleosynthesis. The weak gamma-rays seen in EM170817 are dissimilar to classical short gamma-ray bursts with ultra-relativistic jets. Instead, we suggest that breakout of a wide-angle, mildly-relativistic cocoon engulfing the jet elegantly explains the low-luminosity gamma-rays, the high-luminosity ultraviolet-optical-infrared and the delayed radio/X-ray emission. We posit that all merging neutron stars may lead to a wide-angle cocoon breakout; sometimes accompanied by a successful jet and sometimes a choked jet.

We present an implementation of the Gehrels et al. (2016) galaxy-targeted strategy for gravitational-wave (GW) follow-up using the Las Cumbres Observatory global network of telescopes. We use the Galaxy List for the Advanced Detector Era (GLADE) galaxy catalog, which we show is complete (with respect to a Schechter function) out to ~300 Mpc for galaxies brighter than the median Schechter function galaxy luminosity. We use a prioritization algorithm to select the galaxies with the highest chance of containing the counterpart given their luminosity, their position, and their distance relative to a GW localization, and in which we are most likely to detect a counterpart given its expected brightness compared to the limiting magnitude of our telescopes. This algorithm can be easily adapted to any expected transient parameters and telescopes. We implemented this strategy during the second Advanced Detector Observing Run (O2) and followed the black hole merger GW170814 and the neutron star merger GW170817. For the latter, we identified an optical kilonova/macronova counterpart thanks to our algorithm selecting the correct host galaxy fifth in its ranked list among 182 galaxies we identified in the Laser Interferometer Gravitational-wave Observatory LIGO-Virgo localization. This also allowed us to obtain some of the earliest observations of the first optical transient ever triggered by a GW detection (as presented in a companion paper).

The three gravitational wave events detected by LIGO are opening a new era for high- energy astrophysics. Nevertheless, location of such events remain unknown. A promising solution to the localization problem is to find an electromagnetic (EM) counterpart of GW- generating events, such as binary neutron star mergers (BNS). Indeed, their GW emission will be above sensitivity threshold in the near future. BNS are also considered as short Gamma Ray Bursts (sGRB) progenitors. However, sGRB are highly beamed. In this study, we will therefore focus on another EM counterpart candidate: the cocoon afterglow. The propagation of the GRB jet inside the matter ejected by the BNS produces a cocoon. Then, similarly to the GRB afterglow, a cocoon afterglow is produced, but with a mildly relativistic velocity. Firstly, we propose a model that gives the full hydrodynamic evolution of the cocoon including the mildly relativistic regime. Then we calculate the cocoon afterglow emission in X-ray, optical and radio wavelengths. Finally, we compare the cocoon afterglow emission to the GRB afterglow emission and conclude that the cocoon afterglow is a promising EM counterpart.

The merger of two neutron stars is predicted to give rise to three major detectable phenomena: a short burst of gamma-rays, a gravitational wave signal, and a transient optical/near-infrared source powered by the synthesis of large amounts of very heavy elements via rapid neutron capture (the r-process). Such transients, named "macronovae" or "kilonovae", are believed to be centres of production of rare elements such as gold and platinum. The most compelling evidence so far for a kilonova was a very faint near-infrared rebrightening in the afterglow of a short gamma-ray burst at z = 0.356, although findings indicating bluer events have been reported. Here we report the spectral identification and describe the physical properties of a bright kilonova associated with the gravitational wave source GW 170817 and gamma-ray burst GRB 170817A associated with a galaxy at a distance of 40 Mpc from Earth. Using a series of spectra from ground-based observatories covering the wavelength range from the ultraviolet to the near-infrared, we find that the kilonova is characterized by rapidly expanding ejecta with spectral features similar to those predicted by current models. The ejecta is optically thick early on, with a velocity of about 0.2 times light speed, and reaches a radius of about 50 astronomical units in only 1.5 days. As the ejecta expands, broad absorption-like lines appear on the spectral continuum indicating atomic species produced by nucleosynthesis that occurs in the post-merger fast-moving dynamical ejecta and in two slower (0.05 times light speed) wind regions. Comparison with spectral models suggests that the merger ejected 0.03-0.05 solar masses of material, including high-opacity lanthanides.

We present the first relativistic MHD numerical simulation of a magnetic jet that propagates through and emerges from the dynamical ejecta of a binary neutron star merger. Generated by the magnetized rotation of the merger remnant, the jet propagates through the ejecta and produces an energetic cocoon that expands at mildly relativistic velocities and breaks out of the ejecta. We show that if the ejecta has a low-mass ($\sim10^{-7} M_\odot$) high-velocity ($v\sim0.85$ c) tail, the cocoon shock breakout will generate $\gamma$-ray emission that is comparable to the observed short GRB170817A that accompanied the recent gravitational wave event GW170817. Thus, we propose that this GRB, which is quite different from all other short GRBs observed before, was produced by a different mechanism. We expect, however, that such events are numerous and many will be detected in coming LIGO-Virgo runs.

Gravitational waves have been detected from a binary neutron star merger event, GW170817. The detection of electromagnetic radiation from the same source has shown that the merger occurred in the outskirts of the galaxy NGC 4993, at a distance of 40 megaparsecs from Earth. We report the detection of a counterpart radio source that appears 16 days after the event, allowing us to diagnose the energetics and environment of the merger. The observed radio emission can be explained by either a collimated ultra-relativistic jet viewed off-axis, or a cocoon of mildly relativistic ejecta. Within 100 days of the merger, the radio light curves will distinguish between these models and very long baseline interferometry will have the capability to directly measure the angular velocity and geometry of the debris.

The Telescope Array (TA) shows a 20$^{\circ}$ hotspot as well as an excess of UHECRs above 50~EeV when compared with the Auger spectrum. We consider the possibility that both the TA excess and hotspot are due to a dominant source in the Northern sky. We carry out detailed simulations of UHECR propagation in both the intergalactic medium and the Galaxy, using different values for the intergalactic magnetic field. We consider two general classes of sources: transients and steady, adopting a mixed UHECR composition that is consistent with the one found by Auger. The spatial location of the sources is draw randomly. We generate Auger-like and TA-like data sets from which we determine the spectrum, the sky maps and the level of anisotropy. We find that, while steady sources are favored over transients, it is unlikely to account for all the currently available observational data. While we reproduce fairly well the Auger spectrum for the vast majority of the simulated data sets, most of the simulated data sets with a spectrum compatible with that of TA (at most a few percent depending on density model tested) show a much stronger anisotropy than the one observed. We find that the rare cases in which both the spectrum and the anisotropy are consistent require a steady source within $\sim 10$ Mpc, to account for the flux excess, and a strong extragalactic magnetic field $\sim 10$ nG, to reduce the excessive anisotropy.

One of the puzzles in the recent observations of gravitational waves from binary black hole mergers is the observed low (projected) spins of the progenitor black holes. In two of the four events, GW150914, and the recent GW170104, the observed spins are most likely negative (but consistent with zero). In the third case LVT151012 it is practically zero and only in the forth case, GW151226, the spin is positive but low. These observations are puzzling within the field binary scenario in which positive higher spins are expected. Considering the most favorable Wolfe Rayet (WR) progenitors we estimate the expected spin distribution for different evolution scenarios and compare it to the observations. With typical parameters one expects a significant fraction ($\ge 25\%$) of the mergers to have high effective spin values. However due to uncertainties in the outcome of the common envelope phase (typical separation and whether the stars are rotating or not) and in the late stages of massive star evolution (the strength of the winds) one cannot rule our scenarios in which the expected spins are low. While observations of high effective spin events will support this scenario, further observations of negative spin events would rule it out.

GeV-TeV gamma-ray and PeV-EeV neutrino backgrounds provide a unique window on the nature of the ultra-high-energy cosmic-rays (UHECRs). We discuss the implications of the recent Fermi-LAT data regarding the extragalactic gamma-ray background (EGB) and related estimates of the contribution of point sources as well as IceCube neutrino data on the origin of the UHECRs. We calculate the diffuse flux of cosmogenic $\gamma$-rays and neutrinos produced during the UHECRs propagation and derive constraints on the possible cosmological evolution of UHECR sources. In particular, we show that the mixed-composition scenario which is in agreement with both (i) Auger measurements of the energy spectrum and composition up to the highest energies and (ii) the ankle-like feature in the light component detected by KASCADE-Grande, is compatible with both the Fermi-LAT measurements and with current IceCube limits.

After several decades of extensive research the mechanism driving core-collapse supernovae (CCSNe) is still unclear. A common mechanism is a neutrino driven outflow, but others have been proposed. Among those, a long-standing idea is that jets play an important role in SN explosions. Gamma-ray bursts (GRBs) that accompany rare and powerful CCSNe, sometimes called "hypernovae", provide a clear evidence for a jet activity. The relativistic GRB jet punches a hole in the stellar envelope and produces the observed gamma-rays far outside the progenitor star. While SNe and jets coexist in long GRBs, the relation between the mechanisms driving the hypernova and the jet is unknown. Also unclear is the relation between the rare hypernovae and the more common CCSNe. Here we {present observational evidence that indicates} that choked jets are active in CCSNe types that are not associated with GRBs. A choked jet deposits all its energy in a cocoon. The cocoon eventually breaks out from the star releasing energetic material at very high, yet sub-relativistic, velocities. This fast moving material has a unique signature that can be detected in early time SN spectra. We find a clear evidence for this signature in several CCSNe, all involving progenitors that have lost all, or most, of their hydrogen envelope prior to the explosion. These include CCSNe that don't harbor GRBs or any other relativistic outflows. Our findings suggest a continuum of central engine activities in different types of CCSNe and call for rethinking of the explosion mechanism of regular CCSNe

The plateau in the duration distribution of long Gamma-Ray Bursts (LGRBs) provides a direct observational evidence for the Collapsar model. The plateau reflects the fact that the observed duration satisfies: $T_{90} = t_{e}-t_{b}$ where $t_{e}$ is the time that the central engine operates and $t_{b}$ is a threshold time, interpreted within the Collapsar model as the time it takes for the relativistic jet to penetrate the stellar envelope. Numerical simulation and macronova observations suggest that compact binary mergers involve mass ejection. If short-Gamma Ray Bursts (sGRBs) arise from such mergers, their jets should cross this surrounding ejecta before producing the prompt emission. Like in LGRBs, this should result in a distinct short plateau in the GRBs' duration distribution. We present a new analysis of the duration distribution for the three GRB satellites: BATSE, {\it Swift} and Fermi. We find a clear evidence for a short ($\sim 0.4$ sec) plateau in the duration distribution. This plateau is consistent with the expected jet crossing time, provided that the ejecta is of order of a few percent of solar masses.

We explore the implications of the observed low spins aligned with the orbital axis in Advanced~LIGO~O1 run on binary black hole (BBH)~merger scenarios in which the merging BBHs have evolved from field binaries. The coalescence time determines the initial orbital separation of BBHs. This, in turn, determines whether the stars are synchronized before collapse and hence determines their projected spins. Short coalescence times imply synchronization and large spins. Among known stellar objects, Wolf-Rayet~(WR) stars seem the only progenitors consistent with the low aligned spins observed in LIGO's~O1 provided that the orbital axis maintains its direction during the collapse. We calculate the spin distribution of BBH mergers in the local Universe and its redshift evolution for WR progenitors. Assuming that the BBH formation rate peaks around a redshift of $\sim 2\,$--$\,3$, we show that BBH mergers in the local Universe are dominated by low spin events. The high spin population starts to dominate at a redshift of $\sim0.5\,$--$\,1.5$. WR stars are also progenitors of long Gamma-Ray~Bursts~(LGRBs) that take place at a comparable rate to BBH mergers. We discuss the possible connection between the two phenomena. Additionally, we show that hypothetical Population~III star progenitors are also possible. Although WR and Population~III progenitors are consistent with the current data, both models predict a non-vanishing fraction of high aligned spin BBH mergers. If those are not detected within the coming LIGO/Virgo runs, it will be unlikely that the observed BBHs formed via field binaries.

GeV-TeV gamma-rays and PeV-EeV neutrino backgrounds provide a unique window on the nature of the ultra-high-energy cosmic-rays (UHECRs). We discuss the implications of the recent Fermi-LAT data regarding the extragalactic gamma-ray background (EGB) and related estimates of the contribution of point sources as well as IceCube neutrino data on the origin of the UHECRs. We calculate the diffuse flux of cosmogenic $\gamma$-rays and neutrinos produced by the UHECRs and derive constraints on the possible cosmological evolution of UHECR sources. In particular, we show that the mixed-composition scenario considered in \citet{Globus2015b}, which is in agreement with both (i) Auger measurements of the energy spectrum and composition up to the highest energies and (ii) the ankle-like feature in the light component detected by KASCADE-Grande, is compatible with both the Fermi-LAT measurements and with current IceCube limits. We also discuss the possibility for future experiments to detect associated cosmogenic neutrinos and further constrain the UHECR models, including possible subdominant UHECR proton sources.

The Telescope Array (TA) shows a 20$^{\circ}$ hotspot as well as an excess of UHECRs above 50~EeV when compared with the Auger spectrum. We consider the possibility that both the TA excess and hotspot are due to a dominant source in the Northern sky. We carry out detailed simulations of UHECR propagation in both the intergalactic medium and the Galaxy, using different values for the intergalactic magnetic field. We consider two general classes of sources: transients and steady, adopting a mixed UHECR composition that is consistent with the one found by Auger. The spatial location of the sources is draw randomly. We generate Auger-like and TA-like data sets from which we determine the spectrum, the sky maps and the level of anisotropy. We find that, while steady sources are favored over transients, it is unlikely to account for all the currently available observational data. While we reproduce fairly well the Auger spectrum for the vast majority of the simulated data sets, most of the simulated data sets with a spectrum compatible with that of TA (at most a few percent depending on density model tested) show a much stronger anisotropy than the one observed. We find that the rare cases in which both the spectrum and the anisotropy are consistent require a steady source within $\sim 10$ Mpc, to account for the flux excess, and a strong extragalactic magnetic field $\sim 10$ nG, to reduce the excessive anisotropy.

Models for tidal disruption events (TDEs) in which a supermassive black hole disrupts a star commonly assume that the highly eccentric streams of bound stellar debris promptly form a circular accretion disk at the pericenter scale. However, recent numerical simulations (Shiokawa et al., 2015) demonstrated that dissipation via hydrodynamical shocks is insufficient to circularize debris, and the flow retains its initial semi-major axis scale throughout the first ~10 orbits of the event. The bolometric peak luminosity of most TDE candidates, a few x 10^{44} erg/s, implies that we observe only ~1% of the energy expected from radiatively efficient accretion. Motivated by these results, (Piran et al., 2015) suggested that the observed optical TDE emission is powered by shocks at the apocenter between freshly infalling material and earlier arriving matter. This model explains the small radiated energy, the low temperature, and the large radius implied by the observations as well as the t^{-5/3} light curve. However the question of the system's low {bolometric} efficiency remains unanswered. We suggest that the high orbital energy and low angular momentum of the flow make it possible for magnetic stresses to reduce the matter's already small angular momentum to the point at which it can fall ballistically into the SMBH before circularization. As a result, the efficiency is only ~1--10% of a standard accretion disk's efficiency. Thus, the intrinsically high eccentricity of the tidal debris naturally explains why most TDE candidates are fainter than expected.

Macronovae (kilonovae) that arise in binary neutron star mergers are powered by radioactive beta decay of hundreds of $r$-process nuclides. We derive, using Fermi's theory of beta decay, an analytic estimate of the nuclear heating rate. We show that the heating rate evolves as a power law ranging between $t^{-6/5}$ to $t^{-4/3}$. The overall magnitude of the heating rate is determined by the mean values of nuclear quantities, e.g., the nuclear matrix elements of beta decay. These values are specified by using nuclear experimental data. We discuss the role of higher order beta transitions and the robustness of the power law. The robust and simple form of the heating rate suggests that observations of the late-time bolometric light curve $\propto t^{-\frac{4}{3}}$ would be a direct evidence of a $r$-process driven macronova. Such observations could also enable us to estimate the total amount of $r$-process nuclei produced in the merger.

Most of the diffuse Galactic GeV gamma-ray emission is produced via collisions of cosmic ray (CR) protons with ISM protons. As such the observed spectra of the gamma-rays and the CRs should be strongly linked. Recent observations of Fermi-LAT exhibit a hardening of the gamma-ray spectrum at around a hundred GeV, between the Sagittarius and Carina tangents, and a further hardening at a few degrees above and below the Galactic plane. However, standard CR propagation models that assume a time independent source distribution and a location independent diffusion cannot give rise to a spatially dependent CR (and hence gamma-ray) spectral slopes. Here we consider a dynamic spiral arm model in which the distribution of CR sources is concentrated in the (dynamic) spiral arms, and we study the effects of this model on the $\pi^0$-decay produced gamma-ray spectra. Within this model, near the Galactic arms the observed gamma-ray spectral slope is not trivially related to the CR injection spectrum and energy dependence of the diffusion coefficient. We find unique signatures that agree with the Fermi-LAT observations. This model also provides a physical explanation for the difference between the local CR spectral slope and the CR slope inferred from the average gamma-ray spectrum.

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

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

It is largely recognized that Gamma-Ray Burst (GRB) jets involve ultra-relativistic motion. However, the value of the Lorentz factor Gamma_0 is still not clear and only lower limits are known for most bursts. We suggest here a new method to obtain upper limits on Gamma_0. The early high-energy synchrotron afterglow flux depends strongly on Gamma_0. Upper limits on GeV emission therefore provide uppers limit on Gamma_0. Applying this method to 190 Fermi GRBs that have not been detected by the Fermi-LAT we place upper limits on the high-energy afterglow flux, and in turn on Gamma_0. For bursts at a typical redshift z=2, we find values of the order of 200 (and above) for a homogeneous density medium, and in the range 100-400 for a wind-like medium. These upper limits are consistent with (and are very close to) lower limits and direct estimates inferred using other methods, suggesting that the typical Lorentz factors of GRB jets are of order a few hundred.

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 $\gamma$-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.

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.

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 $\gamma$-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.

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

Supernovae near the galactic center evolve differently from regular galactic supernovae. This is mainly due to the environment into which the supernova remnants propagate. Instead of a static, uniform density medium, SNRs near the galactic center propagate into a wind swept environment with a velocity away from the galactic center, and a graded density profile. This causes these SNRs to be non - spherical, and to evolve faster than their galactic counterparts. We develop an analytic theory for the evolution of explosions within a stellar wind, and verify it using a hydrodynamic code. We show that such explosions can evolve in one of three possible morphologies. Using these results we discuss the association between the two SNRs (SGR East and SGR A's bipolar radio/X-ray Lobes) and the two neutron stars (the cannonball and SGR J1745-2900) near the galactic center. We show that, given the morphologies of the SNR and positions of the neutron stars, the only possible association is between SGR A's bipolar radio/X-ray Lobes and SGR J1745-2900. The compact object created in the explosion of SGR East remains undetected, and the SNR of the supernova that created the cannonball has already disappeared.

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

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.

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

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

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

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.

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.

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

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

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

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.

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.

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

Measurements of the radioactive $^{244}$Pu abundances can break the degeneracy between high-rate/low-yield and low-rate/high-yield scenarios for the production of heavy $r$-process elements. The first corresponds to production by core collapse supernovae (cc-SNe) while the latter corresponds to production by e.g. compact binary mergers. The estimated $^{244}$Pu abundance in the current interstellar medium inferred from deep-sea measurements (Wallner et al. 2015) is significantly lower than that corresponding Early Solar System abundances (Turner et al 2007). We estimate the expected median value of the $^{244}$Pu abundances and fluctuations around this value in both models. We show that while the current and Early Solar System abundances are naturally explained within the low-rate/high-yield (e.g. merger) scenario, they are incompatible with the high-rate/low-yield (cc-SNe) model. The inferred event rate remarkably agrees with compact binary merger rates estimated from Galactic neutron star binaries and from short gamma-ray bursts. Furthermore, the ejected mass of $r$-process elements per event agrees with both theoretical and observational macronova/kilonova estimates.

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

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

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

Evidence has been accumulated on the existence of a thermal-like component during the prompt phase of GRBs. This component, often associated with the GRB jet's photosphere, is usually subdominant compared to a much stronger non-thermal one. The prompt emission of Fermi GRB 131014A provides a unique opportunity to study this thermal-like component. Indeed, the thermal emission in GRB 131014A is much more intense than in other GRBs and a pure thermal episode is observed during the initial 0.16 s. The thermal-like component cools monotonically during the first second while the non-thermal emission kicks off. The intensity of the non-thermal component progressively increases until being energetically dominant at late time. This is a perfect scenario to disentangle the thermal component from the non-thermal one. A low-energy spectral index of +0.6 better fit the thermal component than the typical index value +1 corresponding to a pure Planck function. The non-thermal component is adequately fitted with a Band function whose low and high energy power law indices are ~-0.7 and <~-3, respectively; this is also statistically equivalent to a cutoff power law with a ~-0.7 index. This is in agreement with our previous results. Finally, a strong correlation is observed between the time-resolved luminosity of the non-thermal component, L$_i^{nTh}$, and its corresponding rest frame spectral peak energy, E$_{peak,i}^{rest,nTh}$, with a slope similar to the one reported in our previous articles. Assuming this relation to be universal for all GRBs we estimate a redshift of ~1.55 for GRB 131014A that is a typical value for long GRBs. These observational results are consistent with the models in which the non-thermal emission is produced well above the GRB jet photosphere but they may also be compatible with other scenarios (e.g., dissipative photosphere) that are not discussed in this article.

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

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

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.

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

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

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

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

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.

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

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

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

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

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.

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

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

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

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

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

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

We follow the longterm evolution of the dynamic ejecta of neutron star mergers for up to 100 years and over a density range of roughly 40 orders of magnitude. We include the nuclear energy input from the freshly synthesized, radioactively decaying nuclei in our simulations and study its effects on the remnant dynamics. Although the nuclear heating substantially alters the longterm evolution, we find that running nuclear networks over purely hydrodynamic simulations (i.e. without heating) yields actually acceptable nucleosynthesis results. The main dynamic effect of the radioactive heating is to quickly smooth out inhomogeneities in the initial mass distribution, subsequently the evolution proceeds self-similarly and after 100 years the remnant still carries the memory of the initial binary mass ratio. We also explore the nucleosynthetic yields for two mass ejection channels. The dynamic ejecta very robustly produce "strong" r-process elements with $A > 130$ with a pattern that is essentially independent of the details of the merging system. From a simple model we find that neutrino-driven winds yield "weak" r-process contributions with $50 < A < 130$ whose abundance patterns vary substantially between different merger cases. This is because their electron fraction, set by the ratio of neutrino luminosities, varies considerably from case to case. Such winds do not produce any $^{56}{\rm Ni}$, but a range of radioactive isotopes that are long-lived enough to produce a second, radioactively powered electromagnetic transient in addition to the "macronova" from the dynamic ejecta. While our wind model is very simple, it nevertheless demonstrates the potential of such neutrino-driven winds for electromagnetic transients and it motivates further, more detailed neutrino-hydrodynamic studies. The properties of the mentioned transients are discussed in more detail in a companion paper.

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

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

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

Many of the astrophysical sources and violent phenomena observed in our Universe are potential emitters of gravitational waves (GW) and high-energy neutrinos (HEN). Both GWs and HENs may escape very dense media and travel unaffected over cosmological distances, carrying information from the innermost regions of the astrophysical engines. Such messengers could also reveal new, hidden sources that have not been observed by conventional photon-based astronomy. Coincident observation of GWs and HENs may thus play a critical role in multimessenger astronomy. This is particularly true at the present time owing to the advent of a new generation of dedicated detectors: IceCube, ANTARES, VIRGO and LIGO. Given the complexity of the instruments, a successful joint analysis of this data set will be possible only if the expertise and knowledge of the data is shared between the two communities. This review aims at providing an overview of both theoretical and experimental state-of-the-art and perspectives for such a GW+HEN multimessenger astronomy.

Gamma-ray observations of a stellar tidal disruption event (TDE) detected by the Swift satellite and follow up observations in radio, mm, infrared and x-ray bands have provided a rich data set to study accretion onto massive blackholes, production of relativistic jets and their interaction with the surrounding medium. The radio and x-ray data for TDE Swift J1644+57 provide a conflicting picture regarding the energy in relativistic jet produced in this event: x-ray data suggest jet energy declining with time as t^{-5/3} whereas the nearly flat lightcurves in radio and mm bands lasting for about 100 days have been interpreted as evidence for the total energy output increasing with time. We show in this work that flat lightcurves do not require addition of energy to decelerating external shock (which produced radio and mm emission via synchrotron process), instead the flat behavior is due to inverse-Compton cooling of electrons by x-ray photons streaming through the external shock; the higher x-ray flux at earlier times cools electrons more rapidly thereby reducing the emergent synchrotron flux, and this effect weakens as the x-ray flux declines with time.

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.

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

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

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

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

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

We explore the multi-messenger signatures of encounters between two neutron stars and between a neutron star and a stellar-mass black hole. We focus on the differences between gravitational wave driven binary mergers and dynamical collisions that occur, for example, in globular clusters. For both types of encounters we compare the gravitational wave and neutrino emission properties. We also calculate fallback rates and analyze the properties of the dynamically ejected matter. Last but not least we address the electromagnetic transients that accompany each type of encounter. The canonical nsns merger case ejects more than 1% of a solar mass of extremely neutron-rich ($Y_e\sim 0.03$) material, an amount that is consistent with double neutron star mergers being a major source of r-process in the galaxy. nsbh collisions eject very large amounts of matter ($\sim 0.15$ \msun) which seriously constrains their admissible occurrence rates. The compact object collision rate must therefore be less, likely much less, than 10% of the nsns merger rate. The radioactively decaying ejecta produce optical-UV "macronova" which, for the canonical merger case, peak after $\sim 0.4$ days with a luminosity of $\sim 10^{42}$ erg/s. nsns (nsbh) collisions reach up to 3 (7) times larger peak luminosities. The dynamic ejecta deposit a kinetic energy comparable to a supernova in the ambient medium. The canonical merger case releases approximately $2 \times 10^{50}$ erg, the most extreme (but likely rare) cases deposit kinetic energies of up to $10^{52}$ erg. The deceleration of this mildly relativistic material by the ambient medium produces long lasting radio flares. A canonical ns$^2$ merger at the detection horizon of advanced LIGO/Virgo produces a radio flare that peaks on a time scale of one year with a flux of $\sim$0.1 mJy at 1.4 GHz.

Compact binary mergers are prime sources of gravitational waves, targeted by current and next generation detectors. 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. We present a large set of numerical simulations that focus on the electromagnetic signals that emerge from the dynamically ejected sub-relativistic material. These outflows produce on a time scale of a day macronovae - short-lived IR to UV signals powered by radioactive decay. The interaction of this outflow with the surrounding matter inevitably leads to a long-lasting remnant. The expected radio signals of these remnants last longer than a year, when the sub-relativistic ejecta dominate the emission. We discuss their detectability in 1.4 GHz and 150 MHz and compare it with the detectability of short GRBs' orphan afterglows (which are produced by a different component of this outflow). Mergers with characteristics similar to those of the Galactic neutron star binary population (similar masses and typical circum-merger Galactic disk density of ~1cm^-3) taking place at the detection horizon of advanced GW detectors (300 Mpc) yield 1.4 GHz [150 MHz] signals of ~50 [300] microJy. The signal on time scales of weeks, is dominated by the mildly and/or ultra-relativistic outflow, which is not accounted for by our simulations, and is expected to be even brighter. Upcoming all sky surveys are expected to detect a few dozen, and possibly more, merger remnants at any given time thereby providing robust lower limits to the mergers rate even before the advanced GW detectorsbecome operational. Macronovae from the same distance peak in the IR to UV range at an observed magnitude that may be as bright as 22-23 about 10 hours after the merger but dimmer, redder and longer if the opacity is larger.

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

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

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

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

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.

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.

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

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

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

The recent discovery of the "weak field, old magnetar", the soft gamma repeater SGR 0418+5729, whose dipole magnetic field is less than 7.5 \times 10^{12} G, has raised perplexing questions: How can the neutron star produce SGR-like bursts with such a low magnetic field? What powers the observed X-ray emission when neither the rotational energy nor the magnetic dipole energy are sufficient? These observations, that suggest either a much larger energy reservoir or a much younger true age (or both), have renewed the interest in the evolutionary sequence of magnetars. We examine, here, a phenomenological model for the magnetic field decay: B_dip} \propto (B_dip)^{1+a} and compare its predictions with the observed period, P,the period derivative, \dot{P}, and the X-ray luminosity, L_X, of magnetar candidates. We find a strong evidence for a dipole field decay on a timescale of \sim 10^3 yr for the strongest (\sim 10^{15} G) field objects, with a decay index within the range 1 \leq a < 2 and more likely within 1.5\lesssim a \lesssim 1.8. The decaying field implies a younger age than what is implied by the spinown age. Surprisingly, even with the younger age, the energy released in the dipole field decay is insufficient to power the X-ray emission, suggesting the existence of a stronger internal field, B_int. Examining several models for the internal magnetic field decay we find that it must have a very large (> 10^{16} G) initial value. Our findings suggest two clear distinct evolutionary tracks -- the SGR/AXP branch and the transient branch, with a possible third branch involving high-field radio pulsars that age into low luminosity X-ray dim isolated neutron stars.

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

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

According the Collapsar model long gamma-ray bursts (LGRBs) involve relativistic jets that puncture the envelope of a collapsing star, and produced the \gamma-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.

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.

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

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

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.

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.

The new and extreme population of GRBs detected by Fermi-LAT shows several new features in high energy gamma-rays that are providing interesting and unexpected clues into GRB prompt and afterglow emission mechanisms. Over the last 6 years, it has been Swift that has provided the robust dataset of UV/optical and X-ray afterglow observations that opened many windows into components of GRB emission structure. The relationship between the LAT detected GRBs and the well studied, fainter, less energetic GRBs detected by Swift-BAT is only beginning to be explored by multi-wavelength studies. We explore the large sample of GRBs detected by BAT only, BAT and Fermi-GBM, and GBM and LAT, focusing on these samples separately in order to search for statistically significant differences between the populations, using only those GRBs with measured redshifts in order to physically characterize these objects. We disentangle which differences are instrumental selection effects versus intrinsic properties, in order to better understand the nature of the special characteristics of the LAT bursts.

The Large Area Telescope (LAT) on board of the Fermi satellite detected emission above 30 MeV only in a small fraction of the long gamma-ray bursts (GRBs) detected by the Fermi Gamma-ray Burst Monitor (GBM) at 8 keV - 10 MeV. Those bursts that were detected by the LAT were among the brightest GBM bursts. We examine a sample of the most luminous GBM bursts with no LAT detection and obtain upper limits on their high energy fluence. We find an average upper limit of LAT/GBM fluence ratio of 0.13 for GeV fluence during $T_{90}$ and an average upper limit ratio of 0.45 for GeV fluence during the first 600 seconds after the trigger. These ratios strongly constrain various emission models and in particular rule out SSC models for the prompt emission. In about a third of both LAT detected and LAT non-detected bursts, we find that the extrapolation of the MeV range Band spectrum to the GeV range is larger than the observed GeV fluence (or its upper limit). While this excess is not highly significant for any specific burst, the overall excess in a large fraction of the bursts suggests a decline in the high energy spectral slope in at least some of these bursts. Possibly an evidence for the long sought after pair creation limit.

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

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 $\mu$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.

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

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

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

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.

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

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

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

Gamma-ray bursts (GRBs) are short and intense emission of soft gamma-rays, which have fascinated astronomers and astrophysicists since their unexpected discovery in 1960s. The X-ray/optical/radio afterglow observations confirm the cosmological origin of GRBs, support the fireball model, and imply a long-activity of the central engine. The high energy gamma-ray emission (>20 MeV) from GRBs is particularly important because they shed some lights on the radiation mechanisms and can help us to constrain the physical processes giving rise to the early afterglows. In this work, we review observational and theoretical studies of the high energy emission from GRBs. Special attention is given to the expected high energy emission signatures accompanying the canonical early-time X-ray afterglow that was observed by the Swift X-ray Telescope. We also discuss the detection prospect of the upcoming GLAST satellite and the current ground-based Cerenkov detectors.

Recent measurements of the positron/electron ratio in the cosmic ray (CR) flux exhibits an apparent anomaly, whereby this ratio increases between 10 and 100 GeV. In contrast, this ratio should decrease according to the standard scenario, in which CR positrons are secondaries formed by hadronic interactions between the primary CR protons and the interstellar medium (ISM). The positron excess is therefore interpreted as evidence for either an annihilation/decay of weakly interacting massive particles, or for a direct astrophysical source of pairs. The common feature of all proposed models is that they invoke new physics or new astrophysical sources. However, this line of argumentation relies implicitly on the assumption of a relatively homogeneous CR source distribution. Inhomogeneity of CR sources on a scale of order a kpc, can naturally explain this anomaly. If the nearest major CR source is about a kpc away, then low energy electrons (~GeV) can easily reach us. At higher energies (> 10 GeV), the source electrons cool via synchrotron and inverse-Compton before reaching the solar vicinity. Pairs formed in the local vicinity through the proton/ISM interactions can reach the solar system also at high energies, thus increasing the positron/electron ratio. A natural origin of source inhomogeneity is the strong concentration of supernovae to the galactic spiral arms. Assuming supernova remnants (SNRs) as the sole primary source of CRs, and taking into account their concentration near the galactic spiral arms, we consistently predict the observed positron fraction between 1 and 100 GeV, while abiding to different constraints such as the observed electron spectrum and the CRs cosmogenic age. An ATIC like spectral excess at ~600 GeV can be explained, in this picture, as the contribution of a few known nearby SNRs.

Long-duration gamma-ray bursts (GRBs) are widely believed to be highly-collimated explosions (opening angle theta ~ 1-10 deg). As a result of this beaming factor, the true energy release from a GRB is usually several orders of magnitude smaller than the observed isotropic value. Measuring this opening angle, typically inferred from an achromatic steepening in the afterglow light curve (a "jet" break), has proven exceedingly difficult in the Swift era. Here we undertake a study of five of the brightest (in terms of the isotropic prompt gamma-ray energy release, E(gamma, iso)) GRBs in the Swift era to search for jet breaks and hence constrain the collimation-corrected energy release. We present multi-wavelength (radio through X-ray) observations of GRBs 050820A, 060418, and 080319B, and construct afterglow models to extract the opening angle and beaming-corrected energy release for all three events. Together with results from previous analyses of GRBs 050904 and 070125, we find evidence for an achromatic jet break in all five events, strongly supporting the canonical picture of GRBs as collimated explosions. The most natural explanation for the lack of observed jet breaks from most Swift GRBs is therefore selection effects. However, the opening angles for the events in our sample are larger than would be expected if all GRBs had a canonical energy release of ~ 10e51 erg. The total energy release we measure for those "hyper-energetic" (E(total) >~ 10e52 erg) events in our sample is large enough to start challenging models with a magnetar as the compact central remnant.

Recent measurements of the positron/electron ratio in the cosmic ray (CR) flux exhibits an apparent anomaly, whereby this ratio increases between 10 and 100 GeV. We show that inhomogeneity of CR sources on a scale of order a kpc, can naturally explain this anomaly. If the nearest major CR source is about a kpc away, then low energy electrons ($\sim 1$ GeV) can easily reach us. At higher energies ($\gtrsim 10$ GeV), the source electrons cool via synchrotron and inverse-Compton before reaching Earth. Pairs formed in the local vicinity through the proton/ISM interactions can reach Earth also at high energies, thus increasing the positron/electron ratio. A natural origin of source inhomogeneity is the strong concentration of supernovae in the galactic spiral arms. Assuming supernova remnants (SNRs) as the sole primary source of CRs, and taking into account their concentration near the galactic spiral arms, we consistently recover the observed positron fraction between 1 and 100 GeV. ATIC's electron excess at $\sim 600$ GeV is explained, in this picture, as the contribution of a few known nearby SNRs. The apparent coincident similarity between the cooling time of electrons at 10 GeV (where the positron/electron ratio upturn), $\sim 10$ Myr, and the CRs protons cosmogenic age at the same energy is predicted by this model.

We calculate the high energy (sub-GeV to TeV) prompt and afterglow emission of GRB 080319B that was distinguished by a naked-eye optical flash and by an unusual strong early X-ray afterglow. There are three possible sources for high energy emission: the prompt optical and $\gamma$-ray photons IC scattered by the accelerated electrons, the prompt photons IC scattered by the early external reverse-forward shock electrons, and the higher band of the synchrotron and the synchrotron self-Compton emission of the external shock. There should have been in total {hundreds} high energy photons detectable for the Large Area Telescope (LAT) onboard the Fermi satellite, and {tens} photons of those with energy $> 10$ GeV. The $> 10$ GeV emission had a duration about twice that of the soft $\gamma$-rays. AGILE could have observed these energetic signals if it was not occulted by the Earth at that moment. The physical origins of the high energy emission detected in GRB 080514B, GRB 080916C and GRB 081024B are also discussed. These observations seem to be consistent with the current high energy emission models.

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.

The first structures in the Universe formed at z>7, at higher redshift than all currently known galaxies. Since GRBs are brighter than other cosmological sources at high redshift and exhibit simple power-law afterglow spectra that is ideal for absorption studies, they serve as powerful tools for studying the early universe. New facilities planned for the coming decade will be able to obtain a large sample of high-redshift GRBs. Such a sample would constrain the nature of the first stars, galaxies, and the reionization history of the Universe.

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?

Inverse Compton (IC) scattering is one of two viable mechanisms that can produce the prompt gamma-ray emission in Gamma-Ray Bursts. IC requires low energy seed photons and a population of relativistic electrons that upscatter them. The same electrons can upscatter the gamma-ray photons to even higher energies in the TeV range. Using the current upper limits on the prompt optical emission we show here that under general conservative assumption the IC mechanism suffers from an "energy crisis". Namely, IC will over-produce a very high energy component that would carry much more energy than the observed prompt gamma-rays. Our analysis is general and it makes no assumptions on the specific mechanism that produces the relativistic electrons population.

We examine what the absorbed spectra of extra-galactic TeV gamma-ray sources, such as blazars, would look like in the presence of Lorentz invariance violation (LIV). Pair-production with the extra-galactic background light modifies the observed spectra of such sources, and we show that a violation of Lorentz invariance would generically have a dramatic effect on this absorption feature. Inspecting this effect, an experimental task likely practical in the near future, can provide unique insight on the possibility of LIV.

The extremely bright optical flash that accompanied GRB 080319B suggested, at first glance, that the prompt $\gamma$-rays in this burst were produced by Synchrotron self Compton (SSC). We analyze here the observed optical and $\gamma$ spectrum. We find that the very strong optical emission poses, due to self absorption, very strong constraints on the emission processes and put the origin of the optical emission at a very large radius, almost inconsistent with internal shock. Alternatively it requires a very large random Lorentz factor for the electrons. We find that SSC could not have produced the prompt $\gamma$-rays. We also show that the optical emission and the $\gamma$ rays could not have been produced by synchrotron emission from two populations of electron within the same emitting region. Thus we must conclude that the optical and the $\gamma$-rays were produced in different physical regions. A possible interpretation of the observations is that the $\gamma$-rays arose from internal shocks but the optical flash resulted from external shock emission. This would have been consistent with the few seconds delay observed between the optical and $\gamma$-rays signals.

The only supernovae (SNe) to have shown early gamma-ray or X-ray emission thus far are overenergetic, broad-lined Type Ic SNe (Hypernovae - HNe). Recently, SN 2008D shows several novel features: (i) weak XRF, (ii) an early, narrow optical peak, (iii) disappearance of the broad lines typical of SNIc HNe, (iv) development of He lines as in SNeIb. Detailed analysis shows that SN 2008D was not a normal SN: its explosion energy (KE ~ 6*10^{51} erg) and ejected mass (~7 Msun) are intermediate between normal SNeIbc and HNe. We derive that SN 2008D was originally a ~30Msun star. When it collapsed a black hole formed and a weak, mildly relativistic jet was produced, which caused the XRF. SN 2008D is probably among the weakest explosions that produce relativistic jets. Inner engine activity appears to be present whenever massive stars collapse to black holes.

Long duration gamma-ray bursts (GRBs) release copious amounts of energy across the entire electromagnetic spectrum, and so provide a window into the process of black hole formation from the collapse of a massive star. Over the last forty years, our understanding of the GRB phenomenon has progressed dramatically; nevertheless, fortuitous circumstances occasionally arise that provide access to a regime not yet probed. GRB 080319B presented such an opportunity, with extraordinarily bright prompt optical emission that peaked at a visual magnitude of 5.3, making it briefly visible with the naked eye. It was captured in exquisite detail by wide-field telescopes, imaging the burst location from before the time of the explosion. The combination of these unique optical data with simultaneous gamma-ray observations provides powerful diagnostics of the detailed physics of this explosion within seconds of its formation. Here we show that the prompt optical and gamma-ray emissions from this event likely arise from different spectral components within the same physical region located at a large distance from the source, implying an extremely relativistic outflow. The chromatic behaviour of the broadband afterglow is consistent with viewing the GRB down the very narrow inner core of a two-component jet that is expanding into a wind-like environment consistent with the massive star origin of long GRBs. These circumstances can explain the extreme properties of this GRB.

Bohdan Paczy'nski was one of the pioneers of the cosmological GRB model. His ideas on how GRBs operate and what are their progenitors have dominated the field of GRBs in the hectic nineties during which the distances and the origin of GRBs were revealed. I discuss here Bohdan's contributions in some historical perspective.

We point out that previous studies of possible Lorentz-violating effects in astronomical time-of-flight data did not take into account the entire implications of the universe's cosmological expansion. We present the derivation of the accurate formulation of the problem and show that the resulting correction of the limits on Lorentz violation is significant.

We calculate the very high energy (sub-GeV to TeV) inverse Compton emission of GRB afterglows. We argue that this emission provides a powerful test of the currently accepted afterglow model. We focus on two processes: synchrotron self-Compton (SSC) emission within the afterglow blast wave, and external inverse Compton (EIC) emission which occurs when flare photons (produced by an internal process) pass through the blast wave. We show that if our current interpretations of the Swift XRT data are correct, there should be a canonical high energy afterglow emission light curve. Our predictions can be tested with high energy observatories such as GLAST, Whipple, H.E.S.S. and MAGIC. Under favorable conditions we expect afterglow detections in all these detectors.

Before 2004, nearly all GRB afterglow data could be understood in the context of the external shocks model. This situation has changed in the past two years, when it became clear that some afterglow components should be attributed to the activity of the central engine; i.e., the {\it central engine afterglow}. We review here the afterglow emission that is directly related to the GRB central engine. Such an interpretation proposed by Katz, Piran & Sari, peculiar in pre-{\it Swift} era, has become generally accepted now.

We present self consistent calculations of Synchrotron self Compton (SSC) radiation that takes place within the afterglow blast wave and External inverse Compton (EIC) radiation that takes place when flare photons (produced by an internal process) pass through the blast wave. We show that if our current interpretations of the Swift XRT data are correct, there should be a canonical high energy afterglow emission light curve. We expect that GRBs with a long term X-ray flattening or X-ray flares should show similar high energy features. The EIC emission, however, is long lasting and weak and might be outshined by the SSC emission of the forward shock. The high energy emission could be well detected by the soon to be launched GLAST satellite. Its detection could shed new light on the conditions within the emitting regions of GRBs.

We compare the luminosity function and rate inferred from the BATSE long bursts peak flux distribution with those inferred from the Swift peak flux distribution. We find that both the BATSE and the Swift peak fluxes can be fitted by the same luminosity function and the two samples are compatible with a population that follows the star formation rate. The estimated local long GRB rate (without beaming corrections) varies by a factor of five from 0.05 Gpc^(-3)yr^(-1) for a rate function that has a large fraction of high redshift bursts to 0.27 Gpc^(-3)yr^(-1) for a rate function that has many local ones. We then turn to compare the BeppoSax/HETE2 and the Swift observed redshift distributions and compare them with the predictions of the luminosity function found. We find that the discrepancy between the BeppoSax/HETE2 and Swift observed redshift distributions is only partially explained by the different thresholds of the detectors and it may indicate strong selection effects. After trying different forms of the star formation rate (SFR) we find that the observed Swift redshift distribution, with more observed high redshift bursts than expected, is inconsistent with a GRB rate that simply follows current models for the SFR. We show that this can be explained by GRB evolution beyond the SFR (more high redshift bursts). Alternatively this can also arise if the luminosity function evolves and earlier bursts were more luminous or if strong selection effects affect the redshift determination.

On 2006 May 5, a four second duration, low energy, ~10^49 erg, Gamma-Ray Burst (GRB) was observed, spatially associated with a z=0.0894 galaxy. Here, we report the discovery of the GRB optical afterglow and observations of its environment using Gemini-south, Hubble Space Telescope (HST), Chandra, Swift and the Very Large Array. The optical afterglow of this GRB is spatially associated with a prominent star forming region in the Sc-type galaxy 2dFGRS S173Z112. Its proximity to a star forming region suggests that the progenitor delay time, from birth to explosion, is smaller than about 10 Myr. Our HST deep imaging rules out the presence of a supernova brighter than an absolute magnitude of about -11 (or -12.6 in case of ``maximal'' extinction) at about two weeks after the burst, and limits the ejected mass of radioactive Nickel 56 to be less than about 2x10^-4 solar mass (assuming no extinction). Although it was suggested that GRB 060505 may belong to a new class of long-duration GRBs with no supernova, we argue that the simplest interpretation is that the physical mechanism responsible for this burst is the same as for short-duration GRBs.

Over the last decade, long-duration gamma-ray bursts (GRBs) including the subclass of X-ray flashes (XRFs) have been revealed to be a rare variety of Type Ibc supernova (SN). While all these events result from the death of massive stars, the electromagnetic luminosities of GRBs and XRFs exceed those of ordinary Type Ibc SNe by many orders of magnitude. The essential physical process that causes a dying star to produce a GRB or XRF, and not just an SN, remains the crucial open question. Here we present radio and X-ray observations of XRF 060218 (associated with SN 2006aj), the second nearest GRB identified to-date, which allow us to measure its total energy and place it in the larger context of cosmic explosions. We show that this event is 100 times less energetic but ten times more common than cosmological GRBs. Moreover, it is distinguished from ordinary Type Ibc SNe by the presence of 10^48 erg coupled to mildly-relativistic ejecta, along with a central engine (an accretion-fed, rapidly rotating compact source) which produces X-rays for weeks after the explosion. This suggests that the production of relativistic ejecta is the key physical distinction between GRBs/XRFs and ordinary SNe, while the nature of the central engine (black hole or magnetar) may distinguish typical bursts from low-luminosity, spherical events like XRF 060218.

Over the past decade our physical understanding of gamma-ray bursts (GRBs) has progressed rapidly thanks to the discovery and observation of their long-lived afterglow emission. Long-duration (T < 2 s) GRBs are associated with the explosive deaths of massive stars (``collapsars''), which produce accompanying supernovae, while the short-duration (T > 2 s) GRBs arise from a different origin, which has been argued to be the merger of two compact objects, either neutron stars or black holes. Here we present observations of GRB060614, a 100-s long burst discovered by the Swift satellite, which require the invocation of a new explosive process: either a massive ``collapsar'' that powers a GRB without any associated supernova, or a new type of engine, as long-lived as the collapsar but without any such massive stellar host. We also discuss the properties of this burst's redshift z=0.125 host galaxy, which distinguish it from other long-duration GRBs and suggest that an entirely new type of GRB progenitor may be required.

We present observations of the unusually bright and long gamma-ray burst GRB 050820A, one of the best-sampled broadband data sets in the Swift era. The gamma-ray light curve is marked by a soft precursor pulse some 200 s before the main event; the lack of any intervening emission suggests that it is due to a physical mechanism distinct from the GRB itself. The large time lag between the precursor and the main emission enabled simultaneous observations in the gamma-ray, X-ray, and optical band-passes, something only achieved for a handful of events to date. While the contemporaneous X-rays are the low-energy tail of the prompt emission, the optical does not directly track the gamma-ray flux. Instead, the early-time optical data appear mostly consistent with the forward shock synchrotron peak passing through the optical, and are therefore likely the beginning of the afterglow. On hour time scales after the burst, the X-ray and optical light curves are inconsistent with an adiabatic expansion of the shock into the surrounding region, but rather indicate that there is a period of energy injection. Observations at late times allow us to constrain the collimation angle of the relativistic outflow to theta = 6.8 - 9.3 degrees. Our estimates of both the kinetic energy of the afterglow and the prompt gamma-ray energy release make GRB 050820A one of the most energetic events for which such values could be determined.

The nearby GRB 060216/SN 2006aj was an extremely long, weak and very soft GRB. It was peculiar in many aspects. We show here that the X-ray, ultraviolet/optical and radio afterglow of GRB 060218 have to be attributed to different physical processes arising from different emission regions. From the several components in this burst's afterglow only the radio afterglow can be interpreted in terms of the common external shock model. We infer from the radio that the blast wave's kinetic energy was $\sim 10^{50}$ erg and the circumburst matter had a constant rather than a wind profile. The lack of a "jet break" up to 22 days implies that the outflow was wide $\theta_j >1$. Even though the late X-ray afterglow decays normally it cannot result from an external shock because of its very steep spectrum. Furthermore, the implied kinetic energy would have produced far too much radio. We suggest that this X-ray afterglow could be attributed to a continued activity of the central engine that within the collapsar scenario could arise from fall-back accretion. "Central engine afterglow" may be common in under-luminous GRBs where the kinetic energy of the blast wave is small and the external shock does not dominate over this component. Such under-luminous GRBs might be very common but they are rarely recorded because they can be detected only from short distances.

Lorentz Invariance Violation (LIV) arises in various quantum-gravity theories. As the typical energy for quantum gravity is the Planck mass, $M_{pl}$, LIV will, most likely, be manifested at very high energies that are not accessible on Earth in the foreseeable future. One has to turn to astronomical observations. Time of flight measurement from different astronomical sources set current limits on the energy scale of possible LIV to $> 0.01 M_{pl}$ (for n=1 models) and $> 10^{-9} M_{pl}$ (for n=2). According to current models Gamma-Ray Bursts (GRBs) are accompanied by bursts of high energy ($\gsim 100$TeV) neutrinos. At this energy range the background level of currently constructed neutrino detectors is so low that a detection of a single neutrino from the direction of a GRB months or even years after the burst would imply an association of the neutrino with the burst and will establish a measurement of a time of flight delay. Such time of flight measurements provide the best way to observe (or set limits) on LIV. Detection of a single GRB neutrino would open a new window on LIV and would improve current limits by many orders of magnitude.

Gamma ray bursts are excellent candidates to constrain physical models which break Lorentz symmetry. We consider deformed dispersion relations which break the boost invariance and lead to an energy-dependent speed of light. In these models, simultaneously emitted photons from cosmological sources reach Earth with a spectral time delay that depends on the symmetry breaking scale. We estimate the possible bounds which can be obtained by comparing the spectral time delays with the time resolution of available telescopes. We discuss the best strategy to reach the strongest bounds. We compute the probability of detecting bursts that improve the current bounds. The results are encouraging. Depending on the model, it is possible to build a detector that within several years will improve the present limits of 0.015 m_pl.

Bright optical and X-ray flares have been observed in many Gamma-ray Burst (GRB) afterglows. These flares have been attributed to late activity of the central engine. In most cases the peak energy is not known and it is possible and even likely that there is a significant far-ultraviolet component. These far-ultraviolet photons escape our detection because they are absorbed by the neutral hydrogen before reaching Earth. However, these photons cross the blast wave produced by the ejecta that have powered the initial GRB. They can be inverse Compton upscattered by hot electrons within this blast wave. This process will produce a strong sub-GeV flare that follows the high energy (soft X-ray) tail of the far-UV flare but lasts much longer and can be detected by the upcoming {\em Gamma-Ray Large Area Telescope} (GLAST) satellite. This signature can be used to probe the spectrum of the underlying far-ultraviolet flare. The extra cooling produced by this inverse Compton process can lower the X-ray emissivity of the forward shock and explain the unexpected low early X-ray flux seen in many GRBs.

The discovery by Swift that a good fraction of Gamma Ray Bursts (GRBs) have a slowly decaying X-ray afterglow phase led to the suggestion that energy injection into the blast wave takes place several hundred seconds after the burst. This implies that right after the burst the kinetic energy of the blast wave was very low and in turn the efficiency of production of $\gamma$-rays during the burst was extremely high, rendering the internal shocks model unlikely. We re-examine the estimates of kinetic energy in GRB afterglows and show that the efficiency of converting the kinetic energy into gamma-rays is moderate and does not challenge the standard internal shock model. We also examine several models, including in particular energy injection, suggested to interpret this slow decay phase. We show that with proper parameters, all these models give rise to a slow decline lasting several hours. However, even those models that fit all X-ray observations, and in particular the energy injection model, cannot account self-consistently for both the X-ray and the optical afterglows of well monitored GRBs such as GRB 050319 and GRB 050401. We speculate about a possible alternative resolution of this puzzle.

On December 27, 2004, a giant flare from SGR 1806$-$20 was detected on earth. Its thermal spectrum and temperature suggest that the flare resulted from an energy release of about $10^{47}$ erg/sec close to the surface of a neutron star in the form of radiation and/or pairs. This plasma expanded under its own pressure producing a fireball and the observed gamma-rays escaped once the fireball became optically thin. The giant flare was followed by a bright radio afterglow, with an observable extended size, implying an energetic relativistic outflow. We revisit here the evolution of relativistic fireballs and we calculate the Lorentz factor and energy remaining in relativistic outflow once the radiation escapes. We show that pairs that arise naturally in a pure pairs-radiation fireball do not carry enough energy to account for the observed afterglow. We consider various alternatives and we show that if the relativistic outflow that causes the afterglow is related directly to the prompt flare, then the initial fireball must be loaded by baryons or Poynting flux. While we focus on parameters applicable to the giant flare and the radio afterglow of SGR 1806$-$20 the calculations presented here might be also applicable to GRBs.

The peak energy - isotropic energy (EpEi) relation is among the most intriguing recent discoveries concerning GRBs. It can have numerous implications on our understanding of the emission mechanism of the bursts and on the application of GRBs for cosmological studies. However, this relation was verified only for a small sample of bursts with measured redshifts. We propose here a test whether a burst with an unknown redshift can potentially satisfy the EpEi relation. Applying this test to a large sample of BATSE bursts we find that a significant fraction of those bursts cannot satisfy this relation. Our test is sensitive only to dim and hard bursts and therefore this relation might still hold as an inequality (i.e. there are no intrinsically bright and soft bursts). We conclude that the observed relation seen in the sample of bursts with a known redshift might be influenced by observational biases and from the inability to locate and well localize hard and weak bursts that have only a small number of photons. In particular we point out that the threshold for detection, localization and redshift measurement is essentially higher than the threshold for detection alone. We predict that Swift will detect some hard and weak bursts that would be outliers to the EpEi relation. However, we cannot quantify this prediction. We stress the importance of understanding the detection-localization-redshift threshold for the coming Swift detections.

Shortly after the discovery of the binary system PSR J0737-3039 we predicted that it has a small proper motion and suggested that this implies that the progenitor mass of the younger pulsar (B) must be around 1.45M_sun. This is in contradiction to standard evolutionary scenarios that suggest a He-star progenitor with M >2.1-2.3M_sun and in sharp contrast to a velocity of ~100 km/s predicted by these models. It requires a new previously unseen type of gravitational collapse. Our (generally ignored) prediction of a low peculiar motion was confirmed later by pulsar timing observations that found v_\bot < 30 km/s. Recently Willems et al. put forward an elaborate simulation suggesting that in spite of these observations, the most likely natal kick is large (50-170 km/s) and the corresponding most likely progenitor mass is large as well. We refute this claim and show that Willems et al. implicitly assume the result that they obtain.

We compare the luminosity function and rate inferred from the BATSE short hard bursts (SHBs) peak flux distribution with the redshift and luminosity distributions of SHBs observed by Swift/HETE II. The Swift/HETE II SHB sample is incompatible with SHB population that follows the star formation rate. However, it is compatible with a distribution of delay times after the SFR. This would be the case if SHBs are associated with binary neutron star mergers. The implied SHB rates that we find range from \sim 8 to \sim 30h_{70}^3Gpc^{-3}yr^{-1}. This rate is a much higher than what was previously estimated and, when beaming is taken into account, it is comparable to the rate of neutron star mergers estimated from statistics of binary pulsars. If GRBs are produced in mergers the implied rate practically guarantees detection by LIGO II and possibly even by LIGO I, if we are lucky. Our analysis, which is based on observed short hard burst is limited to bursts with luminosities above 10^{49}erg/sec. Weaker bursts may exist but if so they are hardly detected by BATSE or Swift and hence their rate is very weakly constrained by current observations. Thus the rate of mergers that lead to a detection of a gravitational radiation signal might be even higher.

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

We compare the luminosity function and rate inferred from the BATSE short hard bursts (SHBs) peak flux distribution with the redshift and luminosity distributions of SHBs observed by Swift/HETE II. While the Swift/HETE II SHB sample is incompatible with SHB population that follows the star formation rate, it is compatible with a SHB rate that reflect a distribution of delay times after the SFR. This would be the case if SHBs are associated with binary neutron star mergers. The available data allows, however, different interpretations. For example, a population whose rate is independent of the redshift fits the data very well. The implied SHB rates that we find range from $\sim 8$ to $\sim 30h_{70}^3$Gpc$^{-3}$yr$^{-1}$. This is a much higher rate than what was previously estimated. A detailed (2 dimensional) look at the best fit models shows, however, some discrepancy between the four Swift/HETE II SHBs and the models based on BATSE SHBs. This could be a statistical fluke. It could also arise from wrong estimates of the triggering criteria or from selection effects. If real it may indicate the existence of two SHB populations with different luminosity functions and redshift distributions.

According to current models, gamma-ray bursts (GRBs) are produced when the energy carried by a relativistic outflow is dissipated and converted into radiation. The efficiency of this process, $\epsilon_\gamma$, is one of the critical factors in any GRB model. The X-ray afterglow light curves of Swift GRBs show an early stage of flat decay. This has been interpreted as reflecting energy injection. When combined with previous estimates, which have concluded that the kinetic energy of the late ($\gtrsim 10$ hr) afterglow is comparable to the energy emitted in gamma-rays, this interpretation implies very high values of $\epsilon_\gamma \gtrsim 0.9$. Such a high efficiency is hard to reconcile with most models, including in particular the popular internal-shocks model. We re-analyze the derivation of the kinetic energy from the afterglow X-ray flux and re-examine the resulting estimates of the efficiency. We confirm that, if the flatt decay arises from energy injection and the pre-Swift broad-band estimates of the kinetic energy are correct, then $\epsilon_\gamma \gtrsim 0.9$. We point out, however, that another interpretation of the flattish decay - a variable X-ray afterglow efficiency (e.g., due to a time dependence of afterglow shock microphysical parameters) - is possible. We also show that direct estimates of the kinetic energy from the late X-ray afterglow flux are sensitive to the assumed values of the shock microphysical parameters and suggest that broad-band afterglow fits might have underestimated the kinetic energy (e.g., by overestimating the fraction of electrons that are accelerated to relativistic energies). Either one of these possibilities implies a lower gamma-ray efficiency, and their joint effect could conceivably reduce the estimate of the typical $\epsilon_\gamma$ to $\sim 0.1-0.5$.

Despite a rich phenomenology, gamma-ray bursts (GRBs) are divided into two classes based on their duration and spectral hardness -- the long-soft and the short-hard bursts. The discovery of afterglow emission from long GRBs was a watershed event, pinpointing their origin to star forming galaxies, and hence the death of massive stars, and indicating an energy release of about 10^51 erg. While theoretical arguments suggest that short GRBs are produced in the coalescence of binary compact objects (neutron stars or black holes), the progenitors, energetics, and environments of these events remain elusive despite recent localizations. Here we report the discovery of the first radio afterglow from a short burst, GRB 050724, which unambiguously associates it with an elliptical galaxy at a redshift, z=0.257. We show that the burst is powered by the same relativistic fireball mechanism as long GRBs, with the ejecta possibly collimated in jets, but that the total energy release is 10-1000 times smaller. More importantly, the nature of the host galaxy demonstrates that short GRBs arise from an old (>1 Gyr) stellar population, strengthening earlier suggestions, and providing support for coalescing compact object binaries as the progenitors.

We present an analytic model for the sizes of voids in the galaxy distribution. Peebles and others have recently emphasized the possibility that the observed characteristics of voids may point to a problem in galaxy formation models, but testing these claims has been difficult without any clear predictions for their properties. In order to address such questions, we build a model to describe the distribution of galaxy underdensities. Our model is based on the "excursion set formalism," the same technique used to predict the dark matter halo mass function. We find that, because of bias, galaxy voids are typically significantly larger than dark matter voids and should fill most of the universe. We show that voids selected from catalogs of luminous galaxies should be larger than those selected from faint galaxies: the characteristic radii range from ~5-10h^{-1} Mpc for galaxies with absolute r-band magnitudes M_r-5 log(h) < -16 to -20. These are reasonably close to, though somewhat smaller than, the observed sizes. The discrepancy may result from the void selection algorithm or from their internal structure. We also compute the halo populations inside voids. We expect small haloes (M<10^{11} Msun) to be up to a factor of two less underdense than the haloes of normal galaxies. Within large voids, the mass function is nearly independent of the size of the underdensity, but finite-size effects play a significant role in small voids (R<7h^{-1} Mpc).

We calculate the reverse shock (RS) synchrotron emission in the optical and the radio wavelength bands from electron-positron pair enriched gamma-ray burst ejecta with the goal of determining the pair content of GRBs using early time observations. We take into account an extensive number of physical effects that influence radiation from the reverse-shock heated GRB ejecta. We find that optical/IR flux depends very weakly on the number of pairs in the ejecta, and there is no unique signature of ejecta pair enrichment if observations are confined to a single wavelength band. It may be possible to determine if the number of pairs per proton in the ejecta is > 100 by using observations in optical and radio bands; the ratio of flux in the optical and radio at the peak of each respective reverse-shock light curve is dependent on the number of pairs per proton. We also find that over a large parameter space, RS emission is expected to be very weak; GRB 990123 seems to have been an exceptional burst in that only a very small fraction of the parameter space produces optical flashes this bright. Also, it is often the case that the optical flux from the forward shock is brighter than the reverse shock flux at deceleration. This could be another possible reason for the paucity of prompt optical flashes with a rapidly declining light curve at early times as was seen in 990123 and 021211. Some of these results are a generalization of similar results reported in Nakar & Piran (2004).

It is generally accepted that neutron stars form in core collapse events that are accompanied by a supernovae (types II or Ib or Ic). Typical progenitors are, therefore, larger than $\sim 2.1 \Ms$. We suggest \cite{PS04,PS05} that the binary pulsar J0737-3039 provides evidence for a new formation channel: collapse of a light progenitor. This binary pulsar J0737-3039 has several remarkable features including among others: a very tight orbit with a Keplerian velocity of ~600km/sec, a low eccentricity, and a location ~50pc from the Galactic plane implying that the system has, at high likelihood, a small (compared to Keplerian) center of mass velocity. A significant mass loss during the formation of the second pulsar would have lead either to an eccentric orbit or to a large center of mass velocity or to both. Therefore, we can set a strong upper limit on the progenitor's mass. A progenitor more massive than 1.9$\Ms$ is ruled out (at 97% confidence). The kinematically favored option is of a progenitor mass around 1.45$\Ms$. Recent evidence for a rather low velocity proper motion supports this prediction and decreases the likelihood of the standard (high mass progenitor) scenario. Lack of variations in the pulses' profiles (which indicate no significant geodetic precession) provides further support for the no kick and low progenitor mass formation scenario.

The final chapter in the long-standing mystery of the gamma-ray bursts (GRBs) centres on the origin of the short-hard class, suspected on theoretical grounds to result from the coalescence of neutron star or black hole binary systems. Numerous searches for the afterglows of short-hard bursts have been made, galvanized by the revolution in our understanding of long-duration GRBs that followed the discovery in 1997 of their broadband (X-ray, optical, and radio) afterglow emission. Here we present the discovery of the X-ray afterglow of a short-hard burst whose accurate position allows us to unambiguously associate it with a star-forming galaxy at redshift z=0.160, and whose optical lightcurve definitively excludes a supernova association. Together with results from three other recent short-hard bursts, this suggests that short-hard bursts release much less energy than the long-duration GRBs. Models requiring young stellar populations, such as magnetars and collapsars, are ruled out, while coalescing degenerate binaries remain the most promising progenitor candidates.

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\sigma$ 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.

We suggest that part of the infalling material during the core-collapse of a massive star goes into orbit around the compact core to form a hot, dense, centrifugally-supported accretion disk whose evolution is strongly influenced by neutrino interactions. Under a wide range of conditions, this neutrino-dominated accretion flow will be advection-dominated and will develop a substantial outflowing wind. We estimate the energy carried out in the wind and find that it exceeds $10^{50}$ erg for a wide range of parameters and even exceeds 10^{51} erg for reasonable parameter choices. We propose that the wind energy will revive a stalled shock and will help produce a successful supernova explosion. We discuss the role of the disk wind in both prompt and delayed explosions. While both scenarios are feasible, we suggest that a delayed explosion is more likely, and perhaps even unavoidable. Finally, we suggest that the disk wind may be a natural site for r-process nucleosynthesis.

In Nakar & Piran (2004) we have shown that a large fraction of BATSE bursts cannot satisfy the Amati relation between the isotropic energy - peak energy of GRBs. Our results have been corroborated by a further analysis carried out by Band and Preece (2005). Recently Ghirlanda et al. (2005) claimed that when the exact distribution of the scatter around the Amati relation is taken into account, the BATSE sample is consistent with the Amati relation. We show that the probability that the BATSE sample is consistent with the Amati relation (including the scattering), as found by Ghirlanda et al. (2005) for a sample of 23 bursts with known redshift, is vanishingly small.

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

Magnetic fields play a crucial role in the physics of Gamma-Ray Bursts (GRBs). Strong observational evidence indicates that the observed afterglow and most likely the prompt emission arise from synchrotron emission. It is possible that Poynting flux plays an important or even dominant role in the relativistic outflow from the inner engine, but like in other astronomical relativistic jets this suggestion is controversial. Finally, it is likely that magnetic fields larger than $10^{15}$ G occur within GRBs' inner engines and contribute to the acceleration and collimation of the relativistic jets. I review here the GRB fireball model and discuss the role that magnetic fields play in its various components. I suggest that the early afterglow, that reflects the initial interaction of the relativistic jet with its surrounding matter is the best available tool to explore the nature of relativistic outflow in astronomical relativistic jets.

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 $\gamma$-ray emission results from internal shocks, the early reverse shock emission should resemble a smoothed version of the prompt $\gamma$-ray light curve.

Using the BATSE peak flux distribution we rederive the short GRBs luminosity function and compare it with the observed redshift distribution of long bursts. We show that both distributions are compatible with the assumption that short as well as long bursts follow the star formation rate. In this case the difference between the two observed distributions can be interpreted as arising mostly from differences in the detector's sensitivity to long and short bursts, while the local rate of short bursts is 0.1 h_70^3 Gpc^-3 yr^-1 . We also consider the possibility that short GRBs may be associated with binary neutron star mergers and estimate the effect of the merging time delay on the luminosity function and redshift distribution. We find that in this case the local rate of short GRBs is ~0.8 h_70^3 Gpc^-3 yr^-1. Assuming that all binary merging systems lead to short GRBs, we find a typical jet opening angle of 1.6^o

Evolutionary scenarios suggest that the progenitor of the new binary pulsar J0737-3039B \cite{ref1,ref2} was a He-star with $M > 2.1-2.3~\Ms$ \cite{ref3,ref4}. We show that this case implies that the binary must have a large ($>120$ km/s) center of mass velocity. However, the location, $\sim 50$ pc from the Galactic plane, suggests that the system has, at high likelihood, a significantly smaller center of mass velocity and a progenitor more massive than 2.1~$\Ms$ is ruled out (at 97% c.l.). A progenitor mass around 1.45~$\Ms$, involving a new previously unseen gravitational collapse, is kinematically favored. The low mass progenitor is consistent with the recent scintillations based velocity measurement of 66$\pm 15$ km/s \cite{ref12new} (and which also rules out the high mass solution at 99% c.l.) and inconsistent with the higher earlier estimates of 141$\pm 8.5$ km/s \cite{ref11new}. Direct proper motion measurements, that should be available within a year or so, should better help to distinguish between the two scenarios.

Gamma ray bursts (GRBs) are short and intense pulses of $\gamma$-rays arriving from random directions in the sky. Several years ago Amelino-Camelia et al. pointed out that a comparison of time of arrival of photons at different energies from a GRB could be used to measure (or obtain a limit on) possible deviations from a constant speed of light at high photons energies. I review here our current understanding of GRBs and reconsider the possibility of performing these observations.

The realization that the total energy of GRBs is correlated with their jet break angles motivates the search for a similar relation between the peak luminosity, $L$, and the jet break angles, $L\propto\theta^{-2}$. Such a relation implies that the GRB luminosity function determines the angular distribution. We re-derive the GRB luminosity function using the BATSE peak flux distribution and compare the predicted distribution with the observed redshift distribution. The luminosity function can be approximated by a broken power law with a break peak luminosity of $4.4 \times 10^{51}$ erg/sec, a typical jet angle of 0.12 rad and a local GRB rate of $0.44 h_{65}^{3}$ Gpc$^{-3}$yr$^{-1}$. The angular distribution implied by $L\propto\theta^{-2}$ agrees well with the observed one, and implies a correction factor to the local rate due to beaming of $75\pm25$ (instead of 500 as commonly used). The inferred overall local GRB rate is $33\pm11 h_{65}^{3}$ Gpc$^{-3}$yr$^{-1}$. The luminosity function and angle distribution obtained within the universal structured jet model, where the angular distribution is essentially $\propto \theta$ and hence the luminosity function must be $\propto L^{-2}$, deviate from the observations at low peak fluxes and, correspondingly, at large angles. The corresponding correction factor for the universal structure jet is $\sim 20 \pm 10$.

The Gamma-Ray burst (GRB) - afterglow transition is one of the most interesting and least studied GRB phases. During this phase the relativistic ejecta begins interacting with the surrounding matter. A strong short lived reverse shock propagates into the ejecta (provided that it is baryonic) while the forward shock begins to shape the surrounding matter into a Blandford-McKee profile. We suggest a parametrization of the early afterglow light curve and we calculate (analytically and numerically) the observed parameters that results from a reverse shock emission (in an interstellar medium [ISM] environment). We present a new fingerprint of the reverse shock emission that is added to the well known $t^{-2}$ optical decay. Observation of this signature would indicate that the reverse shock dominates the emission during the early afterglow. The existence of a reverse shock will in turn imply that the relativistic ejecta contains a significant baryonic component. This signature would also imply that the surrounding medium is an ISM. We further show that: (i) The reverse shock optical flash depends strongly on initial conditions of the relativistic ejecta. (ii) Previous calculations have generally overestimated the strength of this optical flash. (iii) If the reverse shock dominates the optical flash then detailed observations of the early afterglow light curve would possibly enable us to determine the initial physical conditions within the relativistic ejecta and specifically to estimate its Lorentz factor and its width.

Taylor et al. (2004) reported recently on the first direct measurement of the apparent size of a GRB afterglow. Here we confront these observations with the predictions of the fireball model. We use a simple model to study numerically the evolution of the fireball and determine its apparent size, starting at the early spherical expansion phase through the jet break and into the Newtonian regime. We perform these calculations on the background of a uniform interstellar medium and a massive stellar wind environment. We find that the calculated apparent size at the time of the measurements taken by Taylor et al. depends only weakly on the chosen parameters (the jet opening angle, the energy, and the external density profile). Thus it is reassuring that the only possible outcome of the model, within the plausible range of parameters, agrees very well with the data. These measurements therefore present us with a strong test of possible GRB models, which the fireball model passes successfully.

Gamma-Ray Bursts (GRBs), short and intense pulses of low energy gamma-rays, have fascinated astronomers and astrophysicists since their unexpected discovery in the late sixties. During the last decade, several space missions: BATSE (Burst and Transient Source Experiment) on Compton Gamma-Ray Observatory, BeppoSAX and now HETE II (High-Energy Transient Explorer), together with ground optical, infrared and radio observatories have revolutionized our understanding of GRBs showing that they are cosmological, that they are accompanied by long lasting afterglows and that they are associated with core collapse Supernovae. At the same time a theoretical understanding has emerged in the form of the fireball internal-external shocks model. According to this model GRBs are produced when the kinetic energy of an ultra-relativistic flow is dissipated in internal collisions. The afterglow arises when the flow is slowed down by shocks with the surrounding circum-burst matter. This model has numerous successful predictions like the prediction of the afterglow itself, the prediction of jet breaks in the afterglow light curve and of an optical flash that accompanies the GRBs themselves. In this review I focus on theoretical aspects and on physical processes believed to take place in GRBs.

Recently we have presented a new theoretical analysis of the reverse shock emission. We use this analysis here to revisit GRB 990123. We find new and compelling evidences that the optical flash and the radio flare of GRB 990123 resulted from a reverse shock. This suggests that a significant fraction of the energy of the relativistic ejecta must have been carried by baryons. It also suggests that the external medium is an ISM and that in this burst the reverse shock emission dominated at early time over other possible processes. We use the early optical emission to constrain the physical parameters of the original ejecta and the microscopic parameters in the emitting reverse shocked region.

The observed orbital parameters of the recently discovered binary pulsar J0737-3039 are used to constrain the progenitor system. In particular, the small observed eccentricity and the small inferred peculiar velocity imply that during the formation of the second NS in the system, a very small mass Delta M <~ 0.15 Msun, was ejected. A progenitor more massive than 1.4 Msun is unlikely and a progenitor more massive than 2.3 Msun is practically ruled out. We therefore argue that the companion removed most of the progenitor envelope prior to its collapse. The second collapse may have been a SN of Type Ib, or more probably, of Type Ic. Alternatively, it could have been a bare collapse, an option which is kinematically favored. Future constraints on the peculiar motion could solidify this constraint on the mass loss. A future determination of the pulsars' spin vectors and proper motion could allow us the complete determination of the conditions during the second collapse.

We consider some recent developments in GRB/afterglow observations: (i) the appearance of a very hard prompt component in GRB 941017, and (ii) variability in the early afterglow light curves of GRB 021004 and GRB 030329. We show that these observations fit nicely within the internal-external shocks model. The observed variability indicates that the activity of the inner engine is more complicated than was thought earlier and that it involves patchy shells and refreshed shocks. We refute the claims of Berger et al. and of Sheth et al. that the radio and mm observations of GRB 030329 are inconsistent with refreshed shocks.

We compare two possible scenarios for the producing of high level of polarization within the prompt emission of a GRB: synchrotron emission from a relativistic jet with a uniform (in space and time) magnetic field and synchrotron emission from a jet with a random magnetic field in the plane of the shock. Somewhat surprisingly we find that both scenarios can produce a comparable level of polarization ($\sim 45-50$% for the uniform field and $\sim 30-35$% for a random field). Uniform time independent field most naturally arises by expansion of the field from the compact object. It requires a $10^{12}$G field at the source and a transport of the field as $\propto R^{-1}$. It {\it does not} imply Poynting flux domination of the energy of the wind. There is a serious difficulty however, within this scenario, accounting for particle acceleration (which requires random magnetic fields) both for Poynting flux and non-Poynting flux domination. Significant polarization can also arise from a random field provided that the observer is located within $1/\Gamma$ orientation from a narrow ($\theta_j \sim 1/\Gamma$) jet. While most jets are wider, the jet of GRB 021206 from which strong polarization was recently observed, was most likely very narrow. GRB 021206 is among the strongest bursts ever. Adopting the energy-angle relation we find an estimated angle of $<1/40$rad or even smaller. Thus, for this particular burst the required geometry is not unusual. We conclude that the RHESSI observations suggest that the prompt emission results from synchrotron radiation. However, in view of the comparable levels of polarizations predicted by both the random field and the homogeneous field scenarios these observations are insufficient to rule out or confirm either one.

There is growing interest in quantum-spacetime models in which small departures from Lorentz symmetry are governed by the Planck scale. In particular, several studies have considered the possibility that these small violations of Lorentz symmetry may affect various astrophysical observations, such as the evaluation of the GZK limit for cosmic rays, the interaction of TeV photons with the Far Infrared Background and the arrival time of photons with different energies from cosmological sources. We show that the same Planck-scale departures from Lorentz symmetry that lead to a modification of the GZK limit which would be consistent with the observations reported by AGASA, also have significant implications for the evaluation of the Bahcall-Waxman bound on the flux of high-energy neutrinos produced by photo-meson interactions in sources of size not much larger than the proton photo-meson mean free path.

GRB 030329 is unique in many aspects. It has a very low redshift for a GRB, $z=0.1685$, and is therefore very bright and easy to monitor, making it the most well studied afterglow to date. It shows a supernova bump in the light curve, with a spectrum very similar to SN 1998bw, thus establishing with much better confidence the connection between GRBs and core collapse SNe. There are also two important physical characteristics that make this burst especially interesting, aside from its remarkably low redshift. First, unlike most GRB afterglows, the light curve of GRB 030329 shows a very large variability a few days after the burst. These fluctuations show a roughly constact amplitude, and a constant duration $\Delta t$, while $\Delta t/t$ decreases with time $t$. Second, its $\gamma$-ray energy output and X-ray luminosity at $10 $hr are a factor of $\sim 20$ and $\sim 30$, respectively, below the average value around which most GRBs are narrowly clustered. We consider several interpretations for the variability in the light curve, in the context of different physical mechanisms, and find that the most likely cause is refreshed shocks, i.e. slow shells that are ejected from the source and catch up with the afterglow shock at late times. In GRB 030329 this happens after the jet break, which implies an approximately constant duration $\Delta t$ of the bumps, in agreement with the observations. This interpretation also explains the anomalously low initial energy of this burst, as the total energy of the afterglow shock is increased by a factor of $\sim 10$ due to the refreshed shocks, thus bringing the total energy output close to the average value for all GRBs.

The fluctuations observed in the light curves of some GRB-afterglows (such as GRB 021004) provide a useful tool to probe the circum-burst density profile and to probe the variations in the energy of blast-wave with time. We present a general formalism that reduces the calculation of the observed light curve from a Blandford-Mackee blast-wave to the evaluation of a one dimensional integral. Using this formalism we obtain a simple approximation to the general light curve that arises in more complex situations where the afterglow's energy or the external density vary. The solution is valid for spherically symmetric profiles and it takes a full consideration of the angular time delay effects. We present the light curves of several external density profiles and demonstrate the effects of density variations on the light curve. We also re-visit the afterglow of GRB021004 and we find that the steep decay after the first bump ($\sim 4000sec$) cannot result from a spherically symmetric density variation or from the passage of the synchrotron frequency through the optical band. This suggests that an angular structure is responsible to some of the observed features in the light curve. This may be the first evidence that an angular structure is important in the early stages of the afterglow.

We present general analytic expressions for GRB afterglow light curves arising from a variable external density profile and/or a variable energy in the blast wave. The former could arise from a clumpy ISM or a variable stellar wind; The latter could arise from refreshed shocks or from an angular dependent jet structure (patchy shell). Both scenarios would lead to a variable light curve. Our formalism enables us to invert the observed light curve and obtain possible density or energy profiles. The optical afterglow of GRB 021004 was detected 537 seconds AB (after the burst) (Fox et al. 2002). Extensive follow up observations revealed a significant temporal variability. We apply our formalism to the R-band light curve of GRB 021004 and we find that several models provide a good fit to the data. We consider the patchy shell model with $p=2.2$ as the most likely explanation. According to this model our line of sight was towards a "cold spot" that has lead to a relativity low $\gamma$-ray flux and an initially weak afterglow (while the X-ray afterglow flux after a day was above average). Observations above the cooling frequency, $\nu_c$, could provide the best way to distinguish between our different models.

We investigate how the angular structure of GRB jets effects the afterglow light curves at different viewing angles, $\theta_{v}$, from the jet symmetry axis. A numerical hydrodynamical modeling for the evolution of a relativistic collimated outflow, as it interacts with the surrounding medium, is carried out, and compared to two simple models that make opposite and extreme assumptions for the degree of lateral energy transfer. The Lorentz factor, $\Gamma$, and kinetic energy per unit solid angle, $\epsilon$, are initially taken to be power laws of the angle $\theta$ from the jet axis. We find that the lateral velocity in the comoving frame, $v'_\theta$, is typically much smaller than the sound speed, $c_s$, as long as $\Gamma\gg 1$, and the dynamics of relativistic structured jets may be reasonably described by a simple analytic model where $\epsilon$ is independent of time, as long as $\Gamma(\theta=0)\ga$ a few. We perform a qualitative comparison between the resulting light curves and afterglow observations. This constrains the jet structure, and poses problems for a `universal' jet model, where all GRB jets are assumed to be intrinsically identical, and differ only by our viewing angle, $\theta_{v}$.

In many GRB inner engine models the highly relativistic GRB jets are engulfed by slower moving matter. This could result in different beaming for the prompt $\gamma$-ray emission and for the lower energy afterglow. In this case we should expect that some observer will see {\it on-axis orphan afterglows}: X-ray, optical and radio afterglows within the initial relativistic ejecta with no preceding GRB; The prompt $\gamma$-ray emission is pointing elsewhere. We show that the observations of the WFC on BeppoSAX constrain with high certainty the prompt X-ray beaming factor to be less than twice the prompt $\gamma$-ray beaming. The results of Ariel 5 are consistent with this interpretation. The RASS from ROTSE and HEAO-1 constrain the X-ray beaming factor at 400 and 20 minutes after the burst respectively to be comparable and certainly not much larger than the $\gamma$-ray beaming factor. There is no direct limit on the optical beaming. However, we show that observations of several months with existing hardware could result in a useful limit on the optical beaming factor of GRB afterglows.

The realization that GRBs release a constant amount of energy implies that the post jet-break afterglow evolution would be largely universal. For a given redshift all afterglows should be detected up to a fixed observer angle. We estimate the observed magnitude and the implied detectability of orphan afterglows. We show that for reasonable limiting magnitudes ($m_{lim}=25$) orphan afterglows would typically be detected from small ($\sim 10^\circ$) angles away from the GRB jet axis. A detected orphan afterglow would generally correspond to a "near-miss" of the GRB whose jet was pointing just slightly away from us. With our most optimistic parameters we expect that 15 orphan afterglows will be recorded in the SDSS and 35 transients will be recorded in a dedicated 2m class telescope operating full time for a year in an orphan afterglow search. The rate is smaller by a factor of 15 for our "canonical" parameters. We show that for a given facility an optimal survey should be shallower, covering larger area rather than deeper. The limiting magnitude should not be, however, lower than $\sim$ 23rd as in this case more transients from on-axis GRBs will be discovered than orphan afterglows. About 15% of the transients could be discovered with a second exposure of the same area provided that it follows after 3, 4 and 8 days for $m_{lim}=23$, 25 and 27.

The realization that GRBs are narrowly beamed implied that the actual rate of GRBs is much larger than the observed one. There are 500 unobserved GRBs for each observed one. The lack of a clear trigger makes it hard to detect these unobserved GRBs as orphan afterglows. At late time, hundreds or thousands of years after a GRB, we expect to observe a GRB remnant (GRBR). These remnants could be distinguished from the more frequent SNRs using their different morphology. While SNRs are spherical, GRBRs that arise from a highly collimated flow, are expected to be initially nonspehrical. Using SPH simulations we follow the evolution of a GRBR and calculate the image of the remnant produced by bremsstrahlung and by synchrotron emission. We find that the GRBR becomes spherical after $ \sim 3000 {\rm yr}(E_{51}/n)^{1/3}$ at $R\sim 12{\rm pc} (E_{51}/n)^{1/3}$, where $E_{51}$ is the initial energy in units of $10^{51}{\rm erg}$ and $n$ is the surrounding ISM number density in ${\rm cm}^{-3}$. We expect $0.5 (E_{51}/n)^{1/3}$ non-spherical GRBs per galaxy. Namely, we expect $\sim ~ 20$ non spherical GRBRs with angular sizes $\sim 1 \mu$arcsec within a distance of 10Mpc. These results are modified if there is an underlying spherical supernova. In this case the GRBR will remain spherical only for $\sim150 {\rm yr} (E_{51}/n)^{1/3}$ and the number of non-spherical GRBRs is smaller by a factor of 10 and their size is smaller by a factor of 3.

The nature of the `inner engine' that accelerate and collimate the relativistic flow at the cores of GRBs is the most interesting current puzzle concerning GRBs. Numerical simulations have shown that the internal shocks' light curve reflects the activity of this inner engine. Using a simple analytic toy model we clarify the relations between the observed $ \gamma $-rays light curve and the inner engine's activity and the dependence of the light curves on the inner engine's parameters. This simple model also explains the observed similarity between the observed distributions of pulses widths and the intervals between pulses and the correlation between the width of a pulse and the length of the preceding interval. Our analysis suggests that the variability in the wind's Lorentz factors arises due to a modulation of the mass injected into a constant energy flow.

According to the fireball model GRBs are produced when a relativistic flow is dissipated by shocks. The questions how does a GRB's "inner engine" accelerates and collimated the relativistic flow is today the most interesting (and most difficult) open question concerning GRBs. We summarize, here, the known constrains on the "inner engines". These constrains arise mostly from the temporal structure seen in the light curves of GRBs. We discuss the well-known NO GO theorem: An explosive source cannot produce a variable GRB. We show that various attempts to find a way around this theorem are rules out by new observational data: similarity between pulse width and pulse intervals and a correlation between a pulse and the preceding interval. We show that a simple toy model explains these features in terms of the internal shocks model. We discuss the implications of these findings to GRBs' inner engines.

We analyze a sample of bright short bursts from the BATSE 4B-catalog and find that many short bursts are highly variable ($\delta t_{min}/T\ll 1$, where $\delta t_{min}$ is the shortest pulse duration and T is the burst duration). This indicates that it is unlikely that short bursts are produced by external shocks. We also analyze the available (first 1-2 seconds) high resolution (TTE) data of some of the long bursts. We find that variability on a 10ms time-scale is common in long bursts. This result shows that some long bursts are even more variable than it was thought before ($\delta t_{min}/T\approx 10^{-4}-10^{-3}$ )

We analyze a sample of bright long bursts and find that the pulses duration have a lognormal distribution while the intervals between pulses have an excess of long intervals (relative to lognormal distribution). This excess can be explained by the existence of quiescent times, long periods with no signal above the background. The lognormal distribution of the intervals (excluding the quiescent times) is similar to the distribution of the pulses width. This result suggests that the quiescent times are made by a different mechanism than the rest of the intervals. It also suggests that the intervals (excluding the quiescent times) and the pulse width are connected to the same parameters of the source. We find that there is a correlation between a pulse width and the duration of the interval preceding it. There is a weaker, but still a significant, correlation between a pulse width and the interval following it. The significance of the correlation drops substantially when the intervals considered are not adjacent to the pulse.

The energy of a Gamma-Ray Burst is one of the most interesting factors that can help determining the origin of these mysterious explosions. After the discovery that GRBs are cosmological it was thought, for a while that they are standard candles releasing $\sim 10^{51}$ergs. Redshift measurements that followed the discovery of GRB afterglow revealed that GRBs are (i) further than what was initially believed and (ii) have a very wide luminosity function. Some bursts revealed an enormous energy release with a record estimate of $ 1.4 \times 10^{54}$ergs for GRB990123. The energy budget was stretched even further when it was suggested that the conversion efficiency of producing gamma-rays is low. The realization that GRBs are beamed changed this perspective and brought the energy budget of GRBs down back to a "modest" $\sim 5 \times 10^{50}$ergs. I discuss here various estimates showing that GRB energy is narrowly distributed and discuss the implications of this conclusion to GRB models.

The energy release in gamma-ray bursts is one of the most interesting clues on the nature of their "inner engines". We show here that the total energy release in GRBs varies by less than one order of magnitude from one burst to another while the energy emitted in $\gamma$-ray photons varies by more than an order of magnitude. This result indicates that the central engine of long duration GRB has a remarkably constant energy output which provides very important constraint on the nature of these enigmatic explosions. The broader distribution of the observed $\gamma$-ray flux, about three orders of magnitude in width, can be attributed, in part, to a variation in the opening angle of the collimated explosion, and in part to the variation of Lorentz factor across the jet as well as a variable efficiency for converting the kinetic energy of explosion to $\gamma$-rays.

The most interesting current open question in the theory of GRB afterglow is the propagation of jetted afterglows during the sideway expansion phase. Recent numerical simulations show hydrodynamic behavior that differs from the one suggested by simple analytic models. Still, somewhat surprisingly, the calculated light curves show a `jet break' at about the expected time. These results suggest that the expected rate of orphan optical afterglows should be smaller than previously estimated.

Many models of gamma-ray bursts (GRBs) involve accretion onto a compact object, usually a black hole, at a mass accretion rate of order a fraction of a solar mass per second. If the accretion disk is larger than a few tens or hundreds of Schwarzschild radii, the accretion will proceed via a convection-dominated accretion flow (CDAF) in which most of the matter escapes to infinity rather than falling onto the black hole. Models involving the mergers of black hole white dwarf binaries and black hole helium star binaries fall in this category. These models are unlikely to produce GRBs since very little mass reaches the black hole. If the accretion disk is smaller, then accretion will proceed via neutrino cooling in a neutrino-dominated accretion disk (NDAF) and most of the mass will reach the center. Models involving the mergers of double neutron star binaries and black hole neutron star binaries fall in this category and are capable of producing bright GRBs. If the viscosity parameter $\alpha$ in the NDAF has a standard value $\sim0.1$, these mergers can explain short GRBs with durations under a second, but they are unlikely to produce long GRBs with durations of tens or hundred of seconds. If the accretion disk is fed by fallback of material after a supernova explosion, as in the collapsar model, then the time scale of the burst is determined by fallback, not accretion. Such a model can produce long GRBs. Fallback models again require that the accretion should proceed via an NDAF rather than a CDAF in order for a significant amount of mass to reach the black hole. This condition imposes an upper limit on the radius of injection of the gas.

We introduce an adaptation of the well known Tree+SPH numerical scheme to Post Newtonian (PN) hydrodynamics and gravity. Our code solves the (0+1+2.5)PN equations. These equations include Newtonian hydrodynamics and gravity (0PN), the first order relativistic corrections to those (1PN) and the lowest order gravitational radiation terms (2.5PN). We test various aspects of our code using analytically solvable test problems. We then proceed to study the 1PN effects on binary neutron star coalescence by comparing calculations with and without the 1PN terms. We find that the effect of the 1PN terms is rather small. The largest effect arises with a stiff equation of state for which the maximum rest mass density increases. This could induce black hole formation. The gravitational wave luminosity is also affected.

If a light gluino exists, the lightest gluino-containing baryon, the \OSO, is a possible candidate for self-interacting dark matter. In this scenario, the simplest explanation for the observed ratio $\Omega_{dm}/\Omega_b \approx 6-10$ is that $m_{S^0} \sim 900$\MeVcs; this is not at present excluded by particle physics. Such an \OSO could be present in neutron stars, with hyperon formation serving as an intermediate stage. We calculate equilibrium compositions and equation of state for high density matter with the \OSO, and find that for a wide range of parameters the properties of neutron stars with the \OSO are consistent with observations. In particular, the maximum mass of a nonrotating star is $1.7-1.8 M_\odot$, and the presence of the \OSO is helpful in reconciling observed cooling rates with hyperon formation.

Gamma-Ray Bursts (GRBs) are the brightest objects observed. They are also the most relativistic objects known so far. GRBs occur when an ultrarelativisitic ejecta is slowed down by internal shocks within the flow. Relativistic particles accelerated within these shocks emit the observed gamma-rays by a combination of synchrotron and inverse Compton emission. External shocks with the circumstellar matter slow down further the ejecta and produce the afterglow, which lasts for months. Comparison of the predictions of this fireball model with observations confirm a relativistic macroscopic motion with a Lorentz factor of $\Gamma \ge 100$. Breaks in the light curves of the afterglow indicate that GRBs are beamed with typical opening angles of a few degrees. The temporal variability of the gamma-rays signal provide us with the best indirect evidence on the nature of the ``internal engine'' that powers the GRBs and accelerates the relativistic ejecta, suggesting accretion of a massive disk onto a newborn black hole: GRBs are the birth cries of these black holes. Two of the most promising models: Neutron Star Mergers and Collapars lead naturally to this scenario.

We determine the distribution of total energy emitted by gamma-ray bursts for bursts with fluences and distance information. Our core sample consists of eight bursts with BATSE spectra and spectroscopic redshifts. We extend this sample by adding four bursts with BATSE spectra and host galaxy R magnitudes. From these R magnitudes we calculate a redshift probability distribution; this method requires a model of the host galaxy population. From a sample of ten bursts with both spectroscopic redshifts and host galaxy R magnitudes (some do not have BATSE spectra) we find that the burst rate is proportional to the galaxy luminosity at the epoch of the burst. Assuming that the total energy emitted has a log-normal distribution, we find that the average emitted energy (assumed to be radiated isotropically) is $<E_{\gamma iso} > = 1.3^{+1.2}_{-1.0} \times 10^{53}$ ergs (for H$_0$ = 65 km s$^{-1}$ Mpc$^{-1}$, $\Omega_m=0.3$ and $\Omega_\Lambda=0.7$); the distribution has a logarithmic width of $\sigma_\gamma=1.7^{+0.7}_{-0.3}$. The corresponding distribution of X-ray afterglow energy (for seven bursts) has $<E_{X iso} > = 4.0^{+1.6}_{-1.8} \times 10^{51}$ergs and $\sigma_X = 1.3^{+0.4}_{-0.3}$. For completeness, we also provide spectral fits for all bursts with BATSE spectra for which there were afterglow searches.

We perform fully relativistic hydrodynamic simulations of the deceleration and lateral expansion of a relativistic jet as it expands into an ambient medium. The hydrodynamic calculations use a 2D adaptive mesh refinement (AMR) code, which provides adequate resolution of the thin shell of matter behind the shock. We find that the sideways propagation is different than predicted by simple analytic models. The physical conditions at the sides of the jet are found to be significantly different than at the front of the jet, and most of the emission occurs within the initial opening angle of the jet. The light curves, as seen by observers at different viewing angles with respect to the jet axis, are then calculated assuming synchrotron emission. For an observer along the jet axis, we find a sharp achromatic `jet break' in the light curve at frequencies above the typical synchrotron frequency, at $t_{jet}\approx 5.8(E_{52}/n_1)^{1/3}(\theta_0/0.2)^{8/3}$ days, while the temporal decay index $\alpha$ ($F_{\nu}\propto t^{\alpha}$) after the break is steeper than $-p$ ($\alpha=-2.85$ for $p=2.5$). At larger viewing angles $t_{jet}$ increases and the jet break becomes smoother.

We model the intermediate time evolution of a "jetted" gamma-ray burst by two blobs of matter colliding with the interstellar medium. We follow the hydrodynamical evolution of this system numerically and calculate the bremsstrahlung and synchrotron images of the remnant. We find that for a burst energy of $10^{51}$ erg the remnant becomes spherical after $\sim 5000$ years when it collects $\sim 50M_\odot$ of interstellar mass. This result is independent of the exact details of the GRB, such as the opening angle. After this time a gamma-ray burst remnant has an expanding sphere morphology. The similarity to a supernova remnant makes it difficult distinguish between the two at this stage. The expected number of non-spherical gamma-ray burst remnants is $\sim0.05$ per galaxy for a beaming factor of 0.01 and a burst energy of $10^{51}$ erg. Our results suggest that that the double-shell object DEM L 316 is not a GRB remnant.

We analyze the temporal structure of long ((T_{90}>2sec)) and short ((T_{90}<2sec)) BATSE bursts. We find that: (i) In many short bursts (\delta t_{min}/T\ll 1) (where (\delta t_{min}) is the shortest pulse). This indicates that short bursts arise, like long ones, in internal shocks. (ii) In long bursts there is an excess of long intervals between pulses (relative to a lognormal distribution). This excess can be explained by the existence of \emph{quiescent times}, long periods with no signal above the background that arise, most likely, from periods with no source activity. The lognormal distribution of the intervals (excluding the \emph{quiescent times}) is similar and correlated with the distribution of the pulses width, in agreement with the predictions of the internal shock model.

Gamma Ray Bursts (GRBs) are the most relativistic objects known so far, involving, on one hand an ultra-relativistic motion with a Lorentz factor $\Gamma > 100$ and on the other hand an accreting newborn black hole. The two main routes leading to this scenario: binary neutron star mergers and Collapsar - the collapse of a rotating star to a black hole, are classical sources for gravitational radiation. Additionally one expect a specific a gravitational radiation pulse associated with the acceleration of the relativistic ejecta. I consider here the implication of the observed rates of GRBs to the possibility of detection of a gravitational radiation signal associated with a GRB. Unfortunately I find that, with currently planned detectors it is impossible to detect the direct gravitational radiation associated with the GRB. It is also quite unlikely to detect gravitational radiation associated with Collapsars. However, the detection of gravitational radiation from a neutron star merger associated with a GRB is likely.

We study the long term evolution of a solar type star that is being disrupted by a super massive (10^6 solar mass) black hole. The evolution is followed from the disruption event, which turns the star into a long thin stream of gas, to the point where some of this gas returns to pericenter and begins its second orbit around the black hole. Following the evolution for this long allows us to determine the amount of mass that is accreted by the black hole. We find that approximately 75% of the returning mass is not accreted but instead becomes unbound, following the large compression characterizing the return to pericenter. The impact of a tidal disruption on the surrounding gas may therefore be like that of two consecutive supernove-type events.

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

The apparent lack of suitable astrophysical sources for cosmic rays with E > 10^{19.7} eV (UHECRs) is the "GZK Paradox". We argue that whatever mechanism produces them must also account for events down to ~10^{18.7} eV, including their isotropy and spectral smoothness. This rules out galactic sources, distributed sources such as topological defects, and Gamma Ray Bursts (GRBs). We are lead to identify the powerful radio galaxy Cen A, at 3.4 Mpc, as the probable source of most UHECRs observed at Earth today, and to estimate the extragalactic magnetic field to be ~0.3 microG.

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

We introduce a simple model for the formation of voids. In this model the underdensity of galaxies in voids is the product of two factors. The first arises from a gravitational expansion of the negative density perturbation. The second is due to biasing: galaxies are less likely to form in an underdense region. One feature of the model is an upper cutoff in void sizes. We calculate the volume filling factor of characteristic voids for different CDM models and find that our formation model points to $\Lambda$CDM models as preferred models of the power spectrum. A natural consequence of our model is that the underdensity of the dark matter inside voids is smaller than the galaxy underdensity.

The hard X-ray (>2 keV) emission of the local and distant Universe as observed with the HEAO1-A2 experiment is reconsidered in the context of large scale cosmic structure. Using all-sky X-ray samples of AGN and galaxy clusters we remove the dominant local X-ray flux from within a redshift of ~ 0.02. We evaluate the dipolar and higher order harmonic structure in 4 X-ray colours. The estimated dipole anisotropy of the unresolved flux appears to be consistent with a combination of the Compton-Getting effect due to the Local Group motion (dipole amplitude Delta = 0.0042) and remaining large scale structure (0.0023 <~ Delta <~ 0.0085), in good agreement with the expectations of Cold Dark Matter models. The observed anisotropy does however also suggest a non-negligible Galactic contribution which is more complex than current, simple models of >2 keV Galactic X-ray emission. Comparison of the soft and hard colour maps with a harmonic analysis of the 1.5 keV ROSAT all-sky data qualitatively suggests that at least a third of the faint, unresolved ~ 18 deg scale structure in the HEAO1-A2 data may be Galactic in origin. However, the effect on measured flux dipoles is small (<~3%). We derive an expression for dipole anisotropy and acceleration and demonstrate how the dipole anisotropy of the distant X-ray frame can constrain the amplitude of bulk motions of the universe. From observed bulk motions over a local ~ 50 Mpc/h radius volume we determine 0.14 <~ Omega^0.6/b_x(0) <~ 0.59.

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

We introduce a set of four new publicly available N-body simulations, the most recent additions to the Texas P^3M Database. Our models probe the less studied parameter space region of moderate volume (100/h Mpc box) combined with fine mass resolution (~10^12 M_Sol, roughly comparable to a L_* galaxy), making these simulations especially suitable for study of major large-scale structure features such as voids, and for comparison with the largest three-dimensional redshift surveys currently available. Our cosmological models (LCDM, TOCDM, OCDM, TCDM) are all COBE-normalized, and when possible (LCDM and TOCDM) also cluster-normalized, based on the X-ray cluster M--T relation. The COBE- and cluster-normalized LCDM model reiterates the attractiveness of this currently favored model which does not require the introduction of tilt in order to fit the constraints imposed by observations of other cosmological parameters.

Results of 3D, Newtonian hydrodynamic calculations of the last stages of the inspiral and the final coalescence of neutron star binary systems are described. The focus is on sligthly asymmetric systems with the asymmetry stemming either from different spins or different masses of the binary components. Implications for nucleosynthesis and gamma ray bursts are discussed. Test calculations for the used viscosity scheme are provided in the appendix.

In the internal shock scenario for GRBs we expect some fraction of the energy of the burst to be carried by slow moving shells that were ejected at late times. These slow shells collide with faster moving outer shells when the outer shells have slowed down as a result of sweeping up material from the ISM. This gives rise to a forward shock that moves into the outer shell producing a bump in the afterglow light curve of amplitude roughly proportional to the ratio of the energy in the inner and the outer shells. In addition, a reverse shock propagates in the inner shell and produces emission at a characteristic frequency that is typically much smaller than the peak of the emission from the outer shell by a factor of $\sim 7 \gamma_{0c}^2 (E_2/E_1)^{1.1}$, and the observed flux at this frequency from the reverse shock is larger compared to the flux from the outer shell by a factor of $\sim 8 (\gamma_{0c} E_2/E_1)^{5/3} $; where $\gamma_{0c}$ is the bulk Lorentz factor of the outer shell at the time of collision, and $E_1 & E_2$ are the total energy in the outer and the inner shells respectively. The Lorentz factor is related to the observer time as $\sim 5 (t/day)^{3/8}$. The shell collision could produce initial temporal variability in the early afterglow signal. The lack of significant deviation from a power-law decline of the optical afterglow from half a dozen bursts suggests that $E_2/E_1$ is small. Future multi-wavelength observations should be able to either detect bumps in the light curve corresponding to both the forward and the reverse shocks or further constrain the late time release of energy in ejecta with small Lorentz factor, which is expected generically in the internal shock models for the gamma-ray bursts.

We discuss the spectrum arising from synchrotron emission by fast cooling (FC) electrons, when fresh electrons are continually accelerated by a strong blast wave, into a power law distribution of energies. The FC spectrum was so far described by four power law segments divided by three break frequencies $\nu_{sa} < \nu_{c} < \nu_{m}$. This is valid for a homogeneous electron distribution. However, hot electrons are located right after the shock, while most electrons are farther down stream and have cooled. This spatial distribution changes the optically thick part of the spectrum, introducing a new break frequency, $\nu_{ac} < \nu_{sa}$, and a new spectral slope, $F_{\nu }\propto \nu^{11/8}$ for $\nu_{ac} < \nu < \nu_{sa}$. The familiar $F_{\nu }\propto \nu ^{2}$ holds only for $\nu < \nu_{ac}$. This ordering of the break frequencies is relevant for typical gamma-ray burst (GRB) afterglows in an ISM environment. Other possibilities arise for internal shocks or afterglows in dense circumstellar winds. We discuss possible implications of this spectrum for GRBs and their afterglows, in the context of the internal-external shock model. Observations of $F_{\nu} \propto \nu^{11/8}$ would enable us to probe scales much smaller than the typical size of the system, and constrain the amount of turbulent mixing behind the shock.

The apparent lack of suitable astrophysical sources for the observed highest energy cosmic rays (UHECRs) within $\approx 20$ Mpc is the "GZK Paradox". We constrain representative models of the extra-galactic magnetic field structure by Faraday Rotation measurements; limits are at the $\mu$G level rather than the nG level usually assumed. In such fields, even the highest energy cosmic rays experience large deflections. This allows nearby AGNs (possibly quiet today) or GRBs to be the source of ultra-high energy cosmic rays without contradicting the GZK distance limit.

We describe fully relativistic three dimensional calculations of the slowing down and spreading of a relativistic jet by an external medium like the ISM. We calculate the synchrotron spectra and light curves using the conditions determined by the hydrodynamic calculations. Preliminary results with a moderate resolution are presented here. Higher resolution calculations are in progress.

Gamma-ray bursts are believed to be some catastrophic event in which material is ejected at a relativistic velocity, and internal collisions within this ejecta produce the observed $\gamma$-ray flash. The angular size of a causally connected region within a relativistic flow is of the order the angular width of the relativistic beaming, $\gamma^{-1}$. Thus, different observers along different lines of sights could see drastically different fluxes from the same burst. Specifically, we propose that the most energetic bursts correspond to exceptionally bright spots along the line of sight on colliding shells, and do not represent much larger energy release in the explosion. We calculate the distribution function of the observed fluence for random angular-fluctuation of ejecta. We find that the width of the distribution function for the observed fluence is about two orders of magnitude if the number of shells ejected along different lines of sight is ten or less. The distribution function becomes narrower if number of shells along typical lines of sight increases. The analysis of the $\gamma$-ray fluence and afterglow emissions for GRBs with known redshifts provides support for our model i.e. the large width of GRB luminosity function is not due to a large spread in the energy release but instead is due to large angular fluctuations in ejected material. We outline several observational tests of this model. In particular, we predict little correlation between the $\gamma$-ray fluence and the afterglow emission as in fact is observed. We predict that the early (minutes to hours) afterglow would depict large temporal fluctuations whose amplitude decreases with time. Finally we predict that there should be many weak bursts with about average afterglow luminosity in this scenario.

We present a comparison between the voids in two nearly all sky redshift surveys: the ORS and the IRAS 1.2 Jy. While the galaxies in these surveys are selected differently and their populations are known to be biased relative to each other, the two void distributions are similar. We compare the spatial distribution of the two void populations and demonstrate the correlation between them. The voids also agree with regard to the overall void statistics -- a filling factor of ~0.45 of the volume, an average void diameter ~45/h Mpc and an average galaxy underdensity in the voids of approximately -0.9. Our measurements of the underdensities of the voids in the two surveys enable us to estimate the relative bias in the voids between optical and IRAS samples. We find b(opt) / b(IRAS) ~ 1, showing that on average there is little -- or no -- biasing between the two void populations.

For a few seconds a gamma-ray burst (GRB) becomes the brightest object in the Universe, over-shining the rest of the Universe combined! Clearly this reflects extreme conditions that are fascinating and worth exploring. The recent discovery of GRB afterglow have demonstrated that we are on the right track towards the resolution of this long standing puzzle. These observations have confirmed the relativistic fireball model (more specifically the internal-external shocks model). The prompt optical emission seen in GRB 990123 have demonstrated that GRBs involve ultra-relativistic motion. The breaks in the light curves of GRB 990123 and GRB 990510 and the peculiar light curves of GRB 980519 and GRB 980326 disclosed that these GRBs are beamed. I examine these recent developments and discuss their implications to the models of the source. I argue that the current understanding implies that GRBs signal the birth of stellar mass black holes.

GRB 980519 had the most rapidly fading of the well-documented GRB afterglows, consistent with t^{-2.05 +/- 0.04} in BVRI as well as in X-rays during the two days in which observations were made. We report VRI observations from the MDM 1.3m and WIYN 3.5m telescopes, and we synthesize an optical spectrum from all of the available photometry. The optical spectrum alone is well fitted by a power law of the form nu^{-1.20 +/- 0.25}, with some of the uncertainty due to the significant Galactic reddening in this direction. The optical and X-ray spectra together are adequately fitted by a single power law nu^{-1.05 +/- 0.10}. This combination of steep temporal decay and flat broad-band spectrum places a severe strain on the simplest afterglow models involving spherical blast waves in a homogeneous medium. Instead, the rapid observed temporal decay is more consistent with models of expansion into a medium of density n(r) proportional to r^{-2}, or with predictions of the evolution of a jet after it slows down and spreads laterally. The jet model would relax the energy requirements on some of the more extreme GRBs, of which GRB 980519 is likely to be an example because of its large gamma-ray fluence and faint host galaxy.

In several GRBs afterglows, rapid temporal decay is observed which is inconsistent with spherical (isotropic) blast-wave models. In particular, GRB 980519 had the most rapidly fading of the well-documented GRB afterglows, with t^{-2.05\pm 0.04} in optical as well as in X-rays. We show that such temporal decay is more consistent with the evolution of a jet after it slows down and spreads laterally, for which t^{-p} decay is expected (where p is the index of the electron energy distribution). Such a beaming model would relax the energy requirements on some of the more extreme GRBs by a factor of several hundreds. It is likely that a large fraction of the weak (or no) afterglow observations are also due to the common occurrence of beaming in GRBs, and that their jets have already transitioned to the spreading phase before the first afterglow observations were made. With this interpretation, a universal value of p~2.5 is consistent with all data.

Afterglow, or long-lived emission, has now been detected from about a dozen well-positioned gamma-ray bursts. Distance determinations made by measuring optical emission lines from the host galaxy, or absorption lines in the afterglow spectrum, place the burst sources at significant cosmological distances, with redshifts ranging from ~1--3. The energy required to produce the bright gamma-ray flashes is enormous: up to ~10^{53} erg or 10 percent of the rest mass energy of a neutron star, if the emission is isotropic. Here we present the discovery of the optical afterglow and the redshift of GRB 990123, the brightest well-localized GRB to date. With our measured redshift of >1.6, the inferred isotropic energy release exceeds the rest mass of a neutron star thereby challenging current theoretical models for the origin of GRBs. We argue that the optical and IR afterglow measurements reported here may provide the first observational evidence of beaming in a GRB, thereby reducing the required energetics to a level where stellar death models are still tenable.

We compare the ongoing observations of the remarkable burst GRB990123, the mother of all bursts, with the predictions of the afterglow theory. We show that the observations agree with the recent prediction that a reverse shock propagating into the ejecta would produce a very strong prompt optical flash. This reverse shock has also produced the 8.46GHz radio signal, observed after one day. The forward shock, which propagates into the ISM is the origin of the classical afterglow. It has produced the prompt X-ray signal as well as the late optical and IR emission. It would most likely produce a radio emission within the next few weeks. The observations suggest that the initial Lorentz factor of the ejecta was $\sim 200$. Within factors of order unity, this crude model explains all current observations of GRB990123.

According to the internal-external shocks model for $\gamma $-ray bursts (GRBs), the GRB is produced by internal shocks within a relativistic flow while the afterglow is produced by external shocks with the ISM. We explore the early afterglow emission. For short GRBs the peak of the afterglow will be delayed, typically, by few dozens of seconds after the burst. For long GRBs the early afterglow emission will overlap the GRB signal. We calculate the expected spectrum and the light curves of the early afterglow in the optical, X-ray and $\gamma $-ray bands. These characteristics provide a way to discriminate between late internal shocks emission (part of the GRB) and the early afterglow signal. If such a delayed emission, with the characteristics of the early afterglow, will be detected it can be used both to prove the internal shock scenario as producing the GRB, as well as to measure the initial Lorentz factor of the relativistic flow. The reverse shock, at its peak, contains energy which is comparable to that of the GRB itself, but has a much lower temperature than that of the forward shock so it radiates at considerably lower frequencies. The reverse shock dominates the early optical emission, and an optical flash brighter than 15th magnitude, is expected together with the forward shock peak at x-rays or $\gamma$-rays. If this optical flash is not observed, strong limitations can be put on the baryonic contents of the relativistic shell deriving the GRBs, leading to a magnetically dominated energy density.

We calculate the expected spectrum and light curves of the early afterglow. For short GRBs the peak of the afterglow will be delayed, typically, by a few dozens of seconds after the burst. The x-ray and gamma-ray characteristics of this delayed emission provide a way to discriminate between late internal shocks emission (part of the GRB) and the early afterglow signal. Detection of this delayed emission will prove the internal shock scenario as producing the GRB, and will pinpoint the initial Lorentz factor gamma_0. In the optical band, the dominant emission arises from the reverse shock. This shock, carries an energy comparable to that of the forward shock. It radiates this energy at much lower frequencies, producing a short optical flash of 15th magnitude or brighter.

We calculate the light curve and spectra near the peak and the self absorption break, for an adiabatic blast wave described by the Blandford-McKee solution, considering the emission from the whole region behind the shock front. The expected light curve and spectra are flat near the peak. This rules out the interpretation of the sharp peak observed in the optical afterglow of GRB970508 as the expected peak of the light curve. The observed image of an afterglow is calculated for a broad range of frequencies. We show that for frequencies below the self absorption frequency the image is rather homogeneous, as opposed to the bright ring at the outer edge and dim center, which appear at higher frequencies. We fit the observed spectra of GRB970508 to the detailed theory and obtain estimates of the physical parameters of this burst.

We present the results of 3D Newtonian SPH simulations of the merger of a neutron star binary. The microscopic properties of matter are described by the physical equation of state of Lattimer and Swesty (LS-EOS). To test for the robustness of our results we check the sensitivity to the approximations of our model as well as to the binary system parameters. The main and new result is that for the realistic LS-EOS, depending on the initial spin, between 4e-3 and 4e-2 solar masses of material become unbound. If, as suggested, large parts of this matter consist of r-process nuclei, neutron star mergers could account for the whole observed r-process material in the Galaxy.

Gamma-ray bursts (GRBs) have puzzled astronomers since their accidental discovery in the late sixties. The BATSE detector on the COMPTON-GRO satellite has been detecting one burst per day for the last six years. Its findings have revolutionized our ideas about the nature of these objects. They have shown that GRBs are at cosmological distances. This idea was accepted with difficulties at first. The recent discovery of an X-ray afterglow by the Italian/Dutch satellite BeppoSAX has led to a detection of high red-shift absorption lines in the optical afterglow of GRB970508 and in several other bursts and to the identification of host galaxies to others. This has confirmed the cosmological origin. Cosmological GRBs release $\sim 10^{51}-10^{53}$ergs in a few seconds making them the most (electromagnetically) luminous objects in the Universe. The simplest, most conventional, and practically inevitable, interpretation of these observations is that GRBs result from the conversion of the kinetic energy of ultra-relativistic particles or possibly the electromagnetic energy of a Poynting flux to radiation in an optically thin region. This generic "fireball" model has also been confirmed by the afterglow observations. The "inner engine" that accelerates the relativistic flow is hidden from direct observations. Consequently it is difficult to infer its structure directly from current observations. Recent studies show, however, that this ``inner engine'' is responsible for the complicated temporal structure observed in GRBs. This temporal structure and energy considerations indicates that the ``inner engine'' is associated with the formation of a compact object - most likely a black hole.

We use general arguments to show that a continuously powered radiative blast wave can behave self similarly if the energy injection and radiation mechanisms are self similar. In that case, the power-law indices of the blast wave evolution are set by only one of the two constituent physical mechanisms. If the luminosity of the energy source drops fast enough, the radiation mechanisms set the power-law indices, otherwise, they are set by the behavior of the energy source itself. We obtain self similar solutions for the Newtonian and the ultra-relativistic limits. Both limits behave self similarly if we assume that the central source supplies energy in the form of a hot wind, and that the radiative mechanism is the semi-radiative mechanism of Cohen, Piran & Sari (1998). We calculate the instantaneous radiative efficiencies for both limits and find that a relativistic blast wave has a higher efficiency than a Newtonian one. The instantaneous radiative efficiency depends strongly on the hydrodynamics and cannot be approximated by an estimate of local microscopic radiative efficiencies, since a fraction of the injected energy is deposited in shocked matter. These solutions can be used to calculate Gamma Ray Bursts afterglows, for cases in which the energy is not supplied instantaneously.

GRB afterglow is reasonably described by synchrotron emission from relativistic blast waves at cosmological distances. We perform detailed calculations taking into account the effect of synchrotron self absorption. We consider emission from the whole region behind the shock front, and use the Blandford McKee self similar solution to describe the fluid behind the shock. We calculate the spectra and the observed image of a GRB afterglow near the self absorption frequency $\nu_a$ and derive an accurate expression for $\nu_a$. We show that the image is rather homogeneous for $\nu<\nu_a$, as opposed to the bright ring at the outer edge and dim center, which appear at higher frequencies. We compare the spectra we obtain to radio observations of GRB970508. We combine the calculations of the spectra near the self absorption frequency with other parts of the spectra and obtain revised estimates for the physical parameters of the burst: $E_{52}=0.53$, $\epsilon_e=0.57$, $\epsilon_B=0.0082$, $n_1=5.3$. These estimates are different by up to two orders of magnitude than the estimates based on an approximate spectrum.

Gamma-ray bursts GRBs are among the most mysterious astronomical phenomenon ever discovered. Unlike most astronomical discoveries which were explained within weeks or months after their initial discovery, GRBs remain a puzzle for more than thirty years. During the last decade our understanding of GRBs has undergone two major revolutions. First, BATSE discovered that GRBs are distributed isotropically over the sky and thereby demonstrated their cosmological origin. The second revolution tool place more recently when BeppoSAX discovered GRB afterglow. This confirmed the fireball model and led to a wealth of observational data, some of which has not been fully understood yet. The emerging picture is that GRBs are the most luminous objects and the most relativistic objects ever discovered: (i) GRBs involve relativistic motion at a velocity of 0.9999c or larger. (ii) Most current GRB models involve the formation of a black hole in one way or another. (iii) If binary neutron star mergers are the sources of GRBs then GRBs are also associated with gravitational radiation signals. Finally, (iv) as cosmological power-houses that are observed to high red-shift GRBs might be used to measure cosmological parameters and to teach us about the epoch of galaxy formation.

We present an attempt to measure the large angular scale fluctuations in the X-Ray Background (XRB) from the HEAO1-A2 data, expressed in terms of spherical harmonics. We model the harmonic coefficients assuming a power spectrum and an epoch-dependent bias parameter, and using a phenomenological scenario describing the evolution of the X-ray sources. From the few low-order multipoles detected above shot noise, we estimate the power-spectrum normalization on scales intermediate between those explored by local galaxy redshift surveys (~ 100 Mpc) and by the COBE Microwave Background measurements (~ 1000 Mpc). We find that the HEAO1 harmonics are consistent with present epoch rms fluctuations of the X-ray sources bx(0)sigma8 ~ 1-2 in 8 Mpc spheres. Therefore the observed fluctuations in the XRB are roughly as expected from interpolating between the local galaxy surveys and the COBE CMB experiment. We predict that an X-ray all-sky surface brightness survey resolving sources a factor of 10 fainter than HEAO1, may reveal fluctuations to significantly larger angular scales and therefore more strongly constrain the large scale structure of the Universe on scales of hundreds of Mpcs.

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

Recent calculations place questions on the ability of supernovae to produce r-process nuclei in the correct amounts. We present results from 3D Newtonian SPH calculations of mergers of equal mass neutron star binaries. We find the amounts of ejected mass to be very promising for an explanation of the observed r-process abundances. Preliminary abundance calculations show that practically all the material is subject to r-process conditions. The calculated abundance patterns reproduce the basic features of the solar r-process abundances very well. Especially the peak around A= 195 is easily reproduced without any tuning of the initial entropies. Thus, eventually all the observed r-process material could be explained by mass ejection during neutron star mergers.

We derive a self similar solution for the propagation of an extreme relativistic (or Newtonian) radiative spherical blast wave into a surrounding cold medium. The solution is obtained under the assumption that the radiation process is fast, it takes place only in the vicinity of the shock and that it radiates away a fixed fraction of the energy generated by the shock. In the Newtonian regime these solutions generalize the Sedov-Taylor adiabatic solution and the pressure-driven fully radiative solution. In the extreme relativistic case these solutions generalize the Blandford-McKee adiabatic solution. They provide a new fully radiative extreme relativistic solution which is different from the Blandford-McKee fully radiative relativistic solution. This new solution develops a hot interior which causes it to cool faster than previous estimates. We find that the energy of the blast wave behaves as a power law of the location of the shock. The power law index depends on the fraction of the energy emitted by the shock. We obtain an analytic solution for the interior of the blast wave. These new solutions might be applicable to the study of GRB afterglow or SNRs.

We study numerically the evolution of an adiabatic relativistic fireball expanding into a cold uniform medium. We follow the stages of initial free expansion and acceleration, coasting and then deceleration and slowing down to a non-relativistic velocity. We compare the numerical results with simplified analytical estimates. We show that the relativistic self similar Blandford-McKee solution describes well the relativistic deceleration epoch. It is an excellent approximation throughout the relativistic deceleration stage, down to $\gamma \sim 5$, and a reasonable approximation even down to $\gamma \sim 2$ though the solution is rigorous only for $\gamma \gg 1$. We examine the transition into the Blandford-McKee solution, and the transition from the solution to the non-relativistic self similar Sedov-Taylor solution. These simulations demonstrate the attractive nature of the Blandford-McKee solution and its stability to radial perturbations.

It is often stated that gamma-ray bursts (GRBs) have typical energies of several hundreds $\keV$, where the typical energy may be characterized by the hardness H, the photon energy corresponding to the peak of $\nu F_{\nu}$. Among the 54 BATSE bursts analyzed by Band et al. (1993), and 136 analyzed by us, more then 60% have 50 keV < H < 300 keV. Is the narrow range of H a real feature of GRBs or is it due to an observational difficulty to detect harder bursts? We consider a population of standard candle bursts with a hardness distribution: rho(H) d log H \propto H^gamma d log H and no luminosity - hardness correlation. We model the detection algorithm of BATSE as a function of H, including cosmological effects, detector characteristics and triggering procedure, and we calculate the expected distribution of H in the observed sample for various values of gamma. Both samples shows a paucity of soft (X-ray) bursts, which may be real. However, we find that the observed samples are consistent with a distribution above H=120 keV with gamma \sim -0.5 (a slowly decreasing numbers of GRBs per decade of hardness). Thus, we suggest that a large population of unobserved hard gamma-ray bursts may exist.

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

The discovery by BeppoSAX and coordinated ground-based observations of persistent X-ray, visible and radio counterparts to GRB has successfully concluded a search begin in 1973. The problem of explaining the mechanisms of GRB and their persistent counterparts remains. There are two classes of models: 1) GRB continue weakly for days at all frequencies; 2) GRB emission shifts to lower frequencies as relativistic debris sweeps up surrounding gas (in an ``external shock'') and slows. In 1) the visible afterglow is accompanied by continuing gamma-ray emission, as hinted by the high energy emission of GRB940217 and the ``Gang of Four'' bursts of October 27--29, 1996; the persistent emission may fluctuate. Behavior of this sort may be found in ``internal shock'' models. Models of class 2) have been the subject of several theoretical studies which disagree in assumptions and details but which predict that at each frequency the flux should rise and then decline, with the maximum coming later at lower frequencies. Some of this behavior has been observed, but data from GRB970508 show that its afterglow cannot be simply extrapolated from its gamma-ray emission. It is likely that both classes of processes occur in most GRB. Comparisons between GRB show that they are not all scaled versions of the same event. These results suggest that most gamma-ray emission is the result of ``internal shocks'' while most afterglow is the result of ``external shocks'', and hint at the presence of collimated outflows. Self-absorption in the radio spectrum of GRB970508 permitted the size of the radiating surface to be estimated, and in future GRB it may be possible to follow the expansion of the shell in detail and to construct an energy budget.

The recently discovered GRB afterglow is believed to be described reasonably well by synchrotron emission from a slowing down relativistic shell that collides with an external medium. To compare theoretical models with afterglow observations we calculate here the broad band spectrum and corresponding light curve of synchrotron radiation from a power-law distribution of electrons in an expanding relativistic shock. Both the spectrum and the light curve consist of several power-law segments. The light curve is constructed under two limiting models for the hydrodynamical evolution of the shock: fully adiabatic and fully radiative. We compare the results with observations of $\gamma$-ray burst afterglows.

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

The recent discoveries of X-ray and optical counterparts for GRBs, and a possible discovery of a host galaxy, implies that a direct measurement of the red-shift of some GRBs host galaxies is eminent. We discuss the implications of such measurements. They could enable us to determine the GRBs luminosity distribution, the variation of the rate of GRBs with cosmic time, and even, under favorable circumstances, to estimate $\Omega$. Using GRB970508 alone, and assuming standard candles, we constrain the intrinsic GRB evolution to $\rho(z)=(1+z)^{-0.5\pm0.7}$.

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

The gamma-ray burst source GRB970228 has been observed after a delay of 8--12 hours in X-rays and after one day in visible and near infrared light. This marks the first detection of emission at lower frequencies following the gamma-ray observation of a GRB and the first detection of any visible counterpart to a GRB. We consider possible delayed visible and X-ray emission mechanisms, and conclude that the intrinsic gamma-ray activity continued at a much reduced intensity for at least a day. There are hints of such continued activity in other GRB, and future observations can decide if this is true of GRB in general. The observed multi-band spectrum of GRB970228 agrees with the predictions of relativistic shock theory when the flux is integrated over a time longer than that required for a radiating electron to lose its energy.

Relativistic shock models of gamma-ray bursts may be tested by comparison of their predicted low energy asymptotic spectral indices s to observations. Synchrotron radiation theory predicts that the instantaneous spectrum has s = - 1/3 and the spectrum integrated over the radiative decay of the electrons' energies has s = 1/2, with other cases lying between these limits. We examine the spectra of 11 bursts obtained by the Large Area Detectors on BATSE. One agrees with the predicted instantaneous spectrum, as does the initial portion of a second, and three are close to the predicted integrated spectrum. All the observed asymptotic spectral slopes lie in the predicted range. This evidence for relativistic shocks is independent of detailed models of bursts and of assumptions about their distances. Radiation observed with the predicted instantaneous spectrum has a comparatively smooth time dependence, consistent with the necessarily long radiation time, while that with the predicted integrated spectrum has a spiky time dependence, consistent with the necessarily short radiation time.

Voids are the most prominent feature of the LSS of the universe. Still, they have been generally ignored in quantitative analysis of it, essentially due to the lack of an objective tool to identify and quantify the voids. To overcome this, we present the Void-Finder algorithm, a novel tool for objectively quantifying galaxy voids. The algorithm classifies galaxies as either wall- or field-galaxies. Then it identifies voids in the wall-galaxy distribution. Voids are defined as continuous volumes that do not contain any wall-galaxies. The voids must be thicker than an adjustable limit, which is refined in successive iterations. We test the algorithm using Voronoi tessellations. By appropriate scaling of the parameters we apply it to the SSRS2 survey and to the IRAS 1.2 Jy. Both surveys show similar properties: ~50% of the volume is filled by the voids, which have a scale of at least 40 Mpc, and a -0.9 under-density. Faint galaxies populate the voids more than bright ones. These results suggest that both optically and IRAS selected galaxies delineate the same LSS. Comparison with the recovered mass distribution further suggests that the observed voids in the galaxy distribution correspond well to under-dense regions in the mass distribution. This confirms the gravitational origin of the voids.

The highly variable temporal structure observed in most GRBs provides us with unexpected clues. We show that variable GRBs cannot be produced by external shocks models and consequently cannot be produced by an ``explosive'' inner engine. The observed temporal structure must reflect the activity of the inner engine that must be producing unsteady and irregular ``wind'' which is converted to radiation via internal shocks.

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

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

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

We have laid down the formalism and techniques necessary for computing the multipole components in a spherical harmonic expansion for bursting sources given any specific power spectrum of density perturbations. Using this formalism we have explicitly computed and tabulated the expected first few multipole components for the GRB distribution using a Cold Dark Matter power spectrum. Unfortunately, our analysis leads us to expect that an anisotropy signal for this model of structure formation will be below the shot noise level for the foreseeable future. Earlier studies by others had claimed that it may be possible within a decade or so to probe structure formation using the angular distribution of GRBs. We find that while it may be possible to probe the dipole due to our motion with respect to the cosmic rest-frame, it does not seem feasible to probe the multipole components due to intrinsic fluctuations if the actual power spectrum is fixed to within an order of magnitude by the Cold Dark Matter model of structure formation.

Gamma ray Bursts (GRBs) - short bursts of few hundred keV $\gamma$-rays - have fascinated astronomers since their accidental discovery in the sixties. GRBs were ignored by most relativists who did not expect that they are associated with any relativistic phenomenon. The recent observations of the BATSE detector on the Compton GRO satellite have revolutionized our ideas on these bursts and the picture that emerges shows that GRBs are the most relativistic objects discovered so far.

We find that Gamma-ray bursts (GRBs) at the 3B catalog are correlated (at 95\% confidence level) with Abell clusters. This is the first known correlation of GRBs with any other astronomical population. It confirms the cosmological origin of GRBs. Comparison of the rich clusters auto-correlation with the cross-correlation found here suggests that \sim 26 \pm 15\% of an accurate (\delta<2.3 deg.) position GRBs sub-sample members are located within 600 \hmpc.

A leading mechanism for producing cosmological gamma-ray bursts (GRBs) is via ultra-relativistic particles in an expanding fireball. The kinetic energy of the particles is converted into thermal energy in a forward shock and a reverse shock, when the outward flowing particles encounter the interstellar medium. The thermal energy is radiated via synchrotron emission and Comptonization. We estimate the synchrotron cooling time scale of the shocked regions for electrons of various Lorentz factors, focusing in particular on those electrons whose radiation falls within the energy detection range of the BATSE detectors. We find that in order to produce the rapid variability observed in most bursts the magnetic energy density in the shocked material must be greater than about 1\% of the thermal energy density. Additionally, the electrons must be nearly in equipartition with the protons, otherwise the radiative efficiency is low. Inverse Compton scattering can increase the cooling rate of the relevant electrons but the Comptonized emission itself is never within the BATSE range. These arguments allow us to pinpoint the conditions within the radiating regions in GRBs and to determine the important radiation processes. In addition, the model predicts that the duty cycle of intensity variations in GRB light curves should be nearly independent of burst duration, and should scale inversely as the square root of the observed photon energy. Both correlations are in agreement with observations. The model also provides a plausible explanation for the bimodal distribution of burst durations. There is no explanation, however, for the characteristic break energy in GRB spectra.

We analyze the anisotropy signature expected if the high energy (above $10^{19}$eV) cosmic ray (CR) sources are extra-Galactic and trace the distribution of luminous matter on large scales. We investigate the dependence of the anisotropy on both the relative bias between the CR sources and the galaxy distribution and on the (unknown) intrinsic CR source density. We find that the expected anisotropy associated with the large scale structure (LSS) should be detected once the number of CR events observed above $10^{19}{\rm eV}$ is increased by a factor of $\sim10$. This would require $\sim30$ observation-years with existing experiments, but less then $1$ year with the proposed $\sim5000\ {\rm km}^2$ Auger detectors. We find that the recently reported concentration of the Haverah Park CR events towards the super-galactic plane is not consistent with the known LSS. If real, the Haverah Park result suggests that the CR sources are much more concentrated towards the super-galactic plane than the known LSS. Our results are not sensitive to the number density of CR sources. We show that once the number of detected events is increased by a factor of $\sim10$, the number density would be strongly constrained by considering the probability for having repeating sources.

We present a new void search algorithm for automated detection of voids in three-dimensional redshift surveys. Based on a model in which the main features of the LSS of the Universe are voids and walls, we classify the galaxies into wall galaxies and field galaxies and we define voids as continuous volumes that are devoid of any wall galaxies. Field galaxies are allowed within the voids. The algorithm makes no assumptions regarding the shapes of the voids and the only constraint that is imposed is that the voids are always thicker than a preset limit, thus eliminating connections between adjacent voids through smal breaches in the walls. By appropriate scaling of the parameters with the selection function this algorithm can be used to analyze flux-limited surveys. We test the algorithm on Voronoi tessellations and apply it to the SSRS2 redshift survey to derive the spectrum of void sizes and other void properties. We find that the average diameter of a void is $37\pm 8 \h$ Mpc. We suggest the usage of this fully automated algorithm to derive a quantitative description of the voids, as another tool in describing the large scale structure of the Universe and for comparison with numerical simulations.

The analysis of the BATSE's count distribution within cosmological models suffers from observational uncertainties due to the variability of the bursts' spectra: when BATSE observes bursts from different redshifts at a fixed energy band it detects photons from different energy bands at the source. This adds a spectral dependence to the count distribution $N(C)$. Similarly variation of the duration as a function of energy at the source complicates the time dilation analysis. We describe here a new statistical formalism that performs the required ``blue shifting" of the count number and the burst duration in a statistical manner. This formalism allows us to perform a combined best fit (maximal likelihood) to the count distribution, $N(C)$, and the duration distribution simultaneously. The outcome of this analysis is a single best fit value for the redshift of the observed bursts.

We calculate the hydrodynamic time scales for a spherical ultra-relativistic shell that is decelerated by the ISM and discuss the possible relations between these time scales and the observed temporal structure in $\gamma$-ray bursts. We suggest that the bursts' duration is related to the deceleration time, the variability is related to the ISM inhomogeneities and precursors are related to internal shocks within the shell. Good agreement can be achieved for these quantities with reasonable, not fined tuned, astrophysical parameters. The difference between Newtonian and relativistic reverse shocks may lead to the observed bimodal distribution of bursts' durations.

It is often stated that gamma-ray bursts (GRBs) have typical energies of several hundred keV. Is this a real feature of GRBs or is it due to an observational bias? We consider the possibility that bursts of a given bolometric luminosity occur with a hardness distribution $p(H)d \log H \propto H^\gamma d \log H$. We model the detection efficiency of BATSE as a function of $H$ and calculate the expected distribution of $H$ in the observed sample for various values of $\gamma$. We show that because the detection efficiency of BATSE falls steeply with increasing $H$, the paucity of hard bursts need not be real. We find that the observed sample is consistent with a distribution above $H = 100$ keV with $\gamma \approx 0$ or even $\gamma =0.5$. Thus, a large population of unobserved hard gamma-ray bursts may exist. It is important to extend the present analysis to a larger sample of BATSE bursts and to include the OSSE and COMPTEL limits. If the full sample is consistent with $\gamma\ \sgreat\ 0$, then it would be interesting to look for MeV bursts in the future.

We calculate the hydrodynamic time scales for a spherical ultra-relativistic shell that is decelerated by the ISM and discuss the possible relations between these time scales and the observed temporal structure in $\gamma$-ray bursts. We suggest that the bursts' duration is related to the deceleration time, the variability is related to the ISM inhomogeneities and precursors are related to internal shocks within the shell. Good agreement can be achieved for these quantities with reasonable, not fined tuned, astrophysical parameters. The difference between Newtonian and relativistic reverse shocks may lead to the observed bimodal distribution of bursts' durations.

$\gamma$-ray bursts (GRBs) have puzzled astronomers since their accidental discovery in the sixties. The BATSE detector on COMPTON-GRO satellite has been detecting GRBs for the last four years at a rate of one burst per day. Its findings has revolutionized our ideas about the nature of these objects. In this lecture I show that the simplest, most conventional and practically inevitable, interpretation of the observations is that GRBs form during the conversion of the kinetic energy of ultra-relativistic particles to radiation. The inner ``engine" that accelerates these particles is well hidden from direct observations and its origin might remain mysterious for a long time.

Following Blumenthal \etl \cite{blu} and Piran \etl \cite{Pi} I discuss a simple model for the gravitational formation of voids. This model enables us to estimate the primordial perturbation required on the scale of voids. A comparison of this condition with the constraints on the spectrum of fluctuations derived from CMBR observations reveals that gravitational growth of large voids require more power on the scale of $\approx 1500$\kms than is generally accepted in CDM models. For a generic power law primordial spectrum $P(k) \propto k^n$ gravitational growth of voids requires $n \gsim 1.25$, which is marginally compatible with the observed CMBR fluctuations. If future measurements will show that the amplitude of fluctuations required to produce the voids is too large they will imply that the voids did not grow gravitationally and that galaxies do not trace matter on very large scale.

It is generally accepted that the observed CMBR dipole arises from the motion of the local group relative to the CMBR frame. An alternative interpretation is that the dipole results from an ultra-large scale ($\lambda > 100 c/H_0)$ isocurvature perturbation. Recently it was argued that this alternative possibility is ruled out. We examine the growth of perturbations on scales larger than the Hubble radius and in view of this analysis, we show that the isocurvature interpretation is still a viable explanation. If the dipole is due to peculiar motion then it should appear in observations of other background sources provided that they are distant enough.

We compare the burst distribution of the new (2B) BATSE catalogue to a cosmological distribution. The observed distribution agrees well with a cosmological one, however, it is insensitive to cosmological parameters such as $\O$ and $\Lambda$. The bursts are not necessarily standard candles and their luminosity can vary by up to a factor of ten. The maximal red shift, $\zm$, of bursts longer than 2 sec is $2.1^{+1.}_{-.7}$ (assuming no evolution). The present data is insufficient to determine maximal red shift, $\zm$, of bursts shorter than 2 sec.

Neutron star binaries, such as the one observed in the famous binary pulsar PSR 1916+13, end their life in a catastrophic merge event (denoted here NS$^2$M). The merger releases $\approx 5 \cdot 10^{53}$ergs, mostly as neutrinos and gravitational radiation. A small fraction of this energy suffices to power $\gamma$-ray bursts (GRBs) at cosmological distances. Cosmological GRBs must pass, however, an optically thick fireball phase and the observed $\gamma$-rays emerge only at the end of this phase. Hence, it is difficult to determine the nature of the source from present observations (the agreement between the rates of GRBs and NS$^2$Ms being only an indirect evidence for this model). In the future a coinciding detection of a GRB and a gravitational radiation signal could confirm this model.

We compare the burst distribution of the new (2B) BATSE catalogue to a cosmological distribution. We find that the distribution is insensitive to cosmological parameters such as Omega and Lambda and to the width of the bursts luminosity function. The maximal red shift of the long bursts is ~2.1 (assuming no evolution) while Zm(long) of the short bursts is significantly lower Zm(short) =~ 0.5 In agreement with this relatively nearby origin of the short burst we find an indication that these bursts are correlated ( >=2 sigma level at 10 degrees) with Abell clusters. This is the first known correlation of the bursts with any other astrophysical population and if confirmed by further observations it will provides additional evidence for the cosmological origin of those bursts.

Fireballs are an essential part of any cosmological $\gamma$-ray burst. We drive here a stability criterion for fireballs and show that fireballs are Rayleigh-Taylor unstable in any region in which the entropy decreases outward. The instability begins to operate when the fireball becomes matter dominated. Among the possible implication of the instability are: (i) Conversion of a fraction of the radiation energy to a convective energy expressed in the motion of bubbles relative to each other. (ii) Penetration of fast bubbles through slower ones and creation of high $\gamma$ regimes which are essential for efficient conversion of the energy to $\gamma$-rays. (iii) Formation of rapid variation (of the scale of the bubbles) in the observed $\gamma$-rays.

$\gamma$-ray bursts have baffled theorists ever since their accidental discovery at the sixties. We suggest that these bursts originate in merger of neutron star binaries, taking place at cosmological distances. These mergers release $\approx 10^{54}ergs$, in what are possibly the strongest explosions in the Universe. If even a small fraction of this energy is channeled to an electromagnetic signal it will be detected as a grbs. We examine the virtues and limitations of this model and compare it with the recent Compton \g-ray observatory results. We also discuss the potential application of grbs to study cosmology and show that these burst might lead to a new and independent determination of both $H_0$ and $\Omega$.

Binary neutron stars merger (NS$^2$M) at cosmological distances is probably the only $\gamma$-ray bursts model based on an independently observed phenomenon which is known to be taking place at a comparable rate. We describe this model, its predictions and some open questions.

We present 3D Newtonian simulations of the coalescence of two neutron stars, using a Smoothed Particle Hydrodynamics (SPH) code. We begin the simulations with the two stars in a hard, circular binary, and have them spiral together as angular momentum is lost through gravitational radiation at the rate predicted by modeling the system as two point masses. We model the neutron stars as hard polytropes ($\gamma=2.4$) of equal mass, and investigate the effect of the initial spin of the two stars on the coalescence. The process of coalescence, from initial contact to the formation of an axially symmetric object, takes only a few orbital periods. Some of the material from the two neutron stars is shed, forming a thick disk around the central, coalesced object. The mass of this disk is dependent on the initial neutron star spins; higher spin rates resulting in greater mass loss, and thus more massive disks. For spin rates that are most likely to be applicable to real systems, the central coalesced object has a mass of $2.4M_\odot$, which is tantalizingly close to the maximum mass allowed by any neutron star equation of state for an object that is supported in part by rotation. Using a realistic nuclear equation of state, we estimate the temperature of the material after the coalescence. We find that the central object is at a temperature of $\sim 10$MeV, whilst the disk is heated by shocks to a temperature of 2-4MeV.

The sudden release of copious $\g$-ray photons into a compact region creates an opaque photon--lepton fireball due to the prolific production of electron--positron pairs. The photons that we observe in the bursts emerge only at the end of the fireball phase after it expanded sufficiently to become optically thin or after it converted its energy to the kinetic energy of relativistic baryons which convert it, in turn, to electromagnetic pulse via the interaction with interstellar matter. It is essential, therefore, to analyze the evolution of a fireball in order to comprehend the observed features of $\gamma$-ray bursts. We discuss various aspects of fireball hydrodynamics and the resulting emitted spectra.

We study the propagator of a non-relativistic, non-interacting particle in any non-relativistic ``time-machine'' spacetime of the type shown in Fig.~1: an external, flat spacetime in which two spatial regions, $V_-$ at time $t_-$ and $V_+$ at time $t_+$, are connected by two temporal wormholes, one leading from the past side of $V_-$ to t the future side of $V_+$ and the other from the past side of $V_+$ to the future side of $V_-$. We express the propagator explicitly in terms of those for ordinary, flat spacetime and for the two wormholes; and from that expression we show that the propagator satisfies completeness and unitarity in the initial and final ``chronal regions'' (regions without closed timelike curves) and its propagation from the initial region to the final region is unitary. However, within the time machine it satisfies neither completeness nor unitarity. We also give an alternative proof of initial-region-to-final-region unitarity based on a conserved current and Gauss's theorem. This proof can be carried over without change to most any non-relativistic time-machine spacetime; it is the non-relativistic version of a theorem by Friedman, Papastamatiou and Simon, which says that for a free scalar field, quantum mechanical unitarity follows from the fact that the classical evolution preserves the Klein-Gordon inner product.

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

We suggest to use the common appearance of voids with a scale of $5000 \kms$ in the galaxy distribution to estimate the power spectrum on this scale. We use a simple model for a gravitational formation of voids and we compare the results with the matter fluctuations as constrained by the CMBR observations of COBE. We find that a power spectrum $P(k)\propto k^n$ with $n \approx 1.25$ is compatible both with COBE and with gravitational growth of large voids in an $\Omega=1$ universe. A Harrison-Zel'dovich spectrum, $n=1$, normalized to produce the observed CMBR fluctuations, does not have enough power for gravitational growth of voids with a diameter of $5000 \kms$ in an $\Omega=1$ Universe. Such a spectrum would be compatible if (i) void diameters are smaller ($3500 \kms$), the voids are shallower or if the voids are rare (we assume that the universe is void filled); (ii) $\Omega < 1$ i.e. the Universe is open; (iii) galaxies do not trace matter on very large scale or if (iv) the voids do not grow gravitationally.

Coalescing binaries in distant galaxies are one of the most promising sources of gravitational waves detectable by the LIGO project.$^{[1-5]}$ They are also a copious source of neutrinos,$^{[1]}$ however these neutrino pulses are far too weak to be detected on earth. Several years ago Eichler \etal$\,$$^{[6]}$ suggested that they are also sources of $\gamma$-ray bursts (GRBs). Recently it was found$^{[7]}$ that GRBs are likely to be cosmological in origin, and coalescing binary systems$^{[6,8-11]}$ are probably the most promising cosmological sources. The current estimates of the burst and LIGO rates from a cosmologically distributed population are based on the systems observed in our galaxy.$^{[12-14]}$ These estimates are based on only three binarys so there are large statistical uncertainties in the coalescence rate. If we accept the cosmological/coalescing binary hypothesis for GRBs, then we get a more accurate estimate of the rate at which binaries coalesce and hence of the predicted LIGO signal. The association between GRBs and binaries can significantly improve the performance of LIGO. The detection of gravitational radiation from a GRB source not only confirms the coalescing binary model, but it will also provide information on the geometry and energy generation mechanism of the burst.

Kouveliotou et al. (1993) recently confirmed that gamma-ray bursts are bimodal in duration. In this paper we compute the statistical properties of the short ($\le 2$~s) and long ($>2$~s) bursts using a method of analysis that makes no assumption regarding the location of the bursts, whether in the Galaxy or at a cosmological distance. We find the 64 ms channel on BATSE to be more sensitive to short bursts and the 1024 ms channel is more sensitive to long bursts. We show that all the currently available data are consistent with the simple hypothesis that both short and long bursts have the same spatial distribution and that within each population the sources are standard candles. The rate of short bursts is $\sim 0.4$ of the rate of long bursts. Although the durations of short and long gamma-ray bursts span several orders of magnitude and the total energy of a typical short burst is smaller than that of a typical long burst by a factor of $\sim 20$, surprisingly the peak luminosities of the two kinds of bursts are equal to within a factor of $\sim 2$.

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

$\gamma$-ray bursts have baffled theorists ever since their accidental discovery at the sixties. We suggest that these bursts originate in merger of neutron star binaries, taking place at cosmological distances. These mergers release $\approx 10^{54}ergs$, in what are possibly the strongest explosions in the Universe. If even a small fraction of this energy is channeled to an electromagnetic signal it will be detected as a grbs. We examine the virtues and limitations of this model and compare it with the recent Compton \g-ray observatory results.

If gamma ray burst sources are in the galactic halo they inevitably involve a creation of an opaque pair plasma fireball, just like in cosmological sources. We find that the typical physical conditions in a galactic halo fireball are: optical depth $\ap 10^8$, thermal energy $\ap 100 KeV$, maximal relativstic expansion $\G \ap 300 $ and a maximal baryonic load of $\ap 10^{-15} M_\odot$. This does not rule out galctic halo models but it poses an addtional severe constraint on all such sources. A comparison of these conditions with the physical conditions at cosmological fireballs reveal that galactic halo fireballs are less favorable than cosmological ones as sources of gamma ray bursts.

It is proposed that gamma-ray bursts are created in the mergers of double neutron star binaries and black hole neutron star binaries at cosmological distances. Bursts with complex profiles and relatively long durations are the result of magnetic flares generated by the Parker instability in a post-merger differentially-rotating disk. Some bursts may also be produced through neutrino-antineutrino annihilation into electrons and positrons. In both cases, an optically thick fireball of size $\sles\ 100$ km is initially created, which expands ultrarelativistically to large radii before radiating. Several previous objections to the cosmological merger model are eliminated. It is predicted that $\gamma$-ray bursts will be accompanied by a burst of gravitational radiation from the spiraling-in binary which could be detected by LIGO.