Spark gaps are likely the source of plasma in active black hole (BH) magnetospheres. In this paper, we present results of 1D general-relativistic particle-in-cell simulations of a starved BH magnetosphere with a realistic treatment of inverse Compton scattering and pair production, for a broad range of conditions, run times longer than in previous studies, and different setups. We find that following the initial discharge, the system undergoes gradual evolution over prolonged time until either, restoring the vacuum state or reaching a state of quasi-periodic oscillations, depending on the spectral shape and luminosity of the ambient radiation field. The oscillations occur near the null charge surface in cases where the global magnetospheric current is negative, and near the boundary of the simulation box when it is positive. Their amplitude and the resultant luminosity of TeV photons emitted from the gap depend sensitively on the conditions; for the cases studied here the ratio of TeV luminosity to the Blandford-Znajek power ranges from $10^{-5}$ to $10^{-2}$, suggesting that strong flares may be generated by moderate changes in disk emission. We also examined the dependence of the solution on the initial number of particles per cell (PPC) and found convergence for PPC of about 50 for the cases examined. At lower PPC values the pair multiplicity is found to be artificially high, affecting the solution considerably.

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.

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.

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.

Black hole - accretion disc systems are the central engines of relativistic jets from stellar to galactic scales. We numerically quantify the unsteady outgoing Poynting flux through the horizon of a rapidly spinning black hole endowed with a rotating accretion disc. The disc supports small-scale, concentric, flux tubes with zero net magnetic flux. Our General Relativistic force-free electrodynamics simulations follow the accretion onto the black hole over several hundred dynamical timescales in 3D. For the case of counter-rotating accretion discs, the average process efficiency reaches up to $\left\langle\epsilon\right\rangle\approx 0.43$, compared to a stationary energy extraction by the Blandford/Znajek process. The process efficiency depends on the cross-sectional area of the loops, i.e. on the product $l\times h$, where $l$ is the radial loop thickness and $h$ its vertical scale height. We identify a strong correlation between efficient electromagnetic energy extraction and the quasi-stationary setting of ideal conditions for the operation of the Blandford/Znajek process (e.g. optimal field line angular velocity and fulfillment of the so-called Znajek condition). Remarkably, the energy extraction operates intermittently (alternating episodes of high and low efficiency) without imposing any large-scale magnetic field embedding the central object. Scaling our results to supermassive black holes, we estimate that the typical variability timescale of the system is of the order of days to months. Such timescales may account for the longest variability scales of TeV emission observed, e.g. in M87.

Radiation mediated shocks (RMS) play a key role in shaping the early emission observed in many transients. In most cases, e.g., shock breakout in supernovae, llGRBs and neutron star mergers, the upstream plasma is devoid of radiation, and the photons that ultimately reach the observer are generated predominantly inside and downstream of the shock. Predicting the observed spectrum requires detailed calculations of the shock structure and thermodynamic state that account properly for the shock microphysics. We present results of self-consistent Monte-Carlo simulations of photon-starved RMS, that yield the shock structure and emission for a broad range of shock velocities, from sub-relativistic ($\beta_{sh} = 0.1$) to highly relativistic ($\Gamma_{sh} = 20$). Our simulations confirm that in relativistic RMS the immediate downstream temperature is regulated by exponential pair creation, ranging from $50$ keV at $\beta_{sh}=0.5$ to $200$ keV at $\Gamma_{sh}=20$. At lower velocities the temperature becomes sensitive to the shock velocity, with $kT\sim 0.5$ keV at $\beta_{sh}=0.1$. We also confirm that in relativistic shocks the opacity is completely dominated by newly created pairs, which has important implications for the breakout physics. We find the transition to pair dominance to occur at $\beta_{sh}=0.5$ roughly. In all cases examined, the spectrum below the $\nu F_\nu$ peak has been found to be substantially softer than the Planck distribution. This has important implications for the optical emission in fast and relativistic breakouts, and their detection. The applications to GRB 060218 and GRB 170817A are discussed.

Multi-messenger emissions from SN1987A and GW170817/GRB170817A suggest a Universe rife with multi-messenger transients associated with black holes and neutron stars. For LIGO-Virgo, soon to be joined by KAGRA, these observations promise unprecedented opportunities to probe the central engines of core-collapse supernovae (CC-SNe) and gamma-ray bursts. Compared to neutron stars, central engines powered by black hole-disk or torus systems may be of particular interest to multi-messenger observations by the relatively large energy reservoir $E_J$ of angular momentum, up to 29\% of total mass in the Kerr metric. These central engines are expected from relatively massive stellar progenitors and compact binary coalescence involving a neutron star. We review prospects of multi-messenger emission by catalytic conversion of $E_J$ by a non-axisymmetric disk or torus. Observational support for this radiation process is found in a recent identification of ${\cal E}\simeq (3.5\pm1)\%M_\odot c^2$ in Extended Emission to GW170817 at a significance of 4.2\,$\sigma$ concurrent with GRB170817A. A prospect on similar emissions from nearby CC-SNe justifies the need for all-sky blind searches of long duration bursts by heterogeneous computing.

GW170817/GRB170817A probably marks a double neutron star coalescence. Extended Emission $t_s\simeq (0.67\pm0.03)$\,s post-merger shows an estimated energy output ${\cal E}\simeq (3.5\pm1)\%M_\odot c^2$ determined by response curves to power-law signal injections, where $c$ is the velocity of light. It provides calorimetric evidence for a rotating black hole of $\sim 3M_\odot$, inheriting the angular momentum $J$ of the merged hyper-massive neutron star in the immediate aftermath of GW170817 following core-collapse about or prior to $t_s$. Core collapse greatly increases the central energy reservoir to $E_J\lesssim 1M_\odot c^2$, accounting for ${\cal E}$ even at modest efficiencies in radiating gravitational waves through a non axisymmetric thick torus. The associated multi-messenger output in ultra-relativistic outflows and sub-relativistic mass-ejecta is consistent with observational constraints from the GRB-afterglow emission of GRB170817A and accompanying kilonova.

Blue-shifted Broad Absorption Lines (BALs) detected in quasar's spectra are indicative of AGN outflows. We show, using 2D hydrodynamical simulations, that disruption of interstellar clouds by a fast AGN wind can lead to formation of cold, dense high speed blobs that give rise to broad absorption features in the transmission spectrum of the AGN continuum source. For a wind velocity of $0.1 c$ and sufficiently high cloud density ($n_{c}$ < $10^4$ cm$^{-3}$, depending on size), absorption troughs with velocities up to about $3000$ km s$^{-1}$ can be produced. For slower winds and/or lower cloud density the anticipated velocity of the absorbing clouds should be smaller.

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

Radio Galaxies have by now emerged as a new $\gamma$-ray emitting source class on the extragalactic sky. Given their remarkable observed characteristics, such as unusual gamma-ray spectra or ultrafast VHE variability, they represent unique examples to probe into the nature and physics of AGN in general. This review provides a compact summary of their observed characteristics at very high $\gamma$-ray energies (VHE; $> 100$ GeV) along with a discussion of their possible physics implications. A particular focus is given to a concise overview of fundamental concepts concerning the origin of variable VHE emission, including recent developments in black hole gap physics.

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.

