Cosmic ray (CR) sources leave signatures in the isotopic abundances of CRs. Current models of Galactic CRs that consider supernovae (SNe) shocks as the main sites of particle acceleration cannot satisfactorily explain the higher ${\rm ^{22}Ne/^{20}Ne}$ ratio in CRs compared to the interstellar medium. Although stellar winds from massive stars have been invoked, their contribution relative to SNe ejecta has been taken as a free parameter. Here we present a theoretical calculation of the relative contributions of wind termination shocks (WTSs) and SNe shocks in superbubbles, based on the hydrodynamics of winds in clusters, the standard stellar mass function, and stellar evolution theory. We find that the contribution of WTSs towards the total CR production is at least $25\%$, which rises to $\gtrsim 50\%$ for young ($\lesssim 10$ Myr) clusters, and explains the observed $^{22}{\rm Ne}/^{20} {\rm Ne}$ ratio. We argue that since the progenitors of apparently isolated supernovae remnants (SNRs) are born in massive star clusters, both WTS and SNe shocks can be integrated into a combined scenario of CRs being accelerated in massive clusters. This scenario is consistent with the observed ratio of SNRs to $\gamma$-ray bright ($L_\gamma \gtrsim 10^{35}$ erg s$^{-1}$) star clusters, as predicted by star cluster mass function. Moreover, WTSs can accelerate CRs to PeV energies, and solve other longstanding problems of the standard supernova paradigm of CR acceleration.

Neutron star mergers, referring to both binary neutron star and neutron star black hole mergers, are the canonical multimessenger events. They have been detected across the electromagnetic spectrum, have recently been detected in gravitational waves, and are likely to produce neutrinos over several decades in energy. The non-thermal prompt and afterglow emission of short gamma-ray bursts and the quasi-thermal emission from the radioactively powered kilonovae provide distinct insights into the physics of neutron star mergers. When combined with direct information on coalescence from gravitational waves and neutrinos these sources may become the best understood astrophysical transients. Multimessenger observations of these cataclysmic events will determine sources of gravitational waves and astrophysical neutrinos, enable precision cosmology, and unique tests of fundamental physics, the origin of heavy elements, the behavior of relativistic jets, and the equation of state of supranuclear matter. In this white paper we present a summary of the science discoveries possible with multimessenger observations of neutron star mergers and provide recommendations to enable them in the new era of time-domain, multimessenger astronomy.

With upcoming (continuum) surveys of high-resolution radio telescopes, detection rates of fast radio bursts (FRBs) might approach $10^5$ per sky per day by future extremely large observatories, such as the possible extension of the Square Kilometer Array (SKA) to a phase 2 array. Depending on the redshift distribution of FRBs and using the repeating FRB121102 as a model, we calculate a detection rate of multiply-imaged FRBs with their multiply-imaged hosts caused by the distribution of galaxy-cluster scale gravitational lenses of the order of $10^{-4}$ per square degree per year for a minimum total flux of the host of 10 $\mu$Jy at 1.4 GHz for SKA phase 2. Our comparison of estimated detection rates for quasars, supernovae, gamma ray bursts, and FRBs shows that multiple images of FRBs could be more numerous than those of gamma ray bursts and supernovae and as numerous as multiple images of quasars. Time delays between the multiple images of an FRB break degeneracies in model-based and model-independent lens reconstructions as other time-varying sources do, yet without a microlensing bias as FRBs are more point-like and have shorter duration times. We estimate the relative imprecision of FRB time-delay measurements to be $10^{-10}$ for time delays on the order of 100 days for galaxy-cluster scale lenses, yielding more precise (local) lens properties than time delays from the other time-varying sources. Using the lens modelling software Grale, we show the increase in accuracy and precision of the reconstructed scaled surface mass density map of a simulated cluster-scale lens when adding time delays for one set of multiple images to the set of observational constraints.

The variation in the intensity of cosmic rays at small angular scales is attributed to the interstellar turbulence in the vicinity of the Solar system. We show that {a turbulent origin of the small-scale structures implies that} the morphology of the observed cosmic-ray intensity skymap varies with our location in the interstellar turbulence. The gyroradius of cosmic rays is shown to be the length scale associated with an observable change in the skymap over a radian angular scale. The extent to which the intensity \mpo{at a certain} angular scale varies is proportional to the change in our location with a maximum change of about the amplitude of intensity variation at that scale in the existing skymap.} We suggest that for TeV cosmic rays a measurable variation could occur over a time scale of a decade due to the Earth's motion through the interstellar medium, if interstellar turbulence persists down to the gyroradius, \mpo{about $00\ \mu\mathrm{pc}$ for TeV-ish cosmic rays}. Observational evidence of the variability, or an absence of it, could provide a useful insight into the physical origin of the small-scale anisotropy.

The recent radio observations (Mooley et al, 2018) of a superluminal radio afterglow following GRB 170817A are interpreted in terms of a jet impacting a baryonic cloak, which is presumably the material caught at the front of the jet as the latter emerges from a denser ejected material. Assuming that we the observers are located {\bf at a viewing angle of $\sim 0.2$ radians from the emitting material (perhaps slightly more from jet axis)}, we suggest that the Lorentz factor of the jet is $\lesssim 20$ at the time of the prompt emission, and that, as suggested previously, it is accelerated to much higher values before finally decelerating during the afterglow phase. A less extreme example of a short GRB being observed off axis may have been GRB 150101b (Fong, et al., 2016). A feature of GRBs viewed from large offset angles is a large afterglow isotropic equivalent energy as compared to prompt emission, as predicted (Eichler, 2017), and this is born out by the observations of these two GRB. It is also shown that the prompt emission of GRB 170817A, if seen way off-axis ($\theta \gg 1/\Gamma $), could not be made by internal shocks in the baryonic material that powers the afterglow.

It is noted that the duration of a fast radio burst (FRB), about $10^{-3}$ s, is a smaller fraction of the time delay between multiple images of a source gravitationally lensed by a galaxy or galaxy cluster than the human lifetime is to the age of the universe. Thus repeating, strongly lensed FRBs may offer an unprecedented opportunity for observing cosmological evolution in "real time". The possibility is discussed of observing cosmic expansion, transverse proper motion, mass accretion and perhaps growth of density perturbations, as a function of redshift.

It is shown that the antiproton spectrum reported by the AMS02 collaboration can be accounted for by dark matter (DM) decay if the residence time in the Galactic halo is of order 90 Myr. The DM lifetime assumed, $ 5 \cdot 10^{27}$ s for DM particle mass of 3 TeV, is shown to be consistent with the constraints of diffuse gamma ray background. Alternative sources of antiprotons, such as those from compact objects, are shown to be strongly constrained by neutrino astronomy and therefore unlikely candidates.

