Strong explosion of a compact star surrounded by a thick stellar wind drives a fast ($>0.1c$) radiation mediated shock (RMS) that propagates in the wind, and ultimately breaks out gradually once photons start escaping from the shock transition layer. In exceptionally strong or aspherical explosions the shock velocity may even be relativistic. The properties of the breakout signal depend on the dynamics and structure of the shock during the breakout phase. Here we present, for the first time, spectra and lightcurves of the breakout emission of fast Newtonian and mildly relativistic shocks, that were calculated using self-consistent Monte-Carlo simulations of finite RMS with radiative losses. We find a strong dependence of the $\nu F_\nu$ peak on shock velocity, ranging from $\sim 1$ keV for $v_s/c=0.1$ to $\sim 100$ keV for $v_s/c=0.5$, with a shift to lower energies as losses increase. For all cases studied the spectrum below the peak exhibits a nearly flat component ($F_\nu \sim \nu^0$) that extends down to the break frequency below which absorption becomes important. This implies much bright optical/UV emission than hitherto expected. The computed lightcurves show a gradual rise over tens to hundreds of seconds for representative conditions. The application to SN 2008D/XRT 080109 and the detectability limits are also discussed. We predict a detection rate of about one per year with eROSITA.

After its launch, a GRB jet propagates through a dense medium prior to its breakout. The interaction of the jet with the medium gives rise to the formation of a complex structured outflow, often referred to as a "structured jet". This structure is essential for our understanding of GRBs as it ultimately dictates their emission signatures. However, to date, the underlying physics which sets the post-breakout jet morphology remained unexplored, and its modeling has been mainly done by assuming ad-hoc functions. Here, using a set of 3D simulations we follow the evolution of hydrodynamic long and short GRB jets after breakout and provide a physically motivated post-breakout outflow structure of GRB jets (the structure of weakly magnetized jets is underway). Our simulations feature Rayleigh-Taylor fingers, which grow from the cocoon into the jet and destabilize it. The mixing of jet-cocoon material gives rise to a previously unidentified region between the two, which we denote the jet-cocoon interface (JCI). In lGRBs the mixing is strong, leading to most of the outflow energy to drift into the JCI. In sGRBs, where the medium is lighter, the mixing is weaker, and the JCI and the jet core hold a comparable amount of energy. Remarkably, the jet structure (core and JCI) in all systems can be characterized by simple angular power-law distributions of power and velocity, with power-law indices that depend solely on the level of mixing. This result supports the commonly used power-law angular distribution, and disfavors a Gaussian jet modeling. At larger angles, where the cocoon dominates, the structure is more complex including both an angular and a radial structure. The mixing shapes the prompt light curve and implies that typical afterglows of lGRBs are different from those of sGRBs. The predictions that we provide can be used to infer jet characteristics from prompt and afterglow observations.

We present the full panchromatic afterglow light curve data of GW170817, including new radio data as well as archival optical and X-ray data, between 0.5 and 940 days post-merger. By compiling all archival data, and reprocessing a subset of it, we have ensured that the panchromatic dataset is uniform and therefore immune to the differences in data processing or flux determination methods used by different groups. Simple power-law fits to the uniform afterglow light curve indicate a $t^{0.86\pm0.04}$ rise, a $t^{-1.90\pm0.12}$ decline, and a peak occurring at $155\pm4$ days. The afterglow is optically thin throughout its evolution, consistent with a single spectral index ($-0.569\pm0.002$) across all epochs. This gives a precise and updated estimate of the electron power-law index, $p=2.138\pm0.004$. By studying the diffuse X-ray emission from the host galaxy, we place a conservative upper limit on the hot ionized ISM density, $<$0.01 cm$^{-3}$, consistent with previous afterglow studies. Using the late-time afterglow data we rule out any long-lived neutron star remnant having magnetic field strength between 10$^{10.4}$ G and 10$^{16}$ G. Our fits to the afterglow data using an analytical model that includes VLBI proper motion from Mooley et al (2018), and a structured jet model that ignores the proper motion, indicates that the proper motion measurement needs to be considered while seeking an accurate estimate of the viewing angle.

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.

Strong variability is a common characteristic of the prompt emission of gamma-ray bursts (GRB). This observed variability is widely attributed to an intermittency of the central engine, through formation of strong internal shocks in the GRB-emitting jet expelled by the engine. In this paper we study numerically the propagation of hydrodynamic jets, injected periodically by a variable engine, through the envelope of a collapsed star. By post-processing the output of 3D numerical simulations, we compute the net radiative efficiency of the outflow. We find that all intermittent jets are subject to heavy baryon contamination that inhibits the emission at and above the photosphere well below detection limits. This is in contrast to continuous jets that, as shown recently, produce a highly variable gamma-ray photospheric emission with high efficiency, owing to the interaction of the jet with the stellar envelope. Our results challenge the variable engine model for hydrodynamic jets, and may impose constraints on the duty cycle of GRB engines. If such systems exist in nature, they are not expected to produce bright gamma-ray emission, but should appear as X-ray, optical and radio transients that resemble a delayed GRB afterglow signal.

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

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

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.

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

Mergers of binaries consisting of two neutron stars, or a black hole and a neutron star, offer a unique opportunity to study a range of physical and astrophysical processes using two different and almost orthogonal probes - gravitational waves (GW) and electromagnetic (EM) emission. The GW signal probes the binary and the physical processes that take place during the last stages of the merger, while the EM emission provides clues to the material that is thrown out following the merger. The accurate localization, which only the EM emission can provide, also indicates the astrophysical setting in which the merger took place. In addition, the combination of the two signals provides constraints on the nature of gravity and on the expansion rate of the Universe. The first detection of a binary neutron star merger by the LIGO-Virgo collaboration, GW170817, initiated the era of multi-messenger GW-EM astrophysics and demonstrated the great promise it holds. The event produced an unprecedented data set, and although it was only a single event, it provided remarkable results that revolutionized our knowledge of neutron star mergers. GW170817 is especially exciting since we know that it is not one of a kind and that many more events will be detected during the next decade. In this review, I summarize, first, the theory of EM emission from compact binary mergers, highlighting the unique information that the combined GW-EM detection provides. I then describe the entire set of GW and EM observations of GW170817, and summarize the range of insights that it offers. This includes clues about the role that similar events play in the r-process elements budget of the Universe, the neutron star equation of state, the properties of the relativistic outflow that followed the merger, and the connection between neutron star binary mergers and short gamma-ray bursts.

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

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.

The primary dissipation mechanism in jets of gamma-ray bursts (GRBs), and the high efficiency of the prompt emission are long standing issues. One possibility is strong collimation of a weakly magnetized relativistic jet by the surrounding medium, which can considerably enhance the efficiency of the photospheric emission. We derive a simple analytic criterion for the radiative efficiency of a collimated jet showing that it depends most strongly on the baryon loading. We confirm this analytic result by 3D numerical simulations, and further find that mixing of jet and cocoon material at the collimation throat leads to a substantial stratification of the outflow as well as sporadic loading, even if the injected jet is uniform and continuous. One consequence of this mixing is a strong angular dependence of the radiative efficiency. Another is large differences in the Lorentz factor of different fluid elements that lead to formation of internal shocks. Our analysis indicates that in both long and short GRBs a prominent photospheric component cannot be avoided when observed within an angle of a few degrees to the axis, unless the asymptotic Lorentz factor is limited by baryon loading at the jet base to $\Gamma_\infty <100$ (with a weak dependence on outflow power). Photon generation by newly created pairs behind the collimation shock regulates the observed temperature at $\sim 50~\theta_0^{-1}$ keV, where $\theta_0$ is the initial jet opening angle, in remarkable agreement with the observed peak energies of prompt emission spectra. Further consequences for the properties of the prompt emission are discussed at the end.

The light curves of type-II supernovae (SNe) are believed to be highly affected by recombination of hydrogen that takes place in their envelopes. In this work, we analytically investigate the transition from a fully ionized envelope to a partially recombined one and its effects on the SN light curve. The motivation is to establish the underlying processes that dominate the evolution at late times when recombination takes place in the envelope, yet early enough so that $^{56}$Ni decay is a negligible source of energy. We consider the diffusion of photons through the envelope while analyzing the ionization fraction and the coupling between radiation and gas, and find that the main effect of recombination is on the evolution of the observed temperature. Before recombination the temperature decreases relatively fast, while after recombination starts it significantly reduces the rate at which the observed temperature drops with time. This behaviour is the main cause for the observed flattening in the optical bands, where for a typical red supergiant explosion, the recombination wave affects the bolometric luminosity only mildly during most of the photospheric phase. Moreover, the plateau phase observed in some type-II SNe is not a generic result of recombination, and it also depends on the density structure of the progenitor. This is one possible explanation to the different light curve decay rates observed in type II (P and L) SNe.

