Fast Radio Bursts (FRBs) are millisecond-long radio pulses of extragalactic origin with peak luminosities far exceeding any Milky Way sources. The prevalent invocation for the FRB origin involves magnetars: young, magnetically powered neutron stars with the strongest magnetic fields in the Universe. A magnetar-defining signature is the emission of bright, hard X-ray bursts of sub-second duration. These occur in isolation or during a burst storm, when several hundred are observed within minutes to hours. On April 27th 2020, the Galactic magnetar SGR J1935+2154 entered an active period, emitting hundreds of X-ray bursts in a few hours. Remarkably, only one of these temporally coincided with an FRB-like radio burst. Here we report on the spectral and temporal analyses of 24 X-ray bursts emitted 13 hours prior to the FRB and seen simultaneously with NASA NICER and Fermi/GBM missions in their combined energy range. We demonstrate that the FRB-associated X-ray burst is very similar temporally, albeit strikingly different spectrally, from the 24 NICER/GBM bursts. If the FRB-associated burst were drawn from this magnetar burst population, its occurrence rate would be at most around 1 in 7000. This rarity combined with the unusual X-ray burst spectrum is perhaps indicative of an uncommon locale for the origin of the FRB-associated burst. We suggest that this unique event originated in quasi-polar open or closed magnetic field lines extending to high altitudes where radio emission can be generated, possibly from a collimated plasma flow.

We report on the timing and spectral properties of the soft X-ray emission from the magnetar 1E 2259+586 from January 2013, $\sim 8$ months after the detection of an anti-glitch, until September 2019, using the Neil Gehrels Swift and NICER observatories. During this time span, we detect two timing discontinuities. The first, occurring around 5 years after the April 2012 anti-glitch, is a relatively large spin-up glitch with a fractional amplitude $\Delta\nu/\nu=1.24(2)\times10^{-6}$. We find no evidence for flux enhancement or change in the spectral or pulse profile shape around the time of this glitch. This is consistent with the picture that a significant number of magnetar spin-up glitches are radiatively-quiet. Approximately 1.5 years later in April 2019, 1E 2259+586 exhibited an anti-glitch with spin-down of a fractional amplitude $\Delta\nu/\nu=-5.8(1)\times10^{-7}$; similar to the fractional change detected in 2012. We do not, however, detect any change to the pulse-profile shape or increase in the rms pulsed flux of the source, nor do we see any possible bursts from its direction around the time of the anti-glitch; all of which occurred during the 2012 event. Hence, similar to spin-up glitches, anti-glitches can occur silently. This may suggest that these phenomena originate in the neutron star interior, and that their locale and triggering mechanism do not necessarily have to be connected to the magnetosphere. Lastly, our observations suggest that the occurrence rate of spin-up and spin-down glitches is about the same in 1E 2259+586, with the former having a larger net fractional change.

During a pointed 2018 NuSTAR observation, we detected a flare with a 2.2 hour duration from the magnetar 1RXS J170849.0$-$400910. The flare, which rose in $\sim25$ seconds to a maximum flux 6 times larger than the persistent emission, is highly pulsed with an rms pulsed fraction of $53\%$. The pulse profile shape consists of two peaks separated by half a rotational cycle, with a peak flux ratio of $\sim$2. The flare spectrum is thermal with an average temperature of 2.1 keV. Phase resolved spectroscopy show that the two peaks possess the same temperature, but differ in size. These observational results along with simple light curve modeling indicate that two identical antipodal spots, likely the magnetic poles, are heated simultaneously at the onset of the flare and for its full duration. Hence, the origin of the flare has to be connected to the global dipolar structure of the magnetar. This might best be achieved externally, via twists to closed magnetospheric dipolar field lines seeding bombardment of polar footpoint locales with energetic pairs. Approximately 1.86 hours following the onset of the flare, a short burst with its own 3-minute thermal tail occurred. The burst tail is also pulsating at the spin period of the source and phase-aligned with the flare profile, implying an intimate connection between the two phenomena. The burst may have been caused by a magnetic reconnection event in the same twisted dipolar field lines anchored to the surface hot spots, with subsequent return currents supplying extra heat to these polar caps.

We explore the burst energy distribution of fast radio bursts (FRBs) in the low-twist magnetar model of Wadiasingh and Timokhin (2019). Motivated by the power-law fluence distributions of FRB 121102, we propose an elementary model for the FRB luminosity function of individual repeaters with an inversion protocol which directly relates the power-law distribution index of magnetar short burst fluences to that for FRBs. The protocol indicates the FRB energy scales virtually linearly with crust/field dislocation amplitude, if magnetar short bursts prevail in the magnetoelastic regime. Charge starvation in the magnetosphere during bursts (required in WT19) for individual repeaters implies the predicted burst fluence distribution is narrow, $\lesssim 3$ decades for yielding strains and oscillation frequencies feasible in magnetar crusts. Requiring magnetic confinement and charge starvation, we obtain a death line for FRBs which segregates magnetars from the normal pulsar population, suggesting only the former will host recurrent FRBs. We convolve the burst energy distribution for individual magnetars to define the distribution of luminosities in evolved magnetar populations. The broken power-law luminosity function's low energy character depends on the population model, while the high energy index traces that of individual repeaters. Independent of the evolved population, the broken power-law isotropic-equivalent energy/luminosity function peaks at $\sim10^{37}-10^{40}$ erg with a low-energy cutoff at $\sim 10^{37}$ erg. Lastly, we consider the local fluence distribution of FRBs, and find that it can constrain the subset of FRB-producing magnetar progenitors. Our model suggests that improvements in sensitivity may reveal flattening of the global FRB fluence distribution and saturation in FRB rates.

The accretion-powered X-ray pulsar GX 301-2 was observed with the balloon-borne X-Calibur hard X-ray polarimeter during late December 2018, with contiguous observations by the NICER X-ray telescope, the Swift X-ray Telescope and Burst Alert Telescope, and the Fermi Gamma-ray Burst Monitor spanning several months. The observations detected the pulsar in a rare apastron flaring state coinciding with a significant spin-up of the pulsar discovered with the Fermi GBM. The X-Calibur, NICER, and Swift observations reveal a pulse profile strongly dominated by one main peak, and the NICER and Swift data show strong variation of the profile from pulse to pulse. The X-Calibur observations constrain for the first time the linear polarization of the 15-35 keV emission from a highly magnetized accreting neutron star, indicating a polarization degree of (27+38-27)% (90% confidence limit) averaged over all pulse phases. We discuss the spin-up and the X-ray spectral and polarimetric results in the context of theoretical predictions. We conclude with a discussion of the scientific potential of future observations of highly magnetized neutron stars with the more sensitive follow-up mission XL-Calibur.

Jets in blazars are an excellent forum for studying acceleration at relativistic shocks using the highly-variable emission seen across the electromagnetic spectrum. Our recent work on combining multi-wavelength leptonic emission models with simulated thermal+non-thermal distributions from shock acceleration theory has resulted in new insights into plasma conditions in blazars. This has demonstrated the ability to infer the cyclotron frequency, the plasma density and thus also the Alfven speed, thereby determining the rapidity of particle energization. An important inference was that MHD turbulence levels decline with remoteness from jet shocks. This paper outlines new results from our recent extension of this program to a two-zone, time-evolving construction, modeling together both extended, enhanced emission states from larger radiative regions, and prompt flare events from compact acceleration zones. These are applied to flares in the FSRQ blazar 3C 279 monitored by Fermi-LAT in gamma-rays in late 2013. With impulsive injection episodes from the shock zone, as the acceleration first proceeds and then abates, the radiative simulations obtain a pronounced spectral hardening in the optical and gamma-ray bands as the flare grows, followed by a softening during the decay phase. For 3C 279, while model radio and X-ray synchrotron flares are temporally correlated, there is a lag in both bands relative to GeV gamma rays and optical emission on timescales of a number of hours. This delay is governed by the short cooling time associated with the bright external Compton signal.

The All-sky Medium Energy Gamma-ray Observatory (AMEGO) is a probe class mission concept that will provide essential contributions to multimessenger astrophysics in the late 2020s and beyond. AMEGO combines high sensitivity in the 200 keV to 10 GeV energy range with a wide field of view, good spectral resolution, and polarization sensitivity. Therefore, AMEGO is key in the study of multimessenger astrophysical objects that have unique signatures in the gamma-ray regime, such as neutron star mergers, supernovae, and flaring active galactic nuclei. The order-of-magnitude improvement compared to previous MeV missions also enables discoveries of a wide range of phenomena whose energy output peaks in the relatively unexplored medium-energy gamma-ray band.

Mildly-relativistic shocks that are embedded in colliding magnetohydrodynamic flows are prime sites for relativistic particle acceleration and production of strongly variable, polarized multi-wavelength emission from relativistic jet sources such as blazars and gamma-ray bursts. The principal energization mechanisms at these shocks are diffusive shock acceleration and shock drift acceleration. In recent work, we had self-consistently coupled shock acceleration and radiation transfer simulations in blazar jets in a basic one-zone scenario. These one-zone models revealed that the observed spectral energy distributions (SEDs) of blazars strongly constrain the nature of the shock layer hydromagnetic turbulence. In this paper, we expand our previous work by including full time dependence and treating two zones, one being the site of acceleration, and the second being a larger emission zone. This construction is applied to multiwavelength flares of the flat spectrum radio quasar 3C 279, fitting snap-shot SEDs and generating light curves that are consistent with observed variability timescales. We also present a generic study for the typical flaring behavior of the BL Lac object Mrk 501. The model predicts correlated variability across all wavebands, but cross-band time lags depending on the type of blazar (FSRQ vs. BL Lac), as well as distinctive spectral hysteresis patterns in all wavelength bands, from mm radio waves to gamma-rays. These evolutionary signatures serve to provide diagnostics on the competition between acceleration and radiative cooling.

MAXI J1621-501 is the first Swift/XRT Deep Galactic Plane Survey transient that was followed up with a multitude of space missions (NuSTAR, Swift, Chandra, NICER, INTEGRAL, and MAXI) and ground-based observatories (Gemini, IRSF, and ATCA). The source was discovered with MAXI on 2017 October 19 as a new, unidentified transient. Further observations with NuSTAR revealed 2 Type I X-ray bursts, identifying MAXI J1621-501 as a Low Mass X-ray Binary (LMXB) with a neutron star primary. Overall, 24 Type I bursts were detected from the source during a 15 month period. At energies below 10 keV, the source spectrum was best fit with three components: an absorbed blackbody with kT = 2.3 keV, a cutoff power law with index $\Gamma{}$ = 0.7, and an emission line centered on 6.3 keV. Timing analysis of the X-ray persistent emission and burst data has not revealed coherent pulsations from the source or an orbital period. We identified, however, a super-orbital period $\sim{}$78 days in the source X-ray light curve. This period agrees very well with the theoretically predicted radiative precession period of $\sim{}$82 days. Thus, MAXI J1621-501 joins a small group of sources characterized with super-orbital periods.

The X-ray Polarization Probe (XPP) is a second generation X-ray polarimeter following up on the Imaging X-ray Polarimetry Explorer (IXPE). The XPP will offer true broadband polarimetery over the wide 0.2-60 keV bandpass in addition to imaging polarimetry from 2-8 keV. The extended energy bandpass and improvements in sensitivity will enable the simultaneous measurement of the polarization of several emission components. These measurements will give qualitatively new information about how compact objects work, and will probe fundamental physics, i.e. strong-field quantum electrodynamics and strong gravity.

This white paper highlights compact object and fundamental physics science opportunities afforded by high-throughput broadband (0.1-60 keV) X-ray polarization observations. X-ray polarimetry gives new observables with geometric information about stellar remnants which are many orders of magnitude too small for direct imaging. The X-ray polarimetric data also reveal details about the emission mechanisms and the structure of the magnetic fields in and around the most extreme objects in the Universe. Whereas the Imaging X-ray Polarimetry Explorer (IXPE) to be launched in 2021 will obtain first results for bright objects, a follow-up mission could be one order of magnitude more sensitive and would be able to use a broader bandpass to perform physics type experiments for representative samples of sources.

Over the last four decades, persistent and flaring emission of magnetars observed by various telescopes has provided us with a suite of light curves and spectra in soft and hard X-rays, with no emission yet detected above around 1 MeV. Attenuation of such high-energy photons by magnetic pair creation and photon splitting is expected to be active in the magnetospheres of magnetars, possibly accounting for the paucity of gamma-rays in their signals. This paper explores polarization-dependent opacities for these two QED processes in static vacuum dipole magnetospheres of highly-magnetized neutron stars, calculating attenuation lengths and determining escape energies, which are the maximum photon energies for transparency out to infinity. The numerical trajectory integral analysis in flat and curved spacetimes provides upper bounds of a few MeV or less to the visible energies for magnetars for locales proximate to the stellar surface. Photon splitting opacity alone puts constraints on the possible emission locales in their magnetospheres: regions within field loops of maximum altitudes 2-4 stellar radii are not commensurate with maximum detected energies of around 250 keV. These constraints apply not only to magnetar flares but also to their quiescent hard X-ray tail emission. An exploration of photon splitting attenuation in the context of a resonant inverse Compton scattering model for the hard X-ray tails derives distinctive phase-resolved spectroscopic and polarimetric signatures, of significant interest for future MeV-band missions such as AMEGO and e-ASTROGAM.

Mildly relativistic, oblique shocks are frequently invoked as possible sites of relativistic particle acceleration and production of strongly variable, polarized multi-wavelength emission from relativistic jet sources such as blazars, via diffusive shock acceleration (DSA). In recent work, we had self-consistently coupled DSA and radiation transfer simulations in blazar jets. These one-zone models determined that the observed spectral energy distributions (SEDs) of blazars strongly constrain the nature of the hydromagnetic turbulence responsible for pitch-angle scattering. In this paper, we expand our previous work by including full time dependence and treating two emission zones, one being the site of acceleration. This modeling is applied to a multiwavelength flare of the flat spectrum radio quasar 3C~279, fitting snap-shot SEDs and light curves. We predict spectral hysteresis patterns in various energy bands as well as cross-band time lags with optical and GeV gamma-rays as well as radio and X-rays tracing each other closely with zero time lag, but radio and X-rays lagging behind the optical and gamma-ray variability by several hours.

MeV blazars are the most luminous persistent sources in the Universe and emit most of their energy in the MeV band. These objects display very large jet powers and accretion luminosities and are known to host black holes with a mass often exceeding $10^9 M_{\odot}$. An MeV survey, performed by a new generation MeV telescope which will bridge the entire energy and sensitivity gap between the current generation of hard X-ray and gamma-ray instruments, will detect $>$1000 MeV blazars up to a redshift of $z=5-6$. Here we show that this would allow us: 1) to probe the formation and growth mechanisms of supermassive black holes at high redshifts, 2) to pinpoint the location of the emission region in powerful blazars, 3) to determine how accretion and black hole spin interplay to power the jet.

A next generation of Compton and pair telescopes that improve MeV-band detection sensitivity by more than a decade beyond current instrumental capabilities will open up new insights into a variety of astrophysical source classes. Among these are magnetars, the most highly magnetic of the neutron star zoo, which will serve as a prime science target for a new mission surveying the MeV window. This paper outlines the core questions pertaining to magnetars that can be addressed by such a technology. These range from global magnetar geometry and population trends, to incisive probes of hard X-ray emission locales, to providing cosmic laboratories for spectral and polarimetric testing of exotic predictions of QED, principally the prediction of the splitting of photons and magnetic pair creation. Such fundamental physics cannot yet be discerned in terrestrial experiments. State of the art modeling of the persistent hard X-ray tail emission in magnetars is presented to outline the case for powerful diagnostics using Compton polarimeters. The case highlights an inter-disciplinary opportunity to seed discovery at the interface between astronomy and physics.

The constantly improving sensitivity of ground-based and space-borne observatories has made possible the detection of high-energy emission (X-rays and gamma-rays) from several thousands of extragalactic sources. Enormous progress has been made in measuring the continuum flux enabling us to perform imaging, spectral and timing studies. An important remaining challenge for high-energy astronomy is measuring polarization. The capability to measure polarization is being realized currently at X-ray energies (e.g. with IXPE), and sensitive gamma-ray telescopes capable of measuring polarization, such as AMEGO, AdEPT, e-ASTROGAM, etc., are being developed. These future gamma-ray telescopes will probe the radiation mechanisms and magnetic fields of relativistic jets from active galactic nuclei at spatial scales much smaller than the angular resolution achieved with continuum observations of the instrument. In this white paper, we discuss the scientific potentials of high-energy polarimetry, especially gamma-ray polarimetry, including the theoretical implications, and observational technology advances being made. In particular, we will explore the primary scientific opportunities and wealth of information expected from synergy of multi-wavelength polarimetry that will be brought to multi-messenger astronomy.

Active galactic nuclei (AGN) with relativistic jets, powered by gas accretion onto their central supermassive black hole (SMBH), are unique laboratories for studying the physics of matter and elementary particles in extreme conditions that cannot be realized on Earth. For a long time since the discovery of AGN, photons were the only way to probe the underlying physical processes. The recent discovery of a very high energy neutrino, IceCube-170922A, coincident with a flaring blazar, TXS 0506+056, provides the first evidence that AGN jets are multi-messenger sources; they are capable of accelerating hadrons to very high energies, while producing non-thermal EM radiation and high-energy neutrinos. This new era of multi-messenger astronomy, which will mature in the next decade, offers us the unprecedented opportunity to combine more than one messenger to solve some long-standing puzzles of AGN jet physics: How do jets dissipate their energy to accelerate particles? What is the jet total kinetic power? Where and how do jets produce high-energy emission and neutrinos? What physical mechanisms drive the particle acceleration?