It has been suggested that short gamma ray bursts (GRBs) have shorter or undetectable spectral lags than long GRB because in the former, the observer's line of sight makes a larger angle with the GRB jet axis than for the latter. It is proposed that simultaneous gravitational wave - short GRB events could provide a simple test of this hypothesis. Multimessenger astronomy eventually may test whether event horizons are a necessary ingredient for GRB.

It is suggested that fast radio bursts can probe gravitational lensing by clumpy dark matter objects that range in mass from $10^{-3}M_{\odot}$ to $10^2 M_{\odot}$. They may provide a more sensitive probe than observations of lensing of objects in the Magellanic Clouds, and could find or rule out clumpy dark matter with an extended mass spectrum.

It is suggested that the decline with energy of the boron to carbon abundance ratio in Galactic cosmic rays is due, in part, to a correlation between the maximum energy attainable by shock acceleration in a given region of the Galactic disk and the grammage traversed before escape. In this case the energy dependence of the escape rate from the Galaxy may be less than previously thought and the spectrum of antiprotons becomes easier to understand.

Using analytic methods and $1$-D two-fluid simulations, we study the effect of cosmic rays (CRs) on the dynamics of interstellar superbubbles (ISBs) driven by multiple supernovae (SNe)/stellar winds in OB associations. In addition to CR advection and diffusion, our models include thermal conduction and radiative cooling. We find that CR injection at the reverse shock or within a central wind-driving region can affect the thermal profiles of ISBs and hence their X-ray properties. Even if a small fraction ($10-20\%$) of the total mechanical power is injected into CRs, a significant fraction of the ram pressure at the reverse shock can be transferred to CRs. The energy transfer becomes efficient if (1) the reverse shock gas Mach number exceeds a critical value ($M_{\rm th}\gtrsim 12$) and (2) the CR acceleration time scale $\tau_{\rm acc}\sim \kappa_{\rm cr}/v^2$ is shorter than the dynamical time, where $\kappa_{\rm cr}$ is CR diffusion constant and $v$ is the upstream velocity. We show that CR affected bubbles can exhibit a volume averaged hot gas temperature $1-5\times10^{6}$ K, lower by a factor of $2-10$ than without CRs. Thus CRs can potentially solve the long-standing problem of the observed low ISB temperatures.

We simulate decaying turbulence in a homogeneous pair plasma using three dimensional electromagnetic particle-in-cell (PIC) method. A uniform background magnetic field permeates the plasma such that the magnetic pressure is three times larger than the thermal pressure and the turbulence is generated by counter-propagating shear Alfv\'en waves. The energy predominately cascades transverse to the background magnetic field, rendering the turbulence anisotropic at smaller scales. We simultaneously move several ion species of varying charge to mass ratios in our simulation and show that the particles of smaller charge to mass ratios are heated and accelerated to non-thermal energies at a faster rate, in accordance with the enhancement of heavy ions and non-thermal tail in their energy spectrum observed in the impulsive solar flares. We further show that the heavy ions are energized mostly in the direction perpendicular to the background magnetic field with a rate consistent with our analytical estimate of the rate of heating due to cyclotron resonance with the Alfv\'en waves of which a large fraction is due to obliquely propagating waves.

Particle acceleration by cascading Alfven wave turbulence was suggested (Eichler, 1979b) as being responsible for energetic particle populations in $^3$He-rich solar flares. In particular, it was noted that the damping of the turbulence by the tail of the particle distribution in rigidity naturally leads to dramatic enhancement of pre-accelerated species - as $^3$He is posited to be - and superheavy elements. The subsequent detection of large enrichment of ultraheavies, relative to iron, has apparently confirmed this prediction, lending support to the original idea. It is shown here that this picture could be somewhat sharpened by progress in understanding the 3-dimensional geometrical details of cascading Alfven turbulence (Sridhar and Goldreich, 1995). The mechanism may be relevant in other astrophysical environments where the source of turbulence is non-magnetic, such as clusters of galaxies.

It is suggested that essentially all of the UHECRs we detect, including those at the highest energy, originate in our Galaxy. It is shown that even if the density of sources decreases with Galactic radius, then the anisotropy and composition can be understood. Inward anisotropy, as recently reported by the Auger collaboration can be understood as drift along the current sheet of UHECRs originating outside the solar circle, as predicted in Kumar and Eichler (2014), while those originating within the solar circle exit the Galaxy at high latitudes.

The Fermi Large Area Telescope (Fermi LAT) collaboration has recently reported on 50 months of measurements of the isotropic Extragalactic Gamma Ray Background (EGRB) spectrum between $100\mbox{MeV}$ and $820\mbox{GeV}$. Ultrahigh Energy Cosmic Ray (UHECR) protons interact with the Cosmic Microwave Background (CMB) photons and produce cascade photons of energies $10\mbox{MeV}\mbox{-}1\mbox{TeV}$ that contribute to the EGRB flux. We examine seven possible evolution models for UHECRs and find that UHECR sources that evolve as the Star Formation Rate (SFR), medium low luminosity active galactic nuclei type-1 ($L = 10^{43.5}\mbox{erg sec}^{-1}$ in the $[0.5\mbox{-}2]\mbox{KeV}$ band), and BL Lacertae objects (BL Lacs) are the most acceptable given the constraints imposed by the observed EGRB. Other possibilities produce too much secondary $\gamma$-radiation. In all cases, the decaying dark matter contribution improves the fit at high energy, but the contribution of still unresolved blazars, which would leave the smallest role for decaying dark matter, may yet provide an alternative improvement. The possibility that the entire EGRB can be fitted with resolvable but not-yet-resolved blazars, as recently claimed by Ajello et al. (2015), would leave little room in the EGRB to accommodate $\gamma$-rays from extragalactic UHECR production, even for many source evolution rates that would otherwise be acceptable. We find that, under the assumption of UHECRs being mostly protons, there is not enough room for producing extragalactic UHECRs with AGN, GRB, or even SFR source evolution. Sources that evolve as BL Lacs on the other hand, would produce much less secondary $\gamma$-radiation and would remain a viable source of UHECRs, provided that they dominate.

The dynamics of two initially unmagnetized relativistic counter-streaming homogeneous ion-electron plasma beams are simulated in two dimensions using the particle-in-cell (PIC) method. It is shown that current filaments, which form due to the Weibel instability, develop a large scale longitudinal electric field in the direction opposite to the current carried by the filaments as predicted by theory. Fast moving ions in the current filaments decelerate due to this longitudinal electric field. The same longitudinal electric field, which is partially inductive and partially electrostatic, is identified as the main source of acceleration of electrons in the current filaments. The transverse electric field, though larger than the longitudinal one, is shown to play a smaller role in heating electrons, contrary to previous claims. It is found that, in 1D, the electrons become strongly magnetized and are \textit{not} accelerated beyond their initial kinetic energy. Rather, the heating of the electrons is enhanced by the bending and break-up of the filaments, which releases electrons that would otherwise be trapped within a single filament and hence slow the development of the Weibel instability (i.e. the magnetic field growth) via induction as per Lenz's law. In 2D simulations electrons are heated to about one quarter of the initial kinetic energy of ions. The magnetic energy at maximum is about 4 percent, decaying to less than 1 percent by the end of the simulation. Most of the heating of electrons takes place while the longitudinal electric field is still growing while only a small portion of the heating is a result of subsequent magnetic field decay. The ions are found to gradually decelerate until the end of the simulation by which time they retain a residual anisotropy less than 10 percent.