We investigate the dynamics and emission of a starved magnetospheric region (gap) formed in the vicinity of a Kerr black hole horizon, using a new, fully general relativistic particle-in-cell code that implements Monte Carlo methods to compute gamma-ray emission and pair production through the interaction of pairs and gamma rays with soft photons emitted by the accretion flow. It is found that when the Thomson length for collision with disk photons exceeds the gap width, screening of the gap occurs through low-amplitude, rapid plasma oscillations that produce self-sustained pair cascades, with quasi-stationary pair and gamma-ray spectra, and with a pair multiplicity that increases in proportion to the pair production opacity. The gamma-ray spectrum emitted from the gap peaks in the TeV band, with a total luminosity that constitutes a fraction of about $10^{-5}$ of the corresponding Blandford-Znajek power. This stage is preceded by a prompt discharge phase of duration $\sim r_g/c$, during which the potential energy initially stored in the gap is released as a flare of curvature TeV photons. We speculate that the TeV emission observed in M87 may be produced by pair discharges in a spark gap.

Under conditions prevailing in certain classes of compact astrophysical systems, the active magnetosphere of a rotating black hole becomes charge-starved, giving rise to formation of a spark gap in which plasma is continuously produced. The plasma production process is accompanied by curvature and inverse Compton emission of gamma rays in the GeV-TeV band, that may be detectable by current and future experiments. The properties of the gap emission have been studied recently using a fully general relativistic model of a local steady gap. However, this model requires artificial adjustment of the electric current which is determined, in reality, by the global properties of the magnetosphere. In this paper we map the parameter regime in which steady gap solutions exist, using a steady-state gap model in Kerr geometry, and show that such solutions are allowed only under restrictive conditions that may not apply to most astrophysical systems. We further argue that even the allowed solutions are inconsistent with the global magnetospheric structure. We conclude that magnetospheric gaps are inherently intermittent, and point out that this may drastically change their emission properties.

We explore the physics of relativistic radiation mediated shocks (RRMSs) in the regime where photon advection dominates over photon generation. For this purpose, a novel iterative method for deriving a self-consistent steady-state structure of RRMS is developed, based on a Monte-Carlo code that solves the transfer of photons subject to Compton scattering and pair production/annihilation. Systematic study is performed by imposing various upstream conditions which are characterized by the following three parameters: the photon-to-baryon inertia ratio $\xi_{u *}$, the photon-to-baryon number ratio $\tilde{n}$, and the shock Lorentz factor $\gamma_u$. We find that the properties of RRMSs vary considerably with these parameters. In particular, while a smooth decline in the velocity, accompanied by a gradual temperature increase is seen for $\xi_{u*} \gg 1$, an efficient bulk Comptonization, that leads to a heating precursor, is found for $\xi_{u*} \lesssim 1$. As a consequence, although particle acceleration is highly inefficient in these shocks, a broad non-thermal spectrum is produced in the latter case. The generation of high energy photons through bulk Comptonization leads, in certain cases, to a copious production of pairs that provide the dominant opacity for Compton scattering. We also find that for certain upstream conditions a weak subshock appears within the flow. For a choice of parameters suitable to gamma-ray bursts, the radiation spectrum within the shock is found to be compatible with that of the prompt emission, suggesting that subphotospheric shocks may give rise to the observed non-thermal features despite the absence of accelerated particles.

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.

We examine the effect of Compton drag on the dynamics of a relativistic shock wave. Similarity solutions describing a radiation-supported shock are obtained for certain profiles of the external radiation intensity and the density of the unshocked ejecta, and are compared with 1D numerical simulations of a blast wave expanding into an ambient medium containing isotropic seed radiation. Both the analytic model and the simulations indicate that under realistic conditions the radiation drag should strongly affect the dynamics and structure of the shocked layer. In particular, our calculations show significant strengthening of the reverse shock and weakening of the forward shock over time of the order of the inverse Compton cooling time. We conclude that the effect of radiation drag on the evolution of the emitting plasma should affect the resultant light curves and, conceivably, spectra of the observed emission during strong blazar flares.

Global linear stability analysis of a self-similar solution describing the interaction of a relativistic shell with an ambient medium is performed. The solution is shown to be unstable to convective Rayleigh-Taylor modes having angular scales smaller than the causality scale. Longer wavelength modes are stable and decay with time. For modes of sufficiently large spherical harmonic degree $l$ the dimensionless growth rate scales as $\sqrt{l/\Gamma}$, where $\Gamma$ is the Lorentz factor of the shell. The instability commences at the contact interface separating the shocked eject a and shocked ambient gas and propagates to the shocks. The reverse shock front responds promptly to the in stability and exhibits rapidly growing distortions at early times. Propagation to the forward shock is slower, and it is anticipated that the region near the contact will become fully turbulent before the instability is communicated to the forward shock. The non-universality of the Blandford-McKee blast wave solution suggests that turbulence generated by the instability in the shocked ambient medium may decay slowly with time and may be the origin of magnetic fields over a long portion of the blast wave evolution. It is also speculated that the instability may affect the emission from the shocked ejecta in the early post-prompt phase of GRBs.

Stationary features are occasionally observed in AGN jets. A notable example is the HST-1 knot in M87. Such features are commonly interpreted as re-confinement shocks in hydrodynamic jets or focusing nozzles in Poynting jets. In this paper we compute the structure and Lorentz factor of a highly magnetized jet confined by external pressure having a profile that flattens abruptly at some radius. We find the development of strong oscillations upon transition from the steeper to the flatter pressure profile medium. Analytic formula is derived for the location of the nodes of these oscillations. We apply the model to the M87 jet and show that if the jet remains magnetically dominated up to sub-kpc scales, then focusing is expected. The location of the HST-1 knot can be reconciled with recent measurements of the pressure profile around the Bondi radius if the jet luminosity satisfies $L_j\simeq10^{43}$ erg/s. However, we find that magnetic domination at the collimation break implies a Lorentz factor in excess of $10^2$, atypical to FRI sources. A much lower value of the asymptotic Lorentz factor would require substantial loading close to the black hole. In that case HST-1 may be associated with a collimation nozzle of a hydrodynamic flow.

The collimation of a Poynting-flux dominated jet by a wind emanating from the surface of an accretion flow is computed using a semi-analytic model. The injection of the disk wind is treated as a boundary condition in the equatorial plane, and its evolution is followed by invoking a prescribed geometry of streamlines. Solutions are obtained for a wide range of disk wind parameters. It is found that jet collimation generally occurs when the total wind power exceeds about ten percents of the jet power. For moderate wind powers we find gradual collimation. For strong winds we find rapid collimation followed by focusing of the jet, after which it remains narrow over many Alfv\'en crossing times before becoming conical. We estimate that in the later case the jet's magnetic field may be dissipated by the current-driven kink instability over a distance of a few hundreds gravitational radii. We apply the model to M87 and show that the observed parabolic shape of the radio jet within the Bondi radius can be reproduced provided that the wind injection zone extends to several hundreds gravitational radii, and that its total power is about one third of the jet power. The radio spectrum can be produced by synchrotron radiation of relativistically hot, thermal electrons in the sheath flow surrounding the inner jet.

The dynamics of a dissipative Poynting dominated flow subject to a radiation drag due to Compton scattering of ambient photons by relativistic electrons accelerated in reconnecting current sheets is studied. It is found that the efficiency at which magnetic energy is converted to radiation is limited to a maximum value of $\epsilon_c=3l_{dis}\,\sigma_0/4(\sigma_0+1)$, where $\sigma_0$ is the initial magnetization of the flow and $l_{dis}\le1$ the fraction of initial Poynting flux that can dissipate. The asymptotic Lorentz factor satisfies $\Gamma_\infty\ge\Gamma_0(1+l_{dis}\,\sigma_0/4)$, where $\Gamma_0$ is the initial Lorentz factor. This limit is approached in cases where the cooling time is shorter than the local dissipation time. A somewhat smaller radiative efficiency is expected if radiative losses are dominated by synchrotron and SSC emissions. It is suggested that under certain conditions magnetic field dissipation may occur in two distinct phases: On small scales, asymmetric magnetic fields that are advected into the polar region and dragged out by the outflow dissipate to a more stable configuration. The dissipated energy is released predominantly as gamma rays. On much larger scales, the outflow encounters a flat density profile medium and re-collimates. This leads to further dissipation and wobbling of the jet head by the kink instability, as found recently in 3D MHD simulations. Within the framework of a model proposed recently to explain the dichotomy of radio loud AGN, this scenario can account for the unification of gamma-ray blazars with FRI and FRII radio sources.

