Jonathan Arons

Relativistic magnetosonic shock waves in synchrotron sources - Shock structure and nonthermal acceleration of positrons

The theoretical properties of relativistic, transverse, magnetosonic collisionless shock waves in electron-positron-heavy ion plasmas of relevance to astrophysical sources of synchrotron radiation are investigated. Both 1D electromagnetic particle-in-cell simulations and quasi-linear theory are used to examine the spatial and kinetic structure of these nonlinear flows. A new process of shock acceleration of nonthermal positrons, in which the gyrating reflected heavy ions dissipate their energy in the form of collectively emitted, left-handed magnetosonic waves which are resonantly absorbed by the positrons immediately behind the ion reflection region, is described. Applications of the results to the termination shocks of pulsar winds and to the termination shocks of jets emanating from the AGN are outlined. 61 refs.

Magnetars in the Metagalaxy: An Origin for Ultra-High-Energy Cosmic Rays in the Nearby Universe

I show that the relativistic winds of newly born magnetars (neutron stars with petagauss surface magnetic fields) with initial spin rates close to the centrifugal breakup limit, occurring in all normal galaxies with massive star formation, can provide a source of ultrarelativistic light ions with an E-1 injection spectrum, steepening to E-2 at higher energies, with an upper cutoff at 1021-1022 eV. Interactions with the cosmic microwave background yield a spectrum at the Earth that compares favorably with the spectrum of ultra-high-energy cosmic rays (UHECRs) observed at energies up to a few times 1020 eV. The fit to the observations suggests that ~5%-10% of the magnetars are born with rotation rates and voltages sufficiently high to allow the acceleration of the UHECR. The form the spectrum incident on the Earth takes depends sensitively on the mechanism and the magnitude of gravitational wave losses during the early spin-down of these neutron stars: pure electromagnetic spin-down (the E-1 injection spectrum) yields a GZK feature [a flattening of the E3J(E) spectrum] below 1020 eV, rather than a cutoff, while a moderate GZK cutoff appears if gravitational wave losses are strong enough to steepen the injection spectrum above 1020 eV. The flux above 1020 eV comes from magnetars in relatively nearby galaxies (D 100 EeV air showers, the model predicts gravitational wave strains ~3 × 10-21. Such bursts of gravitational radiation should correlate with bursts of ultra-high-energy particles. The Auger experiment should see bursts of particles with energy above 100 EeV every few years.

The Maximum Energy of Accelerated Particles in Relativistic Collisionless Shocks

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 ? 10?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 ? 10?3. For electron-ion plasmas, the transition to efficient acceleration occurs for ? 3 ? 10?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 ?maxt 1/2. This scaling is shallower than the so-called (and commonly assumed) Bohm limit ?maxt, 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 ?sat/?0 mic 2 ~ ??1/4, where ?0 mic 2 is the mean energy per particle in the upstream bulk flow. Further energization is prevented by the fact that the self-generated turbulence is confined within a finite region of thickness ??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.

Long-Term Evolution of Magnetic Turbulence in Relativistic Collisionless Shocks: Electron-Positron Plasmas

We study the long term evolution of magnetic fields generated by an initially unmagnetized collisionless relativistic e+e− shock. Our two-dimensional 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. The nonpropagating magnetic fields decay in amplitude and do not merge. The fields start with magnetic energy density ~ 0.1-0.2 of equipartition, but rapid downstream decay drives the fields to much smaller values, below ~10−3 of equipartition after ~103 skin depths. To construct a theory to follow field decay to these smaller values, we hypothesize that the observed damping is a variant of Landau damping. The model is based on the small value of the downstream magnetic energy density, which only weakly perturbs particle orbits, for homogeneous turbulence. Using linear kinetic theory, we find a simple analytic form for the damping rates for small-amplitude, subluminous electromagnetic fields. Our theory predicts that overall magnetic energy decays as (ωpt)−q with q ~ 1, which compares with simulations. However, our theory predicts overly rapid damping of short-wavelength modes. Magnetic trapping of particles within the highly spatially intermittent downstream magnetic structures may be the origin of this discrepancy and may allow for some of this initial magnetic energy to persist. Absent additional physical processes that create longer wavelength, more persistent fields, 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. We also comment on the relevance of these results for relativistic electron-ion shocks.

Pair Production Multiplicities in Rotation-powered Pulsars

We discuss the creation of electron-positron cascades in the context of pulsar polar cap acceleration models and derive several useful analytic and semianalytic results for the spatial extent and energy response of the cascade. Instead of Monte Carlo simulations, we use an integrodifferential equation that 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 1012 G 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 γ-rays emitted through inverse Compton scattering of thermal photons from the neutron stars surface, both resonant and nonresonant, are important for the dynamics of the polar cap region in many pulsars. In these objects, the expected multiplicity of pairs generated by a single input particle is lower than previously found in cascades initiated by curvature emission, frequently being on the order of 10 rather than ~1000 as usually quoted. Such pulsars also are expected to be less luminous in polar cap γ-rays than when curvature emission triggers the cascade, a topic that will be the subject of a subsequent paper.

Pair Multiplicities and Pulsar Death

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 PFF 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 106 K are dominated by curvature at 105 K. Our principle result is a new determination of the theoretical pulsar death line. Proper inclusion of ICS allows us to reach the conclusion that all known radio pulsars are consistent with pair creation above their polar caps, assuming steady acceleration of a space charge-limited particle beam. We identify a new region of the P- diagram where slow pulsars with narrow radiation beams should be found in sufficiently large surveys. We also show that the injection rate of electrons and positrons into the Crab Nebula inferred from the polar cap pair creation model at the present epoch (1039 electrons and positrons s-1) suffices to explain the nebular X-ray and ?-ray emission, the best empirical measure of the instantaneous particle-loss rate from any pulsar, but that the contemporary injection rate is about a factor of 5 below the rate averaged over the nebulas history required to explain the nebular radio emission (assuming that the nebular radio source is homogeneous). It is not clear whether this discrepancy can be resolved by evolutionary effects and by better treatment of nebular inhomogeneity or is an indication of another particle source in the pulsars magnetosphere.

Magnetized relativistic jets and long‐duration GRBs from magnetar spin‐down during core‐collapse supernovae

We use ideal axisymmetric relativistic magnetohydrodynamic simulations to calculate the spin-down 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 stars atmosphere and the winds interaction with the surrounding star. The magnetar drives a relativistic magnetized wind into a cavity created by the outgoing SN shock. For high spin-down powers (∼10 51 ―10 52 erg s ―1 ), the magnetar wind is superfast at almost all latitudes, while for lower spin-down powers (∼10 50 erg s ―1 ), the wind is subfast 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 per cent 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°) relativistic jet out along the rotation axis of the star. Nearly all of the spin-down 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.