It is suggested that many $\gamma$-ray bursts (GRBs) are cloaked by an ultra-relativistic baryonic shell that has high optical depth when the photons are manufactured. Such a shell would not fully block photons reflected or emitted from its inner surface, because the radial velocity of the photons can be less than that of the shell. This avoids the standard problem associated with GRBs that the thermal component should be produced where the flow is still obscured by high optical depth. The radiation that escapes high optical depth obeys the Amati relation. Observational implications may include a) anomalously high ratios of afterglow to prompt emission, such as may have been the case in the recently discovered PTF 11agg, and b) ultrahigh-energy neutrino pulses that are non-coincident with detectable GRB. It is suggested that GRB 090510, a short, very hard GRB with very little afterglow, was an {\it exposed} GRB, in contrast to those cloaked by baryonic shells. \end{abstract}

We study the time dependent propagation of sub-ankle ultra-high energy cosmic rays (UHECRs) originating from point-like Galactic sources. We show that drift in the Galactic magnetic field (GMF) may play an important role in the propagation of UHECRs and their measured anisotropy, particularly when the transport is anisotropic. To fully account for the discreteness of UHECR sources in space and time, a Monte Carlo method is used to randomly place sources in the Galaxy. The low anisotropy measured by Auger is not generally characteristic of the theoretical models, given that the sources are distributed in proportion to the star formation rate, but it can possibly be understood as a) intermittency effects due to the discrete nature of the sources or, with extreme parameters, b) a cancellation of drift current along a current sheet with the outward radial diffusive flux. We conclude that it is possible to interpret the Galactic sub-ankle CR flux as being due entirely to intermittent discrete Galactic sources distributed in proportion to star formation, but only with probability of roughly 35 percent, of which the spectrum is in accord with observations about 30 percent of the time. An alternative explanation for the low anisotropy may be that they are mostly extragalactic and/or heavy.

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

The hypothesis is considered that the present Galactic cosmic ray spectrum is at present softer than its time average due to source intermittency. Measurements of muogenic nuclides underground could provide an independent measurement of the time averaged spectrum. Source intermittency could also account for the surprising low anisotropy reported by the IceCube collaboration. Predictions for Galactic emission of ultrahigh-energy quanta, such as UHE gamma rays and neutrinos, might be higher or lower than previously estimated.

The $^{14}$C production of shock-accelerated particles is calculated in terms of the total energy released in energetic particles. The recently reported 1.2% jump in the $^{14}$C content of the atmosphere in the year C.E. 775, it is found, would require $\gtrsim 10^{34}$ erg in energetic particles, less than first estimates but far more than any known solar flare on record. It is noted that the superflare from a large comet (comparable to C/Hale-Bopp) colliding with the sun could produce shock-accelerated GeV cosmic rays in the solar corona and/or solar wind, and possibly account for the CE 775 event. Several additional predictions of cometary encounters with the sun and other stars may be observable in the future.

We investigate the temporal and spectral correlations between flux and anisotropy fluctuations of TeV-band cosmic rays in the light of recent data taken with IceCube. We find that for a conventional distribution of cosmic-ray sources the dipole anisotropy is higher than observed, even if source discreteness is taken into account. Moreover, even for a shallow distribution of galactic cosmic-ray sources and a reacceleration model, fluctuations arising from source discreteness provide a probability only of the order of 10% that the cosmic-ray anisotropy limits of the recent IceCube analysis are met. This probability estimate is nearly independent of the exact choice of source rate, but generous for a large halo size. The location of the intensity maximum far from the Galactic Center is naturally reproduced.

We show that an alpha effect is driven by the cosmic ray Bell instability exciting left-right asymmetric turbulence. Alfven waves of a preferred polarization have maximally helical motion, because the transverse motion of each mode is parallel to its curl. We show how large-scale Alfven modes, when rendered unstable by cosmic ray streaming, can create new net flux over any finite region, in the direction of the original large-scale field. We perform direct numerical simulations (DNS) of an MHD fluid with a forced cosmic ray current and use the test-field method to determine the alpha effect and the turbulent magnetic diffusivity. As follows from DNS, the dynamics of the instability has the following stages: (i) in the early stage, the small-scale Bell instability that results in a production of small-scale turbulence is excited; (ii) in the intermediate stage, there is formation of larger-scale magnetic structures; (iii) finally, quasi-stationary large-scale turbulence is formed at a growth rate that is comparable to that expected from the dynamo instability, but its amplitude over much longer timescales remains unclear. The results of DNS are in good agreement with the theoretical estimates. It is suggested that this dynamo is what gives weakly magnetized relativistic shocks such as those from gamma ray bursts a macroscopic correlation length. It may also be important for large-scale magnetic field amplification associated with cosmic ray production and diffusive shock acceleration in supernova remnants (SNR) and blast waves from gamma ray bursts. Magnetic field amplification by Bell turbulence in SNR is found to be significant, but it is limited owing to the finite time available to the super-Alfvenicly expanding remnant. The effectiveness of the mechanisms is shown to be dependent on the shock velocity.

We study the propagation of ultra-high-energy cosmic rays (UHECR) in the Galaxy, concentrating on the energy range below the ankle in the spectrum at 4 EeV. A Monte-Carlo method, based on analytical solutions to the time-dependent diffusion problem, is used to account for intermittency by placing sources at random locations. Assuming a source population that scales with baryon mass density or star formation (e.g. long GRB), we derive constraints arising from intermittency and the observational limits on the composition and anisotropy. It is shown that the composition and anisotropy at 1e18 eV are difficult to reproduce and require that either a) the particle mean free path is much smaller than a gyroradius, implying the escape time is very long, b) the composition is heavier than suggested by recent Auger data, c) the ultrahigh-energy sub-ankle component is mostly extragalactic, or d) we are living in a rare lull in the ultrahigh-energy cosmic-ray production, and the current ultrahigh-energy cosmic-ray intensity is far below the Galactic time average. We therefore recommend a strong observational focus on determining the UHECR composition around 1e18 eV.