Black holes in core-collapse of massive stars are expected to surge in mass and angular momentum by hyper-accretion immediately following their formation. We here describe a general framework of extended emission in gravitational waves from non-axisymmetric accretion flows from fallback matter of the progenitor envelope. It shows (a) a maximum efficiency in conversion of accretion energy into gravitational waves at hyper-accretion rates exceeding a critical value set by the ratio of the quadrupole mass inhomogeneity and viscosity with (b) a peak characteristic strain amplitude at the frequency $f_b=\Omega_b/\pi$, where $\Omega_b$ is the Keplerian angular velocity at which viscous torques equal angular momentum loss in gravitational radiation, with $h_{char}\propto f^{1/6}$ at $f<f_b$ and $h_{char}\propto f^{-1/6}$ at $f>f_b$. Upcoming gravitational wave observations may probe this scaling by extracting broadband spectra using time-sliced matched filtering with chirp templates, recently developed for identifying turbulence in noisy time series.

We compute the dynamics and emission of dissipative shells that are subject to a strong Compton drag, under simplifying assumptions about the dissipation mechanism. We show that under conditions prevailing in blazars, substantial deceleration is anticipated on sub-parsec and parsec scales in cases of rapid dissipation. Such episodes may be the origin of some of the flaring activity occasionally observed in gamma ray blazars. The shape of the light curves thereby produced reflects the geometry of the emitting surface if the deceleration is very rapid, or the dynamics of the shell if the deceleration is delayed, or initially more gradual, owing, e.g., to continuous injection of energy and momentum.

The Burst and Transient Source Experiment (BATSE) classifies cosmological gamma-ray bursts (GRBs) into short (less than 2 s) and long (over 2 s) events, commonly attributed to mergers of compact objects and, respectively, peculiar core-collapse supernovae. This standard classification has recently been challenged by the {\em Swift} discovery of short GRBs showing Extended Emission (SGRBEE) and nearby long GRBs without an accompanying supernova (LGRBN). Both show an initial hard pulse, characteristic of SGRBs, followed by a long duration soft tail. We here consider the spectral peak energy ($E_{p,i}$)-radiated energy $(E_{iso})$ correlation and the redshift distributions to probe the astronomical and physical origin of these different classes of GRBs. We consider {\em Swift} events of $15$ SGRBs, $7$ SGRBEEs, 3 LGRBNs and 230 LGRBs. The spectral-energy properties of the initial pulse of both SGRBEE and LGRBNs are found to coincide with those of SGRBs. A Monte Carlo simulation shows that the redshift distributions of SGRBs, SGRBEE and LGRBNs fall outside the distribution of LGRBs at 4.75$\,\sigma$, 4.67$\,\sigma$ and 4.31$\sigma$, respectively. A distinct origin of SGRBEEs with respect to LGRBs is also supported by the elliptical host galaxies of the SGRBEE events 050509B and 050724. This combined evidence supports the hypothesis that SGRBEE and LGRBNs originate in mergers as SGRBs. Moreover, long/soft tail of SGRB and LGRBNs satisfy the same $E_{p,i}-E_{iso}$ Amati-correlation holding for normal LGRBs. This fact points to rapidly rotating black holes as a common long-lived inner engine produced by different astronomical progenitors (mergers and supernovae).

We construct a semi-analytic model for MHD flows in Kerr geometry, that incorporates energy loading via neutrino annihilation on magnetic field lines threading the horizon. We compute the double-flow structure for a wide range of energy injection rates, and identify the different operation regimes. At low injection rates the outflow is powered by the spinning black hole via the Blandford-Znajek mechanism, whereas at high injection rates it is driven by the pressure of the plasma deposited on magnetic field lines. In the intermediate regime both processes contribute to the outflow formation. The parameter that quantifies the load is the ratio of the net power injected below the stagnation radius and the maximum power that can be extracted magnetically from the black hole.

We compute the time integrated, thermal emission produced by a series of radiation mediated shocks that emerge from the photosphere of a GRB outflow. We show that for a sufficiently broad distribution of shock strengths, the overall shape of the time integrated spectral energy distribution below the peak is a power law, $\nu E_\nu\propto \nu^\alpha$, with a slope $1<\alpha<2$. A substructure in the SED can also be produced in this model for certain choices of the shock train distribution. In particular, we demonstrate that our model can reproduce the double-peak SED observed in some bursts, in events whereby a strong shock is followed by a sequence sufficiently weaker ones.

We discuss the role plasma injection plays in the formation of outflows in Kerr spacetime. Using a model for the double flow established in the polar region of a rotating black hole, we study the interplay between the different processes that can power the outflow. In particular, we find two types of flows with distinct properties that depend on the rate at which energy is deposited in the magnetosphere. We discuss the implications of this result for gamma ray bursts outflows.

The effect of mass and energy loading on the efficiency at which energy can be extracted magnetically from a Kerr black hole is explored, using a semi-analytic, ideal MHD model that incorporates plasma injection on magnetic field lines. We find a critical load below which the specific energy of the plasma inflowing into the black hole is negative, and above which it is positive, and identify two types of flows with distinct properties; at sub-critical loads a magnetic outflow is launched from the ergosphere, owing to extraction of the black hole spin energy, as originally proposed by Blandford and Znajek. At super-critical loads the structure of the flow depends on the details of the injection process. In cases where the injected plasma is relativistically hot, a pressure driven, double trans-magnetosonic flow is launched from a stagnation point located outside the ergosphere, between the inner and outer light cylinders. Some fraction of the energy deposited in the magnetosphere is then absorbed by the black hole and the rest emerges at infinity in the form of a relativistic outflow. When the injected plasma is cold an outflow may not form at all. We discuss the implications of our results to gamma ray bursts and active galactic nuclei.

The structure of a hydrodynamic, double transonic flow driven by neutrino annihilation in the polar region of a Schwarzschild black hole is computed for different energy deposition profiles. The requirement that both, the inflow into the black hole and the outflow to infinity pass smoothly through their sonic points fixes the stagnation radius and stagnation pressure. The asymptotic power of the outflow is shown to be the integral of the energy deposition rate above the stagnation radius. The outflow production efficiency depends on the energy deposition profile, and is generally higher for shallower profiles. Using recent calculations of the neutrino annihilation rate, we estimate that over 50 percents of the total energy deposited above the horizon can emerge in the form of a relativistic outflow at infinity. The continuous creation of plasma during the expansion of the outflow leads to generation of a large specific entropy. This has important implications for the prompt photospheric emission in GRBs.

We study the collimation of a highly magnetized jet by a surrounding cocoon that forms as a result of the interaction of the jet with the external medium. We show that in regions where the jet is well confined by the cocoon, current-driven instabilities should develop over timescales shorter than the expansion time of the jet's head. We speculate that these instabilities would give rise to complete magnetic field destruction, whereby the jet undergoes a transition from high-to-low sigma above the collimation zone. Using this assumption, we construct a self-consistent model for the evolution of the jet-cocoon system in an ambient medium of arbitrary density profile. We apply the model to jet breakout in long GRBs, and show that the jet is highly collimated inside the envelope of the progenitor star, and is likely to remain confined well after breakout. We speculate that this strong confinement may provide a channel for magnetic field conversion in GRB outflows, whereby the hot, low-sigma jet section thereby produced is the source of the photospheric emission observed in many bursts.

