Mysteries about the origin of high-energy cosmic neutrinos have deepened by the recent IceCube measurement of a large diffuse flux in the 10-100 TeV range. Based on the standard disk-corona picture of active galactic nuclei (AGN), we present a phenomenological model enabling us to systematically calculate the spectral sequence of multimessenger emission from the AGN coronae. We show that protons in the coronal plasma can be stochastically accelerated up to PeV energies by plasma turbulence, and find that the model explains the large diffuse flux of medium-energy neutrinos if the cosmic rays carry only a few percent of the thermal energy. We find that the Bethe-Heitler process plays a crucial role in connecting these neutrinos and cascaded MeV gamma rays, and point out that the gamma-ray flux can even be enhanced by the reacceleration of secondary pairs. Critical tests of the model are given by its prediction that a significant fraction of the MeV gamma-ray background correlates with about 10 TeV neutrinos, and nearby Seyfert galaxies including NGC 1068 are promising targets for IceCube, KM3Net, IceCube-Gen2, and future MeV gamma-ray telescopes.

The LIGO discoveries have rekindled suggestions that primordial black holes (BHs) may constitute part to all of the dark matter (DM) in the Universe. Such suggestions came from 1) the observed merger rate of the BHs, 2) their unusual masses, 3) their low/zero spins, and 4) also from the independently uncovered cosmic infrared background (CIB) fluctuations signal of high amplitude and coherence with unresolved cosmic X-ray background (CXB). Here we summarize the prospects to resolve this important issue with electromagnetic observations using the instruments and tools expected in the 2020's. These prospects appear promising to make significant, and potentially critical, advances. We demonstrate that in the next decade, new space- and ground-borne electromagnetic instruments, combined with concurrent theoretical efforts, should shed critical light on the long-considered link between primordial BHs and DM. Specifically the new data and methodologies under this program will involve: I) Probing with high precision the spatial spectrum of source-subtracted CIB with Euclid and WFIRST, and its coherence with unresolved cosmic X-ray background using eROSITA and Athena, II) Advanced searches for microlensing of Galactic stars by the intervening Galactic Halo BHs with OGLE, Gaia, LSST and WFIRST, III) Supernovae (SNe) lensing in the upcoming surveys with WFIRST, LSST and also potentially with Euclid and JWST, IV) Advanced theoretical work to understand the details of PBH accretion and evolution and their influence on cosmic microwave background (CMB) anisotropies in light of the next generation CMB experiments, V) Better new samples and theoretical understanding involving stability and properties of ultra faint dwarf galaxies, pulsar timing, and cosmological quasar lensing.

The recent detection of gravitational waves and electromagnetic counterparts from the double neutron star merger event GW+EM170817, supports the standard paradigm of short gamma-ray bursts (SGRBs) and kilonovae/macronovae. It is important to reveal the nature of the compact remnant left after the merger, either a black hole or neutron star, and their physical link to the origin of the long-lasting emission observed in SGRBs. The diversity of the merger remnants may also lead to different kinds of transients that can be detected in future. Here we study the high-energy emission from the long-lasting central engine left after the coalescence, under certain assumptions. In particular, we consider the X-ray emission from a remnant disk and the non-thermal nebular emission from disk-driven outflows or pulsar winds. We demonstrate that late-time X-ray and high-frequency radio emission can provide useful constraints on properties of the hidden compact remnants and their connections to long-lasting SGRB emission, and we discuss the detectability of nearby merger events through late-time observations at $\sim30-100$ d after the coalescence. We also investigate the GeV-TeV gamma-ray emission that occurs in the presence of long-lasting central engines, and show the importance of external inverse-Compton radiation due to up-scattering of X-ray photons by relativistic electrons in the jet. We also search for high-energy gamma-rays from GW170817 in the Fermi-LAT data, and report upper limits on such long-lasting emission. Finally, we consider the implications of GW+EM170817 and discuss the constraints placed by X-ray and high-frequency radio observations.

High energy neutrino astrophysics has come of age with the discovery by IceCube of neutrinos in the TeV to PeV energy range attributable to extragalactic sources at cosmological distances. At such energies, astrophysical neutrinos must have their origin in cosmic ray interactions, providing information about the sources of high energy cosmic rays, as well as leading to the co-production of high energy gamma-rays. The intimate link with these latter two independently observed types of radiation provides important tools in the search for identifying and understanding the nature of their astrophysical sources. These neutrinons can provide important constraints about the cosmic ray acceleration process, and since they travel essentially unimpeded they can probe our Universe out to the farthest cosmological distances.

High-energy neutrinos are expected to originate from different stages in a gamma-ray burst (GRB) event. In this work we revisit the dissipative photospheric scenario, in which the GRB prompt emission is produced around the photospheric radius. Meanwhile, possible dissipation mechanisms (e.g., internal shocks or magnetic reconnection) could accelerate cosmic-rays (CRs) to ultra-high energies and then produce neutrinos via hadronuclear and photohadronic processes, which are referred to as prompt neutrinos. In this paper, we obtain the prompt neutrino spectrum of a single GRB within a self-consistent analytical framework, in which the jet-cocoon structure and possible collimation effects are included. We investigate a possible neutrino signal from the cocoon, which has been ignored in the previous studies. We show that if a GRB event happens at a distance of the order of Mpc, there is a great chance to observe the neutrino emission from the cocoon by IceCube, which is even more promising than jet neutrinos since the opening angle of the cocoon is much larger. We also figure out the diffuse neutrino flux of GRB cocoons and find that it could be comparable with that of the jets. Our results are consistent with the latest result reported by the IceCube collaboration that no significant correlation between neutrino events and observed GRBs is seen in the new data (Aartsen et al. 2017).

X-ray and $\gamma$-ray observations by the Swift satellite revealed that a fraction of tidal disruption events (TDEs) have relativistic jets. Jetted TDEs have been considered as potential sources of very high-energy cosmic-rays and neutrinos. In this work, using semi-analytical methods, we calculate neutrino spectra of X-ray bright TDEs with powerful jets and dark TDEs with possible choked jets, respectively. We estimate their neutrino fluxes and find that the non-detection would give us an upper limit on the cosmic-ray loading factor, $\xi_{\rm cr}\lesssim20$, for Sw J1644+57. We show that X-ray bright TDEs make a sub-dominant ($\lesssim5-10$%) contribution to IceCube's diffuse neutrino flux, and study possible contributions of X-ray dark TDEs given that particles are accelerated in choked jets or disk-winds. We discuss future prospects for multi-messenger searches of the brightest TDEs.

Tenuous wind bubbles, which are formed by the spin-down activity of central compact remnants, are relevant in some models of fast radio bursts (FRBs) and super-luminous supernovae. We study their high-energy signatures, focusing on the role of pair-enriched bubbles produced by young magnetars, rapidly-rotating neutron stars, and magnetized white dwarfs. (i) First, we study the nebular properties and the conditions allowing for escape of high-energy gamma-rays and radio waves, showing that their escape is possible for nebulae with ages of >10-100 yr. In the rapidly-rotating neutron star scenario, we find that radio emission from the quasi-steady nebula itself may be bright enough to be detected especially at sub-mm frequencies, which is relevant as a possible counterpart of pulsar-driven SNe and FRBs. (ii) Second, we consider the fate of bursting emission in the nebulae. We suggest that an impulsive burst may lead to a highly relativistic flow, which would interact with the nebula. If the shocked nebula is still relativistic, pre-existing non-thermal particles in the nebula can be significantly boosted by the forward shock, leading to short-duration (maybe millisecond or longer) high-energy gamma-ray flashes. Possible dissipation at the reverse shock may also lead to gamma-ray emission. (iii) After such flares, interactions with the baryonic ejecta may lead to afterglow emission with a duration of days to weeks. In the magnetar scenario, this burst-in-bubble model leads to the expectation that nearby (<10-100 Mpc) high-energy gamma-ray flashes may be detected by HAWC and CTA, and the subsequent afterglow emission may be seen by radio telescopes such as VLA. (iv) Finally, we discuss several implications specific to FRBs, including constraints on the emission regions and limits on soft gamma-ray counterparts.

We consider some general implications of bright gamma-ray counterparts to fast radio bursts (FRBs). We show that even if these manifest in only a fraction of FRBs, gamma-ray detections with current satellites (including Swift) can provide stringent constraints on cosmological FRB models. If the energy is drawn from the magnetic energy of a compact object such as a magnetized neutron star, the sources should be nearby and be very rare. If the intergalactic medium is responsible for the observed dispersion measure, the required gamma-ray energy is comparable to that of the early afterglow or extended emission of short gamma-ray bursts. While this can be reconciled with the rotation energy of compact objects, as expected in many merger scenarios, the prompt outflow that yields the gamma-rays is too dense for radio waves to escape. Highly relativistic winds launched in a precursor phase, and forming a wind bubble, may avoid the scattering and absorption limits and could yield FRB emission. Largely independent of source models, we show that detectable radio afterglow emission from gamma-ray bright FRBs can reasonably be anticipated. Gravitational wave searches can also be expected to provide useful tests.

The origin of prompt emission from gamma ray bursts remains to be an open question. Correlated prompt optical and gamma-ray emission observed in a handful of GRBs strongly suggests a common emission region, but failure to adequately fit the broadband GRB spectrum prompted the hypothesis of different emission mechanisms for the low- and high-energy radiations. We demonstrate that our multi-component model for GRB gamma-ray prompt emission provides an excellent fit to GRB 110205A from optical to gamma-ray energies. Our results show that the optical and highest gamma-ray emissions have the same spatial and spectral origin, which is different from the bulk of the X- and softest gamma-ray radiation. Finally, our accurate redshift estimate for GRB 110205A demonstrates promise for using GRBs as cosmological standard candles.

We present a targeted search for blazar flux-correlated high-energy ($\varepsilon_\nu > 1$ TeV) neutrinos from six bright northern blazars, using the public database of northern-hemisphere neutrinos detected during "IC40" 40-string operations of the IceCube neutrino observatory (April 2008 to May 2009). Our six targeted blazars are subjects of long-term monitoring campaigns by the VERITAS TeV gamma-ray observatory. We use the publicly-available VERITAS lightcurves to identify periods of excess and flaring emission. These predefined intervals serve as our "active temporal windows" in a search for an excess of neutrinos, relative to Poisson fluctuations of the near-isotropic atmospheric neutrino background which dominates at these energies. After defining the parameters of an optimized search, we confirm the expected Poisson behavior with Monte Carlo simulations prior to testing for excess neutrinos in the actual data. We make two searches: One for excess neutrinos associated with the bright flares of Mrk 421 that occurred during the IC40 run, and one for excess neutrinos associated with the brightest emission periods of five other blazars (Mrk 501, 1ES 0805+524, 1ES 1218+304, 3C66A, and W Comae), all significantly fainter than the Mrk 421 flares. We find no significant excess of neutrinos from the preselected blazar directions during the selected temporal windows. We derive 90%-confidence upper limits on the number of expected flux-associated neutrinos from each search. These limits are consistent with previous point-source searches and Fermi GeV flux-correlated searches. Our upper limits are sufficiently close to the physically-interesting regime that we anticipate future analyses using already-collected data will either constrain models or yield discovery of the first blazar-associated high-energy neutrinos.

We propose a novel model to produce ultrahigh-energy cosmic-rays (UHECRs) in gamma-ray burst jets. After the prompt gamma-ray emission, hydrodynamical turbulence is excited in the GRB jets at or before the afterglow phase. The mildly relativistic turbulence stochastically accelerates protons. The acceleration rate is much slower than the usual first-order shock acceleration rate, but in this case it can be energy-independent. The resultant UHECR spectrum is so hard that the bulk energy is concentrated in the highest energy range, resulting in a moderate requirement for the typical cosmic ray luminosity of $\sim 10^{53.5}~\mbox{erg}~\mbox{s}^{-1}$. In this model, the secondary gamma-ray and neutrino emissions initiated by photopion production are significantly suppressed. Although the UHECR spectrum at injection shows a curved feature, this does not conflict with the observed UHECR spectral shape. The cosmogenic neutrino spectrum in the $10^{17}$--$10^{18}$ eV range becomes distinctively hard in this model, which may be verified by future observations.

High-energy neutrino (HEN) and gravitational wave (GW) can probe astrophysical sources in addition to electromagnetic observations. Multimessenger studies can reveal nature of the sources which may not be discerned from one type of signal alone. We discuss HEN emission in connection with the Advanced Laser Interferometer Gravitational-wave Observatory (ALIGO) event GW150914 which could be associated with a short gamma-ray burst (GRB) detected by the $Fermi$ Gamma-ray Burst Monitor (GBM) 0.4 s after the GW event and within localization uncertainty of the GW event. We calculate HEN flux from this short GRB, GW150914-GBM, and show that non-detection of a high-energy starting event (HESE) by the IceCube Neutrino Observatory can constrain the total isotropic-equivalent jet energy of this short burst to be less than $3\times 10^{52}$ erg.

Recently, the direct detection of gravitational waves from black hole (BH) mergers was announced by the Advanced LIGO Collaboration. Multi-messenger counterparts of stellar-mass BH mergers are of interest, and it had been suggested that a small disk or celestial body may be involved in the binary of two BHs. To test such possibilities, we consider the fate of a wind powered by an active mini-disk in a relatively short, super-Eddington accretion episode onto a BH with ~10-100 solar masses. We show that its thermal emission could be seen as a fast optical transient with the duration from hours to days. We also find that the coasting outflow forms external shocks due to interaction with the interstellar medium, whose synchrotron emission might be expected in the radio band on a time scale of years. Finally, we also discuss a possible jet component and the associated high-energy neutrino emission as well as ultra-high-energy cosmic-ray acceleration.

We consider gamma-ray burst (GRB) jets that are choked by extended material as sources of high-energy cosmic neutrinos. We take into account the jet propagation physics both inside the progenitor star and the surrounding dense medium. Radiation constraints, which are relevant for high-energy neutrino production are considered as well. Efficient shock acceleration of cosmic rays is possible for sufficiently low-power jets and/or jets buried in a dense, extended wind or outer envelope. Such conditions also favor GRB jets to become stalled, and the necessary conditions for stalling are explicitly derived. Such choked jets may explain transrelativistic supernovae (SNe) and low-luminosity (LL) GRBs, giving a unified picture of GRBs and GRB-SNe. Focusing on this unified scenario for GRBs, we calculate the resulting neutrino spectra from choked jets including the relevant microphysical processes such as multipion production in pp and pgamma interactions, as well as the energy losses of mesons and muons. We obtain diffuse neutrino spectra using the latest results for the luminosity function of LL GRBs. Although uncertainties are large, we confirm that LL GRBs can potentially give a significant contribution to the diffuse neutrino flux. Our results are consistent with the present IceCube data and do not violate the stacking limits on classical high-luminosity GRBs. We find that high-energy neutrino production in choked jets is dominated by pgamma interactions. These sources are dark in GeV-TeV gamma-rays, and do not contribute significantly to the Fermi diffuse gamma-ray background. Assuming stalled jets can launch a quasi-spherical shock in the dense medium, precursor TeV neutrinos emerging prior to the shock breakout gamma-ray emission can be used as smoking gun evidence for a choked jet model for LL GRBs. Our results strengthen the relevance of wide field-of-view sky monitors.

The first gravitational-wave (GW) observations will greatly benefit from the detection of coincident electromagnetic counterparts. Electromagnetic follow-ups will nevertheless be challenging for GWs with poorly reconstructed directions. GW source localization can be inefficient (i) if only two GW observatories are in operation; (ii) if the detectors' sensitivities are highly non-uniform; (iii) for events near the detectors' horizon distance. For these events, follow-up observations will need to cover 100-1000 square degrees of the sky over a limited period of time, reducing the list of suitable telescopes. We demonstrate that the Cherenkov Telescope Array will be capable of following up GW event candidates over the required large sky area with sufficient sensitivity to detect short gamma-ray bursts, which are thought to originate from compact binary mergers, out to the horizon distance of advanced LIGO/Virgo. CTA can therefore be invaluable starting with the first multimessenger detections, even with poorly reconstructed GW source directions. This scenario also provides a further scientific incentive for GW observatories to further decrease the delay of their event reconstruction.

Recent observations by the $Fermi$ Gamma-ray Space Telescope have confirmed the existence of thermal and non-thermal components in the prompt photon spectra of some Gamma-ray bursts (GRBs). Through an analysis of six bright Fermi GRBs, we have discovered a correlation between the observed photospheric and non-thermal $\gamma$-ray emission components of several GRBs using a physical model that has previously been shown to be a good fit to the Fermi data. From the spectral parameters of these fits we find that the characteristic energies, $E_{\rm p}$ and $kT$, of these two components are correlated via the relation $E_{\rm p} \propto T^{\alpha}$ which varies from GRB to GRB. We present an interpretation in which the value of index $\alpha$ indicates whether the jet is dominated by kinetic or magnetic energy. To date, this jet composition parameter has been assumed in the modeling of GRB outflows rather than derived from the data.

We revisit the neutrino and ultra high-energy cosmic ray (UHECR) production from gamma-ray bursts (GRBs) with time-dependent simulations for the proton-induced cascades. This method can generate self-consistent photon, neutrino and escaped neutron spectra. To obtain the integrated background spectra, we take into account the distributions of the burst luminosity and pulse duration timescale. A benchmark case with standard GRB luminosity function, a bulk Lorentz factor $\Gamma=300$ and a proton to gamma-ray luminosity fraction $f_{\rm p}=10$, is consistent with both the neutrino upper-limits and the observed UHECR intensity at $\sim 10^{20}$ eV, while requiring a different type of UHECR source at the ankle. For the benchmark case the GRBs in the bright end of the luminosity function, which contribute most of the neutrinos, have their photon spectrum substantially distorted by secondary photons. Such bright GRBs are few in number, and reducing their $f_p$ eliminates the distortion, while reducing the neutrino production. Even if we neglect the contribution of the brightest GRBs, the UHECR production rate at GZK energies is almost unchanged. These nominal GRB models, especially with $L_{\rm iso} \lesssim 10^{53} ~\mbox{erg} ~\mbox{s}^{-1}$, appear to meet the current constraints as far as being candidate UHECR sources above the ankle energy.

High energy astrophysics is one of the most active branches in contemporary astrophysics. It studies astrophysical objects that emit X-ray and gamma-ray photons, such as accreting super-massive and stellar-size black holes, and various species of neutron stars. With the operations of many space-borne and ground-based observational facilities, high energy astrophysics has enjoyed a rapid development in the past decades. It is foreseen that the field will continue to advance rapidly in the coming decade, with possible ground-breaking discoveries of astrophysical sources in the high-energy neutrino and gravitational wave channels. This Special Issue of Frontiers of Physics is dedicated to a systematic survey of the field of high energy astrophysics as it stands in 2013.

Gamma-ray bursts are the most luminous explosions in the Universe. They appear connected to supernova remnants from massive stars or the merger of their remnants, and their brightness makes them temporarily detectable out to the larges distances yet explored in the Universe. After pioneering breakthroughs from space and ground experiments, their study is entering a new phase with observations from recently launched satellites, as well as the prospect of detections or limits from large neutrino and gravitational wave detectors. The interplay between such observations and theoretical models of gamma-ray bursts and related objects is reviewed.

In the internal shock model for gamma-ray bursts (GRBs), the synchrotron spectrum from the fast cooling electrons in a homogeneous downstream magnetic field (MF) is too soft to produce the low-energy slope of GRB spectra. However the magnetic field may decay downstream with distance from the shock front. Here we show that the synchrotron spectrum becomes harder if electrons undergo synchrotron and inverse-Compton cooling in a decaying MF. To reconcile this with the typical GRB spectrum with low energy slope $\nu F_\nu\propto\nu$, it is required that the postshock MF decay time is comparable to the cooling time of the bulk electrons (corresponding to a MF decaying length typically of $\sim10^5$ skin depths); that the inverse-Compton cooling should dominate synchrotron cooling after the MF decay time; and/or that the MF decays with comoving time roughly as $B\propto t^{-1.5}$. An internal shock synchrotron model with a decaying MF can account for the majority of GRBs with low energy slopes not harder than $\nu^{4/3}$.

Relativistic outflows with neutrons inevitably lead to inelastic collisions, and resulting subphotospheric gamma rays may explain prompt emission of gamma-ray bursts. In this model, hadronuclear, quasithermal neutrinos in the 10-100 GeV range should be generated, and they may even have a high-energy tail by neutron-proton-converter or shock acceleration mechanisms. We demonstrate the importance of dedicated searches with DeepCore+IceCube, though such analyses have not been performed. Successful detections enable us to discriminate among prompt emission mechanisms, probe the jet composition, and see roles of relativistic neutrons as well as effects of cosmic-ray acceleration.

