CHIME has now detected 18 repeating fast radio bursts (FRBs). We explore what can be learned about the energy distribution and activity level of the repeaters by constructing a realistic FRB population model, taking into account wait-time clustering and cosmological effects. For a power-law energy distribution dN/dE ~ E^{-gamma} for the repeating bursts, a steep energy distribution means that most repeaters should be found in the local Universe with low dispersion measure (DM), whereas a shallower distribution means some repeaters may be detected at large distances with high DM. It is especially interesting that there are two high-DM repeaters (FRB 181017 and 190417) with DM ~ 1000 pc/cm^3. These can be understood if: (i) the energy distribution is shallow gamma = 1.7 + 0.3 - 0.1 (68% confidence) or (ii) a small fraction of sources are extremely active. In the second scenario, these high-DM sources should be repeating more than 100 times more frequently than FRB 121102, and the energy index is constrained to be gamma = 1.9 + 0.3 - 0.2 (68% confidence). In either case, this power-law index is consistent with the energy dependence of the non-repeating ASKAP sample, which suggests that they are drawn from the same population. Finally, we show that the CHIME repeating fraction can be used to infer the distribution of activity level in the whole population.

Non-thermal particles and high-energy radiation can play a role in the dynamical processes in star-forming regions and provide an important piece of the multiwavelength observational picture of their structure and components. Powerful stellar winds and supernovae in compact clusters of massive stars and OB associations are known to be favourable sites of high-energy particle acceleration and sources of non-thermal radiation and neutrinos. Namely, young massive stellar clusters are likely sources of the PeV (petaelectronvolt) regime cosmic rays (CRs). They can also be responsible for the cosmic ray composition, e.g., 22Ne/20Ne anomalous isotopic ratio in CRs. Efficient particle acceleration can be accompanied by super-adiabatic amplification of the fluctuating magnetic fields in the systems converting a part of kinetic power of the winds and supernovae into the magnetic energy through the CR-driven instabilities. The escape and CR propagation in the vicinity of the sources are affected by the non-linear CR feedback. These effects are expected to be important in starburst galaxies, which produce high-energy neutrinos and gamma-rays. We give a brief review of the theoretical models and observational data on high-energy particle acceleration and their radiation in star-forming regions with young stellar population.

The finite time, $\tau_{\rm dep}$, over which positrons from $\beta^{+}$ decays of $^{56}$Co deposit energy in type Ia supernovae ejecta lead, in case the positrons are trapped, to a slower decay of the bolometric luminosity compared to an exponential decline. Significant light-curve flattening is obtained when the ejecta density drops below the value for which $\tau_{\rm dep}$ equals the $^{56}$Co life-time. We provide a simple method to accurately describe this "delayed deposition" effect, which is straightforward to use for analysis of observed light curves. We find that the ejecta heating is dominated by delayed deposition typically from 600 to 1200~day, and only later by longer lived isotopes $^{57}$Co and $^{55}$Fe decay (assuming solar abundance). For the relatively narrow $^{56}$Ni velocity distributions of commonly studied explosion models, the modification of the light curve depends mainly on the $^{56}$Ni mass-weighted average density, $\langle \rho \rangle t^{3}$. Accurate late-time bolometric light curves, which may be obtained with JWST far-infrared (far-IR) measurements, will thus enable to discriminate between explosion models by determining $\langle \rho \rangle t^3$ (and the $^{57}$Co and $^{55}$Fe abundances). The flattening of light curves inferred from recent observations, which is uncertain due to the lack of far-IR data, is readily explained by delayed deposition in models with $\langle \rho\rangle t^{3} \approx 0.2\,M_{\odot}\,(10^{4}\, \textrm{km}\,\textrm{s}^{-1})^{-3}$, and does not imply supersolar $^{57}$Co and $^{55}$Fe abundances.

High cadence transient surveys are able to capture supernovae closer to their first light than before. Applying analytical models to such early emission, we can constrain the progenitor stars properties. In this paper, we present observations of SN2018fif (ZTF18abokyfk). The supernova was discovered close to first light and monitored by the Zwicky Transient Facility (ZTF) and the Neil Gehrels Swift Observatory. Early spectroscopic observations suggest that the progenitor of SN2018fif was surrounded by relatively small amounts of circumstellar material (CSM) compared to all previous cases. This particularity, coupled with the high cadence multiple-band coverage, makes it a good candidate to investigate using shock-cooling models. We employ the SOPRANOS code, an implementation of the model by Sapir & Waxman (2017). Compared with previous implementations, SOPRANOS has the advantage of including a careful account of the limited temporal validity domain of the shock-cooling model. We find that the progenitor of SN2018fif was a large red supergiant, with a radius of R=1174_{-81}^{+208} solar radii, and an ejected mass of M=5.6_{-1.0}^{+9.1} solar masses. Our model also gives information on the explosion epoch, the progenitor inner structure, the shock velocity and the extinction. The large radius differs from previously modeled objects, and the difference could be either intrinsic or due to the relatively small amount of CSM around SN2018fif, perhaps making it a "cleaner" candidate for applying shock-cooling analytical models.

Dewitt et al. (1973) derived a useful relation between the plasma screening factor for a reaction of two fusing ions and their chemical potentials, based on the plasma pair distribution functions. We show that their result can be derived in a simpler, more straightforward way, by applying the principle of detailed balance, which also enables us to generalize the relation to reactions involving $N$ fusing ions. In order to demonstrate the usefulness of applying the principle of detailed balance, we calculate the screening factor for the $pep$ reaction, $p+e+p\rightarrow ^{2}$D $+\nu_{e}$. For the plasma conditions near the centre of the Sun, the reaction is suppressed by roughly the same amount by which the reaction $p+p\rightarrow ^{2}$D $+$ $e^{+}+\nu_{e}$ is enhanced. This effect may be measured in the near future.

We discuss the late time (tens of days) emission from the radioactive ejecta of mergers involving neutron stars, when the ionization energy loss time of beta-decay electrons and positrons exceeds the expansion time. We show that if the e$^\pm$ are confined to the plasma (by magnetic fields), then the time dependence of the plasma heating rate, $\dot{\varepsilon}_d$, and hence of the bolometric luminosity $L=\dot{\varepsilon}_d$, are given by $d\log L/d\log t\simeq-2.8$, nearly independent of the composition and of the instantaneous radioactive energy release rate, $\dot{\varepsilon}$. This universality of the late time behavior is due to the weak dependence of the ionization loss rate on composition and on e$^\pm$ energy. The late time IR and optical measurements of GW 170817 are consistent with this expected behavior provided that the ionization loss time exceeds the expansion time at $t>t_\varepsilon\approx 7$~d, as predicted based on the early (few day) electromagnetic emission.

A sufficient condition for maser instability in a weakly magnetized relativistic plasma with an isotropic particle distribution function is given. The maser growth rates and polarizations are computed starting from the exact dielectric permittivity tensor of a magnetized plasma. For very weak magnetic fields, our results confirm the approximate validity of the 'standard maser theory', which is based on the Einstein coefficients method, with one significant exception. For inclined propagation and realistic (small but finite) field, the growth rates of the two (nearly circular) polarizations differ significantly, while the standard theory predicts two (nearly circular) polarizations with similar growth rates. We show that this deviation is due to circularly polarized synchrotron emission, which is neglected in the standard theory. The maser is shown to grow slower than Langmuir waves. Nevertheless, significant generation of EM waves is seen in (highly simplified) direct numerical simulations. We study the nonlinear saturation of the maser instability and find that it offers a mechanism for the conversion of a significant fraction of the plasma energy into radio waves. We briefly discuss the conditions under which the maser instability may operate in astrophysical sources, and provide rough estimates that may be used as a guidance when studying particular astrophysical sources/phenomena.

We present a simple analytic model, that captures the key features of the emission of radiation from material ejected by the merger of neutron stars (NS), and construct the multi-band and bolometric luminosity light curves of the transient associated with GW170817, AT\,2017gfo, using all available data. The UV to IR emission is shown to be consistent with a single $\approx0.05$\,M$_\odot$ component ejecta, with a power-law velocity distribution between $\approx 0.1\,c$ and $>0.3\,c$, a low opacity, {$\kappa<1$\,cm$^2$\,g$^{-1}$}, and a radioactive energy release rate consistent with an initial $Y_{\rm e}<0.4$. The late time spectra require an opacity of $\kappa_\nu\approx0.1$\,cm$^2$\,g$^{-1}$ at 1 to $2\mu$m. If this opacity is provided entirely by Lanthanides, their implied mass fraction is $X_{\rm Ln}\approx10^{-3}$, approximately 30 times below the value required to account for the solar abundance. The inferred value of $X_{\rm Ln}$ is uncertain due to uncertainties in the estimates of IR opacities of heavy elements, which also do not allow the exclusion of a significant contribution to the opacity by other elements (the existence of a slower ejecta rich in Lanthanides, that does not contribute significantly to the luminosity, can also not be ruled out). The existence of a relatively massive, $\approx 0.05$\,M$_\odot$, ejecta with high velocity and low opacity is in tension with the results of numerical simulations of NS mergers.

Galaxy clusters, the largest gravitationally bound objects in the Universe, are thought to grow by accreting mass from their surroundings through large-scale virial shocks. Due to electron acceleration in such a shock, it should appear as a $\gamma$-ray, hard X-ray, and radio ring, elongated towards the large-scale filaments feeding the cluster, coincident with a cutoff in the thermal Sunyaev-Zel'dovich (SZ) signal. However, no such signature was found until now, and the very existence of cluster virial shocks has remained a theory. We find preliminary evidence for a large, $\sim 5$ Mpc minor axis $\gamma$-ray ring around the Coma cluster, elongated towards the large scale filament connecting Coma and Abell 1367, detected at the nominal $2.7\sigma$ confidence level ($5.1\sigma$ using control signal simulations). The $\gamma$-ray ring correlates both with a synchrotron signal and with the SZ cutoff, but not with Galactic tracers. The $\gamma$-ray and radio signatures agree with analytic and numerical predictions, if the shock deposits $\sim 1\%$ of the thermal energy in relativistic electrons over a Hubble time, and $\sim 1\%$ in magnetic fields. The implied inverse-Compton and synchrotron cumulative emission from similar shocks can significantly contribute to the diffuse extragalactic $\gamma$-ray and low frequency radio backgrounds. Our results, if confirmed, reveal the prolate structure of the hot gas in Coma, the feeding pattern of the cluster, and properties of the surrounding large scale voids and filaments. The anticipated detection of such shocks around other clusters would provide a powerful new cosmological probe.

In recent years, space-born experiments have delivered new measurements of high energy cosmic-ray (CR) $\bar p$ and $e^+$. In addition, unprecedented sensitivity to CR composite anti-nuclei anti-d and anti-He is expected to be achieved in the near future. We report on the theoretical interpretation of these measurements. While CR antimatter is a promising discovery tool for new physics or exotic astrophysical phenomena, an irreducible background arises from secondary production by primary CR collisions with interstellar matter. Understanding this irreducible background or constraining it from first principles is an interesting challenge. We review the attempt to obtain such understanding and apply it to CR $\bar p,\, e^+,$ anti-d and anti-He. Based on state of the art Galactic cosmic ray measurements, dominated currently by the AMS-02 experiment, we show that: (i) CR $\bar p$ most likely come from CR-gas collisions; (ii) $e^+$ data is consistent with, and suggestive of the same secondary astrophysical production mechanism responsible for $\bar p$ and dominated by proton-proton collisions. In addition, based on recent accelerator analyses we show that the flux of secondary high energy anti-He may be observable with a few years exposure of AMS-02. We highlight key open questions, as well as the role played by recent and upcoming space and accelerator data in clarifying the origins of CR antimatter.

We derive stringent constraints on the persistent source associated with FRB 121102: Size $10^{17}$ cm $<R<10^{18}$ cm, age $<300$ yr, characteristic electron energy $\varepsilon_e\sim0.3$ GeV, total energy $\sim10^{49}$ erg. The hot radiating plasma is confined by a cold plasma of mass $M_c<0.03 (R/10^{17.5}{\rm cm})^4 M_\odot$. The source is nearly resolved, and may be resolved by 10 GHz observations. The fact that $\varepsilon_e\sim m_p c^2$ suggests that the hot plasma was created by the ejection of a mildly relativistic, $M\sim10^{-5} M_\odot$ shell, which propagated into an extended ambient medium or collided with a pre-ejected shell of mass $M_c$. The inferred plasma properties are inconsistent with typical "magnetar wind nebulae" model predictions. We suggest a physical mechanism for the generation of FRBs (independent of the persistent source model): Ejection from an underlying compact object, $R_s\sim10^{6}$ cm, of highly relativistic shells, with energy $E_s=10^{41}$ erg and Lorentz factor $\gamma_s$~$10^3$, into a surrounding e-p plasma with density $n\sim0.1/cm^3$ (consistent with that inferred for the plasma producing the persistent emission associated with FRB 121102). Such shell ejections with energy typical for FRB events lead to plasma conditions appropriate for strong synchrotron maser emission at the GHz range, $\nu_{ coh.}\sim0.5(E/10^{41}erg)^{1/4}$ GHz. In this model, a significant fraction of the deposited energy is converted to an FRB with duration $R_s/c$, accompanied by ~10 MeV photons carrying less energy than the FRB. The inferred energy and mass associated with the source are low compared to those of typical supernova ejecta. This may suggest some type of a "weak stellar explosion", where a neutron star is formed with relatively low mass and energy ejection. However, the current upper limit on R does not allow one to rule out $M_c\sim1M_\odot$.

The earliest supernova (SN) emission is produced when the optical depth of the plasma lying ahead of the shock, which ejects the envelope, drops below c/v, where v is the shock velocity. This "breakout" may occur when the shock reaches the edge of the star, producing a bright X-ray/UV flash on time scales of seconds to a fraction of an hour, followed by UV/optical "cooling" emission from the expanding cooling envelope on a day time-scale. If the optical depth of circumstellar material (CSM) ejected from the progenitor star prior to the explosion is larger than c/v, the breakout will take place at larger radii, within the CSM, extending its duration to days time scale. The properties of the early, breakout and cooling, emission carry unique signatures of the structure of the progenitor star (e.g. its radius and surface composition) and of its mass-loss history. The recent progress of wide-field transient surveys enable SN detections on a day time scale, and are being used to set unique constraints on the progenitors of SNe of all types. This chapter includes a pedagogical description of SN breakout theory, and a concise overview of what we have learned from observations so far, and of advances in observational capabilities that are required in order to make further significant progress.

The early part of a supernova (SN) light-curve is dominated by radiation escaping from the expanding shock-heated progenitor envelope. For polytropic Hydrogen envelopes, the properties of the emitted radiation are described by simple analytic expressions and are nearly independent of the polytropic index, $n$. This analytic description holds at early time, $t<$~few days, during which radiation escapes from shells initially lying near the stellar surface. We use numerical solutions to address two issues. First, we show that the analytic description holds at early time also for non-polytropic density profiles. Second, we extend the solutions to later times, when the emission emerges from deep within the envelope and depends on the progenitor's density profile. Examining the late time behavior of polytropic envelopes with a wide range of core to envelope mass and radius ratios, $0.1\le M_{\rm c}/M_{\rm env}\le10$ and $10^{-3}\le R_{\rm c}/R\le10^{-1}$, we find that the effective temperature is well described by the analytic solution also at late time, while the luminosity $L$ is suppressed by a factor, which may be approximated to better than $20[30]\%$ accuracy up to $t=t_{\rm tr}/a$ by $A\exp[-(at/t_{\rm tr})^\alpha]$ with $t_{\rm tr} = 15 (M_{\rm env}/M_\odot)^{3/4}(E/10^{51}{\rm erg})^{-1/4}~\rm d$, $A=0.9[0.8]$, $a=1.7[4.6]$ and $\alpha=0.8[0.7]$ for $n=3/2[3]$. This description holds as long as the opacity is approximately that of a fully ionized gas, i.e. for $T>0.7~\rm eV$, $t<14(R/10^{13.5}{\rm cm})^{0.55}~\rm d$. The suppression of $L$ at $t_{\rm tr}/a$ obtained for standard polytropic envelopes may account for the first optical peak of double-peaked SN light curves, with first peak at a few days for $M_{\rm env}<1M_\odot$.

We consider limits on the local ($z=0$) density ($n_0$) of extragalactic neutrino sources set by the nondetection of steady high-energy neutrino sources producing $\gtrsim50$ TeV muon multiplets in the present IceCube data, taking into account the redshift evolution, luminosity function and neutrino spectrum of the sources. We show that the lower limit depends moderately on source spectra and strongly on redshift evolution. We find $n_0\gtrsim{10}^{-8}-{10}^{-7}~{\rm Mpc}^{-3}$ for standard candle sources evolving rapidly, $n_s\propto{(1+z)}^3$, and $n_0\gtrsim{10}^{-6}-{10}^{-5}~{\rm Mpc}^{-3}$ for nonevolving sources. The corresponding upper limits on their neutrino luminosity are $L_{{\nu_\mu}}^{\rm eff}\lesssim10^{42}-10^{43}~{\rm erg}~{\rm s}^{-1}$ and $L_{{\nu_\mu}}^{\rm eff}\lesssim10^{41}-10^{42}~{\rm erg}~{\rm s}^{-1}$, respectively. Applying these results to a wide range of classes of potential sources, we show that powerful blazar jets associated with active galactic nuclei are unlikely to be the dominant sources. For almost all other steady candidate source classes (including starbursts, radio galaxies, and galaxy clusters and groups), an order of magnitude increase in the detector sensitivity at $\sim0.1-1$ PeV will enable a detection (as point sources) of the few brightest objects. Such an increase, which may be provided by next-generation detectors like IceCube-Gen2 and an upgraded KM3NET, can improve the limit on $n_0$ by more than two orders of magnitude. Future gamma-ray observations (by Fermi, HAWC and CTA) will play a key role in confirming the association of the neutrinos with their sources.

We compare the isotropic equivalent 15-2000 keV gamma-ray energy, E_gamma, emitted by a sample of 91 swift Gamma-Ray Bursts (GRBs) with known redshifts, with the isotropic equivalent fireball energy, E_fb, as estimated within the fireball model framework from X-ray afterglow observations of these bursts. The uncertainty in E_gamma, which spans the range of ~10^51 erg to ~10^53.5 erg, is approximately 25% on average, due mainly to the extrapolation from the BAT detector band to the 15-2000 keV band. The uncertainty in E_fb is approximately a factor of 2, due mainly to the X-ray measurements' scatter. We find E_gamma and E_fb to be tightly correlated. The average(std) of {\eta}^11hr_gamma is approximately log_10(E_gamma/(3{\epsilon} _eE^11hr_fb)) are -0.34(0.60), and the upper limit on the intrinsic spread of {\eta}_gamma is approximately 0.5 ({\epsilon}_e is the fraction of shocked plasma energy carried by electrons and E^x hr_fb is inferred from the X-ray flux at x hours). We also find that E_fb inferred from X-ray observations at 3 and 11 hours are similar, with an average(std) of log_10(E^3hr_fb/E^11hr_fb) of 0.04(0.28). The small variance of {\eta}_gamma implies that burst-to-burst variations in {\epsilon}_e and in the efficiency of fireball energy conversion to gamma-rays are small, and suggests that both are of order unity. The small variance of {\eta}_gamma and the similarity of E^3hr_fb and E^11hr_fb further imply that {\epsilon}_e does not vary significantly with shock Lorentz factor, and that for most bursts the modification of fireball energy during the afterglow phase, by processes such as radiative losses or extended duration energy injection, are not significant. Finally, our results imply that if fireballs are indeed jets, then the jet opening angle satisfies {\theta}>0.1 for most cases. [abridged]