A shock that form below the photosphere of a GRB outflow is mediated by Compton scattering of radiation advected into the shock by the upstream fluid. The characteristic scale of such a shock, a few Thomson depths, is larger than any kinetic scale involved by several orders of magnitudes, hence, unlike collisionless shocks, radiation mediated shocks cannot accelerate particles to non-thermal energies. The spectrum emitted by a shock that emerges from the photosphere of a GRB jet, reflects the temperature profile downstream of the shock, with a possible contribution at the highest energies from the shock transition layer itself. We study the properties of radiation mediated shocks that form during the prompt phase of GRBs, and compute the time integrated spectrum emitted by the shocked fluid following shock breakout. We show that the time integrated emission from a single shock exhibits a prominent thermal peak, with the location of the peak depending on the shock velocity profile. We also point out that multiple shock emission can produce a spectrum that mimics a Band spectrum.

Some core-collapse supernovae appear to be hyper-energetic, and a subset of these are aspherical and associated with long GRBs. We use observations of electromagnetic emission from core-collapse supernovae and GRBs to impose constraints on their free energy source as a prior to searches for their gravitational wave emission. We review these events based on a finite efficiency for the conversion of spin energy to magnetic winds powering supernovae. We find that some of the hyper-energetic events cannot be powered by the spindown of rapidly rotating proto-neutron stars by virtue of their limited rotational energy. They can, instead, be produced by the spindown of black holes providing a distinct prospect for gravitational-wave emission of interest to LIGO, Virgo, and the LCGT.

It is proposed that the variable TeV emission observed in M87 may be produced in a starved magnetospheric region, above which the outflow associated with the VLBA jet is established. It is shown that annihilation of MeV photons emitted by the radiative inefficient flow in the vicinity of the black hole, can lead to injection of seed charges on open magnetic field lines, with a density that depends sensitively on accretion rate, $n_\pm\propto\dot{m}^{4}$. For an accretion rate that corresponds to the inferred jet power, and to a fit of the observed SED by an ADAF model, the density of injected pairs is found to be smaller than the Goldreich-Julian density by a factor of a few. It is also shown that inverse Compton scattering of ambient photons by electrons (positrons) accelerating in the gap can lead to a large multiplicity, $\sim 10^3$, while still allowing photons at energies of up to a few TeV to freely escape the system. The estimated gap width is not smaller than $0.01 r_s$ if the density of seed charges is below the Goldreich-Julian value. The VHE power radiated by the gap can easily account for the luminosity of the TeV source detected by H.E.S.S. The strong dependence of injected pair density on accretion rate should render the gap emission highly intermittent. We also discuss briefly the application of this mechanism to Sgr A$^\star$.

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.

A recent progress in the study of $\gamma$-ray jets is reviewed, with a focus on some theoretical interpretations of the VHE emission from M87, and possibly other misaligned blazars; the connection betwe en the GeV breaks exhibited by bright LAT blazars and opacity sources in the broad line region; the consequences of the detection of GeV emission from GRBs to models of magnetic outf lows; and the implications of the thermal emission observed is some GRBs to dissipation of the outflow bulk energy.

The interaction of relativistic magnetized ejecta with an ambient medium is studied for a range of structures and magnetization of the unshocked ejecta. We particularly focus on the effect of the ambient medium on the dynamics of an impulsive, high-sigma shell. It is found that for sufficiently high values of the initial magnetization $\sigma_0$ the evolution of the system is significantly altered by the ambient medium well before the shell reaches its coasting phase. The maximum Lorentz factor of the shell is limited to values well below $\sigma_0$; for a shell of initial energy $E=10^{52}E_{52}$ ergs and size $r_0=10^{12}T_{30}$ cm expelled into a medium having a uniform density $n_i$ we obtain $\Gamma_{\rm max}\simeq180(E_{52}/T_{30}^3 n_i)^{1/8}$ in the high sigma limit. The reverse shock and any internal shocks that might form if the source is fluctuating are shown to be very weak. The restriction on the Lorentz factor is more severe for shells propagating in a stellar wind. Intermittent ejection of small sub-shells doesn't seem to help, as the shells merge while still highly magnetized. Lower sigma shells start decelerating after reaching the coasting phase and spreading away. The properties of the reverse shock then depend on the density profiles of the coasting shell and the ambient medium. For a self-similar cold shell the reverse shock becomes strong as it propagates inwards, and the system eventually approaches the self-similar solution derived recently by Nakamura \& Shigeyama.

A brief overview of jets and their central drivers is presented, with a focus on accreting black hole systems. In particular, scaling relations that elucidate some basic properties of the engines are derived, and the implications for the associated outflows are discussed. The kinematics and dynamics of relativistic jets in various systems and the dissipation of their bulk energy is considered, with an emphasis on consequences of recent observations. Also considered is the interaction of the jets with their environment. Comments on multi-messenger probes are made at the end.

Global linear stability analysis of a self-similar solution describing a relativistic shell decelerated by an ambient medium is performed. The system is shown to be subject to the convective Rayleigh-Taylor instability, with a rapid growth of eigenmodes having angular scale much smaller than the causality scale. The growth rate appears to be largest at the interface separating the shocked ejecta and shocked ambient gas. The disturbances produced at the contact interface propagate in the shocked media and cause nonlinear oscillations of the forward and reverse shock fronts. It is speculated that such oscillations may affect the emission from the shocked ejecta in the early afterglow phase of GRBs, and may be the origin of the magnetic field in the shocked circum-burst medium.

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.

The $\gamma$-ray emission observed in several classes of Galactic and extragalactic astrophysical sources appears to be linked to accreting black holes and rotational powered neutron stars. These systems are prodigious cosmic accelerators, and are also potential sources of the UHE cosmic rays detected by several experiments and VHE neutrinos. We review a recent progress in our understanding of these objects, and demonstrate how recent and future observations can be employed to probe the conditions in the sources.

This is a report on the findings of the Galactic compact objects working group for the white paper on the status and future of TeV gamma-ray astronomy. The white paper is an APS commissioned document, and the full version has also been released and can be found on astro-ph. This section of the white paper discusses the potential of future ground-based gamma-ray experiments to advance our understanding of the physics of Galactic compact objects including pulsars, pulsar wind nebulae, and X-ray binaries.

Rapid variability of the TeV emission in several blazars implies a central black hole mass $M_{BH}<10^8M_\odot$, appreciably smaller than the values estimated from the $M_{BH}-L_{bulg}$ relation, and Doppler factors for the $\gamma$-ray emitting fluid much larger than those associated with radio patterns. We discuss the conditions in the central engine required to account for the short timescales and large luminosities observed, and propose some explanations for the inferred kinematics of the source on various scales.

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.

The rapid variability of the VHE emission reported for some TeV blazars implies Doppler factors well in excess of those inferred from superluminal motions and unification schemes. We propose that those extreme flares may result from radiative deceleration of blobs on scales where local dissipation occurs. The minimum jet power estimated from the resolved synchrotron emission on VLBI scales appears to be consistent with this model. It is shown that if the energy distribution of nonthermal electrons accelerated locally in the blob is reasonably flat, then a background radiation field having a luminosity in the range 10$^{41}-10^{42}$ erg s$^{-1}$ can give rise to a substantial deceleration of the blob, but still be transparent enough to allow the TeV $\gamma$-rays thereby produced to escape the system.