Gamma-ray astronomy explores the most energetic photons in nature to address some of the most pressing puzzles in contemporary astrophysics. It encompasses a wide range of objects and phenomena: stars, supernovae, novae, neutron stars, stellar-mass black holes, nucleosynthesis, the interstellar medium, cosmic rays and relativistic-particle acceleration, and the evolution of galaxies. MeV gamma-rays provide a unique probe of nuclear processes in astronomy, directly measuring radioactive decay, nuclear de-excitation, and positron annihilation. The substantial information carried by gamma-ray photons allows us to see deeper into these objects, the bulk of the power is often emitted at gamma-ray energies, and radioactivity provides a natural physical clock that adds unique information. New science will be driven by time-domain population studies at gamma-ray energies. This science is enabled by next-generation gamma-ray instruments with one to two orders of magnitude better sensitivity, larger sky coverage, and faster cadence than all previous gamma-ray instruments. This transformative capability permits: (a) the accurate identification of the gamma-ray emitting objects and correlations with observations taken at other wavelengths and with other messengers; (b) construction of new gamma-ray maps of the Milky Way and other nearby galaxies where extended regions are distinguished from point sources; and (c) considerable serendipitous science of scarce events -- nearby neutron star mergers, for example. Advances in technology push the performance of new gamma-ray instruments to address a wide set of astrophysical questions.

The breakout of a fast ($>0.1 c$), yet sub-relativistic shock from a thick stellar wind is expected to produce a pulse of X-rays with a rise time of seconds to hours. Here, we construct a semi-analytic model for the breakout of a sub-relativistic, radiation-mediated shock from a thick stellar wind, and use it to compute the spectrum of the breakout emission. The model incorporates photon escape through the finite optical depth wind, assuming a diffusion approximation and a quasi-steady evolution of the shock structure during the breakout phase. We find that in sufficiently fast shocks, for which the breakout velocity exceeds about $0.1c$, the time-integrated spectrum of the breakout pulse is non-thermal, and the time-resolved temperature is expected to exhibit substantial decrease (roughly by one order of magnitude) during breakout, when the flux is still rising, because of the photon generation by the shock compression associated with the photon escape. We also derive a closure relation between the breakout duration, peak luminosity, and characteristic temperature that can be used to test whether an observed X-ray flare is consistent with being associated with a sub-relativistic shock breakout from a thick stellar wind or not. We also discuss implications of the spectral softening for a possible breakout event XRT 080109/SN 2008D.

Motivated by GW170817 we examine constraints that observations put on the angular structure of long gamma-ray burst (GRB) jets. First, the relatively narrow observed distribution of $E_{X}/E_{\gamma}$ (the isotropic equivalent early X-ray afterglow to prompt $\gamma$-ray energy ratio) implies that at any angle that $\gamma$-rays are emitted the Lorentz factor must be high. Specifically, the Lorentz factor of $\gamma$-ray emitting material cannot drop rapidly with angle, and must be $\Gamma(\theta)\gtrsim 50$ even if there are angles for which the gamma-ray received energy is lower by three orders of magnitude compared to the jet core. Second, jets with an angular structure of the $\gamma$-ray emission that over-produce events with a $\gamma$-ray luminosity below the peak of the observed luminosity function are ruled-out. This eliminates models in which the $\gamma$-ray energy angular distribution isn't sufficiently steep and the Lorentz factor distribution isn't sufficiently shallow. Finally, models with a steep structure (e.g. Gaussian) which are detected away from the jet core generate afterglow light-curves that were never observed. We conclude that even if the jet kinetic energy distribution drops continuously with latitude, efficient $\gamma$-ray emission seems to be restricted to material with $\Gamma\gtrsim 50$ and is most likely confined to a narrow region around the core. While our study is confined to long GRBs, where the observed sample is larger and more complete, there are indications that similar conclusions may be applicable also to short GRBs. We discuss the possible implications to the $\gamma$-rays observed in GRB 170817A.

In the present study, we systematically explore the effect of the radioactive 56Ni and its mixing properties in the ejecta on the plateau of Type IIP supernovae (SNe). We evaluate the importance of 56Ni in shaping light curves of SNe IIP by simulating light curves for two red supergiant models using different amounts of 56Ni and with different types of mixing: uniform distribution of 56Ni out to different fractions of the envelope and "boxcar" distribution of 56Ni. We find, similarly to previous studies, that 56Ni extends duration of the plateau. We find a formula to estimate the extension based on the observed bolometric light curves and show that for most SNe IIP 56Ni extends the plateau by about 20%. Another effect of 56Ni consists in reduction of the plateau decline rate, i.e. 56Ni presented in the ejecta flattens the plateau. Our simulations suggest that for typical SNe IIP it can reduce the decline rate by about 1 mag/100 day. We find that for the contribution of 56Ni seen in most SNe our simulated bolometric light curves resemble observed ones for various types of 56Ni mixing. We thereby cannot determine the level of 56Ni mixing in these SNe based on the light curve alone. However, for SN2009ib we find that only a model where 56Ni is mixed significantly throughout most of the hydrogen envelope is consistent with the observed light curve. Our light curves are available via link https://wwwmpa.mpa-garching.mpg.de/ccsnarchive/data/Kozyreva2018/.

Several processes may limit the accretion rate onto a super-massive black hole (SMBH). Two processes that are commonly considered (e.g., for sub-grid prescriptions) are Bondi-Hoyle-Lyttleton accretion and the Eddington limit. A third one is AGN wind feedback. It has been long suggested that such a wind feedback regulates the final SMBH mass, however, it has been shown recently that AGN winds can also regulate the average accretion rate at a level consistent with observations of high redshift AGNs. In this paper we study the effect of wind feedback on the accretion rate using 2D, high resolution hydrodynamic simulations, that incorporate a self-consistent wind injection scheme and resolves the SMBH sphere of influence. Two different cases are explored and compared: one in which the initial gas density is uniform, and one in which it has an isothermal sphere profile. We also compare simulations with and without cooling. Our main finding is that for reasonable parameters, AGN feedback always limits the accretion rate to be far below the Bondi-Hoyle-Lyttleton limit. For typical wind parameters and a uniform ISM densities of the order $1 {cm}^{-3}$, the accretion rate is found to be several orders of magnitudes smaller than that inferred in large samples of high redshift AGNs. On the other hand, the accretion rate obtained for initially isothermal density profile is found to be consistent with the observations, particularly when cooling is included. Furthermore, it roughly scales as $\sigma^{5}$ with the velocity dispersion of the bulge, in accord with the $M-\sigma$ relation.

The binary neutron star merger GW170817 was accompanied by radiation across the electromagnetic spectrum and localized to the galaxy NGC 4993 at a distance of 41+/-3 Mpc. The radio and X-ray afterglows of GW170817 exhibited delayed onset, a gradual rise in the emission with time as t^0.8, a peak at about 150 days post-merger, followed by a relatively rapid decline. To date, various models have been proposed to explain the afterglow emission, including a choked-jet cocoon and a successful-jet cocoon (a.k.a. structured jet). However, the observational data have remained inconclusive as to whether GW170817 launched a successful relativistic jet. Here we show, through Very Long Baseline Interferometry, that the compact radio source associated with GW170817 exhibits superluminal motion between two epochs at 75 and 230 days post-merger. This measurement breaks the degeneracy between the models and indicates that, while the early-time radio emission was powered by a wider-angle outflow (cocoon), the late-time emission was most likely dominated by an energetic and narrowly-collimated jet, with an opening angle of <5 degrees, and observed from a viewing angle of about 20 degrees. The imaging of a collimated relativistic outflow emerging from GW170817 adds substantial weight to the growing evidence linking binary neutron star mergers and short gamma-ray bursts.

We present new 0.6-10 GHz observations of the binary neutron star merger GW170817 covering the period up to 300 days post-merger, taken with the Karl G. Jansky Very Large Array, the Australia Telescope Compact Array, the Giant Metrewave Radio Telescope and the MeerKAT telescope. We use these data to precisely characterize the decay phase of the late-time radio light curve. We find that the temporal decay is consistent with a power-law slope of t^-2.2, and that the transition between the power-law rise and decay is relatively sharp. Such a slope cannot be produced by a quasi-isotropic (cocoon-dominated) outflow, but is instead the classic signature of a relativistic jet. This provides strong observational evidence that GW170817 produced a successful jet, and directly demonstrates the link between binary neutron star mergers and short-hard GRBs. Using simple analytical arguments, we derive constraints on the geometry and the jet opening angle of GW170817. These results are consistent with those from our companion Very Long Baseline Interferometry (VLBI) paper, reporting superluminal motion in GW170817.

We present visible-light and ultraviolet (UV) observations of the supernova PTF12glz. The SN was discovered and monitored in near-UV and R bands as part of a joint GALEX and Palomar Transient Factory campaign. It is among the most energetic Type IIn supernovae observed to date (~10^51erg). If the radiated energy mainly came from the thermalization of the shock kinetic energy, we show that PTF12glz was surrounded by ~1 solar mass of circumstellar material (CSM) prior to its explosive death. PTF12glz shows a puzzling peculiarity: at early times, while the freely expanding ejecta are presumably masked by the optically thick CSM, the radius of the blackbody that best fits the observations grows at ~8000km/s. Such a velocity is characteristic of fast moving ejecta rather than optically thick CSM. This phase of radial expansion takes place before any spectroscopic signature of expanding ejecta appears in the spectrum and while both the spectroscopic data and the bolometric luminosity seem to indicate that the CSM is optically thick. We propose a geometrical solution to this puzzle, involving an aspherical structure of the CSM around PTF12glz. By modeling radiative diffusion through a slab of CSM, we show that an aspherical geometry of the CSM can result in a growing effective radius. This simple model also allows us to recover the decreasing blackbody temperature of PTF12glz. SLAB-Diffusion, the code we wrote to model the radiative diffusion of photons through a slab of CSM and evaluate the observed radius and temperature, is made available on-line.

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.