The IceCube collaboration discovery of 28 high-energy neutrinos over the energy range 30 TeV <~ E_nu <~ 1 PeV, a 4.3-sigma excess over expected backgrounds, represents the first high-confidence detection of cosmic neutrinos at these energies. In light of this discovery, we explore the possibility that some of the Sub-PeV cosmic neutrinos might originate in our Galaxy's TeV unidentified (TeV UnID) sources. While typically resolved at TeV energies, these sources lack prominent radio or X-ray counterparts, and so have been considered promising sites for hadron acceleration within our Galaxy. Modeling the TeV UnID sources as Galactic hypernova remnants, we predict Sub-PeV neutrino fluxes and spectra consistent with their contributing a minority of n_nu <~ 2 of the observed events. This is consistent with our analysis of the spatial distribution of the Sub-PeV neutrinos and TeV UnID sources, which finds that a best-fit of one, and maximum of 3.8 (at 90%-confidence), of the ~16 non-atmospheric Sub-PeV neutrinos may originate in TeV UnID sources, with the remaining 75% to 95% of events being drawn from an isotropic background. If our scenario is correct, we expect excess Sub-PeV neutrinos to accumulate along the Galactic plane, within |l| <~ 30 deg of the Galactic center and in the Cygnus region, as observations by IceCube and other high-energy neutrino facilities go forward. Our scenario also has implications for radio, X-ray, and TeV observations of the TeV UnID sources.

The recent gamma-ray burst GRB 130427A has an isotropic electromagnetic energy $E^{iso}\sim10^{54}$ erg, suggesting an ample supply of target photons for photo-hadronic interactions, which at its low redshift of $z\sim0.34$ would appear to make it a promising candidate for neutrino detection. However, the IceCube collaboration has reported a null result based on a search during the prompt emission phase. We show that this neutrino non-detection can provide valuable information about this GRB's key physical parameters such as the emission radius $R_{d}$, the bulk Lorentz factor $\Gamma$ and the energy fraction converted into cosmic rays $\epsilon_{p}$. The results are discussed both in a model-independent way and in the specific scenarios of an internal shock model (IS), a baryonic photospheric model (BPH) and magnetic photospheric model (MPH). We find that the constraints are most stringent for the magnetic photospheric model considered, but the constraints on the internal shock and the baryonic photosphere models are fairly modest.

We report the first measurement of the atmospheric electron neutrino flux in the energy range between approximately 80 GeV and 6 TeV, using data recorded during the first year of operation of IceCube's DeepCore low energy extension. Techniques to identify neutrinos interacting within the DeepCore volume and veto muons originating outside the detector are demonstrated. A sample of 1029 events is observed in 281 days of data, of which 496 $\pm$ 66(stat.) $\pm$ 88(syst.) are estimated to be cascade events, including both electron neutrino and neutral current events. The rest of the sample includes residual backgrounds due to atmospheric muons and charged current interactions of atmospheric muon neutrinos. The flux of the atmospheric electron neutrinos is consistent with models of atmospheric neutrinos in this energy range. This constitutes the first observation of electron neutrinos and neutral current interactions in a very large volume neutrino telescope optimized for the TeV energy range.

We calculate the high energy neutrino spectrum from gamma-ray bursts where the emission arises in a dissipative jet photosphere determined by either baryonically or magnetically dominated dynamics, and compare these neutrino spectra to those obtained in conventional internal shock models. We also calculate the diffuse neutrino spectra based on these models, which appear compatible with the current IceCube 40+59 constraints. While a re-analysis based on the models discussed here and the data from the full array would be needed, it appears that only those models with the most extreme parameters are close to being constrained at present. A multi-year operation of the full IceCube and perhaps a next generation of large volume neutrino detectors may be required in order to distinguish between the various models discussed.

Gamma-ray bursts (GRBs) are bright flashes of gamma-rays coming from the cosmos. They occur roughly once per day, last typically 10s of seconds and are the most luminous events in the universe. More than three decades after their discovery, and after pioneering advances from space and ground experiments, they still remain mysterious. The launch of the Swift and Fermi satellites in 2004 and 2008 brought in a trove of qualitatively new data. In this review we survey the interplay between these recent observations and the theoretical models of the prompt GRB emission and the subsequent afterglows.

Heavy nuclei such as nickel-56 are synthesized in a wide range of core-collapse supernovae (CCSN), including energetic supernovae associated with gamma-ray bursts (GRBs). Recent studies suggest that jet-like outflows are a common feature of CCSN. These outflows may entrain synthesized nuclei at launch or during propagation, and provide interesting multi-messenger signals including heavy ultra-high energy cosmic rays. Here, we investigate the destruction processes of nuclei during crossing from the stellar material into the jet material via a cocoon, and during propagation after being successfully loaded into the jet. We find that nuclei can survive for a range of jet parameters because collisional cooling is faster than spallation. While canonical high-luminosity GRB jets may contain nuclei, magnetic dominated models or low-luminosity jets with small bulk Lorentz factors are more favorable for having a more significant heavy nuclei component.

We investigate the expected neutrino emissivity from nuclear collisions in magnetically dominated collisional models of gamma-ray bursts, motivated by recent observational and theoretical developments. The results indicate that significant multi-GeV neutrino fluxes are expected for model parameter values which are typical of electromagnetically detected bursts. We show that for detecting at least one muon event in Icecube and its Deep Core sub-array, a single burst must be near the high end of the luminosity function and at a redshift $z\lesssim 0.2$. We also calculate the luminosity and distance ranges that can generate $0.01-1$ muon events per GRB in the same detectors, which may be of interest if simultaneously detected electromagnetically, or if measured with future extensions of Icecube or other neutrino detectors with larger effective volume and better sensitivity.

The double burst, GRB 110709B, triggered Swift/BAT twice at 21:32:39 UT and 21:43:45 UT, respectively, on 9 July 2011. This is the first time we observed a GRB with two BAT triggers. In this paper, we present simultaneous Swift and Konus-WIND observations of this unusual GRB and its afterglow. If the two events originated from the same physical progenitor, their different time-dependent spectral evolution suggests they must belong to different episodes of the central engine, which may be a magnetar-to-BH accretion system.

We study effects of particle re-acceleration (or heating) in the post-shock region via magnetohydrodynamic/plasma turbulence, in the context of a mixed hadronic-leptonic model for the prompt emission of gamma-ray bursts (GRBs), using both analytical and numerical methods. We show that stochastically accelerated (or heated) leptons, which are injected via pp and pg reactions and subsequent pair cascades, are plausibly able to reproduce the Band function spectra with alpha~1 and beta~2-3 in the ~MeV range. An additional hard component coming from the proton-induced cascade emission is simultaneously expected, which is compatible with observed extra power-law spectra far above the MeV range. We also discuss the specific implications of hadronic models for ongoing high-energy neutrino observations.

We present calculations of the time evolution of the prompt spectra of gamma-ray burst models involving generic internal dissipation regions, including internal shocks, either by itself or in the presence of an external photon source such as a photosphere. The method uses a newly developed time-dependent code involving synchrotron emission and absorption, inverse Compton scattering and pair formation. The models reproduce the typical observed Band spectra and their generic time evolution, including the appearance of an extra keV-GeV component, whose delay in simple SSC models, however, is only partially able to explain the several seconds observed GeV delays. On the other hand, models involving both a photosphere and an internal dissipation region at a larger radius produce both an extra GeV component and time delays which are in the range of the observations.

The origin of the shallow-decay emission during early X-ray afterglows has been an open issue since the launch of the Swift satellite. One of the appealing models is the late internal dissipation model, where X-ray emission during the shallow-decay phase is attributed to internal dissipation, analogous to the prompt gamma-ray emission. We discuss possible scenarios of the late prompt emission, such as late internal shocks, magnetic reconnection, and photospheric dissipation. We also consider the consequences of late dissipation and a two-component (early and late) jet model for the high-energy (GeV-TeV) emission. We study not only synchrotron self-Compton (SSC) emission from the early and late jets but also external inverse-Compton (EIC) emission, which is naturally predicted in the late dissipation model. For the latter, we perform numerical calculations taking into account the equal-arrival-time surface of EIC photons and show that the EIC component typically has a peak at ~1-100 GeV which may dominate over the SSC components. We demonstrate that very high energy gamma rays from both these components are detectable for nearby and/or energetic gamma-ray bursts, with current and future Cherenkov detectors such as MAGIC, VERITAS, CTA and HAWC, and possibly Fermi. Although the expected event rate would not be large, detections should be useful as a test of the model. Multi-wavelength observations using both the ground-based telescopes and the Swift and/or Fermi satellites are also important to constrain the models.