We discuss some phenomenological aspects of $\gamma$-ray emitting jets. In particular, we present calculations of the $\gamma$-sphere and $\pi$-sphere for various target photon fields, and employ them to demonstrate how $\gamma$-ray observations at very high energies can be used to constraint the Doppler factor of the emitting plasma and the production of VHE neutrinos. We also consider the implications of the rapid TeV variability observed in M87 and the TeV blazars, and propose a model 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. Finally, we briefly discuss the prospects for detecting VHE neutrinos from relativistic jets.

We compute the structure and degree of neutronization of general relativistic magnetohydrodynamic (GRMHD) outflows originating from the inner region of neutrino-cooled disks. We consider both, outflows expelled from a hydrostatic disk corona and outflows driven by disk turbulence. We show that in outflows driven thermally from a static disk the electron fraction quickly evolves to its equilibrium value which is dominated by neutrino capture. Those outflows are generally proton rich and, under certain conditions, can be magnetically dominated. They may also provide sites for effective production of 56Ni. Centrifugally driven outflows and outflows driven by disk turbulence, on the other hand, can preserve the large in-disk neutron excess. Those outflows are, quite generally, subrelativistic by virtue of the large mass flux driven by the additional forces.

Constraints on the dynamics, dissipation, and production of VHE neutrinos in relativistic jets are derived using opacity calculations and VHE $\gamma$-ray observations. In particular, it is demonstrated how rapid variability of the $\gamma$-ray emission at very high energies ($> 100$ GeV) can be used to map the location of the $\gamma$-spheres, to derive lower limits on the Doppler factor of the $\gamma$-ray emission zone, and to constrain the photopion production opacity. The apparent discrepancy between jet Lorentz factors inferred from superluminal motions and source statistics in the TeV blazars and those derived from the $\gamma$-ray emission is discussed. The relation to the high-energy emission from the HST1 knot in M87 is briefly mentioned. Estimates of neutrino yields in upcoming neutrino telescopes are given for various sources. It is shown that for TeV blazars, the rapid variability of the TeV emission implies neutrino yields well below detection limit.

Various aspects of the high-energy emission from relativistic jets associated wi th compact astrophysical systems are reviewed. The main leptonic and hadronic processes responsible for the production of high-energy $\gamma$-rays, very-high energy neutrinos and ultra-high energy cosmic rays are discussed. Relations between the $\gamma\gamma$ pair production and photomeson production opacities are derived, and their consequences for the relative emission of $\gamma$-rays and neutrinos are examined. The scaling of the size and location of the various emission zones and other quantities with black hole mass and dimensionless luminosity is elucidated. The results are applied to individual classes of objects, including blazars, microquasars and gamma-ray bursts. It is concluded that if baryons are present in the jet at sufficient quantities, then under optimal conditions most systems exhibiting relativistic jets may be detectable by upcoming neutrino telescopes. An exception is the class of TeV blazars, for which $\gamma$-ray observations imply neutrino yields well below detection limit.

A GRMHD model of disk outflows with neutrino-driven mass ejection is presented. The model is used to calculate the structure of the outflow in the sub-slow magnetosonic region and the mass loading of the outflow, under conditions anticipated in the central engines of gamma-ray bursts. It is concluded that magnetic launching of ultra-relativistic polar outflows is in principle possible along low inclination field lines (with respect to the symmetry axis), provided the neutrino luminosity is sufficiently low, $L_\nu\simlt10^{52}$ erg s$^{-1}$.

The planets capture model for the eruption of V838 Mon is discussed. We used three methods to estimate the location where the planets were consumed. There is a nice consistency for the results of the three different methods, and we find that the typical stopping / slowing radius for the planets is about 1Ro. The three peaks in the optical light curve of V838 Mon are either explained by the swallowing of three planets at different radii or by three steps in the slowing down process of a single planet. We discuss the other models offered for the outburst of V838 Mon, and conclude that the binary merger model and the planet/s scenario seem to be the most promising. These two models have several similarities, and the main differences are the stellar evolutionary stage, and the mass of the accreted material. We show that the energy emitted in the V838 Mon event is consistent with the planets scenario. We suggest a few explanations for the trigger for the outburst and for the double structure of the optical peaks in the light curve of V838 Mon.

It is noted that the Liang-Zhang correlation can be accounted for with the viewing angle interpretation proposed earlier. The Ghirlanda correlation, recently generalized by Nava et al (2006) to a wind profile, can be accounted for by the viewing angle interpretation accordingly generalized to a wind profile. Most of the scatter in the spectra and time-integrated brightness in $\gamma$-ray bursts (GRB) can thus be accounted for by variation in two parameters, 1) the viewing angle and 2) the jet opening angle, with very little variation in any other intrinsic parameters. The scatter in apparent isotropic equivalent fluence and other parameters is reduced by a factor of order 30 when each of these parameters is considered. Possible difficulties with alternative explanations are briefly discussed. It is also noted that the relative scatter in the Amati and Ghirlanda correlations suggests certain conclusions about the inner engine.

V838 Mon is the prototype of a new class of objects. Understanding the nature of its multi-stage outburst and similar systems is challenging. So far, several scenarios have been invoked to explain this group of stars. In this work, the planets-swallowing model for V838 Mon is further investigated, taking into account the findings that the progenitor is most likely a massive B type star. We find that the super-Eddington luminosity during the eruption can explain the fast rising times of the three peaks in the optical light curve. We used two different methods to estimate the location where the planets were consumed. There is a nice agreement between the values obtained from the luminosities of the peaks and from their rising time scale. We estimate that the planets were stopped at a typical distance of one solar radius from the center of the host giant star. The planets-devouring model seems to give a satisfying explanation to the differences in the luminosities and rising times of the three peaks in the optical light curve of V838 Mon. The peaks may be explained by the consumption of three planets or alternatively by three steps in the terminal falling process of a single planet. We argue that only the binary merger and the planets-swallowing models are consistent with the observations of the new type of stars defined by V838 Mon.

(short version) - A model for GRMHD disk outflows with neutrino-driven mass ejection is developed,and employed to calculate the structure of the outflow in the sub-slow magnetosonic region and the mass loading of the outflow, under conditions anticipated in the central engines of gamma-ray bursts. The dependence of the mass flux on the conditions in the disk, on magnetic field geometry, and on other factors is carefully examined for a range of neutrino luminosities expected in hyperaccreting black holes. The fraction of neutrino luminosity that is being converted to kinetic energy flux is shown to be a sensitive function of the effective neutrino temperature at the flow injection point, and the shape of magnetic field lines in the sub-slow region, but is practically independent of the strength of poloidal and toroidal magnetic fields. We conclude that magnetic launching of ultra-relativistic polar outflows from the innermost parts of the disk is in principle possible provided the neutrino luminosity is sufficiently low, $L_\nu\simlt10^{52}$ erg s$^{-1}$ or so. The conditions found to be optimal for the launching of an ultra-relativistic jet are also the conditions favorable for large neutron-to-proton ratio in the disk.

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

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.

It is widely believed that some classes of high-energy transients may be powered by the rotational energy of a rapidly spinning black hole. The energy extraction mechanism commonly discussed involves macroscopic magnetic fields that are produced by currents flowing in a disk or torus surrounding the black hole. In spite of large efforts that have led to impressive progress in our understanding of the physics of magnetized flows in Kerr geometry, several important issues remain unresolved: the nature of the load in the Blandford-Znajeck model, the causality and stability of a force-free magnetosphere, global current closure, and the role of boundary conditions on the channels through which the extracted energy is released. This chapter provides an account of recent work that addresses some of those open questions.

