We present a model for magnetic energy dissipation in a pulsar wind nebula. Better understanding of this process is required to assess the likelihood that certain astrophysical transients may be powered by the spin-down of a "millisecond magnetar." Examples include superluminous supernovae, gamma-ray bursts, and anticipated electromagnetic counterparts to gravitational wave detections of binary neutron star coalescence. Our model leverages recent progress in the theory of turbulent magnetic relaxation to specify a dissipative closure of the stationary magnetohydrodynamic (MHD) wind equations, yielding predictions of the magnetic energy dissipation rate throughout the nebula. Synchrotron losses are treated self-consistently. To demonstrate the model's efficacy, we show that it can reproduce many features of the Crab Nebula, including its expansion speed, radiative efficiency, peak photon energy, and mean magnetic field strength. Unlike ideal MHD models of the Crab (which lead to the so-called sigma problem) our model accounts for the transition from ultra to weakly magnetized plasma flow, and for the associated heating of relativistic electrons. We discuss how the predicted heating rates may be utilized to improve upon models of particle transport and acceleration in pulsar wind nebulae. We also discuss implications for the Crab Nebula's gamma-ray flares, and point out potential modifications to models of astrophysical transients invoking the spin-down of a millisecond magnetar.

One of the most intriguing results from the gamma-ray instruments in orbit has been the detection of powerful flares from the Crab Nebula. These flares challenge our understanding of pulsar wind nebulae and models for particle acceleration. We report on the portion of a multiwavelength campaign using Keck, HST, and Chandra concentrating on a small emitting region, the Crab's inner knot, located a fraction of an arcsecond from the pulsar. We find that the knot's radial size, tangential size, peak flux, and the ratio of the flux to that of the pulsar are correlated with the projected distance of the knot from the pulsar. A new approach, using singular value decomposition for analyzing time series of images, was introduced yielding results consistent with the more traditional methods while some uncertainties were substantially reduced. We exploit the characterization of the knot to discuss constraints on standard shock-model parameters that may be inferred from our observations assuming the inner knot lies near to the shocked surface. These include inferences as to wind magnetization, shock shape parameters such as incident angle and poloidal radius of curvature, as well as the IR/optical emitting particle enthalpy fraction. We find that while the standard shock model gives good agreement with observation in many respects, there remain two puzzles: (a) The observed angular size of the knot relative to the pulsar--knot separation is much smaller than expected; (b) The variable, yet high degree of polarization reported is difficult to reconcile with a highly relativistic downstream flow.

The afterglow emission from gamma-ray bursts (GRBs) is usually interpreted as synchrotron radiation from electrons accelerated at the GRB external shock, that propagates with relativistic velocities into the magnetized interstellar medium. By means of multi-dimensional particle-in-cell simulations, we investigate the acceleration performance of weakly magnetized relativistic shocks, in the magnetization range 0<sigma<1e-1. The pre-shock magnetic field is orthogonal to the flow, as generically expected for relativistic shocks. We find that relativistic perpendicular shocks propagating in electron-positron plasmas are efficient particle accelerators if the magnetization is sigma<1e-3. For electron-ion plasmas, the transition to efficient acceleration occurs for sigma<3e-5. Here, the acceleration process proceeds similarly for the two species, since the electrons enter the shock nearly in equipartition with the ions, as a result of strong pre-heating in the self-generated upstream turbulence. In both electron-positron and electron-ion shocks, we find that the maximum energy of the accelerated particles scales in time as t^(1/2). This scaling is shallower than the so-called (and commonly assumed) Bohm limit, and it naturally results from the small-scale nature of the Weibel turbulence generated in the shock layer. In magnetized plasmas, the energy of the accelerated particles increases until it reaches a saturation value that scales with the magnetization as sigma^(-1/4). Further energization is prevented by the fact that the self-generated turbulence is confined within a finite region of thickness proportional to sigma^(-1/2) around the shock. Our results can provide physically-grounded inputs for models of non-thermal emission from a variety of astrophysical sources, with particular relevance to GRB afterglows.

(shortened) We report the results of an investigation of particle acceleration and electron-positron plasma generation at low altitude in the polar magnetic flux tubes of Rotation Powered Pulsars, when the stellar surface is free to emit whatever charges and currents are demanded by the force-free magnetosphere. We observe novel behavior. a) When the current density is less than the Goldreich-Julian (GJ) value (0<j/j_{GJ}<1), space charge limited acceleration of the current carrying beam is mild, with the full GJ charge density being comprised of the charge density of the beam, co-existing with a cloud of electrically trapped particles with the same sign of charge as the beam. The voltage drops are on the order of mc^2/e, and pair creation is absent. b) When the current density exceeds the GJ value (j/j_{GJ}>1), the system develops high voltage drops, causing emission of gamma rays and intense bursts of pair creation. The bursts exhibit limit cycle behavior, with characteristic time scales somewhat longer than the relativistic fly-by time over distances comparable to the polar cap diameter (microseconds). c) In return current regions, where j/j_{GJ}<0, the system develops similar bursts of pair creation. In cases b) and c), the intermittently generated pairs allow the system to simultaneously carry the magnetospherically prescribed currents and adjust the charge density and average electric field to force-free conditions. We also elucidate the conditions for pair creating beam flow to be steady, finding that such steady flows can occupy only a small fraction of the current density parameter space of the force-free magnetospheric model. The generic polar flow dynamics and pair creation is strongly time dependent. The model has an essential difference from almost all previous quantitative studies, in that we sought the accelerating voltage as a function of the applied current.