We discuss the high energy neutrino emission from gamma-ray bursts resulting from the earliest generation (`population III') stars forming in the Universe, whose core collapses into a black hole. These gamma-ray bursts are expected to produce a highly relativistic, magnetically dominated jet, where protons can be accelerated to ultra-high energies. These interact with the photons produced by the jet, leading to ultra-high energy photo-meson neutrinos as well as secondary leptons and photons. The photon luminosity and the shock properties, and thus the neutrino spectrum, depend on the mass of the black holes as well as on the density of the surrounding external gas. We calculate the individual source neutrino spectral fluxes and the expected diffuse neutrino flux for various source parameters and evolution scenarios. Both the individual and diffuse signals appear detectable in the 1-300 PeV range with current and planned neutrino detectors such as IceCube and ARIANNA, provided the black hole mass is in excess of 30-100 solar masses. This provides a possible test for the debated mass of the progenitor stellar objects, as well as a probe for the early cosmological environment and the formation rate of the earliest structures.

Radially inhomogeneous gamma-ray burst (GRB) jets release variable photospheric emission and can have internal shocks occurring above the photosphere. We generically formulate a photospheric emission model of GRBs including Compton up-scattered photospheric (UP) emission off the electrons (and positrons) in the internal shocks, and find that the photospheric emission may correspond to the traditional (Band) component at <~1 MeV and the UP emission to the high-energy emission observed by Fermi/LAT for some GRBs at >~ 10 MeV. The two components can be separate in the spectrum in some cases or can mimic a smooth broad Band spectrum in other cases. We apply our formulation to the well-studied long and short LAT GRBs, GRB 080916C, GRB 090902B, and GRB 090510, and typically find reasonable parameters for fitting the time-binned spectra, although fine tuning of several parameters is required. The observed delays of the high-energy emission with respect to the MeV emission which are large compared to the variability times are unlikely to be due to simple kinematic effects of a non-evolving jet. These delays may be attributed to the temporal evolution of the physical parameters of the jet, and thus the delay timescales could provide a potential tool for investigating the structures of GRB jets themselves and their progenitors. The difference of the delay timescales of long and short GRBs inferred from the Fermi data might be due to the differences in the progenitors of long and short GRBs. Some other properties and consequences of this model are discussed, including temporal correlations among the prompt optical, the soft X-ray, and the distinct high-energy component as well as the Band component.

Long-lived high-energy (>100MeV) emission, a common feature of most Fermi-LAT detected gamma-ray burst, is detected up to \sim 10^2 s in the short GRB 090510. We study the origin of this long-lived high-energy emission, using broad-band observations including X-ray and optical data. We confirm that the late > 100 MeV, X-ray and optical emission can be naturally explained via synchrotron emission from an adiabatic forward shock propagating into a homogeneous ambient medium with low number density. The Klein-Nishina effects are found to be significant, and effects due to jet spreading and magnetic field amplification in the shock appear to be required. Under the constraints from the low-energy observations, the adiabatic forward shock synchrotron emission is consistent with the later-time (t>2s) high-energy emission, but falls below the early-time (t < 2s) high energy emission. Thus we argue that an extra high energy component is needed at early times. A standard reverse shock origin is found to be inconsistent with this extra component. Therefore, we attribute the early part of the high-energy emission (t< 2s) to the prompt component, and the long-lived high energy emission (t>2s) to the adiabatic forward shock synchrotron afterglow radiation. This avoids the requirement for an extremely high initial Lorentz factor.

Population III stars are theoretically expected to be prominent around redshifts z ~ 20, consisting of mainly very massive stars with M_* >~ 10 M_sun$, but there is no direct observational evidence for these objects. They may produce collapsar gamma-ray bursts (GRBs), with jets driven by magnetohydrodynamic processes, whose total isotropic-equivalent energy could be as high as E_iso >~ 10^{57} erg over a cosmological-rest-frame duration of t_d >~ 10^4 s, depending on the progenitor mass. Here we calculate the afterglow spectra of such Pop. III GRBs based on the standard external shock model, and show that they will be detectable with the Swift BAT/XRT and Fermi LAT instruments. We find that in some cases a spectral break due to electron-positron pair creation will be observable in the LAT energy range, which can put constraints on the ambient density of the pre-collapse Pop. III star. Thus, high redshift GRB afterglow observations could be unique and powerful probes of the properties of Pop. III stars and their environments. We examine the trigger threshold of the BAT instrument in detail, focusing on the image trigger system, and show that the prompt emission of Pop. III GRBs could also be detected by BAT. Finally we briefly show that the late-time radio afterglows of Pop. III GRBs for typical parameters, despite the large distances, can be very bright: ~ 140 mJy at 1 GHz, which may lead to a constraint on the Pop. III GRB rate from the current radio survey data, and ~2.4 mJy at 70 MHz, which implies that Pop. III GRB radio afterglows could be interesting background source candidates for 21 cm absorption line detections.

A fraction of gamma-ray bursts exhibit distinct spectral features in their prompt emission below few 10s of keV that exceed simple extrapolations of the low-energy power-law portion of the Band spectral model. This is also true for the prompt optical emission observed in several bursts. Through Monte Carlo simulations, we model such low-energy spectral excess components as hadronic cascade emission initiated by photomeson interactions of ultra-high-energy protons accelerated within GRB outflows. Synchrotron radiation from the cascading, secondary electron-positron pairs can naturally reproduce the observed soft spectra in the X-ray band, and in some cases the optical spectra as well. These components can be directly related to the higher energy radiation at GeV energies due to the hadronic cascades. Depending on the spectral shape, the total energy in protons is required to be comparable to or appreciably larger than the observed total photon energy. In particular, we apply our model to the excess X-ray and GeV emission of GRB 090902B, and the bright optical emission of the "naked-eye" GRB 080319B. Besides the hard GeV components detected by {\it Fermi}, such X-ray or optical spectral excesses are further potential signatures of ultra-high-energy cosmic ray production in gamma-ray bursts.

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.

Gamma-ray bursts have been proved to be detectable up to distances much larger than any other astrophysical object, providing the most effective way, complementary to ordinary surveys, to study the high redshift universe. To this end, we present here the results of an observational campaign devoted to the study of the high-z GRB 090205. We carried out optical/NIR spectroscopy and imaging of GRB 090205 with the ESO-VLT starting from hours after the event up to several days later to detect the host galaxy. We compared the results obtained from our optical/NIR observations with the available Swift high-energy data of this burst. Our observational campaign led to the detection of the optical afterglow and host galaxy of GRB 090205 and to the first measure of its redshift, z=4.65. Similar to other, recent high-z GRBs, GRB 090205 has a short duration in the rest-frame with T_{90,rf}=1.6 s, which suggests the possibility that it might belong to the short GRBs class. The X-ray afterglow of GRB 090205 shows a complex and interesting behaviour with a possible rebrightening at 500-1000s from the trigger time and late flaring activity. Photometric observations of the GRB 090205 host galaxy argue in favor of a starburst galaxy with a stellar population younger than ~ 150 Myr. Moreover, the metallicity of Z > 0.27 Z_Sun derived from the GRB afterglow spectrum is among the highest derived from GRB afterglow measurement at high-z, suggesting that the burst occurred in a rather enriched envirorment. Finally, a detailed analysis of the afterglow spectrum shows the existence of a line corresponding to Lyman-alpha emission at the redshift of the burst. GRB 090205 is thus hosted in a typical Lyman-alpha emitter (LAE) at z=4.65. This makes the GRB 090205 host the farthest GRB host galaxy, spectroscopically confirmed, detected to date.

IceCube Collaboration Contributions to the 2009 International Cosmic Ray Conference

We discuss a model of Poynting-dominated gamma-ray bursts from the collapse of very massive first generation (pop. III) stars. From redshifts of order 20, the resulting relativistic jets would radiate in the hard X-ray range around 50 keV and above, followed after roughly a day by an external shock component peaking around a few keV. On the same timescales an inverse Compton component around 75 GeV may be expected, as well as a possible infra-red flash. The fluences of these components would be above the threshold for detectors such as Swift and Fermi, providing potentially valuable information on the formation and properties of what may be the first luminous objects and their black holes in the high redshift Universe.

The Fermi satellite has been reporting the detailed temporal properties of gamma-ray bursts (GRBs) in an extremely broad spectral range, 8 keV - 300 GeV, in particular, the unexpected delays of the GeV emission onsets behind the MeV emission of some GRBs. We focus on GRB 080916C, one of the Fermi-LAT GRBs for which the data of the delayed high-energy emission are quite extensive, and we show that the behavior of the high-energy emission of this burst can be explained by a model in which the prompt emission consists of two components: one is the MeV component due to the synchrotron-self-Compton radiation of electrons accelerated in the internal shock of the jet and the other is the high-energy component due to inverse Compton scattering of the photospheric X-ray emission of the expanding cocoon off the same electrons in the jet. Such an external inverse Compton effect could be important for other Fermi-LAT GRBs, including short GRBs as well. In this model, the delay timescale is directly linked to the physical properties of GRB progenitor.

We investigate hypernova (hyper-energetic supernova) and gamma-ray burst (GRB) remnants in our Galaxy as TeV gamma-ray sources, particularly in the role of potential TeV unidentified sources, which have no clear counterpart at other wavelengths. We show that the observed bright sources in the TeV sky could be dominated by GRB/hypernova remnants, even though they are fewer than supernova remnants (SNRs). If this is the case, TeV SNRs are more extended (and more numerous) than deduced from current observations. In keeping with their role as cosmic ray accelerators, we discuss hadronic gamma-ray emission from pi^0 decay, from beta decay followed by inverse Compton emission, and propose a third, novel process of TeV gamma-ray emission arising from the decay of accelerated radioactive isotopes such as 56Co entrained by relativistic or semi-relativistic jets in GRBs/hypernovae. We discuss the relevant observational signatures which could discriminate between these three mechanisms.

The observations of gamma-ray bursts (GRBs) such as 980425, 031203 and 060218, with luminosities much lower than those of other classic bursts, lead to the definition of a new class of GRBs -- low-luminosity GRBs. The nature of the outflow responsible for them is not clear yet. Two scenarios have been suggested: one is the conventional relativistic outflow with initial Lorentz factor of order of $\Gamma_0\ga 10$ and the other is a trans-relativistic outflow with $\Gamma_0\simeq 1-2$. Here we compare the high energy gamma-ray afterglow emission from these two different models, taking into account both synchrotron self inverse-Compton scattering (SSC) and the external inverse-Compton scattering due to photons from the cooling supernova or hypernova envelope (SNIC). We find that the conventional relativistic outflow model predicts a relatively high gamma-ray flux from SSC at early times ($<10^4 {\rm s}$ for typical parameters) with a rapidly decaying light curve, while in the trans-relativistic outflow model, one would expect a much flatter light curve of high-energy gamma-ray emission at early times, which could be dominated by both the SSC emission and SNIC emission, depending on the properties of the underlying supernova and the shock parameter $\epsilon_e$ and $\epsilon_B$. The Fermi Gamma-ray Space Telescope should be able to distinguish between the two models in the future.

The Fermi Gamma-ray Space Telescope recently detected the most energetic gamma-ray burst so far, GRB 080916C, and reported its detailed temporal properties in an extremely broad spectral range: (i) the time-resolved spectra are well described by broken power-law forms over the energy range of $10 {\rm keV}-10$ GeV, (ii) the high-energy emission (at $\epsilon > 100$ MeV) is delayed by $\approx 5$s with respect to the $\epsilon \lesssim 1$ MeV emission, and (iii) the emission onset times shift towards later times in the higher energy bands. We show that this behavior of the high-energy emission can be explained by a model in which the prompt emission consists of two components: one is the emission component peaking at $\epsilon \sim 1$ MeV due to the synchrotron-self-Compton radiation of electrons accelerated in the internal shock of the jet and the other is the component peaking at $\epsilon \sim 100$ MeV due to up-scattering of the photospheric X-ray emission of the expanding cocoon (i.e., the hot bubble produced by dissipation of the jet energy inside the progenitor star) off the same electrons in the jet. Based on this model, we derive some constraints on the radius of the progenitor star and the total energy and mass of the cocoon of this GRB, which may provide information on the structure of the progenitor star and the physical conditions of the jet propagating in the star. The up-scattered cocoon emission could be important for other Fermi-GRBs as well. We discuss some predictions of this model, including a prompt bright optical emission and a soft X-ray excess.

The existence of a shallow decay phase in the early X-ray afterglows of gamma-ray bursts is a common feature. Here we investigate the possibility that this is connected to the formation of a highly magnetized millisecond pulsar, pumping energy into the fireball on timescales longer than the prompt emission. In this scenario the nascent neutron star could undergo a secular bar-mode instability, leading to gravitational wave losses which would affect the neutron star spin-down. In this case, nearby gamma-ray bursts with isotropic energies of the order of 1e50 ergs would produce a detectable gravitational wave signal emitted in association with an observed X-ray light-curve plateau, over relatively long timescales of minutes to about an hour. The peak amplitude of the gravitational wave signal would be delayed with respect to the gamma-ray burst trigger, offering gravitational wave interferometers such as the advanced LIGO and Virgo the challenging possibility of catching its signature on the fly.

(Abridged) The two high-redshift gamma-ray bursts, GRB 080913 at z=6.7 and GRB 090423 at z=8.3, recently detected by Swift appear as intrinsically short, hard GRBs. They could have been recognized by BATSE as short/hard GRBs should they have occurred at z <= 1. We perform a more thorough investigation on two physically distinct types (Type I/II) of cosmological GRBs and their observational characteristics. We reiterate the definitions of Type I/II GRBs and review the observational criteria and their physical motivations. Contrary to the traditional approach of assigning the physical category based on the gamma-ray properties (duration, hardness, and spectral lag), we take an alternative approach to define the Type I and Type II Gold Samples using several criteria that are more directly related to the GRB progenitors, and study the properties of the two Gold Samples and compare them with the traditional long/soft and short/hard samples. We find that the Type II Gold Sample reasonably tracks the long/soft population, although it includes several intrinsically short (shorter than 1s in the rest frame) GRBs. The Type I Gold Sample only has 5 GRBs, 4 of which are not strictly short but have extended emission. Other short/hard GRBs detected in the Swift era represent the BATSE short/hard sample well, but it is unclear whether all of them belong to Type I. We suggest that some (probably even most) high-luminosity short/hard GRBs instead belong to Type II. We suggest that GRB 080913 and GRB 090423 are more likely Type II events. We re-emphasize the importance of invoking multiple observational criteria, and cautiously propose an operational procedure to infer the physical origin of a given GRB with available multiple observational criteria, with various caveats laid out.

Fermi observations of high-energy gamma-ray emission from GRB 080916C shows that its spectrum is consistent with the Band function from MeV to tens of GeV. Assuming one single emission mechanism dominates in the whole energy range, we show that this spectrum is consistent with synchrotron origin by shock-accelerated electrons. The simple electron inverse-Compton model and the hadronic model are found to be less viable. In the synchrotron scenario, the synchrotron self-Compton scattering is likely to be in the Klein-Nishina regime and therefore the resulting high-energy emission is subdominant, even though the magnetic field energy density is lower than that in relativistic electrons. The Klein-Nishina inverse-Compton cooling may also affect the low-energy electron number distribution and hence results in a low-energy synchrotron photon spectrum $n(\nu)\propto\nu^{-1}$ below the peak energy. Under the framework of the electron synchrotron interpretation, we constrain the shock microphysical parameters and derive a lower limit of the upstream magnetic fields. The detection of synchrotron emission extending to about 70 GeV in the source frame in GRB 080916C favors the Bohm diffusive shock acceleration if the bulk Lorentz factor of the relativistic outflow is not significantly greater than thousands.

We investigate the high-energy neutrino emission expected from newly born magnetars surrounded by their stellar ejecta. Protons might be accelerated up to 0.1-100 EeV energies possibly by, e.g., the wave dissipation in the winds, leading to hadronic interactions in the stellar ejecta. The resulting PeV-EeV neutrinos can be detected by IceCube/KM3Net with a typical peak time scale of a few days after the birth of magnetars, making the characteristic soft-hard-soft behavior. Detections would be important as a clue to the formation mechanism of magnetars, although there are ambiguities coming from uncertainties of several parameters such as velocity of the ejecta. Non-detections would also lead to useful constraints on the scenario.

Gamma-ray bursts (GRBs) have been an enigma since their discovery forty years ago. However, considerable progress unraveling their mysteries has been made in recent years. Developments in observations, theory, and instrumentation have prepared the way so that the next decade can be the one in which we finally answer the question, "What are gamma-ray bursts?" This question encompasses not only what the progenitors are that produce the GRBs, but also how the enormous luminosity of the GRBs, concentrated in gamma rays, is achieved. Observations across the electromagnetic spectrum, from both the ground and space, will be required to fully tackle this important question. This white paper, mostly distilled from a recent study commissioned by the Division of Astrophysics of the American Physical Society, focuses on what very high energy (~100 GeV and above) gamma-ray observations can contribute. Very high energy gamma rays probe the most extreme high energy particle populations in the burst environment, testing models of lepton and proton acceleration in GRBs and constraining the bulk Lorentz factor and opacity of the outflow. Sensitivity improvements of more than an order of magnitude in the very high energy gamma-ray band can be achieved early in the next decade, in order to contribute to this science.

Gamma-ray Bursts (GRBs) are relativistic cosmological beacons of transient high energy radiation whose afterglows span the electromagnetic spectrum. Theoretical expectations of correlated neutrino and/or gravitational wave (GW) emission position GRBs at an astrophysical nexus for a metamorphosis in our understanding of the Cosmos. This new dawn in the era of experimental (particle) astrophysics and cosmology is afforded by current and planned facilities enabling the novel astronomies of high energy neutrinos and gravitational waves, in concert with unprecedented electromagnetic coverage. In this white paper, we motivate GRBs as a compelling scientific theme and highlight key technical advances that may facilitate fundamental breakthroughs in the context of Swift, Fermi, IceCube and LIGO (Laser Interferometer Gravitational Wave Observatory), whose capabilities would be augmented with JANUS (Joint Astrophysics Nascent Universe Satellite), EXIST (Energetic X-ray Imaging Survey Telescope) and LISA (Laser Interferometer Space Antenna). Scientific synergy will be achieved by leveraging the combined sensitivity of contemporaneous ground-based and satellite observatories, thus optimizing their collective discovery potential for: (i) revealing the origin(s) and acceleration mechanism(s) of cosmic rays, (ii) exposing GRB progenitor(s) and (iii) exploring the high-z Cosmos. Hence, the advent of GRB multi-messenger astronomy may cement an explicit connection to fundamental physics, via nascent cosmic windows, throughout the next decade and beyond.

We report on observations of correlated behavior between the prompt gamma-ray and optical emission from GRB 080319B, which confirm that (i) they occurred within the same astrophysical source region and (ii) their respective radiation mechanisms were dynamically coupled. Our results, based upon a new CCF methodology for determining the time-resolved spectral lag, are summarized as follows. First, the evolution in the arrival offset of prompt gamma-ray photon counts between Swift-BAT 15-25 keV and 50-100 keV energy bands (intrinsic gamma-ray spectral lag) appears to be anti-correlated with the arrival offset between prompt 15-350 keV gamma-rays and the optical emission observed by TORTORA (extrinsic optical/gamma-ray spectral lag), thus effectively partitioning the burst into two main episodes at ~T+28+/-2 sec. Second, the rise and decline of prompt optical emission at ~T+10+/-1 sec and ~T+50+/-1 sec, respectively, both coincide with discontinuities in the hard to soft evolution of the photon index for a power law fit to 15-150 keV Swift-BAT data at ~T+8+/-2 sec and ~T+48+/-1 sec. These spectral energy changes also coincide with intervals whose time-resolved spectral lag values are consistent with zero, at ~T+12+/-2 sec and ~T+50+/-2 sec. These results, which are robust across heuristic permutations of Swift-BAT energy channels and varying temporal bin resolution, have also been corroborated via independent analysis of Konus-Wind data. This potential discovery may provide the first observational evidence for an implicit connection between spectral lags and GRB emission mechanisms in the context of canonical fireball phenomenology. Future work includes exploring a subset of bursts with prompt optical emission to probe the unique or ubiquitous nature of this result.

We report on the detection by Swift of GRB 080913, and subsequent optical/near-infrared follow-up observations by GROND which led to the discovery of its optical/NIR afterglow and the recognition of its high-z nature via the detection of a spectral break between the i' and z' bands. Spectroscopy obtained at the ESO-VLT revealed a continuum extending down to lambda = 9400 A, and zero flux for 7500 A < lambda<9400 A, which we interpret as the onset of a Gunn-Peterson trough at z=6.695+-0.025 (95.5% conf. level), making GRB 080913 the highest redshift GRB to date, and more distant than the highest-redshift QSO. We note that many redshift indicators which are based on promptly available burst or afterglow properties have failed for GRB 080913. We report on our follow-up campaign and compare the properties of GRB 080913 with bursts at lower redshift. In particular, since the afterglow of this burst is fainter than typical for GRBs, we show that 2m-class telescopes can identify most high-redshift GRBs.

This is a report on the findings of the SNR/cosmic-ray working group for the white paper on the status and future of ground-based gamma-ray astronomy. The white paper is an APS commissioned document, and the overall version has also been released and can be found on astro-ph. This detailed section of the white paper discusses the status of past and current attempts to observe shell-type supernova remnants and diffuse emission from cosmic rays at GeV-TeV energies. We concentrate on the potential of future ground-based gamma-ray experiments to study the acceleration of relativistic charged particles which is one of the main unsolved, yet fundamental, problems in modern astrophysics. The acceleration of particles relies on interactions between energetic particles and magnetic turbulence. In the case of SNRs we can perform spatially resolved studies in systems with known geometry, and the plasma physics deduced from these observations will help us to understand other systems where rapid particle acceleration is believed to occur and where observations as detailed as those of SNRs are not possible.

This is a report on the findings of the gamma ray burst working group for the white paper on the status and future of TeV gamma-ray astronomy. The white paper is an APS commissioned document, and the overall version has also been released and can be found on astro-ph. This detailed section of the white paper discusses the status of past and current attempts to observe gamma ray bursts at GeV-TeV energies. We concentrate on the potential of future ground-based gamma-ray experiments to observe the highest energy emission ever recorded for GRBs, particularly for those that are nearby and have high Lorentz factors in the GRB jet. It is clear that major advances are possible and that the detection of very high energy emission would have strong implications for GRB models, as well as cosmic ray origin.

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.

A class of very energetic supernovae (hypernovae) is associated with long gamma-ray bursts, in particular with a less energetic but more frequent population of gamma-ray bursts. Hypernovae also appear to be associated with mildly relativistic jets or outflows, even in the absence of gamma-ray bursts. Here we consider radiation from charged particles accelerated in such mildly relativistic outflows with kinetic energies of ~10^{50} erg. The radiation processes of the primarily accelerated electrons considered are synchrotron radiation and inverse-Compton scattering of synchrotron photons (synchrotron self-Compton; SSC) and of supernova photons (external inverse-Compton; EIC). In the soft X-ray regime, both the SSC and EIC flux can be the dominant component, but due to their very different spectral shapes it should be easy to distinguish between them. When the fraction of the kinetic energy going into the electrons (\epsilon_e) is large, the SSC is expected to be important; otherwise the EIC will dominate. The EIC flux is quite high, almost independently of \epsilon_e, providing a good target for X-ray telescopes such as XMM-Newton and Chandra. In the GeV gamma-ray regime, the EIC would be the dominant radiation process and the Gamma-ray Large Area Space Telescope (GLAST) should be able to probe the value of \epsilon_e, the spectrum of the electrons, and their maximum acceleration energy. Accelerated protons also lead to photon radiation through the secondary electrons produced by the photopion and photopair processes. We find that over a significant range of parameters the proton component is generally less prominent than the primary electron component. We discuss the prospects for the detection of the X-ray and GeV signatures of the mildly relativistic outflow of hypernovae.

We discuss the relevance of the curvature radiation of pions in strongly magnetized pulsars or magnetars, and their implications for the production of TeV energy neutrinos detectable by cubic kilometer scale detectors, as well as high energy photons.

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

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.

This is a short report on the preliminary findings of the gamma ray burst (GRB) working group for the white paper on the status and future of very high energy (VHE; >50 GeV) gamma-ray astronomy. The white paper discusses the status of past and current attempts to observe GRBs at GeV-TeV energies, including a handful of low-significance, possible detections. The white paper concentrates on the potential of future ground-based gamma-ray experiments to observe the highest energy emission ever recorded for GRBs, particularly for those that are nearby and have high Lorentz factors in the GRB jet. It is clear that the detection of VHE emission would have strong implications for GRB models, as well as cosmic ray origin. In particular, the extended emission phase (including both afterglow emission and possible flaring) of nearby long GRBs could provide the best possibility for detection. The difficult-to-obtain observations during the prompt phase of nearby long GRBs and short GRBs could also provide particularly strong constraints on the opacity and bulk Lorentz factors surrounding the acceleration site. The synergy with upcoming and existing observatories will, of course, be critical for both identification of GRBs and for multiwavelength/multimessenger studies.

The chemical composition of the ultra-high-energy (UHE) cosmic rays serves as an important clue for their origin. Recent measurements of the elongation rates by the Pierre Auger Observatory hint at the possible presence of heavy or intermediate mass nuclei in the UHE cosmic rays. Gamma-ray bursts (GRBs) and hypernovae have been suggested as possible sources of the UHE cosmic rays. Here we derive the constraints on the physical conditions under which UHE heavy nuclei, if they are accelerated in these sources, can survive in their intense photon fields. We find that in the GRB external shock and in the hypernova scenarios, UHE nuclei can easily survive photo-disintegration. In the GRB internal shock scenario, UHE nuclei can also survive, provided the dissipation radius and/or the bulk Lorentz factor of the relativistic outflow are relatively large, or if the low-energy self-absorption break in the photon spectrum of the prompt emission occurs above several KeV. In internal shocks and in the other scenarios, intermediate-mass UHE nuclei have a higher probability of survival against photo-disintegration than UHE heavy nuclei such as Fe.

The early optical emission of the moderately high redshift ($z=3.08$) GRB 060607A shows a remarkable broad and strong peak with a rapid rise and a relatively slow power-law decay. It is not coincident with the strong early-time flares seen in the X-ray and gamma-ray energy bands. There is weak evidence for variability superposed on this dominant component in several optical bands that can be related to flares in high energy bands. While for a small number of GRBs, well-sampled optical flares have been observed simultaneously with X-ray and gamma ray pulses, GRB 060607A is one of the few cases where the early optical emission shows no significant evidence for correlation with the prompt emission. In this work we first report in detail the broad band observations of this burst by Swift. Then by applying a simple model for the dynamics and the synchrotron radiation of a relativistic shock, we show that the dominant component of the early emissions in optical wavelengths has the same origin as the tail emission produced after the main gamma ray activity. The most plausible explanation for the peak in the optical light curve seems to be the cooling of the prompt after the main collisions, shifting the characteristic synchrotron frequency to the optical bands. It seems that the cooling process requires a steepening of the electron energy distribution and/or a break in this distribution at high energies. The sharp break in the X-ray light curve at few thousands of seconds after the trigger, is not observed in the IR/optical/UV bands, and therefore can not be a jet break. Either the X-ray break is due to a change in the spectrum of the accelerated electrons or the lack of an optical break is due to the presence of a related delayed response component (Abbreviated).

The origin of the ultrahigh-energy (UHE) cosmic rays (CRs) from the second knee ($\sim6\times10^{17}$eV) above in the CR spectrum is still unknown. Recently, there has been growing evidence that a peculiar type of supernovae, called hypernovae, are associated with sub-energetic gamma-ray bursts (GRBs), such as SN1998bw/GRB980425 and SN2003lw/GRB031203. Such hypernovae appear to have high (up to mildly relativistic) velocity ejecta, which may be linked to the sub-energetic GRBs. Assuming a continuous distribution of the kinetic energy of the hypernova ejecta as a function of its velocity $E_k\propto (\Gamma\beta)^{-\alpha}$ with $\alpha\sim 2$, we find that 1) the external shock wave produced by the high velocity ejecta of a hypernova can accelerate protons up to energies as high as $10^{19} {\rm eV}$; 2) the cosmological hypernova rate is sufficient to account for the energy flux above the second knee; and 3) the steeper spectrum of CRs at these energies can arise in these sources. In addition, hypernovae would also give rise to a faint diffuse UHE neutrino flux, due to $p\gamma$ interactions of the UHE CRs with hypernova optical-UV photons.

Precursor emission has been observed in a non-negligible fraction of gamma-ray bursts.The time gap between the precursor and the main burst extends in some case up to hundreds of seconds, such as in GRB041219A, GRB050820A and GRB060124. Both the origin of the precursor and the large value of the time gap are controversial. Here we investigate the maximum possible time gaps arising from the jet propagation inside the progenitor star, in models which assume that the precursor is produced by the jet bow shock or the cocoon breaking out of the progenitor. Due to the pressure drop ahead of the jet head after it reaches the stellar surface, a rarefaction wave propagates back into the jet at the sound speed, which re-accelerates the jet to a relativistic velocity and therefore limits the gap period to within about ten seconds. This scenario therefore cannot explain gaps which are hundreds of seconds long. Instead, we ascribe such long time gaps to the behavior of the central engine, and suggest a fallback collapsar scenario for these bursts. In this scenario, the precursor is produced by a weak jet formed during the initial core collapse, possibly related to MHD processes associated with a short-lived proto-neutron star, while the main burst is produced by a stronger jet fed by fallback accretion onto the black hole resulting from the collapse of the neutron star. We have examined the propagation times of the weak precursor jet through the stellar progenitor. We find that the initial weak jet can break out of the progenitor in a time less than ten seconds (a typical precursor duration) provided that it has a moderately high relativistic Lorentz factor \Gamma>=10 (abridged).

GRB061007 is the brightest gamma-ray burst (GRB) to be detected by Swift and is accompanied by an exceptionally luminous afterglow that had a V-band magnitude <11.1 at 80s after the prompt emission. From the start of the Swift observations the afterglow decayed as a power law with a slope of \alpha_X=1.66+/-0.01 in the X-ray and \alpha_{opt}=1.64+/-0.01 in the UV/optical, up to the point that it was no longer detected above background in the optical or X-ray bands. The brightness of this GRB and the similarity in the decay rate of the X-ray, optical and gamma-ray emission from 100s after the trigger distinguish this burst from others and present a challenge to the fireball model. The lack of a cooling or jet break in the afterglow up to \~10^5s constrains any model that can produce the large luminosity observed in GRB061007, which we found to require either an excessively large kinetic energy or highly collimated outflow. Analysis of the multi-wavelength spectral and high-resolution temporal data taken with Swift suggest an early time jet-break to be a more plausible scenario. This must have occurred within 80s of the prompt emission, which places an upper limit on the jet opening angle of \theta_j=0.8deg. Such a highly collimated outflow resolves the energy budget problem presented in a spherical emission model, reducing the isotropic equivalent energy of this burst to E_{\gamma}^{corr}=10^{50} ergs; consistent with other GRBs.

We calculate the GeV spectra of GRB afterglows produced by inverse Compton scattering of the sub-MeV emission of these objects. We improve on earlier treatments by using refined afterglow parameters and new model developments motivated by recent Swift observations. We present time-dependent GeV spectra for standard, constant parameter models, as well as for models with energy injection and with time-varying parameters, for a range of burst parameters. We evaluate the limiting redshift to which such afterglows can be detected by the GLAST LAT, as well as AGILE.

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.

GRB 050904 at redshift z=6.29, discovered and observed by Swift and with spectroscopic redshift from the Subaru telescope, is the first gamma-ray burst to be identified from beyond the epoch of reionization. Since the progenitors of long gamma-ray bursts have been identified as massive stars, this event offers a unique opportunity to investigate star formation environments at this epoch. Apart from its record redshift, the burst is remarkable in two respects: first, it exhibits fast-evolving X-ray and optical flares that peak simultaneously at t~470 s in the observer frame, and may thus originate in the same emission region; and second, its afterglow exhibits an accelerated decay in the near-infrared (NIR) from t~10^4 s to t~3 10^4 s after the burst, coincident with repeated and energetic X-ray flaring activity. We make a complete analysis of available X-ray, NIR, and radio observations, utilizing afterglow models that incorporate a range of physical effects not previously considered for this or any other GRB afterglow, and quantifying our model uncertainties in detail via Markov Chain Monte Carlo analysis. In the process, we explore the possibility that the early optical and X-ray flare is due to synchrotron and inverse Compton emission from the reverse shock regions of the outflow. We suggest that the period of accelerated decay in the NIR may be due to suppression of synchrotron radiation by inverse Compton interaction of X-ray flare photons with electrons in the forward shock; a subsequent interval of slow decay would then be due to a progressive decline in this suppression. The range of acceptable models demonstrates that the kinetic energy and circumburst density of GRB 050904 are well above the typical values found for low-redshift GRBs.

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.

We present an interpretation of the phenomenological relations between the spectral peak, isotropic luminosity and duration of long gamma ray bursts that have been discovered by Amati et al., Ghirlanda et al., Firmani et al., and Liang & Zhang. In our proposed model, a jet undergoes internal dissipation which prevents its bulk Lorentz factor from exceeding 1/theta (theta being the jet opening angle) until it escapes from the core of its progenitor star, at a radius of order 1e10 cm; dissipation may continue at larger radii. The dissipated radiation will be partially thermalized, and we identify its thermal peak (Doppler boosted by the outflow) with E_pk. The radiation comes, in effect, from within the jet photosphere. The non-thermal, high energy part of the GRB emission arises from Comptonization of this radiation by relativistic electrons and positrons outside the effective photosphere. This model can account naturally not only for the surprisingly small scatter in the various claimed correlations, but also for the normalization, as well as the slopes. It then has further implications for the jet energy, the limiting jet Lorentz factor, and the relation of the energy, opening angle and burst duration to the mass and radius of the stellar stellar progenitor. The observed relation between pulse width and photon frequency can be reproduced by inverse-Compton emission at ~ 1e14 cm from the engine, but there are significant constraints on the energy distribution and isotropy of the radiating particles.

The low energy spectra of some gamma-ray bursts' show excess components beside the power-law dependence. The consequences of such a feature allows to estimate the gamma photometric redshift of the long gamma-ray bursts in the BATSE Catalog. There is good correlation between the measured optical and the estimated gamma photometric redshifts. The estimated redshift values for the long bright gamma-ray bursts are up to z=4, while for the the faint long bursts - which should be up to z=20 - the redshifts cannot be determined unambiguously with this method. The redshift distribution of all the gamma-ray bursts with known optical redshift agrees quite well with the BATSE based gamma photometric redshift distribution.

We report the results of the Swift and XMM observations of the Swift-discovered long Gamma-Ray Burst GRB 060729 ($T_{90}$=115s). The afterglow of this burst was exceptionally bright in X-rays as well as at UV/Optical wavelengths showing an unusually long slow decay phase ($\alpha$=0.14\plm0.02) suggesting a larger energy injection phase at early times than in other bursts. The X-ray light curve displays a break at about 60 ks after the burst. The X-ray decay slope after the break is $\alpha$=1.29\plm0.03. Up to 125 days after the burst we do not detect a jet break, suggesting that the jet opening angle is larger than 28 degrees. In the first 2 minutes after the burst (rest frame) the X-ray spectrum of the burst changed dramatically from a hard X-ray spectrum to a very soft one. We find that the X-ray spectra at this early phase can all be fitted by an absorbed single power law model or alternatively by a blackbody plus power law model. The power law fits show that the X-ray spectrum becomes steeper while the absorption column density decreases. In Swift's UV/Optical telescope the afterglow was clearly detected up to 9 days after the burst in all 6 filters and even longer in some of the UV filters with the latest detection in the UVW1 31 days after the burst. A break at about 50 ks is clearly detected in all 6 UVOT filters from a shallow decay slope of about 0.3 and a steeper decay slope of 1.3. In addition to the \swift observations we also present and discuss the results from a 61 ks ToO observation by XMM. (Abriviated)

The recent detection of X-ray flares during the afterglow phase of gamma-ray bursts (GRBs) suggests an inner-engine origin, at radii inside the forward shock. There must be inverse Compton (IC) emission arising from such flare photons scattered by forward shock afterglow electrons when they are passing through the forward shock. We find that this IC emission produces high energy gamma-ray flares, which may be detected by AGILE, GLAST and ground-based TeV telescopes. The anisotropic IC scattering between flare photons and forward shock electrons does not affect the total IC component intensity, but cause a time delay of the IC component peak relative to the flare peak. The anisotropic scattering effect may also weaken, to some extent, the suppression effect of the afterglow intensity induced by the enhanced electron cooling due to flare photons. We speculate that this IC component may already have been detected by EGRET from a very strong burst--GRB940217. Future observations by GLAST may help to distinguish whether X-ray flares originate from late central engine activity or from external shocks.

We provided separate bivariate log-normal distribution fits to the BATSE short and long burst samples using the durations and fluences. We show that these fits present an evidence for a power-law dependence between the fluence and the duration, with a statistically significant different index for the long and short groups. We argue that the effect is probably real, and the two subgroups are different physical phenomena. This may provide a potentially useful constraint for models of long and short bursts.

The absence of a supernova accompanying the nearby long GRB 060614 poses a great puzzle about the progenitor of this event and challenges the current GRB classification scheme. This burst displays a short-hard emission episode followed by extended soft emission with strong spectral evolution. Noticing that this burst has an isotropic gamma-ray energy only ~8 times that of GRB 050724, a good candidate of merger-type short GRBs, we generate a ``pseudo'' burst that is ~8 times less energetic than GRB 060614 based on the spectral properties of GRB 060614 and the Ep ~ Eiso^{1/2} (Amati) relation. We find that this pseudo-burst would have been detected by BATSE as a marginal short-duration GRB, and would have properties in the Swift BAT and XRT bands similar to GRB 050724. This suggests that GRB 060614 is likely a more intense event in the traditional short-hard GRB category as would be detected by BATSE. Events like GRB 060614 that seem to defy the traditional short vs. long classification of GRBs may require modification of our classification terminology for GRBs. By analogy with supernova classifications, we suggest that GRBs be classified into Type I (typically short and associated with old populations) and Type II (typically long and associated with young populations). We propose that GRB 060614 belongs to Type I, and predict that similar events will be detected in elliptical galaxies.

%auto-ignore This paper has been withdrawn by the authors due to dublicate submission. To download the paper please go to astro-ph/0611081

Swift discovered the high redshift (z=6.29) GRB050904 with the Burst Alert Telescope (BAT) and began observing with its narrow field instruments 161 s after the burst onset. This gamma-ray burst is the most distant cosmic explosion ever observed. Because of its high redshift, the X-ray Telescope (XRT) and BAT simultaneous observations provide 4 orders of magnitude of spectral coverage (0.2-150 keV; 1.4-1090 keV in the source rest frame) at a very early source-frame time (22 s). GRB050904 was a long, multi-peaked, bright GRB with strong variability during its entire evolution. The light curve observed by the XRT is characterized by the presence of a long flaring activity lasting up to 1-2 hours after the burst onset in the burst rest frame, with no evidence of a smooth power-law decay following the prompt emission as seen in other GRBs. However, the BAT tail extrapolated to the XRT band joins the XRT early light curve and the overall behavior resembles that of a very long GRB prompt. The spectral energy distribution softens with time, with the photon index decreasing from -1.2 during the BAT observation to -1.9 at the end of the XRT observation. The dips of the late X-ray flares may be consistent with an underlying X-ray emission arising from the forward shock and with the properties of the optical afterglow reported by Tagliaferri et al. (2005b). We interpret the BAT and XRT data as a single continuous observation of the prompt emission from a very long GRB. The peculiarities observed in GRB050904 could be due to its origin within one of the first star-forming regions in the Universe; very low metallicities of the progenitor at these epochs may provide an explanation.

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

We study synchrotron self-inverse Compton radiation from a reverse-shocked fireball. If the inverse Compton process dominates the cooling of shocked electrons, an X-ray flare produced by the first order Compton scattering would emerge in the very early afterglow phase, with the bulk of the shock energy radiated in the second order scattering component at 10-100 MeV. The dominance of inverse Compton cooling leads to the lack of prompt optical flashes. We show that for plausible parameters this scattering process can produce an X-ray flare with a relative amplitude change by a factor of several. Flares with a larger amplitude and multiple X-ray flares in a single event are likely to be produced by another mechanism (e.g. internal shocks).

We systematically analyze the prompt emission and the early afterglow data of a sample of 31 GRBs detected by {\em Swift} before September 2005, and estimate the GRB radiative efficiency. BAT's narrow band inhibits a precise determination of the GRB spectral parameters, and we have developed a method to estimate these parameters with the hardness ratio information. The shallow decay component commonly existing in early X-ray afterglows, if interpreted as continuous energy injection in the external shock, suggests that the GRB efficiency previously derived from the late-time X-ray data were not reliable. We calculate two radiative efficiencies using the afterglow kinetic energy E_K derived at the putative deceleration time t_{dec}) and at the break time (t_b) when the energy injection phase ends, respectively. At t_b XRFs appear to be less efficient than normal GRBs. However, when we analyze the data at t_{dec} XRFs are found to be as efficient as GRBs. Short GRBs have similar radiative efficiencies to long GRBs despite of their different progenitors. Twenty-two bursts in the sample are identified to have the afterglow cooling frequency below the X-ray band. Assuming \epsilon_e = 0.1, we find \eta_\gamma(t_b) usually <10% and \eta_\gamma (t_{dec}) varying from a few percents to > 90%. Nine GRBs in the sample have the afterglow cooling frequency above the X-ray band for a very long time. This suggests a very small \epsilon_B and/or a very low ambient density n.

We review the fireball shock model of gamma-ray burst prompt and early afterglow emission in light of rapid follow-up measurements made and enabled by the multi-wavelength Swift satellite. These observations are leading to a reappraisal and expansion of the previous standard view of the GRB and its fireball. New information on the behavior of the burst and afterglow on minutes to hour timescales has led, among other results, to the discovery and follow-up of short GRB afterglows, the opening up of the z>6 redshift range, and the first prompt multi-wavelength observations of a long GRB-supernova. We discuss the salient observational results and some associated theoretical issues.

We report the results of the \chandra observations of the \swift-discovered short Gamma-Ray Burst GRB 050724. \chandra observed this burst twice, about two days after the burst and a second time three weeks later. The first \chandra pointing occurred at the end of a strong late-time flare. About 150 photons were detected during this 49.3 ks observation in the 0.4-10.0 keV range. The spectral fit is in good agreement with spectral analysis of earlier \swift XRT data. In the second \chandra pointing the afterglow was clearly detected with 8 background-subtracted photons in 44.6 ks. From the combined \swift XRT and \chandra-ACIS-S light curve we find significant flaring superposed on an underlying power-law decay slope of $\alpha$=0.98$^{+0.11}_{-0.09}$. There is no evidence for a break between about 1 ks after the burst and the last \chandra pointing about three weeks after the burst. The non-detection of a jet break places a lower limit of 25$^{\circ}$ on the jet opening angle, indicating that the outflow is less strongly collimated than most previously-reported long GRBs. This implies that the beaming corrected energy of GRB 050724 is at least $4\times 10^{49}$ ergs.

We report the best evidence to date of a jet break in a short Gamma-Ray Burst (GRB) afterglow, using Chandra and Swift XRT observations of the X-ray afterglow of GRB 051221A. The combined X-ray light curve, which has three breaks, is similar to those commonly observed in Swift observations of long GRBs. A flat segment of the light curve at ~0.1 days after the burst represents the first clear case of strong energy injection in the external shock of a short GRB afterglow. The last break occurs at ~4 days post-burst and breaks to a power-law decay index of ~2. We interpret this as a jet break, with important implications for models of short GRBs, since it requires collimation of the afterglow into a jet with an initial opening angle ~4-8 degrees and implies a total jet kinetic energy of E_jet ~(1-5) x 10^{49} erg. Combined with the lack of a jet break in GRB 050724, this suggests a wide range in jet collimation in short GRBs, with at least some having collimation similar to that found in long GRBs, though with significantly lower jet energies.

We present a physical framework that can account for most of the observed spectral properties of the prompt gamma-ray burst emission. This includes the variety of spectral shapes, shape evolutions, and spectral correlations between flux and spectral peak, within bursts described by Borgonovo & Ryde, and among bursts described by Amati/Ghirlanda. In our proposed model the spectral peak is given by the photospheric emission from a relativistic outflow for which the horizon length is much smaller than the radial width. The observed duration of the thermal flash will be given by the radial light-crossing time. This then gives that the typical emission site is at ~10e11 cm from the center, with a Lorentz factor of ~300. This emission is accompanied by non-thermal emission from dissipation locations outside the photosphere. The relative strength of these two components depend on injection effects at the central engine leading to varying relative location of the saturation and photospheric radii. The total emission can then reproduce the observed variety. The spectral correlations are found by assuming that the amount of energy dissipated depends non-linearly on the averaged particle density. Beside the spectral correlations this also gives a description of how the relative strength of the thermal component varies with temperature within a burst.

Very early observations with the Swift satellite of gamma-ray burst (GRB) afterglows reveal that the optical component is not detected in a large number of cases. This is in contrast to the bright optical flashes previously discovered in some GRBs (e.g. GRB 990123 and GRB 021211). Comparisons of the X-ray afterglow flux to the optical afterglow flux and prompt gamma-ray fluence is used to quantify the seemingly deficient optical, and in some cases X-ray, light at these early epochs. This comparison reveals that some of these bursts appear to have higher than normal gamma-ray efficiencies. We discuss possible mechanisms and their feasibility for explaining the apparent lack of early optical emission. The mechanisms considered include: foreground extinction, circumburst absorption, Ly-alpha blanketing and absorption due to high redshift, low density environments, rapid temporal decay, and intrinsic weakness of the reverse shock. Of these, foreground extinction, circumburst absorption, and high redshift provide the best explanations for most of the non-detections in our sample. There is tentative evidence of suppression of the strong reverse shock emission. This could be because of a Poynting-flux-dominated flow or a pure non-relativistic hydrodynamical reverse shock.

We report the simultaneous observations of the prompt emission in the gamma-ray and hard X-ray bands by the Swift-BAT and the KONUS-Wind instruments of the short-hard burst, GRB 060313. The observations reveal multiple peaks in both the gamma-ray and hard X-ray bands suggesting a highly variable outflow from the central explosion. We also describe the early-time observations of the X-ray and UV/Optical afterglows by the Swift XRT and UVOT instruments. The combination of the X-ray and UV/Optical observations provide the most comprehensive lightcurves to date of a short-hard burst at such an early epoch. The afterglows exhibit complex structure with different decay indices and flaring. This behavior can be explained by the combination of a structured jet, radiative loss of energy, and decreasing microphysics parameters occurring in a circum-burst medium with densities varying by a factor of approximately two on a length scale of 10^17 cm. These density variations are normally associated with the environment of a massive star and inhomogeneities in its windy medium. However, the mean density of the observed medium (n approximately 10^&#8722;4 cm^3) is much less than that expected for a massive star. Although the collapse of a massive star as the origin of GRB 060313 is unlikely, the merger of a compact binary also poses problems for explaining the behavior of this burst. Two possible suggestions for explaining this scenario are: some short bursts may arise from a mechanism that does not invoke the conventional compact binary model, or soft late-time central engine activity is producing UV/optical but no X-ray flaring.

Recent discovery of the afterglow emission from short gamma-ray bursts suggests that binary neutron star or black hole-neutron star binary mergers are the likely progenitors of these short bursts. The accretion of neutron star material and its subsequent ejection by the central engine implies a neutron-rich outflow. We consider here a neutron-rich relativistic jet model of short bursts, and investigate the high energy neutrino and photon emission as neutrons and protons decouple from each other. We find that upcoming neutrino telescopes are unlikley to detect the 50 GeV neutrinos expected in this model. For bursts at z~0.1, we find that GLAST and ground-based Cherenkov telescopes should be able to detect prompt 100 MeV and 100 GeV photon signatures, respectively, which may help test the neutron star merger progenitor identification.

Core collapse of massive stars and binary neutron stars or black hole-neutron star binary mergers are likely progenitors of long and short duration gamma-ray bursts respectively. Neutronized material in the former and neutron star material in the latter are ejected by the central engine implying a neutron-rich jet outflow. A free neutron, however, beta decays to a proton, an electron (beta) and an anti-neutrino in about fifteen minutes in its rest frame. Sudden creation of a relativistic electron is accompanied by radiation with unique temporal and spectral signature. We calculate here this radiation signature collectively emitted by all beta decay electrons from neutron-rich outflow. Detection of this signature may thus provide strong evidence for not only neutron but also for proton content in the relativistic gamma-ray burst jets.

Gamma-ray bursts are the most luminous explosions in the Universe, and their origin and mechanism are the focus of intense research and debate. More than three decades after their discovery, and after pioneering breakthroughs from space and ground experiments, their study is entering a new phase with the recently launched Swift satellite. The interplay between these observations and theoretical models of the prompt gamma ray burst and its afterglow is reviewed.

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.

Neutrinos at energies ranging from sub-TeV to EeV from astrophysical sources can yield interesting physical information about fundamental interactions, about cosmic rays and about the nature of the sources and their environment. Gamma-ray bursts are a leading candidate source, and their expected neutrino emission can address a number of current questions, which may be answered with forthcoming experiments such as IceCube, Auger, ANITA and KM3NeT.

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.

Shock breakout X-ray emission has been reported for the first time from a supernova connected with a gamma-ray burst, namely GRB060218/SN2006aj. The gamma-ray emission and the power-law decaying X-ray afterglow are ascribed to a highly relativistic jet, while the thermal soft X-rays are thought to be produced when the radiation-dominated shock breaks from the optically thick stellar wind. We study the inverse Compton emission of the breakout thermal soft X-rays scattered by relativistic electrons in the jet forward shock, which is expected to be at larger radii than the breakout shock. This IC emission produces sub-GeV to GeV photons, which may be detectable by GLAST. The detection of such GeV photons simultaneously with the supernova shock breakout emission would give evidence for the presence of a GRB jet ahead of the shock while the shock is breaking out. The anisotropic scattering between the X-rays and relativistic electrons may lead to large angle emission outside of the jet opening angle. This has implications for the detection of GeV photons from burstless type Ib/c hypernova shock breakout, which due to its more isotropic emission might be observed with wide-field X-ray cameras such as LOBSTER.

The power-law decay of the X-ray emission of GRB afterglows 050319, 050401, 050607, 050713A, 050802 and 050922C exhibits a steepening at about 1--4 hours after the burst which, surprisingly, is not accompanied by a break in the optical emission. If it is assumed that both the optical and X-ray afterglows arise from the same outflow then, in the framework of the standard forward shock model, the chromaticity of the X-ray light-curve breaks indicates that they do not arise solely from a mechanism related to the outflow dynamics (e.g. energy injection) or the angular distribution of the blast-wave kinetic energy (structured outflows or jets). The lack of a spectral evolution accompanying the X-ray light-curve breaks shows that these breaks do not arise from the passage of a spectral break (e.g. the cooling frequency) either. Under these circumstances, the decoupling of the X-ray and optical decays requires that the microphysical parameters for the electron and magnetic energies in the forward shock evolve in time, whether the X-ray afterglow is synchrotron or inverse-Compton emission. For a steady evolution of these parameters with the Lorentz factor of the forward shock and an X-ray light-curve break arising from cessation of energy injection into the blast-wave, the optical and X-ray properties of the above six Swift afterglows require a circumburst medium with a r^{-2} radial stratification, as expected for a massive star origin for long GRBs. Alternatively, the chromatic X-ray light-curve breaks may indicate that the optical and X-ray emissions arise from different outflows. Neither feature (evolution of microphysical parameters or the different origin of the optical and X-ray emissions) were clearly required by pre-Swift afterglows.

Gamma-ray Bursts (GRBs) are bright flashes of high energy photons that can last from about 10 milliseconds to 10 minutes. Their origin and nature have puzzled the scientific community for about 25 years until 1997, when the first X-ray afterglows of long (> 2 s duration) bursts were detected and the first optical and radio counterparts were found. These measurements established that long GRBs are typically at high redshift (z 1.6) and are in sub-luminous star-forming host galaxies. They are likely produced in core-collapse explosions of a class of massive stars that give rise to highly relativistic jets (collapsar model). Internal inhomogeneities in the velocity field of the relativistic expanding flow lead to collisions between fast moving and slow moving fluid shells and to the formation of internal shock waves. These shocks are believed to produce the observed prompt emission in the form of irregularly shaped and spaced pulses of gamma-rays, each pulse corresponding to a distinct internal collision. The expansion of the jet outward into the circumstellar medium is believed to give rise to ``external'' shocks, responsible for producing the smoothly fading afterglow emission seen in the X-ray, optical and radio bands. Here we report on the gamma-ray and x-ray observation of the most distant gamma-ray burst ever observed: its redshift of 6.29 translates to a distance of 13 billion light-years from Earth, corresponding to a time when the Universe was just 700 million to 750 million years old. The discovery of a gamma-ray burst at such a large redshift implies the presence of massive stars only 700 million years after the Big Bang. The very high redshift bursts represent a good way to study the re-ionization era soon after the Universe came out of the Dark Ages.

We report the results of Swift observations of the Gamma Ray Burst GRB 050603. With a V magnitude V=18.2 about 10 hours after the burst the optical afterglow was the brightest so far detected by Swift and one of the brightest optical afterglows ever seen. The Burst Alert Telescope (BAT) light curves show three fast-rise-exponential-decay spikes with $T_{90}$=12s and a fluence of 7.6$\times 10^{-6}$ ergs cm$^{-2}$ in the 15-150 keV band. With an $E_{\rm \gamma, iso} = 1.26 \times 10^{54}$ ergs it was also one of the most energetic bursts of all times. The Swift spacecraft began observing of the afterglow with the narrow-field instruments about 10 hours after the detection of the burst. The burst was bright enough to be detected by the Swift UV/Optical telescope (UVOT) for almost 3 days and by the X-ray Telescope (XRT) for a week after the burst. The X-ray light curve shows a rapidly fading afterglow with a decay index $\alpha$=1.76$^{+0.15}_{-0.07}$. The X-ray energy spectral index was $\beta_{\rm X}$=0.71\plm0.10 with the column density in agreement with the Galactic value. The spectral analysis does not show an obvious change in the X-ray spectral slope over time. The optical UVOT light curve decays with a slope of $\alpha$=1.8\plm0.2. The steepness and the similarity of the optical and X-ray decay rates suggest that the afterglow was observed after the jet break. We estimate a jet opening angle of about 1-2$^{\circ}$

The recent detection of delayed X-ray flares during the afterglow phase of gamma-ray bursts (GRBs) suggests an inner-engine origin, at radii inside the deceleration radius characterizing the beginning of the forward shock afterglow emission. Given the observed temporal overlapping between the flares and afterglows, there must be inverse Compton (IC) emission arising from such flare photons scattered by forward shock afterglow electrons. We find that this IC emission produces GeV-TeV flares, which may be detected by GLAST and ground-based TeV telescopes. We speculate that this kind of emission may already have been detected by EGRET from a very strong burst--GRB940217. The enhanced cooling of the forward shock electrons by the X-ray flare photons may suppress the synchrotron emission of the afterglows during the flare period. The detection of GeV-TeV flares combined with low energy observations may help to constrain the poorly known magnetic field in afterglow shocks. We also consider the self-IC emission in the context of internal-shock and external-shock models for X-ray flares. The emission above GeV from internal shocks is low, while the external shock model can also produce GeV-TeV flares, but with a different temporal behavior from that caused by IC scattering of flare photons by afterglow electrons. This suggests a useful approach for distinguishing whether X-ray flares originate from late central engine activity or from external shocks.

The short/hard GRB051210 was detected and located by the Swift-BAT instrument and rapidly pointed towards by the narrow field instrumens. The XRT was able to observe a bright X-ray afterglow, one of the few ever observed for this class of bursts. We present the analysis of the prompt and afterglow emission of this event The BAT spectrum is a power-law with photon index 1.1 +/-0.3. The X-ray light curve decays with slope 2.58+/-0.11 and shows a small flare in the early phases. The spectrum can be described with a power law with photon index 1.54+/-0.16 and absorption (7.5 (-3.2, +4.3)*10^20 cm-2 We find that the X-ray emission is consistent with the hypothesis that we are observing the curvature effect of a GRB occurred in a low density medium, with no detectable afterglow. We estimate the density of the circumburst medium to be lower than 4*10^-3 cm^-3. We also discuss different hypothesis on the possible origin of the flare.

Observations of gamma ray bursts (GRBs) with Swift produced the initially surprising result that many bursts have large X-ray flares superimposed on the underlying afterglow. These flares were sometimes intense, rapid, and late relative to the nominal prompt phase. The most intense of these flares was observed by XRT with a flux >500 the afterglow. This burst then surprised observers by flaring again after >10000 s. The intense flare can be most easily understood within the context of the standard fireball model, if the internal engine that powers the prompt GRB emission is still active at late times. Recent observations indicate that X-ray flares are detected in ~1/3 of XRT detected afterglows. By studying the properties of the varieties of flares (such as rise/fall time, onset time, spectral variability, etc.) and relating them to overall burst properties, models of flare production and the GRB internal engine can be constrained.

The rapid follow-up of gamma-ray burst (GRB) afterglows made possible by the multi-wavelength satellite Swift, launched in November 2004, has put under a microscope the GRB early post-burst behavior, This is leading to a significant reappraisal and expansion of the standard view of the GRB early afterglow behavior, and its connection to the prompt gamma-ray emission. In addition to opening up the previously poorly known behavior on minutes to hours timescales, two other new pieces in the GRB puzzle being filled in are the the discovery and follow-up of short GRB afterglows, and the opening up of the $z\simg 6$ redshift range. We review some of the current theoretical interpretations of these new phenomena.

We present observations of XRF 050406, the first burst detected by Swift showing a flare in its X-ray light curve. During this flare, which peaks at t_peak ~210s after the BAT trigger, a flux variation of (delta F)/F~6 in a very short time (delta t)/t_peak<<1 was observed. Its measured fluence in the 0.2-10 keV band was ~1.4x10^-8 erg cm^-2, which corresponds to 1-15% of the prompt fluence. We present indications of spectral variations during the flare. We argue that the producing mechanism is late internal shocks, which implies that the central engine is still active at 210s, though with a reduced power with respect to the prompt emission. The X-ray light curve flattens to a very shallow slope with decay index of ~0.5 after ~4400s, which also supports continued central engine activity at late times. This burst is classified as an X-ray flash, with a relatively low fluence (~10^-7 erg cm^-2 in the 15-350 keV band, E_iso~10^51 erg), a soft spectrum (photon index 2.65), no significant flux above ~50 keV and a peak energy E_p<15 keV. XRF 050406 is one of the first examples of a well-studied X-ray light curve of an XRF. We show that the main afterglow characteristics are qualitatively similar to those of normal GRBs. In particular, X-ray flares superimposed on a power-law light curve have now been seen in both XRFs and GRBs. This indicates that a similar mechanism may be at work for both kinds of events.

Until recently, X-ray flares during the afterglow of gamma ray bursts (GRBs) were a rarely detected phenomenon, thus their nature is unclear. During the afterglow of GRB 050502B, the largest X-ray flare ever recorded rose rapidly above the afterglow lightcurve detected by the Swift X-ray Telescope. The peak flux of the flare was >500 times that of the underlying afterglow, and it occurred at >12 minutes after the nominal prompt burst emission. The fluence of this X-ray flare, (1.0 +/- 0.05) x 10^{-6} erg cm^{-2} in the 0.2-10.0 keV energy band, exceeded the fluence of the nominal prompt burst. The spectra during the flare were significantly harder than those measured before and after the flare. Later in time, there were additional flux increases detected above the underlying afterglow, as well as a break in the afterglow lightcurve. All evidence presented below, including spectral and particularly timing information during and around the giant flare, suggests that this giant flare was the result of internal dissipation of energy due to late central engine activity, rather than an afterglow-related effect. We also find that the data are consistent with a second central engine activity episode, in which the ejecta is moving slower than that of the initial episode, causing the giant flare and then proceeding to overtake and refresh the afterglow shock, thus causing additional activity at even later times in the lightcurve.

(Abridged) The Swift X-Ray Telescope (XRT) reveals some interesting features of early X-ray afterglows, including a distinct rapidly decaying component preceding the conventional afterglow component in many sources, a shallow decay component before the more ``normal'' decay component observed in a good fraction of GRBs (e.g. GRB 050128, GRB 050315, GRB 050319, and GRB 050401), and X-ray flares in nearly half of the afterglows (e.g. GRB 050406, GRB 050502B, GRB 050607, and GRB 050724). In this paper, we systematically analyze the possible physical processes that shape the properties of the early X-ray afterglow lightcurves, and use the data to constrain various models. We suggest that the steep decay component is consistent with the tail emission of the prompt gamma-ray bursts and/or of the X-ray flares. This provides clear evidence that the prompt emission and afterglow emission are two distinct components, supporting the internal origin of the GRB prompt emission. The shallow decay segment observed in a group of GRBs suggests that the forward shock keeps being refreshed for some time. This might be caused either by a long-lived central engine, or by a power law distribution of the shell Lorentz factors, or else by the deceleration of a Poynting flux dominated flow. X-ray flares suggest that the GRB central engine is still active after the prompt gamma-ray emission is over, but with a reduced activity at later times. In some cases, the central engine activity even extends days after the burst trigger. Analyses of early X-ray afterglow data reveal that GRBs are indeed highly relativistic events. Early afterglow data of many bursts, starting from the beginning of the XRT observations, are consistent with the afterglow emission from an interstellar medium (ISM) environment.

We present Swift and XMM observations of GRB 050326, detected by Swift-BAT. The fluence was 7.7x10^-6 erg cm^-2 (20-150 keV), and its spectrum was hard, with a power law photon index 1.25. The afterglow light curve did not show any break nor flares between ~1 hr and ~6 d after the burst, and decayed with a slope 1.70. The afterglow spectrum is well fitted by a power-law model, suffering absorption both in the Milky Way and in the host galaxy. The rest-frame Hydrogen column density is significant, N_H_z > 4x10^21 cm^-2, and the redshift of the absorber is z > 1.5. There was good agreement between the Swift-XRT and XMM results. By comparing the prompt and afterglow fluxes, we found that an early break occurred before the XRT observation. The properties of the GRB 050326 afterglow are well described by a spherical fireball expanding in a uniform external medium, so a further steepening is expected at later times. The lack of such a break constrains the jet angle to be >7 deg. Using the redshift constraints provided by the X-ray analysis, we also estimated that the beaming-corrected gamma-ray energy was >3x10^51 erg, at the high end of GRB energies. Despite the brightness in X rays, only deep limits could be placed by Swift-UVOT at optical/UV wavelengths. Thus, this GRB was "truly dark", with the optical-to-X-ray spectrum violating the synchrotron limit. The optical and X-ray observations are consistent either with an absorbed event or with a high-redshift one. To obey the Ghirlanda relation, a moderate/large redshift z>4.5 is required. (abridged)

The X-ray light-curves of 9 Swift XRT afterglows (050126, 050128, 050219A, 050315, 050318, 050319, 050401, 050408, 050505) display a complex behaviour: a steep t^{-3.0 \pm 0.3} decay until ~400 s, followed by a significantly slower t^{-0.65+/-0.20} fall-off, which at 0.2--2 d after the burst evolves into a t^{-1.7+/-0.5} decay. We consider three possible models for the geometry of relativistic blast-waves (spherical outflows, non-spreading jets, and spreading jets), two possible dynamical regimes for the forward shock (adiabatic and fully radiative), and we take into account a possible angular structure of the outflow and delayed energy injection in the blast-wave, to identify the models which reconcile the X-ray light-curve decay with the slope of the X-ray continuum for each of the above three afterglow phases. By piecing together the various models for each phase in a way that makes physical sense, we identify possible models for the entire X-ray afterglow. The major conclusion of this work is that a long-lived episode of energy injection in the blast-wave, during which the shock energy increases at t^{1.0+/-0.5}, is required for five afterglows and could be at work in the other four as well. Optical observations in conjunction with the X-ray can distinguish among these various models. Our simple tests allow the determination of the location of the cooling frequency relative to the X-ray domain and, thus, of the index of the electron power-law distribution with energy in the blast-wave. The resulting indices are clearly inconsistent with an universal value.

BAT and XRT observations of two recent well-covered GRBs observed by Swift, GRB 050315 and GRB 050319, show clearly a prompt component joining the onset of the afterglow emission. By fitting a power law form to the gamma-ray spectrum, we extrapolate the time dependent fluxes measured by the BAT, in the energy band 15-350 keV, into the spectral regime observed by the XRT 0.2-10 keV, and examine the functional form of the rate of decay of the two light curves. We find that the BAT and XRT light curves merge to form a unified curve. There is a period of steep decay up to ~300 s, followed by a flatter decay. The duration of the steep decay, ~100 s in the source frame after correcting for cosmological time dilation, agrees roughly with a theoretical estimate for the deceleration time of the relativistic ejecta as it interacts with circumstellar material. For GRB 050315, the steep decay can be characterized by an exponential form, where one e-folding decay time t_e(BAT) ~ 24+/-2 s, and t_e(XRT) ~ 35+/-2 s. For GRB 050319, a power law decay -(d ln f/d ln t) = n, where n is about 3, provides a reasonable fit. The early time X-ray fluxes are consistent with representing the lower energy tail of the prompt emission, and provide our first quantitative measure of the decay of the prompt gamma-ray emission over a large dynamic range in flux. The initial steep decay is expected due to the delayed high latitude photons from a curved shell of relativistic plasma illuminated only for a short interval. The overall conclusion is that the prompt phase of GRBs remains observable for hundreds of seconds longer than previously thought.

Two short (<2 s) g-ray bursts (GRBs) have recently been localized and fading afterglow counterparts detected. The combination of these two results left unclear the nature of the host galaxies of the bursts, because one was a star-forming dwarf, while the other was probably an elliptical galaxy. Here we report the X-ray localization of a short burst (GRB 050724) with unusual g-ray and X-ray properties. The X-ray afterglow lies off the centre of an elliptical galaxy at a redshift of z=0.258, coincident with the position determined by ground-based optical and radio observations. The low level of star formation typical for elliptical galaxies makes it unlikely that the burst originated n a supernova explosion. A supernova origin was also ruled out for GRB 050709, even though that burst took place in a galaxy with current star formation. The isotropic energy for the short bursts is 2-3 orders of magnitude lower han that for the long bursts. Our results therefore suggest that an alternative source of bursts -- the of binary systems of neutron stars or a neutron star-black hole pair -- are the progenitors of short bursts.

The Swift XRT has been observing GRB afterglows since December 23, 2004. Three-quarters of these observations begin within 300 s of the burst onset, providing an unprecendented look at the behavior of X-ray emission from GRB afterglows in the first few hours after the burst. While most of the early afterglows have smoothly declining lightcurves, a substantial fraction has large X-ray flares on short time-scales. We suggest that these flares provide support for models with extended central engine activity producing late-time internal shocks.

We consider the spatial clustering of massive black hole (MBH) mergers, and discuss possible ways to use gravitational wave observations in the LISA and DECIGO/BBO range for obtaining cosmological and cosmogonical information. Constraints on large scale structure (LSS) and merger histories may be possible by detections of an alignment of the GW polarization direction with principal axes of the LSS. Constraints on the merger physics and the reionization epoch may be obtained by GW measurements of MBH correlation lengths, in the case when the MBH angular momentum loss occurs through gas drag. Such measurements would provide information about the LSS and the reionization epoch, additional to and independent of any obtained from electromagnetic signals, as well as about the astrophysics of MBH mergers.

This paper discusses Swift observations of the gamma-ray burst GRB 050315 (z=1.949) from 80 s to 10 days after the onset of the burst. The X-ray light curve displayed a steep early decay (t^-5) for ~200 s and several breaks. However, both the prompt hard X-ray/gamma-ray emission (observed by the BAT) and the first ~ 300 s of X-ray emission (observed by the XRT) can be explained by exponential decays, with similar decay constants. Extrapolating the BAT light curve into the XRT band suggests the rapidly decaying, early X-ray emission was simply a continuation of the fading prompt emission; this strong similarity between the prompt gamma-ray and early X-ray emission may be related to the simple temporal and spectral character of this X-ray rich GRB. The prompt (BAT) spectrum was a steep down to 15 keV, and appeared to continue through the XRT bandpass, implying a low peak energy, inconsistent with the Amati relation. Following the initial steep decline the X-ray afterglow did not fade for ~1.2*10^4 s, after which time it decayed with a temporal index of alpha ~ 0.7, followed by a second break at ~2.5*10^5 s to a slope of alpha ~ 2. The apparent `plateau' in the X-ray light curve, after the early rapid decay, makes this one of the most extreme examples of the steep-flat-steep X-ray light curves revealed by Swift. If the second afterglow break is identified with a jet break then the jet opening angle was theta_0 ~ 5 deg, and implying E_gamma > 10^50 erg.

The Swift Gamma-Ray Burst Explorer performed its first autonomous, X-ray follow-up to a newly detected GRB on 2005 January 17, within 193 seconds of the burst trigger by the Swift Burst Alert Telescope. While the burst was still in progress, the X-ray Telescope obtained a position and an image for an un-catalogued X-ray source; simultaneous with the gamma-ray observation. The XRT observed flux during the prompt emission was 1.1 x 10^{-8} ergs cm^{-2} s^{-1} in the 0.5-10 keV energy band. The emission in the X-ray band decreased by three orders of magnitude within 700 seconds, following the prompt emission. This is found to be consistent with the gamma-ray decay when extrapolated into the XRT energy band. During the following 6.3 hours, the XRT observed the afterglow in an automated sequence for an additional 947 seconds, until the burst became fully obscured by the Earth limb. A faint, extremely slowly decaying afterglow, alpha=-0.21$, was detected. Finally, a break in the lightcurve occurred and the flux decayed with alpha<-1.2$. The X-ray position triggered many follow-up observations: no optical afterglow could be confirmed, although a candidate was identified 3 arcsecs from the XRT position.

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.

Swift discovered the high redshift GRB 050319 with the Burst Alert Telescope and began observing with its narrow field instruments only 225 s after the burst onset. The afterglow X-ray emission was monitored by the XRT up to 28 days after the burst. The light curve shows a decay with three different phases, each characterized by a distinct slope: an initial steep decay with a power law index of ~ 5.5, a second phase characterized by a flat decay slope of \~ 0.54, and a third phase with a decay slope of ~ 1.14. During the first phase the spectral energy distribution is softer than in the following two phases and the photon index is consistent with the GRB prompt spectrum. The extrapolation of the BAT light curve to the XRT band suggests that the initial fast decaying phase of the XRT afterglow might be the low energy tail of the prompt emission. The second break in the afterglow light curve occurs about 27000 s after the burst. The spectral energy distribution before and after the second break does not change and it can be tentatively interpreted as a jet break or the end of a delayed or continuous energy injection phase.

We analyze the Swift XRT light curves and spectra of the gamma-ray bursts (GRBs) for which the redshift has been measured. The sample consists of seven GRBs. The soft X-ray light curves of all these GRBs are separated into at least two morphological classes: 1) those starting off with a very steep light curve decay and 2) those showing a rather mild initial decay. This initial decay is followed by a flattening and by a further steepening. During these transitions the soft X-ray spectrum of these GRBs remains constant within the observational errors (except for one case, GRB050319). For the first time we are able to exploit the early light curve of GRB afterglows in the comoving frame. Besides the temporal and spectral behavior we find that the energy of the afterglow emitted in the (rest frame) time interval 20-200 s and 1300-12600 s after the trigger correlates with the mean energy of the prompt emission, hinting at a close link between the two. Averaging over the bursts that have a rather well sampled light curve and starting immediately after the prompt phase, the energy emitted in the soft X-rays is about 3% of the energy emitted during the prompt phase, except for GRB050315, where the soft X-ray emission is as high as 14% of the prompt emission.