It is pointed out that the beaming correction commonly inferred from the achromatic breaks in the afterglow light curve is biased in situations where the isotropic equivalent energy is affected by factors other than the spread in opening angles. In particular, it underestimates the beaming factor of sources observed off-axis. Here we show that both the slopes and scatters in the Amati and Ghirlanda relations, and the difference between them, are quantitatively consistent with a model proposed recently, in which the $E_{iso}-\nu_{peak}$ relation, as originally derived by Amati et al., is due to viewing angle effects. The quantitative difference between them confirms the relations between opening angle and break time suggested by Frail et al.

A time-dependent model for pair creation in a pulsar magnetosphere is developed. It is argued that the parallel electric field that develops in a charge-starved region (a gap) of a pulsar magnetosphere oscillates with large amplitude. Electrons and positrons are accelerated periodically and the amplitude of the oscillations is assumed large enough to cause creation of upgoing and downgoing pairs at different phases of the oscillation. With a charge-starved initial condition, we find that the oscillations result in bursts of pair creation in which the pair density rises exponentially with time. The pair density saturates at $N_\pm\simeq E_{0}^2/(8\pi m_ec^2\Gamma_{\rm thr})$, where $E_0$ is the parallel electric field in the charge-starved initial state, and $\Gamma_{\rm thr}$ is the Lorentz factor for effec tive pair creation. The frequency of oscillations following the pair creation burst is given roughly by $\omega_{\rm osc}=eE_0/(8m_ec\Gamma_{\rm thr})$. A positive feedback keeps the system stable, such that the average pair creation rate balances the loss rate due to pairs escaping the magnetosphere.

The $h\nu_{peak}/100 KeV \sim [E_{iso}/10^{52} \rm{erg}]^{1/2}$ correlation reported recently by Amati et al. (2002) and Lamb et al. (2003) for long GRB's can be interpreted as a viewing angle effect if the emitting region of the fireball is ring-shaped.

The polarization of gamma ray emission scattered off the baryon rich material that collimates a GRB fireball and the coasting speed of the irradiated matter are calculated numerically for different geometries of the radiation source and the collimating wall. It is shown that when the scattering material is Compton sailing, the direction of maximum polarization is quite generally well within the beaming cone of scattered radiation. As a result, the probability for observing bright highly polarized GRB's increases considerably, provided the Lorentz factor of the coasting matter is not well below 30, and the scattered radiation is highly polarized even when the beam thickness is large compared to $1/\Gamma$. It is suggested that correlation between polarization and intensity could provide clues as to whether energy flows from matter to photons or the reverse.

The rotational energy of a Kerr black hole is often invoked as the free energy source that powers compact astrophysical systems An important question concerning the energy extraction mechanism is whether a Kerr black hole can communicate stresses to a distant load via a surrounding, force-free magnetosphere, and whether such structures are stable. In this paper we address this question. We first re-examine the properties of short wavelength force-free waves, and show that contrary to earlier claims, fast magnetosonic disturbances do affect the global electric current flowing in the system, even in flat spacetime; beyond the light cylinder the perturbed charge density is on the order of the perturbed Goldreich-Julian charge density. This implies that the fast magnetosonic surface is a causal wind boundary. We then go on to study the evolution of a force-free system driven by a spinning-up black hole, by solving the time dependent Maxwell equations near the horizon using perturbation method. We find that the electromagnetic field, current and charge density grow exponentially in the linear regime. We conclude that after time $t=2M\ln(1/2\alpha^2)$, where $\alpha$ is the lapse function, the evolution will reach an adiabatic stage, whereby the solution is to a good approximation the static one, with the angular momentum of the hole being the adiabatic parameter. This implies that black hole magnetospheres are stable. A brief discussion on the relation between global current closure and causality in force-free systems is given at the end.

Gamma-ray bursts are believed to originate in core-collapse of massive stars. This produces an active nucleus containing a rapidly rotating Kerr black hole surrounded by a uniformly magnetized torus represented by two counter-oriented current rings. We quantify black hole spin-interactions with the torus and charged particles along open magnetic flux-tubes subtended by the event horizon. A major output of Egw=4e53 erg is radiated in gravitational waves of frequency fgw=500 Hz by a quadrupole mass-moment in the torus. Consistent with GRB-SNe, we find (i) Ts=90s (tens of s, Kouveliotou et al. 1993), (ii) aspherical SNe of kinetic energy Esn=2e51 erg (2e51 erg in SN1998bw, Hoeflich et al. 1999) and (iii) GRB-energies Egamma=2e50 erg (3e50erg in Frail et al. 2001). GRB-SNe occur perhaps about once a year within D=100Mpc. Correlating LIGO/Virgo detectors enables searches for nearby events and their spectral closure density 6e-9 around 250Hz in the stochastic background radiation in gravitational waves. At current sensitivity, LIGO-Hanford may place an upper bound around 150MSolar in GRB030329. Detection of Egw thus provides a method for identifying Kerr black holes by calorimetry.

GRB fireballs could be ensheathed by a relativistic, baryon-enriched outflow that would in selected directions scatter a polarized beam of gamma rays. The strong polarization recently announced by RHESSI is attainable by the scattering mechanism, though not most of the time. The scattering mechanism, when acting in addition to the synchrotron mechanism, could circumvent some of the predictions implied by a purely synchrotron origin for the polarization. The time profile of the scattered component could differ from that of the primary gamma rays.

It is proposed that the baryons in gamma ray burst (GRB) fireballs originate as "pick-up" neutrons that leak in sideways from surrounding baryonic matter and convert to protons in a collision avalanche. The asymptotic Lorentz factor is estimated, and, in the absence of collimation, is shown to be angle dependent. Reasonable agreement is obtained with existing limits on the GRB baryonic component. The charged decay and collision products of the neutrons become ultrarelativistic immediately, and a UHE neutrino burst is produced with an efficiency that can exceed 0.5. Other signatures may include lithium, beryllium and/or boron lines in the supernova remnants associated with GRB's and high polarization of the gamma rays.

(Abbrev.) We analyze the topology, lifetime, and emissions of a torus around a black hole formed in hypernovae and black hole-neutron star coalescence. The torus is ab initio uniformly magnetized, represented by two counter oriented current-rings, and develops a state of suspended accretion against a "magnetic wall" around the black hole. Magnetic stability of the torus gives rise to a new fundamental limit EB/Ek<0.1 for the ratio of poloidal magnetic field energy-to-kinetic energy. The lifetime of rapid spin of the black hole is effectively defined by the timescale of dissipation of black hole-spin energy in the horizon, and satisfies T= 40s (MH/7MSun)(R/6MH)^4(0.03MH/MT) for a black hole of mass MH surrounded by a torus of mass MT and radius R. The torus converts a major fraction Egw/Erot=0.1 into gravitational radiation through a finite number of multipole mass-moments, and a smaller fraction into MeV neutrinos and baryon-rich winds. At a source distance of 100Mpc, these emissions over N=2e4 periods give rise to a characteristic strain amplitude \sqrt{N}hchar=6e-21. We argue that torus winds create an open magnetic flux-tube on the black hole, which carries a minor and standard fraction Ej/Erot=1e-3 in baryon-poor outflows to infinity. We identify this baryon poor output of tens of seconds with GRBs with contemporaneous and strongly correlated emissions in gravitational radiation, conceivably at multiple frequencies. Ultimately, this leaves a black hole binary surrounded by a supernova remnant.