The IceCube collaboration reports a detection of extra-terrestrial neutrinos. The isotropy and flavor content of the signal, and the coincidence, within current uncertainties, of the 50 TeV to 2 PeV flux and the spectrum with the Waxman-Bahcall bound, suggest a cosmological origin of the neutrinos, related to the sources of ultra-high energy, $>10^{10}$ GeV, cosmic-rays (UHECR). The most natural explanation of the UHECR and neutrino signals is that both are produced by the same population of cosmological sources, producing CRs (likely protons) at a similar rate, $E^2d\dot{n}/dE\propto E^{0}$, over the [$1$ PeV,$10^{11}$ GeV] energy range, and residing in "calorimetric" environments, like galaxies with high star formation rate, in which $E/Z<100$ PeV CRs lose much of their energy to pion production. A tenfold increase in the effective mass of the detector at $\gtrsim100$ TeV is required in order to significantly improve the accuracy of current measurements, to enable the detection of a few bright nearby starburst "calorimeters", and to open the possibility of identifying the CR sources embedded within the calorimeters, by associating neutrinos with photons accompanying transient events responsible for their generation. Source identification and a large neutrino sample may enable one to use astrophysical neutrinos to constrain new physics models.

The radius and surface composition of an exploding massive star,as well as the explosion energy per unit mass, can be measured using early UV observations of core collapse supernovae (SNe). We present the first results from a simultaneous GALEX/PTF search for early UV emission from SNe. Six Type II SNe and one Type II superluminous SN (SLSN-II) are clearly detected in the GALEX NUV data. We compare our detection rate with theoretical estimates based on early, shock-cooling UV light curves calculated from models that fit existing Swift and GALEX observations well, combined with volumetric SN rates. We find that our observations are in good agreement with calculated rates assuming that red supergiants (RSGs) explode with fiducial radii of 500 solar, explosion energies of 10^51 erg, and ejecta masses of 10 solar masses. Exploding blue supergiants and Wolf-Rayet stars are poorly constrained. We describe how such observations can be used to derive the progenitor radius, surface composition and explosion energy per unit mass of such SN events, and we demonstrate why UV observations are critical for such measurements. We use the fiducial RSG parameters to estimate the detection rate of SNe during the shock-cooling phase (<1d after explosion) for several ground-based surveys (PTF, ZTF, and LSST). We show that the proposed wide-field UV explorer ULTRASAT mission, is expected to find >100 SNe per year (~0.5 SN per deg^2), independent of host galaxy extinction, down to an NUV detection limit of 21.5 mag AB. Our pilot GALEX/PTF project thus convincingly demonstrates that a dedicated, systematic SN survey at the NUV band is a compelling method to study how massive stars end their life.

The flux, spectrum and angular distribution of the excess neutrino signal detected by IceCube between 50TeV and 2PeV are inconsistent with those expected for Galactic sources. The coincidence of the excess, $E_\nu^2\Phi_\nu=3.6\pm1.2\times10^{-8}(GeV/ cm^2 sr s)$, with the Waxman-Bahcall (WB) bound, $E_\nu^2\Phi_{WB}=3.4\times10^{-8}(GeV/cm^2 sr s)$, is probably a clue to the origin of IceCube's neutrinos. The most natural explanation of this coincidence is that both the neutrino excess and the ultra-high energy, $>10^{19}$ eV, cosmic-ray (UHECR) flux are produced by the same population of cosmologically distributed sources, producing CRs, likely protons, at a similar rate, $E^2 dQ/dE=0.5\times10^{44}(erg/Mpc^3yr)$ (at z=0), across a wide range of energies, from $10^{15}$ eV to $>10^{20}$ eV, and residing in environments (such as starburst galaxies) in which CRs of rigidity $E/Z< 10^{17}$ eV lose much of their energy to pion production. Identification of the neutrino sources will allow one to identify the UHECR accelerators, to resolve open questions related to the accelerator models, and to study neutrino properties (related e.g. to flavor oscillations and coupling to gravity) with an accuracy many orders of magnitude better than is currently possible. The most promising method for identifying the sources is by association of a neutrino with an electromagnetic signal accompanying a transient event responsible for its generation. The neutrino flux that is produced within the sources, and that may thus be directly associated with transient events, may be significantly lower than the total observed neutrino flux, which may be dominated by neutrino production at the environment in which the sources reside.

We show that the recent AMS02 positron fraction measurement is consistent with a secondary origin for positrons, and does not require additional primary sources such as pulsars or dark matter. The measured positron fraction at high energy saturates the previously predicted upper bound for secondary production (Katz et al 2009), obtained by neglecting radiative losses. This coincidence, which will be further tested by upcoming AMS02 data at higher energy, is a compelling indication for a secondary source. Within the secondary model the AMS02 data imply a cosmic ray propagation time in the Galaxy of < Myr and an average traversed interstellar matter density of order 1/cc, comparable to the density of the Milky Way gaseous disk, at a rigidity of 300 GV.

The energy output (per logarithmic interval of particle energies) of Cosmic Rays (CRs) with energies $10{\rm GeV}\lesssim\varepsilon_p\lesssim100{\rm GeV}$ is $\sim 10^{47}\rm erg$ per solar mass of star$-$formation, based on the CR production rate in the Milky Way and in starburst galaxies, implying a generation rate of $\varepsilon_p^2Q\sim 10^{45}\rm erg~Mpc^{-3}~yr^{-1}$ in the local universe. It is only $\sim 10$ times larger than the output, $\varepsilon_p^2 Q=0.5\pm0.2\times 10^{44}\rm erg~Mpc^{-3}~yr^{-1}$, of Ultra High Energy CRs (UHECRs) at energies $10^{10.5}{\rm GeV}<\varepsilon_p<10^{12}\rm GeV$ (obtained assuming they are mostly protons), which in turn is comparable to the lower limit of $\varepsilon_p^2 Q\ge 0.5\times 10^{44}\rm erg~Mpc^{-3}~yr^{-1}$ of high energy CRs with $10^6{\rm GeV}\lesssim\varepsilon_p\lesssim 10^{8}\rm GeV$ implied by the saturation of the Waxman-Bahcall bound by the neutrino excess recently discovered by IceCube. These similarities are consistent with a flat production spectrum, $\varepsilon_p^2 Q\sim \text{const}$ for CRs at all observed energies. If a flat production spectrum is generated by our galaxy, the observed CR flux in the range $10^{6.5}-10^{9.5}\rm GeV$, above the "knee", is suppressed compared to lower energies due to propagation effects rather than acceleration upper limits. As suggested by Parizot and Aublin, the most exciting possibility is that cosmic rays at all energies are emitted from a single type of (unknown) sources, which can not be supernova remnants.

The joint analysis of anisotropy signals and chemical composition of ultra-high energy cosmic rays offers strong potential for shedding light on the sources of these particles. Following up on an earlier idea, this paper studies the anisotropies produced by protons of energy >E/Z, assuming that anisotropies at energy >E have been produced by nuclei of charge Z, which share the same magnetic rigidity. We calculate the number of secondary protons produced through photodisintegration of the primary heavy nuclei. Making the extreme assumption that the source does not inject any proton, we find that the source(s) responsible for anisotropies such as reported by the Pierre Auger Observatory should lie closer than ~20-30, 80-100 and 180-200 Mpc if the anisotropy signal is mainly composed of oxygen, silicon and iron nuclei respectively. A violation of this constraint would otherwise result in the secondary protons forming a more significant anisotropy signal at lower energies. Even if the source were located closer than this distance, it would require an extraordinary metallicity >120, 1600, 1100 times solar metallicity in the acceleration zone of the source, for oxygen, silicon and iron respectively, to ensure that the concomitantly injected protons not to produce a more significant low energy anisotropy. This offers interesting prospects for constraining the nature and the source of ultra-high energy cosmic rays with the increase in statistics expected from next generation detectors.

We construct a simple and robust approach for deriving constraints on magnetic fields in galaxy clusters from rotation measure (RM) maps. Relaxing the commonly used assumptions of a correlation between the magnetic field strength and the plasma density and of a power-law (in wave number) magnetic field power spectrum, and using an efficient numerical analysis method, we test the consistency of a wide range of magnetic field models with RM maps of 11 extended sources in 5 clusters, for which the data were made available to us. We show that the data reveal no indication for a radial dependence of the average magnetic field strength, and in particular no indication for a correlation between the gas density and the field strength. The RM maps of a considerable fraction of the sources either require or are consistent with the presence of a spatially uniform magnetic field of a relatively small strength, 0.02-0.3 muG, which contributes significantly to the RM. The RM maps of all but one source do not require a power-law magnetic field power spectrum, and most are consistent with a power spectrum dominated by a single wave length. The uncertainties in the magnetic field strengths (and spatial correlation lengths) derived from RM maps exceed an order of magnitude (and often more). These uncertainties imply, in particular, that there is no indication in current RM data for a systematic difference between the magnetic field strengths in radio-halo clusters and in radio-quiet clusters. With the improvement expected in the near future of the quality and quantity of RM data, our analysis method will enable one to derive more accurate constraints on magnetic fields in galaxy clusters.

The spectrum of radiation emitted following shock breakout from a star's surface with a power-law density profile $\rho \propto x^n$ is investigated. Assuming planar geometry, local Compton equilibrium and bremsstrahlung emission as the dominant photon production mechanism, numerical solutions are obtained for the photon number density and temperature profiles as a function of time, for hydrogen-helium envelopes. The temperature solutions are determined by the breakout shock velocity $v_0$ and the pre-shock breakout density $\rho_0$, and depend weakly on the value of n. Fitting formulas for the peak surface temperature at breakout as a function of $v_0$ and $\rho_0$ are provided, with $T_{peak}\approx 9.44\exp{[12.63(v_0/c)^{1/2}]}$ eV, and the time dependence of the surface temperature is tabulated. The time integrated emitted spectrum is a robust prediction of the model, determined by $\mathcal{T}_{\rm peak}$ and $v_0$ alone and insensitive to details of light travel time or slight deviations from spherical symmetry. Adopting commonly assumed progenitor parameters, breakout luminosities of ~10^45 erg/s and ~10^44 erg/s in the 0.3-10 keV band are expected for BSG and RSG/He-WR progenitors respectively (T_{peak} is well below the band for RSGs, unless their radius is ~10^13 cm). >30 detections of SN1987A-like (BSG) breakouts are expected over the lifetime of ROSAT and XMM-Newton. An absence of such detections would imply that either the typical parameters assumed for BSG progenitors are grossly incorrect or that their envelopes are not hydrostatic. The observed spectrum and duration of XRF 080109/SN2008D are in tension with a non-relativistic breakout from a stellar surface interpretation.

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

We study the arrival directions of 69 ultra-high energy cosmic rays (UHECRs) observed at the Pierre Auger Observatory (PAO) with energies exceeding 55 EeV. We investigate whether the UHECRs exhibit the anisotropy signal expected if the primary particles are protons that originate in galaxies in the local universe, or in sources correlated with these galaxies. We cross-correlate the UHECR arrival directions with the positions of IRAS-PSCz and 2MASS-6dF galaxies taking into account particle energy losses during propagation. This is the first time that the 6dF survey is used in a search for the sources of UHECRs and the first time that the PSCz survey is used with the full 69 PAO events. The observed cross-correlation signal is larger for the PAO UHECRs than for 94% (98%) of realisations from an isotropic distribution when cross-correlated with the PSCz (6dF). On the other hand the observed cross-correlation signal is lower than that expected from 85% of realisations, had the UHECRs originated in galaxies in either survey. The observed cross-correlation signal does exceed that expected by 50% of the realisations if the UHECRs are randomly deflected by intervening magnetic fields by 5 degrees or more. We propose a new method of analysing the expected anisotropy signal, by dividing the predicted UHECR source distribution into equal predicted flux radial shells, which can help localise and constrain the properties of UHECR sources. We find that the 69 PAO events are consistent with isotropy in the nearest of three shells we define, whereas there is weak evidence for correlation with the predicted source distribution in the two more distant shells in which the galaxy distribution is less anisotropic.

We derive a simple approximate model describing the early, hours to days, UV/optical supernova emission, which is produced by the expansion of the outer <~0.01 solar mass part of the shock-heated envelope, and precedes the optical emission driven by radioactive decay. Our model includes an approximate description of the time dependence of the opacity (due mainly to recombination), and of the deviation of the emitted spectrum from a black body spectrum. We show that the characteristics of the early UV/O emission constrain the radius of the progenitor star, its envelope composition, and the ratio of the ejecta energy to its mass, E/M. For He envelopes, neglecting the effect of recombination may lead to an over estimate of progenitor radius by more than an order of magnitude. We also show that the relative extinction at different wavelengths may be inferred from the light-curves at these wave-lengths, removing the uncertainty in the estimate of progenitor radius due to reddening (but not the uncertainty in E/M due to uncertainty in absolute extinction). The early UV/O observations of the type Ib SN2008D and of the type IIp SNLS-04D2dc are consistent with our model predictions. For SN2008D we find progenitor radius to be approx. 10^11 cm, and an indication that the He envelope contains a significant C/O fraction.

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

We show that converging spherical and cylindrical shock waves may ignite a detonation wave in a combustible medium, provided the radius at which the shocks become strong exceeds a critical radius, R_c. An approximate analytic expression for R_c is derived for an ideal gas equation of state and a simple (power-law-Arrhenius) reaction law, and shown to reproduce the results of numerical solutions. For typical acetylene-air experiments we find R_c~0.1 mm (spherical) and R_c~1 mm (cylindrical). We suggest that the deflagration to detonation transition (DDT) observed in these systems may be due to converging shocks produced by the turbulent deflagration flow, which reaches sub (but near) sonic velocities on scales >>R_c. Our suggested mechanism differs from that proposed by Zel'dovich et al., in which a fine-tuned spatial gradient in the chemical induction time is required to be maintained within the turbulent deflagration flow. Our analysis may be readily extended to more complicated equations of state and reaction laws. An order of magnitude estimate of R_c within a white dwarf at the pre-detonation conditions believed to lead to Type Ia supernova explosions is 0.1 km, suggesting that our proposed mechanism may be relevant for DDT initiation in these systems. The relevance of our proposed ignition mechanism to DDT initiation may be tested by both experiments and numerical simulations.

A unique feature of deflagration-to-detonation (DDT) white dwarf explosion models of SNe of type Ia is the presence of a strong shock wave propagating through the outer envelope. We consider the early emission expected in such models, which is produced by the expanding shock-heated outer part of the ejecta and precedes the emission driven by radioactive decay. We expand on earlier analyses by considering the modification of the pre-detonation density profile by the weak-shocks generated during the deflagration phase, the time evolution of the opacity, and the deviation of the post-shock equation of state from that obtained for radiation pressure domination. A simple analytic model is presented and shown to provide an acceptable approximation to the results of 1D numerical DDT simulations. Our analysis predicts a thousand second long UV/optical flash with a luminosity of ~1 to 3*1e39 erg/s. Lower luminosity corresponds to faster (turbulent) deflagration velocity. The predicted luminosity of the UV flash is an order of magnitude lower than that of earlier estimates, and is expected to be strongly suppressed at times longer than an hour due to the deviation from pure radiation domination.

On May 31, 2011 UT a supernova (SN) exploded in the nearby galaxy M51 (the Whirlpool Galaxy). We discovered this event using small telescopes equipped with CCD cameras, as well as by the Palomar Transient Factory (PTF) survey, and rapidly confirmed it to be a Type II supernova. Our early light curve and spectroscopy indicates that PTF11eon resulted from the explosion of a relatively compact progenitor star as evidenced by the rapid shock-breakout cooling seen in the light curve, the relatively low temperature in early-time spectra and the prompt appearance of low-ionization spectral features. The spectra of PTF11eon are dominated by H lines out to day 10 after explosion, but initial signs of He appear to be present. Assuming that He lines continue to develop in the near future, this SN is likely a member of the cIIb (compact IIb; Chevalier and Soderberg 2010) class, with progenitor radius larger than that of SN 2008ax and smaller than the eIIb (extended IIb) SN 1993J progenitor. Our data imply that the object identified in pre-explosion Hubble Space Telescope images at the SN location is possibly a companion to the progenitor or a blended source, and not the progenitor star itself, as its radius (~10^13 cm) would be highly inconsistent with constraints from our post-explosion photometric and spectroscopic data.

We show that a collisionless shock necessarily forms during the shock breakout of a supernova (SN) surrounded by an optically thick wind. An intense non-thermal flash of <~ MeV gamma rays, hard X-rays and multi-TeV neutrinos is produced simultaneously with and following the soft X-ray breakout emission, carrying similar or larger energy than the soft emission. The non-thermal flash is detectable by current X-ray telescopes and may be detectable out to 10's of Mpc by km-scale neutrino telescopes.

The GeV and TeV emission from M82 and NGC 253 observed by Fermi, HESS, and VERITAS constrains the physics of cosmic rays (CRs) in these dense starbursts. We argue that the gamma rays are predominantly hadronic in origin, as expected by previous studies. The measured fluxes imply that pionic losses are efficient for CR protons in both galaxies: we show that a fraction F_cal ~ 0.2 - 0.4 of the energy injected in high energy primary CR protons is lost to inelastic proton-proton collisions (pion production) before escape, producing gamma rays, neutrinos, and secondary electrons and positrons. We discuss the factor ~2 uncertainties in this estimate, including supernova rate and leptonic contributions to the GeV-TeV emission. We argue that gamma-ray data on ULIRGs like Arp 220 can test whether M82 and NGC 253 are truly calorimetric, and we present upper limits on Arp 220 from the Fermi data. We show that the observed ratio of the GeV to GHz fluxes of the starbursts suggests that non-synchrotron cooling processes are important for cooling the CR electron/positron population. We briefly reconsider previous predictions in light of the gamma-ray detections, including the starburst contribution to the gamma-ray background and CR energy densities. Finally, as a guide for future studies, we list the brightest star-forming galaxies on the sky and present updated predictions for their gamma-ray and neutrino fluxes.