We present the Spectroscopic Time-Resolving Observatory for Broadband Energy X-rays (STROBE-X), a probe-class mission concept selected for study by NASA. It combines huge collecting area, high throughput, broad energy coverage, and excellent spectral and temporal resolution in a single facility. STROBE-X offers an enormous increase in sensitivity for X-ray spectral timing, extending these techniques to extragalactic targets for the first time. It is also an agile mission capable of rapid response to transient events, making it an essential X-ray partner facility in the era of time-domain, multi-wavelength, and multi-messenger astronomy. Optimized for study of the most extreme conditions found in the Universe, its key science objectives include: (1) Robustly measuring mass and spin and mapping inner accretion flows across the black hole mass spectrum, from compact stars to intermediate-mass objects to active galactic nuclei. (2) Mapping out the full mass-radius relation of neutron stars using an ensemble of nearly two dozen rotation-powered pulsars and accreting neutron stars, and hence measuring the equation of state for ultradense matter over a much wider range of densities than explored by NICER. (3) Identifying and studying X-ray counterparts (in the post-Swift era) for multiwavelength and multi-messenger transients in the dynamic sky through cross-correlation with gravitational wave interferometers, neutrino observatories, and high-cadence time-domain surveys in other electromagnetic bands. (4) Continuously surveying the dynamic X-ray sky with a large duty cycle and high time resolution to characterize the behavior of X-ray sources over an unprecedentedly vast range of time scales. STROBE-X's formidable capabilities will also enable a broad portfolio of additional science.

A multiwavelength campaign was organized to take place between March and July of 2012. Excellent temporal coverage was obtained with more than 25 instruments, including the MAGIC, FACT and VERITAS Cherenkov telescopes, the instruments on board the Swift and Fermi spacecraft, and the telescopes operated by the GASP-WEBT collaboration. Mrk 501 showed a very high energy (VHE) gamma-ray flux above 0.2 TeV of $\sim$0.5 times the Crab Nebula flux (CU) for most of the campaign. The highest activity occurred on 2012 June 9, when the VHE flux was $\sim$3 CU, and the peak of the high-energy spectral component was found to be at $\sim$2 TeV. This study reports very hard X-ray spectra, and the hardest VHE spectra measured to date for Mrk 501. The fractional variability was found to increase with energy, with the highest variability occurring at VHE, and a significant correlation between the X-ray and VHE bands. The unprecedentedly hard X-ray and VHE spectra measured imply that their low- and high-energy components peaked above 5 keV and 0.5 TeV, respectively, during a large fraction of the observing campaign, and hence that Mrk 501 behaved like an extreme high-frequency- peaked blazar (EHBL) throughout the 2012 observing season. This suggests that being an EHBL may not be a permanent characteristic of a blazar, but rather a state which may change over time. The one-zone synchrotron self-Compton (SSC) scenario can successfully describe the segments of the SED where most energy is emitted, with a significant correlation between the electron energy density and the VHE gamma-ray activity, suggesting that most of the variability may be explained by the injection of high-energy electrons. The one-zone SSC scenario used reproduces the behaviour seen between the measured X-ray and VHE gamma-ray fluxes, and predicts that the correlation becomes stronger with increasing energy of the X-rays.

Hard, non-thermal, persistent pulsed X-ray emission extending between 10 keV and $\sim 150$ keV has been observed in nearly ten magnetars. For inner-magnetospheric models of such emission, resonant inverse Compton scattering of soft thermal photons by ultra-relativistic charges is the most efficient production mechanism. We present angle-dependent upscattering spectra and pulsed intensity maps for uncooled, relativistic electrons injected in inner regions of magnetar magnetospheres, calculated using collisional integrals over field loops. Our computations employ a new formulation of the QED Compton scattering cross section in strong magnetic fields that is physically correct for treating important spin-dependent effects in the cyclotron resonance, thereby producing correct photon spectra. The spectral cut-off energies are sensitive to the choices of observer viewing geometry, electron Lorentz factor, and scattering kinematics. We find that electrons with energies $\lesssim 15$ MeV will emit most of their radiation below 250 keV, consistent with inferred turnovers for magnetar hard X-ray tails. More energetic electrons still emit mostly below 1 MeV, except for viewing perspectives sampling field line tangents. Pulse profiles may be singly- or doubly-peaked dependent upon viewing geometry, emission locale, and observed energy band. Magnetic pair production and photon splitting will attenuate spectra to hard X-ray energies, suppressing signals in the Fermi-LAT band. The resonant Compton spectra are strongly polarized, suggesting that hard X-ray polarimetry instruments such as X-Calibur, or a future Compton telescope, can prove central to constraining model geometry and physics.

We report the analysis of 5 NuSTAR observations of SGR 1806-20 spread over a year from April 2015 to April 2016, more than 11 years following its Giant Flare (GF) of 2004. The source spin frequency during the NuSTAR observations follows a linear trend with a frequency derivative $\dot{\nu}=(-1.25\pm0.03)\times10^{-12}$ Hz s$^{-1}$, implying a surface dipole equatorial magnetic field $B\approx7.7\times10^{14}$ G. Thus, SGR 1806-20 has finally returned to its historical minimum torque level measured between 1993 and 1998. The source showed strong timing noise for at least 12 years starting in 2000, with $\dot{\nu}$ increasing one order of magnitude between 2005 and 2011, following its 2004 major bursting episode and GF. SGR 1806-20 has not shown strong transient activity since 2009 and we do not find short bursts in the NuSTAR data. The pulse profile is complex with a pulsed fraction of $\sim8\%$ with no indication of energy dependence. The NuSTAR spectra are well fit with an absorbed blackbody, $kT=0.62\pm0.06$ keV, plus a power-law, $\Gamma=1.33\pm0.03$. We find no evidence for variability among the 5 observations, indicating that SGR 1806-20 has reached a persistent and potentially its quiescent X-ray flux level after its 2004 major bursting episode. Extrapolating the NuSTAR model to lower energies, we find that the 0.5-10 keV flux decay follows an exponential form with a characteristic timescale $\tau=543\pm75$ days. Interestingly, the NuSTAR flux in this energy range is a factor of $\sim2$ weaker than the long-term average measured between 1993 and 2003, a behavior also exhibited in SGR $1900+14$. We discuss our findings in the context of the magnetar model.

As supernova remnants (SNRs) age, they become efficient cosmic ray accelerators at their outer shell shocks. The current paradigm for shock acceleration theory favors turbulent field environs in the proximity of these shocks, turbulence driven by current instabilities involving energetic ions. With the imminent prospect of dedicated X-ray polarimeters becoming a reality, the possibility looms of probing turbulence on scales that couple to the super-TeV electrons that emit X-rays. This paper presents model X-ray polarization signatures from energetic electrons moving in simulated MHD turbulence of varying levels of "chaos." The emission volumes are finite slabs that represent the active regions of young SNR shells. We find that the turbulent field energy must be quite limited relative to that of the total field in order for the X-ray polarization degree to be as strong as the radio measures obtained in some remnants. Results presented are pertinent to the planned IXPE and XIPE polarimeters.

Pulsed non-thermal quiescent emission between 10 keV and around 150 keV has been observed in $\sim10$ magnetars. For inner magnetospheric models of such hard X-ray signals, resonant Compton upscattering of soft thermal photons from the neutron star surface is the most efficient radiative process. We present angle-dependent hard X-ray upscattering model spectra for uncooled monoenergetic relativistic electrons. The spectral cut-off energies are critically dependent on the observer viewing angles and electron Lorentz factor. We find that electrons with energies less than around 15 MeV will emit most of their radiation below 250 keV, consistent with the observed turnovers in magnetar hard X-ray tails. Moreover, electrons of higher energy still emit most of the radiation below around 1 MeV, except for quasi-equatorial emission locales for select pulses phases. Our spectral computations use new state-of-the-art, spin-dependent formalism for the QED Compton scattering cross section in strong magnetic fields.

Non-thermal quiescent X-ray emission extending between 10 keV and around 150 keV has been seen in about 10 magnetars by RXTE, INTEGRAL, Suzaku, NuSTAR and Fermi-GBM. For inner magnetospheric models of such hard X-ray signals, inverse Compton scattering is anticipated to be the most efficient process for generating the continuum radiation, because the scattering cross section is resonant at the cyclotron frequency. We present hard X-ray upscattering spectra for uncooled monoenergetic relativistic electrons injected in inner regions of pulsar magnetospheres. These model spectra are integrated over bundles of closed field lines and obtained for different observing perspectives. The spectral turnover energies are critically dependent on the observer viewing angles and electron Lorentz factor. We find that electrons with energies less than around 15 MeV will emit most of their radiation below 250 keV, consistent with the turnovers inferred in magnetar hard X-ray tails. Electrons of higher energy still emit most of the radiation below around 1 MeV, except for quasi-equatorial emission locales for select pulse phases. Our spectral computations use a new state-of-the-art, spin-dependent formalism for the QED Compton scattering cross section in strong magnetic fields.

Multiwavelength followup of unidentified Fermi sources has vastly expanded the number of known galactic-field "black widow" and "redback" millisecond pulsar binaries. Focusing on their rotation-powered state, we interpret the radio to X-ray phenomenology in a consistent framework. We advocate the existence of two distinct modes differing in their intrabinary shock orientation, distinguished by the phase-centering of the double-peaked X-ray orbital modulation originating from mildly-relativistic Doppler boosting. By constructing a geometric model for radio eclipses, we constrain the shock geometry as functions of binary inclination and shock stand-off $R_0$. We develop synthetic X-ray synchrotron orbital light curves and explore the model parameter space allowed by radio eclipse constraints applied on archetypal systems B1957+20 and J1023+0038. For B1957+20, from radio eclipses the stand-off is $R_0 \sim 0.15$ -- $0.3$ fraction of binary separation from the companion center, depending on the orbit inclination. Constructed X-ray light curves for B1957+20 using these values are qualitatively consistent with those observed, and we find occultation of the shock by the companion as a minor influence, demanding significant Doppler factors to yield double peaks. For J1023+0038, radio eclipses imply $R_0 \lesssim 0.4$ while X-ray light curves suggest $0.1\lesssim R_0 \lesssim 0.3$ (from the pulsar). Degeneracies in the model parameter space encourage further development to include transport considerations. Generically, the spatial variation along the shock of the underlying electron power-law index should yield energy-dependence in the shape of light curves motivating future X-ray phase-resolved spectroscopic studies to probe the unknown physics of pulsar winds and relativistic shock acceleration therein.

We analyzed broad-band X-ray and radio data of the magnetar SGR J1935+2154 taken in the aftermath of its 2014, 2015, and 2016 outbursts. The source soft X-ray spectrum <10 keV is well described with a BB+PL or 2BB model during all three outbursts. NuSTAR observations revealed a hard X-ray tail, $\Gamma=0.9$, extending up to 79 keV, with flux larger than the one detected <10 keV. Imaging analysis of Chandra data did not reveal small-scale extended emission around the source. Following the outbursts, the total 0.5-10 keV flux from SGR J1935+2154 increased in concordance to its bursting activity, with the flux at activation onset increasing by a factor of $\sim7$ following its strongest June 2016 outburst. A Swift/XRT observation taken 1.5 days prior to the onset of this outburst showed a flux level consistent with quiescence. We show that the flux increase is due to the PL or hot BB component, which increased by a factor of $25$ compared to quiescence, while the cold BB component $kT=0.47$ keV remained more or less constant. The 2014 and 2015 outbursts decayed quasi-exponentially with time-scales of $\sim40$ days, while the stronger May and June 2016 outbursts showed a quick short-term decay with time-scales of $\sim4$ days. Our Arecibo radio observations set the deepest limits on the radio emission from a magnetar, with a maximum flux density limit of 14 $\mu$Jy for the 4.6 GHz observations and 7 $\mu$Jy for the 1.4 GHz observations. We discuss these results in the framework of the current magnetar theoretical models.

We have compiled the most comprehensive burst sample from magnetar 4U 0142+61, comprising 27 bursts from its three burst-active episodes in 2011, 2012 and the latest one in 2015 observed with Swift/BAT and Fermi/GBM. Bursts from 4U 0142+61 morphologically resemble typical short bursts from other magnetars. However, 4U 0142+61 bursts are less energetic compared to the bulk of magnetar bursts. We uncovered an extended tail emission following a burst on 2015 February 28, with a thermal nature, cooling over a time-scale of several minutes. During this tail emission, we also uncovered pulse peak phase aligned X-ray bursts , which could originate from the same underlying mechanism as that of the extended burst tail, or an associated and spatially coincident but different mechanism.

The first wind nebula around a magnetar was recently discovered in X-rays around Swift~J1834.9$-$0846. We study this magnetar's global energetics and the properties of its particle wind or outflows. At a distance of $\sim4\;$kpc, Swift~J1834.9$-$0846 is located at the center of the supernova remnant (SNR) W41 whose radius is $\sim 19\;$pc, an order of magnitude larger than that of the X-ray nebula ($\sim2\;$pc). The association with SNR W41 suggests a common age of $\sim5-100\;$kyr, while its spin-down age is $4.9$~kyr. A small natal kick velocity may partly explain why a wind nebula was detected around this magnetar but not around other magnetars, most of which appear to have larger kick velocities and may have exited their birth SNR. We find that the GeV and TeV source detected by Fermi/LAT and H.E.S.S., respectively, of radius $\sim11\;$pc is most likely of hadronic origin. The dynamics and internal structure of the nebula are examined analytically to explain the nebula's current properties. Its size may naturally correspond to the diffusion-dominated cooling length of the X-ray emitting $e^+e^-$ pairs. This may also account for the spectral softening of the X-ray emission from the nebula's inner to outer parts. Analysis of the X-ray synchrotron nebula implies that (i) the nebular magnetic field is $\gtrsim 11\;\mu$G (and likely $\lesssim30\;\mu$G), and (ii) the nebula is not powered predominantly by the magnetar's quiescent spin-down-powered MHD wind, but by other outflows that contribute most of its energy. The latter are most likely associated with the magnetar's bursting activity, and possibly dominated by outflows associated with its past giant flares. The energy source for the required outflows cannot be the decay of the magnetar's dipole field alone, and is most likely the decay of its much stronger internal magnetic field.

Diffusive shock acceleration (DSA) at relativistic shocks is widely thought to be an important acceleration mechanism in various astrophysical jet sources, including radio-loud active galactic nuclei such as blazars. Such acceleration can produce the non-thermal particles that emit the broadband continuum radiation that is detected from extragalactic jets. An important recent development for blazar science is the ability of Fermi-LAT spectroscopy to pin down the shape of the distribution of the underlying non-thermal particle population. This paper highlights how multi-wavelength spectra spanning optical to X-ray to gamma-ray bands can be used to probe diffusive acceleration in relativistic, oblique, magnetohydrodynamic (MHD) shocks in blazar jets. Diagnostics on the MHD turbulence near such shocks are obtained using thermal and non-thermal particle distributions resulting from detailed Monte Carlo simulations of DSA. These probes are afforded by the characteristic property that the synchrotron $\nu F_{\nu}$ peak energy does not appear in the gamma-ray band above 100 MeV. We investigate self-consistently the radiative synchrotron and inverse Compton signatures of the simulated particle distributions. Important constraints on the diffusive mean free paths of electrons, and the level of electromagnetic field turbulence are identified for three different case study blazars, Mrk 501, BL Lacertae and AO 0235+164. The X-ray excess of AO 0235+164 in a flare state can be modelled as the signature of bulk Compton scattering of external radiation fields, thereby tightly constraining the energy-dependence of the diffusion coefficient for electrons. The concomitant interpretations that turbulence levels decline with remoteness from jet shocks, and the probable significant role for non-gyroresonant diffusion, are posited.

The very high energy (VHE; E $>$ 100 GeV) blazar Markarian 501 was observed between April 17 and May 5 (MJD 54938--54956), 2009, as part of an extensive multi-wavelength campaign from radio to VHE. Strong VHE $\gamma$-ray activity was detected on May 1st with Whipple and VERITAS, when the flux (E $>$ 400 GeV) increased to 10 times the pre-flare baseline flux ($3.9{\times 10^{-11}}~{\rm ph~cm^{-2}~s^{-1}}$), reaching five times the flux of the Crab Nebula. This coincided with a decrease in the optical polarization and a rotation of the polarization angle by 15$^{\circ}$. This VHE flare showed a fast flux variation with an increase of a factor $\sim$4 in 25 minutes, and a falling time of $\sim$50 minutes. We present the observations of the quiescent state previous to the flare and of the high state after the flare, focusing on the flux and spectral variability from Whipple, VERITAS, Fermi-LAT, RXTE, and Swift combined with optical and radio data.

We present the Fermi Gamma-ray Burst Monitor (GBM) and Large Area Telescope (LAT) observations of the LIGO binary black hole merger event GW151226 and candi- date LVT151012. No candidate electromagnetic counterparts were detected by either the GBM or LAT. We present a detailed analysis of the GBM and LAT data over a range of timescales from seconds to years, using automated pipelines and new techniques for char- acterizing the upper limits across a large area of the sky. Due to the partial GBM and LAT coverage of the large LIGO localization regions at the trigger times for both events, dif- ferences in source distances and masses, as well as the uncertain degree to which emission from these sources could be beamed, these non-detections cannot be used to constrain the variety of theoretical models recently applied to explain the candidate GBM counterpart to GW150914.

Magnetar giant flares show oscillatory modulations in the tails of their light curves, which can only be explained via some form of beaming. The fireball model for magnetar bursts has been used successfully to fit the phase-averaged light curves of the tails of giant flares, but so far no attempts have been made to fit the pulsations. We present a relatively simple numerical model to simulate beaming of magnetar flare emission. In our simulations, radiation escapes from the base of a fireball trapped in a dipolar magnetic field, and is scattered through the optically thick magnetosphere of the magnetar until it escapes. Beaming is provided by the presence of a relativistic outflow, as well as by the geometry of the system. We find that a simple picture for the relativistic outflow is enough to create the pulse fraction and sharp peaks observed in pulse profiles of magnetar flares, while without a relativistic outflow the beaming is insufficient to explain giant flare rotational modulations.