We discuss the possibility that microquasar jets may be powerful emitters of TeV neutrinos. We estimate the neutrino fluxes produced by photopion production in the jets of a sample of identified microquasars and microquasar candidates, for which available data enables rough determination of the jet parameters. We demonstrate that in several of the sources considered, the neutrino flux at Earth, produced in events similar to those observed, can exceed the detection threshold of a km^2 neutrino detector. Sources with bulk Lorentz factors larger than those characteristic of the sample considered here, directed along our line of sight may be very difficult to resolve at radio wavelengths and hence may be difficult to identify as microqusar candidates. However these sources can be identified through their neutrino and gamma-ray emission.

The number of orphan radio afterglows associated with gamma-ray bursts (GRBs) that should be detected by a flux limited radio survey, is calculated. It is shown that for jetted GRBs this number is smaller for smaller jet opening angle (theta), contrary to naive expectation. For a beaming factor f_b^{-1}=(theta^2/2)^{-1} = 500, roughly the value inferred by Frail et al. (2001) from analysis of afterglow light curves, we predict that between several hundreds to several thousands orphan radio afterglows should be detectable (over all sky) above 1 mJy at GHz frequencies at any given time. This orphan population is dominated by sources lying at distances of a few hundred Mpc, and having an age of ~1 yr. A search for point-like radio transients with flux densities greater than 6 mJy was conducted using the FIRST and NVSS surveys, yielding a list of 25 orphan candidates. We argue that most of the candidates are unlikely to be radio supernovae. However, the possibility that they are radio loud AGNs cannot be ruled out without further observations. Our analysis sets an upper limit for the all sky number of radio orphans, which corresponds to a lower limit f_b^{-1}>10 on the beaming factor. Rejection of all candidates found in our search would imply f_b^{-1}>100. This, and the fact that some candidates may indeed be radio afterglows, strongly motivate further observations of these transients.

It has been proposed recently that Galactic microquasars may be prodigious emitters of TeV neutrinos that can be detected by upcoming km^2 neutrino telescopes. In this paper we consider a sample of identified microquasars and microquasar candiates, for which available data enables rough determination of the jet parameters. By employing the parameters inferred from radio observations of various jet ejection events, we determine the neutrino fluxes that should have been produced during these events by photopion production in the jet. Despite the large uncertainties in our analysis, we demonstrate that in several of the sources considered, the neutrino flux at Earth, produced in events similar to those observed, would exceed the detection threshold of a km^2 neutrino detector. The class of microquasars may contain also sources with bulk Lorentz factors larger than those characteristic of the sample considered here, directed along our line of sight. Such sources, which may be very difficult to resolve at radio wavelengths and hence may be difficult to identify as microqusar candidates, may emit neutrinos with fluxes significantly larger than typically obtained in the present analysis. These sources may eventually be identified through their neutrino and gamma-ray emission.

The jets associated with Galactic micro-quasars are believed to be ejected by accreting stellar mass black-holes or neutron stars. We show that if the energy content of the jets in the transient sources is dominated by electron-proton plasma, then a several hour outburst of 1--100 TeV neutrinos produced by photo- meson interactions should precede the radio flares associated with major ejection events. Several neutrinos may be detected during a single outburst by a 1km^2 detector, thereby providing a powerful probe of micro-quasars jet physics.

Black holes surrounded by a disk or torus may drive the enigmatic cosmological gamma-ray bursts (GRBs). Equivalence in poloidal topology to pulsar magnetospheres shows a high incidence of the black hole-luminosity $L_H$ into the surrounding magnetized matter. We argue that this emission is re-radiated into gravitational waves at $L_{GW}\simeq L_H/3$ in frequencies of order 1kHz, winds and, potentially, MeV neutrinos. The total energy budget and input to the GRB from baryon poor jets are expected to be standard in this scenario, consistent with recent analysis of afterglow data. Collimation of these outflows by baryon rich disk or torus winds may account for the observed spread in opening angles up to about $35^o$. This model may be tested by future LIGO/VIRGO observations.

The possibility that the peak of the spectral energy distribution (SED) of prompt GRB emission is produced on compact scales, prior to the acceleration of the fireball to its terminal Lorentz factor is considered. It is argued that if dissipation on these scales is associated with pair creation processes, then the observed clustering of SED peaks can be accounted for quite naturally by this mechanism. Further dissipation and processing of pairs and gamma rays at larger radii, can give rise to a nonthermal extension of the spectrum up to a rather high observed energy. If the nonthermal processes on these scales conserve the total flux of fireball quanta, then the SED peak, which is established on much smaller scales, will not be altered significantly. The asymptotic bulk Lorentz factor in this model should exceed several hundred.

The possibility that the excess cosmic ray (CR) flux near 10$^{18}$ eV, reported recently by the AGASA group, is due to a compact black hole dynamo associated with the Sgr A* source is considered. Under the assumption that the Galactic center black hole rotates with nearly maximal spin, and that the magnetic field threading the horizon is in rough equipartition with matter accreted by the hole, the spectra and total fluxes of accelerated CRs and their associated curvature emission depend only on accretion rate. For the accretion rate estimated on the basis of observations of stellar winds near Sgr A* ($\sim 10^{-3}$ Eddington), the maximum proton energy achievable by this mechanism is of the order of $10^{18}$ eV, and the corresponding CR power is $\sim 10^{41}\alpha_{CR}$ erg s$^{-1}$, where $\alpha_{CR}<<1$ is the CR production efficiency. The corresponding spectrum of curvature photons peaks at around a hundred GeV, with a total luminosity comparable to the CR power released. For much lower accretion rates both the CR and gamma-ray fluxes predicted are completely insignificant. Upcoming gamma-ray experiments, such as CANGAROO, HESS and GLAST, can be used to probe the parameters of this system.

The possibility that rapidly rotating supermassive black holes associated with quasar remnants may provide sites for the acceleration of ultra-high energy cosmic rays (UHECRs) is considered. It is shown that prodigious TeV emission through curvature losses is an important consequence of this mechanism. Given the measured UHECR flux, it is expected that nearby dormant AGNs will be detectable by current Tev experiments. The implications for the Sgr A* source are also briefly considered.

The idea that the TeV air showers, thought to be produced by >10 TeV gamma rays from Mrk 501, can be mimicked by coherent bunches of sub-TeV photons is reexamined, focusing on fundamental considerations. In particular, it is shown that the minimum spot size of the beam of pulsed TeV photons arriving at Earth is on the order of a few kilometers, unless a lens with certain characteristics is placed between the TeV laser and Earth. The viability of the laser production mechanism proposed by Harwit et al. (2000) is also reassessed.

The curvature radiation produced by particles accelerating near the event horizon of a spinning supermassive black hole, threaded by externally supported magnetic field lines is considered. It is shown that light nuclei suffer catastrophic curvature losses that limit the maximum energy they can attain to values well below that imposed by the maximum potential difference induced by the black hole dynamo, unless the curvature radius of magnetic field lines largely exceeds the gravitational radius. It is further shown that the dominant fraction of the rotational energy extracted from the black hole is radiated in the TeV band. Given the observed flux of ultra-high energy cosmic rays, and the estimated number of nearby supermassive black holes, it is expected that if dormant AGNs are the sources of UHECRs, as proposed recently by Boldt & Ghosh, then they should also be detectable at TeV energies by present TeV experiments.

Analytic solutions are presented for the hydrodynamic collimation of a relativistic fireball by a surrounding baryonic wind emanating from a torus. The opening angle is shown to be the ratio of the power output of the inner fireball to that of the exterior baryonic wind. The gamma ray burst 990123 might thus be interpreted as a baryon-pure jet with an energy output of order 10$^{50}$ erg or less, collimated by a baryonic wind from a torus with an energy output of order $10^{52}$ erg.