The maximum cosmic ray energy achievable by acceleration by a relativistic blast wave is derived. It is shown that forward shocks from long GRB in the interstellar medium are powerful enough to produce the Galactic cosmic-ray component up to the ankle at $4\times 10^{18}$eV, as per an earlier suggestion (Levinson and Eichler, 1993). It is further argued that, were extragalactic long GRB responsible for the component {\it above} the ankle as well, the contribution from an occasional Galactic GRB within the solar circle would yield more than the observational limits on the outward flux from the solar circle, unless intermittency and/or beaming causes the present-day contribution to be less than $10^{-3}$ times the time average, and difficulties with these avoidance scenarios are also noted.

The instabilities of relativistic ion beams in a relativistically hot electron background are derived for general propagation angles. It is shown that the Weibel instability in the direction perpendicular to the streaming direction is the fastest growing mode, and probably the first to appear, consistent with the aligned filaments that are seen in PIC simulations. Oblique, quasiperpendicular modes grow almost as fast, as the growth rate varies only moderately with angle, and they may distort or corrugate the filaments after the perpendicular mode saturates.

It is shown that the contribution of faint gamma-ray bursts (GRB) to the total allsky GRB energy flux is small. The allsky flux of GRB appears to be $5.3 \times 10^{-3}$ erg/$\rm cm^2 y$, with little additional component hidden within weak or otherwise undetectable GRB. This significantly constrains physical models of GRB and dirty fireballs, suggesting a rather sharp dichotomy between them should the latter exist.

The prolific magnetar SGR 1900+14 showed two outbursts in the last decade and has been closely monitored in the X-rays to track the changes in its radiative properties. We use archival Chandra and XMM-Newton observations of SGR 1900+14 to construct a history of its spectrum and persistent X-ray flux spanning a period of about seven years. We show that the decline of its X-ray flux in these two outburst episodes follows the same trend. The flux begins to decline promptly and rapidly subsequent to the flares, then decreases gradually for about 600 days, at which point it resumes a more rapid decline. Utilizing the high quality spectral data in each epoch, we also study the spectral coevolution of the source with its flux. We find that neither the magnetic field strength nor the magnetospheric properties change over the period spanned by the observations, while the surface temperature as well as the inferred emitting area both decline with time following both outbursts. We also show that the source reached the same minimum flux level in its decline from these two subsequent outbursts, suggesting that this flux level may be its steady quiescent flux.

Light curves are calculated for an off-axis observer due the scattering of primary radiation by extended baryonic material. The unusually long duration and the chromaticity of the light curves above several KeV of XRF 060218 can be explained as a result of the acceleration of the baryonic scattering material by the primary radiation. The observed light curves by our model and detailed fits to the data are presented. The model predicts that $\sim 4\times 10^{48}$ ergs are put into accelerated, mildly relativistic baryons by the radiation pressure at large radii from the central engine. It is suggested that the emission below 3 KeV, which lies {\it below} the Amati relation, is a baryon contaminated fireball.

The statistical distribution of energies among particles responsible for long Gamma Ray Burst (GRB) emission is analyzed in light of recent results of the Fermi Observatory. The allsky flux, $F_{\gamma}$, recorded by the Gamma Ray Burst Monitor (GBM) is shown, despite its larger energy range, to be not significantly larger than that reported by the Burst and Transient Explorer (BATSE), suggesting a relatively small flux in the 3 - 30 MeV energy range. The present-day energy input rate in $\gamma$-rays recorded by the GBM from long GRB is found, assuming star-formation rates in the literature, to be $\dot W(0)=0.5 F_{\gamma} H/c = 5 \times 10^{42}\ \rm{erg/Mpc^3 yr}$. The Large Area Telescope (LAT) fluence, when observed, is about 5-10\% per decade of the total, in good agreement with the predictions of saturated, non-linear shock acceleration. The high-energy component of long GRBs, as measured by Fermi, is found to contain only $\sim 10^{-2.5}$ of the energy needed to produce ultrahigh-energy cosmic rays (UHECR) above 4 Eev, assuming the latter to be extragalactic, when various numerical factors are carefully included, if the cosmic ray source spectrum has a spectral index of -2. The observed $\gamma$-ray fraction of the required UHECR energy is even smaller if the source spectrum is softer than $E^{-2}$. The AMANDA II limits rule out such a GRB origin for UHECR if much more than $10^{-2}$ of the cosmic ray energy goes into neutrinos that are within, and simultaneous with, the $\gamma$-ray beam. It is suggested that "orphan" neutrinos out of the $\gamma$-ray beam might be identifiable via orphan afterglow { or other wide angle signatures of GRB in lieu of coincidence with prompt $\gamma$-rays}, and it is recommended that feasible single neutrino trigger criteria be established to search for such coincidences.

Fermi measurements of the high-latitude gamma-ray background strongly constrain a decaying-dark-matter origin for the 1--100 GeV Galactic positron anomaly measured with PAMELA. Inverse-Compton scattering of the microwave background by the emergent positrons produces a bump in the diffuse 100-200 MeV gamma-ray background that would protrude from the observed background at these energies. The positrons are thus constrained to emerge from the decay process at a typical energy between ~100 GeV and ~250 GeV. By considering only gamma-ray emission of the excess positrons and electrons, we derive a minimum diffuse gamma-ray flux that, apart from the positron spectrum assumed, is independent of the actual decay modes. Any gamma-rays produced directly by the dark-matter decay leads to an additional signal that make the observational limits more severe. A similar constraint on the energy of emergent positrons from annihilation in dark-matter substructures is argued to exist, according to recent estimates of enhancement in low-mass dark-matter substructures, and improved simulations of such substructure will further sharpen this constraint.

Recently reported observations of magnetar glitches and coincident X-ray long term brightening events establish enough of a data base to indicate that the brightening events are accompanied by glitches at their outset, and that they are probably triggered by the same event that triggers the glitch. We suggest, on the basis of various observational clues, that a) these are caused by energy releases at depths below 100 meters, and b) the unpinning is due to global mechanical motion triggered by the energy release, not by heat. Because mechanical triggering of a glitch requires less energy than a detectable long term X-ray brightening, the latter does not necessarily accompany every glitch, but it is predicted, when it does occur, to be heralded by a glitch. Crustal oscillation associated with the mechanical energy release may cause short ($\ll 10^3$s) flares.

Gravitational waves (GW) can be emitted from coalescing neutron star (NS) and black hole-neutron star (BH-NS) binaries, which are thought to be the sources of short hard gamma ray bursts (SHBs). The gamma ray fireballs seem to be beamed into a small solid angle and therefore only a fraction of detectable GW events is expected to be observationally coincident with SHBs. Similarly ultrahigh energy (UHE) neutrino signals associated with gamma ray bursts (GRBs) could fail to be corroborated by prompt gamma-ray emission if the latter is beamed in a narrower cone than the neutrinos. Alternative ways to corroborate non-electromagnetic signals from coalescing neutron stars are therefore all the more desirable. It is noted here that the extended X-ray tails (XRT) of SHBs are similar to X-ray flashes (XRFs), and that both can be attributed to an off-axis line of sight and thus span a larger solid angle than the hard emission. It is proposed that a higher fraction of detectable GW events may be coincident with XRF/XRT than with hard gamma-rays, thereby enhancing the possibility to detect it as a GW or neutrino source. Scattered gamma-rays, which may subtend a much larger solid angle that the primary gamma ray jet, are also candidates for corroborating non-electromagnetic signals.