The Hubble constant ($H_0$) measures the current expansion rate of the Universe, and plays a fundamental role in cosmology. Tremendous effort has been dedicated over the past decades to measure $H_0$. Notably, Planck cosmic microwave background (CMB) and the local Cepheid-supernovae distance ladder measurements determine $H_0$ with a precision of $\sim 1\%$ and $\sim 2\%$ respectively. A $3$-$\sigma$ level of discrepancy exists between the two measurements, for reasons that have yet to be understood. Gravitational wave (GW) sources accompanied by electromagnetic (EM) counterparts offer a completely independent standard siren (the GW analogue of an astronomical standard candle) measurement of $H_0$, as demonstrated following the discovery of the neutron star merger, GW170817. This measurement does not assume a cosmological model and is independent of a cosmic distance ladder. The first joint analysis of the GW signal from GW170817 and its EM localization led to a measurement of $H_0=74^{+16}_{-8}$ km/s/Mpc (median and symmetric $68\%$ credible interval). In this analysis, the degeneracy in the GW signal between the source distance and the weakly constrained viewing angle dominated the $H_0$ measurement uncertainty. Recently, Mooley et al. (2018) obtained tight constraints on the viewing angle using high angular resolution imaging of the radio counterpart of GW170817. Here we obtain a significantly improved measurement $H_0=68.9^{+4.7}_{-4.6}$ km/s/Mpc by using these new radio observations, combined with the previous GW and EM data. We estimate that 15 more localized GW170817-like events (comparable signal-to-noise ratio, favorable orientation), having radio images and light curve data, will potentially bring resolution to the tension between the Planck and Cepheid-supernova measurements, as compared to 50-100 GW events without such data.

The gravitational wave event GW170817 from a binary neutron star (NS) merger is accompanied by electromagnetic counterparts, and the optical and near-infared emission is called a macronova (or kilonova). Although the radioactivity of synthesized r-process elements is widely discussed as an energy source, its decisive evidence is not clearly shown yet. We discuss a macronova powered by the central engine activities such as jet activities and X-rays from the matter fallback, and show that the engine model allows much broader parameter spaces, in particular smaller ejecta mass ($\sim10^{-4}-0.01\,\Msun$) than the r-process model. The blue and red macronovae are naturally explained by various combinations of the ejecta such as a cocoon and merger ejecta with the energy sources of jets and X-rays. The required energy injection is very similar to the X-ray excess observed in GRB 130603B with the power-law slope of $\sim-1.3$. The required lanthanoid fraction for the opacity can be also consistent with the Galactic one. Early or late multi-wavelength observations are crucial for revealing the central engine of short gamma-ray bursts and the r-process nucleosynthesis.

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.

Relativistic jets reside in high-energy astrophysical systems of all scales. Their interaction with the surrounding media is critical as it determines the jet evolution, observable signature, and feedback on the environment. During its motion the interaction of the jet with the ambient media inflates a highly pressurized cocoon, which under certain conditions collimates the jet and strongly affects its propagation. Recently, Bromberg et al. (2011) derived a general simplified (semi)analytic solution for the evolution of the jet and the cocoon in case of an unmagnetized jet that propagates in a medium with a range of density profiles. In this work we use a large suite of 2D and 3D relativistic hydrodynamic simulations in order to test the validity and accuracy of this model. We discuss the similarities and differences between the analytic model and numerical simulations and also, to some extent, between 2D and 3D simulations. Our main finding is that although the analytic model is highly simplified, it properly predicts the evolution of the main ingredients of the jet-cocoon system, including its temporal evolution and the transition between various regimes (e.g., collimated to uncollimated). The analytic solution predicts a jet head velocity that is faster by a factor of about 3 compared to the simulations, as long as the head velocity is Newtonian. We use the results of the simulations to calibrate the analytic model which significantly increases its accuracy. We provide an applet that calculates semi-analytically the propagation of a jet in an arbitrary density profile defined by the user at http://www.astro.tau.ac.il/~ore/propagation.html.

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.

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.

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 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.

Every supernova hitherto observed has been considered to be the terminal explosion of a star. Moreover, all supernovae with absorption lines in their spectra show those lines decreasing in velocity over time, as the ejecta expand and thin, revealing slower moving material that was previously hidden. In addition, every supernova that exhibits the absorption lines of hydrogen has one main light-curve peak, or a plateau in luminosity, lasting approximately 100 days before declining. Here we report observations of iPTF14hls, an event that has spectra identical to a hydrogen-rich core-collapse supernova, but characteristics that differ extensively from those of known supernovae. The light curve has at least five peaks and remains bright for more than 600 days; the absorption lines show little to no decrease in velocity; and the radius of the line-forming region is more than an order of magnitude bigger than the radius of the photosphere derived from the continuum emission. These characteristics are consistent with a shell of several tens of solar masses ejected by the star at supernova-level energies a few hundred days before a terminal explosion. Another possible eruption was recorded at the same position in 1954. Multiple energetic pre-supernova eruptions are expected to occur in stars of 95-130 solar masses, which experience the pulsational pair instability. That model, however, does not account for the continued presence of hydrogen, or the energetics observed here. Another mechanism for the violent ejection of mass in massive stars may be required.

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.

With the first direct detection of merging black holes in 2015, the era of gravitational wave (GW) astrophysics began. A complete picture of compact object mergers, however, requires the detection of an electromagnetic (EM) counterpart. We report ultraviolet (UV) and X-ray observations by Swift and the Nuclear Spectroscopic Telescope ARray (NuSTAR) of the EM counterpart of the binary neutron star merger GW170817. The bright, rapidly fading ultraviolet emission indicates a high mass ($\approx0.03$ solar masses) wind-driven outflow with moderate electron fraction ($Y_{e}\approx0.27$). Combined with the X-ray limits, we favor an observer viewing angle of $\approx 30^{\circ}$ away from the orbital rotation axis, which avoids both obscuration from the heaviest elements in the orbital plane and a direct view of any ultra-relativistic, highly collimated ejecta (a gamma-ray burst afterglow).

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.

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.

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.

The effect of AGN wind feedback on the accretion rate and mass evolution of supermassive black holes (SMBH) is considered. It is shown, under reasonable assumptions, that the rate at which gas can be supplied to a SMBH at the center of a young galaxy, is limited to $\sim 20 (\sigma/350~ {\rm km ~s^{-1}})^4~M_\odot/yr$ (where $\sigma$ is the velocity dispersion of the host bulge) by interaction of winds expelled from the innermost regions of the accretion flow with the gas in the bulge. This rate is independent of the black hole mass but is sensitive to the properties of the host bulge. It is further argued that the interaction of the wind and the inflowing gas in the bulge can strongly affect the structure of the accretion flow in the super-Eddington regime, potentially leading to highly super-Eddington accretion into the SMBH. About 300 Myr after the onset of the accretion phase, the AGN wind expels all the gas from the bulge and the accretion rate is strongly suppressed. This scenario is in remarkable agreement with recent observations, as it (i) reveals the origin of the maximal observed accretion rates, (ii) accounts for the inferred growth rate of SMBHs at high redshift (independent of their initial seeds masses), (iii) elucidates the decline in the accretion rate at lower redshifts, and (iv) explains the relation between $\sigma$ and the mass $M_{BH}$ of the central SMBH, measured in the local Universe.

We construct an analytic model for the breakout of a relativistic radiation mediated shock from a stellar wind, and exploit it to calculate the observational diagnostics of the breakout signal. The model accounts for photon escape through the finite optical depth wind, and treats the fraction of downstream photons escaping to infinity as an adiabatic parameter that evolves in a quasi-steady manner. It is shown that the shock is mediated by radiation even when a large fraction of the downstream photons escape, owing to self-generation and adjustment of opacity through accelerated pair creation. Relativistic breakout occurs at radii at which the total optical depth of the wind ahead of the shock is $\sim (m_e/m_p)\Gamma_{sh}$, provided that the local shock Lorentz factor $\Gamma_{sh}$ exceeds unity at this location. Otherwise the breakout occurs in the Newtonian regime. A relativistic breakout is expected in a highly energetic spherical explosion ($10^{52}-10^{53}$ erg) of a Wolf-Rayet star, or in cases where a smaller amount of energy ($\sim 10^{51}$ erg) is deposited by a jet in the outer layers of the star. The properties of the emission observed in such explosions during the relativistic breakout are derived. We find that for typical parameters about $10^{48}$ ergs are radiated in the form of MeV gamma-rays over a duration that can range from a fraction of a second to an hour. Such a signal may be detectable out to 10-100 Mpc by current gamma-ray satellites.