Exact bolometric light curves of supernova shock breakouts are derived based on the universal, non relativistic, planar breakout solutions (Sapir et al. 2011), assuming spherical symmetry, constant Thomson scattering opacity, \kappa, and angular intensity corresponding to the steady state planar limit. These approximations are accurate for progenitors with a scale height much smaller than the radius. The light curves are insensitive to the density profile and are determined by the progenitor radius R, and the breakout velocity and density, v_0 and \rho_0 respectively, and \kappa. The total breakout energy, E_BO, and the maximal ejecta velocity, v_max, are shown to be E_BO=8.0\pi R^2\kappa^-1cv_0 and v_max=2.0v_0 respectively, to an accuracy of about 10%. The calculated light curves are valid up to the time of transition to spherical expansion, t_sph\approx R/4v_0. Approximate analytic expressions for the light curves are provided for breakouts in which the shock crossing time at breakout, t_0=c/\kappa\rho_0v_0^2, is << R/c (valid for R<10^14 cm). Modifications of the flux angular intensity distribution and differences in shock arrival times to the surface, \Delta t_asym, due to moderately asymmetric explosions, affect the early light curve but do not affect v_max and E_BO. For 4v_0<<c, valid for large (RSG) progenitors, L\propto t^{-4/3} at max(\Delta t_asym,R/c)< t<t_sph and R may be accurately estimated from R\approx 2*10^13 (L/10^43 erg s^-1)^{2/5}(t/1 hr)^{8/15}.

The problem of a non-steady planar radiation mediated shock (RMS) breaking out from a surface with a power-law density profile, \rho\propto x^n, is numerically solved in the approximation of diffusion with constant opacity. For an appropriate choice of time, length and energy scales, determined by the breakout opacity, velocity and density, the solution is universal, i.e. depends only on the density power law index n. The resulting luminosity depends weakly on the value of n. An approximate analytic solution, based on the self-similar hydrodynamic solutions and on the steady RMS solutions, is constructed and shown to agree with the numerical solutions as long as the shock is far from the surface, \tau>> c/v_{sh}. Approximate analytic expressions, calibrated based on the exact solutions, are provided, that describe the escaping luminosity as a function of time. These results can be used to calculate the bolometric properties of the bursts of radiation produced during supernova (SN) shock breakouts. For completeness, we also use the exact breakout solutions to provide an analytic approximation for the maximum surface temperature for fast (v_{sh}>~0.1) non-thermal breakouts, and show that it is few times smaller than inferred based on steady-state RMS solutions.

We present numerical solutions of the 2D relativistic hydrodynamics equations describing the deceleration and expansion of highly relativistic conical jets, of opening angles 0.05<theta<0.2, propagating into a medium of uniform density. Jet evolution is followed from a collimated relativistic outflow through to the quasi-spherical non-relativistic phase. We show that relativistic sideways expansion becomes significant beyond the radius R at which the expansion Lorentz factor drops to 1/theta. This is consistent with simple analytic estimates, which predict faster sideways expansion than has been claimed based on earlier numerical modeling. For t>R/c the emission of radiation from the jet blast wave is similar to that of a spherical blast wave carrying the same energy. Thus, the total (calorimetric) energy of GRB blast waves may be estimated with only a small fractional error based on t>R/c observations.

Cosmic-rays with energies exceeding 10^{19} eV are referred to as Ultra High Energy Cosmic Rays (UHECRs). The sources of these particles and their acceleration mechanism are unknown, and for many years have been the issue of much debate. The first part of this review describes the main constraints, that are implied by UHECR observations on the properties of candidate UHECR sources, the candidate sources, and the related main open questions. In order to address the challenges of identifying the UHECR sources and of probing the physical mechanisms driving them, a "multi-messenger" approach will most likely be required, combining electromagnetic, cosmic-ray and neutrino observations. The second part of the review is devoted to a discussion of high energy neutrino astronomy. It is shown that detectors, which are currently under construction, are expected to reach the effective mass required for the detection of high energy extra-Galactic neutrino sources, and may therefore play a key role in the near future in resolving the main open questions. The detection of high energy neutrinos from extra-Galactic sources will not only provide constraints on the identity and underlying physics of UHECR sources, but may furthermore provide information on fundamental neutrino properties.

In a recent analysis it was found that the local (z=0) rate at which gamma-ray bursts (GRBs) produce energy in 1 MeV photons, Q_GRB(z=0), is 300 times lower than the local energy production rate in ultra-high energy cosmic-rays. This may appear to be in contradiction with earlier results, according to which Q_GRB(z=0) is similar to the local energy production rate in >10^{19} eV cosmic-rays, Q_{10EeV}(z=0). This short (1 page) note identifies the origin of the apparent discrepancy and shows that Q_GRB(z=0) \sim Q_{10EeV}(z=0) holds.

Type-IIn supernovae (SNe), which are characterized by strong interaction of their ejecta with the surrounding circumstellar matter (CSM), provide a unique opportunity to study the mass-loss history of massive stars shortly before their explosive death. We present the discovery and follow-up observations of a Type IIn SN, PTF 09uj, detected by the Palomar Transient Factory (PTF). Serendipitous observations by GALEX at ultraviolet (UV) wavelengths detected the rise of the SN light curve prior to the PTF discovery. The UV light curve of the SN rose fast, with a time scale of a few days, to a UV absolute AB magnitude of about -19.5. Modeling our observations, we suggest that the fast rise of the UV light curve is due to the breakout of the SN shock through the dense CSM (n~10^10 cm^-3). Furthermore, we find that prior to the explosion the progenitor went through a phase of high mass-loss rate (~0.1 solar mass per year) that lasted for a few years. The decay rate of this SN was fast relative to that of other SNe IIn.

We derive expressions for the time integrated spectrum of Cosmic Rays (CRs) that are accelerated in a decelerating relativistic shock wave and escape ahead of the shock. It is assumed that at any given time the CRs have a power law form, carry a constant fraction of the energy E_tot of the shocked plasma, and escape continuously at the maximal energy attainable. The spectrum of escaping particles is highly sensitive to the instantaneous spectral index due to the fact that the minimal energy, E_min ~ \Gamma^2 m_pc^2 where \Gamma is the shock Lorentz factor, changes with time. In particular, the escaping spectrum may be considerably harder than the canonical N(E)\propto E^-2 spectrum. For a shock expanding into a plasma of density n, a spectral break is expected at the maximal energy attainable at the transition to non relativistic velocities, E ~ 10^19 (\epsilon_B/0.1)(n/1 cm^-3)^(1/6)(E_tot/10^51 erg)^(1/3) eV where \epsilon_B is the fraction of the energy flux carried by the magnetic field. If ultra-high energy CRs are generated in decelerating relativistic blast waves arising from the explosion of stellar mass objects, their generation spectrum may therefore be different than the canonical N(E)\propto E^-2.

Shock waves driven by the release of energy at the center of a cold ideal gas sphere of initial density rho\propto r^{-omega} approach a self-similar (SLS) behavior, with velocity \dot{R}\propto R^delta, as R->\infty. For omega>3 the solutions are of the second-type, i.e., delta is determined by the requirement that the flow should include a sonic point. No solution satisfying this requirement exists, however, in the 3\leq omega\leq omega_{g}(gamma) ``gap'' (\omega_{g}=3.26 for adiabatic index gamma=5/3). We argue that second-type solutions should not be required in general to include a sonic point. Rather, it is sufficient to require the existence of a characteristic line r_c(t), such that the energy in the region r_c(t)<r<R approaches a constant as R->\infty, and an asymptotic solution given by the SLS solution at r_c(t)<r<R and deviating from it at r<r_c may be constructed. The two requirements coincide for omega>omega_g and the latter identifies delta=0 solutions as the asymptotic solutions for 3\leq omega\leq omega_{g} (as suggested by Gruzinov03). In these solutions, r_c is a C_0 characteristic. It is difficult to check, using numerical solutions of the hydrodynamic equations, whether the flow indeed approaches a delta=0 SLS behavior as R->\infty, due to the slow convergence to SLS for omega~3. We show that in this case the flow may be described by a modified SLS solution, d\ln\dot{R}/d\ln R=delta with slowly varying delta(R), eta\equiv d delta/d\ln R<<1, and spatial profiles given by a sum of the SLS solution corresponding to the instantaneous value of delta and a SLS correction linear in eta. The modified SLS solutions provide an excellent approximation to numerical solutions obtained for omega~3 at large R, with delta->0 (and eta\neq0) for 3\leq omega\leq omega_{g}. (abridged)

For a wide variety of initial and boundary conditions, adiabatic one dimensional flows of an ideal gas approach self-similar behavior when the characteristic length scale over which the flow takes place, $R$, diverges or tends to zero. It is commonly assumed that self-similarity is approached since in the $R\to\infty(0)$ limit the flow becomes independent of any characteristic length or time scales. In this case the flow fields $f(r,t)$ must be of the form $f(r,t)=t^{\alpha_f}F(r/R)$ with $R\propto(\pm t)^\alpha$. We show that requiring the asymptotic flow to be independent only of characteristic length scales imply a more general form of self-similar solutions, $f(r,t)=R^{\delta_f}F(r/R)$ with $\dot{R}\propto R^\delta$, which includes the exponential ($\delta=1$) solutions, $R\propto e^{t/\tau}$. We demonstrate that the latter, less restrictive, requirement is the physically relevant one by showing that the asymptotic behavior of accelerating blast-waves, driven by the release of energy at the center of a cold gas sphere of initial density $\rho\propto r^{-\omega}$, changes its character at large $\omega$: The flow is described by $0\le\delta<1$, $R\propto t^{1/(1-\delta)}$, solutions for $\omega<\omega_c$, by $\delta>1$ solutions with $R\propto (-t)^{1/(\delta-1)}$ diverging at finite time ($t=0$) for $\omega>\omega_c$, and by exponential solutions for $\omega=\omega_c$ ($\omega_c$ depends on the adiabatic index of the gas, $\omega_c\sim8$ for $4/3<\gamma<5/3$). The properties of the new solutions obtained here for $\omega\ge\omega_c$ are analyzed, and self-similar solutions describing the $t>0$ behavior for $\omega>\omega_c$ are also derived.

We present a simple analytic model for the structure of non-relativistic and relativistic radiation mediated shocks. At shock velocities \beta_s\equiv v_s/c\gtrsim 0.1, the shock transition region is far from thermal equilibrium, since the transition crossing time is too short for the production of a black-body photon density (by Bremsstrahlung emission). In this region, electrons and photons (and positrons) are in Compton (pair) equilibrium at temperatures T_s significantly exceeding the far downstream temperature, T_s\gg T_d\approx 2(\varepsilon n_u \hbar^3c^3)^{1/4}. T_s\gtrsim 10 keV is reached at shock velocities \beta_s\approx 0.2. At higher velocities, \beta_s\gtrsim0.6, the plasma is dominated in the transition region by e^\pm pairs and 60 keV\lesssim T_s \lesssim 200 keV. We argue that the spectrum emitted during the breaking out of supernova shocks from the stellar envelopes (or the surrounding winds) of Blue Super Giants and Wolf-Rayet stars, which reach \beta_s>0.1 for reasonable stellar parameters, may include a hard component with photon energies reaching tens or even hundreds of keV. This may account for the X-ray outburst associated with SN2008D, and possibly for other SN-associated outbursts with spectra not extending beyond few 100 keV (e.g. XRF060218/SN2006aj).

The growth rate of long wavelength kinetic instabilities arising due to the interaction of a collimated beam of relativistic particles and a cold unmagnetized plasma are calculated in the ultra relativistic limit. For sufficiently culminated beams, all long wave-length modes are shown to be Weibel-unstable, and a simple analytic expression for their growth rate is derived. For large transverse velocity spreads, these modes become stable. An analytic condition for stability is given. These analytic results, which generalize earlier ones given in the literature, are shown to be in agreement with numerical solutions of the dispersion equation and with the results of novel PIC simulations in which the electro-magnetic fields are restricted to a given k-mode. The results may describe the interaction of energetic cosmic rays, propagating into the far upstream of a relativistic collisionless shock, with a cold unmagnetized upstream. The long wavelength modes considered may be efficient in deflecting particles and could be important for diffusive shock acceleration. It is shown that while these modes grow in relativistic shocks propagating into electron-positron pair plasmas, they are damped in relativistic shocks propagating into electron-proton plasmas with moderate Lorenz factors \Gamma_{sh}\lesssim 100. If these modes dominate the deflection of energetic cosmic rays in electron-positron shocks, it is argued that particle acceleration is suppressed at shock frame energies that are larger than the downstream thermal energy by a factor greater than the shock Lorentz factor.

The structure of relativistic radiation mediated shocks (RRMS) propagating into a cold electron-proton plasma is calculated and analyzed. A qualitative discussion of the physics of relativistic and non relativistic shocks, including order of magnitude estimates for the relevant temperature and length scales, is presented. Detailed numerical solutions are derived for shock Lorentz factors $\Gamma_u$ in the range $6\le\Gamma_u\le30$, using a novel iteration technique solving the hydrodynamics and radiation transport equations (the protons, electrons and positrons are argued to be coupled by collective plasma processes and are treated as a fluid). The shock transition (deceleration) region, where the Lorentz factor $ \Gamma $ drops from $ \Gamma_u $ to $ \sim 1 $, is characterized by high plasma temperatures $ T\sim \Gamma m_ec^2 $ and highly anisotropic radiation, with characteristic shock-frame energy of upstream and downstream going photons of a few~$\times\, m_ec^2$ and $\sim \Gamma^2 m_ec^2$, respectively.Photon scattering is dominated by e$^\pm$ pairs, with pair to proton density ratio reaching $\approx10^2\Gamma_u$. The width of the deceleration region, in terms of Thomson optical depths for upstream going photons, is large, $\Delta\tau\sim\Gamma_u^2$ ($\Delta\tau\sim1$ neglecting the contribution of pairs) due to Klein Nishina suppression of the scattering cross section. A high energy photon component, narrowly beamed in the downstream direction, with a nearly flat power-law like spectrum, $\nu I_\nu\propto\nu^0$, and an energy cutoff at $ \sim \Gamma_u^2 m_ec^2 $ carries a fair fraction of the energy flux at the end of the deceleration region. An approximate analytic model of RRMS, reproducing the main features of the numerical results, is provided.

We present a critical analysis of the observational constraints on, and of the theoretical modeling of, aspects of cosmic ray (CR) generation and propagation in the Galaxy, which are relevant for the interpretation of recent positron and anti-proton measurements. We give simple, analytic, model independent expressions for the secondary pbar flux, and an upper limit for the secondary e+ flux, obtained by neglecting e+ radiative losses, e+/(e+ + e-)<0.2\pm0.1 up to ~300 GeV. These expressions are completely determined by the rigidity dependent grammage, which is measured from stable CR secondaries up to ~150 GeV/nuc, and by nuclear cross sections measured in the laboratory. pbar and e+ measurements, available up to ~100 GeV, are consistent with these estimates, implying that there is no need for new, non-secondary, pbar or e+ sources. The radiative loss suppression factor f_{s,e+} of the e+ flux depends on the e+ propagation in the Galaxy, which is not understood theoretically. A rough, model independent estimate of f_{s,e+} 1/3 can be obtained at a single energy, E\sim20 GeV, from unstable secondary decay and is found to be consistent with e+ measurements, including the positron fraction measured by PAMELA. We show that specific detailed models, that agree with compositional CR data, agree with our simple expressions for the e+ and pbar flux, and that the claims that the positron fraction measured by PAMELA requires new primary e+ sources are based on assumptions, that are not supported by observations. If PAMELA results are correct, they suggest that f_{s,e+} is slightly increasing with energy, which provides an interesting constraint on CR propagation models. We argue that measurements of the e+ to pbar ratio are more useful for challenging secondary production models than the positron fraction.

We investigate the constraints imposed on the luminosity function (LF) of long duration Gamma Ray Bursts (LGRBs) by the flux distribution of bursts detected by the GBM at ~1 MeV, and the implications of the non detection of the vast majority, ~95%, of the LGRBs at higher energy, ~1 GeV, by the LAT detector. We find a LF that is consistent with those determined by BATSE and Swift. The non detections by LAT set upper limits on the ratio R of the prompt fluence at ~1 GeV to that at ~1 MeV. The upper limits are more stringent for brighter bursts, with R<{0.1,0.3,1} for {5,30,60}% of the bursts. This implies that for most bursts the prompt ~1 GeV emission may be comparable to the ~1 MeV emission, but can not dominate it. The value of R is not universal, with a spread of (at least) an order of magnitude around R~10^(-1). For several bright bursts with reliable determination of the photon spectral index at ~1 MeV, the LAT non detection implies an upper limit to the ~100 MeV flux which is <0.1 of the flux obtained by extrapolating the ~1 MeV flux to high energy. For the widely accepted models, in which the ~1 MeV power-law photon spectrum reflects the power-law energy distribution of fast cooling electrons, this suggests that either the electron energy distribution does not follow a power-law over a wide energy range, or that the high energy photons are absorbed. Requiring an order unity pair production optical depth at ~100 MeV sets an upper limit for the Lorentz factor, Gamma<=10^(2.5).

We show that the hard X-ray (HXR) emission observed from several galaxy clusters is naturally explained by a simple model, in which the nonthermal emission is produced by inverse Compton scattering of cosmic microwave background photons by electrons accelerated in cluster accretion shocks: The dependence of HXR surface brightness on cluster temperature is consistent with that predicted by the model, and the observed HXR luminosity is consistent with the fraction of shock thermal energy deposited in relativistic electrons being \lesssim 0.1. Alternative models, where the HXR emission is predicted to be correlated with the cluster thermal emission, are disfavored by the data. The implications of our predictions to future HXR observations (e.g. by NuStar, Simbol-X) and to (space/ground based) gamma-ray observations (e.g. by Fermi, HESS, MAGIC, VERITAS) are discussed.

We argue that the observed correlation between the radio luminosity and the X-ray luminosity in radio emitting galaxy clusters implies that the radio emission is due to secondary electrons that are produced by p-p interactions and lose their energy by emitting synchrotron radiation in a strong magnetic field, B>(8\pi a T_{CMB}^4)^{1/2}\simeq 3\muG. We construct a simple model that naturally explains the correlation, and show that the observations provide stringent constraints on cluster magnetic fields and cosmic rays (CRs): Within the cores of clusters, the ratio beta_{core} between the CR energy (per logarithmic particle energy interval) and the thermal energy is beta_{core}\sim 2*10^{-4}; The source of these CRs is most likely the cluster accretion shock, which is inferred to deposit in CRs ~ 0.1 of the thermal energy it generates; The diffusion time of 100 GeV CRs over scales \gtrsim100 kpc is not short compared to the Hubble time; Cluster magnetic fields are enhanced by mergers to \gtrsim 1 % of equipartition, and decay (to <1 muG) on 1 Gyr time scales. The inferred value of beta_{core} implies that high energy gamma-ray emission from secondaries at cluster cores will be difficult to detect with existing and planned instruments.

We show that the spectral and radial distribution of the nonthermal emission of massive, M>10^{14.5}M_sun, galaxy clusters (GCs) may be approximately described by simple analytic expressions, which depend on the GC thermal X-ray properties and on two model parameter, beta_{core} and eta_e. beta_{core} is the ratio of CR energy density (within a logarithmic CR energy interval) and the thermal energy density at the GC core, and eta_{e(p)} is the fraction of the thermal energy generated in strong collisionless shocks, which is deposited in CR electrons (protons). Using a simple analytic model for the evolution of ICM CRs, which are produced by accretion shocks (primary CRs), we find that beta_{core} ~ eta_{p}/200, nearly independent of GC mass and with a scatter Delta ln(beta_{core}) ~ 1 between GCs of given mass. We show that the HXR and gamma-ray luminosities produced by IC scattering of CMB photons by primary electrons exceed the luminosities produced by secondary particles (generated in hadronic interactions within the GC) by factors ~500(eta_e/eta_p)(T/10 keV)^{-1/2} and ~150(eta_e/eta_p)(T/10 keV)^{-1/2} respectively, where T is the GC temperature. Secondary particle emission may dominate at the radio and VHE (> 1 TeV) gamma-ray bands. Our model predicts, in contrast with some earlier work, that the HXR and gamma-ray emission from GCs are extended, since the emission is dominated at these energies by primary electrons. Our predictions are consistent with the observed nonthermal emission of the Coma cluster for eta_peta_e ~ 0.1. The implications of our predictions to future HXR observations (e.g. by NuStar, Simbol-X) and to (space/ground based) gamma-ray observations (e.g. by Fermi, HESS, MAGIC, VERITAS) are discussed. Finally, we show that our model's results agree with results of detailed numerical calculations.