The effect of magnetic fields on the frequencies of toroidal oscillations of neutron stars is derived to lowest order. Interpreting the fine structure in the QPO power spectrum of magnetars following giant flares reported by Strohmayer and Watts (2006) to be "Zeeman splitting" of degenerate toroidal modes, we estimate a crustal magnetic field of order 10^{15} Gauss or more. We suggest that residual m, -m symmetry following such splitting might allow beating of individual frequency components that is slow enough to be observed.

It is noted that X-ray tails (XRTs) of short, hard $\gamma$-ray bursts (SHBs) are similar to X-ray flashes (XRFs). We suggest a universal central engine hypothesis, as a way of accounting for this curiosity, in which SHBs differ from long $\gamma$-ray bursts (GRBs) in prompt emission because of the differences in the host star and attendant differences in the environment they present to the compact central engine (as opposed to differences in the central engine itself). Observational constraints and implications are discussed, especially for confirming putative detections of gravitational waves from merging compact objects.

A quantitative theory of spectral lags for $\gamma$-ray bursts (GRBs) is given. The underlying hypothesis is that GRB subpulses are photons that are scattered into our line of sight by local concentrations of baryons that are accelerated by radiation pressure. For primary spectra that are power laws with exponential cutoffs, the width of the pulse and its fast rise, slow decay asymmetry is found to increase with decreasing photon energy, and the width near the exponential cutoff scales approximately as $E_{ph}^{-\eta}$, with $\eta\sim 0.4$, as observed. The spectral lag time is naturally inversely proportional to luminosity, all else being equal, also as observed.

Electrons in shocks are efficiently energized due to the cross-shock potential, which develops because of differential deflection of electrons and ions by the magnetic field in the shock front. The electron energization is necessarily accompanied by scattering and thermalization. The mechanism is efficient in both magnetized and non-magnetized relativistic electron-ion shocks. It is proposed that the synchrotron emission from the heated electrons in a layer of strongly enhanced magnetic field is responsible for gamma ray burst afterglows.