Dwarf spheroidal galaxies (dSphs) are promising targets for the indirect detection of dark matter through gamma-ray emission due to their proximity, lack of astrophysical backgrounds and high dark matter density. They are often used to place restrictive bounds on the dark matter annihilation cross section. In this paper, we analyze six years of {\it Fermi}-LAT gamma-ray data from 19 dSphs that are satellites of the Milky Way, and derive from a stacked analysis of 15 dSphs, robust 95\% confidence level lower limits on the dark matter lifetime for several decay channels and dark matter masses between $\sim 1$GeV and $10$TeV. Our findings are based on a bin-by-bin maximum likelihood analysis treating the J-factor as a nuisance parameter using PASS 8 event-class. Our constraints from this ensemble are among the most stringent and solid in the literature, and competitive with existing ones coming from the extragalactic gamma-ray background, galaxy clusters, AMS-02 cosmic ray data, Super-K and ICECUBE neutrino data, while rather insensitive to systematic uncertainties. In particular, among gamma-ray searches, we improve existing limits for dark matter decaying into $\bar{b}b$, ($\mu^+\mu^-$) for DM masses below $\sim 30\, (200)$~GeV, demonstrating that dSphs are compelling targets for constraining dark matter decay lifetimes.

We report on the analysis of two deep XMM-Newton observations of the magnetar Swift J1834.9-0846 and its surrounding extended emission taken in March 2014 and October 2014, 2.5 and 3.1 years after the source went into outburst. The magnetar is only weakly detected in the first observation with an absorption corrected flux $F_{\rm 0.5-10 keV}\approx4\times10^{-14}$ erg s$^{-1}$ cm$^{-2}$, and a $3\sigma$ upper limit during the second observation of about $3\times10^{-14}$ erg s$^{-1}$ cm$^{-2}$. This flux level is more than 3 orders of magnitude lower than the flux measured at the outburst onset on September 2011. The extended emission, centered at the magnetar position and elongated towards the south-west, is clearly seen in both observations; it is best fit by a highly absorbed power-law (PL), with a hydrogen column density of $N_{\rm H}=8.0\times10^{22}$ cm$^{-2}$ and PL photon index $\Gamma=2.2\pm0.2$. Its flux is constant between the two observations at $F_{\rm 0.5-10 keV}=1.3\times10^{-12}$ erg s$^{-1}$ cm$^{-2}$. We find no statistically significant changes in the spectral shape or the flux of this extended emission over a period of 9 years from 2005 to 2014. These new results strongly support the extended emission nature as a wind nebula and firmly establish Swift J1834.9-0846 as the first magnetar to show a surrounding wind nebula. Further, our results imply that such nebulae are no longer exclusive to rotation-powered pulsars and narrow the gap between these two sub-populations of isolated neutron stars. The size and spectrum of the nebula are compatible with those of pulsar-wind nebulae but its radiative efficiency $\eta_{\rm X}=L_{\rm X}/\dot{E}\approx0.1$ is markedly high, possibly pointing to an additional wind component in Swift J1834.9-0846.

This paper describes the Polarization Spectroscopic Telescope Array (PolSTAR), a mission proposed to NASA's 2014 Small Explorer (SMEX) announcement of opportunity. PolSTAR measures the linear polarization of 3-50 keV (requirement; goal: 2.5-70 keV) X-rays probing the behavior of matter, radiation and the very fabric of spacetime under the extreme conditions close to the event horizons of black holes, as well as in and around magnetars and neutron stars. The PolSTAR design is based on the technology developed for the Nuclear Spectroscopic Telescope Array (NuSTAR) mission launched in June 2012. In particular, it uses the same X-ray optics, extendable telescope boom, optical bench, and CdZnTe detectors as NuSTAR. The mission has the sensitivity to measure ~1% linear polarization fractions for X-ray sources with fluxes down to ~5 mCrab. This paper describes the PolSTAR design as well as the science drivers and the potential science return.

Since launch in 2008, the Fermi Gamma-ray Burst Monitor (GBM) has detected many hundreds of bursts from magnetar sources. While the vast majority of these bursts have been attributed to several known magnetars, there is also a small sample of magnetar-like bursts of unknown origin. Here we present the Fermi/GBM magnetar catalog, giving the results of the temporal and spectral analyses of 440 magnetar bursts with high temporal and spectral resolution. This catalog covers the first five years of GBM magnetar observations, from July 2008 to June 2013. We provide durations, spectral parameters for various models, fluences and peak fluxes for all the bursts, as well as a detailed temporal analysis for SGR J1550-5418 bursts. Finally, we suggest that some of the bursts of unknown origin are associated with the newly discovered magnetar 3XMM J185246.6+0033.7.

X-ray polarimetry promises to give qualitatively new information about high-energy astrophysical sources, such as binary black hole systems, micro-quasars, active galactic nuclei, neutron stars, and gamma-ray bursts. We designed, built and tested a X-ray polarimeter, X-Calibur, to be used in the focal plane of the balloon-borne InFOCuS grazing incidence X-ray telescope. X-Calibur combines a low-Z scatterer with a CZT detector assembly to measure the polarization of 20-80keV X-rays making use of the fact that polarized photons scatter preferentially perpendicular to the electric field orientation. X-Calibur achieves a high detection efficiency of ~80%. The X-Calibur detector assembly is completed, tested, and fully calibrated. The response to a polarized X-ray beam was measured successfully at the Cornell High Energy Synchrotron Source. This paper describes the design, calibration and performance of the X-Calibur polarimeter. In principle, a similar space-borne scattering polarimeter could operate over the broader 2-100keV energy band.

The quantum electrodynamical (QED) process of Compton scattering in strong magnetic fields is commonly invoked in atmospheric and inner magnetospheric models of x-ray and soft gamma-ray emission in high-field pulsars and magnetars. A major influence of the field is to introduce resonances at the cyclotron frequency and its harmonics, where the incoming photon accesses thresholds for the creation of virtual electrons or positrons in intermediate states with excited Landau levels. At these resonances, the effective cross section typically exceeds the classical Thomson value by over 2 orders of magnitude. Near and above the quantum critical magnetic field of 44.13 TeraGauss, relativistic corrections must be incorporated when computing this cross section. This paper presents formalism for the QED magnetic Compton differential cross section valid for both subcritical and supercritical fields, yet restricted to scattered photons that are below pair creation threshold. Calculations are developed for the particular case of photons initially propagating along the field, mathematically simple specializations that are germane to interactions involving relativistic electrons frequently found in neutron star magnetospheres. This exposition of relativistic, quantum, magnetic Compton cross sections treats electron spin dependence fully, since this is a critical feature for describing the finite decay lifetimes of the intermediate states. The formalism employs both the Johnson and Lippmann (JL) wave functions and the Sokolov and Ternov (ST) electron eigenfunctions of the magnetic Dirac equation. The ST states are formally correct for self-consistently treating spin-dependent effects that are so important in the resonances. Relatively compact analytic forms for the cross sections are presented that will prove useful for astrophysical modelers.

Magnetic pair creation $\gamma \to e^+e^-$ has been at the core of radio pulsar paradigms and central to polar cap models of gamma-ray pulsars for over three decades. The Fermi gamma-ray pulsar population now exceeds 140 sources and has defined an important part of Fermi's science legacy. Among the population characteristics well established is the common occurrence of exponential turnovers in their spectra in the 1--10 GeV range. These turnovers are too gradual to arise from magnetic pair creation in the strong magnetic fields of pulsar inner magnetospheres. By demanding insignificant photon attenuation precipitated by such single-photon pair creation, the energies of these turnovers for Fermi pulsars can be used to compute lower bounds for the typical altitude of GeV band emission. This paper explores such pair transparency constraints below the turnover energy, and updates earlier altitude bound determinations of that have been deployed in various Fermi pulsar papers. For low altitude emission locales, general relativistic influences are found to be important. Rotational aberration influences are also explored, and are found to be small at low altitudes, except near the magnetic pole. The altitude bounds are typically in the range of 2-7 stellar radii for the young Fermi pulsar population, and provide key information on the emission altitude in radio quiet pulsars that do not possess double-peaked pulse profiles. The bound for the Crab pulsar is at a much higher altitude, with the putative detection by MAGIC out to 350-400 GeV implying a lower bound of 310km to the emission region, i.e., approximately 20% of the light cylinder radius. These results are also extended to the super-critical field domain, where it is found that emission in magnetars originating below around 10 stellar radii will not appear in the Fermi-LAT band.

The discovery of quasi-periodic oscillations (QPOs) in magnetar giant flares has opened up prospects for neutron star asteroseismology. The scarcity of giant flares makes a search for QPOs in the shorter, far more numerous bursts from Soft Gamma Repeaters (SGRs) desirable. In Huppenkothen et al (2013), we developed a Bayesian method for searching for QPOs in short magnetar bursts, taking into account the effects of the complicated burst structure, and have shown its feasibility on a small sample of bursts. Here, we apply the same method to a much larger sample from a burst storm of 286 bursts from SGR J1550-5418. We report a candidate signal at 260 Hz in a search of the individual bursts, which is fairly broad. We also find two QPOs at 93 Hz and one at 127 Hz, when averaging periodograms from a number of bursts in individual triggers, at frequencies close to QPOs previously observed in magnetar giant flares. Finally, for the first time, we explore the overall burst variability in the sample, and report a weak anti-correlation between the power-law index of the broadband model characterising aperiodic burst variability, and the burst duration: shorter bursts have steeper power law indices than longer bursts. This indicates that longer bursts vary over a broader range of time scales, and are not simply longer versions of the short bursts.

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 report on time-resolved spectroscopy of the 63 brightest bursts of SGR J1550-5418, detected with Fermi/Gamma-ray Burst Monitor during its 2008-2009 intense bursting episode. We performed spectral analysis down to 4 ms time-scales, to characterize the spectral evolution of the bursts. Using a Comptonized model, we find that the peak energy, E_peak, anti-correlates with flux, while the low-energy photon index remains constant at -0.8 up to a flux limit F~10^-5 erg s-1 cm-2. Above this flux value the E_peak-flux correlation changes sign, and the index positively correlates with flux reaching 1 at the highest fluxes. Using a two black-body model, we find that the areas and fluxes of the two emitting regions correlate positively. Further, we study here for the first time, the evolution of the temperatures and areas as a function of flux. We find that the area-kT relation follows lines of constant luminosity at the lowest fluxes, R^2 \propto kT^-4, with a break at higher fluxes ($F>10^-5.5 erg s-1 cm-2). The area of the high-kT component increases with flux while its temperature decreases, which we interpret as due to an adiabatic cooling process. The area of the low-kT component, on the other hand, appears to saturate at the highest fluxes, towards R_max~30 km. Assuming that crust quakes are responsible for SGR bursts and considering R_max as the maximum radius of the emitting photon-pair plasma fireball, we relate this saturation radius to a minimum excitation radius of the magnetosphere, and put a lower limit on the internal magnetic field of SGR J1550-5418, B_int>~4.5x10^15 G.

Diffusive shock acceleration (DSA) at relativistic shocks is likely to be an important acceleration mechanism in various astrophysical jet sources, including radio-loud AGN. An important recent development for blazar science is the ability of Fermi-LAT data to pin down the power-law index of the high energy portion of emission in these sources, and therefore also the index of the underlying non-thermal particle population. This diagnostic potential was not possible prior to Fermi launch, when gamma-ray information was dominated by the highly-absorbed TeV band. This paper highlights how multiwavelength spectra including X-ray band and Fermi data can be used to probe diffusive acceleration in relativistic, oblique, MHD shocks in blazar jets. The spectral index of the non-thermal particle distributions resulting from Monte Carlo simulations of DSA, and the fraction of thermal particles accelerated to non-thermal energies, depend sensitively on the particles' mean free path scale, and also on the magnetic field obliquity to the shock normal. We investigate self-consistently the radiative synchrotron/Compton signatures of the resulting thermal and non-thermal particle distributions. Important constraints on the frequency of particle scattering and the level of field turbulence are identified for the blazar AO 0235+164. The possible interpretation that turbulence levels decline with remoteness from jet shocks, and a significant role for non-gyroresonant diffusion, are discussed.

Time-resolved spectroscopy is performed on eight bright, long gamma-ray bursts (GRBs) dominated by single emission pulses that were observed with the {\it Fermi Gamma-ray Space Telescope}. Fitting the prompt radiation of GRBs by empirical spectral forms such as the Band function leads to ambiguous conclusions about the physical model for the prompt radiation. Moreover, the Band function is often inadequate to fit the data. The GRB spectrum is therefore modeled with two emission components consisting of optically thin nonthermal synchrotron radiation from relativistic electrons and, when significant, thermal emission from a jet photosphere, which is represented by a blackbody spectrum. To produce an acceptable fit, the addition of a blackbody component is required in 5 out of the 8 cases. We also find that the low-energy spectral index \alpha is consistent with a synchrotron component with \alpha = -0.81\pm 0.1. This value lies between the limiting values of \alpha = -2/3 and \alpha = -3/2 for electrons in the slow and fast-cooling regimes, respectively, suggesting ongoing acceleration at the emission site. The blackbody component can be more significant when using a physical synchrotron model instead of the Band function, illustrating that the Band function does not serve as a good proxy for a nonthermal synchrotron emission component. The temperature and characteristic emission-region size of the blackbody component are found to, respectively, decrease and increase as power laws with time during the prompt phase. In addition, we find that the blackbody and nonthermal components have separate temporal behaviors.

Gamma-ray burst (GRB) 130427A is one of the most energetic GRBs ever observed. The initial pulse up to 2.5 s is possibly the brightest well-isolated pulse observed to date. A fine time resolution spectral analysis shows power-law decays of the peak energy from the onset of the pulse, consistent with models of internal synchrotron shock pulses. However, a strongly correlated power-law behavior is observed between the luminosity and the spectral peak energy that is inconsistent with curvature effects arising in the relativistic outflow. It is difficult for any of the existing models to account for all of the observed spectral and temporal behaviors simultaneously.

We investigate the possibility of Photospheric Radius Expansion (PRE) during magnetar bursts. Identification of PRE would enable a determination of the magnetic Eddington limit (which depends on field strength and neutron star mass and radius), and shed light on the burst mechanism. To do this we model hydrostatic atmospheres in a strong radial magnetic field, determining both their maximum extent and photospheric temperatures. We find that spatially-extended atmospheres cannot exist in such a field configuration: typical maximum extent for magnetar-strength fields is ~10 m (as compared to 200 km in the non-magnetic case). Achieving balance of gravitational and radiative forces over a large range of radii, which is critical to the existence of extended atmospheres, is rendered impossible in strong fields due to the dependence of opacities on temperature and field strength. We conclude that high luminosity bursts in magnetars do not lead to expansion and cooling of the photosphere, as in the non-magnetic case. We also find the maximum luminosity that can propagate through a hydrostatic magnetar atmosphere to be lower than previous estimates. The proximity and small extent of the photospheres associated with the two different polarization modes also calls into question the interpretation of two blackbody fits to magnetar burst spectra as being due to extended photospheres.

In this document, we describe the scientific potential of blazar observations with a X-ray polarimetry mission like GEMS (Gravity and Extreme Magnetism SMEX). We describe five blazar science investigations that such a mission would enable: (i) the structure and the role of magnetic fields in AGN jets, (ii) analysis of the polarization of the synchrotron X-ray emission from AGN jets, (iii) discrimination between synchrotron self-Compton and external Compton models for blazars with inverse Compton emission in the X-ray band, (iv) a precision study of the polarization properties of the X-ray emission from Cen-A, (v) tests of Lorentz Invariance based on X-ray polarimetric observations of blazars. We conclude with a discussion of a straw man observation program and recommended accompanying multiwavelength observations.

Cosmic rays are particles (mostly protons) accelerated to relativistic speeds. Despite wide agreement that supernova remnants (SNRs) are the sources of galactic cosmic rays, unequivocal evidence for the acceleration of protons in these objects is still lacking. When accelerated protons encounter interstellar material, they produce neutral pions, which in turn decay into gamma rays. This offers a compelling way to detect the acceleration sites of protons. The identification of pion-decay gamma rays has been difficult because high-energy electrons also produce gamma rays via bremsstrahlung and inverse Compton scattering. We detected the characteristic pion-decay feature in the gamma-ray spectra of two SNRs, IC 443 and W44, with the Fermi Large Area Telescope. This detection provides direct evidence that cosmic-ray protons are accelerated in SNRs.

We examine the expected X-ray polarization properties of neutron-star X-ray sources of various types, e.g., accretion and rotation powered pulsars, magnetars, and low-mass X-ray binaries. We summarize the model calculations leading to these expected properties. We describe how a comparison of these with their observed properties, as inferred from GEMS data, will probe the essential dynamical, electromagnetic, plasma, and emission processes in neutron-star binaries, discriminate between models of these processes, and constrain model parameters. An exciting goal is the first observational demonstration in this context of the existence of vacuum resonance, a fundamental quantum electrodynamical phenomenon first described in the 1930s.

We present here a summary of the scientific goals behind the Gravity and Extreme Magnetism SMEX (GEMS) X-ray polarimetry mission's black hole (BH) observing program. The primary targets can be divided into two classes: stellar-mass galactic BHs in accreting binaries, and super-massive BHs in the centers of active galactic nuclei (AGN). The stellar-mass BHs can in turn be divided into various X-ray spectral states: thermal-dominant (disk), hard (radio jet), and steep power-law (hot corona). These different spectral states are thought to be generated by different accretion geometries and emission mechanisms. X-ray polarization is an ideal tool for probing the geometry around these BHs and revealing the specific properties of the accreting gas.