In this work we present a uniform analysis of the temperature evolution and bolometric luminosity of a sample of 29 type-II supernovae (SNe), by fitting a black body model to their multi-band photometry. Our sample includes only SNe with high quality multi-band data and relatively well sampled time coverage. Most of the SNe in our sample were detected less than a week after explosion so their light curves cover the evolution both before and after recombination starts playing a role. We use this sample to study the signature of hydrogen recombination, which is expected to appear once the observed temperature drops to $\approx 7,000$K. Theory predicts that before recombination starts affecting the light curve, both the luminosity and the temperature should drop relatively fast, following a power-law in time. Once the recombination front reaches inner parts of the outflow, it sets the observed temperature to be nearly constant, and slows the decline of the luminosity (or even leads to a re-brightening). We compare our data to analytic studies and find strong evidence for the signature of recombination. We also find that the onset of the optical plateau in a given filter, is effectively the time at which the black body peak reaches the central wavelength of the filter, as it cools, and it does not correspond to the time at which recombination starts affecting the emission.

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

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

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.

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.

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.

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

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.

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

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.

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

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

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.

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.

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

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

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

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

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

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

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

(Abridged) Some supernovae (SNe) may be powered by the interaction of the SN ejecta with a large amount of circumstellar matter (CSM). Here we outline a method to measure the mass of the optically thick CSM around such SNe. We present observations of SN2010jl, including the first detection of a SN using NuSTAR. The total radiated luminosity of SN2010jl is extreme, at least 9e50 erg. By modeling the visible-light data, we robustly show that the mass of the circumstellar material within ~1e16 cm of the progenitor was in excess of 10 solar masses, likely ejected tens of years prior to the SN explosion. Our modeling suggests that the shock velocity during shock breakout was ~6000 km/s, decelerating to ~2600 km/s about two years after maximum light. Our late-time NuSTAR+XMM spectra of the SN presumably provide the first direct measurement of SN shock velocity two years after the SN maximum light -- measured to be in the range of 2000 to 4500 km/s if the ions and electrons are in equilibrium, and >~2000 km/s if they are not in equilibrium. This measurement is in agreement with the shock velocity predicted by our modeling of the optical data. We also show that the mean radial density distribution of the CSM roughly follows an r^-2 law. A possible explanation for the massive CSM with a wind-like profile is that they are the result of multiple pulsational pair instability events prior to the SN explosion, separated from each other by years.

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.

The time-variable electromagnetic sky has been well-explored at a wide range of wavelengths. Numerous high-energy space missions take advantage of the dark Gamma-ray and X-ray sky and utilize very wide field detectors to provide almost continuous monitoring of the entire celestial sphere. In visible light, new wide-field ground-based surveys cover wide patches of sky with ever decreasing cadence, progressing from monthly-weekly time scale surveys to sub-night sampling. In the radio, new powerful instrumentation offers unprecedented sensitivity over wide fields of view, with pathfinder experiments for even more ambitious programs underway. In contrast, the ultra-violet (UV) variable sky is relatively poorly explored, even though it offers exciting scientific prospects. Here, we review the potential scientific impact of a wide-field UV survey on the study of explosive and other transient events, as well as known classes of variable objects, such as active galactic nuclei and variable stars. We quantify our predictions using a fiducial set of observational parameters which are similar to those envisaged for the proposed ULTRASAT mission. We show that such a mission would be able to revolutionize our knowledge about massive star explosions by measuring the early UV emission from hundreds of events, revealing key physical parameters of the exploding progenitor stars. Such a mission would also detect the UV emission from many tens of tidal-disruption events of stars by super massive black holes at galactic nuclei and enable a measurement of the rate of such events. The overlap of such a wide-field UV mission with existing and planned gravitational-wave and high-energy neutrino telescopes makes it especially timely.

Various lines of evidence suggest that very massive stars experience extreme mass-loss episodes shortly before they explode as a supernova. Interestingly, several models predict such pre-explosion outbursts. Establishing a causal connection between these mass-loss episodes and the final supernova explosion will provide a novel way to study pre-supernova massive-star evolution. Here we report on observations of a remarkable mass-loss event detected 40 days prior to the explosion of the Type IIn supernova SN 2010mc (PTF 10tel). Our photometric and spectroscopic data suggest that this event is a result of an energetic outburst, radiating at least 6x10^47 erg of energy, and releasing about 0.01 Solar mass at typical velocities of 2000 km/s. We show that the temporal proximity of the mass-loss outburst and the supernova explosion implies a causal connection between them. Moreover, we find that the outburst luminosity and velocity are consistent with the predictions of the wave-driven pulsation model and disfavor alternative suggestions.

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.

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

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

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.

We present X-ray, optical, near-infrared, and radio observations of GRBs 110709B and 111215A, as well as optical and near-IR observations of their host galaxies. The combination of X-ray detections and deep optical/near-infrared limits establish both bursts as "dark". Sub-arcsecond positions enabled by radio detections lead to robust host galaxy associations, with optical detections that indicate z < 4 (110709B) and 1.8 < z < 2.7 (111215A). Using the radio and X-ray data for each burst we find that GRB 110709B requires A_V > 5.3 mag and GRB 111215A requires A_V > 8.5 mag (z=2), among the largest extinction values inferred for dark bursts to-date. The two bursts also exhibit large neutral hydrogen column densities (N_H > 10^22/cm^2; z=2) as inferred from their X-ray spectra, in agreement with the trend for dark GRBs. Finally, we find that for both bursts the afterglow emission is best explained by a collimated outflow with a total beaming-corrected energy of E_gamma+E_K ~ 7-9 x 10^51 erg (z=2) expanding into a wind medium with a high density (n~100-350 cm^-3 at 10^17 cm). While the energy release is typical of long GRBs, the inferred density may be indicative of larger mass loss rates for GRB progenitors in dusty (and hence metal rich) environments. This study establishes the critical role of radio observations in determining the origin and properties of dark GRBs. Observations with the JVLA and ALMA will provide a sample with sub-arcsecond positions and robust host associations that will help shed light on obscured star formation and the role of metallicity in GRB progenitors.

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.

We report the discovery and detailed monitoring of X-ray emission associated with the Type IIb SN 2011dh using data from the Swift and Chandra satellites, placing it among the best studied X-ray supernovae to date. We further present millimeter and radio data obtained with the SMA, CARMA, and EVLA during the first three weeks after explosion. Combining these observations with early optical photometry, we show that the panchromatic dataset is well-described by non-thermal synchrotron emission (radio/mm) with inverse Compton scattering (X-ray) of a thermal population of optical photons. In this scenario, the shock partition fractions deviate from equipartition by a factor, (e_e/e_B) ~ 30. We derive the properties of the shockwave and the circumstellar environment and find a shock velocity, v~0.1c, and a progenitor mass loss rate of ~6e-5 M_sun/yr. These properties are consistent with the sub-class of Type IIb SNe characterized by compact progenitors (Type cIIb) and dissimilar from those with extended progenitors (Type eIIb). Furthermore, we consider the early optical emission in the context of a cooling envelope model to estimate a progenitor radius of ~1e+11 cm, in line with the expectations for a Type cIIb SN. Together, these diagnostics are difficult to reconcile with the extended radius of the putative yellow supergiant progenitor star identified in archival HST observations, unless the stellar density profile is unusual. Finally, we searched for the high energy shock breakout pulse using X-ray and gamma-ray observations obtained during the purported explosion date range. Based on the compact radius of the progenitor, we estimate that the breakout pulse was detectable with current instruments but likely missed due to their limited temporal/spatial coverage. [Abridged]

We present the discovery of PTF 10vgv, a Type Ic supernova detected by the Palomar Transient Factory, using the Palomar 48-inch telescope (P48). R-band observations of the PTF 10vgv field with P48 probe the supernova emission from its very early phases (about two weeks before R-band maximum), and set limits on its flux in the week prior to the discovery. Our sensitive upper limits and early detections constrain the post-shock-breakout luminosity of this event. Via comparison to numerical (analytical) models, we derive an upper-limit of R \lesssim 4.5 Rsun (R \lesssim 1 Rsun) on the radius of the progenitor star, a direct indication in favor of a compact Wolf-Rayet star. Applying a similar analysis to the historical observations of SN 1994I, yields R \lesssim 1/4 Rsun for the progenitor radius of this supernova.

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

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 report on a re-analysis of archival data from the Very Large Array for a sample of ten long duration radio transients reported by Bower and others. These transients have an implied all-sky rate that would make them the most common radio transient in the sky and yet most have no quiescent counterparts at other wavelengths and therefore no known progenitor (other than Galactic neutron stars). We find that more than half of these transients are due to rare data artifacts. The remaining sources have lower signal-to-noise ratio (SNR) than initially reported by 1 to 1.5-sigma. This lowering of SNR matters greatly since the sources are at the threshold. We are unable to decisively account for the differences. By two orthogonal criteria one source appears to be a good detection. Thus the rate of long duration radio transients without optical counterparts is, at best, comparable to that of the class of recently discovered Swift J1644+57 nuclear radio transients. We revisit the known and expected classes of long duration radio transients and conclude that the dynamic radio sky remains a rich area for further exploration. Informed by the experience of past searches for radio transients, we suggest that future surveys pay closer attention to rare data errors and ensure that a wealth of sensitive multi-wavelength data be available in advance of the radio observations and that the radio searches should have assured follow-up resources.

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.

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.