Gamma ray bursts (GRB's) often feature subpulses that have a distinctively asymmetric profile -- they rise quickly and decay much more slowly, while their spectrum softens slightly with observer time. It is suggested that these subpulses are caused by slow baryonic clouds embedded within a primary $\gamma$-ray beam, which scatter the $\gamma$-radiation into our line of sight as they accelerate. Good quantitative agreement is obtained with observed light curves and spectral evolution. The kinetic energy that the baryonic component of GRB jets receives from the primary $\gamma$-radiation is predicted to be about equal to the amount of $\gamma$-radiation that is scattered, consistent with observations of afterglow. Several other observational consequences are briefly discussed. The possibility is raised that the time scale of short GRB is established by radiative acceleration and/or baryon injection rather than the time scale of the central engine.

Slow dissipation of non-potential magnetic fields in the magnetosphere of the magnetar is assumed to accelerate particles to hundreds MeV along the magnetic field lines. We consider interaction of fast particles with the surface of the magnetar. We argue that the collisionless dissipation does not work in the atmosphere of the neutron star because the two-stream instability is stabilized by the inhomogeneity of the atmosphere. Rather, the dominant dissipation mechanism is collisional Landau levelexcitations followed by pair production via the deexcitation gamma-rays ultimately leading to electrons with the energy below the Landau energy. We show that, because of the effects of the superstrong magnetic field, these electrons could emerge from the surface carrying most of the original energy so that a hot corona arises with the temperature of 1 - 2 MeV. This extended corona is better suited than a thin atmosphere to convert most of the primary beam energy to non-thermal radiation and, as we show, most of the coronal energy release is radiated away in the hard X-ray and the soft gamma-ray bands by Comptonization and bremsstrahlung. The radiation spectrum is a power-law with the photon index $1<\alpha<2$. The model may account for the persistent hard X-ray emission discovered recently from the soft gamma-ray repeaters and anomalous X-ray pulsars and predicts that the radiation spectrum is extended into the MeV band.

A model is presented for the quasiperiodic component of magnetar emission during the tail phase of giant flares. The model invokes modulation of the particle number density in the magnetosphere. The magnetospheric currents are modulated by torsional motion of the surface and we calculate that the amplitude of neutron star surface oscillation should be ~1% of the NS radius in order to produce the observed features in the power spectrum. Using an axisymmetric analytical model for structure of the magnetosphere of an oscillating NS, we calculate the angular distribution of the optical depth to the resonant Compton scattering. The anisotropy of the optical depth may be why QPO are observed only at particular rotational phases.

This paper has been withdrawn by the authors due to missing references to earlier work and an error in the interpretation of Figures 2 and 3. A corrected version is in preparation.

X-ray emission following giant flares of magnetars can be categorized into three categories of time scales (a) short term afterglow (b) medium term afterglow and (c) long term afterglow. Short term afterglow, which declines over several hours, seems to correspond to gravitational resettling of uplifted material. Medium term afterglow, which declines over several weeks or months, appears to be the cooling of the heated outer crust, and long term afterglow, which declines over a period of many years, can be understood to be the cooling of the inner crust. The long term afterglow profile may be a very sensitive indicator of neutron star mass.

We present detailed spectral and timing analysis of the hard x-ray transient IGR J16358-4726 using multi-satellite archival observations. A study of the source flux time history over 6 years, suggests that lower luminosity transient outbursts can be occuring in intervals of at most 1 year. Joint spectral fits of the higher luminosity outburst using simultaneous Chandra/ACIS and INTEGRAL/ISGRI data reveal a spectrum well described by an absorbed power law model with a high energy cut-off plus an Fe line. We detected the 1.6 hour pulsations initially reported using Chandra/ACIS also in the INTEGRAL/ISGRI light curve and in subsequent XMM-Newton observations. Using the INTEGRAL data we identified a spin up of 94 s (dP/dt = 1.6E-4), which strongly points to a neutron star nature for IGR J16358-4726. Assuming that the spin up is due to disc accretion, we estimate that the source magnetic field ranges between 10^13 - 10^15 G, depending on its distance, possibly supporting a magnetar nature for IGR J16358-4726.

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

Early X-ray afterglows recently detected by {\it Swift} frequently show a phase of very shallow flux decay lasting from $\sim 10^{2.5} $s up to $\sim 10^4 $s, followed by a steeper, more familiar decay. We suggest that the flat early part of the light curve may be a combination of the decaying tail of the prompt emission and the delayed onset of the afterglow emission observed from viewing angles slightly outside the edge of the region within the jet with prominent afterglow emission, as predicted previously. This would imply that a significant fraction of viewers get very little external shock energy along their line of sight, and, therefore, see a very high $\gamma$-ray to kinetic energy ratio at early times. The early flat phase in the afterglow light curve implies, in a rather robust and model independent manner, a very large $\gamma$-ray efficiency, typically $\gtrsim 90%$, which is very difficult to extract from baryons by internal shocks.

It is estimated that the Weibel instability is not generally an effective mechanism for generating ultrarelativistic astrophysical shocks. Even if the upstream magnetic field is as low as in the interstellar medium, the shock is mediated not by the Weibel instability but by the Larmor rotation of protons in the background magnetic field. Future simulations should be able to verify or falsify our conclusion.

The 2004 Dec. 27 giant Gamma-ray flare detected from the magnetar SGR 1806-20 created an expanding radio nebula which we have monitored with the Australia Telescope Compact Array and the Very Large Array. These data indicate that there was an increase in the observed flux ~25 days after the initial flare that lasted for ~8 days, which we believe is the result of ambient material swept-up and shocked by this radio nebula. For a distance to SGR 1806-20 of 15 kpc, using the properties of this rebrightening we infer that the initial blast wave was dominated by baryonic material of mass M>10^{24.5} g. For an initial expansion velocity v~0.7c (as derived in an accompanying paper), we infer this material had an initial kinetic energy E>10^{44.5} ergs. If this material originated from the magnetar itself, it may have emitted a burst of ultra-high energy (E > 1 TeV) neutrinos far brighter than that expected from other astrophysical sources.

On 2004 Dec. 27, the soft gamma repeater (SGR) 1806-20 emitted the brightest giant flare (GF) ever detected from an SGR, with an (isotropic) energy release $\sim 100$ times greater than the only two other known SGR GFs. It was followed by a very bright, fading radio afterglow. Extensive follow-up radio observations provided a wealth of information with unprecedented astrometric precision, revealing the temporal evolution of the source size, along with densely sampled light curves and spectra. Here we expand on our previous work on this source, by explaining these observations within one self-consistent dynamical model. In this scenario, the early radio emission is due to the outflow ejected during the GF energizing a thin shell surrounding a pre-existing cavity, where the observed steep temporal decay of the radio emission seen beginning on day 9 is attributed to the adiabatic cooling of the shocked shell. The shocked ejecta and external shell move outward together, driving a forward shock into the ambient medium, and are eventually decelerated by a reverse shock. As we show in Gelfand et al. (2005), the radio emission from the shocked external medium naturally peaks when significant deceleration occurs, and then decays relatively slowly. The dynamical modeling of the collision between the ejecta and the external shell together with the observed evolution of the source size (which is nicely reproduced in our model) suggest that most of the energy in the outflow was in mildly relativistic material, with an initial expansion velocity $v/c \lesssim 0.7d_{15}$, for a distance of $15d_{15}$ kpc to SGR 1806-20. An initially highly relativistic outflow would not have produced a long coasting phase at a mildly relativistic expansion velocity, as was observed.

The extraordinary giant flare (GF) of 2004 December 27 from the soft gamma repeater (SGR) 1806-20 was followed by a bright radio afterglow. We present an analysis of VLA observations of this radio afterglow from SGR 1806-20, consisting of previously reported 8.5 GHz data covering days 7 to 20 after the GF, plus new observations at 8.5 and 22 GHz from day 24 to 81. We detect motion in the flux centroid of the afterglow, at an average velocity of 0.26 +/- 0.03 c (assuming a distance of 15 kpc) at a position angle of -45 degrees. This motion, in combination with the growth and polarization measurements, suggests an asymmetric outflow, mainly from one side of the magnetar. We find a deceleration in the expansion, from ~9 mas/day to <5 mas/day. The time of deceleration is roughly coincident with the rebrightening in the radio light curve, as expected to result when the ejecta from the GF sweeps up enough of the external medium, and transitions from a coasting phase to the Sedov-Taylor regime. The radio afterglow is elongated and maintains a 2:1 axis ratio with an average position angle of -40 degrees (north through east), oriented perpendicular to the average intrinsic linear polarization angle.

An astronomical complex intended to detect optical transients (OTs) in a wide field and follow them up with high time resolution investigation is described.

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

Soft gamma repeaters (SGRs) are "magnetars", a small class of slowly spinning neutron stars with extreme surface magnetic fields, ~10^15 gauss. On 2004 December 27, a giant flare was detected from the magnetar SGR 1806-20, the third such event ever recorded. This burst of energy was detected by a variety of instruments and even caused an ionospheric disturbance in the Earth's upper atmosphere recorded around the globe. Here we report the detection of a fading radio afterglow produced by this outburst, with a luminosity 500 times larger than the only other detection of a similar source. From day 6 to day 19 after the flare from SGR 1806-20, a resolved, linearly polarized, radio nebula was seen, expanding at approximately a quarter the speed of light. To create this nebula, at least 4x10^43 ergs of energy must have been emitted by the giant flare in the form of magnetic fields and relativistic particles. The combination of spatially resolved structure and rapid time evolution allows a study in unprecedented detail of a nearby analog to supernovae and gamma-ray bursts.

The recent superflare of 27 December 2004 from the magnetar SGR 1806-20 was the brightest extrasolar flash ever recorded in the modern era. The chances for seeing exotic ultrahigh energy (UHE) radiation - neutrons, neutrinos, gamma rays and charged cosmic rays - from it are far better from an energetic point of view than from cosmological gamma ray bursts (GRBs). The chances for detecting the various components are discussed in light of recent data from the 27 December event.

We have suggested previously that baryons in GRB fireballs infiltrate from the surrounding walls that collimate the fireball. The efficiency $\epsilon_b$ for generating blast energy can then be angle dependent. Delayed onset of afterglow can be interpreted as being due to a baryon-poor viewing angle.

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

I review several ongoing connections between space physics and astrophysics: a) Measurements of energetic particle spectra have confirmed theoretical prediction of the highest energy to which shocks can accelerate particles, and this has direct bearing on the origin of the highest energy cosmic rays. b) Mass ejection in solar flares may help us understand photon ejection in the giant flares of magnetar outbursts. c) Measurements of electron heat fluxes in the solar wind can help us understand whether heat flux in tenuous astrophysical plasma is in accordance with the classical Spitzer-Harm formula or whether it is reduced well below this value by plasma instabilities.