We report the discovery of PSR J1838-0537, a gamma-ray pulsar found through a blind search of data from the Fermi Large Area Telescope (LAT). The pulsar has a spin frequency of 6.9 Hz and a frequency derivative of -2.2e-11 Hz/s, implying a young characteristic age of 4970 years and a large spin-down power of 5.9e36 erg/s. Follow-up observations with radio telescopes detected no pulsations, thus PSR J1838-0537 appears radio-quiet as viewed from Earth. In September 2009 the pulsar suffered the largest glitch so far seen in any gamma-ray-only pulsar, causing a relative increase in spin frequency of about 5.5e-6. After the glitch, during a putative recovery period, the timing analysis is complicated by the sparsity of the LAT photon data, the weakness of the pulsations, and the reduction in average exposure from a coincidental, contemporaneous change in the LAT's sky-survey observing pattern. The pulsar's sky position is coincident with the spatially extended TeV source HESS J1841-055 detected by the High Energy Stereoscopic System (H.E.S.S.). The inferred energetics suggest that HESS J1841-055 contains a pulsar wind nebula powered by the pulsar.

We present the results of our broadband spectral analysis of 42 SGR J1550-5418 bursts simultaneously detected with the Swift/X-ray Telescope (XRT) and the Fermi/Gamma-ray Burst Monitor (GBM), during the 2009 January active episode of the source. The unique spectral and temporal capabilities of the XRT Windowed Timing mode have allowed us to extend the GBM spectral coverage for these events down to the X-ray domain (0.5-10 keV). Our earlier analysis of the GBM data found that the SGR J1550-5418 burst spectra were described equally well with a Comptonized model or with two blackbody functions; the two models were statistically indistinguishable. Our new broadband (0.5 - 200 keV) spectral fits show that, on average, the burst spectra are better described with two blackbody functions than with the Comptonized model. Thus, our joint XRT/GBM analysis clearly shows for the first time that the SGR J1550-5418 burst spectra might naturally be expected to exhibit a more truly thermalized character, such as a two-blackbody or even a multi-blackbody signal. Using the Swift and RXTE timing ephemeris for SGR J1550-5418 we construct the distribution of the XRT burst counts with spin phase and find that it is not correlated with the persistent X-ray emission pulse phase from SGR J1550-5418. These results indicate that the burst emitting sites on the neutron star need not be co-located with hot spots emitting the bulk of the persistent X-ray emission. Finally, we show that there is a significant pulse phase dependence of the XRT burst counts, likely demonstrating that the surface magnetic field of SGR J1550-5418 is not uniform over the emission zone, since it is anticipated that regions with stronger surface magnetic field could trigger bursts more efficiently.

The extremely efficient process of resonant Compton upscattering by relativistic electrons in high magnetic fields is believed to be a leading emission mechanism of high field pulsars and magnetars in the production of intense X-ray radiation. New analytic developments for the Compton scattering cross section using Sokolov & Ternov (S&T) states with spin-dependent resonant widths are presented. These new results display significant numerical departures from both the traditional cross section using spin-averaged widths, and also from the spin-dependent cross section that employs the Johnson & Lippmann (J&L) basis states, thereby motivating the astrophysical deployment of this updated resonant Compton formulation. Useful approximate analytic forms for the cross section in the cyclotron resonance are developed for S&T basis states. These calculations are applied to an inner magnetospheric model of the hard X-ray spectral tails in magnetars, recently detected by RXTE and INTEGRAL. Relativistic electrons cool rapidly near the stellar surface in the presence of intense baths of thermal X-ray photons. We present resonant Compton cooling rates for electrons, and the resulting photon spectra at various magnetospheric locales, for magnetic fields above the quantum critical value. These demonstrate how this scattering mechanism has the potential to produce the characteristically flat spectral tails observed in magnetars.

In early October 2008, the Soft Gamma Repeater SGRJ1550 - 5418 (1E 1547.0 - 5408, AXJ155052 - 5418, PSR J1550 - 5418) became active, emitting a series of bursts which triggered the Fermi Gamma-ray Burst Monitor (GBM) after which a second especially intense activity period commenced in 2009 January and a third, less active period was detected in 2009 March-April. Here we analyze the GBM data all the bursts from the first and last active episodes. We performed temporal and spectral analysis for all events and found that their temporal characteristics are very similar to the ones of other SGR bursts, as well the ones reported for the bursts of the main episode (average burst durations \sim 170 ms). In addition, we used our sample of bursts to quantify the systematic uncertainties of the GBM location algorithm for soft gamma-ray transients to < 8 deg. Our spectral analysis indicates significant spectral evolution between the first and last set of events. Although the 2008 October events are best fit with a single blackbody function, for the 2009 bursts an Optically Thin Thermal Bremsstrahlung (OTTB) is clearly preferred. We attribute this evolution to changes in the magnetic field topology of the source, possibly due to effects following the very energetic main bursting episode.

In this white paper, we discuss the concept of a next-generation X-ray mission called BEST (Black hole Evolution and Space Time). The mission concept uses a 3000 square centimeter effective area mirror (at 6 keV) to achieve unprecedented sensitivities for hard X-ray imaging spectrometry (5-70 keV) and for broadband X-ray polarimetry (2-70 keV). BEST can make substantial contributions to our understanding of the inner workings of accreting black holes, our knowledge about the fabric of extremely curved spacetime, and the evolution of supermassive black holes. BEST will allow for time resolved studies of accretion disks. With a more than seven times larger mirror area and a seven times wider bandpass than GEMS, BEST will take X-ray polarimetry to a new level: it will probe the time variability of the X-ray polarization from stellar mass and supermassive black holes, and it will measure the polarization properties in 30 independent energy bins. These capabilities will allow BEST to conduct tests of accretion disk models and the underlying spacetimes. With three times larger mirror area and ten times better angular resolution than NuSTAR, BEST will be able to make deep field observations with a more than 15 times better sensitivity than NuSTAR. The mission will be able to trace the evolution of obscured and unobscured black holes in the redshift range from zero to six, covering the most important epoch of supermassive black hole growth. The hard X-ray sensitivity of BEST will enable a deep census of non-thermal particle populations. BEST will give us insights into AGN feedback by measuring the particle luminosity injected by AGNs into the interstellar medium (ISM) of their hosts, and will map the emission from particles accelerated at large scale structure shocks. Finally, BEST has the potential to constrain the equation of state of neutron stars (NS).

We have performed detailed temporal and time-integrated spectral analysis of 286 bursts from SGR J1550-5418 detected with the Fermi Gamma-ray Burst Monitor (GBM) in January 2009, resulting in the largest uniform sample of temporal and spectral properties of SGR J1550-5418 bursts. We have used the combination of broadband and high time-resolution data provided with GBM to perform statistical studies for the source properties. We determine the durations, emission times, duty cycles and rise times for all bursts, and find that they are typical of SGR bursts. We explore various models in our spectral analysis, and conclude that the spectra of SGR J1550-5418 bursts in the 8-200 keV band are equally well described by optically thin thermal bremsstrahlung (OTTB), a power law with an exponential cutoff (Comptonized model), and two black-body functions (BB+BB). In the spectral fits with the Comptonized model we find a mean power-law index of -0.92, close to the OTTB index of -1. We show that there is an anti-correlation between the Comptonized Epeak and the burst fluence and average flux. For the BB+BB fits we find that the fluences and emission areas of the two blackbody functions are correlated. The low-temperature BB has an emission area comparable to the neutron star surface area, independent of the temperature, while the high-temperature blackbody has a much smaller area and shows an anti-correlation between emission area and temperature. We compare the properties of these bursts with bursts observed from other SGR sources during extreme activations, and discuss the implications of our results in the context of magnetar burst models.

The Large Area Telescope on board the Fermi satellite observed a gamma-ray flare in the Crab nebula lasting for approximately nine days in April of 2011. The source, which at optical wavelengths has a size of ~11 ly across, doubled its gamma-ray flux within eight hours. The peak photon flux was (186 +- 6) 10-7 cm-2 s-1 above 100 MeV, which corresponds to a 30-fold increase compared to the average value. During the flare, a new component emerged in the spectral energy distribution, which peaked at an energy of (375 +- 26) MeV at flare maximum. The observations imply that the emission region was likely relativistically beamed toward us and that variations in its motion are responsible for the observed spectral variability.

Diffusive shock acceleration (DSA) at relativistic shocks is expected to be an important acceleration mechanism in a variety of astrophysical objects including extragalactic jets in active galactic nuclei and gamma ray bursts. These sources remain good candidate sites for the generation of ultra-high energy cosmic rays. In this paper, key predictions of DSA at relativistic shocks that are germane to production of relativistic electrons and ions are outlined. The technique employed to identify these characteristics is a Monte Carlo simulation of such diffusive acceleration in test-particle, relativistic, oblique, magnetohydrodynamic (MHD) shocks. Using a compact prescription for diffusion of charges in MHD turbulence, this approach generates particle angular and momentum distributions at any position upstream or downstream of the shock. Simulation output is presented for both small angle and large angle scattering scenarios, and a variety of shock obliquities including superluminal regimes when the de Hoffmann-Teller frame does not exist. The distribution function power-law indices compare favorably with results from other techniques. They are found to depend sensitively on the mean magnetic field orientation in the shock, and the nature of MHD turbulence that propagates along fields in shock environs. An interesting regime of flat spectrum generation is addressed; we provide evidence for it being due to shock drift acceleration, a phenomenon well-known in heliospheric shock studies. The impact of these theoretical results on blazar science is outlined. Specifically, Fermi-LAT gamma-ray observations of these relativistic jet sources are providing significant constraints on important environmental quantities for relativistic shocks, namely the field obliquity, the frequency of scattering and the level of field turbulence.

X-ray polarimetry promises to give new information about high-energy astrophysical sources, such as binary black hole systems, micro-quasars, active galactic nuclei, and gamma-ray bursts. We designed, built and tested a hard X-ray polarimeter X-Calibur to be used in the focal plane of the InFOCuS grazing incidence hard X-ray telescope. X-Calibur combines a low-Z Compton scatterer with a CZT detector assembly to measure the polarization of 10-80 keV X-rays making use of the fact that polarized photons Compton scatter preferentially perpendicular to the electric field orientation. X-Calibur achieves a high detection efficiency of order unity.

X-ray polarimetry promises to give qualitatively new information about high-energy astrophysical sources, such as binary black hole systems, micro-quasars, active galactic nuclei, and gamma-ray bursts. We designed, built and tested a hard X-ray polarimeter X-Calibur to be used in the focal plane of the InFOCuS grazing incidence hard X-ray telescope. X-Calibur combines a low-Z Compton scatterer with a CZT detector assembly to measure the polarization of 10-80 keV X-rays making use of the fact that polarized photons Compton scatter preferentially perpendicular to the electric field orientation. X-Calibur achieves a high detection efficiency of order unity.

We report the detection of gamma-ray pulsations from the high-magnetic-field rotation-powered pulsar PSR J1119-6127 using data from the Fermi Large Area Telescope. The gamma-ray light curve of PSR J1119-6127 shows a single, wide peak offset from the radio peak by 0.43 pm 0.02 in phase. Spectral analysis suggests a power law of index 1.0 pm 0.3 with an energy cut-off at 0.8 pm 0.2 GeV. The first uncertainty is statistical and the second is systematic. We discuss the emission models of PSR J1119-6127 and demonstrate that despite the object's high surface magnetic field---near that of magnetars---the field strength and structure in the gamma-ray emitting zone are apparently similar to those of typical young pulsars. Additionally, we present upper limits on the \gam-ray pulsed emission for the magnetically active PSR J1846-0258 in the supernova remnant Kesteven 75 and two other energetic high-B pulsars, PSRs J1718-3718 and J1734-3333. We explore possible explanations for the non-detection of these three objects, including peculiarities in their emission geometry.

Swift/BAT detected the first burst from 1E 1841-045 in May 2010 with intermittent burst activity recorded through at least July 2011. Here we present Swift and Fermi/GBM observations of this burst activity and search for correlated changes to the persistent X-ray emission of the source. The T90 durations of the bursts range between 18-140 ms, comparable to other magnetar burst durations, while the energy released in each burst ranges between (0.8 - 25)E38 erg, which is in the low side of SGR bursts. We find that the bursting activity did not have a significant effect on the persistent flux level of the source. We argue that the mechanism leading to this sporadic burst activity in 1E 1841-045 might not involve large scale restructuring (either crustal or magnetospheric) as seen in other magnetar sources.

Discerning the radiative dissipation mechanism for prompt emission in Gamma-Ray Bursts (GRBs) requires detailed spectroscopic modeling that straddles the $\nu F_{\nu}$ peak in the 100 keV - 1 MeV range. Historically, empirical fits such as the popular Band function have been employed with considerable success in interpreting the observations. While extrapolations of the Band parameters can provide some physical insight into the emission mechanisms responsible for GRBs, these inferences do not provide a unique way of discerning between models. By fitting physical models directly this degeneracy can be broken, eliminating the need for empirical functions; our analysis here offers a first step in this direction. One of the oldest, and leading, theoretical ideas for the production of the prompt signal is the synchrotron shock model (SSM). Here we explore the applicability of this model to a bright {\it Fermi} GBM burst with a simple temporal structure, GRB {\it 090820}A. Our investigation implements, for the first time, thermal and non-thermal synchrotron emissivities in the RMFIT forward-folding spectral analysis software often used in GBM burst studies. We find that these synchrotron emissivities, together with a blackbody shape, provide at least as good a match with the data as the Band GRB spectral fitting function. This success is achieved in both time-integrated and time-resolved spectral fits.

We present our temporal and spectral analyses of 29 bursts from SGR J0501+4516, detected with the Gamma-ray Burst Monitor onboard the Fermi Gamma-ray Space Telescope during the 13 days of the source activation in 2008 (August 22 to September 3). We find that the T90 durations of the bursts can be fit with a log-normal distribution with a mean value of ~ 123 ms. We also estimate for the first time event durations of Soft Gamma Repeater (SGR) bursts in photon space (i.e., using their deconvolved spectra) and find that these are very similar to the T90s estimated in count space (following a log-normal distribution with a mean value of ~ 124 ms). We fit the time-integrated spectra for each burst and the time-resolved spectra of the five brightest bursts with several models. We find that a single power law with an exponential cutoff model fits all 29 bursts well, while 18 of the events can also be fit with two black body functions. We expand on the physical interpretation of these two models and we compare their parameters and discuss their evolution. We show that the time-integrated and time-resolved spectra reveal that Epeak decreases with energy flux (and fluence) to a minimum of ~30 keV at F=8.7e-6 erg/cm2/s, increasing steadily afterwards. Two more sources exhibit a similar trend: SGRs J1550-5418 and 1806-20. The isotropic luminosity corresponding to these flux values is roughly similar for all sources (0.4-1.5 e40 erg/s).

For inner magnetospheric models of hard X-ray and gamma-ray emission in high-field pulsars and magnetars, resonant Compton upscattering is anticipated to be the most efficient process for generating continuum radiation. This is due in part to the proximity of a hot soft photon bath from the stellar surface to putative radiation dissipation regions in the inner magnetosphere. Moreover, because the scattering process becomes resonant at the cyclotron frequency, the effective cross section exceeds the classical Thomson value by over two orders of magnitude, thereby enhancing the efficiency of continuum production and the cooling of relativistic electrons. This paper presents computations of the electron cooling rates for this process, which are needed for resonant Compton models of non-thermal radiation from such highly-magnetized pulsars. The computed rates extend previous calculations of magnetic Thomson cooling to the domain of relativistic quantum effects, sampled near and above the quantum critical magnetic field of 44.13 TeraGauss. This is the first exposition of fully relativistic, quantum magnetic Compton cooling rates for electrons, and it employs both the traditional Johnson and Lippman cross section, and a newer Sokolov and Ternov (ST) formulation of Compton scattering in strong magnetic fields. Such ST formalism is formally correct for treating spin-dependent effects that are important in the cyclotron resonance, and has not been addressed before in the context of cooling by Compton scattering. The QED effects are observed to profoundly lower the rates below extrapolations of the familiar magnetic Thomson results, as expected, when recoil and Klein-Nishina reductions become important.

We report on the first year of Fermi gamma-ray observations of pulsed high-energy emission from the old PSR J2043+2740. The study of the gamma-ray efficiency of such old pulsars gives us an insight into the evolution of pulsars' ability to emit in gammma rays as they age. The gamma-ray lightcurve of this pulsar above 0.1 GeV is clearly defined by two sharp peaks, 0.353+/-0.035 periods apart. We have combined the gamma-ray profile characteristics of PSR J2043+2740 with the geometrical properties of the pulsar's radio emission, derived from radio polarization data, and constrained the pulsar-beam geometry in the framework of a Two Pole Caustic and an Outer Gap model. The ranges of magnetic inclination and viewing angle were determined to be {alpha,zeta}~{52-57,61-68} for the Two Pole Caustic model, and {alpha,zeta}~{62-73,74-81} and {alpha,zeta}~{72-83,60-75} for the Outer Gap model. Based on this geometry, we assess possible birth locations for this pulsar and derive a likely proper motion, sufficiently high to be measurable with VLBI. At a characteristic age of 1.2 Myr, PSR J2043+2740 is the third oldest of all discovered, non-recycled, gamma-ray pulsars: it is twice as old as the next oldest, PSR J0357+32, and younger only than the recently discovered PSR J1836+5925 and PSR J2055+25, both of which are at least 5 and 10 times less energetic, respectively.