Internal shocks propagating through an ambient radiation field, are subject to a radiative drag that, under certain conditions, can significantly affect their dynamics and, consequently, the evolution of the beaming cone of emission produced behind the shocks. The resultant change of the Doppler factor combined with opacity effects leads to a strong dependence of the variability pattern produced by such systems, specifically, the shape of the light curves and the characteristics of correlated emission, on viewing angle. One implication is that objects oriented at relatively large viewing angles to the observer should exhibit a higher level of activity at high synchrotron frequencies (above the self-absorption frequency) and at gamma-ray energies below the threshold energy to pair production, than at lower (radio/millimeter) frequencies.

It is proposed that the gamma ray burst photons near the peak of the spectrum at several hundred KeV are produced on very compact scales, where photon production is limited by blackbody effects and/or the requirement of energetic quanta ($E>2m_e c^2$) for efficient further production. The fast variation of order milliseconds in the time profile is then a natural expectation, given the other observed GRB parameters. Analytic calculations are presented to show that the escape of non-thermal, energetic gamma rays can emerge within a second of the thermal photons from a gammasphere of below $10^{12}$ cm. The minimum asymptotic bulk Lorentz factor in this model is found to be of order several hundred if the photosphere is of order $3 \times 10^{11}$ cm and greater for larger or smaller photospheric radii. It is suggested that prompt UHE gamma rays might provide a new constraint on the asymptotic Lorentz factor of the outflow.

The time dependent synchrotron emission from relativistic jets, and the relation between the synchrotron and ERC emission is considered within the framework of the radiative front model. The timescale and profile of the optically thin emission are shown to be determined, in this model, by the shock formation radius, the thickness of expelled fluid slab and the variation of the front's parameters due to its transverse expansion. For a range of reasonable conditions, a variety of flare shapes can be produced, varying from roughly symmetric with exponential rises and decays, as often seen in blazars, to highly asymmetric with a fast rise and a much slower, power law decay, as seen in GRB afterglows. The onset, duration, and fluence of low-frequency (below the initial turnover frequency) and hard gamma-ray (above the initial gamma-spheric energy) outbursts are limited by opacity effects; the emission at these energies is quite generally delayed and, in the case of sufficiently short length outbursts, severely attenuated. The observational consequences are discussed. One distinctive prediction of this model is that in a single, powerful source, the upper cutoff of the gamma-ray spectrum should be correlated with the timescale of the outburst and with the amplitude of variations at long wavelengths (typically radio to millimeter).

The gamma ray burst (GRB) 980425 is distinctive in that it seems to be associated with supernova (SN) 1998bw, has no X-ray afterglow, and has a single peak light curve and a soft spectrum. The supernova is itself unusual in that its expansion velocity exceeds c/6. We suggest that many of these features can be accounted for with the hypothesis that we observe the GRB along a penumbral line of sight that contains mainly photons that have scattered off ejected baryons. The hypothesis suggests a baryon poor jet (BPJ) existing within a baryon rich outflow. The sharp distinction can be attributed to whether or not the magnetic field lines thread an event horizon. Such a configuration suggests that there will be some non-thermal acceleration of pick-up ex-neutrons within the BPJ. This scenario might produce observable spallation products and neutrinos.

A large fraction of the AGNs detected by EGRET exhibit strong variability in many wavebands- from radio to the highest energy gamma-rays observed. The origin of the transient emission in blazars is presumably associated with relativistic jets. In this talk plausible variability mechanisms that may operate in blazars are discussed.

The transient emission produced behind internal shocks that are driven by overtaking collisions of a magnetized, relativistic outflow is considered. A self-consistent model capable of describing the structure and dynamics of the shocks and the time evolution of the pair-and gamma-ray distribution functions is developed and applied to gamma-ray flares in blazars, in the case in which gamma-ray production is dominated by inverse Compton scattering of external radiation (ERC). The dependence of the flare properties on magnetic field dissipation rate, intensity of ambient radiation, and the thickness of expelled fluid slabs is analyzed. It is shown that i) the type of gamma-ray flare produced by the model is determined by the ratio of the thickness of ejected fluid slab and the gradient length scale of ambient radiation intensity, ii) the radiative efficiency depends sensitively on the opacity contributed by the background radiation, owing to a radiative feedback, and is typically very high for parameters characteristic to the powerful blazars, and iii) the emitted flux is strongly suppressed at energies for which the pair-production optical depth is initially larger than unity; the time lag and flare duration in this energy range increase with increasing gamma-ray energy. At lower energies, flaring at different gamma-ray bands occurs roughly simultaneously, but with possibly different amplitudes. Some observational consequences are discussed.

The theory of gamma-ray emission from e$^{\pm}$ jets and the implications for jet formation, dynamics and structure are reviewed. In particular, possible carriers of the jet's thrust on small scales, the transition from electromagnetic to particle dominance in Poynting flux jets, formation of pair cascades, synchrotron emission by cascading pairs, and formation of shocks due to unsteadiness in the jet parameters are considered, with emphasis on the observational consequences. Some recent progress in modeling transient emission from blazars is also briefly discussed.

We study molecular hydrogen formation in and on solids. We construct a model with surface sites and bulk sites capable of describing (1) the motion and exchange of H and H_2 between surface and bulk, (2) the recombination of H and dissociation of H_2 in and on the solid, and (3) the injection of H from the gas phase and the loss of H and H_2 from the solid. The basic physical processes include thermally activated reactions, collisionally induced reactions, and tunneling reactions. Our main application is to the astronomical problem of H_2 formation on grains in space but the model has more general applicability. We investigate the steady-state H and H_2 concentrations in and on the solid when gas phase atoms or ions stick to the surface or penetrate the body of the grain. The model identifies ranges of physical parameters for which the solid becomes saturated (surface and/or bulk) with particles (H and/or H_2) and facilitates the calculation of the efficiency of molecule formation (the fraction of the gas phase atoms that leave as molecules). These solutions are highly degenerate in the sense that they depend only on a small number of dimensionless parameters. We find that a variety of recombination pathways operate under a broad range of conditions. As an example we study H_2 formation in and on carbon grains.

We study the formation, evolution and structure of dissipative fronts produced by overtaking collisions of relativistic streams, with emphasis on strongly magnetized flows. The evolution of the system is followed using analytical approach in the simple wave regime and numerical simulations in the non-simple wave regime, until a steady-state is reached. The steady-state structure of the front is then examined by solving the appropriate jump conditions. The conversion of magnetic energy into kinetic energy is parametrized in terms of the Alfv\'en 4-velocity inside the front. The implications for gamma-ray jets are briefly discussed.

Recent observations of GRS 1915+105 and GRO J1655+40 reveal superluminal motions in Galactic sources. This letter examines the physical conditions within these Galactic sources, their interaction with their environment, their possible formation, and contrasts them with their extragalactic counterparts. In particular, e$^{+}$-e$^{-}$ and e-p jets are contrasted, constraints on particle acceleration in the jets are imposed using X-ray and radio observations, the $\gamma$-ray flux from e$^+$-e$^-$ jets expected at EGRET energies and the flux in infrared lines from an e-p jet are estimated. It is also suggested that these sources may exhibit low frequency radio lobes extending up to several hundred parsecs in size, strong, soft X-ray absorption during the birth of the radio components and emission line strengths anti-correlated with the X-ray flux. The implications for other X-ray transients are briefly discussed.