Active galactic nuclei (AGN), powered by long-term accretion onto central supermassive black holes, produce relativistic jets with lifetimes of greater than one million yr that preclude observations at birth. Transient accretion onto a supermassive black hole, for example through the tidal disruption of a stray star, may therefore offer a unique opportunity to observe and study the birth of a relativistic jet. On 2011 March 25, the Swift {\gamma}-ray satellite discovered an unusual transient source (Swift J164449.3+573451) potentially representing such an event. Here we present the discovery of a luminous radio transient associated with Swift J164449.3+573451, and an extensive set of observations spanning centimeter to millimeter wavelengths and covering the first month of evolution. These observations lead to a positional coincidence with the nucleus of an inactive galaxy, and provide direct evidence for a newly-formed relativistic outflow, launched by transient accretion onto a million solar mass black hole. While a relativistic outflow was not predicted in this scenario, we show that the tidal disruption of a star naturally explains the high-energy properties, radio luminosity, and the inferred rate of such events. The weaker beaming in the radio compared to {\gamma}-rays/X-rays, suggests that radio searches may uncover similar events out to redshifts of z ~ 6.

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

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.

Supernovae (SNe) are stellar explosions driven by gravitational or thermonuclear energy, observed as electromagnetic radiation emitted over weeks or more. In all known SNe, this radiation comes from internal energy deposited in the outflowing ejecta by either radioactive decay of freshly-synthesized elements (typically 56Ni), stored heat deposited by the explosion shock in the envelope of a supergiant star, or interaction between the SN debris and slowly-moving, hydrogen-rich circumstellar material. Here we report on a new class of luminous SNe whose observed properties cannot be explained by any of these known processes. These include four new SNe we have discovered, and two previously unexplained events (SN 2005ap; SCP 06F6) that we can now identify as members. These SNe are all ~10 times brighter than SNe Ia, do not show any trace of hydrogen, emit significant ultra-violet (UV) flux for extended periods of time, and have late-time decay rates which are inconsistent with radioactivity. Our data require that the observed radiation is emitted by hydrogen-free material distributed over a large radius (~10^15 cm) and expanding at high velocities (>10^4 km s^-1). These long-lived, UV-luminous events can be observed out to redshifts z>4 and offer an excellent opportunity to study star formation in, and the interstellar medium of, primitive distant galaxies.

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

We use a large sample of GRB afterglow and prompt-emission data (adding further GRB afterglow observations in this work) to compare the optical afterglows (or the lack thereof) of Type I GRBs with those of Type II GRBs. In comparison to the afterglows of Type II GRBs, we find that those of Type I GRBs have a lower average luminosity and show an intrinsic spread of luminosities at least as wide. From late and deep upper limits on the optical transients, we establish limits on the maximum optical luminosity of any associated supernova, confirming older works and adding new results. We use deep upper limits on Type I GRB optical afterglows to constrain the parameter space of possible mini-SN emission associated with a compact-object merger. Using the prompt emission data, we search for correlations between the parameters of the prompt emission and the late optical afterglow luminosities. We find tentative correlations between the bolometric isotropic energy release and the optical afterglow luminosity at a fixed time after trigger (positive), and between the host offset and the luminosity (negative), but no significant correlation between the isotropic energy release and the duration of the GRBs. We also discuss three anomalous GRBs, GRB 060505, GRB 060614, and GRB 060121, in the light of their optical afterglow luminosities. (Abridged)

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

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.

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

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

We have gathered optical photometry data from the literature on a large sample of Swift-era gamma-ray burst (GRB) afterglows including GRBs up to September 2009, for a total of 76 GRBs, and present an additional three pre-Swift GRBs not included in an earlier sample. Furthermore, we publish 840 additional new photometry data points on a total of 42 GRB afterglows, including large data sets for GRBs 050319, 050408, 050802, 050820A, 050922C, 060418, 080413A and 080810. We analyzed the light curves of all GRBs in the sample and derived spectral energy distributions for the sample with the best data quality, allowing us to estimate the host galaxy extinction. We transformed the afterglow light curves into an extinction-corrected z=1 system and compared their luminosities with a sample of pre-Swift afterglows. The results of a former study, which showed that GRB afterglows clustered and exhibited a bimodal distribution in luminosity space, is weakened by the larger sample. We found that the luminosity distribution of the two afterglow samples (Swift-era and pre-Swift) are very similar, and that a subsample for which we were not able to estimate the extinction, which is fainter than the main sample, can be explained by assuming a moderate amount of line-of-sight host extinction. We derived bolometric isotropic energies for all GRBs in our sample, and found only a tentative correlation between the prompt energy release and the optical afterglow luminosity at one day after the GRB in the z=1 system. A comparative study of the optical luminosities of GRB afterglows with echelle spectra (which show a high number of foreground absorbing systems) and those without reveals no indication that the former are statistically significantly more luminous. (abridged)

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

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.

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

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

We present analytic approximations to the optically thin synchrotron and synchrotron self-Compton (SSC) spectra when Klein-Nishina (KN) effects are important and pair production and external radiation fields can be neglected. This theory is useful for analytical treatment of radiation from astrophysical sources, such as gamma-ray bursts (GRBs), active galactic nuclei and pulsar wind nebula, where KN effects may be important. We consider a source with a continuous injection of relativistic electrons with a power-law energy distribution above some typical injection energy. We find that the synchrotron-SSC spectra can be described by a broken power-law, and provide analytic estimates for the break frequencies and power-law indices. In general, we show that the dependence of the KN cross-section on the energy of the upscattering electron results in a hardening of the energy distribution of fast cooling electrons and therefore in a hardening of the observed synchrotron spectrum. As a result the synchrotron spectrum of fast cooling electrons, below the typical injection energy, can be as hard as $F_\nu \propto \nu^0$, instead of the classical $\nu^{-1/2}$ when KN effects are neglected. The synchrotron energy output can be dominated by electrons with energy above the typical injection energy. We solve self-consistently for the cooling frequency and find that the transition between synchrotron and SSC cooling can result in a discontinuous variations of the cooling frequency and the synchrotron and SSC spectra. We demonstrate the application of our results to theory by applying them to prompt and afterglow emission models of GRBs.

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.

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.

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?

We investigate the GeV emission from gamma-ray bursts (GRBs), using the results from the Energetic Gamma Ray Experimental Telescope (EGRET), and in view of the Gamma-ray Large Area Space Telescope (GLAST). Assuming that the conventional prompt and afterglow photons originate from synchrotron radiation, we compare an accompanying inverse-Compton component with EGRET measurements and upper limits on GeV fluence, taking Klein-Nishina feedback into account. We find that EGRET constraints are consistent with the theoretical framework of the synchrotron self-Compton model for both prompt and afterglow phases, and discuss constraints on microphysical parameters in both phases. Based on the inverse-Compton model and using EGRET results, we predict that GLAST would detect GRBs with GeV photons at a rate >~20 yr^{-1} from each of the prompt and afterglow phases. This rate applies to the high-energy tail of the prompt synchrotron emission and to the inverse-Compton component of the afterglow. Theory predicts that in a large fraction of the cases where synchrotron GeV prompt emission would be detected by GLAST, inverse-Compton photons should be detected as well at high energies >~10 GeV. Therefore GLAST will enable a more precise test of the high-energy emission mechanism. Finally, we show that the contribution of GRBs to the flux of the extragalactic gamma-ray background measured with EGRET is at least 0.01% and likely around 0.1%.

We investigate the interaction of nonthermal ions (protons and nuclei) accelerated in an ultrarelativistic blastwave with the pre-existing magnetic field of the medium into which the blastwave propagates. While particle acceleration processes such as diffusive shock acceleration can accelerate ions and electrons, the accelerated electrons suffer larger radiative losses. Under certain conditions, the ions can attain higher energies and reach farther ahead of the shock than the electrons, and so the nonthermal particles can be partially charge-separated. To compensate for the charge separation, the upstream plasma develops a return current, which, as it flows across the magnetic field, drives transverse acceleration of the upstream plasma and a growth of density contrast in the shock upstream. If the density contrast is strong by the time the fluid is shocked, vorticity is generated at the shock transition. The resulting turbulence can amplify the post-shock magnetic field to the levels inferred from gamma-ray burst afterglow spectra and light curves. Therefore, since the upstream inhomogeneities are induced by the ions accelerated in the shock, they are generic even if the blastwave propagates into a medium of uniform density. We speculate about the global structure of the shock precursor, and delineate several distinct physical regimes that are classified by an increasing distance from the shock and, correspondingly, a decreasing density of nonthermal particles that reach that distance.

Massive stars end their short lives in spectacular explosions, supernovae, that synthesize new elements and drive galaxy evolution. Throughout history supernovae were discovered chiefly through their delayed optical light, preventing observations in the first moments (hours to days) following the explosion. As a result, the progenitors of some supernovae and the events leading up to their violent demise remain intensely debated. Here we report the serendipitous discovery of a supernova at the time of explosion, marked by an extremely luminous X-ray outburst. We attribute the outburst to the break-out of the supernova shock-wave from the progenitor, and show that the inferred rate of such events agrees with that of all core-collapse supernovae. We forecast that future wide-field X-ray surveys will catch hundreds of supernovae each year in the act of explosion, and thereby enable crucial neutrino and gravitational wave detections that may ultimately unravel the explosion mechanism.