I outline, from a theoretical and somewhat personal perspective, significant features of Pulsar Wind Nebulae as Cosmic Accelerators. I discuss recent studies of Pulsar Wind Nebulae (PWNe). I pay special attention to the recently discovered gamma ray flares in the Crab Nebula's emission, focusing on the possibility, raised by the observations, that the accelerating electric field exceeds the magnetic field, suggesting that reconnection in the persistent current layer (a current sheet) plays a significant role in the behavior of this well studied Pevatron. I address the present status of the termination shock model for the particle accelerator that converts the wind flow energy to the observed non thermal particle spectra, concluding that it has a number of major difficulties related to the transverse magnetic geometry of the shock wave. I discuss recent work on the inferred pair outflow rates, which are in excess of those predicted by existing theories of pair creation, and use those results to point out that the consequent mass loading of the wind reduces the wind's bulk flow 4-velocity to the point that dissipation of the magnetic field in a pulsar's wind upstream of the termination shock is restored to life as a viable model for the solution of the "sigma problem". I discuss some suggestions that current starvation in the current flow supporting the structured ("striped") upstream magnetic field induces a transition to superluminal wave propagation. I show that current starvation probably does not occur, because those currents are carried in the current sheet separating there stripes rather than in the stripes themselves.

I outline a new model of particle acceleration in the current sheet separating the closed from the open field lines in the force-free model of pulsar magnetospheres, based on reconnection at the light cylinder and "auroral" acceleration occurring in the return current channel that connects the light cylinder to the neutron star surface. I discuss recent studies of Pulsar Wind Nebulae, which find that pair outflow rates in excess of those predicted by existing theories of pair creation occur, and use those results to point out that dissipation of the magnetic field in a pulsar's wind upstream of the termination shock is restored to life as a viable model for the solution of the "$\sigma$" problem as a consequence of the lower wind 4-velocity implied by the larger mass loading.

We present a new model for the spectral evolution of Pulsar Wind Nebulae inside Supernova Remnants. The model couples the long-term dynamics of these systems, as derived in the 1-D approximation, with a 1-zone description of the spectral evolution of the emitting plasma. Our goal is to provide a simplified theoretical description that can be used as a tool to put constraints on unknown properties of PWN-SNR systems: a piece of work that is preliminary to any more accurate and sophisticated modeling. In the present paper we apply the newly developed model to a few objects of different ages and luminosities. We find that an injection spectrum in the form of a broken-power law gives a satisfactory description of the emission for all the systems we consider. More surprisingly, we also find that the intrinsic spectral break turns out to be at a similar energy for all sources, in spite of the differences mentioned above. We discuss the implications of our findings on the workings of pulsar magnetospheres, pair multiplicity and on the particle acceleration mechanism(s) that might be at work at the pulsar wind termination shock.

I discuss the scientific rationale and opportunities in the study of high energy particle accelerators away from the Earth; mostly, those outside the Solar System. I also briefly outline the features to be desired in telescopes used to probe accelerators studied by remote sensing.

We highlight recent theoretical and observational progress in several areas of neutron star astrophysics, and discuss the prospect for advances in the next decade.

This is a white paper submitted to the Stars and Stellar Evolution (SSE) Science Frontier Panel (SFP) of the NRC's Astronomy and Astrophysics 2010 Decadal Survey. The white paper is endorsed by the American Physical Society's (APS) Topical Group on Plasma Astrophysics (GPAP).

In this White Paper, we emphasize the need for and the important role of plasma astrophysics in the studies of formation, evolution of, and feedback by Active Galaxies. We make three specific recommendations: 1) We need to significantly increase the resolution of VLA, perhaps by building an EVLA-II at a modest cost. This will provide the angular resolution to study jets at kpc scales, where, for example, detailed Faraday rotation diagnosis can be done at 1GHz transverse to jets; 2) We need to build coordinated programs among NSF, NASA, and DOE to support laboratory plasma experiments (including liquid metal) that are designed to study key astrophysical processes, such as magneto-rotational instability (origin of angular momentum transport), dynamo (origin of magnetic fields), jet launching and stability. Experiments allowing access to relativistic plasma regime (perhaps by intense lasers and magnetic fields) will be very helpful for understanding the stability and dissipation physics of jets from Supermassive Black Holes; 3) Again through the coordinated support among the three Agencies, we need to invest in developing comprehensive theory and advanced simulation tools to study the accretion disks and jets in relativistic plasma physics regime, especially in connecting large scale fluid scale phenomena with relativistic kinetic dissipation physics through which multi-wavelength radiation is produced.

We use ideal axisymmetric relativistic magnetohydrodynamic simulations to calculate the spindown of a newly formed millisecond, B ~ 10^{15} G, magnetar and its interaction with the surrounding stellar envelope during a core-collapse supernova (SN) explosion. The mass, angular momentum, and rotational energy lost by the neutron star are determined self-consistently given the thermal properties of the cooling neutron star's atmosphere and the wind's interaction with the surrounding star. The magnetar drives a relativistic magnetized wind into a cavity created by the outgoing SN shock. For high spindown powers (~ 10^{51}-10^{52} ergs/s), the magnetar wind is super-fast at almost all latitudes, while for lower spindown powers (~ 10^{50} erg/s), the wind is sub-fast but still super-Alfvenic. In all cases, the rates at which the neutron star loses mass, angular momentum, and energy are very similar to the corresponding free wind values (<~ 30% differences), in spite of the causal contact between the neutron star and the stellar envelope. In addition, in all cases that we consider, the magnetar drives a collimated (~5-10 deg.) relativistic jet out along the rotation axis of the star. Nearly all of the spindown power of the neutron star escapes via this polar jet, rather than being transferred to the more spherical SN explosion. The properties of this relativistic jet and its expected late-time evolution in the magnetar model are broadly consistent with observations of long duration gamma-ray bursts (GRBs) and their associated broad-lined Type Ic SN.

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.