Some problems of pulsar physics or I'm madly in love with electricity

Some current topics in the theory of pulsar magnetospheres and their emission are reviewed. The mode of plasma supply and its consequences for structure of planetary and stellar magnetospheres is discussed. In the pulsar case, the plasma is supplied by electrical forces, in contrast to all other known examples. The resulting theories of particle acceleration along polar field lines are then reviewed, and the total energization of the charge separated plasma is summarized, when pair creation is absent. The effects of pair creation are reviewed using models of the resulting steady and unsteady flows, when the polar zones of the pulsar emit either electrons or ions. The application of these theories of acceleration and plasma supply to pulsars is discussed, with particular attention paid to the total amount of electron-positron plasma created and its momentum distribution. Qualitative agreement is shown between the spatial structure of the relativistically outflowing plasma described in one version of these models and the morphology of pulsar wave forms. Various aspects of radiation emission and transport are summarized, based on the polar current flow model with pair creation, and the phenomenon of marching subpulses is discussed. The corotation beaming and the relativistically expanding current sheet models for pulsar emission are also discussed briefly, and the paper concludes with a brief discussion of the relation between the theories of polar flow with pair plasma and the problem of the energization of the Crab Nebula.

High-energy emission from the eclipsing millisecond pulsar PSR 1957+20

The properties of the high-energy emission expected from the eclipsing millisecond pulsar system PSR 1957+20 are investigated. Emission is considered by both the relativistic shock produced by the pulsar wind in the nebula surrounding the binary and by the shock constraining the mass outflow from the companion star of PSR 1957+20. On the basis of the results of microscopic plasma physical models of relativistic shocks it is suggested that the high-energy radiation is produced in the range from X-rays to MeV gamma rays in the binary and in the range from 0.01 eV to about 40 keV in the nebula. Doppler boost of the emission in the radiating wind suggests the flux should vary on the orbital time scale, with the largest flux observed roughly coincident with the pulsars radio eclipse.

Relativistic, perpendicular shocks in electron-positron plasmas

One-dimensional particle-in-cell plasma simulations are used to examine the mechanical structure and thermalization properties of collisionless relativistic shock waves in electron-positron plasmas. Shocks propagating perpendicularly to the magnetic field direction are considered. It is shown that these shock waves exist, and that they are completely parameterized by the ratio of the upstream Poynting flux to the upstream kinetic energy flux. The way in which the Rankine-Hugoniot shock jump conditions are modified by the presence of wave fluctuations is shown, and they are used to provide a macroscopic description of these collisionless shock flows. The results of a 2D simulation that demonstrates the generality of these results beyond the assumption of the 1D case are discussed. It is suggested that the thermalization mechanism is the formation of a synchrotron maser by the coherently reflected particles in the shock front. Because the downstream medium is thermalized, it is argued that perpendicular shocks in pure electron-positron plasmas are not candidates as nonthermal particle accelerators. 40 refs.