The merging galaxy cluster 1E 0657-56, known as the "bullet cluster," is one of the hottest clusters known. The X-ray emitting plasma exhibits bow-shock-like temperature and density jumps. The segregation of this plasma from the peaks of the mass distribution determined by gravitational lensing has been interpreted as a direct proof of collisionless dark matter. If the high shock speed inferred from the shock jump conditions equals the relative speed of the merging CDM halos, however, this merger is predicted to be such a rare event in a LCDM universe that observing it presents a possible conflict with the LCDM model. We examined this question using high resolution, 2D simulations of gas dynamics in cluster collisions to analyze the relative motion of the clusters, the bow shock, and the contact discontinuity, and relate these to the X-ray data for the bullet cluster. We find that the velocity of the fluid shock need not equal the relative velocity of the CDM components. An illustrative simulation finds that the present relative velocity of the CDM halos is 16% lower than that of the shock. While this conclusion is sensitive to the detailed initial mass and gas density profiles of the colliding clusters, such a decrease of the inferred halo relative velocity would significantly increase the likelihood of finding 1E 0657-56 in a LCDM universe. (Conference proceedings based on a poster at Bash Symposium 2007)

We present optical spectra of the host galaxies of the short-duration gamma-ray burst GRB 070429B and the likely short-duration with extended emission GRB 070714B. In both cases, we find a single emission line that we identify as O II lambda 3727 at z ~ 0.9. Both events are more distant than any previous short-duration GRB with a secure host association from the sub-arcsecond position of an optical afterglow. GRBs 070429B and 070714B provide strong evidence in support of our previous claims in Berger et al. that a significant fraction of short-duration hosts (>~ 33%) reside at z > 0.7. We discuss the implications of the existence this population on the energetics of short-duration GRBs, as well as on progenitor models. In the context of the degenerate binary merger scenario, such events require progenitor systems with a range of lifetimes and disfavor progenitor models with a long, narrow lifetime distribution.

The gamma-ray burst (GRB) 070201 was a bright short-duration hard-spectrum GRB detected by the Inter-Planetary Network (IPN). Its error quadrilateral, which has an area of 0.124 sq. deg, intersects some prominent spiral arms of the nearby M31 (Andromeda) galaxy. Given the properties of this GRB, along with the fact that LIGO data argues against a compact binary merger origin in M31, this GRB is an excellent candidate for an extragalactic Soft Gamma-ray Repeater (SGR) giant flare, with energy of 1.4x10^45 erg. Analysis of ROTSE-IIIb visible light observations of M31, taken 10.6 hours after the burst and covering 42% of the GRB error region, did not reveal any optical transient down to a limiting magnitude of 17.1. We inspected archival and proprietary XMM-Newton X-ray observations of the intersection of the GRB error quadrilateral and M31, obtained about four weeks prior to the outburst, in order to look for periodic variable X-ray sources. No SGR or Anomalous X-ray Pulsar (AXP) candidates (periods in range 1 to 20 s) were detected. We discuss the possibility of detecting extragalactic SGRs/AXPs by identifying their periodic X-ray light curves. Our simulations suggest that the probability of detecting the periodic X-ray signal of one of the known Galactic SGRs/AXPs, if placed in M31, is about 10% (50%), using 50 ks (2 Ms) XMM-Newton exposures.

GRB 070610 is a typical high-energy event with a duration of 5s.Yet within the burst localization we detect a highly unusual X-ray and optical transient, SwiftJ195509.6+261406. We see high amplitude X-ray and optical variability on very short time scales even at late times. Using near-infrared imaging assisted by a laser guide star and adaptive optics, we identified the counterpart of SwiftJ195509.6+261406. Late-time optical and near-infrared imaging constrain the spectral type of the counterpart to be fainter than a K-dwarf assuming it is of Galactic origin. It is possible that GRB 070610 and Swift J195509.6+261406 are unrelated sources. However, the absence of a typical X-ray afterglow from GRB 070610 in conjunction with the spatial and temporal coincidence of the two motivate us to suggest that the sources are related. The closest (imperfect) analog to Swift J195509.6+261406 is V4641 Sgr, an unusual black hole binary. We suggest that Swift J195509.6+261406 along with V4641 Sgr define a sub-class of stellar black hole binaries -- the fast X-ray novae. We further suggest that fast X-ray novae are associated with bursts of gamma-rays. If so, GRB 070610 defines a new class of celestial gamma-ray bursts and these bursts dominate the long-duration GRB demographics

Because conical segments of quasispherical ultrarelativistic blastwaves are causally disconnected on angular scales larger than the blastwave inverse Lorentz factor, astrophysical blastwaves can sustain initial anisotropy, imprinted by the process that drives the explosion, while they remain relativistic. We show that initial angular energy fluctuations in ultrarelativistic blastwaves imply a production of vorticity in the blastwave, and calculate the vortical energy production rate. In gamma-ray burst (GRB) afterglows, the number of vortical eddy turnovers as the shocked fluid crosses the blastwave shell is about unity for marginally nonlinear anisotropy. Thus the anisotropy must be nonlinear to explain the magnetic energy density inferred in measured GRB spectra.

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

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

We present detailed optical, near-infrared, and radio observations of the X-ray flash 050416a obtained with Palomar and Siding Springs Observatories as well as HST and the VLA, placing this event among the best-studied X-ray flashes to date. In addition, we present an optical spectrum from Keck LRIS from which we measure the redshift of the burst, z=0.6528. At this redshift the isotropic-equivalent prompt energy release was about 10^51 erg, and using a standard afterglow synchrotron model we find that the blastwave kinetic energy is a factor of 10 larger, E_K,iso ~ 10^52 erg. The lack of an observed jet break to t ~ 20 days indicates that the opening angle is larger than 7 deg and the total beaming-corrected relativistic energy is larger than 10^50 erg. We further show that the burst produced a strong radio flare at t ~ 40 days accompanied by an observed flattening in the X-ray band which we attribute to an abrupt circumburst density jump or an episode of energy injection (either from a refreshed shock or off-axis ejecta). Late-time observations with HST show evidence for an associated supernova with peak optical luminosity roughly comparable to that of SN 1998bw. Next, we show that the host galaxy of XRF 050416a is actively forming stars at a rate of at least 2 M_solar per year with a luminosity of L_B ~ 0.5L* and metallicity of Z ~ 0.2-0.8 Z_solar. Finally, we discuss the nature of XRF 050416a in the context of short-hard gamma-ray bursts and under the framework of off-axis and dirty fireball models for X-ray flashes.

The redshift distribution of the short-duration GRBs is a crucial, but currently fragmentary, clue to the nature of their progenitors. Here we present optical observations of nine short GRBs obtained with Gemini, Magellan, and the Hubble Space Telescope. We detect the afterglows and host galaxies of two short bursts, and host galaxies for two additional bursts with known optical afterglow positions, and five with X-ray positions (<6'' radius). In eight of the nine cases we find that the most probable host galaxies are faint, R~23-26.5 mag, and are therefore starkly different from the first few short GRB hosts with R~17-22 mag and z<0.5. Indeed, we measure spectroscopic redshifts of z~0.4-1.1 for the four brightest hosts. A comparison to large field galaxy samples, as well as the hosts of long GRBs and previous short GRBs, indicates that the fainter hosts likely reside at z>1. Our most conservative limit is that at least half of the five hosts without a known redshift reside at z>0.7 (97% confidence level), suggesting that about 1/3-2/3 of all short GRBs originate at higher redshifts than previously determined. This has two important implications: (i) We constrain the acceptable age distributions to a wide lognormal (sigma>1) with tau~4-8 Gyr, or to a power law, P(tau)~tau^n, with -1<n<0; and (ii) the inferred isotropic energies, E_{gamma,iso}~10^50-10^52 erg, are significantly larger than ~10^48-10^49 erg for the low redshift short GRBs, indicating a large spread in energy release or jet opening angles. Finally, we re-iterate the importance of short GRBs as potential gravitational wave sources and find a conservative Advanced LIGO detection rate of ~2-6 yr^-1.

Shock waves driven in the intergalactic medium during the merging of galaxy clusters have been observed in X-ray imaging and spectroscopy. Fluid motions inferred from the shock strength and morphology can be compared to the cold dark matter (CDM) distribution inferred from gravitational lensing. A detailed reconstruction of the CDM kinematics, however, must take into account the nontrivial response of the fluid intracluster medium to the collisionless CDM motions. We have carried out two-dimensional simulations of gas dynamics in cluster collisions. We analyze the relative motion of the clusters, the bow shock wave, and the contact discontinuity and relate these to X-ray data. We focus on the "bullet cluster," 1E 0657-56, a near head-on collision of unequal-mass clusters, for which the gas density and temperature jumps across the prominent bow shock imply a high shock velocity 4,700 km/s. The velocity of the fluid shock has been widely interpreted as the relative velocity of the CDM components. This need not be the case, however. An illustrative simulation finds that the present relative velocity of the CDM halos is 16% lower than that of the shock. While this conclusion is sensitive to the detailed initial mass and gas density profile of the colliding clusters, such a decrease of the inferred halo relative velocity would increase the likelihood of finding 1E 0657-56 in a LambdaCDM universe.

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.