Pulsars are powerful sources of radiation across the electromagnetic spectrum. This paper highlights some theoretical insights into non-thermal, magnetospheric pulsar gamma-ray radiation. These advances have been driven by NASA's Fermi mission, launched in mid-2008. The Large Area Telescope (LAT) instrument on Fermi has afforded the discrimination between polar cap and slot gap/outer gap acceleration zones in young and middle-aged pulsars. Altitude discernment using the highest energy pulsar photons will be addressed, as will spectroscopic interpretation of the primary radiation mechanism in the LAT band, connecting to both polar cap/slot gap and outer gap scenarios. Focuses will mostly be on curvature radiation and magnetic pair creation, including population trends that may afford probes of the magnetospheric accelerating potential.

X-ray polarimetry promises to give qualitatively new information about high-energy sources. Examples of interesting source classes are binary black hole systems, rotation and accretion powered neutron stars, Microquasars, Active Galactic Nuclei and Gamma-Ray Bursts. Furthermore, X-ray polarimetry affords the possibility for testing fundamental physics, e.g. to observe signatures of light bending in the strong gravitational field of a black hole, to detect third order Quantum Electrodynamic effects in the magnetosphere of Magnetars, and to perform sensitive tests of Lorentz Invariance. In this paper we discuss scientific drivers of hard (>10 keV) X-ray polarimetry emphasizing how observations in the hard band can complement observations at lower energies (0.1 - 10 keV). Subsequently, we describe four different technical realizations of hard X-ray polarimeters suitable for small to medium sized space borne missions, and study their performance in the signal-dominated case based on Monte Carlo simulations. We end with confronting the instrument requirements for accomplishing the science goals with the capabilities of the four polarimeters.

We report on the findings of a 364 ksec observation of the BL LAC object Mrk 421 with the X-ray observatory Suzaku. The analysis in this paper uses fluxes and hardness ratios in the broad energy range from 0.5 keV to 30 keV. During the course of the observation, the 0.5 keV - 30 keV flux decreased by a factor of $\sim$2 and was accompanied by several large flares occurring on timescales of a few hours. We find that fitting a broken power model to spectra from isolated epochs during the observation describes the data well. Different flares exhibit different spectral and hardness ratio evolutions. The cumulative observational evidence indicates that the particle acceleration mechanism in the Mrk 421 jet produces electron energy distributions with a modest range of spectral indices and maximum energies. We argue that the short-timescale X-ray spectral variability in the flares can be attributed mostly to intrinsic changes in the acceleration process, dominating other influences such as fluctuations in the Doppler beaming factor, or radiative cooling in or outside the acceleration zone.

Pulsar Wind Nebulae (PWNe) act as calorimeters for the relativistic pair winds emanating from within the pulsar light cylinder. Their radiative dissipation in various wavebands is significantly different from that of their pulsar central engines: the broadband spectra of PWNe possess characteristics distinct from those of pulsars, thereby demanding a site of lepton acceleration remote from the pulsar magnetosphere. A principal candidate for this locale is the pulsar wind termination shock, a putatively highly-oblique, ultra-relativistic MHD discontinuity. This paper summarizes key characteristics of relativistic shock acceleration germane to PWNe, using predominantly Monte Carlo simulation techniques that compare well with semi-analytic solutions of the diffusion-convection equation. The array of potential spectral indices for the pair distribution function is explored, defining how these depend critically on the parameters of the turbulent plasma in the shock environs. Injection efficiencies into the acceleration process are also addressed. Informative constraints on the frequency of particle scattering and the level of field turbulence are identified using the multiwavelength observations of selected PWNe. These suggest that the termination shock can be comfortably invoked as a principal injector of energetic leptons into PWNe without resorting to unrealistic properties for the shock layer turbulence or MHD structure.

On August 24th 2008 the new magnetar SGR 0501+4516 (discovered by SWIFT) emitted a bright burst with a pronounced double-peak structure in hard X-rays, reminiscent of the double-peak temporal structure seen in some bright thermonuclear bursts on accreting neutron stars. In the latter case this is due to Photospheric Radius Expansion (PRE): when the flux reaches the Eddington limit, the photosphere expands and cools so that emission becomes softer and drops temporarily out of the X-ray band, re-appearing as the photosphere settles back down. We consider the factors necessary to generate double-peaked PRE events, and show that such a mechanism could plausibly operate in magnetar bursts, despite the vastly different emission process. Identification of the magnetic Eddington limit in a magnetar would constrain magnetic field and distance and could, in principle, enable a measurement of gravitational redshift. It would also locate the emitting region at the neutron star surface, constraining the burst trigger mechanism. Conclusive confirmation of PRE events will require more detailed radiative models for bursts. However for SGR 0501+4516 the predicted critical flux (using the magnetic field strength inferred from timing and the distance suggested by its probable location in the Perseus arm of our Galaxy) is consistent with that observed in the August 24th burst.

Gamma-ray astronomy presents an extraordinary scientific potential for the study of the most powerful sources and the most violent events in the Universe. In order to take full advantage of this potential, the next generation of instrumentation for this domain will have to achieve an improvement in sensitivity over present technologies of at least an order of magnitude. The DUAL mission concept takes up this challenge in two complementary ways: a very long observation of the entire sky, combined with a large collection area for simultaneous observations of Type Ia SNe. While the Wide-Field Compton Telescope (WCT) accumulates data from the full gamma-ray sky (0.1-10 MeV) over the entire mission lifetime, the Laue-Lens Telescope (LLT) focuses on 56Co emission from SNe Ia (0.8-0.9 MeV), collecting gamma-rays from its large area crystal lens onto the WCT. Two separated spacecraft flying in formation will maintain the DUAL payloads at the lens' focal distance.

We report the discovery of eight gamma-ray pulsars in blind frequency searches using the LAT, onboard the Fermi Gamma-ray Space Telescope. Five of the eight pulsars are young (tau_c<100 kyr), energetic (Edot>10^36 erg/s), and located within the Galactic plane (|b|<3 deg). The remaining three are older, less energetic, and located off the plane. Five pulsars are associated with sources included in the LAT bright gamma-ray source list, but only one, PSR J1413-6205, is clearly associated with an EGRET source. PSR J1023-5746 has the smallest characteristic age (tau_c=4.6 kyr) and is the most energetic (Edot=1.1E37 erg/s) of all gamma-ray pulsars discovered so far in blind searches. PSRs J1957+5033 and J2055+25 have the largest characteristic ages (tau_c~1 Myr) and are the least energetic (Edot~5E33 erg/s) of the newly-discovered pulsars. We present the timing models, light curves, and detailed spectral parameters of the new pulsars. We used recent XMM observations to identify the counterpart of PSR J2055+25 as XMMU J205549.4+253959. In addition, publicly available archival Chandra X-ray data allowed us to identify the likely counterpart of PSR J1023-5746 as a faint, highly absorbed source, CXOU J102302.8-574606. The large X-ray absorption indicates that this could be among the most distant gamma-ray pulsars detected so far. PSR J1023-5746 is positionally coincident with the TeV source HESS J1023-575, located near the young stellar cluster Westerlund 2, while PSR J1954+2836 is coincident with a 4.3 sigma excess reported by Milagro at a median energy of 35 TeV. Deep radio follow-up observations of the eight pulsars resulted in no detections of pulsations and upper limits comparable to the faintest known radio pulsars, indicating that these can be included among the growing population of radio-quiet pulsars in our Galaxy being uncovered by the LAT, and currently numbering more than 20.

We report on the design and estimated performance of a balloon-borne hard X-ray polarimeter called HX-POL. The experiment uses a combination of Si and Cadmium Zinc Telluride detectors to measure the polarization of 50 keV-400 keV X-rays from cosmic sources through the dependence of the angular distribution of Compton scattered photons on the polarization direction. On a one-day balloon flight, HX-POL would allow us to measure the polarization of bright Crab-like sources for polarization degrees well below 10%. On a longer (15-30 day) flight from Australia or Antarctica, HX-POL would be be able to measure the polarization of bright galactic X-ray sources down to polarization degrees of a few percent. Hard X-ray polarization measurements provide unique venues for the study of particle acceleration processes by compact objects and relativistic outflows. In this paper, we discuss the overall instrument design and performance. Furthermore, we present results from laboratory tests of the Si and CZT detectors.

It is widely accepted that the prompt transient signal in the 10 keV - 10 GeV band from gamma-ray bursts (GRBs) arises from multiple shocks internal to the ultra-relativistic expansion. The detailed understanding of the dissipation and accompanying acceleration at these shocks is a currently topical subject. This paper explores the relationship between GRB prompt emission spectra and the electron (or ion) acceleration properties at the relativistic shocks that pertain to GRB models. The focus is on the array of possible high-energy power-law indices in accelerated populations, highlighting how spectra above 1 MeV can probe the field obliquity in GRB internal shocks, and the character of hydromagnetic turbulence in their environs. It is emphasized that diffusive shock acceleration theory generates no canonical spectrum at relativistic MHD discontinuities. This diversity is commensurate with the significant range of spectral indices discerned in prompt burst emission. Such system diagnostics are now being enhanced by the broadband spectral coverage of bursts by the Fermi Gamma-Ray Space Telescope; while the Gamma-Ray Burst Monitor (GBM) provides key diagnostics on the lower energy portions of the particle population, the focus here is on constraints in the non-thermal, power-law regime of the particle distribution that are provided by the Large Area Telescope (LAT).

Diffusive shock acceleration (DSA) at relativistic shocks is expected to be an important acceleration mechanism in a variety of astrophysical objects including extragalactic jets in active galactic nuclei and gamma ray bursts. These sources remain strong and interesting candidate sites for the generation of ultra-high energy cosmic rays. In this paper, key predictions of DSA at relativistic shocks that are salient to the issue of cosmic ray ion and electron production are outlined. Results from a Monte Carlo simulation of such diffusive acceleration in test-particle, relativistic, oblique, MHD shocks are presented. Simulation output is described for both large angle and small angle scattering scenarios, and a variety of shock obliquities including superluminal regimes when the de Hoffman-Teller frame does not exist. The distribution function power-law indices compare favorably with results from other techniques. They are found to depend sensitively on the mean magnetic field orientation in the shock, and the nature of MHD turbulence that propagates along fields in shock environs. An interesting regime of flat spectrum generation is addressed, providing evidence for its origin being due to shock drift acceleration. The impact of these theoretical results on gamma-ray burst and blazar science is outlined. Specifically, Fermi gamma-ray observations of these cosmic sources are already providing significant constraints on important environmental quantities for relativistic shocks, namely the frequency of scattering and the level of field turbulence.

The LAT instrument on the Fermi mission will reveal the rich spectral and temporal gamma-ray burst phenomena in the > 100 MeV band. The synergy with Fermi's GBM detectors will link these observations to those in the well explored 10-1000 keV range; the addition of the > 100 MeV band observations will resolve theoretical uncertainties about burst emission in both the prompt and afterglow phases. Trigger algorithms will be applied to the LAT data both onboard the spacecraft and on the ground. The sensitivity of these triggers will differ because of the available computing resources onboard and on the ground. Here we present the LAT's burst detection methodologies and the instrument's GRB capabilities.

The cores of neutron stars harbor the highest matter densities known to occur in nature, up to several times the densities in atomic nuclei. Similarly, magnetic field strengths can exceed the strongest fields generated in terrestrial laboratories by ten orders of magnitude. Hyperon-dominated matter, deconfined quark matter, superfluidity, even superconductivity are predicted in neutron stars. Similarly, quantum electrodynamics predicts that in strong magnetic fields the vacuum becomes birefringent. The properties of matter under such conditions is governed by Quantum Chromodynamics (QCD) and Quantum Electrodynamics (QED), and the close study of the properties of neutron stars offers the unique opportunity to test and explore the richness of QCD and QED in a regime that is utterly beyond the reach of terrestrial experiments. Experimentally, this is almost virgin territory.

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.

The principal paradigm for gamma-ray bursts (GRBs) suggests that the prompt transient gamma-ray signal arises from multiple shocks internal to the relativistic expansion. This paper explores how GRB prompt emission spectra can constrain electron (or ion) acceleration properties at the relativistic shocks that pertain to GRB models. The array of possible high-energy power-law indices in accelerated populations is highlighted, focusing on how spectra above 1 MeV can probe the field obliquity in GRB internal shocks, and the character of hydromagnetic turbulence in their environs. When encompassing the MeV-band spectral break, fits to BATSE/EGRET burst data indicate that the preponderance of electrons responsible for the prompt emission reside in an intrinsically non-thermal population. This differs markedly from typical populations generated in acceleration simulations; potential resolutions of this conflict such as the action of self-absorption are mentioned. Spectral modeling also suggests that the synchrotron mechanism is favored over synchrotron self-Compton scenarios due to the latter's typically broad curvature near the peak. Such diagnostics will be enhanced by the broadband spectral coverage of bursts by the Fermi Gamma-Ray Space Telescope; the GBM will provide key information on the lower energy portions of the non-thermal particle population, while the LAT will constrain the power-law regime of particle acceleration.

The Be/X-ray binary system A0535+262 underwent a giant outburst in May-June 2005, followed by a dimmer outburst in August-September 2005. This increased intensity provided an opportunity to search for redshifted neutron-capture lines from the surface of the neutron star. If discovered, such lines would constrain the neutron star equation of state, providing the motivation of this search. The spectrometer (SPI) on board the INTEGRAL satellite observed the dimmer outburst and provided the data for this research. We have not detected a line with enough significance, with the width-dependent upper limits on the broadened and redshifted neutron capture line in the range of (2 - 11) x 10^(-4) photons cm^(-2) s^(-1). To our knowledge, these are the strongest upper limits on the redshifted 2.2 MeV emission from an accreting neutron star. Our analysis of the transparency of the neutron star surface for 2.2 MeV photons shows that photons have a small but finite chance of leaving the atmosphere unscattered, which diminishes the possibility of detection.

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

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.

The acceleration of interstellar pick-up ions as well as solar wind species has been observed at a multitude of interplanetary (IP) shocks by different spacecraft. This paper expands upon previous work modeling the phase space distributions of accelerated ions associated with the shock event encountered on day 292 of 1991 by the Ulysses mission at 4.5 AU. A kinetic Monte Carlo simulation is employed here to model the diffusive acceleration process. This exposition presents recent developments pertaining to the incorporation into the simulation of the diffusive characteristics incurred by field line wandering (FLW), according to the work of Giacalone and Jokipii. For a pure field-line wandering construct, it is determined that the upstream spatial ramp scales are too short to accommodate the HI-SCALE flux increases for 200 keV protons, and that the distribution function for H+ somewhat underpopulates the combined SWICS/HI-SCALE spectra at the shock. This contrasts our earlier theory/data comparison where it was demonstrated that diffusive transport in highly turbulent fields according to kinetic theory can successfully account for both the proton distributions and upstream ramp scales, using a single turbulence parameter. The principal conclusion here is that, in a FLW scenario, the transport of ions across the mean magnetic field is slightly less efficient than is required to effectively trap energetic ions within a few Larmor radii of the shock layer and thereby precipitate efficient acceleration. This highlights the contrast between ion transport in highly turbulent shock environs and remote, less-disturbed interplanetary regions.

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 magnetospheres of neutron stars provide a valuable testing ground for as-yet unverified theoretical predictions of quantum electrodynamics (QED) in strong electromagnetic fields. Exhibiting magnetic field strengths well in excess of a TeraGauss, such compact astrophysical environments permit the action of exotic mechanisms that are forbidden by symmetries in field-free regions. Foremost among these processes are single-photon pair creation, where a photon converts to an electron-positron pair, and magnetic photon splitting, where a single photon divides into two of lesser energy via the coupling to the external field. The pair conversion process is exponentially small in weak fields, and provides the leading order contribution to vacuum polarization. In contrast, photon splitting possesses no energy threshold and can operate in kinematic regimes where the lower order pair conversion is energetically forbidden. This paper outlines some of the key physical aspects of these processes, and highlights their manifestation in neutron star magnetospheres. Anticipated observational signatures include profound absorption turnovers in pulsar spectra at gamma-ray wavelengths. The shapes of these turnovers provide diagnostics on the possible action of pair creation and the geometrical locale of the photon emission region. There is real potential for the first confirmation of strong field QED with the new GLAST mission, to be launched by NASA in 2008. Suppression of pair creation by photon splitting and its implications for pulsars is also discussed.

The paradigm for Anomalous X-ray Pulsars (AXPs) has evolved recently with the discovery by INTEGRAL and RXTE of flat, hard X-ray components in three AXPs. These non-thermal spectral components differ dramatically from the steeper quasi-power-law tails seen in the classic X-ray band in these sources, and can naturally be attributed to activity in the magnetosphere. Resonant, magnetic Compton upscattering is a candidate mechanism for generating this new component, since it is very efficient in the strong fields present near AXP surfaces. In this paper, results from an inner magnetospheric model for upscattering of surface thermal X-rays in AXPs are presented, using a kinetic equation formalism and employing a QED magnetic scattering cross section. Characteristically flat and strongly-polarized emission spectra are produced by non-thermal electrons injected in the emission region. Spectral results depend strongly on the observer's orientation and the magnetospheric locale of the scattering, which couple directly to the angular distributions of photons sampled. Constraints imposed by the Comptel upper bounds for these AXPs are mentioned.

The observed spectra of blazars, their intrinsic emission, and the underlying populations of radiating particles are intimately related. The use of these sources as probes of the extragalactic infrared background, a prospect propelled by recent advances in TeV-band telescopes, soon to be augmented by observations by NASA's upcoming Gamma-Ray Large Area Space Telescope (GLAST), has been a topic of great recent interest. Here, it is demonstrated that if particles in blazar jets are accelerated at relativistic shocks, then gamma-ray spectra with indices less than 1.5 can be produced. This, in turn, loosens the upper limits on the near infrared extragalactic background radiation previously proposed. We also show evidence hinting that TeV blazars with flatter spectra have higher intrinsic TeV gamma-ray luminosities and we indicate that there may be a correlation of flatness and luminosity with redshift.