(abridged) Recently, a new class of radio transients in the 5-GHz band was detected by Bower et al. We present new deep near-Infrared (IR) observations of the field containing these transients, and find no counterparts down to a limiting magnitude of K=20.4 mag. We argue that the bright (>1 Jy) radio transients recently reported by Kida et al. are consistent with being additional examples of the Bower et al. transients. We refer to these groups of events as "long-duration radio transients". The main characteristics of this population are: time scales longer than 30 minute but shorter than several days; rate, ~10^3 deg^-2 yr^-1; progenitors sky surface density of >60 deg^-2 (95% C.L.) at Galactic latitude ~40 deg; 1.4-5 GHz spectral slopes, f_\nu ~ \nu^alpha, with alpha>0; and most notably the lack of any counterparts in quiescence in any wavelength. We rule out an association with many types of objects. Galactic brown-dwarfs or some sort of exotic explosions remain plausible options. We argue that an attractive progenitor candidate for these radio transients is the class of Galactic isolated old neutron stars (NS). We confront this hypothesis with Monte-Carlo simulations of the space distribution of old NSs, and find satisfactory agreement for the large areal density. Furthermore, the lack of quiescent counterparts is explained quite naturally. In this framework we find: the mean distance to events in the Bower et al. sample is of order kpc; the typical distance to the Kida et al. transients are constrained to be between 30 pc and 900 pc (95% C.L.); these events should repeat with a time scale of order several months; and sub-mJy level bursts should exhibit Galactic latitude dependence. We discuss possible mechanisms giving rise to the observed radio emission.

This paper proposes and discusses a test of the chemical composition of ultra-high energy cosmic rays that relies on the anisotropy patterns measured as a function of energy. In particular, we show that if one records an anisotropy signal produced by heavy nuclei of charge Z above an energy E_{thr}, one should record an even stronger (possibly much stronger) anisotropy at energies >E_{thr}/Z due to the proton component that is expected to be associated with the sources of the heavy nuclei. This conclusion remains robust with respect to the parameters characterizing the sources and it does not depend at all on the modelling of astrophysical magnetic fields. As a concrete example, we apply this test to the most recent data of the Pierre Auger Observatory. Assuming that the anisotropy reported above 55EeV is not a statistical accident, and that no significant anisotropy has been observed at energies <10EeV, we show that the apparent clustering toward Cen A cannot be attributed to heavy nuclei. Similar conclusions are drawn regarding the apparent excess correlation with nearby active galactic nuclei. We then discuss a robust lower bound to the magnetic luminosity that a source must possess in order to be able to accelerate particles of charge Z up to 100EeV, L_B>10^{45}Z^{-2}erg/s. Using this bound in conjunction with the above conclusions, we argue that the current PAO data does not support the model of cosmic ray origin in active radio-quiet or even radio-loud galaxies. Finally, we demonstrate that the apparent clustering in the direction of Cen A can be explained by the contribution of the last few gamma-ray bursts or magnetars in the host galaxy thanks to the scattering of the cosmic rays on the magnetized lobes.

We derive simple analytic expressions for the flux and spectrum of ultra-high energy cosmic-rays (UHECRs) predicted in models where the CRs are protons produced by extra-Galactic sources. For a power-law scaling of the CR production rate with redshift and energy, d\dot{n} /dE\propto E^-\alpha (1+z)^m, our results are accurate at high energy, E>10^18.7 eV, to better than 15%, providing a simple and straightforward method for inferring d\dot{n}/dE from the observed flux at E. We show that current measurements of the UHECR spectrum, including the latest Auger data, imply E^2d\dot{n}/dE(z=0)=(0.45\pm0.15)(\alpha-1) 10^44 erg Mpc^-3 yr^-1 at E<10^19.5 eV with \alpha roughly confined to 2\lesseq\alpha<2.7. The uncertainty is dominated by the systematic and statistic errors in the experimental determination of individual CR event energy, (\Delta E/E)_{sys} (\Delta E/E)_{stat} ~20%. At lower energy, d\dot{n}/dE is uncertain due to the unknown Galactic contribution. Simple models in which \alpha\simeq 2 and the transition from Galactic to extra-Galactic sources takes place at the "ankle", E ~10^19 eV, are consistent with the data. Models in which the transition occurs at lower energies require a high degree of fine tuning and a steep spectrum, \alpha\simeq 2.7, which is disfavored by the data. We point out that in the absence of accurate composition measurements, the (all particle) energy spectrum alone cannot be used to infer the detailed spectral shapes of the Galactic and extra-Galactic contributions.

We study relativistic unmagnetized collisionless shocks using unprecedentedly large particle-in-cell simulations of two-dimensional pair plasma. High energy particles accelerated by the shock are found to drive magnetic field evolution on a timescale >10^4 plasma times. Progressively stronger magnetic fields are generated on larger scales in a growing region around the shock. Shock-generated magnetic fields and accelerated particles carry >1% and >10% of the downstream energy flux, respectively. Our results suggest limits on the magnetization of relativistic astrophysical flows.

Ultra high-energy cosmic rays (UHECRs) are believed to be protons accelerated in magnetized plasma outflows of extra-Galactic sources. The acceleration of protons to ~10^{20} eV requires a source power L>10^{47} erg/s. The absence of steady sources of sufficient power within the GZK horizon of 100 Mpc, implies that UHECR sources are transient. We show that UHECR "flares" should be accompanied by strong X-ray and gamma-ray emission, and that X-ray and gamma-ray surveys constrain flares which last less than a decade to satisfy at least one of the following conditions: (i) L>10^{50} erg/s; (ii) the power carried by accelerated electrons is lower by a factor >10^2 than the power carried by magnetic fields or by >10^3 than the power in accelerated protons; or (iii) the sources exist only at low redshifts, z<<1. The implausibility of requirements (ii) and (iii) argue in favor of transient sources with L>10^{50} erg/s.

We suggest that non-thermal emission can be produced by multiple scatterings of the photons between the supernova ejecta and pre-shock material in supernova shock breakout. Such bulk-Comptonization process may significantly change the original thermal photon spectrum, forming a power-law non-thermal component at higher energies. We then show that the luminous X-ray outburst XRO081009 associated with SN2008D is likely to be such shock breakout emission from an ordinary type Ib/c supernova.

We study the anisotropy signature which is expected if the sources of ultra high energy, >10^{19} eV, cosmic-rays (UHECRs) are extragalactic and trace the large scale distribution of luminous matter. Using the PSCz galaxy catalog as a tracer of the large scale structure (LSS), we derive the expected all sky angular distribution of the UHECR intensity. We define a statistic, that measures the correlation between the predicted and observed UHECR arrival direction distributions, and show that it is more sensitive to the expected anisotropy signature than the power spectrum and the two point correlation function. The distribution of the correlation statistic is not sensitive to the unknown redshift evolution of UHECR source density and to the unknown strength and structure of inter-galactic magnetic fields. We show, using this statistic, that recently published >5.7x10^{19} eV Auger data are inconsistent with isotropy at ~98% CL, and consistent with a source distribution that traces LSS, with some preference to a source distribution that is biased with respect to the galaxy distribution. The anisotropy signature should be detectable also at lower energy, >4x10^{19} eV. A few fold increase of the Auger exposure is likely to increase the significance to >99% CL, but not to >99.9% CL (unless the UHECR source density is comparable or larger than that of galaxies). In order to distinguish between different bias models, the systematic uncertainty in the absolute energy calibration of the experiments should be reduced to well below the current ~25%.

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

For more than fifty years, it has been believed that cosmic ray (CR) nuclei are accelerated to high energies in the rapidly expanding shockwaves created by powerful supernova explosions. Yet observational proof of this conjecture is still lacking. Recently, Uchiyama and collaborators reported the detection of small-scale X-ray flares in one such supernova remnant, dubbed 'RX J1713-3946' (a.k.a. G347.3-0.5), which also emits very energetic, TeV (10^12 eV) range, gamma-rays. They contend that the variability of these X-ray 'hotspots' implies that the magnetic field in the remnant is about a hundred times larger than normally assumed; and this, they say, means that the detected TeV range photons were produced in energetic nuclear interactions, providing 'a strong argument for acceleration of protons and nuclei to energies of 1 PeV (10^15 eV) and beyond in young supernova remnants.' We point out here that the existing multiwavelength data on this object certainly do not support such conclusions. Though intriguing, the small-scale X-ray flares are not the long sought-after 'smoking gun' of nucleonic CR acceleration in SNRs.

The prompt gamma-ray emission in gamma-ray bursts is believed to be produced by internal shocks within a relativistic unsteady outflow. The recent detection of prompt optical emission accompanying the prompt gamma-ray emission appears to be inconsistent with this model since the out flowing plasma is expected to be highly optically thick to optical photons. We show here that fluctuations in flow properties on short, ~ 1 ms, time scale, which drive the gamma-ray producing collisions at small radii, are expected to lead to "residual" collisions at much larger radii, where the optical depth to optical photons is low. The late residual collisions naturally account for the relatively bright optical emission. The apparent simultaneity of gamma-ray and optical emission is due to the highly relativistic speed with which the plasma expands. Residual collisions may also account for the X-ray emission during the early "steep decline" phase, where the radius is inferred to be larger than the gamma-ray emission radius. Finally, we point out that inverse-Compton emission from residual collisions at large radii is expected to contribute significantly to the emission at high energy, and may therefore "smear" the pair production spectral cut-off.

We present a simple analytic model for the various contributions to the non-thermal emission from shell type SNRs, and show that this model's results reproduce well the results of previous detailed calculations. We show that the \geq 1 TeV gamma ray emission from the shell type SNRs RX J1713.7-3946 and RX J0852.0-4622 is dominated by inverse-Compton scattering of CMB photons (and possibly infra-red ambient photons) by accelerated electrons. Pion decay (due to proton-proton collisions) is shown to account for only a small fraction, \lesssim10^-2, of the observed flux, as assuming a larger fractional contribution would imply nonthermal radio and X-ray synchrotron emission and thermal X-ray Bremsstrahlung emission that far exceed the observed radio and X-ray fluxes. Models where pion decay dominates the \geq 1 TeV flux avoid the implied excessive synchrotron emission (but not the implied excessive thermal X-ray Bremsstrahlung emission) by assuming an extremely low efficiency of electron acceleration, K_ep \lesssim 10^-4 (K_ep is the ratio of the number of accelerated electrons and the number of accelerated protons at a given energy). We argue that observations of SNRs in nearby galaxies imply a lower limit of K_ep \gtrsim 10^-3, and thus rule out K_ep values \lesssim 10^-4 (assuming that SNRs share a common typical value of K_ep). It is suggested that SNRs with strong thermal X-ray emission, rather than strong non-thermal X-ray emission, are more suitable candidates for searches of gamma rays and neutrinos resulting from proton-proton collisions. In particular, it is shown that the neutrino flux from the SNRs above is probably too low to be detected by current and planned neutrino observatories (Abridged).

Electromagnetic energy losses of charged pions and muons suppress the expected high energy, >1E18 eV, neutrino emission from sources of ultrahigh energy, >1E19 eV, cosmic-rays. We show here that >1E19 eV photons produced in such sources by neutral pion decay may escape the sources, thanks to the Klein-Nishina suppression of the pair production cross section, and produce muon pairs in interactions with the cosmic microwave background. The flux of muon decay neutrinos, which are expected to be associated in time and direction with the electromagnetic emission from the sources, may reach a few percent of the Waxman-Bahcall bound. Their detection may allow one to directly identify the sources of >1E19 eV cosmic-rays, and will provide the most stringent constraints on quantum-gravity-induced Lorentz violation.

We derive constraints that must be satisfied by the sources of ~10^{15} to ~10^{18} eV cosmic rays, under the assumption that the sources are Galactic. We show that while these constraints are not satisfied by ordinary supernovae (SNe), which are believed to be the sources of <10^{15} eV cosmic rays, they may be satisfied by the recently discovered class of trans-relativistic supernovae (TRSNe), which were observed in association with gamma-ray bursts. We define TRSNe as SNe that deposit a large fraction, f_R>10^{-2}, of their kinetic energy in mildly relativistic, \gamma\beta>1, ejecta. The high velocity ejecta enable particle acceleration to ~10^{18} eV, and the large value of f_R (compared to f_R~10^{-7} for ordinary SNe) ensures that if TRSNe produce the observed ~10^{18} eV cosmic ray flux, they do not overproduce the flux at lower energies. This, combined with the estimated rate and energy production of TRSNe, imply that Galactic TRSNe may be the sources of cosmic rays with energies up to ~10^{18}eV .

Measurements of flavor ratios of astrophysical neutrino fluxes are sensitive to the two yet unknown mixing parameters $\theta_{13}$ and $\delta$ through the combination $\sin\theta_{13}\cos\delta$. We extend previous studies by considering the possibility that neutrino fluxes from more than a single type of sources will be measured. We point out that, if reactor experiments establish a lower bound on $\theta_{13}$, then neutrino telescopes might establish an upper bound on $|\cos\delta|$ that is smaller than one, and by that prove that CP is violated in neutrino oscillations. Such a measurement requires several favorable ingredients to occur: (i) $\theta_{13}$ is not far below the present upper bound; (ii) The uncertainties in $\theta_{12}$ and $\theta_{23}$ are reduced by a factor of about two; (iii) Neutrino fluxes from muon-damped sources are identified, and their flavor ratios measured with accuracy of order 10% or better. For the last condition to be achieved with the planned km^3 detectors, the neutrino flux should be close to the Waxman-Bahcall bound. It motivates neutrino telescopes that are effectively about 10 times larger than IceCube for energies of O(100 TeV), even at the expense of a higher energy threshold.

We show that the prompt and afterglow X-ray emission of GRB060218, as well as its early (t<=1 d) optical-UV emission, can be explained by a model in which a radiation- mediated shock propagates through a compact progenitor star into a dense wind. The prompt thermal X-ray emission is produced in this model as the mildly relativistic shock, v/c=0.85 carrying few x 10^49 erg, reaches the wind (Thomson) photosphere, where the post-shock thermal radiation is released and the shock becomes collisionless. Adopting this interpretation of the thermal X-ray emission, a subsequent X-ray afterglow is predicted, due to synchrotron emission and inverse-Compton scattering of SN UV photons by electrons accelerated in the collisionless shock. Early optical-UV emission is also predicted, due to the cooling of the outer \delta M ~10^{-3} M_sun envelope of the star, which was heated to high temperature during shock passage. The observed X-ray afterglow and the early optical-UV emission are both consistent with those expected in this model. Detailed analysis of the early optical-UV emission may provide detailed constraints on the density distribution near the stellar surface.

The predicted thermal flash from SN shock breakout might have been detected for the first time by Swift in GRB 060218/SN 2006aj. The detected thermal X-ray emission in this event implies emergence of a trans-relativistic (TR) SN shock with kinetic energy of E_k>1E49 erg. During TRSN shock breakout, the thermal photons could be "accelerated" by the shock through repeated bulk Compton scattering, forming a nonthermal gamma/X-ray component with dominant energy over thermal one. This mechanism of "photon acceleration" at TRSN shock breakout might also account for gamma-rays in the other similar low-luminosity GRBs, implying that they are atypical GRBs with only TR outflows. TRSNe form a peculiar type of SNe with large kinetic energy, >1E49 erg, in TR ejecta, \Gamma\beta ~2.

Thermal X-ray emission which is simultaneous with the prompt gamma-rays has been detected for the first time from a supernova connected with a gamma-ray burst (GRB), namely GRB060218/SN2006aj. It has been interpreted as arising from the breakout of a mildly relativistic, radiation-dominated shock from a dense stellar wind surrounding the progenitor star. There is also evidence for the presence of a mildly relativistic ejecta in GRB980425/SN1998bw, based on its X-ray and radio afterglow. Here we study the process of repeated bulk Compton scatterings of shock breakout thermal photons by the mildly relativistic ejecta. During the shock breakout process, a fraction of the thermal photons would be repeatedly scattered between the pre-shock material and the shocked material as well as the mildly relativistic ejecta and, as a result, the thermal photons get boosted to increasingly higher energies. This bulk motion Comptonization mechanism will produce nonthermal gamma-ray and X-ray flashes, which could account for the prompt gamma-ray burst emission in low-luminosity supernova-connected GRBs, such as GRB060218. A Monte Carlo code has been developed to simulate this repeated scattering process, which confirms that a significant fraction of the thermal photons get "accelerated" to form a nonthermal component, with a dominant luminosity. This interpretation for the prompt nonthermal emission of GRB060218 may imply that either the usual internal shock emission from highly relativistic jets in these low-luminosity GRBs is weak, or alternatively, that there are no highly relativistic jets in this peculiar class of bursts.