Short-hard and long-soft gamma-ray bursts (GRBs) are two distinct phenomena, but their prompt and afterglow emission show many similarities. This suggests that two different progenitor systems lead to similar physical processes and that the prompt and afterglow observations of short-hard GRBs (SHBs) can be examined using models of long GRBs. Here, I discuss three conclusions that can be drawn from SHB observations. I show that the lower limit on the Lorentz factor of SHBs is typically ``only'' 10-50, significantly lower than that of long GRBs. SHBs with observed X-ray afterglow after 1 day are found to be roughly as efficient as long GRBs in converting the outflow energy into prompt gamma-rays. Finally, I examine the origin of SHBs with X-ray dark afterglows and find that the most plausible explanation is that these SHBs exploded in extremely low density environment ($n \lesssim 10^{-5} cm^{-3}$)

Collisionless shocks are commonly argued to be the sites of cosmic ray (CR) acceleration. We study the influence of CRs on weakly magnetized relativistic collisionless shocks and apply our results to external shocks in gamma-ray burst (GRB) afterglows. The common view is that the transverse Weibel instability (TWI) generates a small-scale magnetic field that facilitates collisional coupling and thermalization in the shock transition. The TWI field is expected to decay rapidly, over a finite number of proton plasma skin depths from the transition. However, the synchrotron emission in GRB afterglows suggests that a strong and persistent magnetic field is present in the plasma that crosses the shock; the origin of this field is a key open question. Here we suggest that the common picture involving TWI demands revision. Namely, the CRs drive turbulence in the shock upstream on scales much larger than the skin depth. This turbulence generates a large-scale magnetic field that quenches TWI and produces a magnetized shock. The new field efficiently confines CRs and enhances the acceleration efficiency. The CRs modify the shocks in GRB afterglows at least while they remain relativistic. The origin of the magnetic field that gives rise to the synchrotron emission is plausibly in the CR-driven turbulence. We do not expect ultrahigh energy cosmic ray production in external GRB shocks.

Gamma ray bursts (GRBs) are known to come in two duration classes, separated at ~2 s. Long bursts originate from star forming regions in galaxies, have accompanying supernovae (SNe) when near enough to observe and are likely caused by massive-star collapsars. Recent observations show that short bursts originate in regions within their host galaxies with lower star formation rates consistent with binary neutron star (NS) or NS - black hole (BH) mergers. Moreover, although their hosts are predominantly nearby galaxies, no SNe have been so far associated with short GRBs. We report here on the bright, nearby GRB 060614 that does not fit in either class. Its ~102 s duration groups it with long GRBs, while its temporal lag and peak luminosity fall entirely within the short GRB subclass. Moreover, very deep optical observations exclude an accompanying supernova, similar to short GRBs. This combination of a long duration event without accompanying SN poses a challenge to both a collapsar and merging NS interpretation and opens the door on a new GRB classification scheme that straddles both long and short bursts.

The bright short-hard GRB 051103 was triangulated by the inter-planetary network and found to occur in the direction of the nearby M81/M82 galaxy group. Given its possible local-Universe nature, we searched for an afterglow associated with this burst. We observed the entire 3-sigma error quadrilateral using the Palomar 60-inch robotic telescope and the Very Large Array (VLA) about three days after the burst. We used the optical and radio observations to constrain the flux of any afterglow related to this burst, and to show that this burst is not associated with a typical supernova out to z~0.15. Our optical and radio observations, along with the Konus/Wind gamma-ray energy and light curve are consistent with this burst being a giant flare of a Soft Gamma-ray Repeater (SGR) within the M81 galaxy group. Furthermore, we find a star forming region associated with M81 within the error quadrilateral of this burst which supports the SGR hypothesis. If confirmed, this will be the first case of a soft gamma-ray repeater outside the local group.

We present detailed optical, X-ray and radio observations of the bright afterglow of the short gamma-ray burst 051221a obtained with Gemini, Swift/XRT, and the Very Large Array, as well as optical spectra from which we measure the redshift of the burst, z=0.5464. At this redshift the isotropic-equivalent prompt energy release was about 1.5 x 10^51 erg, and using the standard afterglow synchrotron model we find that the blastwave kinetic energy is similar, E_K,iso ~ 8.4 x 10^51 erg. An observed jet break at t ~ 5 days indicates that the opening angle is ~ 7 degrees and the total beaming-corrected energy is therefore ~ 2.5 x 10^49 erg, comparable to the values inferred for previous short GRBs. We further show that the burst experienced an episode of energy injection by a factor of 3.4 between t=1.4 and 3.4 hours, which was accompanied by reverse shock emission in the radio band. This result provides continued evidence that the central engines of short GRBs may be active significantly longer than the duration of the burst and/or produce a wide range of Lorentz factors. Finally, we show that the host galaxy of GRB051221a is actively forming stars at a rate of about 1.6 M_solar/yr, but at the same time exhibits evidence for an appreciable population of old stars (~ 1 Gyr) and near solar metallicity. The lack of bright supernova emission and the low circumburst density (n ~ 10^-3 cm^-3) continue to support the idea that short bursts are not related to the death of massive stars and are instead consistent with a compact object merger. Given that the total energy release is a factor of ~ 10 larger than the predicted yield for a neutrino annihilation mechanism, this suggests that magnetohydrodynamic processes may be required to power the burst.

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.

Some GRB afterglow light curves show significant variability, which often includes episodes of rebrightening. This was attributed in several cases to large fluctuations in the external density. Here we examine the effect of a sharp increase in the external density on the light curve by considering, for the first time, a full treatment of both the hydrodynamic evolution and the radiation in this scenario. To this end we develop a semi-analytic model for the light curve and carry out several elaborate numerical simulations using a 1D hydrodynamic code together with a synchrotron radiation code.Two spherically symmetric cases are explored in detail: a density jump in a uniform external medium, and a wind termination shock. Contrary to previous works, we find that even a sharp (step function) and large (by a factor of $a>>1$) increase in the external density does not produce sharp features in the light curve, and cannot account for significant temporal variability. For a wind termination shock, the light curve smoothly transitions between the asymptotic power-laws over about one decade in time, and there is no rebrightening in the optical or X-rays to serve as a clear observational signature. For a sharp jump in a uniform density we find that the maximal deviation $\Delta\alpha_max$ of the temporal decay index $\alpha$ from its asymptotic value, is bounded (e.g., $\Delta\alpha_max <0.4$ for $a=10$); $\Delta\alpha_max$ slowly increases with $a$, converging to $\Delta\alpha_{max} \approx 1$ at very large $a$ values. Therefore, no optical rebrightening is expected here as well. In the X-rays, while the asymptotic flux is unaffected by the density jump, the fluctuations in $\alpha$ are found to be comparable to those in the optical. Finally, we discuss the implications of our results for the origin of the fluctuations in several afterglows.

The recent discovery of the first four afterglows of short-hard gamma-ray bursts (SHBs), suggests that they typically result from long-lived progenitor systems. The most popular progenitor model invokes the merger of either double neutron star (DNS) binaries or neutron star-black hole (NS-BH) systems. Such events are strong sources of gravitational-waves (GWs) and might be detected by ground-based GW observatories. In this work we combine the census of SHB observations with refined theoretical analysis to perform a critical evaluation of the compact binary model. We then explore the implications for GW detection of these events. Beginning from the measured star-formation rate through cosmic time, we consider what intrinsic luminosity and lifetime distributions can reproduce the known SHB redshifts and luminosities as well as the peak flux distribution of the large BATSE SHB sample. We find that: (1) The typical progenitor lifetime is long. Assuming lognormal lifetime distribution the typical lifetime is >4[1] Gyr (2[3]\sigma c.l.). If the lifetime distribution is a power-law with index \eta then \eta>-0.5 [-1] (2[3] \sigma c.l.). This result is difficult to reconcile with the properties of the observed galactic DNS population, suggesting that if SHBs do result from DNS mergers then the observed galactic binaries do not represent the cosmic one. (2) We find that the local rate of SHBs is larger than 10Gpc^-3/yr and may be higher by several orders of magnitude, significantly above previous estimates. (3) Assuming that SHBs do result from compact binaries, our predictions for the LIGO and VIRGO event rates are encouraging: The chance for detection by current facilities is not negligible, while a coincident detection of GW and electromagnetic radiation from an SHB is guaranteed for next-generation observatories.

We present the discovery of radio afterglow emission from the high redshift (z=6.295) burst GRB 050904. The peak flux density for this burst is similar to typical low-redshift gamma-ray bursts (GRB). We further show that beyond a redshift of order unity, the flux density of radio afterglows are largely insensitive to redshift, consistent with predictions. By combining the existing X-ray, near-infrared and radio measurements, we derive estimates for the kinetic energy and opening angle of the blast wave, and for the density of the circumburst medium into which it expands. Both the kinetic and radiated energy indicate that GRB 050904 was an unusally energetic burst (10**52 erg). More importantly, we are able to make an in situ measurement of the density structure of the circumburst medium. We conclude that GRB 050904 exploded into a constant density medium with n=680 cm**-3, which is two orders of magnitude above the nominal value for low-redshift GRBs. The next generation of centimeter (EVLA) and millimeter radio instuments (ALMA) will be able to routinely detect events like GRB 050904 and use them to study magnetic fields, and the atomic and molecular gas in the high redshift Universe.

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.