Particle acceleration at plasma shocks appears to be ubiquitous in the universe, spanning systems in the heliosphere, supernova remnants, and relativistic jets in distant active galaxies and gamma-ray bursts. This review addresses some of the key issues for shock acceleration theory that require resolution in order to propel our understanding of particle energization in astrophysical environments. These include magnetic field amplification in shock ramps, the non-linear hydrodynamic interplay between thermal ions and their extremely energetic counterparts possessing ultrarelativistic energies, and the ability to inject and accelerate electrons in both non-relativistic and relativistic shocks. Recent observational developments that impact these issues are summarized. While these topics are currently being probed by astrophysicists using numerical simulations, they are also ripe for investigation in laboratory experiments, which potentially can provide valuable insights into the physics of cosmic shocks.

Measuring the polarization of the prompt gamma-ray emission from GRBs can significantly improve our understanding of both the GRB emission mechanisms, as well as the underlying engine driving the explosion. We searched for polarization in the prompt gamma-ray emission of GRB 041219a with the SPI instrument on INTEGRAL. Using multiple-detector coincidence events in the 100--350 keV energy band, our analysis yields a polarization fraction from this GRB of 99 +- 33 %. Statistically, we cannot claim a polarization detection from this source. Moreover, different event selection criteria lead to even less significant polarization fractions, e.g. lower polarization fractions are obtained when higher energies are included in the analysis. We cannot strongly rule out the possibility that the measured modulation is dominated by instrumental systematics. Therefore, SPI observations of GRB 041219a do not significantly constrain GRB models. However, this measurement demonstrates the capability of SPI to measure polarization, and the techniques developed for this analysis.

A significant new development in the study of Anomalous X-ray Pulsars (AXPs) has been the recent discovery by INTEGRAL and RXTE of flat, hard X-ray components in three AXPs. These non-thermal spectral components differ dramatically from the steeper quasi-power-law tails seen in the classic X-ray band in these sources. A prime candidate mechanism for generating this new component is resonant, magnetic Compton upscattering. This process is very efficient in the strong magnetic fields present in AXPs. Here an introductory exploration of an inner magnetospheric model for upscattering of surface thermal X-rays in AXPs is offered, preparing the way for an investigation of whether such resonant upscattering can explain the 20-150 keV spectra seen by INTEGRAL. Characteristically flat emission spectra produced by non-thermal electrons injected in the emission region are computed using collision integrals. A relativistic QED scattering cross section is employed so that Klein-Nishina reductions are influential in determining the photon spectra and fluxes. Spectral results depend strongly on the magnetospheric locale of the scattering and the observer's orientation, which couple directly to the angular distributions of photons sampled.

Recent advances in the understanding of the properties of supernova remnant shocks have been precipitated by the Chandra and XMM X-ray Observatories, and the HESS Atmospheric Cerenkov Telescope in the TeV band. A critical problem for this field is the understanding of the relative degree of dissipative heating/energization of electrons and ions in the shock layer. This impacts the interpretation of X-ray observations, and moreover influences the efficiency of injection into the acceleration process, which in turn feeds back into the thermal shock layer energetics and dynamics. This paper outlines the first stages of our exploration of the role of charge separation potentials in non-relativistic electron-ion shocks where the inertial gyro-scales are widely disparate, using results from a Monte Carlo simulation. Charge density spatial profiles were obtained in the linear regime, sampling the inertial scales for both ions and electrons, for different magnetic field obliquities. These were readily integrated to acquire electric field profiles in the absence of self-consistent, spatial readjustments between the electrons and the ions. It was found that while diffusion plays little role in modulating the linear field structure in highly oblique and perpendicular shocks, in quasi-parallel shocks, where charge separations induced by gyrations are small, and shock-layer electric fields are predominantly generated on diffusive scales.

The Advanced Compton Telescope (ACT), the next major step in gamma-ray astronomy, will probe the fires where chemical elements are formed by enabling high-resolution spectroscopy of nuclear emission from supernova explosions. During the past two years, our collaboration has been undertaking a NASA mission concept study for ACT. This study was designed to (1) transform the key scientific objectives into specific instrument requirements, (2) to identify the most promising technologies to meet those requirements, and (3) to design a viable mission concept for this instrument. We present the results of this study, including scientific goals and expected performance, mission design, and technology recommendations.

The spectra obtained above 100 MeV by the EGRET experiment aboard the Compton Gamma-Ray Observatory for a handful of gamma-ray bursts has given no indication of any spectral attenuation that might preclude detection of bursts at higher energies. With the discovery of optical afterglows and counterparts to bursts in the last few years, enabling the determination of significant redshifts for these sources, it is anticipated that profound spectral attenuation will arise in the Gamma-Ray Large Area Space Telescope (GLAST) energy band of 30 MeV-300 GeV for many if not most bursts. This paper explores time-dependent expectations for burst spectral properties in the EGRET/GLAST band, focusing on how attenuation of photons by pair creation internal to the source generates distinctive spectral signatures. The energy of spectral breaks and the associated spectral indices provide valuable information that constrains the bulk Lorentz factor of the GRB outflow at a given time. Moreover, the distinct temporal behavior that is present for internal attenuation is easily distinguished from extrinsic absorption due to intervening cosmic background fields. These characteristics define palpable observational goals for both spaced-based hard gamma-ray experiments such as GLAST, and ground-based Cherenkov telescopes, and strongly impact the observability of bursts above 300 MeV.

Heliospheric shocks are excellent systems for testing theories of particle acceleration in their environs. These generally fall into two classes: (1) interplanetary shocks that are linear in their ion acceleration characteristics, with the non-thermal ions serving as test particles, and (2) non-linear systems such as the Earth's bow shock and the solar wind termination shock, where the accelerated ions strongly influence the magnetohydrodynamic structure of the shock. This paper explores the modelling of diffusive acceleration at a particular interplanetary shock, with an emphasis on explaining in situ measurements of ion distribution functions. The observational data for this event was acquired on day 292 of 1991 by the Ulysses mission. The modeling is performed using a well-known kinetic Monte Carlo simulation, which has yielded good agreement with observations at several heliospheric shocks, as have other theoretical techniques, namely hybrid plasma simulations, and numerical solution of the diffusion-convection equation. In this theory/data comparison, it is demonstrated that diffusive acceleration theory can, to first order, successfully account for both the proton distribution data near the shock, and the observation of energetic protons farther upstream of this interplanetary shock than lower energy pick-up protons, using a single turbulence parameter. The principal conclusion is that diffusive acceleration of inflowing upstream ions can model this pick-up ion-rich event without the invoking any seed pre-acceleration mechanism, though this investigation does not rule out the action of such pre-acceleration.

The acceleration of interstellar pick-up ions as well as solar wind species has been observed at a multitude of interplanetary (IP) shocks by different spacecraft. The efficiency of injection of the pick-up ion component differs from that of the solar wind, and is expected to be strongly enhanced at highly oblique and quasi-perpendicular shock events, in accord with inferences from {\it in situ} observations. This paper explores theoretical modeling of the phase space distributions of accelerated ions obtained by the Ulysses mission for the Day 292, 1991 shock associated with a corotating interaction region, encountered before Ulysses' fly-by of Jupiter. A Monte Carlo simulation is used to model the acceleration process, adapting a technique that has been successfully tested on earlier IP shocks possessing minimal pick-up ion presence. Phase space distributions from the simulation technique for various low mass ions are compared with SWICS and HI-SCALE data to deduce values of a ``turbulence parameter'' that controls the efficiency of injection, and the degree of cross-field diffusion. Acceptable fits are obtained for the $H^+$ and $He^+$ populations using standard prescriptions for the pick-up ion distribution; $He^{++}$ spectral data was only fit well for scenarios very close to the Bohm diffusion limit. It is also found that the simulation successfully accounts for the observation of energetic protons farther upstream of the forward shock than lower energy pick-up protons, using the same turbulence parameter that is required to achieve reasonable spectral fits.

Cyclotron decay and absorption rates have been well studied in the literature, focusing primarily on spectral, angular and polarization dependence. Astrophysical applications usually do not require retention of information on the electron spin state, and these are normally averaged in obtaining the requisite rates. In magnetic fields, higher order quantum processes such as Compton scattering become resonant at the cyclotron frequency and its harmonics, with the resonances being formally divergent. Such divergences are usually eliminated by accounting for the finite lifetimes of excited Landau states. This practice requires the use of spin-dependent cyclotron rates in order to obtain accurate determinations of process rates very near cyclotronic resonances, the phase space domain most relevant for certain applications to pulsar models. This paper develops previous results in the literature to obtain compact analytic expressions for cyclotron decay rates/widths in terms of a series of Legendre functions of the second kind; these expressions can be expediently used in astrophysical models. The rates are derived using two popular eigenstate formalisms, namely that due to Sokolov and Ternov, and that due to Johnson and Lippmann. These constitute two sets of eigenfunctions of the Dirac equation that diagonalize different operators, and accordingly yield different spin-dependent cyclotron rates. This paper illustrates the attractive Lorentz transformation characteristics of the Sokolov and Ternov formulation, which is another reason why it is preferable when electron spin information must be explicitly retained.

The remnant G347.3-0.5 exhibits strong shell emission in the radio and X-ray bands, and has a purported detection in the TeV gamma-ray band by the CANGAROO-II telescope. The CANGAROO results were touted as evidence for the production of cosmic ray ions, a claim that has proven controversial due to constraining fluxes associated with a proximate unidentified EGRET source 3EG J1714-3857. HESS has now seen this source in the TeV band. The complex environment of the remnant renders modeling of its broadband spectrum sensitive to assumptions concerning the nature and parameters of the circumremnant medium. This paper explores a sampling of reasonable possibilities for multiwavelength spectral predictions from this source, using a non-linear model of diffusive particle acceleration at the shocked shell. The magnetic field strength, shell size and degree of particle cross-field diffusion act as variables to which the radio to X-ray to gamma-ray signal is sensitive. The modeling of the extant data constrains these variables, and the potential impact of the recent HESS detection on such parameters is addressed. Putative pion decay signals in hard gamma-rays resulting from hadronic interactions in dense molecular clouds are briefly discussed; the requisite suppression of the GeV component needed to accommodate the 3EG J1714-3857 EGRET data provides potential bounds on the diffusive distance from the shell to the proximate clouds.

The principal paradigm for gamma-ray bursts suggest that the prompt transient gamma-ray signal arises from multiple shocks internal to the relativistic expansion. This paper illustrates some properties of diffusive acceleration at relativistic shocks that pertain to GRB models, providing interpretation of the BATSE/EGRET data. Using a standard Monte Carlo simulation, computations of the spectral shape, and the range of spectral indices are presented, as functions of the shock speed, magnetic field obliquity and the type of particle scattering. It is clear that while quasi-parallel, relativistic shocks with fields approximately normal to the shock plane can efficiently accelerate particles, highly oblique and perpendicular ones cannot unless the particle diffusion is almost isotropic, i.e. extremely close to the Bohm limit. Accordingly, an array of distribution indices should be present in the burst population, as is exhibited by the EGRET data, and even mildly-relativistic internal shocks require strong field turbulence in order to model the >100 MeV observations. In addition, recent spectral fitting to burst data in the BATSE band is discussed, providing probes of the efficiency of injection of non-thermal electrons.

The process of diffusive acceleration of charged particles in shocked plasmas is widely invoked in astrophysics to account for the ubiquitous presence of signatures of non-thermal relativistic electrons and ions in the universe. A key characteristic of this statistical energization mechanism is the absence of a momentum scale; astrophysical systems generally only impose scales at the injection (low energy) and loss (high energy) ends of the particle spectrum. The existence of structure in the cosmic ray spectrum (the "knee") at around 3000 TeV has promoted contentions that there are at least two origins for cosmic rays, a galactic one supplying those up to the knee, and even beyond, and perhaps an extragalactic one that can explain even the ultra-high energy cosmic rays (UHECRs) seen at 1-300 EeV. Accounting for the UHECRs with familiar astrophysical sites of acceleration has historically proven difficult due to the need to assume high magnetic fields in order to reduce the shortest diffusive acceleration timescale, the ion gyroperiod, to meaningful values. Yet active galaxies and gamma-ray bursts remain strong and interesting candidate sources for UHECRs, turning the theoretical focus to relativistic shocks. This review summarizes properties of diffusive shock acceleration that are salient to the issue of UHECR generation. These include spectral indices, acceleration efficencies and timescales, as functions of the shock speed and mean field orientation, and also the nature of the field turbulence. The interpretation of these characteristics in the context of gamma-ray burst models for the production of UHECRs is also examined.

The principal paradigm for the prompt emission of gamma-ray bursts invokes diffusive acceleration at shocks internal to the dynamic ultrarelativistic outflow. This paper explores expectations for burst emission spectra arising from shock acceleration theory in the limit of particles cooling much slower than their acceleration. Parametric fits to burst spectra obtained by the Compton Gamma-Ray Observatory (CGRO) are explored for the cases of the synchrotron, inverse Compton and synchrotron self-Compton radiation mechanisms, using a linear combination of thermal and non-thermal electron populations. These fits demand that the preponderance of electrons that are responsible for the prompt emission reside in an intrinsically non-thermal population, strongly contrasting particle distributions obtained from acceleration simulations. This implies a potential conflict for acceleration scenarios where the non-thermal electrons are drawn directly from a thermal gas. It is also found that the CGRO data preclude effective spectroscopic discrimination between the synchrotron and inverse Compton mechanisms. This situation may be resolved with future missions probing gamma-ray bursts, namely Swift and GLAST. However, the synchrotron self-Compton (SSC) spectrum is characteristically too broad near the $\nu F_{\nu}$ peak to viably account for bursts such as GRB 910601, GRB 910814 and GRB 990123. It is concluded that the SSC mechanism may be generally incompatible with differential burst spectra steeper than around $E^{-2.5}$ above the peak, unless the synchrotron component is strongly self-absorbed.

A principal candidate for quiescent non-thermal gamma-ray emission from magnetars is resonant inverse Compton scattering in the strong fields of their magnetospheres. This paper outlines expectations for such emission, formed from non-thermal electrons accelerated in a pulsar-like polar cap potential upscattering thermal X-rays from the hot stellar surface. The resultant spectra are found to be strikingly flat, with fluxes and strong pulsation that could be detectable by GLAST.

The study of pulsars in the three and a half decades since their discovery has highlighted a handful of issues critical to their understanding. To date there is no consensus on the physical mechanism for their radio radio emission, despite a rapid increase in the observed population due to the Parkes Multi-Beam survey and prospects for similar growth in the radio population database in the near future. The small subset of pulsars that emit at X-ray to gamma-ray wavelengths are critical to refining the pulsar paradigm since this energy band (i) is where the vast majority of radiative luminosity is observed, and (ii) is intimately connected to the pair winds that form the dominant mode of energy deposition in the circum-pulsar environment. The most crucial point of contention pertaining to the high energy astrophysics of pulsars is the location of the acceleration region in their magnetospheres: is an outer gap model or a polar cap scenario (or both) the most appropriate picture. Radiative signatures provide the clues to this current enigma. This review focuses on salient characteristics of the polar cap scenario; these form the basis for discriminating observational diagnostics that should drive pulsars studies in the GLAST era just three years away.

Shock jump conditions are obtained for steady-state, plane shocks with oblique magnetic fields and arbitrary flow speeds. For ultrarelativistic and nonrelativistic shocks, the jump conditions may be solved analytically. For mildly relativistic shocks, analytic solutions are obtained for isotropic pressure using an approximation for the adiabatic index that is valid in high sonic Mach number cases. In the more general case of gyrotropic pressure, the jump conditions cannot be solved analytically without additional assumptions, and the effects of gyrotropic pressure are investigated by parameterizing the distribution of pressure parallel and perpendicular to the magnetic field. Our numerical solutions reveal that relatively small departures from isotropy (e.g., about 20%) produce significant changes in the shock compression ratio, r, at all shock Lorentz factors, including ultrarelativistic ones, where an analytic solution with gyrotropic pressure is obtained. In particular, either dynamically important fields or significant pressure anisotropies can incur marked departures from the canonical gas dynamic value of r=3 for a shocked ultrarelativistic flow and this may impact models of particle acceleration in gamma-ray bursts and other environments where relativistic shocks are inferred. The jump conditions presented apply directly to test-particle acceleration, and will facilitate future self-consistent numerical modeling of particle acceleration at oblique, relativistic shocks.

The advent of the next generation of gamma-ray experiments, led by GLAST, AGILE, INTEGRAL and a host of atmospheric \v{C}erenkov telescopes coming on line in the next few years, will enable ground-breaking discoveries relating to the presently enigmatic set of EGRET/CGRO UID galactic sources that have yet to find definitive identifications. Pulsars are principal candidates for such sources, and many are expected to be detected by GLAST, some that are radio-selected, like most of the present EGRET/Comptel pulsars, and perhaps even more that are detected via independent pulsation searches. At this juncture, it is salient to outline the principal predictions of pulsar models that might aid identification of gamma-ray sources, and moreover propel subsequent interpretation of their properties. This review summarizes relevant characteristics of the polar cap model, emphasizing where possible distinctions from the competing outer gap model. Foremost among these considerations are the hard X-ray to gamma-ray spectral shape, high energy cutoffs and pulse profiles, and how these characteristics generally depend on pulsar period and period derivative, as well as observational viewing angle. The polar cap model exhibits definitive signatures that will be readily tested by the detections of GLAST and other experiments, thereby establishing cogent observational diagnostics. The paper focuses on different classes of pulsars that might define agendas and parameter regimes for blind gamma-ray pulsation searches; examples include the highly-magnetized ones that are currently quite topical in astrophysics.