We study the long term evolution of magnetic fields generated by an initially unmagnetized collisionless relativistic $e^+e^-$ shock. Our 2D particle-in-cell numerical simulations show that downstream of such a Weibel-mediated shock, particle distributions are approximately isotropic, relativistic Maxwellians, and the magnetic turbulence is highly intermittent spatially, nonpropagating, and decaying. Using linear kinetic theory, we find a simple analytic form for these damping rates. Our theory predicts that overall magnetic energy decays like $(\omega_p t)^{-q}$ with $q \sim 1$, which compares favorably with simulations, but predicts overly rapid damping of short wavelength modes. Magnetic trapping of particles within the magnetic structures may be the origin of this discrepancy. We conclude that initially unmagnetized relativistic shocks in electron-positron plasmas are unable to form persistent downstream magnetic fields. These results put interesting constraints on synchrotron models for the prompt and afterglow emission from GRBs.

I outline a mechanism, akin to Weibel instabilities of interpenetrating beams, in which the neighboring current sheets in a striped wind from an oblique rotator interact through a two stream-like mechanism (a Weibel instability in flatland), to create an anomalous resistivity that heats the sheets and causes the magnetic field to diffusively annihilate in the wind upstream of the termination shock. The heating has consequences for observable unpulsed emission from pulsars.

I share a few reminiscences and observations of 40 years of Pulsars.

We study the long term evolution of magnetic fields generated by a collisionless relativistic $e^+e^-$ shock which is initially unmagnetized. Our 2D particle-in-cell numerical simulations show that downstream of such a Weibel-mediated shock, particle distributions are close to isotropic, relativistic Maxwellians, and the magnetic turbulence is highly intermittent spatially, with the non-propagating magnetic fields forming relatively isolated regions with transverse dimension $\sim 10-20$ skin depths. These structures decay in amplitude, with little sign of downstream merging. The fields start with magnetic energy density $\sim (0.1-0.2)$ of the upstream kinetic energy within the shock transition, but rapid downstream decay drives the fields to much smaller values, below $10^{-3}$ of equipartition after $10^3$ skin depths. In an attempt to construct a theory that follows field decay to these smaller values, we explore the hypothesis that the observed damping is a variant of Landau damping in an unmagnetized plasma. The model is based on the small value of the downstream magnetic energy density, which suggests that particle orbits are only weakly perturbed from straight line motion, if the turbulence is homogeneous. Using linear kinetic theory applied to electromagnetic fields in an isotropic, relativistic Maxwellian plasma, we find a simple analytic form for the damping rates, $\gamma_k$, in two and three dimensions for small amplitude, subluminous electromagnetic fields. We find that magnetic energy does damp due to phase mixing of current carrying particles as $(\omega_p t)^{-q}$ with $q \sim 1$. (abridged)

We model the interaction between the wind from a newly formed rapidly rotating magnetar and the surrounding progenitor. In the first few seconds after core collapse the magnetar inflates a bubble of plasma and magnetic fields behind the supernova shock, which expands asymmetrically because of the pinching effect of the toroidal magnetic field, as in PWNe, even if the host star is spherically symmetric. The degree of asymmetry depends on the ratio of the magnetic energy to the total energy in the bubble. We assume that the wind by newly formed magnetars inflating these bubbles is more magnetized than for PWNe. We show that for a magnetic to total power supplied by the central magnetar $\sim 0.1$ the bubble expands relatively spherically while for values greater than 0.3, most of the pressure in the bubble is exerted close to the rotation axis, driving a collimated outflow out through the host star. This can account for the collimation inferred from observations of long-duration gamma-ray bursts (GRBs). Given that the wind magnetization increases in time, we thus suggest that the magnetar-driven bubble initially expands relatively spherically (enhancing the energy of the associated supernova) while at late times it becomes progressivelymore collimated (producing the GRB). Similar processes may operate in more modestly rotating neutron stars to produce asymmetric supernovae and lower energy transients such as X-ray flashes.

We present time-dependent axisymmetric magnetohydrodynamic simulations of the interaction of a relativistic magnetized wind produced by a proto-magnetar with a surrounding stellar envelope, in the first $\sim 10$ seconds after core collapse. We inject a super-magnetosonic wind with $\dot E = 10^{51}$ ergs s$^{-1}$ into a cavity created by an outgoing supernova shock. A strong toroidal magnetic field builds up in the bubble of plasma and magnetic field that is at first inertially confined by the progenitor star. This drives a jet out along the polar axis of the star, even though the star and the magnetar wind are each spherically symmetric. The jet has the properties needed to produce a long-duration gamma-ray burst (GRB). At $\sim 5$ s after core bounce, the jet has escaped the host star and the Lorentz factor of the material in the jet at large radii $\sim 10^{11}$ cm is similar to that in the magnetar wind near the source. Most of the spindown power of the central magnetar escapes via the relativistic jet. There are fluctuations in the Lorentz factor and energy flux in the jet on $\sim 0.01-0.1$ second timescale. These may contribute to variability in GRB emission (e.g., via internal shocks).

I survey recent successes in the application of relativistic MHD and force-free electrodynamics to the modeling of the pulsars' rotational energy loss mechanism as well as to the structure and emission characteristics of Pulsar Wind Nebulae. I suggest that unsteady reconnection in the current sheet separating the closed from the open zones of the magnetosphere is responsible for the torque fluctuations observed in some pulsars, as well as for departures of the braking index from the canonical value of 3. I emphasize the significance of the boundary layer between the closed and open zones as the active site in the outer magnetopshere. I elaborate on the conflict between the models currently in use to interpret the gamma ray and X-ray pulses from these systems with the electric current flows found in the spin down models. Because the polar cap ``gap'' is the essential component in the supply of plasma to pulsar magnetospheres and to pulsar wind nebulae, I emphasize the importance of high sensitivity gamma ray observations of pulsars with core components of radio emission and high magnetospheric voltage, since these observations will look directly into the polar plasma production region. I also discuss the shock conversion of flow energy into the spectra of the synchrotron emitting particles in the Nebulae. I comment on the prospects for future developments and improvements in all these areas.