Magnetars comprise two classes of rotating neutron stars (Soft Gamma Repeaters (SGRs) and Anomalous X-ray Pulsars), whose X-ray emission is powered by an ultrastrong magnetic field, B ~ 10^15 G. Occasionally SGRs enter into active episodes producing many short X-ray bursts; extremely rarely (about once per 50 years per source), SGRs emit a giant flare, an event with total energy at least 1000 times higher than their typical bursts. Here we report that, on 2004 December 27, SGR 1806-20 emitted the brightest extra-solar transient event ever recorded, even surpassing the full moon brightness for 0.2 seconds. The total (isotropic) flare energy is 2x10^46 erg, 100 times higher than the only two previous events, making this flare a once in a century event. This colossal energy release likely occurred during a catastrophic reconfiguration of the magnetar's magnetic field. Such an event would have resembled a short, hard Gamma Ray Burst (GRB) if it had occurred within 40 Mpc, suggesting that extragalactic SGR flares may indeed form a subclass of GRBs.

The anomalous bump in the cosmic ray positron to electron ratio at $10 GeV$ can be explained as being a component from a point source that was originally harder than the primary electron background and degrades due to synchrotron and inverse Compton losses in the Galaxy while propagating to the Earth's vicinity. The fit is better than can be obtained with homogeneous injection and is attributed to a minimum age threshold. Annihilating neutralinos can provide a fair fit to the data if they have a mass just above 1/2 the mass of the $Z^o$ and if they annihilate primarily in distant density concentrations in the Galaxy. A possible observational consequence of this scenario would be intense inverse Comptonization of starlight at the Galactic center, with a sharp energy cutoff in the emergent photons as a possible signature of the neutralino mass.

It is noted that a surface layer of matter in contact with a sufficiently super-Eddington, radially combed photon flux typically attains a relativistic coasting state whereby the radiation does not accelerate the matter. Radiation that scatters off this layer is most likely to be observed along the velocity vector of the matter, where it would be most strongly polarized.

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

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

It is suggested that a Poynting flux-dominated outflow with a sufficiently strong magnetic field can pick up hard X-ray photons when the magnetic field is of sufficient strength. The zeroth generation X-rays pair produce, and the pairs radiate extremely energetic first generation photons that could be detected by extensive air shower arrays and/or MILAGRO if they escape the production site intact. Giant flares from magnetars may thus yield bursts of UHE photons. GRB-associated Poynting flows may be unstable to pair production near their source, and their energy rapidly converted to pairs.

SGR 1627-41 was discovered in 1998 after a single active episode which lasted \~6 weeks. We report here our monitoring results of the decay trend of the persistent X-ray luminosity of the source during the last 5 years. We find an initial temporal power law decay with index 0.47, reaching a plateau which is followed by a sharp (factor of ten) flux decline ~800 days after the source activation. The source spectrum is best described during the entire period by a single power law with high absorption (N_H=9.0(7)x10^(22) cm^(-2)); the spectral index, however, varies dramatically between 2.2-3.8 spanning the entire range for all known SGR sources. We discuss the cooling behavior of the neutron star assuming a deep crustal heating initiated by the burst activity of the source during 1998.

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

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

It is shown that bulk acceleration during reconnection of extremely strong magnetic fields near compact objects can accelerate ions to Lorentz factors of $\sim 10^2 \sigma^{3/5}$ under general conditions, where $\sigma$, the magnetic energy per current-carrying proton rest energy, can approach $10^{15}$. For magnetar-type fields, neutrons and neutrinos can be generated at potentially detectable levels via hadron polarization. Ultrahigh energy photons can also be emitted and escorted from the high field region by Poynting flux.

The X-ray afterglow that is observed following large flares on magnetars can be accurately fit by simple and quantitative theoretical models: The long term afterglow, lasting of order weeks, can be understood as thermal radiation of a heated neutron star crust that is re-scattered in the magnetosphere. Short term afterglow is well fit by the cooling of a non-degenerate, pair-rich layer, which gradually shrinks and releases heat to a pair-free zone above it. Measurements of persistent optical and infrared emission directly probe long lived currents in the magnetosphere which are a likely source of collective plasma emission. The superstrong magnetic field plays an important role in generating these various emissions, and previous inference of its strength in magnetars is supported by the good fits with observation.

Transient X-ray emission, with an approximate t^{-0.7} decay, was observed from SGR 1900+14 over 40 days following the the giant flare of 27 Aug 1998. We calculate in detail the diffusion of heat to the surface of a neutron star through an intense 10^{14}-10^{15} G magnetic field, following the release of magnetic energy in its outer layers. We show that the power law index, the fraction of burst energy in the afterglow, and the return to persistent emission can all be understood if the star is composed of normal baryonic material.

It is proposed that magnetospheric currents above the surfaces of magnetars radiate coherent emission in analogy to pulsars. Scaling the magnetospheric parameters suggests that the coherent emission from magnetars would emerge in the infra-red or optical.

It is suggested that magnetars eventually become unstable to a dynamic overturning instability that destroys most of their dipole moment in a single event. It is further suggested that such an instability would produce a supergiant version of the largest soft gamma ray bursts that have thus far been observed, and that they could be observed out to tens or even $10^2$ Mpc. Persistent emission and afterglows from the fireball might confirm positive identification with nearby galaxies and perhaps provide other distinguishing features as well.

The recent Chandra image of the Crab nebula shows a striking, axisymmetric polar jet. It is shown that jets are formed in axisymmetric, magnetized pulsar winds and that the jet luminosity scales relative to the total as (\gamma_0\sigma_{eq})^{-4/3}, where \sigma_{eq} is the ratio of Poynting flux to particle kinetic energy output at the equator at the base of the flow and \gamma_0 the initial Lorentz factor of the flow. The results are applied to the image of the Crab nebula, and the limit is set for the Crab pulsar of \sigma_{eq} \leq 100. It is argued that conventional pulsar theory needs to be reexamined in light of these limits.

Quantum gravity "foam", among its various generic Lorentz non-invariant effects, would cause neutrino mixing. It is shown here that, if the foam is manifested as a nonrenormalizable effect at scale M, the oscillation length generically decreases with energy $E$ as (E/M)^(-2). Neutrino observatories and long-baseline experiments should have therefore already observed foam-induced oscillations, even if M is as high as the Planck energy scale. The null results, which can be further strengthened by better analysis of current data and future experiments, can be taken as experimental evidence that Lorentz invariance is fully preserved at the Planck scale, as is the case in critical string theory.