The absence of radio pulsars with long periods has lead to the popular notion of a high P ``death line.'' In the standard picture, beyond this boundary, pulsars with low spin rates cannot accelerate particles above the stellar surface to high enough energies to initiate pair cascades, and the pair creation needed for radio emission is strongly suppressed. In this paper we explore the possibility of another pulsar ``death line'' in the context of polar cap models, corresponding to high magnetic fields B in the upper portion of the period-period derivative diagram, a domain where few radio pulsars are observed. The origin of this high B boundary, which may occur when B becomes comparable to or exceeds $B_{\rm cr} = 4.4 \times 10^{13}$ Gauss, is also due to the suppression of magnetic pair creation, but primarily because of ineffective competition with magnetic photon splitting. Threshold pair creation also plays a prominent role in the suppression of cascades. We present Monte Carlo calculations of the pair yields in photon splitting/pair cascades which show that, in the absence of scattering effects, pair production is effectively suppressed, but only if all three modes of photon splitting allowed by QED are operating in high fields. This paper describes the probable shape and position of the new ``death line,'' above which pulsars are expected to be radio quiet, but perhaps still X-ray and gamma-ray bright. The hypothesized existence of radio-quiet sources finds dramatic support in the recent discovery of ultra-strong fields in Soft Gamma-ray Repeaters and Anomalous X-ray Pulsars. Guidelines for moderate to high B pulsar searches at radio wavelengths and also in the soft and hard gamma-ray bands are presented.

This paper explores the effects of strong magnetic fields on the Compton scattering of relativistic electrons. Recent studies of upscattering and energy loss by relativistic electrons that have used the non-relativistic, magnetic Thomson cross section for resonant scattering or the Klein-Nishina cross section for non-resonant scattering do not account for the relativistic quantum effects of strong fields ($ > 4 \times 10^{12}$ G). We have derived a simplified expression for the exact QED scattering cross section for the broadly-applicable case where relativistic electrons move along the magnetic field. To facilitate applications to astrophysical models, we have also developed compact approximate expressions for both the differential and total polarization-dependent cross sections, with the latter representing well the exact total QED cross section even at the high fields believed to be present in environments near the stellar surfaces of Soft Gamma-Ray Repeaters and Anomalous X-Ray Pulsars. We find that strong magnetic fields significantly lower the Compton scattering cross section below and at the resonance, when the incident photon energy exceeds $m_ec^2$ in the electron rest frame. The cross section is strongly dependent on the polarization of the final scattered photon. Below the cyclotron fundamental, mostly photons of perpendicular polarization are produced in scatterings, a situation that also arises above this resonance for sub-critical fields. However, an interesting discovery is that for super-critical fields, a preponderance of photons of parallel polarization results from scatterings above the cyclotron fundamental. This characteristic is both a relativistic and magnetic effect not present in the Thomson or Klein-Nishina limits.

Calculations of reaction rates for the third-order QED process of photon splitting in strong magnetic fields traditionally have employed either the effective Lagrangian method or variants of Schwinger's proper-time technique. Recently, Mentzel, Berg and Wunner (1994) presented an alternative derivation via an S-matrix formulation in the Landau representation. Advantages of such a formulation include the ability to compute rates near pair resonances above pair threshold. This paper presents new developments of the Landau representation formalism as applied to photon splitting, providing significant advances beyond the work of Mentzel et al. by summing over the spin quantum numbers of the electron propagators, and analytically integrating over the component of momentum of the intermediate states that is parallel to field. The ensuing tractable expressions for the scattering amplitudes are satisfyingly compact, and of an appearance familiar to S-matrix theory applications. Such developments can facilitate numerical computations of splitting considerably both below and above pair threshold. Specializations to two regimes of interest are obtained, namely the limit of highly supercritical fields and the domain where photon energies are far inferior to that for the threshold of single-photon pair creation. In particular, for the first time the low-frequency amplitudes are simply expressed in terms of the Gamma function, its integral and its derivatives. In addition, the equivalence of the asymptotic forms in these two domains to extant results from effective Lagrangian/proper-time formulations is demonstrated.

We have extended a simple model of nonlinear diffusive shock acceleration (Berezhko and Ellison 1999) to include the injection and acceleration of electrons and the production of photons from bremsstrahlung, synchrotron, inverse-Compton, and pion-decay processes. We argue that the results of this model, which is simpler to use than more elaborate ones, offer a significant improvement over test-particle, power-law spectra which are often used in astrophysical applications of diffusive shock acceleration. With an evolutionary supernova remnant (SNR) model to obtain shock parameters as functions of ambient interstellar medium parameters and time, we predict broad-band continuum photon emission from supernova remnants in general, and SN1006 in particular, showing that our results compare well with the more complete time-dependent and spherically symmetric nonlinear model of Berezhko, Ksenofontov, and Petukhov (1999). We discuss the implications nonlinear shock acceleration has for X-ray line emission, and use our model to describe how ambient conditions determine the TeV/radio flux ratio, an important parameter for gamma-ray observations of radio SNRs.

This paper summarizes highlights of the OG3.1, 3.2 and 3.3 sessions of the XXVIth International Cosmic Ray Conference in Salt Lake City, which were devoted to issues of origin/composition, acceleration and propagation.

The observation by the CANGAROO experiment of TeV emission from SN 1006, in conjunction with several instances of non-thermal X-ray emission from supernova remnants, has led to inferences of super-TeV electrons in these extended sources. While this is sufficient to propel the theoretical community in their modelling of particle acceleration and associated radiation, the anticipated emergence in the next decade of a number of new experiments probing the TeV and sub-TeV bands provides further substantial motivation for modellers. In particular, the quest for obtaining unambiguous gamma-ray signatures of cosmic ray ion acceleration defines a ``Holy Grail'' for observers and theorists alike. This review summarizes theoretical developments in the prediction of MeV-TeV gamma-rays from supernova remnants over the last five years, focusing on how global properties of models can impact, and be impacted by, hard gamma-ray observational programs, thereby probing the supernova remnant environment. Properties of central consideration include the maximum energy of accelerated particles, the density of the unshocked interstellar medium, the ambient magnetic field, and the relativistic electron-to-proton ratio. Criteria for determining good candidate remnants for observability in the TeV band are identified.

One of the principal results obtained by the EGRET experiment aboard the Compton Gamma-Ray Observatory (CGRO) was the detection of several Gamma-ray bursts (GRBs) above 100 MeV. The broad-band spectra obtained for these bursts gave no indication of any high energy spectral attenuation that might preclude detection of bursts by ground-based Cerenkov telescopes (ACTs), thus motivating several TeV observational programs. This paper explores the expectations for the spectral properties in the TeV and sub-TeV bands for bursts, in particular how attenuation of photons by pair creation internal to the source modifies the spectrum to produce distinctive spectral signatures. The energy of spectral breaks and the associated spectral indices provide valuable information that can constrain the bulk Lorentz factor of the GRB outflow at a given time. These characteristics define palpable observational goals for ACT programs, and strongly impact the observability of bursts in the TeV band.

Most recent extragalactic models of gamma-ray bursts consider the expansion of a relativistic blast wave, emanating from a solar-mass type progenitor, into the surrounding interstellar medium as the site for their activity. The popular perception is that the optical afterglows result from the external shock interface, while the prompt transient gamma-ray signal arises from multiple shocks internal to the expansion. This paper illustrates a number of acceleration properties of relativistic and ultrarelativistic shocks that pertain to GRB models, by way of a standard Monte Carlo simulation. Computations of the spectral shape, the range of spectral indices, and the energy gain per shock crossing are presented, as functions of the shock speed and the type of particle scattering.

Developments over the last couple of years have supported the interpretation that anomalous X-ray pulsars (AXPs) and soft gamma repeaters (SGRs) possess unusually high magnetic fields, and furthermore may represent a class or classes of neutron stars distinct from the population of conventional radio pulsars. We have recently suggested that such a dichotomization of the pulsar population may naturally arise due to the inherently different conditions in subcritical and supercritical magnetic fields. In this paper, we summarize, within the polar gap model, expectations for observable properties of highly magnetized pulsars, conventional or anomalous. This includes a discussion of the potential suppression of pair production and cascade generation in very strong fields by photon splitting and by threshold pair creation, which might explain radio quiescence in AXPs and SGRs. X-ray and hard gamma-ray spectral properties and trends are identified, with a view to establishing goals for future high energy experimental programs.

We have modeled the injection and acceleration of pickup ions at the solar wind termination shock and investigated the parameters needed to produce the observed Anomalous Cosmic Ray (ACR) fluxes. A non-linear Monte Carlo technique was employed which in effect solves the Boltzmann equation and is not restricted to near isotropic particle distribution functions. This technique models the injection of thermal and pickup ions, the acceleration of these ions, and the determination of the shock structure under the influence of the accelerated ions. The essential effects of injection are treated in a mostly self-consistent manner, including effects from shock obliquity, cross-field diffusion, and pitch-angle scattering. Using recent determinations of pickup ion densities, we are able to match the absolute flux of hydrogen in the ACRs by assuming that pickup ion scattering mean free paths, at the termination shock, are much less than an AU and that moderately strong cross-field diffusion occurs. Simultaneously, we match the flux ratios He$^{+}$/H$^+$ or O$^{+}$/H$^+$ to within a factor ~ 3. If the conditions of strong scattering apply, no pre-termination-shock injection phase is required and the injection and acceleration of pickup ions at the termination shock is totally analogous to the injection and acceleration of ions at highly oblique interplanetary shocks recently observed by the Ulysses spacecraft.

We calculate particle spectra and continuum photon emission from the Cassiopeia A supernova remnant (SNR). The particle spectra, ion and electron, result from diffusive shock acceleration at the forward SNR shock and are determined with a nonlinear Monte Carlo calculation. The calculation self-consistently determines the shock structure under the influence of ion pressure, and includes a simple parameterized treatment of electron injection and acceleration. Our results are compared to photon observations, concentrating on the connection between the Radio and GeV-TeV gamma-ray range, and to cosmic ray ion observations. We include new upper limits from the Cherenkov Array at Themis (CAT) imaging Cherenkov telescope and the Whipple 10m gamma-ray telescope at > 400 GeV. These new limits support the suggestion (e.g. Cowsik & Sarkar 1980; Allen et. al. 1997) that energetic electrons are emitting synchrotron radiation in an extremely high magnetic field (~ 1000 microGauss), far greater than values routinely assigned to the ISM, and help to constrain our model. The large magnetic field allows acceleration of cosmic ray ions to well above $10^{15}$ eV per nucleon in the ~ 300 yr lifetime of Cas A.

There has recently been interest in the role of inverse bremsstrahlung, the emission of photons by fast suprathermal ions in collisions with ambient electrons possessing relatively low velocities, in tenuous plasmas in various astrophysical contexts. This follows a long hiatus in the application of suprathermal ion bremsstrahlung to astrophysical models since the early 1970s. The potential importance of inverse bremsstrahlung relative to normal bremsstrahlung, i.e. where ions are at rest, hinges upon the underlying velocity distributions of the interacting species. In this paper, we identify the conditions under which the inverse bremsstrahlung emissivity is significant relative to that for normal bremsstrahlung in shocked astrophysical plasmas. We determine that, since both observational and theoretical evidence favors electron temperatures almost comparable to, and certainly not very deficient relative to proton temperatures in shocked plasmas, these environments generally render inverse bremsstrahlung at best a minor contributor to the overall emission. Hence inverse bremsstrahlung can be safely neglected in most models invoking shock acceleration in discrete sources such as supernova remnants. However, on scales > 100pc distant from these sources, Coulomb collisional losses can deplete the cosmic ray electrons, rendering inverse bremsstrahlung, and perhaps bremsstrahlung from knock-on electrons, possibly detectable.

The complete absence of radio pulsars with periods exceeding a few seconds has lead to the popular notion of the existence of a high period death line. We have recently postulated the existence of another radio quiescence boundary at high magnetic fields ($B\gtrsim 4\times 10^{13}$G) in the upper portion of the period-period derivative diagram, a domain where no radio pulsars are observed. The origin of this high B boundary is also due to the suppression of magnetic pair creation, $\gamma\to e^{\pm}$, but mainly because of competition with the exotic QED process of magnetic photon splitting, $\gamma\to\gamma\gamma$, coupled with ground state pair creation. This mechanism could also explain the low spectral cutoff energy of the gamma-ray pulsar PSR1509-58, which lies near the high B death-line. In this paper, we summarize the hypothesis of this new ``death line,'' and discuss some subtleties of pair suppression that relate to photon polarization and positronium formation. We identify several ways in which Integral will serve as a useful tool in probing high field pulsars, given that the radio quiescence boundary is transparent to the high energy bands.

Supernova remnants (SNRs) are widely believed to be the principal source of galactic cosmic rays. Such energetic particles can produce gamma-rays and lower energy photons via interactions with the ambient plasma. In this paper, we present results from a Monte Carlo simulation of non-linear shock structure and acceleration coupled with photon emission in shell-like SNRs. These non-linearities are a by-product of the dynamical influence of the accelerated cosmic rays on the shocked plasma and result in distributions of cosmic rays which deviate from pure power-laws. Such deviations are crucial to acceleration efficiency and spectral considerations, producing GeV/TeV intensity ratios that are quite different from test particle predictions. The Sedov scaling solution for SNR expansions is used to estimate important shock parameters for input into the Monte Carlo simulation. We calculate ion and electron distributions that spawn neutral pion decay, bremsstrahlung, inverse Compton, and synchrotron emission, yielding complete photon spectra from radio frequencies to gamma-ray energies. The cessation of acceleration caused by the spatial and temporal limitations of the expanding SNR shell in moderately dense interstellar regions can yield spectral cutoffs in the TeV energy range; these are consistent with Whipple's TeV upper limits on unidentified EGRET sources. Supernova remnants in lower density environments generate higher energy cosmic rays that produce predominantly inverse Compton emission at super-TeV energies; such sources will generally be gamma-ray dim at GeV energies.

We have modeled the injection and acceleration of pickup ions at the solar wind termination shock and investigated the parameters needed to produce the observed Anomalous Cosmic Ray (ACR) fluxes. A non-linear Monte Carlo technique was employed, which in effect solves the Boltzmann equation and is not restricted to near-isotropic particle distribution functions. This technique models the injection of thermal and pickup ions, the acceleration of these ions, and the determination of the shock structure under the influence of the accelerated ions. The essential effects of injection are treated in a mostly self-consistent manner, including effects from shock obliquity, cross-field diffusion, and pitch-angle scattering. Using recent determinations of pickup ion densities, we are able to match the absolute flux of hydrogen in the ACRs by assuming that pickup ion scattering mean free paths, at the termination shock, are much less than an AU and that modestly strong cross-field diffusion occurs. Simultaneously, we match the flux ratios He$^{+}$/H$^{+}$ or O$^{+}$/H$^{+}$ to within a factor $\sim 5$. If the conditions of strong scattering apply, no pre-termination-shock injection phase is required and the injection and acceleration of pickup ions at the termination shock is totally analogous to the injection and acceleration of ions at highly oblique interplanetary shocks recently observed by the Ulysses spacecraft. The fact that ACR fluxes can be modeled with standard shock assumptions suggests that the much-discussed "injection problem" for highly oblique shocks stems from incomplete (either mathematical or computer) modeling of these shocks rather than from any actual difficulty shocks may have in injecting and accelerating thermal or quasi-thermal particles.

The notable absence of radio pulsars having measured magnetic dipole surface field strengths above $B_0\sim 3\times 10^{13}$ Gauss naturally raises the question of whether this forms an upper limit to pulsar magnetization. Recently there has been increasing evidence that neutron stars possessing higher dipole spin-down fields do in fact exist, including a growing list of anomalous X-ray pulsars (AXPs) with long periods and spinning down with high period derivatives, implying surface fields of $10^{14}$--$10^{15}$ Gauss. Furthermore, the recently reported X-ray period and period derivative for the Soft Gamma-ray Repeater (SGR) source SGR1806-20 suggest a surface field around $10^{15}$ Gauss. None of these high-field pulsars have yet been detected as radio pulsars. We propose that high-field pulsars should be radio-quiet because electron-positron pair production in their magnetospheres, thought to be essential for radio emission, is efficiently suppressed in ultra-strong fields ($B_0\gtrsim 4\times 10^{13}$ Gauss) by the action of photon splitting, a quantum electrodynamical process in which a photon splits into two. Our computed radio quiescence boundary in the radio pulsar $P-\dot P$ diagram, where photon splitting overtakes pair creation, is located just above the boundary of the known radio pulsar population, neatly dividing them from the AXPs. We thus identify a physical mechanism that defines a new class of high-field radio-quiet neutron stars that should be detectable by their pulsed emission at X-ray and perhaps $\gamma$-ray energies.

We give a rigorous derivation of a theorem showing that charged particles in an arbitrary electromagnetic field with at least one ignorable spatial coordinate remain forever tied to a given magnetic-field line. Such a situation contrasts the significant motions normal to the magnetic field that are expected in most real three-dimensional systems. It is pointed out that, while the significance of the theorem has not been widely appreciated, it has important consequences for a number of problems and is of particular relevance for the acceleration of cosmic rays by shocks.

One of the principal results obtained by the Compton Gamma Ray Observatory relating to the study of gamma-ray bursts was the detection by the EGRET instrument of energetic ($>$100 MeV) photons from a handful of bright bursts. The most extreme of these was the single 18 GeV photon from the GRB940217 source. Given EGRET's sensitivity and limited field of view, the detection rate implies that such high energy emission may be ubiquitous in bursts. Hence expectations that bursts emit out to at least TeV energies are quite realistic, and the associated target-of-opportunity activity of the TeV gamma-ray community is well-founded. This review summarizes the observations and a handful of theoretical models for generating GeV--TeV emission in bursts sources, outlining possible ways that future positive detections could discriminate between different scenarios. The power of observations in the GeV--TeV range to distinguish between spectral structure intrinsic to bursts and that due to the intervening medium between source and observer is also discussed.