We present an analysis indicating that there is a correlation between the fluences and the durations of gamma-ray bursts, and provide arguments that this reflects a correlation between the total emitted energies and the intrinsic durations. For the short (long) bursts the total emitted energies are roughly proportional to the first (second) power of the intrinsic duration. This difference in the energy-duration relationship is statistically significant, and may provide an interesting constraint on models aiming to explain the short and long gamma-ray bursts.

Gamma-ray burst (GRB) afterglows have provided important clues to the nature of these massive explosive events, providing direct information on the nearby environment and indirect information on the central engine that powers the burst. We report the discovery of two bright X-ray flares in GRB afterglows, including a giant flare comparable in total energy to the burst itself, each peaking minutes after the burst. These strong, rapid X-ray flares imply that the central engines of the bursts have long periods of activity, with strong internal shocks continuing for hundreds of seconds after the gamma-ray emission has ended.

We estimate the high energy neutrino flux from the giant flare of SGR 1806-20 on December 27, 2004, which irradiated Earth with a gamma-ray flux ~10^4 times larger than the most luminous gamma-ray bursts (GRBs) ever detected. The Antarctic Cherenkov neutrino detector AMANDA was on-line during the flare, and may either have detected high energy neutrinos for the first time from a cosmic point source, or put constraints on the flare mechanism of magnetars. If TeV neutrinos are detected, one would expect also detectable EeV cosmic rays and possibly TeV gamma-ray emission in coincidence.