We model the interaction between the wind from a newly formed rapidly rotating magnetar and the surrounding supernova shock and host star. The dynamics is modeled using the two-dimensional, axisymmetric thin-shell equations. In the first ~10-100 seconds after core collapse the magnetar inflates a bubble of plasma and magnetic fields behind the supernova shock. The bubble expands asymmetrically because of the pinching effect of the toroidal magnetic field, just as in the analogous problem of the evolution of pulsar wind nebulae. The degree of asymmetry depends on E_mag/E_tot. The correct value of E_mag/E_tot is uncertain because of uncertainties in the conversion of magnetic energy into kinetic energy at large radii in relativistic winds; we argue, however, that bubbles inflated by newly formed magnetars are likely to be significantly more magnetized than their pulsar counterparts. We show that for a ratio of magnetic to total power supplied by the central magnetar L_mag/L_tot ~ 0.1 the bubble expands relatively spherically. For L_mag/L_tot ~ 0.3, however, most of the pressure in the bubble is exerted close to the rotation axis, driving a collimated outflow out through the host star. This can account for the collimation inferred from observations of long-duration gamma-ray bursts (GRBs). Outflows from magnetars become increasingly magnetically dominated at late times, due to the decrease in neutrino-driven mass loss as the young neutron star cools. We thus suggest that the magnetar-driven bubble initially expands relatively spherically, enhancing the energy of the associated supernova, while at late times it becomes progressively more collimated, producing the GRB.

We solve the time-dependent dynamics of axisymmetric, general relativistic MHD winds from rotating neutron stars. The mass loss rate is obtained self consistently as a solution of the MHD equations, subject to a finite thermal pressure at the stellar surface. Conditions are chosen to be representative of the neutrino driven phase in newly born magnetars, which have been considered as a possible engine for GRBs. We compute the angular momentum and energy losses as a function of $\sigma$ and compare them with the analytic expectation from the classical theory of pulsar winds. We observe the convergence to the force-free limit in the energy loss and we study the evolution of the closed zone for increasing magnetization. Results also show that the dipolar magnetic field and the presence of a closed zone do not modify significantly the acceleration and collimation properties of the wind.

We report the results of 1D particle-in-cell simulations of ultrarelativistic shock waves in proton-electron-positron plasmas. We consider magnetized shock waves, in which the upstream medium carries a large scale magnetic field, directed transverse to the flow. Relativistic cyclotron instability of each species as the incoming particles encounter the increasing magnetic field within the shock front provides the basic plasma heating mechanism. The most significant new results come from simulations with mass ratio $m_p/m_\pm = 100$. We show that if the protons provide a sufficiently large fraction of the upstream flow energy density (including particle kinetic energy and Poynting flux), a substantial fraction of the shock heating goes into the formation of suprathermal power-law spectra of pairs. Cyclotron absorption by the pairs of the high harmonic ion cyclotron waves, emitted by the protons, provides the non-thermal acceleration mechanism. As the proton fraction increases, the non-thermal efficiency increases and the pairs' power-law spectra harden. We suggest that the varying power law spectra observed in synchrotron sources powered by magnetized winds and jets might reflect the correlation of the proton to pair content enforced by the underlying electrodynamics of these sources' outflows, and that the observed correlation between the X-ray spectra of rotation powered pulsars with the X-ray spectra of their nebulae might reflect the same correlation.

We solve for the time-dependent dynamics of axisymmetric, general relativistic magnetohydrodynamic winds from rotating neutron stars. The mass loss rate is obtained self-consistently as a solution to the MHD equations, subject to a finite thermal pressure at the stellar surface. We consider both monopole and dipole magnetic field geometries and we explore the parameter regime extending from low magnetization (low-sigma_o), almost thermally-driven winds to high magnetization (high-sigma_o), relativistic Poynting-flux dominated outflows. We compute the angular momentum and rotational energy loss rates as a function of sigma_o and compare with analytic expectations from the classical theory of pulsars and magnetized stellar winds. In the case of the monopole, our high-sigma_o calculations asymptotically approach the analytic force-free limit. If we define the spindown rate in terms of the open magnetic flux, we similarly reproduce the spindown rate from recent force-free calculations of the aligned dipole. However, even for sigma_o as high as ~20, we find that the location of the Y-type point (r_Y), which specifies the radius of the last closed field line in the equatorial plane, is not the radius of the light cylinder R_L = c/omega (R = cylindrical radius), as has previously been assumed in most estimates and force-free calculations. Instead, although the Alfven radius at intermediate latitudes quickly approaches R_L as sigma_o exceeds unity, r_Y remains significantly less than R_L. Because r_Y < R_L, our calculated spindown rates thus exceed the classic ``vacuum dipole'' rate. We discussthe implications of our results for models of rotation-powered pulsars and magnetars, both in their observed states and in their hypothesized rapidly rotating initial state.

We use new and archival Chandra and ROSAT data to study the time variability of the X-ray emission from the pulsar wind nebula (PWN) powered by PSR B1509-58 on timescales of one week to twelve years. There is variability in the size, number, and brightness of compact knots appearing within 20" of the pulsar, with at least one knot showing a possible outflow velocity of ~0.6c (assuming a distance to the source of 5.2 kpc). The transient nature of these knots may indicate that they are produced by turbulence in the flows surrounding the pulsar. A previously identified prominent jet extending 12 pc to the southeast of the pulsar increased in brightness by 30% over 9 years; apparent outflow of material along this jet is observed with a velocity of ~0.5c. However, outflow alone cannot account for the changes in the jet on such short timescales. Magnetohydrodynamic sausage or kink instabilities are feasible explanations for the jet variability with timescale of ~1.3-2 years. An arc structure, located 30"-45" north of the pulsar, shows transverse structural variations and appears to have moved inward with a velocity of ~0.03c over three years. The overall structure and brightness of the diffuse PWN exterior to this arc and excluding the jet has remained the same over the twelve year span. The photon indices of the diffuse PWN and possibly the jet steepen with increasing radius, likely indicating synchrotron cooling at X-ray energies.

