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.

Magnetized jets in GRBs and AGNs are thought to be efficient accelerators of particles, however, the process responsible for the acceleration is still a matter of active debate. In this work, we study the kink-instability in non-rotating force-free jets using first-principle particle-in-cell simulations. We obtain similar overall evolution of the instability as found in MHD simulations. The instability first generates large scale current sheets, which at later times break up into small-scale turbulence. Reconnection in these sheets proceeds in the strong guide field regime, which results in a formation of steep power laws in the particle spectra. Later evolution shows heating of the plasma, which is driven by small-amplitude turbulence induced by the kink instability. These two processes energize particles due to a combination of ideal and non-ideal electric fields.

We study the evolution of kink instability in a force-free, non-rotating plasma column of high magnetization. The main dissipation mechanism is identified as reconnection of magnetic field-lines with various intersection angles, driven by the compression of the growing kink lobes. We measure dissipation rates ${\rm d} U_{B\phi}/{{\rm d}t} \approx -0.1 U_{B\phi}/\tau$, where $\tau$ is the linear growth time of the kink instability. This value is consistent with the expansion velocity of the kink mode, which drives the reconnection. The relaxed state is close to a force-free Taylor state. We constraint the energy of that state using considerations from linear stability analysis. Our results are important for understanding magnetic field dissipation in various extreme astrophysical objects, most notably in relativistic jets. We outline the evolution of the kink instability in such jets and derive constrains on the conditions that allow for the kink instability to grow in these systems.

Globular clusters (GCs) may harbour intermediate-mass black holes (IMBHs) at their centres. In these dynamically active environments stellar-mass black holes (SBHs) sink to the center soon after formation, due to dynamical friction and start interacting among themselves and with the central IMBH. Likely, some of the SBHs will form bound systems with the IMBH. A fraction of those will be triple systems composed of binary SBHs and the IMBH acting as a third distant perturber. If the SBH binary orbit is sufficiently inclined it can develop Lidov-Kozai (LK) oscillations, which can drive the system to high eccentricities and eventually to a merger due to gravitational wave (GW) emission on short timescales. In this work, we focus on the dynamics of the IMBH-SBH-SBH triples and illustrate that these systems can be possible sources of GWs. A distinctive signature of this scenario is that a considerable fraction of these mergers are highly eccentric when entering the LIGO band ($10$ Hz). Assuming that $\sim 20\%$ of GCs host IMBHs and a GC density in the range $n_{_{\rm GC}}=0.32$-$2.31\,\mathrm{Mpc}^{-3}$, we have estimated a rate $\Gamma=0.06$-$0.46\,\mathrm{Gpc}^{-3}\,\mathrm{yr}^{-1}$ of these events. This suggests that dynamically-driven binary SBH mergers in this scenario could contribute to the merger events observed by LIGO/VIRGO. Full $N$-body simulations of GCs harbouring IMBHs are highly desirable to give a more precise constrain on this scenario.

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.

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.

Long-duration, spectrally-soft Gamma-Ray Bursts (GRBs) are associated with Type Ic Core Collapse (CC) Supernovae (SNe), and thus arise from the death of massive stars. In the collapsar model, the jet launched by the central engine must bore its way out of the progenitor star before it can produce a GRB. Most of these jets do not break out, and are instead "choked" inside the star, as the central-engine activity time, $t_{\rm e}$, is not long enough. Modelling the long-soft GRB duration distribution assuming a power-law distribution for their central-engine activity times, $\propto t_{\rm e}^{-\alpha}$ for $t_{\rm e}>t_{\rm b}$, we find a steep distribution ($\alpha\sim4$) and a typical GRB jet breakout time of $t_{\rm b}\sim 60\text{ s}$ in the star's frame. The latter suggests the presence of a low-density, extended envelope surrounding the progenitor star, similar to that previously inferred for low-luminosity GRBs. Extrapolating the range of validity of this power law below what is directly observable, to $t_{\rm e}<t_{\rm b}$, by only a factor of $\sim$4-5 produces enough events to account for all Type Ib/c SNe. Such extrapolation is necessary to avoid fine-tuning the distribution of central engine activity times with the breakout time, which are presumably unrelated. We speculate that central engines launching relativistic jets may operate in all Type Ib/c SNe. In this case, the existence of a common central engine would imply that (i) the jet may significantly contribute to the energy of the SN; (ii) various observational signatures, like the asphericity of the explosion, could be directly related to jet's interaction with the star.