[Abridged] Very recently, relations between the peak energy of Gamma-Ray Burst burst spectra, the isotropic-equivalent energy of the burst, and the radiated energy of the burst have been found. In a way that is exactly analogous to the way in which the relation between the peak luminosity and the rate of decline of the light curve of Type Ia supernovae can be used to make Type Ia supernovae excellent standard candles for cosmology, so too, the relations between Epk, Eiso, and Egamma point toward a methodology for using GRBs as excellent standard candles for cosmology. In addition, GRBs occur over the broad redshift range from z=0.1 to at least z=4.5, and both they and their afterglows are easily detectable out to z > 8. Thus GRBs show great promise as cosmological ``yardsticks'' to measure the rate of expansion of the universe over time, and therefore the properties of dark energy (i.e., Omega_M, Omega_Lambda, w_0, and w_a). We describe a concept for a possible MIDEX-class mission dedicated to using GRBs to constrain the properties of dark energy that would obtain these quantities for > 800 bursts in the redshift range 0.1 >~ z <~ 10 during a 2-year mission. This burst sample would enable both Omega_M and w_0 to be determined to +/- 0.07 and +/- 0.06 (68% CL), respectively, and w_a to be significantly constrained.

Long gamma-ray bursts (GRBs) are commonly accepted to originate in the explosion of particularly massive stars, which gives rise to a highly relativistic jet. Internal inhomogeneities in the expanding flow give rise to internal shock waves that are believed to produce the gamma-rays we see. As the jet travels further outward into the surrounding circumstellar medium `external' shocks give rise to the afterglow emission seen in the X-ray, optical and radio bands. Here we report on the early phases of the X-ray emission of five GRBs. Their X-ray light curves are characterised by a rapid fall-off for the first few hundred seconds, followed by a less rapid decline lasting several hours. This steep decline, together with detailed spectral properties of two particular bursts, shows that violent shock interactions take place in the early jet outflows.

Swift discovered GRB050128 with the Burst Alert Telescope and promptly pointed its narrow field instruments to monitor the afterglow. X-ray observations started 108 s after the trigger time. The early decay of the afterglow is relatively flat with a temporal decay modeled with a power law with index ~ -0.3. A steepening occurs at later times (~ 1500 s) with a power law index of ~ -1.3. During this transition, the observed X-ray spectrum does not change. We interpret this behaviour as either an early jet break or evidence for a transition from the fast cooling regime to the slow cooling regime in a wind environment.

We consider dissipative effects occurring in the optically thick inner parts of the relativistic outflows producing gamma-ray bursts and X-ray flashes, emphasizing specially the Comptonization of the thermal radiation flux that is advected from the base of the outflow. Such dissipative effects --e.g. from magnetic reconnection, neutron decay or shocks -- would boost the energy density of the thermal radiation. The dissipation can lead to pair production, in which case the pairs create an effective photosphere further out than the usual baryonic one. In a slow dissipation scenario, pair creation can be suppressed, and the effects are most important when dissipation occurs below the baryonic photosphere. In both cases an increased photospheric luminosity is obtained. We suggest that the spectral peak in gamma ray bursts is essentially due to the Comptonized thermal component from the photosphere, where the comoving optical depth in the outflow falls to unity. Typical peak photon energies range between those of classical bursts and X-ray flashes. The relationship between the observed photon peak energy and the luminosity depends on the details of the dissipation, but under plausible assumptions can resemble the observed correlations.

Short-duration Gamma-ray bursts (GRBs) (\leq 2{\rm s}) have remained a mystery due to the lack of afterglow detection until recently. The models to interpret short GRBs invoke distinct progenitor scenarios. Here we present a generic analysis of short GRB afterglows, and calculate the optical lightcurves of short GRBs within the framework of different progenitor models. We show that all these optical afterglows are bright enough to be detected by the Ultraviolet and Optical Telescope (UVOT) on board the {\em Swift} observatory, and that different models could be distinguished with a well-monitored lightcurve. We also model the afterglow data of the recently discovered short burst GRB 040924. We find that the limited data are consistent with a low medium density environment which is consistent with the pre-concept of the compact-star merger progenitor model, although the models with a collapsar progenitor are not ruled out.

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 Swift Gamma-Ray Burst Explorer, launched on 2004 November 20, is a multiwavelength, autonomous, rapid-slewing observatory for gamma-ray burst (GRB) astronomy. On 2004 December 23, during the activation phase of the mission, the Swift X-Ray Telescope (XRT) was pointed at a burst discovered earlier that day by the Swift Burst Alert Telescope. A fading, uncataloged X-ray source was discovered by the XRT and was observed over a period of about 3 hours, beginning 4.6 hours after the burst. The X-ray detection triggered a VLT observation of the optical/NIR counterpart, located about 1.1 arcseconds from the XRT position. The X-ray counterpart faded rapidly, with a power law index of -1.72 +/- 0.20. The average unabsorbed X-ray flux 4.6-7.9 hours after the burst was 6.5 x 10^{-12} erg cm^{-2} s^{-1} in the 0.5-10 keV band, for a power-law spectrum of photon index 2.02 +/- 0.13 with Galactic absorption. The NIR counterpart was observed at three epochs between 16 and 87 hours after the burst, and faded with a power-law index of -1.14 +/- 0.08 with a reddening-corrected SED power-law slope of -0.40 +/- 0.03. We find that the X-ray and NIR data are consistent with a two-component jet in a wind medium, with an early jet break in the narrow component and an underlying electron index of 1.8-2.0.

The rapid acquisition of positions by the upcoming Swift satellite will allow the monitoring for X-ray lines in GRB afterglows at much earlier epochs than was previously feasible. We calculate the possible significance levels of iron line detections as a function of source redshift and observing time after the trigger, for the Swift XRT, Chandra ACIS and XMM Epic detectors. For bursts with standard luminosities, decay rates and equivalent widths of 1 keV assumed constant starting at early source-frame epochs, Swift may be able to detect lines up to z~1.5 with a significance of better than 3 sigma for times up to 10^4 s. The same lines would be detectable with better than 4 sigma significance at z up to 6 by Chandra, and up to 8 by XMM, for times of up to 10^5 s. For similar bursts with a variable equivalent width peaking at 1 keV between 0.5 and 1 days in the source frame, Swift achieves the same significance level for z~1 at t~1 day, while Chandra reaches the previous detection significances around t~ 1-2 days for z~ 2-4, i.e. the line is detectable near the peak equivalent width times, and undetectable at earlier or later times. For afterglows in the upper range of initial X-ray luminosites afterglows, which may also be typical of pop. III bursts, similar significance levels are obtained out to substantially higher redshifts. A distinction between broad and narrow lines to better than 3 sigma is possible with Chandra and XMM out to z~2 and ~6.5, respectively, while Swift can do so up to z~1, for standard burst parameters. A distinction between different energy centroid lines of 6.4 keV vs. 6.7 KeV (or 6.7 keV vs. Cobalt 7.2 keV) is possible up to z~0.6, 1.2, and 2 (z~ 1, 5, 7.5), with Swift,Chandra, and XMM respectively.

In 2002, Ray Davis and Masatoshi Koshiba were awarded the Nobel Prize in Physics ``for pioneering contributions to astrophysics, in particular for the detection of cosmic neutrinos.'' However, while astronomy has undergone a revolution in understanding by synthesizing data taken at many wavelengths, the universe has only barely been glimpsed in neutrinos, just the Sun and the nearby SN 1987A. An entire universe awaits, and since neutrinos can probe astrophysical objects at densities, energies, and distances that are otherwise inaccessible, the results are expected to be particularly exciting. Similarly, the revolution in quantitative cosmology has heightened the need for very precise tests that depend on the effects of neutrinos, and prominent among them is the search for the effects of neutrino mass, since neutrinos are a small but known component of the dark matter. In this report, we highlight some of the key opportunties for progress in neutrino astrophysics and cosmology, and the implications for other areas of physics.

We investigate the radio afterglows of gamma-ray bursts (GRBs) and hypernovae (HNe) at high redshifts and quantify their detectability, as well as their potential usefulness for 21 cm absorption line studies of the intergalactic medium (IGM) and intervening structures. We examine several sets of source and environment model parameters that are physically plausible at high redshifts. The radio afterglows of GRBs would be detectable out to z ~ 30, while the energetic HNe could be detectable out to z ~ 20 even by the current Very Large Array (VLA). We find that the 21 cm absorption line due to the diffuse neutral IGM is difficult to detect even by the proposed Square Kilometer Array (SKA), except for highly energetic sources. We also find that the 21 cm line due to collapsed gas clouds with high optical depth may be detected on rare occasions.

Most active galactic nuclei (AGN) lack prominent jets, and show modest radio emission and significant X-ray emission which arises mainly from the galactic core, very near from the central black hole. We use a quantitative scenario of such core-dominated radio-quiet AGN, which attributes a substantial fraction of the X-ray emission to the presence of abortive jets involving the collision of gas blobs in the core. Here we investigate the consequences of the acceleration of protons in the shocks from such collisions. We find that protons will be accelerated up to energies above the pion photoproduction threshold on both the X-rays and the UV photons from the accretion disk. The secondary charged pions decay, producing neutrinos. We predict significant fluxes of TeV-PeV neutrinos, and show that the AMANDA II detector is already constraining several important astrophysical parameters of these sources. Larger cubic kilometer detectors such as IceCube will be able to detect such neutrinos in less than one year of operation, or otherwise rule out this scenario.

We investigate the high-energy photon emission around Gamma-Ray Burst (GRB) remnants caused by ultrahigh-energy cosmic rays (UHECRs) from the GRBs. We apply the results to the recent report that the supernova remnant W49B is a GRB remnant in our Galaxy. If this is correct, and if GRBs are sources of UHECRs, a natural consequence of this identification would be a detectable TeV photon emission around the GRB remnant. The imaging of the surrounding emission could provide new constraints on the jet structure of the GRB.

Gravitational waves and X-ray flares are expected from tidal disruption of stars by a massive black hole. Using a relativistic smoothed particle hydrodynamics code, we investigate the fate of main sequence and Helium stars in plunge orbits passing near Schwarzschild or a Kerr black holes of 10^5-6 Msun. We show that quadrupole gravitational waves emitted during the tidal disruption process are described reasonable well by a point particle approximation even in the strong encounter case. An additional hydrodynamic calculation based on the Godunov method indicates that shocks develop for sufficiently high tidal compressions. The shock-heating results in an X-ray flare, which for solar-type stars disrupted by 10^6Msun black holes is in the keV range, associated with the gravitational wave signal. The hardness and duration of the X-ray flare may serve as a diagnostic of the mass of the central black hole.

We have calculated the opacities and secondary production mechanisms of high energy photons arising in gamma-ray burst internal shocks, using exact cross-sections for the relevant processes. We find that for reasonable choices of parameters, photons in the range of 10's to 100's of GeV may be emitted in the prompt phase. Photons above this range are subject to electron-positron pair production with fireball photons and would be absent from the spectrum escaping the gamma-ray burst. We find that, in such cases, the fireball becomes optically thin again at ultra-high energies ($\gtrsim$ PeV). On the other hand, for sufficiently large fireball bulk Lorentz factors, the fireball is optically thin at all energies. Both for $\gamma\gamma$ self-absorbed and optically thin cases, the escaping high energy photons can interact with infra-red and microwave background photons to produce delayed secondary photons in the GeV-TeV range. These may be observable with GLAST, or at low redshifts with ground-based air Cherenkov telescopes. Detection of the primary prompt spectrum constrains the bulk Lorentz factor, while detection of delayed secondary gamma-rays would provide a consistency check for the primary spectrum and the bulk Lorentz factor as well as constraints on the intergalactic magnetic field strength.

We discuss the effects of particle acceleration at the bow shocks expected in the binary pulsar system J0737-3039, due to the wind from pulsar A interacting with both the interstellar medium (ISM), and with the magnetosphere of pulsar B. In this model, we find that the likeliest source for the X-rays observed by Chandra is the emission from the shocked wind of pulsar A as it interacts with the ISM. In this case, for favorable model parameter values, better statistics might help Chandra marginally resolve the source. A consequence of the model is a power law high energy spectrum extending up to $\lesssim$ 60 keV, at a level of $\sim 2 \times 10^{-13}$ erg cm$^{-2}$ s$^{-1}$.

The Swift mission, scheduled for launch in early 2004, is a multiwavelength observatory for gamma-ray burst (GRB) astronomy. It is the first-of-its-kind autonomous rapid-slewing satellite for transient astronomy and pioneers the way for future rapid-reaction and multiwavelength missions. It will be far more powerful than any previous GRB mission, observing more than 100 bursts per year and performing detailed X-ray and UV/optical afterglow observations spanning timescales from 1 minute to several days after the burst. The objectives are to determine the origin of GRBs; classify GRBs and search for new types; study the interaction of the ultra-relativistic outflows of GRBs with their surrounding medium; and use GRBs to study the early universe out to z>10. The mission is being developed by a NASA-led international collaboration. It will carry three instruments: a new-generation wide-field gamma-ray (15-150 keV) detector; a narrow-field X-ray telescope; and a narrow-field UV/optical telescope. Redshift determinations will be made for most bursts. In addition to the primary GRB science, the mission will perform a hard X-ray survey to a sensitivity of ~1 mCrab (~2x10^{-11} erg cm^{-2} s^{-1} in the 15-150 keV band), more than an order of magnitude better than HEAO A-4. A flexible data and operations system will allow rapid follow-up observations of all types of high-energy transients, with rapid data downlink and uplink available through the NASA TDRSS system. The mission is currently funded for 2 years of operations and the spacecraft will have a lifetime to orbital decay of ~8 years. [ABRIDGED]

Two separate topics are discussed. (1) We describe the classifications of the long GRB early afterglow lightcurves within the framework of the fireball shock model, focusing on the interplay between the reverse and forward shock emission components. We will also provide evidence that the central engine of at least two bursts are entrained with strong magnetic fields, and discuss the implications of this result for our understanding of the GRB phenomenon; (2) We argue that the current gamma-ray burst (GRB) and X-ray flash (XRF) data are consistent with a picture that all GRB-XRF jets are structured and quasi-universal, with a typical Gaussian-like jet structure.

An observed correlation Ep \propto (E_{iso})^{1/2} extending from Gamma-ray Bursts (GRBs) to X-ray flashes (XRFs) poses problems both for a power-law universal jet model where the energy per solid angle decreases as the inverse square of the angle respect to the jet axis, and for a conical jet model with a uniform energy density within the jet beam and a sharp energy cut-off at the jet edge. Here we show that the current GRB-XRF prompt emission/afterglow data can be understood in terms of a picture in which the GRB-XRF jets are quasi-universal and structured, with a Gaussian-like or similar structure, i.e., one where the jet has a characteristic angle, with a mild variation of energy inside and a rapid (e.g. exponential) decrease of energy outside of it. A Monte Carlo simulation shows that the current data is compatible with such a quasi-universal Gaussian jet with a typical opening angle of 5.7(+3.4)(-2.1) degrees, and with a standard jet energy of about log (E_j/1 erg)= 51.1+-0.3. According to this model, the true-to-observed number ratio of the whole GRB-XRF population is about 14 with the current instrumental sensitivity.

Gamma-ray bursts are promising tools for tracing the formation of high redshift stars, including the first generation. At very high redshifts the reverse shock emission lasts longer in the observer frame, and its importance for detection and analysis purposes relative to the forward shock increases. We consider two different models for the GRB environment, based on current ideas about the redshift dependence of gas properties in galaxies and primordial star formation. We calculate the observed flux as a function of the redshift and observer time for typical GRB afterglows, taking into account intergalactic photoionization and Lyman-$\alpha$ absorption opacity as well as extinction by the Milky Way Galaxy. The fluxes in the X-ray and near IR bands are compared with the sensitivity of different detectors such as Chandra, XMM, Swift XRT and JWST. Using standard assumptions, we find that Chandra, XMM and Swift XRT can potentially detect GRBs in the X-ray band out to very high redshifts $z\gtrsim$ 30. In the K and M bands, the JWST and ground-based telescopes are potentially able to detect GRBs even one day after the trigger out to $z\sim$ 16 and 33, if present. While the X-ray band is insensitive to the external density and to reverse shocks, the near IR bands provide a sensitive tool for diagnosing both the environment and the reverse shock component.

We discuss the effects of synchrotron self-absorption in the prompt emission from the reverse shock of GRB afterglows occurring in a dense environment, such as the wind of a massive stellar progenitor or a dense ISM in early galaxies. We point out that, when synchrotron losses dominate over inverse Compton losses, the higher self-absorption frequency in a dense environment implies a bump in the reverse shock emission spectrum, which can result in a more complex optical/IR light curve than previously thought. This bump is prominent especially if the burst ejecta is highly magnetized. In the opposite case of low magnetization, inverse Compton losses lead to a prompt X-ray flare. These effects give a possible new diagnostic for the magnetic energy density in the fireball, and for the presence of a dense environment.

The cosmological gamma-ray burst (GRB) phenomenon is reviewed. The broad observational facts and empirical phenomenological relations of the GRB prompt emission and afterglow are outlined. A well-tested, successful fireball shock model is introduced in a pedagogical manner. Several important uncertainties in the current understanding of the phenomenon are reviewed, and prospects of how future experiments and extensive observational and theoretical efforts may address these problems are discussed.

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.

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.

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.

Early optical afterglows have been observed from GRB 990123, GRB 021004, and GRB 021211, which reveal rich emission features attributed to reverse shocks. It is expected that Swift will discover many more early afterglows. Here we investigate in a unified manner both the forward and the reverse external shock emission components, and introduce a straightforward recipe for directly constraining the initial Lorentz factor of the fireball using early optical afterglow data. The scheme is largely independent of the shock microphysics. We identify two types of combinations of the reverse and forward shock emission, and explore their parameter regimes. We also discuss a possible diagnostic for magnetized ejecta. There is evidence that the central engine of GRB 990123 is strongly magnetized.

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.

We discuss the X-ray flux of gamma-ray burst afterglows at redshifts in the range 3-30, including the effects of the intergalactic He II absorption. We point out that strong X-ray lines may form locally in burst afterglows starting minutes after the trigger. This can provide distinctive X-ray distance indicators out to the redshifts where the first generation of massive stars form.

We present models for reprocessing of an intense flux of X-rays and gamma rays expected in the vicinity of gamma ray burst sources. We consider the transfer and reprocessing of the energetic photons into observable features in the X-ray band, notably the K lines of iron. Our models are based on the assumption that the gas is sufficiently dense to allow the microphysical processes to be in a steady state, thus allowing efficient line emission with modest reprocessing mass and elemental abundances ranging from solar to moderately enriched. We show that the reprocessing is enhanced by down-Comptonization of photons whose energy would otherwise be too high to absorb on iron, and that pair production can have an effect on enhancing the line production. Both "distant" reprocessors such as supernova or wind remnants and "nearby" reprocessors such as outer stellar envelopes can reproduce the observed line fluxes with Fe abundances 30-100 times above solar, depending on the incidence angle. The high incidence angles required arise naturally only in nearby models, which for plausible values can reach Fe line to continuum ratios close to the reported values.

Most of the current knowldege about GRB is based on electromagnetic observations at MeV and lower energies. Here we focus on some recent theoretical work on GRB, in particular the higher energy (GeV-TeV) photon emission, and two potentially important non-electromagnetic channels, the TeV and higher energy neutrino signals, and the gravitational wave signals expected from GRB.

We study the physics behind the neutrino-driven mechanism for gamma-ray bursts and hypernovae, deriving the critical density at which these outbursts occur in the collapsar model. The agreement between this derivation and results from past collapsar simulations (MacFadyen & Woosley 2000) is excellent, implying that we have captured the essential physics. We then use this derivation to study a range of progenitors for collapsar gamma-ray bursts. We derive how much of the star will accrete onto the black hole core before the infall density drops below this critical density, leading to an estimate of the remnant black hole mass for GRBs and hypernovae. We also estimate the time delays between gravity wave or neutrino signals and the onset of the explosion or burst event. This derivation, combined with future observational constraints, provides a physical insight into the structure of the GRB progenitor.

We present spectra of the optical transient of GRB021004 obtained with the Hobby-Eberly telescope starting 15.48, 20.31 hours, and 4.84 days after the burst and a spectrum obtained with the H. J. Smith 2.7 m Telescope starting 14.31 hours after the burst. GRB021004 is the first afterglow whose spectrum is dominated by absorption lines from high ionization species with multiple velocity components separated by up to 3000 km/s. We argue that these lines are likely to come from shells around a massive progenitor star. The high velocities and high ionizations arise from a combination of acceleration and flash-ionization by the burst photons and the wind velocity and steady ionization by the progenitor. We also analyze the broad-band spectrum and the light curve. We distinguish six components along the line of sight: (1) The z~2.293 absorption lines arise from the wind of a massive star. For a mass loss rate of ~6 x 10^{-5} solar masses per year, this component also provides the external medium to create the afterglow light. (2) A second shell produces absorption lines with a relative velocity of 560 km/s, and this is associated with the shell created by the fast massive star wind blowing a bubble in the preceding slow wind at a radial distance of order 10 pc. (3) More distant clouds within the host galaxy lie between 30-2500 pc, where they have been ionized by the burst. (4-6) The massive star wind has clumps with radii and over-densities of 0.022, 0.063, and 0.12 parsecs and 50%, 10%, and 10% respectively. The immediate progenitor of the burster could either be a WC-type Wolf-Rayet star or a highly evolved star whose original mass was just too small for it to become a WN-type Wolf-Rayet star.

We study gravitational radiation from various proposed gamma-ray burst (GRB) progenitor models, in particular compact mergers and massive stellar collapses. These models have in common a high angular rotation rate, and the final stage involves a rotating black hole and accretion disk system. We consider the in-spiral, merger and ringing phases, and for massive collapses we consider the possible effects of asymmetric collapse and break-up, as well bar-mode instabilities in the disks. We calculate the strain and frequency of the gravitational waves expected from various progenitors, at distances based on occurrence rate estimates. Based on simplifying assumptions, we give estimates of the probability of detection of gravitational waves by the advanced LIGO system from the different GRB scenarios.

Significant gravitational wave emission is expected from gamma-ray bursts arising from compact stellar mergers, and possibly also from bursts associated with fast-rotating massive stellar core collapses. These models have in common a high angular rotation rate, and observations provide evidence for jet collimation of the photon emission, with properties depending on the polar angle, which may also be of relevance for X-ray flashes. Here we consider the gravitational wave emission and its polarization as a function of angle which is expected from such sources. We discuss possible correlations between the burst photon luminosity, or the delay between gravitational wave bursts and X-ray flashes, and the polarization degree of the gravitational waves.

We argue that the distributions of both the intrinsic fluence and the intrinsic duration of the gamma-ray emission in gamma-ray bursts from the BATSE sample are well represented by log-normal distributions, in which the intrinsic dispersion is much larger than the cosmological time dilatation and redshift effects. We perform separate bivariate log-normal distribution fits to the BATSE short and long burst samples. The bivariate log-normal behaviour results in an ellipsoidal distribution, whose major axis determines an overall statistical relation between the fluence and the duration. We show that this fit provides evidence for a power-law dependence between the fluence and the duration, with a statistically significant different index for the long and short groups. We discuss possible biases, which might affect this result, and argue that the effect is probably real. This may provide a potentially useful constraint for models of long and short bursts.

Extensive observational campaigns of afterglow hunting have greatly enriched our understanding of the gamma-ray burst (GRB) phenomenon. Efforts have been made recently to explore some afterglow properties or signatures that will be tested by the on-going or the future observational campaigns yet come. These include the properties of GRB early afterglows in the temporal domain; the GeV-TeV afterglow signatures in the spectral domain; as well as a global view about the GRB universal structured jet configuration. These recent efforts are reviewed. Within the standard cosmological fireball model, the very model(s) responsible for the GRB prompt emission is (are) not identified. These models are critically reviewed and confronted with the current data.

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.

Typical gamma-ray burst spectra are characterized by a spectral break, Ep, which for bright BATSE bursts is found to be narrowly clustered around 300 keV. Recently identified X-ray flashes, which may account for a significant portion of the whole GRB population, seem to extend the Ep distribution to a broader range below 40 keV. Positive correlations among Ep and some other observed parameters have been noticed. On the other hand, within the cosmological fireball model, the issues concerning the dominant energy ingredient of the fireball as well as the location of the GRB emission site are still unsettled, leading to several variants of the fireball model. Here we analyze these models within a unified framework, and critically reexamine the Ep predictions in the various model variants. Attention is focused on the predictions of the narrowness of the Ep distribution in different models, and the correlations among Ep and some other measurable observables. These model properties may be tested against the current and upcoming GRB data, through which the nature of the fireball as well as the mechanism and site of the GRB emission will be identified. In view of the current data, various models are appraised through a simple Monte-Carlo simulation, and a tentative discussion about the possible nature of X-ray flashes is presented.

We investigate the relationship between the quasi-thermal baryon-related photosphere in relativistic outflows, and the internal shocks arising outside them, which out to a limiting radius may be able to create enough pairs to extend the optically thick region. Variable gamma-ray light curves are likely to arise outside this limiting pair-forming shock radius, while X-ray excess bursts may arise from shocks occurring below it; a possible relation to X-ray flashes is discussed. This model leads to a simple physical interpretation of the observational gamma-ray variability-luminosity relation.

We consider a gamma-ray burst (GRB) model based on an anisotropic fireball with an axisymmetric power-law distribution of the energy per solid angle with index -k, and allow for the observer's viewing direction being at an arbitrary angle with respect to the jet axis. This model can reproduce the key features expected from the conventional on-axis uniform jet models, with the novelty that the achromatic break time in the broadband afterglow lightcurves corresponds to the epoch when the relativistic beaming angle is equal to the viewing angle rather than to the jet half opening angle. If all the GRB fireballs have such a similar energy distribution form with 1.5 < k < (or \sim) 2, GRBs may be modeled by a quasi-universal beaming configuration, and an approximately standard energy reservoir. The conclusion also holds for some other forms of angular energy distributions, such as the Gaussian function.

The gamma ray burst phenomenon is reviewed from a theoretical point of view, with emphasis on the fireball shock scenario of the prompt emission and the longer wavelenght afterglow. Recent progress and issues are discussed, including spectral-temporal evolution, localizations, jets, spectral lines, environmental and cosmological aspects, as well as some prospects for future experiments in both electromagnetic and non-electromagnetic channels.

(Abridged) We discuss a Gamma-ray burst (GRB) fireball with an additional energy injection, either in the form of a Poynting-flux-dominated outflow or a kinetic-energy-dominated matter shell injected after the burst. A total injection energy comparable to that of the impulsive energy in the initial fireball is required to make a detectable signature in the afterglow lightcurves. For the case of a Poynting-flux-dominated injection, this leads to a gradual achromatic bump feature. For the case when the injection is kinetic-energy-dominated, the collision between the rear (injected) and the forward shell mild if the relative velocity between the colliding shells does not exceed a critical value defined by their energy ratio. Otherwise, the injection is violent, and an additional pair of shocks will form at the discontinuity between two colliding shells, so that there are altogether three shock-heated regions from which the emission contributes to the final lightcurves. We describe the shell-merging process in detail and present broadband lightcurves with the injection signatures. A violent collision results in very complicated features, which differs from the gradual bump signature found in the Poynting-flux or mild matter injection cases. In all the cases, the energetics of the fireball as well as the absolute afterglow flux level after the injection are boosted after the injection. Identifying the different injection signatures from future early afterglow observations may provide diagnostics about the nature of the fireball and of the central engine.

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.

In the collapsar scenario, gamma ray bursts are caused by relativistic jets expelled along the rotation axis of a collapsing stellar core. We discuss how the structure and time-dependence of such jets depends on the stellar envelope and central engine properties, assuming a steady jet injection. It takes a few seconds for the jet to bore its way through the stellar core; most of the energy output during that period goes into a cocoon of relativistic plasma surrounding the jet. This material subsequently forms a bubble of magnetized plasma that takes several hours to expand, subrelativistically, through the envelope of a high-mass supergiant. Jet break-through and a conventional burst would be expected not only in He stars but possibly also in blue supergiants. Shock waves and magnetic dissipation in the escaping bubble can contribute a non thermal UV/X-ray afterglow, and also excite Fe line emission from thermal gas, in addition to the standard jet deceleration power-law afterglow.

We investigate two high energy radiation mechanisms, the proton synchrotron and the electron inverse Compton emission, and explore their possible signatures in the broad-band spectra and in the keV to GeV light curves of gamma-ray burst afterglows. We develop a simple analytical approach, allowing also for the effects of photon-photon pair production, and explore the conditions under which one or the other of these components dominates. We identify three parameter space regions where different spectral components dominate: (I) a region where the proton synchrotron and other hadron-related emission components dominate, which is small; (II) a region where the electron inverse Compton component dominates, which is substantial; (III) a third substantial region where electron synchrotron dominates. We discuss the prospects and astrophysical implications of directly detecting the inverse Compton and the proton high energy components in various bands, in particular in the GeV band with future missions such as GLAST, and in the X-ray band with Chandra. We find that regime II parameter space is the most favorable regime for high energy emission. The inverse Compton component is detectable by GLAST within hours for bursts at typical cosmological distances, and by Chandra in days if the ambient density is high.

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 argue that the distributions of both the intrinsic fluence and the intrinsic duration of the gamma-ray emission in gamma-ray bursts from the BATSE sample are well represented by log-normal distributions, in which the intrinsic dispersion is much larger than the cosmological time dilatation and redshift effects. We perform separate bivariate log-normal distribution fits to the BATSE short and long burst samples. The bivariate log-normal behavior results in an ellipsoidal distribution, whose major axis determines an overall statistical relation between the fluence and the duration. We show that this fit provides evidence for a power-law dependence between the fluence and the duration, with a statistically significant different index for the long and short groups. We discuss possible biases which might affect this result, and argue that the effect is probably real. This may provide a potentially useful constraint for models of long and short bursts.

Gamma-ray burst sources with a high luminosity can produce electron-positron pair cascades in their environment as a result of back-scattering of a seed fraction of their original hard spectrum. The resulting spectral modifications offer the possibility of diagnosing not only the compactness parameter of the gamma-ray emitting region but also the baryonic density of the environment external to the burst.

We discuss three aspects of the interaction between GRB and their surroundings. The illumination of the progenitor remnant and/or the surroundings by the X-ray afterglow continuum can produce substantial Fe K-alpha line and edge emission, with implications for the progenitor model. The presence of large dust column densities, capable of obscuring the GRB optical afterglow, will lead to characteristic delayed X-ray and far-IR light curve signatures. Pair production induced by the initial gamma-rays in the nearby environment will modify the initial spectrum and the afterglow light curve, and the magnitude of these changes provides a diagnostic for the external density.

We investigate the consequences of a continuously injecting central engine on the gamma-ray burst afterglow emission, focusing more specifically on a highly-magnetized millisecond pulsar engine. For initial pulsar parameters within a certain region of the parameter space, the afterglow lightcurves are predicted to show a distinctive achromatic bump feature, the onset and duration of which range from minutes to months, depending on the pulsar and the fireball parameters. The detection of or upper limits on such features would provide constraints on the burst progenitor and on magnetar-like central engine models. An achromatic bump such as that in GRB 000301C afterglow may be caused by a millisecond pulsar with P0=3.4 millisecond and Bp=2.7e14 Gauss.

Gamma-ray bursts (GRB) are sudden, intense flashes of gamma-rays which, for a few blinding seconds, light up in an otherwise fairly dark gamma-ray sky. They are detected at the rate of about once a day, and while they are on, they outshine every other gamma-ray source in the sky, including the sun. Major advances have been made in the last three or four years, including the discovery of slowly fading x-ray, optical and radio afterglows of GRBs, the identification of host galaxies at cosmological distances, and finding evidence for many of them being associated with star forming regions and possibly supernovae. Progress has been made in understanding how the GRB and afterglow radiation arises in terms of a relativistic fireball shock model. These advances have opened new vistas and questions on the nature of the central engine, the identity of their progenitors, the effects of the environment, and their possible gravitational wave, cosmic ray and neutrino luminosity. The debates on these issues indicate that GRB remain among the most mysterious puzzles in astrophysics.

We report on the discovery of two emission features observed in the X-ray spectrum of the afterglow of the gamma-ray burst (GRB) of 16 Dec. 1999 by the Chandra X-Ray Observatory. These features are identified with the Ly$_{\alpha}$ line and the narrow recombination continuum by hydrogenic ions of iron at a redshift $z=1.00\pm0.02$, providing an unambiguous measurement of the distance of a GRB. Line width and intensity imply that the progenitor of the GRB was a massive star system that ejected, before the GRB event, $\approx 0.01 \Ms$ of iron at a velocity $\approx 0.1 c$, probably by a supernova explosion.

The recent report of X-ray Fe features in the afterglow of the gamma-ray burst GRB 991216 may provide important clues for identifying the nature of its progenitor and constraining the burst mechanism. We argue that the strong line emission can be attributed to the interaction of a continuing (but decaying) post-burst relativistic outflow from the central engine with the progenitor stellar envelope at distances less than a light-hour. Only a small mass of Fe is then required, which could have been readily produced by the star itself.

We calculate the synchrotron self-Compton emission from internal shocks occurring in relativistic winds as a source of gamma-ray bursts, with allowance for self-absorption. For plausible model parameters most pulses within a Gamma-Ray Burst (GRB) are optically thick to synchrotron self-absorption at the frequency at which most electrons radiate. Up-scattering of photon number spectra harder than $\nu^0$ (such as the self-absorbed emission) yields inverse Compton photon number spectra that are flat, therefore our model has the potential of explaining the low-energy indices harder than $\nu^{-2/3}$ (the optically thin synchrotron limit) that have been observed in some bursts. The optical counterparts of the model bursts are sufficiently bright to be detected by such experiments as LOTIS, unless the magnetic field is well below equipartition.

For Gamma-Ray Bursts occurring in dense star-forming regions, the X-ray afterglow behavior minutes to days after the trigger may be dominated by the small-angle scattering of the prompt X-ray emission off dust grains. We give a simple illustrative model for the X-ray light curves at different X-ray energies, and discuss possible implications. A bump followed by a steeper decay in soft X-rays is predicted for bursts which are heavily obscured in the optical.

Sub-photospheric internal shocks and transverse differences of the bulk Lorentz factor in relativistic fireball models of GRB lead to neutron diffusion relative to protons, resulting in inelastic nuclear collisions. This produces significant fluxes of ~3 GeV muon neutrinos (antineutrinos) and ~2 GeV electron neutrinos (antineutrinos), scaling with the Lorentz factor eta < 400. This extends significantly the parameter space for which neutrinos from inelastic collision are expected, which in the absence of the above effects requires eta > 400. A model with sideways diffusion of neutrons from a slower wind into a fast jet can lead to production of muon and electron neutrinos (antineutrinos) in the 2-25 GeV or higher range, depending on the value of eta. The emission from either of these mechanisms at z~1 may be detectable in suitably densely spaced detectors.

A gamma-ray burst fireball is likely to contain an admixture of neutrons, in addition to protons, in essentially all progenitor scenarios. Inelastic collisions between differentially streaming protons and neutrons in the fireball produce muon neutrinos (antineutrinos) of ~ 10 GeV as well as electron neutrinos (antineutrinos) of ~ 5 GeV, which could produce ~ 7 events/year in kilometer cube detectors, if the neutron abundance is comparable to that of protons. Photons of ~ 10 GeV from pi-zero decay and ~ 100 MeV electron antineutrinos from neutron decay are also produced, but will be difficult to detect. Photons with energies < 1 MeV from shocks following neutron decay produce a characteristic signal which may be distinguishable from the proton-related MeV photons.

MeV seed photons produced in shocks in a variable ultra-relativistic outflow gain energy by the Fermi mechanism, because the photons Compton scatter off relativistically colliding shells. The Fermi-modified high-energy photon spectrum has a non-universal slope and a universal cutoff. A significant increase in the total radiative efficiency is possible. In some gamma ray bursts, most of the power might be emitted at the high-energy cutoff for this mechanism, which would be close to 100 MeV for outflows with a mean bulk Lorentz factor of 100.

The standard model of Gamma-Ray Bursts afterglows is based on synchrotron radiation from a blast wave produced when the relativistic ejecta encounters the surrounding medium. We reanalyze the refreshed shock scenario, in which slower material catches up with the decelerating ejecta and reenergizes it. This energization can be done either continuously or in discrete episodes. We show that such scenario has two important implications. First there is an additional component coming from the reverse shock that goes into the energizing ejecta. This persists for as long as the re-energization itself, which could extend for up to days or longer. We find that during this time the overall spectral peak is found at the characteristic frequency of the reverse shock. Second, if the injection is continuous, the dynamics will be different from that in constant energy evolution, and will cause a slower decline of the observed fluxes. A simple test of the continuously refreshed scenario is that it predicts a spectral maximum in the far IR or mm range after a few days.

In a recent paper we have calculated the power density spectrum of Gamma-Ray Bursts arising from multiple shocks in a relativistic wind. The wind optical thickness is one of the factors to which the power spectrum is most sensitive, therefore we have further developed our model by taking into account the photon down-scattering on the cold electrons in the wind. For an almost optically thick wind we identify a combination of ejection features and wind parameters that yield bursts with an average power spectrum in agreement with the observations, and with an efficiency of converting the wind kinetic energy in 50-300 keV emission of order 1%. For the same set of model features the interval time between peaks and pulse fluences have distributions consistent with the log-normal distribution observed in real bursts.

The successful discovery of X-ray, optical and radio afterglows of GRB has made possible the identification of host galaxies at cosmological distances. The energy release inferred in these outbursts place them among the most energetic and violent events in the Universe. They are thought to be the outcome of a cataclysmic stellar collapse or compact stellar merger, leading to a relativistically expanding fireball, in which particles are accelerated at shocks and produce nonthermal radiation. The substantial agreement between observations and the theoretical predictions of the fireball shock model provide confirmation of the basic aspects of this scenario. Among recent issues are the collimation of the outflow and its implications for the energetics, the production of prompt bright flashes at wavelenghts much longer than gamma-rays, the time structure of the afterglow, its dependence on the central engine or progenitor system behavior, and the role of the environment on the afterglow.

Gamma-ray bursts are thought to be the outcome of a cataclysmic event leading to a relativistically expanding fireball, in which particles are accelerated at shocks and produce nonthermal radiation. We discuss the theoretical predictions of the fireball shock model and its general agreement with observations. Some of the recent work deals with the collimation of the outflow and its implications for the energe tics, the production of prompt bright flashes at wavelenghts much longer than gamma-rays, the time structure of the afterglow, its dependence on the central engine or progenitor system behavior, and the role of the environment on the evolution of the afterglow.

Steep X-ray spectral slopes, X-ray excesses and preferred spectral energy breaks in the 01-0.3 MeV range are discussed as the possible consequences of the photospheric component of the GRB relativistic outflow, and of pair breakdown in internal shocks leading to comptonization on semirelativistic electrons or MHD waves. We also discuss the X-ray and UV spectra of GRB afterglows occuring in a dense environment characteristic of massive stellar progenitors, including their ability to produce detectable Fe or other metallic line features.

We calculate the X-ray/UV spectral line signatures expected from the interaction of a gamma-ray burst afterglow and a dense pre-burst environment produced by the progenitor. We explore the conditions under which Fe line and edge equivalent widths of $\sim$ 1 keV can arise, and discuss the possibility of gaining information about possible progenitor scenarios using X-ray metal line spectra in the first few days of a burst. A wind or supernova shell around the burst produces an X-ray absorption line spectrum and later emission lines, while a hypernova funnel model produces mainly emission lines. The Fe \ked can in some cases be more prominent than the Fe \kal line. Under simple assumptions for the input continuum luminosity, current reports of observed Fe line luminosities are compatible with an Fe-enriched funnel model, while lower values are expected in shell models.

The role of a photospheric component and of pair breakdown is examined in the internal shock model of gamma-ray bursts. We discuss some of the mechanisms by which they would produce anomalously steep low energy slopes, X-ray excesses and preferred energy breaks. Sub-relativistic comptonization should dominate in high comoving luminosity bursts with high baryon load, while synchrotron radiation dominates the power law component in bursts which have lower comoving luminosity or have moderate to low baryon loads. A photosphere leading to steep low energy spectral slopes should be prominent in the lowest baryon load

We present an analytical approach to the dynamical evolution of fireballs or axially symmetric jets expanding into an external medium, with application to Gamma-Ray Burst remnants. This method leads to numerical calculations of fireball dynamics that are computationally faster than hydrodynamic simulations. It is also a very flexible approach, that can be easily extended to include more complex situations, such as a continuous injection of energy at the reverse shock, and the sideways expansion in non-spherical ejecta. Some features of the numerical results for the remnant dynamical evolution are discussed and compared to the analytical results. We find that the ratio of the observer times when the jet half-angle reaches twice its initial value and when the light cone becomes wider than the jet is substantially smaller than predicted analytically. The effects arising from the remnant curvature and thickness further reduce this ratio, such that the afterglow light-curve break due to the remnant finite angular extent overlaps the weaker break due to the jet's sideways expansion. An analysis of the effects of the curvature and thickness of the remnant on the afterglow light-curves shows that these effects are important and should be taken into account for accurate calculations of the afterglow emission.

We simulate Gamma-Ray Bursts arising from internal shocks in relativistic winds, calculate their power density spectrum (PDS), and identify the factors to which the PDS is most sensitive: the wind ejection features, which determine the wind dynamics and its optical thickness, and the energy release parameters, which give the pulse 50-300 keV radiative efficiency. For certain combinations of ejection features and wind parameters the resulting PDS exhibits the features observed in real bursts. We found that the upper limit on the efficiency of conversion of wind kinetic energy into 50-300 keV photons is $\sim$ 1%. Winds with a modulated Lorentz factor distribution of the ejecta yield PDSs in accord with current observations and have efficiencies closer to $10^{-3}$, while winds with a random, uniform Lorentz factor ejection must be optically thick to the short duration pulses to produce correct PDSs, and have an overall efficiency around $10^{-4}$.

Major advances have been made in the field of gamma-ray bursts in the last two years. The successful discovery of X-ray, optical and radio afterglows, which were predicted by theory, has made possible the identification of host galaxies at cosmological distances. The energy release inferred in these outbursts place them among the most energetic and violent events in the Universe. Current models envisage this to be the outcome of a cataclysmic event leading to a relativistically expanding fireball, in which particles are accelerated at shocks and produce nonthermal radiation. The substantial agreement between observations and the theoretical predictions of the standard fireball shock model provide confirmation of the basic aspects of this scenario. The continued observations show a diversity of behavior, providing valuable constraints for more detailed, post-standard models which incorporate more realistic physical features. Crucial questions being now addressed are the beaming at different energies and its implications for the energetics, the time structure of the afterglow, its dependence on the central engine or progenitor system behavior, and the role of the environment on the evolution of the afterglow.

The prompt $(t \siml 0.16$ days) light curve and initial 9-th magnitude optical flash from GRB 990123 can be attributed to a reverse external shock, or possibly to internal shocks. We discuss the time decay laws and spectral slopes expected under various dynamical regimes, and discuss the constraints imposed on the model by the observations, arguing that they provide strongly suggestive evidence for features beyond those in the simple standard model. The longer term afterglow behavior is discussed in the context of the forward shock, and it is argued that, if the steepening after three days is due to a jet geometry, this is likely to be due to jet-edge effects, rather than sideways expansion.

The discovery of X-ray, optical and radio afterglows of GRBs provides an important tool for understanding these sources. Most current models envisage GRB as arising in a cataclysmic stellar event leading to a relativistically expanding fireball, where particle acceleration at shocks lead to nonthermal radiation. The predictions of this scenario are in substantial agreement with the bulk of the observations. In addition, the data show a diversity of finer structure behavior, which is providing constraints for more detailed models. Current issues of interest are the implications of the beaming for the energetics, the afterglow time structure, dependence on progenitor system, and the role of the environment.

The successful discovery of X-ray, optical and radio afterglows of gamma-ray bursts has significantly helped our understanding of these sources, and made possible the identification of host galaxies at cosmological distances. The energy release inferred rivals that of supernovae, while its photon energy output may considerably exceed it. Current models envisage this to be the outcome of a cataclysmic stellar event leading to a relativistically expanding fireball, in which particles are accelerated at shocks and produce nonthermal radiation. The substantial agreement between observations and the theoretical predictions of the fireball shock model provide confirmation of the basic aspects of this scenario. The continued observations show a diversity of behavior, providing valuable constraints for more detailed models.

We investigate the spectrum of photohadronically produced neutrinos at very high energies (VHE, >10^14 eV) in astrophysical sources whose physical properties are constrained by their variability, in particular jets in Active Galactic Nuclei (blazars) and Gamma-Ray Bursts (GRBs). We discuss in detail the various competing cooling processes for energetic protons, as well as the cooling of pions and muons in the hadronic cascade, which impose limits on both the efficiency of neutrino production and the maximum neutrino energy. If the proton acceleration process is of the Fermi type, we can derive a model independent upper limit on the neutrino energy from the observed properties of any cosmic transient, which depends only on the assumed total energy of the transient. For standard energetic constraints, we can rule out major contributions above 10^19 eV from current models of both blazars and GRBs; and in most models much stronger limits apply in order to produce measurable neutrino fluxes. For GRBs, we show that the cooling of pions and muons in the hadronic cascade imposes the strongest limit on the neutrino energy, leading to cutoff energies of the electron and muon neutrino spectrum at the source differing by about one order of magnitude. We also discuss the relation of maximum cosmic ray energies to maximum neutrino energies and fluxes in GRBs, and find that the production of both the highest energy cosmic rays and observable neutrino fluxes at the same site can only be realized under extreme conditions; a test implication of this joint scenario would be the existence of strong fluxes of GRB correlated muon neutrinos up to ultra high energies, >10^17 eV. Secondary particle cooling also leads to slightly revised estimates for the neutrino fluxes from (non-transient) AGN cores.

We review the hypothesis that the acceleration of protons at internal shocks in Gamma Ray Bursts (GRB) could be the origin of the ultra-high energy cosmic rays (UHECR) observed at earth, E_max > 10^19 eV. We find that, even though protons may be accelerated to such energies, their ejection into the interstellar/intergalactic medium is problematic because it is likely to be accompanied by considerable adiabatic losses in the expanding shell. The problem is circumvented by neutrons produced in photohadronic interactions, which are not magnetically bound and thus effectively ejected in the moment of their production. They can be both produced in sufficient number and be able to leave the emission region if the optical depth of the emission region to photohadronic interactions is of order 1. We show that this requirement can be fulfilled under the same conditions which allow acceleration of protons to the highest energies. The production of neutrinos in this process correlates the fluxes of cosmic rays and neutrinos, and makes the hypothesis of UHECR origin in GRBs testable.

We simulate the dynamics of the interaction between relativistically expanding spherically symmetric shells using a 1-dimensional hydrodynamic code and calculate spectra and light-curves arising from such collisions by integrating the synchrotron and inverse Compton emission of the shocked gas. The numerical results reflect the most important features observed in Gamma-Ray Bursts: the spectrum exhibits a progressive softening (its break energy decaying exponentially with the 50--300 keV photon fluence), and the pulses that form the burst appear narrower in higher energy bands. Analytical results for the most important physical parameters of the burst are obtained by solving the shock jump conditions for a pair of interacting shells, in the case when both the forward and reverse shocks are relativistic.

A wide range of mechanisms have been proposed to supply the energy for gamma-ray bursts (GRB) at cosmological distances. It is a common misconception that some of these, notably NS-NS mergers, cannot meet the energy requirements suggested by recent observations. We show here that GRB energies, even at the most distant redshifts detected, are compatible with current binary merger or collapse scenarios involving compact objects. This is especially so if, as expected, there is a moderate amount of beaming, since current observations constrain the energy per solid angle much more strongly and directly than the total energy. All plausible progenitors, ranging from NS-NS mergers to various hypernova-like scenarios, eventually lead to the formation of a black hole with a debris torus around it, so that the extractable energy is of the same order, 1E+54 ergs, in all cases. MHD conversion of gravitational into kinetic and radiation energy can significantly increase the probability of observing large photon fluxes, although significant collimation may achieve the same effect with neutrino annihilation in short bursts. The lifetime of the debris torus is dictated by a variety of physical processes, such as viscous accretion and various instabilities; these mechanisms dominate at different stages in the evolution of the torus and provide for a range of gamma-ray burst lifetimes.

We discuss the formation of spectral features in the decelerating ejecta of gamma-ray bursts, including the possible effect of inhomogeneities. These should lead to blueshifted and broadened absorption edges and resonant features, especially from H and He. An external neutral ISM could produce detectable H and He, as well as Fe X-ray absorption edges and lines. Hypernova scenarios may be diagnosed by Fe K-$\alpha$ and H Ly-$\alpha$ emission lines.

Gamma ray burst outflows may entrain small blobs or filaments of dense, highly ionized metal rich material. Such inhomogeneities, accelerated by the flow to Lorentz factors in the range 10-100, could have a significant coverage factor, and give rise to broad features, especially due to Fe K-edges, which influence the spectrum below the MeV range, leading to a progressively decreasing hardness ratio.

We present an analytical approach to calculate the hydrodynamics of the interaction between a relativistic ejecta and a surrounding medium, whose evolution serves as a model for Gamma-Ray Burst afterglows. We investigate the effect of the relevant model parameters on the X-ray, optical and radio fluxes, and the effect of a refreshed shock energy input and anisotropy in the ejecta on the shape of the light curves. We compare our numerical results to observed afterglows and give a quantitative description of the conditions (geometry and physical parameters) in the ejecta that are compatible with the light curves of the 970508 afterglow, for which a large number of accurate flux measurements are available.

We present numerical simulations of Gamma-Ray Bursts arising from external shocks in the impulsive and wind models, including a weak or a strong coupling between electrons and protons plus magnetic fields, and analyze the burst features in each scenario. The dynamics of the ejecta and external medium are followed into the late stages of deceleration, in order to study the hydrodynamics of the remnant and the temporal and spectral evolution of the afterglow. A brief comparison with the optical and radio afterglows of GRB 970228 and GRB 970508 is made.

We consider fireball models where the ejecta have a range of bulk Lorentz factors, so that the inner (lower $\Gamma$) parts may carry most of the mass, or even most of the energy. The outer shock and contact discontinuity decelerate as the fireball sweeps up external matter. This deceleration allows slower ejecta to catch up, replenishing and reenergizing the reverse shock and boosting the momentum in the blast wave. In consequence, the energy available to power the afterglow may substantially exceed that of the burst itself. Such models allow a wide range of possibilities for the afterglow evolution, even in the case of spherically symmetric expansion.

We have carried out a principal component analysis for 625 gamma-ray bursts in the BATSE 3B catalog for which non-zero values exist for the nine measured variables. This shows that only two out of the three basic quantities of duration, peak flux and fluence are independent, even if this relation is strongly affected by instrumental effects, and these two account for 91.6% of the total information content. The next most important variable is the fluence in the fourth energy channel (at energies above 320 keV). This has a larger variance and is less correlated with the fluences in the remaining three channels than the latter correlate among themselves. Thus a separate consideration of the fourth channel, and increased attention on the related hardness ratio $H43$ appears useful for future studies. The analysis gives the weights for the individual measurements needed to define a single duration, peak flux and fluence. It also shows that, in logarithmic variables, the hardness ratio $H32$ is significantly correlated with peak flux, while $H43$ is significantly anticorrelated with peak flux. The principal component analysis provides a potentially useful tool for estimating the improvement in information content to be achieved by considering alternative variables or performing various corrections on available measurements

We discuss the afterglows from the evolution of both spherical and anisotropic fireballs decelerating in an inhomogeneous external medium. We consider both the radiative and adiabatic evolution regimes, and analyze the physical conditions under which these regimes can be used. Afterglows may be expected to differ widely among themselves, depending on the angular anisotropy of the fireball and the properties of the environment. They may be entirely absent, or may be detected without a corresponding $\gamma$-ray event. A tabulation of different representative light curves is presented, covering a wide range of behaviors that resemble what is currently observed in GRB 970228, GRB 970508 and other objects.

Models of gamma ray bursts are reviewed in the light of recent observations of afterglows which point towards a cosmological origin. The physics of fireball shock models is discussed, with attention to the type of light histories and spectra during the gamma-ray phase. The evolution of the remnants and their afterglows is considered, as well as their implications for our current understanding of the mechanisms giving rise to the bursts.

We derive the equation for the surface of equal arrival times of radiation from a thin relativistic shell interacting with an external medium, representing the afterglow of a gamma-ray burst produced by a fireball. Due to the deceleration, these surfaces become distorted ellipsoids and, at sufficiently late times, most of the light (either bolometric or in a given band) comes from a ring-like region whose width depends only on age. We analyze the shape of these surfaces and the radiation received from different angles for different dynamic and radiative regimes and homogeneous or power-law external densities. We calculate angle-integrated bolometric and fixed frequency fluxes, and we tabulate the most relevant parameters that describe the equal arrival time surfaces and the source brightness distribution, quantities that are useful for more accurate analytic estimates of the afterglow evolution.

We compute burst spectra and time structures arising from synchrotron and inverse Compton scattering by non-thermal electrons accelerated in shocks which form during the interaction between a thin ultra-relativistic fireball and a stationary external medium. We investigate the effect of varying the most important model parameters on the resulting burst spectra, and we present a set of correlations among the spectral and temporal features of the bursts. The spectral hardness, various spectral-temporal correlations and the spectral evolution of the simulated bursts are compared to those of observed bursts for a representative set of model parameters. Multi-pulse structures are simulated using a variable magnetic field and anisotropic emission, and the most important spectral and temporal properties of these pulses are compared with observations.

The location accuracy of the BeppoSAX Wide Field Cameras and acute ground-based followup have led to the detection of a decaying afterglow in X rays and optical light following the classical gamma-ray burst GRB 970228. The afterglow in X rays and optical light fades as a power law at all wavelengths. This behaviour was predicted for a relativistic blast wave that radiates its energy when it decelerates by ploughing into the surrounding medium. Because the afterglow has continued with unchanged behaviour for more than a month, its total energy must be of order 10**51 erg, placing it firmly at a redshift of order 1. Further tests of the model are discussed, some of which can be done with available data, and implications for future observing strategies are pointed out. We discuss how the afterglow can provide a probe for the nature of the burst sources.

We numerically model the interaction between an expanding fireball and a stationary external medium whose density is either homogeneous or varies with distance as a power-law. The evolution is followed until most of the fireball kinetic energy is converted into internal energy. The density, pressure and flow Lorentz factor profiles are shown at different stages, including shock and rarefaction wave reflections, for a fireball of initial bulk Lorentz factor Gamma = 100, both in the adiabatic and non-adiabatic (radiative) regimes. For cooling times shorter than the dynamic time, bolometric light-curves are computed for values of Gamma = 50, 100 and 200. We compare the numerical light-curves with analytic results, and find that for a homogeneous external medium there is a simple scaling relationship among light-curves obtained for different parameters. The light-curves for power-law external densities are similar in shape to those in the homogeneous case. We discuss the implications of a comparison of the results with observed Gamma-Ray Burst time histories.

We discuss the properties of magnetically dominated jet-like outflows from stellar mass black holes surrounded by debris tori resulting from neutron star disruption. These jets may have narrow cores (along the rotation axis) which are almost free of baryons and attain very high bulk Lorentz factors $\simg 10^6$. The jets give rise to a characteristic MeV to TeV emission as well as to relativistic shocks producing the usual MeV bursts. Because the outflow is highly directional the properties of the observed gamma-rays would depend on the viewing angle relative to the rotation axis. Even for the most intense bursts, which under the assumption of isotropic emission and substantial redshifts would be inferred to emit $10^{52}-10^{53}$ erg, the efficiencies required are only $10^{-2}-10^{-4}$.

Two models of the gamma ray burst population, one with a standard candle luminosity and one with a power law luminosity distribution, are chi^2-fitted to the union of two data sets: the differential number versus peak flux distribution of BATSE's long duration bursts, and the time dilation and energy shifting versus peak flux information of pulse duration time dilation factors, interpulse duration time dilation factors, and peak energy shifting factors. The differential peak flux distribution is corrected for threshold effects at low peak fluxes and at short burst durations, and the pulse duration time dilation factors are also corrected for energy stretching and similar effects. Within an Einstein-de Sitter cosmology, we place strong bounds on the evolution of the bursts, and these bounds are incompatible with a homogeneous population, assuming a power law spectrum and no luminosity evolution. Additionally, under the implied conditions of moderate evolution, the 90% width of the observed luminosity distribution is shown to be < 10^2, which is less constrained than others have demonstrated it to be assuming no evolution. Finally, redshift considerations indicate that if the redshifts of BATSE's faintest bursts are to be compatible with that which is currently known for galaxies, a standard candle luminosity is unacceptable, and in the case of the power law luminosity distribution, a mean luminosity < 10^57 ph s^-1 is favored.

We calculate the spectra expected from unsteady relativistic wind models of gamma-ray bursts, suitable for events of arbitrary duration. The spectral energy distribution of the burst is calculated over photon energies spanning from eV to TeV, for a range of event durations and variability timescales. The relative strength of the emission at different wavelengths can provide valuable information on the particle acceleration, radiation mechanisms and the possible types of models.

We discuss the evolution of cosmological gamma-ray burst remnants, consisting of the cooling and expanding fireball ejecta together with any swept-up external matter, after the gamma-ray event. We show that significant optical emission is predicted which should be measurable for timescales of hours after the event, and in some cases radio emission may be expected days to weeks after the event. The flux at optical, X-ray and other long wavelengths decays as a power of time, and the initial value of the flux or magnitude, as well as the value of the time-decay exponent, should help to distinguish between possible types of dissipative fireball models.

In energetic nonthermal sources such as gamma-ray bursts, AGN or galactic jet sources, etc., one expects both relativistic and transrelativistic shocks acompanied by violent motions of moderately relativistic plasma. We present general considerations indicating that these sites are electron and positron accelerators leading to a modified power law spectrum. The electron (or $e^\pm$) energy index is very hard, $\propto \gamma^{-1}$ or flatter up to a comoving frame break energy $\gamma_\ast$, and becomes steeper above that. In the example of gamma-ray bursts the Lorentz factor reaches $\gamma_\ast\sim 10^3$ for $e^{\pm}$ accelerated by the internal shock ensemble on subhydrodynamical time scales. For pairs accelerated on hydrodynamical timescales in the external shocks similarly hard spectra are obtained, and the break Lorentz factor can be as high as $\gamma_\star \siml 10^5$. Radiation from the nonthermal electrons produces photon spectra with shape and characteristic energies in qualitative agreement with observed generic gamma-ray burst and blazar spectra. The scenario described here provides a plausible way to solve one of the crucial problems of nonthermal high energy sources, namely the efficient transfer of energy from the proton flow to an apropriate nonthermal lepton component.

We present analytic formulae for the integral number count distribution of cosmological bursting or steady sources valid over the entire range of fluxes, including density evolution and either standard candle or a power law luminosity function. These are used to derive analytic formulae for the mean redshift, the time dilations and the dispersion of these quantities for sources within a given flux range for Friedmann models with $\Omega=1,~ \Lambda=0$ without K-corrections, and we discuss the extension to cases with $\Omega <1$ and inclusion of K-corrections. Applications to the spatial distribution of cosmological gamma ray burst sources are discussed, both with and without an intrinsic energy stretching of the burst time profiles, and the implied ranges of redshift $z$ are considered for a specific time dilation signal value. The simultaneous consideration of time dilation information and of fits of the number distribution versus peak flux breaks the degeneracy inherent in the latter alone, allowing a unique determination of the density evolution index and the characteristic luminosity of the sources. For a reported time dilation signal of 2.25 and neglecting [including] energy stretching we find that the proper density should evolve more steeply with redshift than comoving constant, and the redshifts of the dimmest sources with stretching would be very large. However, the expected statistical dispersion in the redshifts is large, especially for power law luminosity functions, and remains compatible with that of distant quasars. For smaller time dilation values of 1.75 and 1.35 the redshifts are more compatible with conventional ideas about galaxy formation, and the evolution is closer to a comoving constant or a slower evolution. More generally, we have considered a wide range of possible measured time dilation

We discuss detailed fits of the BATSE and PVO gamma-ray burst peak-flux distributions with Friedman models taking into account possible density evolution and standard candle or power law luminosity functions. A chi-square analysis is used to estimate the goodness of the fits and we derive the significance level of limits on the density evolution and luminosity function parameters. Cosmological models provide a good fit over a range of parameter space which is physically reasonable.

We present analytical solutions for the integral distribution of arbitrary bursting or steady source counts as a function of peak photon count rate within Friedmann cosmological models. We discuss both the standard candle and truncated power-law luminosity function cases with a power-law density evolution. While the analysis is quite general, the specific example discussed here is that of a cosmological gamma-ray burst distribution. These solutions show quantitatively the degree of dependence of the counts on the density and luminosity function parameters, as well as the the weak dependence on the closure parameter and the maximum redshift. An approximate comparison with the publicly available Compton Gamma Ray Observatory data gives an estimate of the maximum source luminosity and an upper limit to the minimum luminosity. We discuss possible ways of further constraining the various parameters.

Whatever the ultimate energy source of gamma-ray bursts turns out to be, the resulting sequence of physical events is likely to lead to a fairly generic, almost unavoidable scenario: a relativistic fireball that dissipates its energy after it has become optically thin. This is expected both for cosmological and halo distances. Here we explore the observational motivation of this scenario, and the consequences of the resulting models for the photon production in different wavebands, the energetics and the time structure of classical gamma-ray bursters.

Sudden collapse of a compact object, or coalescence of a compact binary, can generate an unsteady relativistic wind that lasts for a few seconds. The wind is likely to carry a high magnetic field; and its Lorentz factor depends on the extent to which it is 'loaded' with baryons. If the Lorentz factor is $\sim 100$, internal dissipation and shocks in this wind produce a non-thermal gamma-ray burst, detectable in the range $0.1\MeV \siml E_\gamma \siml 0.1-1\GeV$ out to cosmological distances. The cooled wind ejecta would subsequently be decelerated by the external medium. The resultant blast wave and reverse shock can then give rise to a second burst component, mainly detectable in the GeV range, with a time delay relative to the MeV burst ranging from minutes to hours.

The 'event' that triggers a gamma ray burst cannot last for more than a few seconds. This is, however, long compared with the dynamical timescale of a compact stellar-mass object ($\sim 10^{-3}$ seconds). Energy is assumed to be released as an outflow with high mean lorentz factor $\Gamma$. But a compact stellar-mass collapse or merger is, realistically, likely to generate a mass (or energy) flux that is unsteady on some timescales in the range $10^{-3}$ - 10 seconds. If $\Gamma$ fluctuates by a factor of $\sim 2$ around its mean value, relative motions within the outflowing material will themselves (in the comoving frame) be relativistic, and can give rise to internal shocks. For $\Gamma \sim 10^2$, the resultant dissipation occurs outside the 'photosphere' and can convert a substantial fraction of the overall outflow energy into non-thermal radiation. This suggests a mechanism for cosmological bursts that demands less extreme assumptions (in respect of $\Gamma$-values, freedom from baryonic contamination, etc) than earlier proposals.

We calculate the spectrum of blast wave models of gamma-ray burst sources, for various assumptions about the magnetic field density and the relativistic particle acceleration efficiency. For a range of physically plausible models we find that the radiation efficiency is high, and leads to nonthermal spectra with breaks at various energies comparable to those observed in the gamma-ray range. Radiation is also predicted at other wavebands, in particular at X-ray, optical/UV and GeV/TeV energies. We discuss the spectra as a function of duration for three basic types of models, and for cosmological, halo and galactic disk distances. We also evaluate the gamma-ray fluences and the spectral characteristics for a range of external densities. Impulsive burst models at cosmological distances can satisfy the conventional X-ray paucity constraint $S_x/S_\gamma \siml$ few percent over a wide range of durations, but galactic models can do so only for bursts shorter than a few seconds, unless additional assumptions are made. The emissivity is generally larger for bursts in a denser external environment, with the efficiency increasing up to the point where all the energy input is radiated away.

If gamma-ray bursts originate in our Galaxy, they probably involve violent disturbances in the magnetospheres of neutron stars. Any event of this kind is likely to trigger the sudden expulsion of magnetic flux and plasma at relativistic speed. We show that such ejecta would be braked by the interstellar medium (ISM), and that a gamma-ray flash with duration $\sim 0.1 - 100 $ s may result from this interaction. We suggest that the stronger,``smooth"class of bursts, recently identified from analyses of BATSE data, could result from the deceleration of such a blast wave by the ISM. The radiative efficiency,and hence the detectability, of such bursts would depend on the density of the circumstellar ISM. Therefore, even if the neutron star ``triggers" were uniformly distributed in space (at least within 1--2 kpc of the Sun), the bserved locations of bursts would correlate with regions of above-average ISM density.

In almost any scenario for 'cosmological' gamma-ray bursts (and in many models where they originate in our own Galaxy), the initial energy density is so large that the resulting relativistic plasma expands with $v\sim c$ producing a blast wave ahead of it and a reverse shock moving into the ejecta, as it ploughs into the external medium. We evaluate the radiation expected from these shocks,for both cosmological and galactic bursts, for various assumptions about the strength of the magnetic field and the particle acceleration mechanisms in the shocks. The spectra are evaluated over the whole range from the IR to $>$ GeV, and are compared with the variety of spectral behavior reported by BATSE, and with the X-ray and optical constraints. For bursts of duration $\simg 1\s$ acceptable $\gamma$-ray spectra and $L_x/L_\gamma$ ratios are readily obtained for 'cosmological' models. Blast waves in galactic models can produce bursts of similar gamma-ray fluence and duration, but they violate the X-ray paucity constraint, except for the shorter bursts ($\siml 1\s$). We discuss the prospects for using O/UV and X-ray observations to discriminate among alternative models.