We propose using large Air Cerenkov Telescopes (ACT's) to search for optical, pulsed signals from extra-terrestrial intelligence. Such dishes collect tens of photons from a nanosecond-scale pulse of isotropic equivalent power of tens of solar luminosities at a distance of 100 pc. The field of view for giant ACT's can be on the order of ten square degrees, and they will be able to monitor 10 to 10$^2$ stars simultaneously for nanosecond pulses of about 6th mag or brighter. Using the Earth's diameter as a baseline, orbital motion of the planet could be detected by timing the pulse arrival times.

An improved linear stability theory of small-amplitude oscillations of a self-gravitating, infinitesimally thin gaseous disk of spiral galaxies has been developed. It was shown that in the differentially rotating disks for nonaxisymmetric perturbations Toomre's modified hydrodynamical simulations to test the validity of the modified local criterion.

The time evolution of models for an isolated disk of highly flattened galaxies of stars is investigated by direct integration of the Newtonian equations of motion of N=30,000 identical stars over a time span of many galactic rotations. Certain astronomical implications of the simulations to actual disk-shaped (i.e. rapidly rotating) galaxies are explored as well.

The kinetic theory is used to study the evolution of the self-gravitating disk of planetesimals. The effects of frequent collisions between planetesimals are taken into account by using a Krook integral in the Boltzmann kinetic equation. It is shown that as a result of an aperiodic collision-dissipative instability of small gravity disturbances the disk is subdivided into numerous dense fragments. These can eventually condense into the planetary sequence.

Linear theory is used to determine the stability of the self-gravitating, rapidly (and nonuniformly) rotating, two-dimensional, and collisional particulate disk against small-amplitude gravity perturbations. A gas-kinetic theory approach is used by exploring the combined system of the Boltzmann and the Poisson equations. The effects of physical collisions between particles are taken into account by using in the Boltzmann kinetic equation a Krook model integral of collisions modified to allow collisions to be inelastic. It is shown that as a direct result of the classical Jeans instability and a secular dissipative-type instability of small-amplitude gravity disturbances (e.g. those produced by a spontaneous perturbation and/or a companion system) the disk is subdivided into numerous irregular ringlets, with size and spacing of the order of the Jeans length (= (4-6) mean epicyclic radius). The present research is aimed above all at explaining the origin of various structures in highly flattened, rapidly rotating systems of mutually gravitating particles. In particular, it is suggested that forthcoming Cassini spacecraft high-resolution images may reveal this kind of hyperfine of the order of 100 m structure in the main rings A, B, and C of the Saturnian ring system.

The reaction of collective oscillations excited in the interaction between aperiodically growing Jeans-type gravity perturbations and stars of a rapidly rotating disk of flat galaxies is considered. An equation is derived which describes the change in the main body of equilibrium distribution of stars in the framework of the nonresonant weakly nonlinear theory. Certain applications of the theory to the problem of relaxation of the Milky Way at radii where two-body relaxation is not effective are explored. The theory, as applied to the Solar neighborhood, accounts for observed features, such as the shape for the velocity ellipsoid of stars and the increase in star random velocities with age.

The time evolution of initially balanced, rapidly rotating models for an isolated disk of highly flattened galaxies of stars is calculated. The method of direct integration of the Newtonian equations of motion of stars over a time span of many galactic rotations is applied. Use of modern concurrent supercomputers has enabled us to make long simulation runs using a sufficiently large number of particles N=30,000. One of the goals of the present simulation is to test the validities of a modified local criterion for stability of Jeans-type gravity perturbations (e.g. those produced by a barlike structure, a spontaneous perturbation and/or a companion galaxy) in a self-gravitating, infinitesimally thin and collisionless disk. In addition to the local criterion we are interested in how model particles diffuse in velocity. This is of considerable interest in the kinetic theory of stellar disks. Certain astronomical implications of the simulations to actual disk-shaped (i.e. rapidly rotating) galaxies are explored. The weakly nonlinear, or quasi-linear kinetic theory of the Jeans instability in disk galaxies of stars is described as well.

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

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

The recent detection of high energy gamma-rays coming from supernova remnants and active galactic nuclei has revived interest in the diffusive shock acceleration of electrons. In the present paper we examine the basis of the so-called ``box'' model for particle acceleration and present a more physical version of it. Using this we determine simple criteria for the conditions under which ``pile-ups'' can occur in shock accelerated electron spectra subject to synchrotron or inverse Compton losses (the latter in the Thompson limit). An extension to include nonlinear effects is proposed.

The hypothesis is considered that the vacuum is a Lorentz non-invariant foam in which translational symmetry is spontaneously broken at the Planck or grand unification scale. This could possibly be observed via Rayleigh scattering of ultrahigh energy quanta, and it appears to be ruled out for spin 1/2 leptons and independently testable for photons. A weaker version of the hypothesis predicts (for otherwise massless neutrinos) irreversible neutrino mixing over a length scale of $l_{Planck}(E_{\nu}/E_{Planck})^{-2}$ in comoving cosmic coordinates for a purely four dimensional universe, and, for neutrinos that have spin 3/2 in a higher dimensional manifold above the energy scale $\eta E_{Planck}$, the mixing is predicted to take place over a scale of $\eta^{-4}$ times this length. Such neutrino mixing might also be observable with atmospheric and solar neutrino experiments.

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

It is hypotesized that neutrinos are massles and mixed among three (muon, electron and tau) flavors with a mixing length proportional to (E_nu/E_Planck)^-2. It is suggested that this hypothesis mostly reconciles the standard solar models with the observations of the solar Boron and pp neutrinos, and the atmospheric muon-to-electron neutrinos ratio.

We calculate a grid of spherically symmetric OB stellar atmospheres at low metallicities, including both non-local thermodynamic equilibrium (NLTE) and metal line blanketing effects. This is done to assess the uncertainties in helium abundance determination by nebular codes due to input stellar atmosphere models. The more sophisticated stellar atmosphere models we use can differ from LTE models by as much as 40 percent in the ratio of He to H-ionizing photons.

Heat-transfer through weakly magnetized diffuse astrophysical plasmas excites whistlers. This leads to electron whistler resonant scattering, a reduction of the electron mean-free path, and heat-flux inhibition. However, only whistlers propagating at a finite angle to the magnetic field (off-axis) can scatter the heat-flux carrying electrons. Thus, the level of heat flux-inhibition along the magnetic field lines depends on the presence of off-axis whistlers.

Companion stars exposed to high energy radiation from a primary neutron star or accreting black hole can experience significant spallation of their heavy elements, so that their atmospheres would be extremely rich in lithium, beryllium, and especially boron. In this paper we note that the detection or non-detection of these elements, and their relative abundances if detected, would provide a diagnostic of the high energy output of the primary, and possibly the shock acceleration of particles at the companion's bow shock in a pulsar wind.

The possibility of detecting proto-galaxies in the UV band is pointed out, assuming that galaxy formation occured at z $\sim 5-6$. It is shown that the diffuse gas in collapsing galaxy sized objects with temperatures $\sim 10^{6\pm0.5}$ K, and with a modest amount of metallicity, should copiously produce emission lines from highly ionized Iron atoms. The expected luminosity from models of galaxy formation is compared with the sensitivity of HST.