gamma-ray bursts (GRBs) are likely sources of ultra-high energy, >10^{19} eV, protons and high energy, >1 TeV, neutrinos. Large volume detectors of ultra high energy cosmic rays (UHECRs) and high energy neutrinos, which are already operating and are being expanded, may allow to test in the coming few years the predictions of the GRB model for high energy proton and neutrino production. Detection of the predicted signals will allow to identify the sources of UHECRs and will provide a unique probe, which may allow to resolve open questions related to the underlying physics of GRB models. Moreover, detection of GRB neutrinos will allow to test for neutrino properties (e.g., flavor oscillations for which tau's would be a unique signature, and coupling to gravity) with an accuracy many orders of magnitude better than is currently possible.

We show that the radial profiles of the temperature and density of the electrons as well as the magnetic field strength around the massive black hole at the Galactic center, SgrA*, may be constrained directly from existing radio data without any need to make prior assumptions about the dynamics of the emitting gas. The observed spectrum and wavelength-dependent angular size of SgrA* indicate that the synchrotron emission originates from an optically-thick plasma of quasi-thermal electrons. We find that the electron temperature rises above the virial temperature within tens of Schwarzschild radii from the black hole, suggesting that the emitting plasma may be outflowing. Constraints on the electron density profile are derived from polarization measurements. Our best-fit results differ from expectations based on existing theoretical models. However, these models cannot be ruled out as of yet due to uncertainties in the source size measurements. Our constraints could tighten considerably with future improvements in the size determination and simultaneous polarization measurements at multiple wavelengths.

In the past four decades a new type of astronomy has emerged, where instead of looking up into the sky "telescopes" are buried miles underground or deep under water or ice and search not for photons (that is, light), but rather for particles called neutrinos. Neutrinos are nearly massless particles that interact very weakly with matter. The detection of neutrinos emitted by the Sun and by a nearby supernova provided direct tests of the theory of stellar evolution and led to modifications of the standard model describing the properties of elementary particles. At present, several very large neutrino detectors are being constructed, aiming at the detection of the most powerful sources of energy and particles in the universe. The hope is that the detection of neutrinos from these sources, which are extra-Galactic and are most likely powered by mass accretion onto black-holes, will not only allow study of the sources, but, much like solar neutrinos, will also provide new information about fundamental properties of matter.

The construction of large volume detectors of high energy, >1 TeV, neutrinos is mainly driven by the search for extra-Galactic neutrino sources. The existence of such sources is implied by observations of ultra-high energy, >10^{19} eV, cosmic-rays, the origin of which is a mystery. In this lecture I briefly discuss the expected extra-Galactic neutrino signal and the current state of the experimental efforts. Neutrino emission from gamma-ray bursts (GRBs), which are likely sources of both high energy protons and neutrinos, is discussed in some detail. The detection of the predicted GRB neutrino signal, which may become possible in the coming few years, will allow one to identify the sources of ultra-high energy cosmic-rays and to resolve open questions related to the underlying physics of GRB models. Moreover, detection of GRB neutrinos will allow one to test for neutrino properties (e.g. flavor oscillations and coupling to gravity) with an accuracy many orders of magnitude better than is currently possible.

Particle acceleration in collisionless shocks is believed to be responsible for the production of cosmic-rays over a wide range of energies, from few GeV to >10^{20}eV, as well as for the non-thermal emission of radiation from a wide variety of high energy astrophysical sources. A theory of collisionless shocks based on first principles does not, however, exist. Observations of gamma-ray burst (GRB) "afterglows" provide a unique opportunity for diagnosing the physics of relativistic collisionless shocks. Most GRBs are believed to be associated with explosions of massive stars. Their "afterglows," delayed low energy emission following the prompt burst of gamma-rays, are well accounted for by a model in which afterglow radiation is due to synchrotron emission of electrons accelerated in relativistic collisionless shock waves driven by the explosion into the surrounding plasma. Within the framework of this model, some striking characteristics of collisionless relativistic shocks are implied. These include the generation of downstream magnetic fields with energy density exceeding that of the upstream field by ~8 orders of magnitude, the survival of this strong field at distances ~10^{10} skin-depths downstream of the shock, and the acceleration of particles to a power-law energy spectrum, d\log n/d\log E ~ -2, possibly extending to 10^{20}eV. I review in this talk the phenomenological considerations, based on which these characteristics are inferred, and the challenges posed to our current models of particle acceleration and magnetic field generation in collisionless shocks. Some recent theoretical results derived based on the assumption of a self-similar shock structure are briefly discussed.

Observations of gamma-ray burst afterglows suggest that the correlation length of magnetic field fluctuations downstream of relativistic non-magnetized collisionless shocks grows with distance from the shock to scales much larger than the plasma skin depth. We argue that this indicates that the plasma properties are described by a self-similar solution, and derive constraints on the scaling properties of the solution. For example, we find that the scaling of the characteristic magnetic field amplitude with distance from the shock is B \propto D^{s_B} with -1<s_B<=0, that the spectrum of accelerated particles is dn/dE \propto E^{-2/(s_B+1)}, and that the scaling of the magnetic correlation function is <B_i(x)B_j(x+dx)> \propto x^{2s_B} (for x>>D). We show that the plasma may be approximated as a combination of two self-similar components: a kinetic component of energetic particles and an MHD-like component representing "thermal" particles. We argue that the latter may be considered as infinitely conducting, in which case s_B=0 and the scalings are completely determined (e.g. dn/dE \propto E^{-2} and B \propto D^0). Similar claims apply to non- relativistic shocks such as in supernova remnants, if the upstream magnetic field can be neglected. Self-similarity has important implications for any model of particle acceleration and/or field generation. For example, we show that the diffusion function in the angle \mu of momentum p in diffusive shock acceleration models must satisfy D_{\mu\mu}(p,D) = D^{-1}D'_{\mu\mu}(p/D), and that a previously suggested model for the generation of large scale magnetic fields through a hierarchical merger of current-filaments should be generalized. A numerical experiment testing our analysis is outlined (Abridged).

Cosmic ray protons interacting with gas at the mean density of the interstellar medium in starburst galaxies lose energy rapidly via inelastic collisions with ambient nuclei. The resulting pions produce secondary electrons and positrons, high-energy neutrinos, and gamma-ray photons. We estimate the cumulative gamma-ray emission from starburst galaxies. We find a total integrated background above 100 MeV of F_gamma ~ 10^{-6} GeV/cm^2/s/sr and a corresponding specific intensity at GeV energies of nuI_nu ~ 10^{-7} GeV/cm^2/s/sr. Starbursts may thus account for a significant fraction of the extra-galactic $\gamma$-ray background. We show that the FIR-radio correlation provides a strong constraint on the gamma-ray emission from starburst galaxies because pions decay into both gamma-rays and radio-emitting electron/positron pairs. We identify several nearby systems where the potential for observing gamma-ray emission is the most favorable (M82, NGC 253, and IC 342), predict their fluxes, and predict a linear FIR-gamma-ray correlation for the densest starbursts. If established, the FIR-gamma-ray correlation would provide strong evidence for the ``calorimeter'' theory of the FIR-radio correlation and would imply that cosmic rays in starburst galaxies interact with gas at approximately the mean density of the interstellar medium (ISM), thereby providing an important constraint on the physics of the ISM in starbursts.

If cosmic ray protons interact with gas at roughly the mean density of the interstellar medium in starburst galaxies, then pion decay in starbursts is likely to contribute significantly to the diffuse extra-galactic background in both gamma-rays and high energy neutrinos. We describe the assumptions that lead to this conclusion and clarify the difference between our estimates and those of Stecker (2006). Detection of a single starburst by GLAST would confirm the significant contribution of starburst galaxies to the extra-galactic neutrino and gamma-ray backgrounds.

Gamma-ray burst (GRB) afterglow emission is believed to be produced by synchrotron emission of electrons accelerated to high energy by a relativistic collisionless shock propagating into a weakly magnetized plasma. Afterglow observations have been used to constrain the post-shock magnetic field and structure, as well as the accelerated electron energy distribution. Here we show that X-ray afterglow observations on day time scale constrain the pre-shock magnetic field to satisfy B>0.2[n/(1/cc)]^{5/8} mG, where n is the pre-shock density. This suggests that either the shock propagates into a highly magnetized fast, v~10^3 km/s, wind, or that the pre-shock magnetic field is strongly amplified, most likely by the streaming of high energy shock accelerated particles. More stringent constraints may be obtained by afterglow observations at high photon energy at late, >1 d, times.

Although the link between long Gamma Ray Bursts (GRBs) and supernovae (SNe) has been established, hitherto there have been no observations of the beginning of a supernova explosion and its intimate link to a GRB. In particular, we do not know however how a GRB jet emerges from the star surface nor how a GRB progenitor explodes. Here we report on observations of the close GRB060218 and its connection to SN2006aj. In addition to the classical non-thermal emission, GRB060218 shows a thermal component in its X-ray spectrum, which cools and shifts into the optical/UV band as time passes. We interpret these features as arising from the break out of a shock driven by a mildly relativistic shell into the dense wind surrounding the progenitor. Our observations allow us for the first time to catch a SN in the act of exploding, to directly observe the shock break-out and to provide strong evidence that the GRB progenitor was a Wolf-Rayet star.

Short-hard gamma-ray bursts (SHBs) may arise from gravitational wave (GW) driven mergers of double neutron star (DNS) systems. DNSs may be "primordial" or can form dynamically by binary exchange interactions in globular clusters during core-collapse. For primordial binaries, the time delay between formation and merger is expected to be short, tau~0.1 Gyr, implying that the redshift distribution of merger events should follow that of star-formation. We point out here that for dynamically formed DNSs, the time delay between star-formation and merger is dominated by the cluster core-collapse time, rather than by the GW inspiral time, yielding delays comparable to the Hubble time. We derive the redshift distribution of merger events of dynamically formed DNSs, and find it to differ significantly from that typically expected for primordial binaries. The observed redshift distribution of SHBs favors dynamical formation, although a primordial origin cannot be ruled out due to possible detection biases. Future red-shift observations of SHBs may allow to determine whether they are dominated by primordial or dynamically formed DNSs.

We show that starburst galaxies convert efficiently cosmic-rays into pions, which in turn decay into high-energy neutrinos and photons. The cumulative background of GeV neutrinos is 10^{-7}GeV/cm^2/s/sr. Its extrapolation to higher neutrino energies depends on the energy spectrum of the injected cosmic-rays and is proportional to E^{-0.15+-0.1} up to E~0.3PeV and possibly higher neutrino energies. This flux, which constitutes a lower limit to the high energy extra-Galactic neutrino flux, is likely to be detectable by forthcoming km-scale neutrino telescopes.

The ratio of anti-electron to total neutrino flux, $\Phi_{\bar{\nu}_e}:\Phi_\nu$, expected from $p\gamma$ interactions in astrophysical sources is $\le1:15$. We point out that this ratio is enhanced by the decay of $\mu^+\mu^-$ pairs, created by the annihilation of secondary high energy photons from the decay of the neutral pions produced in $p\gamma$ interactions. We show that, under certain conditions, the $\Phi_{\bar{\nu}_e}:\Phi_\nu$ ratio may be significantly enhanced in gamma-ray burst (GRB) fireballs, and that detection at the Glashow resonance of $\bar{\nu}_e$ in kilometer scale neutrino detectors may constrain GRB fireball model parameters, such as the magnetic field and energy dissipation radius.

We estimate minimum energy magnetic fields (B_min) for a sample of galaxies with measured gas surface densities, spanning from normal spirals to starbursts. We show that the ratio of the minimum energy magnetic pressure to the total pressure in the ISM decreases substantially with increasing surface density; for Arp 220 this ratio is ~10^-4. Therefore, if the minimum energy estimate is applicable, magnetic fields in starbursts are dynamically weak compared to gravity, in contrast to normal spiral galaxies. We argue, however, that rapid cooling of relativistic electrons in starbursts invalidates the minimum energy estimate. We critically assess a number of independent constraints on the magnetic field strength in starbursts. In particular, we argue that the existence of the FIR-radio correlation implies that the synchrotron cooling timescale for cosmic ray electrons is much shorter than their escape time from the galactic disk; this in turn implies that the true magnetic field in starbursts is significantly larger than B_min. The strongest argument against such large fields is that one might expect starbursts to have steep radio spectra indicative of strong synchrotron cooling, which is not observed. We show, however, that ionization and bremsstrahlung losses can flatten the nonthermal spectra of starburst galaxies even in the presence of rapid cooling, providing much better agreement with observed spectra. We further demonstrate that ionization and bremsstrahlung losses are likely to be important in shaping the radio spectra of most starbursts at GHz frequencies, thereby preserving the linearity of the FIR-radio correlation. We thus conclude that magnetic fields in starbursts are significantly larger than B_min. We highlight several observations that can test this conclusion.

Electromagnetic (and adiabatic) energy losses of pions and muons modify the flavor ratio (measured at Earth) of neutrinos produced by pion decay in astrophysical sources, $\Phi_{\nu_e}:\Phi_{\nu_\mu}:\Phi_{\nu_\tau}$, from 1:1:1 at low energy to 1:1.8:1.8 at high energy. The transition occurs over 1-2 decades of nuetrino energy, and is correlated with a modification of the neutrino spectrum. For gamma-ray bursts, e.g., the transition is expected at \~100 TeV, and may be detected by km-scale neutrino telescopes. Measurements of the transition energy and energy-width will provide unique probes of the physics of the sources. pion and muon energy losses also affect the ratio of $\bar\nu_e$ flux to total neutrino flux, which may be measured at the W-resonance (6.3 PeV): It is modified from 1/6 (1/15) at low energy to 1/9 (practically 0) at high energy for neutrinos produced in pp ($p\gamma$) interactions.

It has been hypothesized recently that core collapse supernovae are triggered by mildly relativistic jets following observations of radio properties of these explosions. Association of a jet, similar to a gamma-ray burst jet but only slower, allows shock acceleration of particles to high energy and non-thermal neutrino emission from a supernova. Detection of these high energy neutrinos in upcoming kilometer scale Cherenkov detectors may be the only direct way to probe inside these astrophysical phenomena as electromagnetic radiation is thermal and contains little information. Calculation of high energy neutrino signal from a simple and slow jet model buried inside the pre-supernova star is reviewed here. The detection prospect of these neutrinos in water or ice detector is also discussed in this brief review. Jetted core collapse supernovae in nearby galaxies may provide the strongest high energy neutrino signal from point sources.

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

We present a new method for analyzing the global stability of the Sedov-von Neumann-Taylor self-similar solutions, describing the asymptotic behavior of spherical decelerating shock waves, expanding into ideal gas with density \propto r^{-\omega}. Our method allows to overcome the difficulties associated with the non-physical divergences of the solutions at the origin. We show that while the growth rates of global modes derived by previous analyses are accurate in the large wave number (small wavelength) limit, they do not correctly describe the small wave number behavior for small values of the adiabatic index \gamma. Our method furthermore allows to analyze the stability properties of the flow at early times, when the flow deviates significantly from the asymptotic self-similar behavior. We find that at this stage the perturbation growth rates are larger than those obtained for unstable asymptotic solutions at similar [\gamma,\omega]. Our results reduce the discrepancy that exists between theoretical predictions and experimental results.

We consider the emission within the fireball model framework of very high energy, ~1 GeV to >1 TeV photons, on a minute time scale, during the onset of fireball deceleration due to interaction with surrounding medium. Our time dependent numerical model includes exact treatment of electron synchrotron emission, inverse-Compton scattering, pair production, and evolution of electromagnetic cascades (initiated by pair production or photo-production of pions). We find that (i) The 1 GeV--10 GeV flux is not sensitive to model parameters and is ~10^{-7} erg cm^-2 s^-1 for z=1 bursts, well within the detection capabilities of GLAST; (ii) The sub-TeV flux depends on the surrounding medium density and on the fraction of thermal energy carried by the magnetic field, epsilon_B: It ranges from ~10^{-7} erg cm^-2 s^-1 in the case of typical ISM density and epsilon_B <~ 10^{-4} to 10^{-10} erg cm^-2 s^-1 in the case of a source surrounded by a wind and epsilon_B ~ 10^{-0.5}; (iii) The sub-TeV flux is detectable by high energy gamma-ray experiments such as HESS, MAGIC, MILAGRO, and VERITAS; (iv) Combined ~1 keV, ~1 GeV and sub-TeV observations will allow to determine both epsilon_B and the ambient medium density; (v) The spectra depend only weakly on the spectral index of the energy distribution of the accelerated electrons. Pion production energy loss of high energy protons may contribute significantly in the wind case to the luminosity of high energy photons. However, it is difficult to distinguish in this case between the electron and proton contributions since the spectral shape is determined primarily by the energy dependence of the pair production optical depth.

Intergalactic type-Ia supernovae (SNe-Ia) have been discovered recently in rich galaxy clusters, likely the descendants of an intergalactic stellar population found in recent years through a variety of tracers. We estimate the observational signatures of the associated SN remnants (SNRs) in the unusual intracluster medium (ICM) environment. If SNe-Ia still have a circumstellar medium (CSM) at the time of explosion, then their remnants are visible in the optical for ~100-1000 years, with properties similar to young galactic SNRs. In contrast with galactic SNRs, in which the ejecta from the explosion interacts with the ISM, intracluster SNRs become undetectable in the optical once their ejecta passes beyond the CSM and enters the hot and tenuous ICM. If SNe-Ia have a CSM, there should be ~150 young SNRs in the Virgo cluster, with L(H-alpha)~10^{35} erg/s and angular size ~0.1''. We investigate the possibility that members of this SNR population may have recently been detected, but incorrectly identified as intergalactic HII regions. Alternatively, if optical intergalactic SNRs do not exist in Virgo, this will be evidence that SNe-Ia are devoid of a CSM, with implications for progenitor scenarios. Regardless of the presence of a CSM, about 10 older radio SNRs per square degree should be detectable in Virgo, with fluxes of ~0.1 mJy at 1 GHz. Their angular sizes, morphologies, and lack of optical association with distant galaxies can distinguish them from the much more numerous background population. Their detection would provide a measurement of the intracluster SN rate. Observations toward the site of SN1980I, a possibly intergalactic Virgo SN-Ia, can test the existence of a CSM by comparison to our early-time predictions for intergalactic SNR development.

We describe a numerical model constructed for the study of the emission of radiation from relativistic plasma under conditions characteristic, e.g., to gamma-ray bursts (GRB's) and active galactic nuclei (AGN's). The model solves self consistently the kinetic equations for e^\pm and photons, describing cyclo-synchrotron emission, direct Compton and inverse Compton scattering, pair production and annihilation, including the evolution of high energy electromagnetic cascades. The code allows calculations over a wide range of particle energies, spanning more than 15 orders of magnitude in energy and time scales. Our unique algorithm, which enables to follow the particle distributions over a wide energy range, allows to accurately derive spectra at high energies, >100 \TeV. We present the kinetic equations that are being solved, detailed description of the equations describing the various physical processes, the solution method, and several examples of numerical results. Excellent agreement with analytical results of the synchrotron-SSC model is found for parameter space regions in which this approximation is valid, and several examples are presented of calculations for parameter space regions where analytic results are not available.

Afterglow observations are commonly used to determine the parameters of GRB explosions, the energy E, surrounding density n, post-shock magnetic field equipartition fraction \epsilon_B and electron equipartition fraction \epsilon_e, under the frequently made assumption that the efficiency of electron "injection" into relativistic shock acceleration is high, i.e. that the fraction f of electrons which undergo acceleration is f~1. We show that the value of f can not be determined by current observations, since currently testable model predictions for a parameter choice {E'=E/f,n'=n/f,\epsilon'_B=f\epsilon_B,\epsilon'_e=f\epsilon_e} are independent of the value of f for m_e/m_p< f< 1. Current observations imply that the efficiency f is similar for highly relativistic and for sub relativistic shocks, and plausibly suggest that f~1, quite unlike the situation in the Crab Nebula. However, values m_e/m_p< f<< 1 can not be ruled out, implying a factor m_e/m_p uncertainty in determination of model parameters. We show that early, <10 hr, radio afterglow observations, which will be far more accessible in the SWIFT era, may provide constraints on f. Such observations will therefore provide a powerful diagnostic of GRB explosions and of the physics of particle acceleration in collisionless shocks.

A fraction of core collapse supernovae of type Ib/c are associated with Gamma-ray bursts, which are thought to produce highly relativistic jets. Recently, it has been hypothesized that a larger fraction of core collapse supernovae produce slower jets, which may contribute to the disruption and ejection of the supernova envelope, and explain the unusually energetic hypernovae. We explore the TeV neutrino signatures expected from such slower jets, and calculate the expected detection rates with upcoming Gigaton Cherenkov experiments. We conclude that individual jetted SNe may be detectable from nearby galaxies.

The main goal of the construction of large volume, high energy neutrino telescopes is the detection of extra-Galactic neutrino sources. The existence of such sources is implied by observations of ultra-high energy, >10^{19} eV, cosmic-rays (UHECRs), the origin of which is a mystery. The observed UHECR flux sets an upper bound to the extra-Galactic high energy neutrino intensity, which implies that the detector size required to detect the signal in the energy range of 1 TeV to 1 PeV is >=1 giga-ton, and much larger at higher energy. Optical Cerenkov neutrino detectors, currently being constructed under ice and water, are expected to achieve 1 giga-ton effective volume for 1 TeV to 1 PeV neutrinos. Coherent radio Cerenkov detectors (and possibly large air-shower detectors) will provide the >> 1 giga-ton effective volume required for detection at ~10^{19} eV. Detection of high energy neutrinos associated with electromagnetically identified sources will allow to identify the sources of UHECRs, will provide a unique probe of the sources, which may allow to resolve open questions related to the underlying physics of models describing these powerful accelerators, and will provide information on fundamental neutrino properties.

The arguments suggesting an association between the sources of cosmological gamma-ray bursts (GRBs) and the sources of ultra-high energy cosmic rays (UHECRs) are presented. Recent GRB and UHECR observations are shown to strengthen these arguments. Predictions of the GRB model for UHECR production, that may be tested with large area high energy cosmic-ray detectors which are either operating or under construction, are outlined.

We analytically study diffusive particle acceleration in relativistic, collisionless shocks. We find a simple relation between the spectral index s and the anisotropy of the momentum distribution along the shock front. Based on this relation, we obtain s = (3beta_u - 2beta_u*beta_d^2 + beta_d^3) / (beta_u - beta_d) for isotropic diffusion, where beta_u (beta_d) is the upstream (downstream) fluid velocity normalized to the speed of light. This result is in agreement with previous numerical determinations of s for all (beta_u,beta_d), and yields s=38/9 in the ultra-relativistic limit. The spectrum-anisotropy connection is useful for testing numerical studies and for constraining non-isotropic diffusion results. It implies that the spectrum is highly sensitive to the form of the diffusion function for particles travelling along the shock front.

Strong intergalactic shocks are a natural consequence of structure formation in the universe. They are expected to deposit large fractions of their energy in relativistic electrons (xi_e~0.05 according to SNR observations) and magnetic fields (xi_B~0.01 according to cluster halo observations). We calculate the synchrotron emission from such shocks using an analytical model, calibrated with a hydrodynamical LCDM simulation. The resulting signal composes a large fraction of the extragalactic radio background (ERB) below 500 MHz. The associated angular fluctuations dominate the sky for frequencies nu<10 GHz and angular scales arcmin-deg (after a modest removal of point sources), provided that xi_e*xi_B>3*10^-4. The fluctuating signal is most pronounced for nu<500 MHz, dominating the sky even for xi_e*xi_B=5*10^-5. The signal will be easily observable by next generation radio telescopes such as LOFAR and SKA, and is marginally observable with present telescopes. It may be identified using cross-correlations with tracers of large scale structure, possibly even in existing <10 GHz CMB anisotropy maps and high resolution ~1 GHz radio surveys. Detection of the signal will provide the first identification of intergalactic shocks and of the WHIM, and gauge the unknown intergalactic magnetic field. We show that existing observations of the diffuse <500 MHz radio background are well fit by a simple, double-disk Galactic model, precluding a direct identification of the diffuse ERB. Modelling the frequency-dependent anisotropy pattern observed at very low (1-10 MHz) frequencies can disentangle the distributions of Galactic cosmic-rays, ionized gas and magnetic fields. Space missions such as ALFA will thus provide important insight into the structure and composition of our Galaxy (abridged).

Strong intergalactic shocks are a natural consequence of structure formation in the universe. These shocks are expected to deposit large fractions of their energy in relativistic electrons (xi_e~0.05 of the thermal energy according to supernova remnant observations) and magnetic fields (xi_B~0.01 according to cluster halo observations). We discuss the expected synchrotron emission from such shocks, and the observational consequences for next generation radio telescopes such as the Square Kilometer Array. We present an analytical model, calibrated and verified based on a hydrodynamical LCDM simulation. The resulting signal composes a large fraction (up to a few 10%) of the extragalactic radio background below 500 MHz. The associated angular fluctuations, e.g. delta T_l>260(xi_e*xi_B/5*10^-4)(nu/100 MHz)^-3 K for multipoles 400<l<2000, dominate the radio sky for frequencies<10 GHz and angular scales 1 arcmin-1 deg (after a modest removal of point sources), provided that xi_e*xi_B>3*10^-4. The fluctuating signal is most pronounced for nu<500 MHz, dominating the sky there even for xi_e*xi_B=5*10^-5. We find that the signal is easily observable by next generation telescopes such as the SKA, and marginally observable with present telescopes. The signal may be identified using cross-correlations with tracers of large-scale structure (e.g. gamma-ray emission from intergalactic shocks). Detection of the signal will provide the first identification of intergalactic shocks and of the warm-hot intergalactic medium (believed to contain most of the baryons in the low redshift universe), and gauge the unknown intergalactic magnetic field.

(Abr) The origin and structure of magnetic fields in Gamma-Ray Burst (GRB) fireball plasmas are two of the most important open questions in all GRB models. We show that the structure and strength of the magnetic field may be constrained by radio and IR observations of the early afterglow, where plasma effects on the polarization of propagating radiation are significant. We calculate these propagation effects for cold and relativistic plasmas, and find that in the presence of a uniform equipartition field the degree of linear polarization is suppressed, and circular polarization prevails at low frequencies, nu < 1-3 GHz, (2x10^11 Hz < nu < few x 10^14 Hz) in the forward (reverse) shock. At higher frequencies linear polarization dominates. At the frequency of the transition between circular and linear polarization, the net level of polarization is minimal, ~10-20%. These features are nearly independent of the circumburst density. The transition frequency is smaller by a factor of ~10 when the uniform field is much weaker than equipartition. The dependence of these results on viewing geometry, outflow collimation and magnetic field orientation is discussed. When the configuration of the field is entangled over length scales much smaller than the extent of the emitting plasma, the aforementioned effects should not be observed and a linear polarization at the few % level is expected. Polarimetric observations during the early afterglow, and particularly of the reverse shock emission, may therefore place strong constraints on the structure and strength of the magnetic field within the fireball plasma.

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

We present detailed calculations of the prompt spectrum of gamma-ray bursts (GRBs) predicted within the fireball model framework, where emission is due to internal shocks in an expanding relativistic wind. Our time dependent numerical model describes cyclo-synchrotron emission and absorption, inverse and direct Compton scattering, and e^\pm pair production and annihilation (including the evolution of high energy electro-magnetic cascades). It allows, in particular, a self-consistent calculation of the energy distribution of e^\pm pairs produced by photon annihilation, and hence a calculation of the spectra resulting when the scattering optical depth due to pairs, \tau_\pm, is high. We show that emission peaks at ~1MeV for moderate to large \tau_\pm, reaching \tau_\pm ~ 10^2. In this regime of large compactness we find that (i) A large fraction of shock energy can escape as radiation even for large \tau_\pm; (ii) The spectrum depends only weakly on the magnetic field energy fraction; (iii) The spectrum is hard, epsilon^2 dN/depsilon \propto epsilon^\alpha with 0.5<\alpha<1, between the self absorption (epsilon_{ssa}= 10^{0.5\pm0.5} keV) and peak (epsilon_{peak}= 10^{0.5\pm0.5} MeV) photon energy, (iv) and shows a sharp cutoff at ~10 MeV; (v) Thermal Comptonization leads to emission peaking at epsilon_peak>=30 MeV, and can not therefore account for observed GRB spectra. For small compactness, spectra extend to >10 GeV with flux detectable by GLAST, and the spectrum at low energy depends on the magnetic field energy fraction. Comparison of the flux at ~1 GeV and ~ 100 keV may therefore allow to determine the magnetic field strength. For both small and large compactness, the spectra depend only weakly on the spectral index of the energy distribution of accelerated electrons.

Measurements of the diffuse extragalactic gamma-ray background (EGRB) are complicated by a strong Galactic foreground. Estimates of the EGRB flux and spectrum, obtained by modeling the Galactic emission, have produced a variety of (sometimes conflicting) results. The latest analysis of the EGRET data found an isotropic flux I_x=1.45+-0.05 above 100 MeV, in units of 10^-5 s^-1 cm^-2 sr^-1. We analyze the EGRET data in search for robust constraints on the EGRB flux, finding the gamma-ray sky strongly dominated by Galactic foreground even at high latitudes, with no conclusive evidence for an additional isotropic component. The gamma-ray intensity measured towards the Galactic poles is similar to or lower than previous estimates of I_x. The high latitude profile of the gamma-ray data is disk-like for 40<|b[deg]|<70, and even steeper for |b|>70; overall it exhibits strong Galactic features and is well fit by a simple Galactic model. Based on the |b|>40 data we find that I_x<0.5 at a 99% confidence level, with evidence for a much lower flux. We show that correlations with Galactic tracers, previously used to identify the Galactic foreground and estimate I_x, are not satisfactory; the results depend on the tracers used and on the part of the sky examined, because the Galactic emission is not linear in the Galactic tracers and exhibits spectral variations across the sky. The low EGRB flux favored by our analysis places stringent limits on extragalactic scenarios involving gamma-ray emission, such as radiation from blazars, intergalactic shocks and production of ultra-high energy cosmic rays and neutrinos. We suggest methods by which future gamma-ray missions such as GLAST and AGILE could indirectly identify the EGRB.

We show that the recent identification of X-ray emission from SN1998bw is naturally explained as synchrotron emission from a shock driven into the wind surrounding the progenitor by a mildly relativistic shell ejected by the supernova, the existence of which was inferred earlier from radio observations. X-ray observations imply a shell energy E~10^{49.7}erg, and constrain the initial shell velocity \beta*c and normalized wind mass loss rate, \dot{m}=(\dot{M}/10^{-5}M_sun/yr)/(v_w/10^3 km/s), to satisfy \beta^3*\dot{m}~10^{-1.5}. The inferred energy is consistent with energy estimates based on radio observations provided \dot{m}~0.04, in which case radio observations imply \beta~0.8, consistent with the X-ray constraint \beta^3*\dot{m}~10^{-1.5}. While X-ray observations allow to determine the parameters characterizing the pre-explosion wind and the mildly relativistic shell ejected by SN1998bw, they do not provide evidence for existence of an off-axis "standard" GRB jet associated with SN1998bw, that may have produced GRB980425. However, as recently pointed out in (astro-ph/0310320), the lack of observational signatures typically expected to be produced by such an off-axis jet on a 1yr time scale, may be due to a low \dot{m}<0.1, which implies that an off-axis jet will become observable only on >10yr time scale.

The recent detection of an unusually hard spectral component in GRB941017 extending to $\ge200$ MeV is hard to explain as a synchrotron emission from shock-accelerated electrons. It was argued to imply acceleration of protons to ultra-high energy. We show here that the "high energy tail" can be explained as emission from shock-accelerated electrons in the early afterglow epoch, taking into account the effect of synchrotron self-absorption. High energy observations set in this case stringent constraints on model parameters: A lower limit to the total explosion energy $E\gsim5 \times 10^{53}$ erg (assuming spherical symmetry); An upper limit to the density of gas surrounding the explosion, $n\lsim10^{-2}(E/10^{54}{\rm erg}){\rm cm}^{-3}$; A lower limit to the expansion Lorentz factor $\Gamma_i\gsim 200$; and An upper limit to the fraction of thermal energy carried by the magnetic field behind the shock driven into the surrounding medium, $\epsilon_{B,f} <= 10^{-4}$. Such constraints can not be inferred from keV--MeV data alone. The unusually low value of $\epsilon_{B,f}$ and the unusually high ratio $E/n$ may account for the rareness of GRB941017-type high energy tails. Tighter constraints on model parameters may be obtained in the future from optical and sub-TeV observations.

The suggested association between the sources of gamma-ray bursts (GRB's) and the sources of ultra-high energy cosmic rays (UHECR's) is based on two arguments: (i) The average energy generation rate of UHECR's is similar to the gamma-ray generation rate of GRB's, and (ii) The constraints that UHECR sources must satisfy to allow proton acceleration to >10^{20} eV are similar to those inferred for GRB sources from gamma-ray observations. We show that recent GRB and UHECR observations strengthen both arguments, and hence strengthen the suggested association.

The identification of type Ib/c supernovae as GRB progenitors is motivated by the association of GRB980425 with SN1998bw and of GRB030329 with SN2003dh. While the gamma-ray luminosity of GRB030329 was typical to cosmological GRBs, the luminosity of the nearby (40 Mpc) GRB980425 was ~5 orders of magnitude lower. The large luminosity difference is commonly explained by hypothesizing that either SNe Ib/c produce two different classes of GRBs, or that GRB980425 was a typical cosmological GRB jet viewed off-axis. In the latter scenario, strong radio emission, ~10^{30} erg/s Hz, is expected at ~1 yr delay due to jet deceleration to sub-relativistic speed, as observed from GRB970508. The radio luminosity of SN1998bw was 3 orders of magnitude lower than this value. We show that the low radio flux may be consistent with the off-axis jet interpretation, if the density of the wind surrounding the progenitor is lower than typically expected, \dot{m}=(\dot{M}/10^{-5} M_sun/yr)/(v/10^3 km/s)~0.1 instead of \dot{m}>=1. The lower value of \dot{m} is consistent with the observed radio emission from the supernova shock driven into the wind. This interpretation predicts transition to sub-relativistic expansion at ~10 yr delay, with current \~1 mJy 10GHz flux and m_V~23 optical flux, and with ~10 mas angular source size. It also implies that in order to search for the signature of off-axis GRBs associated with nearby Ib/c supernovae, follow up observations should be carried on a multi-yr time scale.

Neutrinos at energies above TeV can serve as probes of the stellar progenitor and jet dynamics of gamma ray bursts arising from stellar core collapses. They can also probe collapses which do not lead to gamma-rays, which may be much more numerous. We calculate detailed neutrino spectra from shock accelerated protons in jets just below the outer stellar envelope, before their emergence. We present neutrino flux estimates from such pre-burst jets for two different massive stellar progenitor models. These should be distinguishable by IceCube, and we discuss the implications.

We calculate the TeV-PeV neutrino fluxes of gamma-ray bursts associated with supernovae, based on the observed association between GRB 030329 and supernova SN 2003dh. The neutrino spectral flux distributions can test for possible delays between the supernova and the gamma-ray burst events down to much shorter timescales than what can be resolved with photons. As an illustrative example, we calculate the probability of neutrino induced muon and electron cascade events in a km scale under-ice detector at the South Pole, from the GRB 030329. Our calculations demonstrate that km scale neutrino telescopes are expected to detect signals that will allow to constrain supernova-GRB models.

The existence of cosmic rays of energies exceeding 10^20 eV is one of the mysteries of high energy astrophysics. The spectrum and the high energy to which it extends rule out almost all suggested source models. The challenges posed by observations to models for the origin of high energy cosmic rays are reviewed, and the implications of recent new experimental results are discussed. Large area high energy cosmic ray detectors and large volume high energy neutrino detectors currently under construction may resolve the high energy cosmic ray puzzle, and shed light on the identity and physics of the most powerful accelerators in the universe.

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

The high energy neutrino signature from proton-proton and photo-meson interactions in a supernova remnant shell ejected prior to a gamma-ray burst provides a test for the precursor supernova, or supranova, model of gamma-ray bursts. Protons in the supernova remnant shell, and photons entrapped from a supernova explosion or a pulsar wind from a fast-rotating neutron star remnant provide ample targets for protons escaping the internal shocks of the gamma-ray burst to interact and produce high energy neutrinos. We calculate the expected neutrino fluxes, which can be detected by current and future experiments.

Magnetars can accelerate cosmic rays to high energies through the unipolar effect, and are also copious soft photon emitters. We show that young, fast-rotating magnetars whose spin and magnetic moment point in opposite directions emit high energy neutrinos from their polar caps through photomeson interactions. We identify a neutrino cut-off band in the magnetar period-magnetic field strength phase diagram, corresponding to the photomeson interaction threshold. Within uncertainties, we point out four possible neutrino emission candidates among the currently known magnetars, the brightest of which may be detectable for a chance on-beam alignment. Young magnetars in the universe would also contribute to a weak diffuse neutrino background, whose detectability is marginal, depending on the typical neutrino energy.

The origin of energetic gamma-ray bursts is still unknown. But the detection of polarization of gamma-rays provides fresh insight into the mechanism driving these powerful explosions.

A pedagogical derivation is presented of the ``fireball'' model of gamma-ray bursts, according to which the observable effects are due to the dissipation of the kinetic energy of a relativistically expanding wind, a ``fireball.'' The main open questions are emphasized, and key afterglow observations, that provide support for this model, are briefly discussed. The relativistic outflow is, most likely, driven by the accretion of a fraction of a solar mass onto a newly born (few) solar mass black hole. The observed radiation is produced once the plasma has expanded to a scale much larger than that of the underlying ``engine,'' and is therefore largely independent of the details of the progenitor, whose gravitational collapse leads to fireball formation. Several progenitor scenarios, and the prospects for discrimination among them using future observations, are discussed. The production in gamma- ray burst fireballs of high energy protons and neutrinos, and the implications of burst neutrino detection by kilometer-scale telescopes under construction, are briefly discussed.

The energy spectra of ultra high energy cosmic rays reported by the AGASA, Fly's Eye, Haverah Park, HiRes, and Yakutsk experiments are all shown to be in agreement with each other for energies below 10^{20} eV (after small adjustments, within the known uncertainties, of the absolute energy scales). The data from HiRes, Fly's Eye, and Yakutsk are consistent with the expected flux suppression above 5\times 10^{19} eV due to interactions of cosmic rays with the cosmic microwave background, the Greisen-Zatsepin-Kuzmin (GZK) "supression," and are inconsistent with a smooth extrapolation of the observed cosmic ray energy spectrum to energies > 5\times 10^{19}$ eV. AGASA data show an excess of events above 10^{20} eV, compared to the predicted GZK suppression and to the flux measured by the other experiments.

We examine the stability of self-similar solutions for an accelerating relativistic blast wave which is generated by a point explosion in an external medium with a steep radial density profile of a power-law index > 4.134. These accelerating solutions apply, for example, to the breakout of a gamma-ray burst outflow from the boundary of a massive star, as assumed in the popular collapsar model. We show that short wavelength perturbations may grow but only by a modest factor <~ 10.

Several high energy, >100 GeV, neutrino telescopes are currently operating or under construction. Their main motivation is the extension of the horizon of neutrino astronomy to cosmological scales. We show that general, model independent, arguments imply that ~1 Gton detectors are required to detect cosmic high energy neutrino sources. Predictions of models of some of the leading candidate sources, gamma-ray bursts and micro-quasars, are discussed, and the question of what can be learned from neutrino observations is addressed.

Structure formation in the intergalactic medium (IGM) produces large-scale, collisionless shock waves, where electrons can accelerate to highly relativistic energies. Such electrons can Compton scatter cosmic microwave background photons up to gamma-ray energies. We study the radiation emitted in this process using a hydrodynamic cosmological simulation of a LCDM universe. This radiation, extending beyond TeV energies, has roughly constant energy flux per decade in photon energy, in agreement with the predictions of Loeb & Waxman (2000). Assuming that a fraction xi_e=0.05 of the shock thermal energy is transferred to relativistic electrons, as inferred from collisionless non-relativistic shocks in the interstellar medium, we find that the radiation energy flux, e^2(dJ/de)~ 50-160 eV cm^-2 s^-1 sr^-1, constitutes ~10% of the extragalactic gamma-ray background. The associated point-sources are too faint to account for the ~60 unidentified EGRET gamma-ray sources, but GLAST should resolve several sources associated with large-scale IGM structures for xi_e~0.03, and many more sources for larger xi_e. The intergalactic origin of the radiation can be verified through a cross-correlation with, e.g., the galaxy distribution that traces the same structure. Its shock-origin may be tested by a cross-correlation with radio synchrotron radiation, emitted as the same electrons gyrate in post-shock magnetic fields. We predict that GLAST and Cherenkov telescopes such as MAGIC, VERITAS and HESS should resolve gamma-rays from nearby (redshifts z < 0.01) rich galaxy clusters, perhaps in the form of a \~5-10 Mpc diameter ring-like emission tracing the cluster accretion shock, with luminous peaks at its intersections with galaxy filaments detectable even at z~0.025.

We show that the collapsar model of gamma-ray bursts results in a series of successive shocks and rarefaction waves propagating in the "cork" of stellar material being pushed ahead of the jet, as it emerges from the massive stellar progenitor. These shocks result in a series of characteristic, increasingly shorter and harder thermal X-ray pulses, as well as a non-thermal gamma-ray pulse, which precede the usual non-thermal MeV gamma-rays. We consider jets escaping from both compact (CO or He dwarf) and blue supergiant stellar progenitors. The period, fluence and hardness of the pulses serves as a diagnostic for the size and density of the outer envelope of the progenitor star.

Several high-frequency peaked BL Lac objects such as Mrk 501 are strong TeV emitters. However, a significant fraction of the TeV gamma rays emitted are likely to be absorbed in interactions with the diffuse IR background, yielding electron-positron pairs. Hence, the observed TeV spectrum must be steeper than the intrinsic one. Using the recently derived intrinsic $\gamma$-ray spectrum of Mrk 501 during its 1997 high state, we study the inverse-Compton scattering of cosmic microwave photons by the resulting electron-positron pairs, which implies the existence of a hitherto undiscovered GeV emission. The typical duration of the GeV emission is determined by the flaring activity time and the energy-dependent magnetic deflection time. We numerically calculate the scattered photon spectrum for different intergalactic magnetic field (IGMF) strengths, and find a spectral turnover and flare duration at GeV energies which are dependent on the field strength. We also estimate the scattered photon flux in the quiescent state of Mrk 501. The GeV flux levels predicted are consistent with existing EGRET upper limits, and should be detectable above the synchrotron -- self Compton (SSC) component with the {\em Gamma-Ray Large Area Space Telescope} ({\em GLAST}) for IGMFs $\lesssim 10^{-16}$ G, as expected in voids. Such detections would provide constraints on the strength of weak IGMFs.

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.

We show that if the sources of ultra-high energy cosmic rays (UHECRs) with energies E>10^{19}eV are associated with galaxies, then the production of UHECRs must occur in transient events. Our galaxy is currently at a dim state in between transients. The time interval between transients is >3x10^4 years if the minimum Lorentz factor to which protons are accelerated in the transients is between 10^3 and 10^7. This constraint is satisfied by gamma-ray bursts.

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.

We consider the effects of collective plasma processes on synchrotron emission from highly relativistic electrons. We find, in agreement with Sazonov (1970), that strong effects are possible also in the absence of a non-relativistic plasma component, due to the relativistic electrons (and protons) themselves. In contrast with Sazonov, who infers strong effects only in cases where the ratio of plasma frequency to cyclotron frequency is much larger than the square of the characteristic electron Lorentz factor, nu_p/nu_B >> gamma^2, we find strong effects also for 1 << nu_p/nu_B << gamma^2. The modification of the spectrum is prominent at frequencies nu < nu_{R*} = nu_p min[gamma, (nu_p/nu_B)^(1/2)], where nu_{R*} generalizes the Razin-Tsytovich frequency, nu_R = gamma nu_p, to the regime nu_p/nu_B << gamma^2. Applying our results to gamma-ray burst (GRB) plasmas, we predict a strong modification of the radio spectrum on minute time scale following the GRB, at the onset of fireball interaction with its surrounding medium, in cases where the ratio of the energy carried by the relativistic electrons to the energy carried by the magnetic field exceeds ~ 10^5. Plausible electron distribution functions may lead to negative synchrotron reabsorption, i.e to coherent radio emission, which is characterized by a low degree of circular polarization. Detection of these effects would constrain the fraction of energy in the magnetic field, which is currently poorly determined by observations, and, moreover, would provide a novel handle on the properties of the environment into which the fireball expands.

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.

We analyze the stability of a relativistic double (forward/reverse) shock system which forms when the fireball of a Gamma-Ray Burst (GRB) impacts on the surrounding medium. We find this shock system to be stable to linear global perturbations for either a uniform or a wind (r^{-2}) density profile of the ambient medium. For the wind case, we calculate analytically the frequencies of the normal modes which could modulate the early short-term variability of GRB afterglows. We find that perturbations in the double shock system could induce oscillatory fluctuations in the observed flux on short (down to seconds) time scales during the early phase of an afterglow.

Core collapse of massive stars resulting in a relativistic fireball jet which breaks through the stellar envelope is a widely discussed scenario for gamma-ray burst production. For very extended or slow rotating stars, the fireball may be unable to break through the envelope. Both penetrating and choked jets will produce, by photo-meson interactions of accelerated protons, a burst of neutrinos with energies in excess of 5 TeV while propagating in the envelope. The predicted flux, from both penetrating and chocked fireballs, should be easily detectable by planned cubic kilometer neutrino telescopes.

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.

We present X-ray observations of the afterglow of GRB000926 performed around and after the break observed in the optical light curve two days after the burst. The steep X-ray light curve observed around the break confirms the presence of this feature in X-rays. However, the spectral and temporal properties are not consistent with a standard jet scenario based on synchrotron emission, requiring a more complicated model. We find that X-ray and optical data are compatible with a moderately collimated fireball (with opening angle $\theta\approx25 \deg$) expanding in a dense medium ($n\approx 4\times 10^4$ cm$^{-3}$). This produces two breaks in the light curve. The first, at $t\approx$ 2 days, is due to jet behavior. The second, around 5 days, is attributed to the transition of the fireball to a non-relativistic expansion. This transition predicts a{\it flattening} of the light curve, that explains the late X-ray measurement in excess above the extrapolation of the simple jet scenario, and is also consistent with optical data.

We present an analysis of the Internal Shock Model of GRBs, where gamma-rays are produced by internal shocks within a relativistic wind. We show that observed GRB characteristics impose stringent constraints on wind and source parameters. We find that a significant fraction, of order 20 %, of the wind kinetic energy can be converted to radiation, provided the distribution of Lorentz factors within the wind has a large variance and provided the minimum Lorentz factor is higher than 10^(2.5)L_(52)^(2/9), where L=10^(52)L_(52)erg/s is the wind luminosity. For a high, >10 %, efficiency wind, spectral energy breaks in the 0.1 to 1 MeV range are obtained for sources with dynamical time R/c < 1 ms, suggesting a possible explanation for the observed clustering of spectral break energies in this range. The lower limit to wind Lorenz factor and the upper limit, around (R/10^7 cm)^(-5/6) MeV to observed break energies are set by Thomson optical depth due to electron positron pairs produced by synchrotron photons. Natural consequences of the model are absence of bursts with peak emission energy significantly exceeding 1 MeV, and existence of low luminosity bursts with low, 1 keV to 10 keV, break energies.

A review is presented of the fireball model of gamma-ray bursts (GRBs), and of the production in GRB fireballs of high energy protons and neutrinos. Constraints imposed on the model by recent afterglow observations, which support the association of GRB and ultra-high energy cosmic-ray (UHECR) sources, are discussed. Predictions of the GRB model for UHECR production, which can be tested with planned large area UHECR detectors and with planned high energy neutrino telescopes, are reviewed.

Observations imply that gamma-ray bursts (GRBs) are produced by the dissipation of the kinetic energy of a highly relativistic fireball. Photo-meson interactions of protons with gamma-rays within the fireball dissipation region are expected to convert a significant fraction of fireball energy to >10^14 eV neutrinos. We present an analysis of the internal shock model of GRBs, where production of synchrotron photons and photo-meson neutrinos are self-consistently calculated, and show that the fraction of fireball energy converted to high energy neutrinos is not sensitive to uncertainties in fireball model parameters, such as the expansion Lorentz factor and characteristic variability time. This is due in part to the constraints imposed on fireball parameters by observed GRB characteristics, and in part to the fact that for parameter values for which the photo-meson optical depth is high (implying high proton energy loss to pion production) neutrino production is suppressed by pion and muon synchrotron losses. The neutrino flux is therefore expected to be correlated mainly with the observed gamma-ray flux. The time averaged neutrino intensity predicted by the model, ~10^-8.5 GeV/cm^2 s sr, is consistent with the flux predicted by the assumption that GRBs are the sources of >10^19 eV cosmic-rays.

We show that as a Type II supernova shock breaks out of its progenitor star, it becomes collisionless and may accelerate protons to energies >10 TeV. Inelastic nuclear collisions of these protons produce a ~1 hr long flash of TeV neutrinos and 10 GeV photons, about 10 hr after the thermal (10 MeV) neutrino burst from the cooling neutron star. A Galactic supernova in a red supergiant star would produce a photon and neutrino flux of ~10^{-4} erg/cm^2 s. A km^2 neutrino detector will detect ~100 muons, thus allowing to constrain both supernova models and neutrino properties.

We elucidate the physical basis for the upper bound on high energy neutrino fluxes implied by the observed cosmic ray flux. We stress that the bound is valid for neutrinos produced either by p,gamma reactions or by p-p(n) reactions in sources which are optically thin for high energy protons to photo-meson and nucleon-meson interactions. We show that the upper bound is robust and conservative. The Waxman-Bahcall bound overestimates the most likely neutrino flux by a factor ~ 5/tau, for small optical depths tau. The upper limit cannot be plausibly evaded by invoking magnetic fields, optically thick AGNs, or large hidden fluxes of extragalactic protons. We describe the implications of the bound for future experiments including the AMANDA, ANTARES, Auger, ICECUBE, NESTOR, and OWL/AIRWATCH detectors.

Observations suggest that gamma-ray bursts (GRBs) are produced by the dissipation of the kinetic energy of a relativistic fireball. In this talk, recent work on the production of high energy neutrinos by GRB fireballs is reviewed. A significant fraction of GRB energy is expected to be converted to an accompanying burst of high energy neutrinos. Photomeson interactions produce a burst of ~100 TeV neutrinos in coincidence with the GRB, and a burst of \~10^18 eV neutrinos following the GRB on a time scale of 10 s. Inelastic p-n nuclear collisions result in the production of a burst of ~10 GeV neutrinos in coincidence with the GRB. Planned 1 km^3 neutrino telescopes are expected to detect tens of 100 TeV neutrino events, and several 10^18 eV events, correlated with GRBs per year. A suitably densely spaced detector may allow the detection of several 10 GeV events per year. The detection of high-energy neutrino events correlated with GRBs will allow to constrain GRB progenitor models and to test the suggestion that GRBs accelerate protons to >10^20 eV. Moreover, such detection will allow to test for neutrino properties, e.g. flavor oscillations (for which upward moving tau's would be a unique signature) and coupling to gravity, with an accuracy many orders of magnitude better than is currently possible.

The shocks produced in the intergalactic medium during large-scale structure formation accelerate a population of highly relativistic electrons which emit synchrotron radiation due to intergalactic magnetic fields. In a previous paper (Loeb & Waxman 2000) we have shown that these electrons cool primarily by inverse-Compton scattering of the microwave background photons and can thereby produce the observed intensity and spectrum of the diffuse gamma-ray background. Here we calculate the intensity and angular fluctuations of the radio synchrotron background that results from the same high-energy electrons, as well as the expected angular fluctuations in the gamma-ray background. On angular scales smaller than a degree, the predicted fluctuations in the microwave background temperature are of order 40micro-K*(xi_B/0.01)(nu/10 GHz)^{-3}, where xi_B is the magnetic fraction of the post-shock energy density. This foreground might have already dominated the anisotropy signal detected in existing low-frequency CMB experiments, and can be identified with confidence through multi-frequency observations. Detection of the synchrotron fluctuations would allow to determine the strength of the intergalactic magnetic field. We predict a strong correlation between high-resolution maps taken at low-frequency radio waves and at high-energy gamma-rays. Young X-ray clusters may also appear as radio or gamma-ray clusters. The detailed study of this correlation will become easily accessible with the future launch of GLAST.

Clustering of cosmic-ray sources affects the flux observed beyond the cutoff imposed by the cosmic microwave background and may be important in interpreting the AGASA, Fly's Eye, and HiRes data. The standard deviation, sigma, in the predicted number, N, of events above 10^{20} eV is sigma/N = 0.9(r_0/10 Mpc)^{0.9}, where r_0 is the unknown scale length of the correlation function (r_0 = 10 Mpc for field galaxies). Future experiments will allow the determination of r_0 through the detection of anisotropies in arrival directions of ~ 10^{20} eV cosmic-rays over angular scales of Theta ~ r_0/30 Mpc.

We show that a significant fraction of the energy of a gamma-ray burst(GRB) is probably converted to a burst of 10^{17}-10^{19} eV neutrinos and multiple GeV gammas that follow the GRB by > 10 s . If, as previously suggested, GRB's accelerate protons to ~10^{20} eV, then both the neutrinos and the gammas may be detectable.

Ultra-high-energy, >10^19 eV, cosmic-ray and high energy, ~10^14 eV, neutrino production in GRBs is discussed in the light of recent GRB and cosmic-ray observations. Emphasis is put on model predictions that can be tested with operating and planned cosmic-ray and neutrino detectors, and on the prospects of testing for neutrino properties.

The universe is filled with a diffuse and isotropic extragalactic background of gamma-ray radiation, containing roughly equal energy flux per decade in photon energy between 3 MeV-100 GeV. The origin of this background is one of the unsolved puzzles in cosmology. Less than a quarter of the gamma-ray flux can be attributed to unresolved discrete sources, but the remainder appears to constitute a truly diffuse background whose origin has hitherto been mysterious. Here we show that the shock waves induced by gravity during the formation of large-scale structure in the intergalactic medium, produce a population of highly-relativistic electrons with a maximum Lorentz factor above 10^7. These electrons scatter a small fraction of the microwave background photons in the present-day universe up to gamma-ray energies, thereby providing the gamma-ray background. The predicted diffuse flux agrees with the observed background over more than four decades in photon energy, and implies a mean cosmological density of baryons which is consistent with Big-Bang nucleosynthesis.

The gamma-ray burst (GRB) model for production of ultra-high-energy cosmic-rays (UHECRs) is based on the hypothesis that GRBs arise from the dissipation of the kinetic energy of relativistic fireballs at cosmological distances. Recent GRB afterglow observations support the validity of this hypothesis and provide quantitative support for the model. The inferred physical fireball parameters imply that protons may be accelerated to >10^20 eV, and the inferred GRB energy generation rate is similar to that required to account for the observed flux of UHECRs. Strong suppression of UHECR flux is expected in this model above 10^19.7 eV, due to proton interaction with CMB photons. Strong deviations from model flux derived under the assumption of uniform source distribution is expected above 10^20 eV, due to source clustering. The flux above 10^20.5 eV is expected to be dominated by few, narrow spectrum sources. While model predictions can not be tested (with high confidence level) at present, the predicted signatures should be observed with the planned Auger and Telescope-Array UHECR detectors. A natural consequence of the GRB model of UHECR production is the conversion of a large fraction, ~10%, of the fireball energy to accompanying burst of \~10^14 eV and ~10^18 eV neutrinos. A km^2 neutrino detector would observe several tens of events per year correlated with GRBs, and test for neutrino properties (e.g. flavor oscillations, for which upward moving tau's would be a unique signature, and coupling to gravity) with an accuracy many orders of magnitude better than is currently possible.

We show that gamma-ray burst (GRB) afterglow observations strongly suggest, within the fireball model framework, that radiating electrons are shock accelerated to a power-law energy distribution, with universal index p \approx 2.2, and that the fraction of shock energy carried by electrons, \xi_e, is universal and close to equipartition, \xi_e ~ 1/3. For universal p and \xi_e, a single measurement of the X-ray afterglow flux on the time scale of a day provides a robust estimate of the fireball energy per unit solid angle, \epsilon, averaged over a conical section of the fireball of opening angle \theta ~ 0.1. Applying our analysis to BeppoSAX afterglow data we find that: (i) Fireball energies are in the range of 4\pi\epsilon=10^{51.5} to 10^{53.5} erg; (ii) The ratio of observed $\gamma$-ray to total fireball energy per unit solid angle, \epsilon_\gamma / \epsilon, is of order unity, satisfying abs[log10(\epsilon_\gamma/\epsilon)]<0.5; (iii) If fireballs are jet like, their opening angle should satisfy \theta>=0.1. Our results imply that if typical opening angles are \theta ~ 0.1, a value consistent with our analysis, the total energy associated with a GRB event is in the range of 10^{50} erg to 10^{51.5} erg.

We report on the results of an extensive monitoring campaign of the radio afterglow of GRB 970508, lasting 450 days after the burst. The spectral and temporal radio behavior indicate that the fireball has undergone a transition to sub-relativistic expansion at t~100 days. This allows us to perform "calorimetry" of the explosion. The derived total energy, ~5\times 10^{50} erg, is well below the ~5\times 10^{51} erg inferred under the assumption of spherical symmetry from gamma-ray and early afterglow observations. A natural consequence of this result, which can also account for deviations at t<100 days from the spherical relativistic fireball model predictions, is that the fireball was initially a wide-angle jet of opening angle ~30 degrees. Our analysis also allows to determine the energy fractions carried by electrons and magnetic field, and the density of ambient medium surrounding the fireball. We find that during the sub-relativistic expansion electrons and magnetic field are close to equipartition, and that the density of the ambient medium is ~1/cm^3. The inferred density rules out the possibility that the fireball expands into a strongly non-uniform medium, as would be expected, e.g., in the case of a massive star progenitor.

The recent detection of delayed, low energy emission from Gamma-Ray Burst (GRB) sources confirmed the cosmological origin of the bursts and provided support for models where GRBs are produced by the dissipation of the kinetic energy of relativistic fireballs. In this review, ultra high energy, >10^{19} eV, cosmic-ray and high energy, 100 TeV, neutrino production in GRBs is discussed in the light of recent GRB and cosmic-ray observations. Emphasis is put on model predictions that can be tested with operating and planned cosmic-ray and neutrino detectors. The predicted neutrino intensity, E^2 dN/dE=3\times 10^{-9} GeV/(cm^2 s sr) for 10^{14} eV<E<10^{16} eV, implies that a km^2 neutrino detector would observe tens of events per year correlated with GRBs, and will be able to test for neutrino properties with an accuracy many orders of magnitude better than is currently possible. The predicted production rate of high-energy protons, which is consistent with that required to account for the observed ultra-high-energy cosmic-ray (UHECR) flux, implies that operating and planned cosmic-ray detectors can test the GRB model for UHECR production. If the predicted sources are found, cosmic-ray detectors will provide us with a technique to investigate the inter-galactic magnetic field.

We consider models for gamma-ray bursts in which a collimated jet expands either into a homogeneous medium or into a stellar wind environment, and calculate the expected afterglow temporal behavior. We show that (i) following a break and a faster decay, afterglows should exhibit a flattening, which may be detectable in both the radio and optical bands; (ii) Only observations at times much shorter than a day can clearly distinguish between a fireball interacting with a homogeneous medium and one interacting with a stellar wind. Using our results we demonstrate that constraints can be placed on progenitor models. In particular, existing data imply that while some long duration bursts may be produced by collapses of massive stars, it is almost certain that not all long duration bursts are produced by such progenitors.

The prompt optical flash recently detected accompanying GRB990123 suggests that, for at least some GRBs, gamma-ray emission is accompanied by prompt optical-UV emission with luminosity L(1-7.5eV)=10^{49}(\Delta\Omega/4\pi)erg/s, where \Delta\Omega is the solid angle into which gamma-ray and optical-UV emission is beamed. Such an optical-UV flash can destroy dust in the beam by sublimation out to an appreciable distance, approximately 10 pc, and may clear the dust out of as much as 10^7(\Delta\Omega/4\pi)M_sun of molecular cloud material on an apparent time scale of 10 seconds. Detection of time dependent extinction on this time scale would therefore provide strong constraints on the GRB source environment. Dust destruction implies that existing, or future, observations of not-heavily-reddened fireballs are not inconsistent with GRBs being associated with star forming regions. In this case, however, if gamma-ray emission is highly beamed, the expanding fireball would become reddened on a 1 week time scale. If the optical depth due to dust beyond approximately 8 pc from the GRB is 0.2<\tau_V<2, most of the UV flash energy is converted to infra-red, \lambda \sim 1 micron, radiation with luminosity \sim 10^{41} erg/s extending over an apparent duration of \sim 20(1+z)(\Delta\Omega/0.01) day. Dust infra-red emission may already have been observed in GRB970228 and GRB980326, and may possibly explain their unusual late time behavior.

We show that cosmic-ray observations set a model-independent upper bound to the flux of high-energy, > 10^14 eV, neutrinos produced by photo-meson (or p-p) interactions in sources of size not much larger than the proton photo-meson (or pp) mean-free-path. The bound applies, in particular, to neutrino production by either AGN jets or GRBs. This upper limit is two orders of magnitude below the flux predicted in some popular AGN jet models, but is consistent with our predictions from GRB models. We discuss the implications of these results for future km^2 high-energy neutrino detectors.

SN1998bw is the most luminous radio supernova ever observed. Previous discussions argued that its exceptional radio luminosity, 4e38 erg/s, must originate from a highly relativistic shock which is fully decoupled from the supernova ejecta. Here we present an alternative model in which the radio emission originates from a sub-relativistic shock with a velocity of 0.3c, generated in the surrounding gas by the expanding ejecta. In this model, thermal electrons heated by the shock to a relativistic temperature of 60 MeV emit synchrotron self-absorbed radiation in the post-shock magnetic field. This model provides an excellent fit to the observed spectra provided that the thermal electrons are in equipartition with the ions behind the shock. The required magnetic field is much weaker than its equipartition value and could have been carried out by the progenitor's wind prior to the supernova explosion. According to this model, the radio emission from SN1998bw is unrelated to the highly relativistic blast wave that produced the gamma-ray burst GRB980425.

GRB afterglow polarization is discussed. We find an observable, up to 10%, polarization, if the magnetic field coherence length grows at about the speed of light after the field is generated at the shock front. Detection of a polarized afterglow would show that collisionless ultrarelativistic shocks can generate strong large scale magnetic fields and confirm the synchrotron afterglow model. Non-detection, at a 1% level, would imply that either the synchrotron emission model is incorrect, or that strong magnetic fields, after they are generated in the shock, somehow manage to stay un-dissipated at ``microscopic'', skin depth, scales. Analytic lightcurves of synchrotron emission from an ultrarelativistic self-similar blast wave are obtained for an arbitrary electron distribution function, taking into account the effects of synchrotron cooling. The peak synchrotron flux and the flux at frequencies much smaller than the peak frequency are insensitive to the details of the electron distribution function; hence their observational determination would provide strong constraints on blast wave parameters.

We present radio observations of the gamma-ray burster GRB 970228 made with the Very Large Array (VLA) and the Owens Valley Radio Observatory (OVRO) spanning a range of postburst timescales from one to 300 days. A search for a time-variable radio source was conducted covering an area which included a fading X-ray source and an optical transient, both of which are thought to be the long wavelength counterparts to the gamma-ray burst. At the position of the optical transient sensitive limits between 10 uJy and 1 mJy can be placed on the absence of a radio counterpart to GRB 970228 between 1.4 and 240 GHz. We apply a simple formulation of a fireball model which has been used with some success to reproduce the behavior of the optical and X-ray light curves. Using this model we conclude that the radio non-detections are consistent with the peak flux density of the afterglow lying between 20-40 uJy and it requires that the optical flux peaked between 4 and 16 hours after the burst.

Annihilation of high energy, $\sim 10^{21}$eV, neutrinos on big bang relic neutrinos of $\sim 1$eV mass, clustered in the Galactic halo or in a nearby galaxy cluster halo, has been suggested to generate, through hadronic Z decay, high energy nucleons and photons which may account for the detected flux of >10^{20}eV cosmic-rays. We show that the flux of high energy nucleons and photons produced by this process is dominated by annihilation on the uniform, non-clustered, neutrino background, and that the energy generation rate of 10^{21}eV neutrinos required to account for the detected flux of >10^{20}eV particles is >10^{48} erg/Mpc^3 yr. This energy generation rate, comparable to the total luminosity of the universe, is 4 orders of magnitude larger than the rate of production of high energy nucleons required to account for the flux of >10^{19}eV cosmic-rays. Thus, in order for neutrino annihilation to contribute significantly to the detected flux of >10^{20}eV cosmic-rays, the existence of a new class of high-energy neutrino sources, likely unrelated to the sources of >10^{19}eV cosmic-rays, must be invoked.

The recent detection of delayed X-ray and optical emission, ``afterglow,'' associated with gamma-ray bursts (GRBs) supports models, where the bursts are produced by relativistic expanding blastwaves, ``fireballs,'' at cosmological distances. The detection of absorption lines in the optical afterglow of the GRB of 8 May 1997 confirms that the sources lie at cosmological distance. We show here that the new features detected in GRB970508 afterglow, radio emission one week following the burst and a 2 day increase in optical flux, are consistent with the blastwave model. The fireball optical depth at radio frequencies is much smaller than previously estimated, which accounts for the observed radio emission. The initial suppression of optical flux is consistent with that predicted due to electron cooling. The combined radio and optical data imply that the fireball energy is $\sim10^{52}$erg, and that the density of the medium into which the blastwave expands is $\sim 1cm^{-3}$, a value typical for gas within galaxies. We predict the time dependence of the radio flux and the absorption frequency, which constitute tests of the fireball model as described in this paper.

Radio observations of the afterglow of the gamma-ray burst GRB970508 provide unique new constraints on afterglow models. The quenching of diffractive scintillation at 4 week delay provides the first direct estimate of source size and expansion rate. It implies an apparent size $R\sim10^{17}$cm and expansion at a speed comparable to that of light at 4 weeks delay, in agreement with the fireball model prediction $R=10^{17}(t/week)^{5/8}$cm. The radio flux and its dependence on time and frequency at 1--5 week delay are in agreement with the model and imply a fireball energy (assuming spherical symmetry) of order 10^{52}erg, consistent with the value inferred from observations at shorter delay. The observed radio behavior deviates from model predictions at delays larger than 5 weeks. This is expected, since at this delay the fireball is in transition from highly-relativistic to sub-relativistic expansion, with Lorentz factor \gamma<2. Deviation may be due to a change in the physical processes associated with the shock wave as it becomes sub-relativistic (e.g. a decrease in the fraction of energy carried by magnetic field), or to the fireball being a cone of opening angle $1/\gamma\sim0.5$. We predict the future behavior of the radio flux assuming that the latter interpretation is valid. These predictions may be tested by radio observations in the frequency range 0.1--10GHz on time scale of months.

The detection of delayed X-ray, optical and radio emission, ``afterglow,'' associated with gamma-ray bursts (GRBs) is consistent with models, where the bursts are produced by relativistic expanding blast waves, driven by expanding fireballs at cosmological distances. In particular, the time scales over which radiation is observed at different wave bands agree with model predictions. It had recently been claimed that the commonly used relation between observation time t and blast wave radius r, t=r/2\gamma^{2}c where \gamma(r) is the fluid Lorentz factor, should be replaced with t=r/16\gamma^{2}c due to blast wave deceleration. Applying the suggested deceleration modification would make it difficult to reconcile observed time scales with model predictions. It would also imply an apparent source size which is too large to allow attributing observed radio variability to diffractive scintillation. We present a detailed analysis of the implications of the relativistic hydrodynamics of expanding blast waves to the observed afterglow. We find that modifications due to shock deceleration are small, therefore allowing for both the observed afterglow time scales and for diffractive scintillation. We show that at time t the fireball appears on the sky as a narrow ring of radius h=r/\gamma and width 0.1h, where r and t are related by t=r/2\gamma^{2}c.

Cosmological fireball models of gamma-ray bursts (GRBs) predict delayed emission, ``after-glow,'' at longer wavelengths. We present several new results regarding the model predictions, and show that X-ray to optical observations of GRB970228 and GRB970402 are naturally explained by the model: The scaling of flux with time and frequency agrees with model predictions and requires a power law distribution of shock accelerated electrons, dlog N/dlog E=2.3+-0.1 (implying, and consistent with the observed, 1/t decline of flux observed at a given frequency); The absolute flux value agrees with that inferred through the model from observed gamma-ray fluence. The future after-glow emission of these bursts is predicted. The observations indicate that the ratio of magnetic field to equipartition value and the fraction F of dissipated kinetic energy carried by electrons are not much smaller than 1. More frequent observations at a fixed wavelength, or a wide spectrum at a fixed time, would put strong constraints on these parameters. We show that inverse-Compton scattering suppresses X-ray/ optical emission at delays t<T=10(F/0.3)^4 hr. Observations therefore imply F<0.3. The strong dependence of T on F implies that T may vary widely from burst to burst, and that frequent X-ray/optical observations at t~T would determine F. For F~0.2, inverse-Compton emission dominates during the first 2hr, producing >1GeV photons and providing a natural explanation to the delayed GeV emission observed in several strong bursts.

Observations suggest that $\gamma$-ray bursts (GRBs) are produced by the dissipation of the kinetic energy of a relativistic fireball. We show that a large fraction, $\ge 10%$, of the fireball energy is expected to be converted by photo-meson production to a burst of $\sim10^{14} eV$ neutrinos. A km^2 neutrino detector would observe at least several tens of events per year correlated with GRBs, and test for neutrino properties (e.g. flavor oscillations, for which upward moving $\tau$'s would be a unique signature, and coupling to gravity) with an accuracy many orders of magnitude better than is currently possible.

The cosmological gamma-ray burst (GRB) model for the production of ultra- high energy cosmic rays is described, and the current observational evidence which support it discussed. Several predictions of the model are presented, which would allow it to be tested by future high energy cosmic ray and gamma- ray experiments. If the predicted signatures of the GRB model are observed, they will not only corroborate the model, but will also provide information about the source population, and will allow to investigate the unknown structure of the inter-galactic magnetic field.

It has recently been shown that the highest energy cosmic rays (CRs) may originate in the same cosmological objects producing $\gamma$-ray bursts. This model requires the presence of intergalactic magnetic fields (IGMF) to delay the arrival times of $\sim 10^{20}$ eV CRs by 50 years or longer relative to the $\gamma$-rays, of an amplitude that is consistent with other observational constraints. Sources of CRs coming from individual bursts should be resolved with the planned ``Auger'' experiment, with as many as hundreds of CRs for the brightest sources. We analyze here the apparent angular and energy distribution of CRs from bright sources below the pion production threshold (in the energy range $10^{19}{\rm eV} < E < 4\times10^{19}{\rm eV}$) expected in this model. This observable distribution depends on the structure of the IGMF: the apparent spectral width $\Delta E$ is small, $\Delta E/E\lesssim1\%$, if the intergalactic field correlation length $\lambda$ is much larger than $1{\rm Mpc}$, and large, $\Delta E/E=0.3$, in the opposite limit $\lambda\ll 1{\rm Mpc}$. The apparent angular size is also larger for smaller $\lambda$. If the sources of CRs we predict are found, they will corroborate the bursting model and they will provide us with a technique to investigate the structure of the IGMF.

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

A scenario in which cosmic rays (CRs) above $10^{20}{\rm eV}$ are produced by cosmological gamma-ray bursts (GRBs) is consistent with observations provided that deflections by the inter-galactic magnetic field (IGMF) delay and spread the arrival time of the CRs over $\geq50{\rm yr}$. The energy lost by the CRs as they propagate and interact with the microwave background is transformed by cascading into secondary GeV-TeV photons. We show that a significant fraction of these photons can arrive with delays much smaller than the CR delay if much of inter-galactic space is occupied by large-scale magnetic ``voids'', regions of size $\gtrsim5{\rm Mpc}$ and field weaker than $10^{-15}{\rm G}$. Such voids might be expected, for example, in models where a weak primordial field is amplified in shocked, turbulent regions of the intergalactic medium during the formation of large-scale structure. For a field strength $\sim4\times10^{-11}{\rm G}$ in the high field regions, the value required to account for observed galactic fields if the IGMF were frozen in the protogalactic plasma, the delay of CRs produced by a burst at a distance of $100{\rm Mpc}$ is $\sim100{\rm yr}$, and the fluence of secondary photons above $10{\rm GeV}$ on hour--day time scales is $I(>E)\sim10^{-6}E_{\rm TeV}^{-1}{\rm cm}^{-2}$. This fluence is close to the detection threshold of current high-energy $\gamma$-ray experiments. Detection of the delayed flux would support the GRB-CR association and would also provide information on the IGMF structure.

We derive observational consequences of the hypothesis that cosmic rays (CR's) of energy $>10^{19}eV$ originate in the same cosmological objects producing gamma-ray bursts (GRB's). Inter-galactic magnetic fields $\gtrsim 10^{-12} G$ are required in this model to allow CR's to be observed continuously in time by producing energy dependent delays in the CR arrival times. This results in individual CR sources having very narrow observed spectra, since at any given time only those CR's having a fixed time delay are observed. Thus, the brightest CR sources should be different at different energies. The average number of sources contributing to the total CR flux decreases with energy much more rapidly than in a model of steady CR sources, dropping to one at $E_{crit} \simeq2\times10^{20}$~eV with very weak sensitivity to the inter-galactic magnetic field strength. Below $E_{crit}$, a very large number of sources is expected, consistent with observations. Above $E_{crit}$, a source may be observed with a flux considerably higher than the time-averaged CR flux from all sources, if a nearby GRB occurred recently. If such a source is present, its narrow spectrum may produce a ``gap'' in the total spectrum. These signatures should be detectable by the planned ``Auger'' CR experiment.

The X-ray emission by hot gas at the centers of clusters of galaxies is commonly modeled assuming the existence of steady-state, inhomogeneous cooling flows. We derive the metallicity profiles of the intracluster medium expected from such models. The inflowing gas is chemically enriched by type Ia supernovae and stellar mass loss in the outer parts of the central galaxy, which may give rise to a substantial metallicity gradient. The amplitude of the expected metallicity enhancement towards the cluster center is proportional to the ratio of the central galaxy luminosity to the mass inflow rate. The metallicity of the hotter phases is expected to be higher than that of the colder, denser phases. The metallicity profile expected for the Centaurus cluster is in good agreement with the metallicity gradient recently inferred from ASCA measurements (Fukazawa et al. 1994). However, current data do not rule out alternative models where cooling is balanced by some heat source. The metallicity gradient does not need to be present in all clusters, depending on the recent merging history of the gas around the central cluster galaxy, and on the ratio of the stellar mass in the central galaxy to the gas mass in the cooling flow.

The cosmic ray spectrum at $10^{19}{\rm eV}-10^{20}{\rm eV}$, reported by the Fly's Eye and the AGASA experiments, is shown to be consistent with a cosmological distribution of sources of protons, with a power law generation spectrum ${\rm d}\ln N/{\rm d}\ln E=-2.3\pm0.5$ and energy production rate of $4.5\pm1.5\times10^{44}{\rm erg}\ {\rm Mpc}^{-3}\ {\rm yr}^{-1}$. The two events measured above $10^{20}{\rm eV}$ are not inconsistent with this model. Verifying the existence of a ``black-body cutoff'', currently observed with low significance, would require $\sim30$ observation-years with existing experiments, but only $\sim1$ year with the proposed $\sim5000\ {\rm km}^2$ detectors. For a cosmological source distribution, no anisotropy is expected in the angular distribution of events with energies up to $\sim5\times10^{19}{\rm eV}$.

We discuss a scenario in which the highest energy cosmic rays (CR's) and cosmological $\gamma$-ray bursts (GRB's) have a common origin. This scenario is consistent with the observed CR flux above $10^{20}\text{eV}$, provided that each burst produces similar energies in $\gamma$-rays and in CR's above $10^{20}\text{eV}$. Protons may be accelerated by Fermi's mechanism to energies $\sim10^{20}\text{eV}$ in a dissipative, ultra-relativistic wind, with luminosity and Lorentz factor high enough to produce a GRB. For a homogeneous GRB distribution, this scenario predicts an isotropic, time-independent CR flux.

Observations of gravitational lensing indicate that the mass distribution in clusters of galaxies (where most of the mass is dark matter) is highly peaked towards the center, while X-ray observations imply that the gas is more extended. The additional fact that the gas is cooling in the center has often led one to expect that the gas temperature should be lower near the center, and therefore the gas should be more concentrated than the dark matter. We show that such expectation is not correct, and that the gas temperature must remain approximately constant within the cooling region in order to have consistency with the observed X-ray profiles and lenses. A multiphase cooling flow naturally produces an approximately constant temperature profile, and a more extended distribution for the gas compared to the mass. Cool phases are deposited at relatively large radius, while hot phases are adiabatically heated as they flow inwards and can keep the average temperature constant. Thus, cooling flows result in an {\it increase} of the central temperature, relative to a case where there is no cooling and the gas follows the mass distribution. The increased central temperatures caused by cooling flows give a characteristic core radius to the gas profiles, which is of order the cooling radius. This provides a natural explanation for the typical cores observed in X-ray clusters. It also brings into better agreement with observations the rate of cluster evolution expected in self-similar hierarchical models. We propose that clusters having core radii much larger than their cooling radii are in the process of merging and are not in dynamical equilibrium.

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