Core-collapse supernovae (SNe) of Types Ib and Ic arise from hydrogen-stripped stars, while the latter are also stripped of their helium. Both SN types have a similar temporal evolution, suggesting broadly similar progenitors. However, while some Ic SNe are associated with relativistic jets of gamma-ray bursts (GRBs), no GRB has yet been found in association with Type Ib SNe. Here we find that, even if GRB-like central engines operate in both SNe Ib and Ic, different properties of their envelopes may accommodate this potential tension. In particular, we focus on the case of a low-mass, extended envelope surrounding the progenitor star (as produced, for example, by strong mass losses prior to explosion). If the envelopes of Type Ib SNe are sufficiently massive ($M_{\rm ext}\sim(0.3-1)M_\odot$) and extended ($R_{\rm ext}\sim 10^{13}\text{ cm}$), we show that (i) GRB-like jets cannot break out of the star; (ii) the SN light curve is compatible with current observations. Different envelope properties of Type Ib SNe with respect to (at least some) Type Ic SNe may be connected to the presence of a helium layer surrounding the star.

Energy deposition by active galactic nuclei jets into the ambient medium can affect galaxy formation and evolution, the cooling of gas flows at the centres of galaxy clusters, and the growth of the supermassive black holes. However, the processes that couple jet power to the ambient medium and determine jet morphology are poorly understood. For instance, there is no agreement on the cause of the well-known Fanaroff-Riley (FR) morphological dichotomy of jets, with FRI jets being shorter and less stable than FRII jets. We carry out global 3D magnetohydrodynamic simulations of relativistic jets propagating through the ambient medium. We show that the flat density profiles of galactic cores slow down and collimate the jets, making them susceptible to the 3D magnetic kink instability. We obtain a critical power, which depends on the galaxy core mass and radius, below which jets become kink-unstable within the core, stall, and inflate cavities filled with relativistically-hot plasma. Jets above the critical power stably escape the core and form powerful backflows. Thus, the kink instability controls the jet morphology and can lead to the FR dichotomy. The model-predicted dependence of the critical power on the galaxy optical luminosity agrees well with observations.

Relativistic jets naturally occur in astrophysical systems that involve accretion onto compact objects, such as core collapse of massive stars in gamma-ray bursts (GRBs) and accretion onto supermassive black holes in active galactic nuclei (AGN). It is generally accepted that these jets are powered electromagnetically, by the magnetised rotation of a central compact object. However, how they produce the observed emission and survive the propagation for many orders of magnitude in distance without being disrupted by current-driven non-axisymmetric instabilities is the subject of active debate. We carry out time-dependent 3D relativistic magnetohydrodynamic simulations of relativistic, Poynting flux dominated jets. The jets are launched self-consistently by the rotation of a strongly magnetised central compact object. This determines the natural degree of azimuthal magnetic field winding, a crucial factor that controls jet stability. We find that the jets are susceptible to two types of instability: (i) a global, external kink mode that grows on long time scales and causes the jets to bodily bend sideways. Whereas this mode does not cause jet disruption over the simulated distances, it substantially reduces jet propagation speed. We show, via an analytic model, that the growth of the external kink mode depends on the slope of the ambient medium density profile. In flat density distributions characteristic of galactic cores, an AGN jet may stall, whereas in stellar envelopes the external kink weakens as the jet propagates outward; (ii) a local, internal kink mode that grows over short time scales and causes small-angle magnetic reconnection and conversion of about half of jet electromagnetic energy flux into heat. Based on the robustness and energetics of the internal kink mode, we suggest that this instability is the main dissipation mechanism responsible for powering GRB prompt emission.

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

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

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

Observations of the FR I radio galaxy Centaurus A in radio, X-ray and gamma-ray bands provide evidence for lepton acceleration up to several TeV and clues about hadron acceleration to tens of EeV. Synthesising the available observational constraints on the physical conditions and particle content in the jets, inner lobes and giant lobes of Centaurus A, we aim to evaluate its feasibility as an ultra-high-energy cosmic-ray source. We apply several methods of determining jet power and affirm the consistency of various power estimates of ~ 1 x 10^43 erg s^-1. Employing scaling relations based on previous results for 3C 31, we estimate particle number densities in the jets, encompassing available radio through X-ray observations. Our model is compatible with the jets ingesting ~ 3 x 10^21 g s^-1 of matter via external entrainment from hot gas and ~ 7 x 10^22 g s^-1 via internal entrainment from jet-contained stars. This leads to an imbalance between the internal lobe pressure available from radiating particles and magnetic field, and our derived external pressure. Based on knowledge of the external environments of other FR I sources, we estimate the thermal pressure in the giant lobes as 1.5 x 10^-12 dyn cm^-2, from which we deduce a lower limit to the temperature of ~ 1.6 x 10^8 K. Using dynamical and buoyancy arguments, we infer ~ 440-645 Myr and ~ 560 Myr as the sound-crossing and buoyancy ages of the giant lobes respectively, inconsistent with their spectral ages. We re-investigate the feasibility of particle acceleration via stochastic processes in the lobes, placing new constraints on the energetics and on turbulent input to the lobes. The same 'very hot' temperatures that allow self-consistency between the entrainment calculations and the missing pressure also allow stochastic UHECR acceleration models to work.

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

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

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

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

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.

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

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