We present a collation of the available data on the opening angles of jets in X-ray binaries, which in most cases are small (less than 10 degrees). Under the assumption of no confinement, we calculate the Lorentz factors required to produce such small opening angles via the transverse relativistic Doppler effect. The derived Lorentz factors, which are in most cases lower limits, are found to be large, with a mean greater than 10, comparable to those estimated for AGN and much higher than the commonly-assumed values for X-ray binaries of 2 to 5. Jet power constraints do not in most cases rule out such high Lorentz factors. The upper limits on the opening angles show no evidence for smaller Lorentz factors in the steady jets of Cygnus X-1 and GRS 1915+105. In those sources in which deceleration has been observed (notably XTE J1550-564 and Cygnus X-3), some confinement of the jets must be occurring, and we briefly discuss possible confinement mechanisms. It is however possible that all the jets could be confined, in which case the requirement for high bulk Lorentz factors can be relaxed.

We present late-time radio observations of 68 local Type Ibc supernovae, including six events with broad optical absorption lines (``hypernovae''). None of these objects exhibit radio emission attributable to off-axis gamma-ray burst jets spreading into our line-of-sight. Comparison with our afterglow models reveals the following conclusions: (1) Less than $\sim 10%$ of Type Ibc supernovae are associated with typical gamma-ray bursts initially directed away from our line-of-sight; this places an empirical constraint on the GRB beaming factor of $<f_b^{-1}> \lesssim 10^4$ corresponding to an average jet opening angle, $\theta_j \gtrsim 0.8$ degrees. (2) This holds in particular for the broad-lined supernovae (SNe 1997dq, 1997ef, 1998ey, 2002ap, 2002bl and 2003jd) which have been argued to host GRB jets. Our observations reveal no evidence for typical (or even sub-energetic) GRBs and rule out the scenario in which every broad-lined SN harbors a GRB at the 84% confidence level. Their large photospheric velocities and asymmetric ejecta (inferred from spectropolarimetry and nebular spectroscopy) appear to be characteristic of the non-relativistic supernova explosion and do not necessarily imply the existence of associated GRB jets.

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

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

Models for the synchrotron emission of gamma-ray burst afterglows suggest that the magnetic field is generated in the shock wave that forms as relativistic ejecta plow through the circum-burst medium. Transverse Weibel instability efficiently generates magnetic fields near equipartition with the post-shock energy density. The detailed saturated state of the instability, as seen in particle-in-cell simulations, consists of magnetically self-pinched current filaments. The filaments are parallel to the direction of propagation of the shock and are about a plasma skin depth in radius, forming a quasi--two-dimensional structure. We use a rudimentary analytical model to argue that the Weibel filaments are unstable to a kink-like mode, which destroys their quasi--two-dimensional structure. For wavelengths longer than than the skin depth, the instability grows at the rate equal to the speed of light divided by the wavelength. We calculate the transport of collisionless test particles in the filaments experiencing the instability and show that the particles diffuse in energy. This diffusion marks the beginning of thermalization in the shock transition layer, and causes initial magnetic field decay as particles escape from the filaments. We discuss the implications of these results for the structure of the shock and the polarization of the afterglow.

It is generally accepted that magnetic fields generated in the nonlinear development of the transverse Weibel instability provide effective collisionality in unmagnetized collisionless shocks. Recently, extensive two and three dimensional simulations improved our understanding of the growth and saturation of the instability in colliding plasma shells. However, the steady-state structure of the shock wave transition layers remains poorly understood. We use basic physical considerations and order-of-magnitude arguments to study the steady state structure in relativistic unmagnetized collisionless shocks in pair plasmas. The shock contains an electrostatic layer resulting from the formation of stationary, magnetically-focused current filaments. The filaments form where the cold upstream plasma and the counterstreaming thermal plasma interpenetrate. The filaments are not entirely neutral and strong electrostatic fields are present. Most of the downstream particles cannot cross this layer into the upstream because they are trapped by the electrostatic field. We identify the critical location in the shock transition layer where the electromagnetic field ceases to be static. At this location, the degree of charge separation in the filaments reaches a maximum value, the current inside the filaments comes close to the Alfv\'en limit, and the phase space distribution function starts to isotropize. We argue that the radius of the current filaments upstream of the critical location is about twice the upstream plasma skin depth. Finally, we show that some downstream particles cross the electrostatic layer and run ahead of the shock into the preshock medium without causing instability. These particles may play an important role in particle acceleration.

We present optical, near-IR, and radio follow up of sixteen Swift bursts, including our discovery of nine afterglows and a redshift determination for three. These observations, supplemented by data from the literature, provide an afterglow recovery rate of 60% in the optical/near-IR, much higher than in previous missions (BeppoSAX, HETE-2, INTEGRAL, and IPN). The optical/near-IR afterglows of Swift events are on average 1.7 mag fainter at t=12 hr than those of previous missions. The X-ray afterglows are similarly fainter compared to those of pre-Swift bursts. In the radio the limiting factor is the VLA threshold and the detection rate for Swift bursts is similar to that for past missions. The redshift distribution of pre-Swift bursts peaked at z~1, whereas the five Swift bursts with measured redshifts are distributed evenly between 1.3 and 3.2. From these results we conclude that (i) the pre-Swift distributions were biased in favor of bright events and low redshift events, (ii) the higher sensitivity and accurate positions of Swift result in a better representation of the true burst redshift and brightness distributions (which are higher and dimmer, respectively), and (iii) as many as 1/3 of the bursts can be optically dark, as a result of a high redshift and/or dust extinction. We remark that the apparent lack of low redshift, low luminosity Swift bursts, and the lower event rate compared to pre-launch estimates (90 vs. 150 per year), are the result of a threshold that is similar to that of BATSE. In view of these inferences, afterglow observers may find it advisable to make significant changes in follow up strategies of Swift events. [abridged]

Over a decade ago it was established that the remarkable high energy transients, known as soft gamma-ray repeaters (SGRs), are a Galactic population and originate from neutron stars with intense (<~ 10^15 G) magnetic fields ("magnetars"). On 27 December 2004 a giant flare (fluence >~ 0.3 erg/cm^2) was detected from SGR 1806-20. Here we report the discovery of a fading radio counterpart. We began a monitoring program from 0.2GHz to 250GHz and obtained a high resolution 21-cm radio spectrum which traces the intervening interstellar neutral Hydrogen clouds. Analysis of the spectrum yields the first direct distance measurement of SGR 1806-20. The source is located at a distance greater than 6.4 kpc and we argue that it is nearer than 9.8 kpc. If true, our distance estimate lowers the total energy of the explosion and relaxes the demands on theoretical models. The energetics and the rapid decay of the radio source are not compatible with the afterglow model that is usually invoked for gamma-ray bursts. Instead we suggest that the rapidly decaying radio emission arises from the debris ejected during the explosion.

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.

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.

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.

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 two leading models for the structure of GRB jets are the uniform jet model and the universal structured jet (USJ) model. In the latter, all GRB jets are intrinsically identical and the energy per solid angle drops as the inverse square of the angle from the jet axis. The simplicity of the USJ model gives it a strong predictive power, including a specific prediction for the observed GRB distribution as a function of both the redshift $z$ and the viewing angle $\theta$. We show that the current sample of GRBs with known $z$ and estimated $\theta$ does not agree with the predictions of the USJ model. This can be best seen for a relatively narrow range in $z$, in which the USJ model predicts that most GRBs should be near the upper end of the observed range in $\theta$, while in the observed sample most GRBs are near the lower end of that range. Since the current sample is very inhomogeneous (i.e. involves many different detectors), it should be taken with care and cannot be used to rule out the USJ model. Nevertheless, this sample strongly disfavors the USJ model. Comparing the prediction for the observed GRB distribution both in $\theta$ and in $z$, with a larger and more homogeneous GRB sample, like the one expected from Swift, would either clearly rule out the USJ model, or alternatively, provide a strong support for it. The test presented here is general, and can be used to test any model that predicts both a luminosity function and a luminosity-angle relation.

Recent advances in early detection and detailed monitoring of GRB afterglows have revealed variability in some afterglow light curves. One of the leading models for this behavior is the patchy shell model. This model attributes the variability to random angular fluctuations in the relativistic jet energy. These an-axisymmetric fluctuations should also impose variations in the degree and angle of polarization that are correlated to the light curve variability. In this letter we present a solution of the light curve and polarization resulting from a given spectrum of energy fluctuations. We compare light curves produced using this solution to the variable light curve of GRB 021004 and we show that the main features in both the light curve and the polarization fluctuations are very well reproduced by this model. We use our results to draw constraints on the characteristics of the energy fluctuations that might have been present in GRB 021004.

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.

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.

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.

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.

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.

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.

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.

Complex adaptive systems have been the subject of much recent attention. It is by now well-established that members (`agents') tend to self-segregate into opposing groups characterized by extreme behavior. However, while different social and biological systems manifest different payoffs, the study of such adaptive systems has mostly been restricted to simple situations in which the prize-to-fine ratio, $R$, equals unity. In this Letter we explore the dynamics of evolving populations with various different values of the ratio $R$, and demonstrate that extreme behavior is in fact {\it not} a generic feature of adaptive systems. In particular, we show that ``confusion'' and ``indecisiveness'' take over in times of depression, in which case cautious agents perform better than extreme ones.

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.

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.