Diffusive shock acceleration in the environs of a remnant's expanding shell is a popular candidate for the origin of SNR gamma-rays. In this paper, results from our study of non-linear effects in shock acceleration theory and their impact on the gamma-ray spectra of SNRs are presented. These effects describe the dynamical influence of the accelerated cosmic rays on the shocked plasma at the same time as addressing how the non-uniformities in the fluid flow force the distribution of the cosmic rays to deviate from pure power-laws. Such deviations are crucial to gamma-ray spectral determination. Our self-consistent Monte Carlo approach to shock acceleration is used to predict ion and electron distributions that spawn neutral pion decay, bremsstrahlung and inverse Compton emission components for SNRs. We demonstrate how the spatial and temporal limitations imposed by the expanding SNR shell quench acceleration above critical energies in the 500 GeV - 10 TeV range, thereby spawning gamma-ray spectral cutoffs that are quite consistent with Whipple's TeV upper limits to the EGRET unidentified sources that have SNR associations. We also discuss the role of electron injection in shocks and its impact on the significance of electromagnetic components to GeV--TeV spectral formation.

Eleven bright gamma-ray bursts (GRBs) detected by BATSE have also been seen at much higher energies by EGRET, six at energies above 10 MeV. Such observations imply that these bursts are optically thin to photon-photon pair production at all observed energies. For bursts more than about 30pc away, internal transparency can be achieved only if the source is moving with a relativistic bulk Lorentz factor $\Gamma\gg 1$, or if the radiation is highly beamed. Early calculations of $\gamma\gamma\to e^+e^-$ considerations for GRBs were limited to cases of a beam with opening half-angle $\Thetab\sim 1/\Gamma$, or expansions of infinitely thin spherical shells. This paper presents our extension of pair production optical depth calculations in relativistically expanding sources to more general geometries, including shells of finite thickness and arbitrary opening angle. The problem is reduced analytically to a single integral in the broadly applicable case of observing photons along the axis of the expansion. We find that the minimum bulk Lorentz factor for the EGRET sources to be optically thin is only moderately dependent on the shell thickness and virtually independent of its opening solid angle if $\Thetab\gtrsim 1/\Gamma$. This insensitivity to $\Thetab$ relieves the commonly-perceived number problem for non-repeating sources at cosmological distances, i.e. it is not necessary to invoke small $\Thetab$ to effect photon escape.

Due to observations made by the Compton Gamma-Ray Observatory over the last six years, it appears that a number of galactic supernova remnants may be candidates for sources of cosmic gamma-rays. These include shell-type remnants such as IC443 and $\gamma$ Cygni, which have no known parent pulsars, but have significant associations with unidentified EGRET sources, and others that appear to be composite, where a pulsar is embedded in a shell (e.g. W44 and Vela), or are purely pulsar-driven, such as the Crab Nebula. This review discusses our present understanding of gamma-ray production in plerionic and non-plerionic supernova remnants, and explores the relationship between such emission and that in other wavebands. Focuses include models of the Crab and Vela nebulae, the composite nature of W44, the relationship of shell-type remnants to cosmic ray production, the relative importance of shock-accelerated protons and electrons, constraints on models placed by TeV, X-ray and radio observations, and the role of electrons injected directly into the remnants by parent pulsars.

The complete absence of radio pulsars with periods exceeding a few seconds has lead to the popular notion of the existence of a high $P$ death line. In the standard picture, beyond this boundary, pulsars with low spin rates cannot accelerate particles above the stellar surface to high enough energies to initiated pair cascades through curvature radiation, and the pair creation needed for radio emission is strongly suppressed. In this paper we postulate the existence of another pulsar ``death line,'' corresponding to high magnetic fields $B$ in the upper portion of the $\dot{P}$--$P$ diagram, a domain where few radio pulsars are observed. The origin of this high $B$ boundary, which occurs when $B$ becomes comparable to or exceeds $10^{13}$ Gauss, is again due to the suppression of magnetic pair creation $\gamma\to e^+e^-$, but in this instance, primarily because of ineffective competition with the exotic QED process of magnetic photon splitting. This paper describes the origin, shape and position of the new ``death line,'' above which pulsars are expected to be radio quiet, but perhaps still X-ray and $\gamma$-ray bright.

Supernova remnants are widely believed to be a principal source of galactic cosmic rays, produced by diffusive shock acceleration in the environs of the remnant's expanding shock. This review discusses recent modelling of how such energetic particles can produce gamma-rays via interactions with the remnants' ambient interstellar medium, specifically via neutral pion decay, bremsstrahlung and inverse Compton emission. Predictions that relate to the handful of associations between EGRET unidentified sources and known radio/optical/X-ray emitting remnants are summarized. The cessation of acceleration above 1 TeV - 10 TeV energies in young shell-type remnants is critical to model consistency with Whipple's TeV upper limits; these observations provide important diagnostics for theoretical models.

The mechanism of diffusive Fermi acceleration at collisionless plasma shock waves is widely invoked in astrophysics to explain the appearance of non-thermal particle populations in a variety of environments, including sites of cosmic ray production, and is observed to operate at several sites in the heliosphere. This review outlines the principal results from the theory of diffusive shock acceleration, focusing first on how it produces power-law distributions in test-particle regimes, where the shock dynamics are dominated by the thermal populations that provide the seed particles for the acceleration process. Then the importance of non-linear modifications to the shock hydrodynamics by the accelerated particles is addressed, emphasizing how these subsequently influence non-thermal spectral formation.

The splitting of photons in the presence of an intense magnetic field has recently found astrophysical applications in polar cap models of gamma-ray pulsars and in magnetar scenarios for soft gamma repeaters. Numerical computation of the polarization-dependent rates of this third order QED process for arbitrary field strengths and energies below pair creation threshold is difficult: thus early analyses focused on analytic developments and simpler asymptotic forms. The recent astrophysical interest spurred the use of the S-matrix approach by Mentzel, Berg and Wunner to determine splitting rates. In this paper, we present numerical computations of a full proper-time expression for the rate of splitting that was obtained by Stoneham, and is exact up to the pair creation threshold. While the numerical results derived here are in accord with the earlier asymptotic forms due to Adler, our computed rates still differ by as much as factors of 3 from the S-matrix re-evaluation of Wilke and Wunner, reflecting the extreme difficulty of generating accurate S-matrix numerics for fields below about \teq{4.4\times 10^{13}}Gauss. We find that our proper-time rates appear very accurate, and exceed Adler's asymptotic specializations significantly only for photon energies just below pair threshold and for supercritical fields, but always by less than a factor of around 2.6. We also provide a useful analytic series expansion for the scattering amplitude valid at low energies.

While the definitive detection of gamma-rays from known supernova remnants (SNRs) remains elusive, the collection of unidentified EGRET sources that may be associated with SNRs has motivated recent modelling of TeV emission from these sources. Current theoretical models use power-law shock-accelerated protons and electrons in their predictions of expected gamma-ray TeV fluxes from those unidentified EGRET sources with remnant associations. In this paper, we explore a more detailed non-linear shock acceleration model, which generates non-thermal proton distributions and includes a self-consistent determination of shock hydrodynamics. We obtain gamma-ray spectra for SNRs allowing for the cessation of acceleration to high energies that is due to the finite ages and sizes of remnants. Gamma-ray spectral cutoffs can be observed in the TeV range for reasonable remnant parameters, and deviations from power-law behaviour are found at all energies ranging from 1 MeV up to the cutoff. Correlated observations by INTEGRAL, Whipple and other instruments may provide stringent constraints to our understanding of supernova remnants.

Gamma-rays from pulsars can be efficiently attenuated in their magnetospheres via the mechanism of single-photon pair production and also the exotic QED process of photon splitting, which become prolific in fields approaching the quantum critical value of $B_{cr}=4.41\times 10^{13}$ Gauss. Recently we have published results of our modelling of strongly-magnetized $\gamma$-ray pulsars, which focused on the escape or attenuation of photons emitted near the pole at the neutron star surface in dipole fields, in a Schwarzschild metric. We found that pair production and splitting totally inhibit emission above around 10--30 MeV in PSR1509-58, whose surface field is inferred to be as high as $0.7B_{cr}$. Our model pulsar spectra are consistent with the EGRET upper limits for PSR1509-58 for a wide range of polar cap sizes. Here we review the principal predictions of our attenuation analysis, and identify how its powerful observational diagnostic capabilities relate to current and future gamma-ray experiments. Diagnostics include the energy of the gamma-ray turnover and the spectral polarization, which constrain the estimated polar cap size and field strength, and can determine the relative strength of splitting and pair creation.

Relativistic bulk motion with large Lorentz factors has recently been inferred for gamma-ray bursts regardless of whether they are of galactic or cosmological origin. This conclusion results from calculations of internal pair production transparency in bursts that usually assume an infinite power-law source spectrum for simplicity, an approximation that is quite adequate for some bursts detected by EGRET. However, for a given bulk Lorentz factor \teq{\Gamma}, photons above the EGRET range can potentially interact with sub-MeV photons in such calculations. Hence it is essential to accurately address the spectral curvature in bursts seen by BATSE. In this paper we present the major properties induced in photon-photon opacity considerations by such spectral curvature. The observed spectral breaks around 1 MeV turn out to be irrelevant to opacity in cosmological bursts, but are crucial to estimates of source transparency in the 1 GeV -- 1 TeV range for sources located in the galactic halo. We find that broad absorption troughs can arise at these energies for suitable bulk motion parameters \teq{\Gamma}. Such troughs are probably an unambiguous signature of a galactic halo population, and if observed by experiments such as Whipple, MILAGRO and GLAST, would provide powerful evidence that such bursts are not at cosmological distances.

The most stringent test of theoretical models of the first-order Fermi mechanism at collisionless astrophysical shocks is a comparison of the theoretical predictions with observational data on particle populations. Such comparisons have yielded good agreement between observations at the quasi-parallel portion of the Earth's bow shock and three theoretical approaches, including Monte Carlo kinetic simulations. This paper extends such model testing to the realm of oblique interplanetary shocks: here observations of proton and alpha particle distributions made by the SWICS ion mass spectrometer on Ulysses at nearby interplanetary shocks are compared with test particle Monte Carlo simulation predictions of accelerated populations. The plasma parameters used in the simulation are obtained from measurements of solar wind particles and the magnetic field upstream of individual shocks. Good agreement between downstream spectral measurements and the simulation predictions are obtained for two shocks by allowing the the ratio of the mean-free scattering length to the ionic gyroradius, to vary in an optimization of the fit to the data. Generally small values of this ratio are obtained, corresponding to the case of strong scattering. The acceleration process appears to be roughly independent of the mass or charge of the species.

We have developed a Monte Carlo technique for self-consistently calculating the hydrodynamic structure of oblique, steady-state shocks, together with the first-order Fermi acceleration process and associated non-thermal particle distributions. This is the first internally consistent treatment of modified shocks that includes cross-field diffusion of particles. Our method overcomes the injection problem faced by analytic descriptions of shock acceleration, and the lack of adequate dynamic range and artificial suppression of cross-field diffusion faced by plasma simulations; it currently provides the most broad and versatile description of collisionless shocks undergoing efficient particle acceleration. We present solutions for plasma quantities and particle distributions upstream and downstream of shocks, illustrating the strong differences observed between non-linear and test-particle cases. It is found that there are only marginal differences in the injection efficiency and resultant spectra for two extreme scattering modes, namely large-angle scattering and pitch-angle diffusion, for a wide range of shock parameters, i.e., for subluminal shocks with field obliquities less than or equal to 75 degrees and de Hoffmann-Teller frame speeds much less than the speed of light.

Magnetic photon splitting, a QED process that becomes important only in magnetic fields approaching the quantum critical value, B_cr = 4.41 X 10^13 Gauss, is investigated as a mechanism for attenuation of gamma-rays emitted near the surface of strongly-magnetized pulsars. We model photon splitting attenuation and subsequent splitting cascades in gamma-ray pulsars, including the dipole field and curved spacetime geometry of the neutron star magnetosphere. We focus specifically on PSR1509-58, which has the highest surface magnetic field of all the gamma-ray pulsars (B_0 = 3 X 10^13 Gauss). We find that splitting will not be important for most gamma-ray pulsars, i.e. those with B_0 <~ 0.2 B_cr, but will be important for gamma-ray pulsars having B_0 >~ 0.3 B_cr, where the splitting attenuation lengths and escape energies become comparable to or less than those for pair production. We compute Monte Carlo spectral models for PSR1509-58. We find that photon splitting, or combined splitting and pair production, can explain the unusually low cutoff energy (between 2 and 30 MeV) of PSR1509-58, and that the model cascade spectra, which display strong polarization, are consistent with the observed spectral points and upper limits for polar cap emission at a range of magnetic colatitudes up to ~ 25 degrees.

The spectra of soft gamma-ray repeaters (SGRs), with the exception of the March 5, 1979 main burst, are characterized by high-energy cutoffs around 30 keV and low-energy turnovers that are much steeper than a Wien spectrum. Baring (1995) found that the spectra of cascades due to photon splitting in a very strong, homogeneous magnetic field can soften spectra and produce good fits to the soft spectra of SGRs. Magnetic field strengths somewhat above the QED critical field strength $B_{\rm cr}$, where $B_{\rm cr} = 4.413 \times 10^{13}$ G, is required to produce cutoffs at 30-40 keV. We have improved upon this model by computing Monte Carlo photon splitting cascade spectra in a neutron star dipole magnetic field, including effects of curved space-time in a Schwarzschild metric. We investigate spectra produced by photons emitted at different locations and observer angles. We find that the general results of Baring hold for surface emission throughout most of the magnetosphere, but that emission in equatorial regions can best reproduce the constancy of SGR spectra observed from different bursts.

The modelling of gamma-ray burst (GRB) spectra has considerable potential for increasing the understanding of these enigmatic sources. A diversity of ideas and analyses has been generated over the last two decades to explain line features and continuum shapes, encompassing both older galactic neutron star and ``new age'' cosmological source models. This paper reviews some of the highlights of these studies, discussing the merits and limitations of various ideas, and in particular their compatibility with the observational data. The first focus will be on continuum models for GRBs, which include optically thin synchrotron emission and resonant Compton upscattering near galactic neutron stars, while the synchrotron and non-magnetic inverse Compton scattering mechanisms are prominent in the less well-developed cosmological scenarios. Line formation scenarios will then be discussed, in particular the scattering model for producing cyclotron features, which remains the only viable explanation for the Ginga observations of double lines. Absorption-like line production in cosmological burst models is generally difficult, though interesting notions such as femtolensing interference patterns have been proffered.

The exotic quantum process of photon splitting has great potential to explain the softness of emission in soft gamma repeaters (SGRs) if they originate in neutron stars with surface fields above the quantum critical field $B_{\rm cr}=4.413\times 10^{13}$Gauss. Splitting becomes prolific at such field strengths: its principal effect is to degrade photon energies, initiating a cascade that softens gamma-ray spectra. Uniform field cascade calculations have demonstrated that emission could be softened to the observed SGR energies for fields exceeding about $10^{14}$Gauss. Recently, we have determined splitting attenuation lengths and maximum energies for photon escape in neutron star environments including the effects of magnetospheric dipole field geometry. Such escape energies $\erg_{esc}$ suitably approximate the peak energy of the emergent spectrum, and in this paper we present results for $\erg_{esc}$ as a function of photon emission angles for polar cap and equatorial emission regions. The escape energy is extremely insensitive to viewing perspective for equatorial emission, arguing in favour of such a site for the origin of SGR activity.

The recent association of soft gamma repeaters (SGRs) with counterparts in other wavebands has sparked much interest in these sources. One of the recent models for these objects is that they originate in the environs of neutron stars with fields much stronger than the quantum critical field \teq{B_{cr}=4.413\times 10^{13}} Gauss. Near such neutron stars, dubbed magnetars, the exotic quantum process of magnetic photon splitting becomes prolific. Its principal effect is to degrade photon energies and thereby soften gamma-ray spectra from neutron stars; it has recently been suggested that splitting may be responsible for limiting the hardness of emission in SGRs, if these sources originate in neutron stars with supercritical surface fields. Seed photons in supercritical fields efficiently generate soft gamma-ray spectra, typical of repeaters. In this paper, the influence of the curved dipole field geometry of a neutron star magnetosphere on the photon splitting rate is investigated. The dependence of the attenuation length on the location and angular direction of the seed photons is explored.

The rate at which particles are accelerated by the first-order Fermi mechanism in shocks depends on the angle, \teq{\Tbone}, that the upstream magnetic field makes with the shock normal. The greater the obliquity the greater the rate, and in quasi-perpendicular shocks rates can be hundreds of times higher than those seen in parallel shocks. In many circumstances pertaining to evolving shocks (\eg, supernova blast waves and interplanetary traveling shocks), high acceleration rates imply high maximum particle energies and obliquity effects may have important astrophysical consequences. However, as is demonstrated here, the efficiency for injecting thermal particles into the acceleration mechanism also depends strongly on obliquity and, in general, varies inversely with \teq{\Tbone}. The degree of turbulence and the resulting cross-field diffusion strongly influences both injection efficiency and acceleration rates. The test particle \mc simulation of shock acceleration used here assumes large-angle scattering, computes particle orbits exactly in shocked, laminar, non-relativistic flows, and calculates the injection efficiency as a function of obliquity, Mach number, and degree of turbulence. We find that turbulence must be quite strong for high Mach number, highly oblique shocks to inject significant numbers of thermal particles and that only modest gains in acceleration rates can be expected for strong oblique shocks over parallel ones if the only source of seed particles is the thermal background.