It is proposed that the prompt emission observed in bursts that exhibit a thermal component originates from relativistic radiation mediated shocks that form below the photosphere of the GRB outflow. It is argue that such shocks are expected to form in luminous bursts via collisions of shells that propagate with moderate Lorentz factors $\Gamma\lesssim 500$. Faster shells will collide above the photosphere to form collisionless shocks. We demonstrate that in events like GRB 090902B a substantial fraction of the explosion energy is dissipated below the photosphere, in a region of moderate optical depth $\tau\lesssim300$, whereas in GRB 080916C the major fraction of the energy dissipates above the photosphere. We show that under conditions anticipated in many GRBs, such relativistic radiation mediated shocks convect enough radiation upstream to render photon production in the shock transition negligible, unlike the case of shock breakout in supernovae. The resulting spectrum, as measured in the shock frame, has a relatively low thermal peak, followed by a broad, nonthermal component extending up to the KN limit.

Recent observations of M87 and some blazars reveal violent activity in small regions located at relatively large distances from the central engine. Motivated by these considerations, we study the hydrodynamic collimation of a relativistic cooling outflow using a semi-analytical model developed earlier. We first demonstrate that radiative cooling of the shocked outflow layer can lead to a focusing of the outflow and its reconfinement in a region having a very small cross-sectional radius. Such a configuration can produce rapid variability at large distances from the central engine via reflections of the converging recollimation shock. Possible applications of this model to TeV blazars are discussed. We then apply our model to M87. The low radiative efficiency of the M87 jet renders focusing unlikely. However, the shallow profile of the ambient medium pressure inferred from observations, results in extremely good collimation that can explain the reported variability of the X-ray flux emitted from the HST-1 knot.

A system of equations governing the structure of a steady, relativistic radiation dominated shock is derived, starting from the general form of the transfer equation obeyed by the photon distribution function. Closure is obtained by truncating the system of moment equations at some order. The anisotropy of the photon distribution function inside the shock is shown to increase with increasing shock velocity, approaching nearly perfect beaming at upstream Lorentz factors $\Gamma_{-}>>1$. Solutions of the shock equations are presented for some range of upstream conditions. These solutions are shown to converge as the truncation order is increased.

We consider some implications of the rapid X-ray and TeV variability observed in M87 and the TeV blazars. We outline a model for jet focusing and demonstrate that modest radiative cooling can lead to recollimation of a relativistic jet in a nozzle having a very small cross-sectional radius. Such a configuration can produce rapid variability at large distances from the central engine and may explain recent observations of the HST-1 knot in M87. Possible applications of this model to TeV blazars are discussed. We also discuss a scenario for the very rapid TeV flares observed with HESS and MAGIC in some blazars, that accommodates the relatively small Doppler factors inferred from radio observations.

A model is developed for the confinement and collimation of a baryon poor outflow by its surrounding medium. Both, confinement by kinetic pressure of a static corona, and confinement by the ram pressure of a supersonic wind emanating from a disk surrounding the inner source are considered. Solutions are presented for the structure of the shocked layers of a deflected baryon poor jet (BPJ) and exterior wind. The dependence of the opening angle of the BPJ on the parameters of the confining medium are carefully examined. It is found that the BPJ shock may either converge to the symmetry axis or diverge away from it, depending on the opening angle of the BPJ injection cone. In the latter case the inner flow exhibits a non-uniform structure, consisting of an ultra-relativistic core containing the unshocked BPJ enveloped by the slower, shocked BPJ layer. The implications of our results to the prompt GRB emission are briefly discussed.

We predict the redshift of the first observable (i.e., in our past light cone) Gamma Ray Burst (GRB) and calculate the GRB-rate redshift distribution of the Population III stars at very early times (z=20-60). Using the last 2 years of data from Swift we place an upper limit on the efficiency (\eta_{GRB}) of GRB production per solar mass from the first generation of stars. We find that the first observable GRB is most likely to have formed at redshift 60. The observed rate of extremely high redshift GRBs (XRGs) is a subset of a group of 15 long GRBs per year, with no associated redshift and no optical afterglow counterparts, detected by Swift. Taking this maximal rate we get that \eta_{GRB}<1.1~10^{-4} GRBs per solar mass in stars. A more realistic evaluation, e.g., taking a subgroup of 5% of the total sample of Swift gives an upper limit of \eta_{GRB}<3.2~10^{-5} GRBs per solar mass.

The gravitational-wave spectrum emitted by a non-axisymmetric torus rotating at an angular velocity $\Omega_T$, is derived in terms of a structure function representing a combination of sausage-tilt modes in the torus in the limit of an incompressible fluid. The analysis of the gravitational-wave spectrum is then applied to a model proposed recently, in which a highly magnetized torus interacts with a stellar mass, Kerr black hole via poloidal field lines that connect the torus and the horizon. The deformation of the torus results from global magnetic instabilities when the magnetic field strength inside the torus exceeds a few times $10^{15}$ Gauss. The dynamics of the system is illustrated using a non-MHD toy model. It is found that, quite generally, most of the gravitational-wave energy emitted is in the frequency range of sensitivity of LIGO and Virgo.