We present a time-lapse sequence of 20 near-infrared (J- and K'-band) snapshots of the central 20x20 arcsec of the Crab pulsar wind nebula, taken at subarcsecond resolution with the Hokupa'a/QUIRC adaptive optics camera on the Gemini North Telescope, and sampled at intervals of 10 min and 24 hr. It is observed that the equatorial wisps and polar knots in the termination shock of the pulsar wind appear to fluctuate in brightness on kilosecond time-scales. Maximum flux variations of +/-24 +/-4 and +/-14 +/-4 per cent relative to the mean (in 1.2 ks) are measured for the wisps and knots respectively, with greatest statistical significance in J band where the nebula background is less prominent. The J and K' flux densities imply different near-infrared spectra for the nonthermal continuum emission from the wisps and outermost polar knot (sprite), giving F_nu ~ nu^{-0.56+/-0.12} and F_nu ~ nu^{-0.21+/-0.13) respectively. The data are compared with existing optical and UV photometry and applied to constrain theories of the variability of the wisps (relativistic ion-cyclotron instability) and knots (relativistic fire hose instability).

We report on the first Green Bank Telescope observations at 427, 820 and 1400 MHz of the newly discovered, highly inclined and relativistic double pulsar binary. We focus on the brief eclipse of PSR J0737-3039A, the faster pulsar, when it passes behind PSR J0737-3039B. We measure a frequency-averaged eclipse duration of 26.6 +/- 0.6 s, or 0.00301 +/- 0.00008 in orbital phase. The eclipse duration is found to be significantly dependent on radio frequency, with eclipses longer at lower frequencies. Specifically, eclipse duration is well fit by a linear function having slope (-4.52 +/- 0.03) x 10^{-7} orbits/MHz. We also detect significant asymmetry in the eclipse. Eclipse ingress takes 3.51 +/- 0.99 times longer than egress, independent of radio frequency. Additionally, the eclipse lasts (40 +/- 7) x 10^{-5} in orbital phase longer after conjunction, also independent of frequency. We detect significant emission from the pulsar on short time scales during eclipse in some orbits. We discuss these results in the context of a model in which the eclipsing material is a shock-heated plasma layer within the slower PSR J0737-3039B's light cylinder, where the relativistic pressure of the faster pulsar's wind confines the magnetosphere of the slower pulsar.

The paper was withdrawn because the effect is due to the data reduction methodology. See the text file for details.

We propose synchrotron absorption in a magnetosheath forming a cocoon around the magnetosphere of pulsar B to be the origin of the eclipse phenomena seen in the recently discovered double pulsar system PSRJ07370-3039 A & B. The model predicts the eclipses will clear at frequencies higher than those of the observations reported to date (nominally, above $\nu \sim 5$ GHz.) The model also predicts synchrotron emission at the level of a few to 10 $\mu$Jy, peaking at $\nu \sim 2-5$ GHz with possible orbital modulation. We use simplified semi-analytic models and particle-in-cell simulations to elucidate the structure of the B magnetosphere, showing that B's magnetic field is confined to within a radius less than 50,000 km from B,on the side facing A. We estimate a ``propellor'' spindown torque on B, which is the dominant torque on this star, which yields a polar dipole field $\sim 7 \times 10^{11}$ Gauss (magnetic moment $\mu_B \sim 3.5 \times 10^{29}$ cgs). The model can explain the known eclipses only if the A wind's density is at least 4 orders of magnitude greater than is expected from popular pair creation models, and discuss the implications of this result for our general understanding of pulsar physics. Our proposal was outlined in Kaspi et al. (2004) and Demorest et al. (2004); a similar proposal has been made by Lyutikov (2004).

We use the two-dimensional Kolmogorov-Smirnov (KS) test to study the correlation between the 60 cosmic ray events above 4x10^19 eV from the AGASA experiment and the positions of infrared luminous galaxies from the IRAS PSCz catalog. These galaxies are expected to be hosts to gamma ray bursts (GRB) and magnetars, both of which are associated with core collapse supernovae and have been proposed as possible acceleration sites for ultra high energy cosmic rays. We find consistency between the models and the AGASA events to have been drawn from the same underlying distribution of positions on the sky with KS probabilities ~50%. Application of the same test to the 11 highest AGASA events above 10^20 eV, however, yields a KS probability of < 0.5%, rejecting the models at >99.5% significance level. Taken at face value, these highest energy results suggest that the existing cosmic ray events above 10^20 eV do not owe their origin to long burst GRBs, rapidly rotating magnetars, or any other events associated with core collapse supernovae. The larger data set expected from the AUGER experiment will test whether this conclusion is real or is a statistical fluke that we estimate to be at the 2 sigma level.

We describe results from time-dependent numerical modeling of the collisionless reverse shock terminating the pulsar wind in the Crab Nebula. We treat the upstream relativistic wind as composed of ions and electron-positron plasma embedded in a toroidal magnetic field, flowing radially outward from the pulsar in a sector around the rotational equator. The relativistic cyclotron instability of the ion gyrational orbit downstream of the leading shock in the electron-positron pairs launches outward propagating magnetosonic waves. Because of the fresh supply of ions crossing the shock, this time-dependent process achieves a limit-cycle, in which the waves are launched with periodicity on the order of the ion Larmor time. Compressions in the magnetic field and pair density associated with these waves, as well as their propagation speed, semi-quantitatively reproduce the behavior of the wisp and ring features described in recent observations obtained using the Hubble Space Telescope and the Chandra X-Ray Observatory. By selecting the parameters of the ion orbits to fit the spatial separation of the wisps, we predict the period of time variability of the wisps that is consistent with the data. When coupled with a mechanism for non-thermal acceleration of the pairs, the compressions in the magnetic field and plasma density associated with the optical wisp structure naturally account for the location of X-ray features in the Crab. We also discuss the origin of the high energy ions and their acceleration in the equatorial current sheet of the pulsar wind.

We report here the first polarimetric measurements of the pulsars in the J0737-3039 binary neutron star system using the Green Bank Telescope. We conclude both that the primary star (A) has a wide hollow cone of emission, which is an expected characteristic of the relatively open magnetosphere given its short spin period, and that A has a small angle between its spin and magnetic dipole axes, $4\pm 3$ degrees. This near alignment of axes suggests that A's wind pressure on B's magnetosphere will depend on orbital phase. This variable pressure is one mechanism for the variation of flux and profile shape of B with respect to the orbital phase that has been reported. The response of B to the A wind pressure will also depend on the particular side of its magnetosphere facing the wind at the spin phase when B is visible. This is a second possible mechanism for variability. We suggest that B may have its spin axis aligned with the orbital angular momentum owing to A's wind torque that contributes to its spindown. Monitoring the pulsars while geodetic precession changes spin orientations will provide essential evidence to test detailed theoretical models. We determine the Rotation Measures of the two stars to be $-112.3\pm 1.5$ and $-118\pm 12$ rad m$^{-2}$.

I argue that ideal MHD relativistic winds are always limited in practice to asymptotic 4-velocity $\gamma_\infty \approx \sigma_0^{1/3}$ and asymptotic magnetization $\sigma \sim \sigma_0^{2/3} \gg 1$, where $\sigma_0$ is the wind magnetization with respect to the rest energy density, evaluated at the light cylinder of the rotating, magnetized compact object that drives the flow. This suggests that the observed low value of he asymptotic $\sigma$ in the equatorial sectors of the winds driving Pulsar Wind Nebulae and the associated high values of the asymptotic 4-velocity are a consequence of magentic dissipation in the wind zone.

I show that the relativistic winds of newly born magnetars with khz initial spin rates, occurring in all normal galaxies, can accelerate ultrarelativistic light ions with an E^{-1} injection spectrum, steepening to E^{-2} at higher energies, with an upper cutoff above 10^{21} eV. Interactions with the CMB yield a spectrum in good accord with the observed spectrum of Ultra-High Energy Cosmic Rays (UHECR), if ~ 5-10% of the magnetars are born with voltages sufficiently high to accelerate the UHECR. The form the spectrum spectrum takes depends on the gravitational wave losses during the magnetars' early spindown - pure electromagnetic spindown yields a flattening of the E^3 J(E) spectrum below 10^{20} eV, while a moderate GZK ``cutoff'' appears if gravitational wave losses are strong enough. I outline the physics such that the high energy particles escape with small energy losses from a magnetar's natal supernova, including Rayleigh-Taylor ``shredding'' of the supernova envelope, expansion of a relativistic blast wave into the interstellar medium, acceleration of the UHE ions through surf-riding in the electromgnetic fields of the wind, and escape of the UHE ions in the rotational equator with negligible radiation loss. The abundance of interstellar supershells and unusually large supernova remnants suggests that most of the initial spindown energy is radiated in khz gravitational waves for several hours after each supernova, with effective strains from sources at typical distances ~ 3 x 10^{-21}. Such bursts of gravitational radiation should correlate with bursts of ultra-high energy particles. The Auger experiment should see such bursts every few years.

I discuss recent work on the nature of relativistic winds from Rotation Powered Pulsars, on the physics of how they transport energy from the central compact object to the surrrounding world, and how that energy gets converted into observable synchrotron emission.

We present initial results of the first self-consistent numerical model of the outer magnetosphere of a pulsar. By using the relativistic ``particle-in-cell'' method with special boundary conditions to represent plasma dynamics in 3D, we are able to follow magnetospheric plasma through the light cylinder into the wind zone for arbitrary magnetic inclination angles. For aligned rotators we confirm the ``disk-dome'' charge-separated structure of the magnetosphere and find that this configuration is unstable to a 3D nonaxisymmetric diocotron instability. This instability allows plasma to move across the field lines and approach the corotating Goldreich-Julian solution within several rotation periods. For oblique rotators formation of the spiral ``striped wind'' in the equatorial direction is demonstrated and the acceleration of the wind and its magnetization is discussed. We find that the wind properties vary with stellar latitude; however, whether injection conditions at the pulsar are responsible for the observed jet-equator geometry of Crab and Vela is currently under investigation. We also comment on the electrodynamics of the simulated magnetospheres, their current closure, and future simulations.

We present observations with the Chandra X-ray Observatory of the pulsar wind nebula (PWN) powered by the energetic young pulsar B1509-58. These data confirm the complicated morphology of the system indicated by previous observations, and in addition reveal several new components to the nebula. The overall PWN shows a clear symmetry axis oriented at a position angle 150 degrees (N through E), which we argue corresponds to the pulsar spin axis. We show that a previously identified radio feature matches well with the overall extent of the X-ray PWN, and propose the former as the long-sought radio nebula powered by the pulsar. We further identify a bright collimated feature, at least 4' long, lying along the nebula's main symmetry axis; we interpret this feature as a physical outflow from the pulsar, and infer a velocity for this jet >0.2c. The lack of any observed counter-jet implies that the pulsar spin axis is inclined at approx 30 deg to the line-of-sight, contrary to previous estimates made from lower-resolution data. We also identify a variety of compact features close to the pulsar. A pair of semi-circular X-ray arcs lie 17 and 30 arcsec to the north of the pulsar; the latter arc shows a highly-polarized radio counterpart. We show that these features can be interpreted as ion-compression wisps in a particle-dominated equatorial flow, and use their properties to infer a ratio of electromagnetic to particle energy in pairs at the wind shock sigma approx 0.005, similar to that seen in the Crab Nebula. We further identify several compact knots seen very close to the pulsar; we use these to infer sigma < 0.003 at a separation from the pulsar of 0.1 pc.

The young and energetic pulsar B1509-58 powers a bright X-ray synchrotron nebula, embedded in the unusual supernova remnant G320.4-1.2. We present observations of this system with the Chandra X-ray Observatory, which show a spectacularly complicated source. The nebula is dominated by a bright collimated feature which we interpret as a relativistic jet directed along the pulsar spin axis. Several compact knots can be seen in the immediate vicinity of the pulsar. While many of these features are similar to those seen around the Crab pulsar, the nebula surrounding PSR B1509-58 shows important differences which are possibly a result of the the latter's low nebular magnetic field and low density environment.

We discuss the creation of electron-positron cascades in the context of pulsar polar cap acceleration models and derive several useful analytic and semi-analytic results for the spatial extent and energy response of the cascade. Instead of Monte Carlo simulations, we use an integro- differential equation which describes the development of the cascade energy spectrum in one space dimension quite well, when it is compared to existing Monte Carlo models. We reduce this full equation to a single integral equation, from which we can derive useful results, such as the energy loss between successive generations of photons and the spectral index of the response. We find that a simple analytic formula represents the pair cascade multiplicity quite well, provided that the magnetic field is below 10^12 Gauss, and that an only slightly more complex formula matches the numerically-calculated cascade at all other field strengths. Using these results, we find that cascades triggered by gamma rays emitted through inverse Compton scattering of thermal photons from the neutron star's surface, both resonant and non-resonant, are important for the dynamics of the polar cap region in many pulsars.

Through a simple model of particle acceleration and pair creation above the polar caps of rotation-powered pulsars, we calculate the height of the pair-formation front (PFF) and the dominant photon emission mechanism for the pulsars in the Princeton catalog. We find that for most low- and moderate-field pulsars, the height of the pair formation front and the final Lorentz factor of the primary beam is set by nonresonant inverse Compton scattering (NRICS), in the Klein-Nishina limit. NRICS is capable of creating pairs over a wide range of pulsar parameters without invoking a magnetic field more complicated than a centered dipole, although we still require a reduced radius of curvature for most millisecond pulsars. For short-period pulsars, the dominant process is curvature radiation, while for extremely high-field pulsars, it is resonant inverse Compton scattering (RICS). The dividing point between NRICS dominance and curvature dominance is very temperature-dependent; large numbers of pulsars dominated by NRICS at a stellar temperature of $10^6$ K are dominated by curvature at $10^5$ K. We apply these results to pulsar death-line calculations and to the issue of particle injection into the Crab Nebula.

We re-examine the characteristic polarization angle sweep of rotation-powered pulsars and calculate the expected deviations from this sweep caused by aberrational effects and by polar-cap current flow. We find that in addition to the previously known phase shift of the entire sweep by $\Delta \Phi = -4 r/R_L$, aberration shifts the polarization angle itself by $\Delta \Psi = -(10/3) (r/R_L) \cos \alpha$. Similarly, current flow above the polar cap shifts the polarization sweep by $\Delta \Psi = (10/3) (r/R_L) (J/J_{GJ}) \cos \alpha$, potentially providing a method of directly measuring the magnitude of the current. The competition between these two effects produces a potentially observable signature in the polarization angle sweep. Although these effects may appear similar to orthogonal mode shifts, they are an independent phenomenon with distinct observational characteristics.

I summarize the theory of acceleration of non-neutral particle beams by starvation electric fields along the polar magnetic field lines of rotation powered pulsars, including the effect of dragging of inertial frames which dominates the acceleration of a space charge limited beam. I apply these results to a new calculation of the radio pulsar death line, under the hypotheses that pulsar ``death'' corresponds to cessation of pair creation over the magnetic poles {\it and} that the magnetic field has a locally dipolar topology. The frame dragging effect in star centered dipole geometry does improve comparison of the theory with observation, but an unacceptably large conflict between observation and theory still persists. Offsetting the dipole improves the comparison, but a fully satisfactory theory requires incorporating magnetic conversion of inverse Compton gamma rays, created by scattering thermal photons from the surface of old neutron stars ($t > 10^8 $ years) kept warm ($T \geq 10^5$ K) by friction between the rotating core and the crust. The result is a ``death valley'' for pulsars; offsets of the dipole center from the stellar center in the oldest stars $\sim (0.7-0.8) R_*$ suffice. The resulting theory predicted the existence of rotation powered pulsars with these advanced ages, a prediction confirmed by the recent discovery that PSR J2144-3933 actually has a rotation period of 8.5 seconds.

We report the discovery of kHz fluctuations, including quasi-periodic oscillations (QPO) at ~330 Hz and ~760 Hz and a broadband kHz continuum in the power density spectrum of the high mass X-ray binary pulsar Centaurus X-3. These observations of Cen X-3 were carried out with the Rossi X-ray Timing Explorer (RXTE). The fluctuation spectrum is flat from mHz to a few Hz, then steepens to $f^{-2}$ behavior between a few Hz and ~100 Hz. Above a hundred Hz, the spectrum shows the QPO features, plus a flat continuum extending to ~1200 Hz and then falling out to ~1800 Hz. These results, which required the co-adding three days of observations of Cen X-3, are at least as fast as the fastest known variations in X-ray emission from an accreting compact object (kHz QPO in LMXB sources) and probably faster since extension to ~1800 Hz is indicated by the most likely parameterization of the data. Multi-dimensional radiation hydrodynamics simulations of optically thick plasma flow onto the magnetic poles of an accreting neutron star show that the fluctuations at frequencies above 100 Hz are consistent with photon bubble turbulence and oscillations (PBO) previously predicted to be observable in this source. For a polar cap opening angle of 0.25 radians, we show that the spectral form above 100 Hz is reproduced by the simulations, including the frequencies of the QPO and the relative power in the QPO and the kHz continuum. This has resulted in the first model-dependent measurement of the polar cap size of an X-ray pulsar.

After discussion of observational constraints on the nature of the MHD wind coupling between the Crab Pulsar and the Crab Nebula, the theory of transverse relativistic shock structure is reviewed and applied to the interpretation of the wisps in the Nebula as the manifestation of the distributed wind termination shock structure, energetically dominated by heavy ions, accelerated in the rotational equator of the pulsar to energies comparable to the total voltage across the pulsar's open field lines and carrying a current comparable to the Goldreich- Julian current. New results on the variability of the shock structure are presented, which show that the gyrating ion bunches emit outwardly traveling finite amplitude compressional waves, in agreement with recent ground based observations. The implications of the theory for X-ray, $\gamma$-ray and high energy neutrino emission are briefly discussed, as are the problems of low magnetic energy density in the upstream wind and the origin of the Nebular radio emission. A brief discussion of other plerions leads to the conclusion that much more detailed observations are needed before these systems can be modeled with the same sophistication as can be done for the Crab Nebula.

I summarize the theory of acceleration of non-neutral particle beams by starvation electric fields along the polar magnetic field lines of rotation powered pulsars, including the effect of dragging of inertial frames which dominates the acceleration of a space charge limited beam. I apply these acceleration results to a new calculation of the radio pulsar death line, under the hypotheses that pulsar ``death'' corresponds to cessation of pair creation over the magnetic poles and that the magnetic field has a locally dipolar topology. While the frame dragging effect in star centered dipole geometry does improve comparison of the theory with observation, an unacceptably large fraction of the observed stars outside the bounds of pair creation theory still persists. Offsetting the dipole improves the correspondence between theory and observation. The result is a ``death valley'' for pulsars; acceptable comparison of observation and theory occurs if the boundary of death valley corresponds to offsets of the dipole center from the stellar center $\sim (0.7-0.8) R_*$. I also point out that pulsars are absent for magnetic moments corresponding to star centered polar fields in excess of $\sim 4 \times 10^{13}$ Gauss, and I suggest that this absence is due to pairs forming as bound positronium atoms in such strong fields, creating a neutral, relativistically outflowing gas which cannot participate in low altitude collective radio emission processes in such strongly magnetized objects.

We report the results of a CGRO 3-week observation of the binary system containing the 47 ms pulsar PSR B1259-63 orbiting around a Be star companion in a very eccentric orbit. The PSR B1259-63 binary is a unique system for the study of the interaction of a rapidly rotating pulsar with time-variable nebular surroundings. CGRO observed the PSR B1259-63 system in coincidence with its most recent periastron passage (January 3-23, 1994). Unpulsed and non-thermal hard X-ray emission was detected up to 200 keV, with a photon index $1.8 \pm 0.2$ and a flux of ~4 mCrab, corresponding to a luminosity of a few 10^{34} erg/s at the distance of 2 kpc. The hard X-ray flux and spectrum detected by CGRO agrees with the X-ray emission measured by simultaneous ASCA observations. EGRET upper limits are significant, and exclude strong inverse Compton cooling in the PSR B1259-63 system. We interpret the observed non-thermal emission as synchrotron radiation of shocked electron/positron pairs of the relativistic pulsar wind interacting with the mass outflow from the Be star. Our results clearly indicate, for the first time in a binary pulsar, that high energy emission can be shock-powered rather than caused by accretion. The lack of X-ray/gamma-ray pulsations constrains models of high-energy emission from rapidly rotating pulsars.

We study the physical processes of the PSR B1259-63 system containing a 47 ms pulsar orbiting around a Be star in a highly eccentric orbit. Motivated by the results of a multiwavelength campaign during the January 1994 periastron passage of PSR B1259-63, we discuss several issues regarding the mechanism of high-energy emission. Unpulsed power law emission from the this system was detected near periastron in the energy range 1-200 keV. We find that the observed high energy emission from the PSR B1259-63 system is not compatible with accretion or propeller-powered emission. Shock-powered high-energy emission produced by the pulsar/outflow interaction is consistent with all high energy observations. By studying the evolution of the pulsar cavity we constrain the magnitude and geometry of the mass outflow outflow of the Be star. The pulsar/outflow interaction is most likely mediated by a collisionless shock at the internal boundary of the pulsar cavity. The system shows all the characteristics of a {\it binary plerion} being {\it diffuse} and {\it compact} near apastron and periastron, respectively. The PSR B1259-63 cavity is subject to different radiative regimes depending on whether synchrotron or inverse Compton (IC) cooling dominates the radiation of electron/positron pairs advected away from the inner boundary of the pulsar cavity. The highly non-thermal nature of the observed X-ray/gamma-ray emission near periastron establishes the existence of an efficient particle acceleration mechanism within a timescale shown to be less than $\sim 10^2-10^3$ s. A synchrotron/IC model of emission of e\pm-pairs accelerated at the inner shock front of the pulsar cavity and adiabatically expanding in the MHD flow provides an excellent explanation of the observed time variableX-ray flux and spectrum from the PSR