We study diffusive shock acceleration (DSA) of electrons in non-relativistic quasi-perpendicular shocks using self-consistent one-dimensional particle-in-cell (PIC) simulations. By exploring the parameter space of sonic and Alfv\'{e}nic Mach numbers we find that high Mach number quasi-perpendicular shocks can efficiently accelerate electrons to power-law downstream spectra with slopes consistent with DSA prediction. Electrons are reflected by magnetic mirroring at the shock and drive non-resonant waves in the upstream. Reflected electrons are trapped between the shock front and upstream waves and undergo multiple cycles of shock drift acceleration before the injection into DSA. Strong current-driven waves also temporarily change the shock obliquity and cause mild proton pre-acceleration even in quasi-perpendicular shocks, which otherwise do not accelerate protons. These results can be used to understand nonthermal emission in supernova remnants and intracluster medium in galaxy clusters.

Since pulsars were discovered as emitters of bright coherent radio emission more than half a century ago, the cause of the emission has remained a mystery. In this Letter we demonstrate that coherent radiation can be directly generated in non-stationary pair plasma discharges which are responsible for filling the pulsar magnetosphere with plasma. By means of large-scale two-dimensional kinetic plasma simulations, we show that if pair creation is non-uniform across magnetic field lines, the screening of electric field by freshly produced pair plasma is accompanied by the emission of waves which are electromagnetic in nature. Using localized simulations of the screening process, we identify these waves as superluminal ordinary (O) modes, which should freely escape from the magnetosphere as the plasma density drops along the wave path. The spectrum of the waves is broadband and the frequency range is comparable to that of observed pulsar radio emission.

Plasmoids -- magnetized quasi-circular structures formed self-consistently in reconnecting current sheets -- were previously considered to be the graveyards of energetic particles. In this paper, we demonstrate the important role of plasmoids in shaping the particle energy spectrum in relativistic reconnection (i.e., with upstream magnetization $\sigma_{\rm up} \gg 1$). Using two dimensional particle-in-cell simulations in pair plasmas with $\sigma_{\rm up}=10$ and $100$, we study a secondary particle energization process that takes place inside compressing plasmoids. We demonstrate that plasmoids grow in time, while their interiors compress, amplifying the internal magnetic field. The magnetic field felt by particles injected in an isolated plasmoid increases linearly with time, which leads to particle energization as a result of magnetic moment conservation. For particles injected with a power-law distribution function, this energization process acts in such a way that the shape of the injected power law is conserved, while producing an additional non-thermal tail $f(E)\propto E^{-3}$ at higher energies followed by an exponential cutoff. The cutoff energy, which increases with time as $E_{\rm cut}\propto\sqrt{t}$, can greatly exceed $\sigma_{\rm up} m_e c^2$. We analytically predict the secondary acceleration timescale and the shape of the emerging particle energy spectrum, which can be of major importance in certain astrophysical systems, such as blazar jets.

Recent advances in numerical algorithms and computational power have enabled first-principles simulations of pulsar magnetospheres using Particle-in-Cell (PIC) techniques. These ab-initio simulations seem to indicate that pair creation through photon-photon collision at the light cylinder is required to sustain the pulsar engine. However, for many rotation-powered pulsars pair creation operates effectively only near the stellar surface where magnetic field is high. How these "weak pulsars" fill their magnetospheres without efficient photon-photon pair conversion in the outer magnetosphere is still an open question. In this paper, we present a range of self-consistent solutions to the pulsar magnetosphere that do not require pair production near the light cylinder. When pair production is very efficient near the star, the pulsar magnetosphere converges to previously-reported solutions. However, in the intermediate regime where pair supply is barely enough to sustain the magnetospheric current, we observe a time-dependent solution with quasi-period about half of the rotation period. This new quasi-periodic solution may explain the observed pulsar death line without invoking multipolar components near the star, and can potentially explain the core vs. conal emission patterns observed in pulsar radio signals.

The self-regulation of cosmic ray (CR) transport in the interstellar and intracluster media has long been viewed through the lenses of linear and quasilinear kinetic plasma physics. Such theories are believed to capture the essence of CR behavior in the presence of self-generated turbulence, but cannot describe potentially critical details arising from the nonlinearities of the problem. We utilize the particle-in-cell numerical method to study the time-dependent nonlinear behavior of the gyroresonant streaming instabilities, self-consistently following the combined evolution of particle distributions and self-generated wave spectra in one-dimensional periodic simulations. We demonstrate that the early growth of instability conforms to the predictions from linear physics, but that the late-time behavior can vary depending on the properties of the initial CR distribution. We emphasize that the nonlinear stages of instability depend strongly on the initial anisotropy of CRs -- highly anisotropic CR distributions do not efficiently reduce to Alfvenic drift velocities, owing to reduced production of left-handed resonant modes. We derive estimates for the wave amplitudes at saturation and the time scales for nonlinear relaxation of the CR distribution, then demonstrate the applicability of these estimates to our simulations. Bulk flows of the background plasma due to the presence of resonant waves are observed in our simulations, confirming the microphysical basis of CR-driven winds.

This is a science white paper submitted to the Astro-2020 and Plasma-2020 Decadal Surveys. The paper describes the present status and emerging opportunities in Extreme Plasma Astrophysics -- a study of astrophysically-relevant plasma processes taking place under extreme conditions that necessitate taking into account relativistic, radiation, and QED effects.

Magnetic reconnection is often invoked to explain the non-thermal radiation of relativistic outflows, including jets of active galactic nuclei (AGN). Motivated by the largely unknown plasma composition of AGN jets, we study reconnection in the unexplored regime of electron-positron-proton (pair-proton) plasmas with large-scale two-dimensional particle-in-cell simulations. We cover a wide range of pair multiplicities (lepton-to-proton number ratio $\kappa=1-199$) for different values of the all-species plasma magnetization ($\sigma=1,3$ and 10) and electron temperature ($\Theta_e\equiv kT_e/m_ec^2=0.1-100$). We focus on the dependence of the post-reconnection energy partition and lepton energy spectra on the hot pair plasma magnetization $\sigma_{e,h}$ (i.e., the ratio of magnetic to pair enthalpy densities). We find that the post-reconnection energy is shared roughly equally between magnetic fields, pairs, and protons for $\sigma_{e,h}\gtrsim 3$. We empirically find that the mean lepton Lorentz factor in the post-reconnection region depends on $\sigma, \Theta_e$, and $\sigma_{e,h}$ as $\langle \gamma_e-1\rangle \approx \sqrt{\sigma}(1+4\Theta_e)\left(1+\sigma_{e,h}/30\right)$, for $\sigma\ge1$. The high-energy part of the post-reconnection lepton energy distributions can be described by a power law, whose slope is mainly controlled by $\sigma_{e,h}$ for $\kappa \gtrsim 3-6$, with harder power laws obtained for higher magnetizations. We finally show that reconnection in pair-proton plasmas with multiplicities $\kappa \sim 1-20$, magnetizations $\sigma \sim 1-10$, and temperatures $\Theta_e \sim 1-10$ results in particle power law slopes and average electron Lorentz factors that are consistent with those inferred in leptonic models of AGN jet emission.

In some Seyfert Galaxies, the hard X-rays that produce fluorescent emission lines are thought to be generated in a hot corona that is compact and located at only a few gravitational radii above the supermassive black hole. We consider the possibility that this X-ray source may be powered by small scale magnetic flux tubes attached to the accretion disk near the black hole. We use three dimensional, time dependent force-free simulations in a simplified setting to study the dynamics of such flux tubes as they get continuously twisted by the central compact star/black hole. We find that, the dynamical evolution of the flux tubes connecting the central compact object and the accretion disk is strongly influenced by the confinement of the surrounding field. Although differential rotation between the central object and the disk tends to inflate the flux tubes, strong confinement from surrounding field quenches the formation of a jet-like outflow, as the inflated flux tube becomes kink unstable and dissipates most of the extracted rotational energy relatively close to the central object. Such a process may be able to heat up the plasma and produce strong X-ray emission. We estimate the energy dissipation rate and discuss its astrophysical implications.

We present the first realistic 3D simulations of flame front instabilities during type I X-ray bursts. The unperturbed front is characterised by the balance between the pressure gradient and the Coriolis force of a spinning neutron star ({\nu} = 450 Hz in our case). This balance leads to a fast horizontal velocity field parallel to the flame front. This flow is strongly sheared in the vertical direction. When we perturb the front an instability quickly corrugates the front. We identify this instability as the baroclinic instability. Most importantly, the flame is not disrupted by the instability and there are two major consequences: the overall flame propagation speed is {\sim} 10 times faster than in the unperturbed case and distinct flame vortices appear. The speedup is due to the corrugation of the front and the dynamics of the vortices. These vortices may also be linked to the oscillations observed in the lightcurves of the bursts.

This white paper, submitted for the Plasma 2020 Decadal Survey, concerns the physics of weakly collisional, high-beta plasmas -- plasmas in which the thermal pressure dominates over the magnetic pressure and in which the inter-particle collision time is comparable to the characteristic timescales of bulk motions. This state of matter, although widespread in the Universe, remains poorly understood: we lack a predictive theory for how it responds to perturbations, how it transports momentum and energy, and how it generates and amplifies magnetic fields. Such topics are foundational to the scientific study of plasmas, and are of intrinsic interest to those who regard plasma physics as a fundamental physics discipline. But these topics are also of extrinsic interest: addressing them directly informs upon our understanding of a wide variety of space and astrophysical systems, including accretion flows around supermassive black holes, the intracluster medium (ICM) between galaxies in clusters, and regions of the near-Earth solar wind. Specific recommendations to advance this field of study are discussed.

In this Letter we propose that coherent radio emission of Crab, other young energetic pulsars, and millisecond pulsars is produced in the magnetospheric current sheet beyond the light cylinder. We carry out global and local two-dimensional kinetic plasma simulations of reconnection to illustrate the coherent emission mechanism. Reconnection in the current sheet beyond the light cylinder proceeds in the very efficient plasmoid-dominated regime, and current layer gets fragmented into a dynamic chain of plasmoids which undergo successive coalescence. Mergers of sufficiently large plasmoids produce secondary perpendicular current sheets, which are also plasmoid-unstable. Collisions of plasmoids with each other and with the upstream magnetic field eject fast-magnetosonic waves, which propagate upstream across the background field and successfully escape from the plasma as electromagnetic waves that fall in the radio band. This model successfully explains many important features of the observed radio emission from Crab and other pulsars with high magnetic field at the light cylinder: phase coincidence with the high-energy emission, nano-second duration (nanoshots), and extreme instantaneous brightness of individual pulses.

We use fully kinetic particle-in-cell simulations with unprecedentedly large transverse box sizes to study particle acceleration in weakly-magnetized mildly relativistic shocks traveling at a velocity $\approx 0.75c$ and a Mach number of 15. We examine both subluminal (quasi-parallel) and superluminal (quasi-perpendicular) magnetic field orientations. We find that quasi-parallel shocks are mediated by a filamentary non-resonant (Bell) instability driven by non-thermal ions, producing magnetic fluctuations on scales comparable to the ion gyro-radius. In quasi-parallel shocks, both electrons and ions are accelerated into non-thermal power-laws whose maximum energy grows linearly with time. The upstream heating of electrons is small, and the two species enter the shock front in rough thermal equilibrium. The shock's structure is complex; the current of reflected non-thermal ions evacuates cavities in the upstream which form filaments of amplified magnetic fields once advected downstream. At late times, $10\%$ of the shock's energy goes into non-thermal protons and $\gtrsim10\%$ into magnetic fields. We find that properly capturing the magnetic turbulence driven by the non-thermal ions is important for properly measuring the energy fraction of non-thermal electrons, $\epsilon_e$. We find $\epsilon_e\sim 5\times10^{-4}$ for quasi-parallel shocks with $v=0.75c$, slightly larger than what was measured in simulations of non-relativistic shocks. In quasi-perpendicular shocks, no non-thermal power-law develops in ions or electrons. The ion acceleration efficiency in quasi-parallel shocks suggests that astrophysical objects that could host mildly relativistic quasi-parallel shocks -- for example, the jets of active galactic nuclei or microquasars -- may be important sources of cosmic rays and their secondaries, such as gamma-rays and neutrinos.

High energy radiation from nonthermal particles accelerated in relativistic magnetic reconnection is thought to be important in many astrophysical systems, ranging from blazar jets and black hole accretion disk coronae to pulsars and magnetar flares. The presence of a substantial density of high energy photons ($>$MeV) in these systems can make two-photon pair production ($\gamma\gamma\to e^-e^+$) an additional source of plasma particles and can affect the radiative properties of these objects. We present the results of novel particle-in-cell simulations that track both the radiated synchrotron photons and the created pairs, with which we study the evolution of a two-dimensional reconnecting current sheet in pair plasma. Synchrotron radiation from accelerated particles in the current sheet produces hot secondary pairs in the upstream which are later advected into the current sheet where they are reaccelerated and produce more photons. In the optically thin regime, when most of the radiation is leaving the upstream unaffected, this process is self-regulating and depends only on the background magnetic field and the optical depth of photons to pair production. The extra plasma loading also affects the properties of reconnection. We study how the inflow of the secondary plasma, with multiplicities up to several hundred, reduces the effective magnetization of the plasma, suppressing the acceleration and thus decreasing the high energy photon spectrum cutoff. This offers an explanation for the weak dependence of the observed gamma-ray cutoff in pulsars on the magnetic field at the light cylinder.

We have performed 2D hybrid simulations of non-relativistic collisionless shocks in the presence of pre-existing energetic particles ("seeds"); such a study applies, for instance, to the re-acceleration of Galactic cosmic rays (CRs) in supernova remnant (SNR) shocks and solar wind energetic particles in heliospheric shocks. Energetic particles can be effectively reflected and accelerated regardless of shock inclination via a process that we call \emph{diffusive shock re-acceleration}. We find that reaccelerated seeds can drive the streaming instability in the shock upstream and produce effective magnetic field amplification. This can eventually trigger the injection of thermal protons even at oblique shocks that ordinarily cannot inject thermal particles. We characterize the current in reflected seeds, finding that it tends to a universal value $J\simeq e n_{cr} v_{sh}$, where $e n_{cr}$ is the seed charge density and $v_{sh}$ is the shock velocity. When applying our results to SNRs, we find that the re-acceleration of Galactic CRs can excite the Bell instability to non-linear levels in less than $\sim 10$ yr, thereby providing a minimum level of magnetic field amplification for any SNR shock. Finally, we discuss the relevance of diffusive shock re-acceleration also for other environments, such as heliospheric shocks, Galactic superbubbles, and clusters of galaxies.

A heat flux in a high-$\beta$ plasma with low collisionality triggers the whistler instability. Quasilinear theory predicts saturation of the instability in a marginal state characterized by a heat flux that is fully controlled by electron scattering off magnetic perturbations. This marginal heat flux does not depend on the temperature gradient and scales as $1/\beta$. We confirm this theoretical prediction by performing numerical particle-in-cell simulations of the instability. We further calculate the saturation level of magnetic perturbations and the electron scattering rate as functions of $\beta$ and the temperature gradient to identify the saturation mechanism as quasilinear. Suppression of the heat flux is caused by oblique whistlers with magnetic-energy density distributed over a wide range of propagation angles. This result can be applied to high-$\beta$ astrophysical plasmas, such as the intracluster medium, where thermal conduction at sharp temperature gradients along magnetic-field lines can be significantly suppressed. We provide a convenient expression for the amount of suppression of the heat flux relative to the classical Spitzer value as a function of the temperature gradient and $\beta$. For a turbulent plasma, the additional independent suppression by the mirror instability is capable of producing large total suppression factors (several tens in galaxy clusters) in regions with strong temperature gradients.

We perform global particle-in-cell simulations of pulsar magnetospheres including pair production, ion extraction from the surface, frame dragging corrections, and high energy photon emission and propagation. In the case of oblique rotators, effects of general relativity increase the fraction of open field lines which support active pair discharge. We find that the plasma density and particle energy flux in the pulsar wind are highly non-uniform with latitude. Significant fraction of the outgoing particle energy flux is carried by energetic ions, which are extracted from the stellar surface. Their energies may extend up to a large fraction of the open field line voltage, making them interesting candidates for ultra-high-energy cosmic rays. We show that pulsar gamma-ray radiation is dominated by synchrotron emission, produced by particles that are energized by relativistic magnetic reconnection close to the Y-point and in the equatorial current sheet. In most cases, calculated light curves contain two strong peaks, in general agreement with Fermi observations. The radiative efficiency decreases with increasing pulsar inclination and increasing efficiency of pair production in the current sheet, explaining the observed scatter in $L_{\gamma}$ vs $\dot{E}$. We find that the high-frequency cutoff in the spectra is regulated by the pair loading of the current sheet. Our findings lay the foundation for quantitative interpretation of Fermi observations of gamma-ray pulsars.

Accreting pulsars power relativistic jets, and display a complex spin phenomenology. These behaviours may be closely related to the large-scale configuration of the star's magnetic field, shaped by its interaction with the surrounding accretion disc. Here we present the first relativistic simulations of the interaction of a pulsar magnetosphere with an accretion flow. Our axisymmetric simulations treat the magnetospheric, or coronal, regions using a resistive extension of force-free electrodynamics. The magnetic field is also evolved inside the disc, which is a defined volume with a specified velocity field and conductivity profile, found using an alpha-disc model. We study a range of disc alpha-parameters, thicknesses, magnetic Prandtl numbers, and inner truncation radii. We find that a large fraction of the magnetic flux in the pulsar's closed zone is opened by the intrusion of the disc, leading to an enhancement of the power extracted by the pulsar wind and the spin-down torque applied to the pulsar. In our simulations, most of the spin-down contribution to the stellar torque acts on open field lines. The efficiency of field-line opening is high in the simulations' long-term quasi-steady states, which implies that a millisecond pulsar's electromagnetic wind could be strong enough to power the observed neutron-star radio jets, and may significantly affect the pulsar's spin evolution.

We study the thermalization, injection, and acceleration of ions with different mass/charge ratios, $A/Z$, in non-relativistic collisionless shocks via hybrid (kinetic ions-fluid electrons) simulations. In general, ions thermalize to a post-shock temperature proportional to $A$. When diffusive shock acceleration is efficient, ions develop a non-thermal tail whose extent scales with $Z$ and whose normalization is enhanced as $(A/Z)^2$, so that incompletely-ionized heavy ions are preferentially accelerated. We discuss how these findings can explain observed heavy-ion enhancements in Galactic cosmic rays.

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

The role of the non-resonant firehose instability in conditions relevant to the precursors of supernova remnant shocks is considered. Using a second order tensor expansion of the Vlasov-Fokker-Planck equation we illustrate the necessary conditions for the firehose to operate. It is found that for very fast shocks, the diffusion approximation predicts that the linear firehose growth rate is marginally faster than its resonant counterpart. Preliminary hybrid MHD-Vlasov-Fokker-Planck simulation results using young supernova relevant parameters are presented.

The interaction of a rotating star's magnetic field with a surrounding plasma disk lies at the heart of many questions posed by neutron stars in X-ray binaries. We consider the opening of stellar magnetic flux due to differential rotation along field lines coupling the star and disk, using a simple model for the disk-opened flux, the torques exerted on the star by the magnetosphere, and the power extracted by the electromagnetic wind. We examine the conditions under which the system enters an equilibrium spin state, in which the accretion torque is instantaneously balanced by the pulsar wind torque alone. For magnetic moments, spin frequencies, and accretion rates relevant to accreting millisecond pulsars, the spin-down torque from this enhanced pulsar wind can be substantially larger than that predicted by existing models of the disk-magnetosphere interaction, and is in principle capable of maintaining spin equilibrium at frequencies less than 1 kHz. We speculate that this mechanism may account for the non-detection of frequency increases during outbursts of SAX J1808.4-3658 and XTE J1814-338, and may be generally responsible for preventing spin-up to sub-millisecond periods. If the pulsar wind is collimated by the surrounding environment, the resulting jet can satisfy the power requirements of the highly relativistic outflows from Cir X-1 and Sco X-1. In this framework, the jet power scales relatively weakly with accretion rate, L_j ~ Mdot^{4/7}, and would be suppressed at high accretion rates only if the stellar magnetic moment is sufficiently low.

Current models of gamma-ray lightcurves in pulsars suffer from large uncertainties on the precise location of particle acceleration and radiation. Here, we present an attempt to alleviate these difficulties by solving for the electromagnetic structure of the oblique magnetosphere, particle acceleration, and the emission of radiation self-consistently, using 3D spherical particle-in-cell simulations. We find that the low-energy radiation is synchro-curvature radiation from the polar-cap regions within the light cylinder. In contrast, the high-energy emission is synchrotron radiation that originates exclusively from the Y-point and the equatorial current sheet where relativistic magnetic reconnection accelerates particles. In most cases, synthetic high-energy lightcurves contain two peaks that form when the current sheet sweeps across the observer's line of sight. We find clear evidence of caustics in the emission pattern from the current sheet. High-obliquity solutions can present up to two additional secondary peaks from energetic particles in the wind region accelerated by the reconnection-induced flow near the current sheet. The high-energy radiative efficiency depends sensitively on the viewing angle, and decreases with increasing pulsar inclination. The high-energy emission is concentrated in the equatorial regions where most of the pulsar spindown is released and dissipated. These results have important implications for the interpretation of gamma-ray pulsar data.

We perform first-principles relativistic particle-in-cell simulations of aligned pulsar magnetosphere. We allow free escape of particles from the surface of a neutron star and continuously populate the magnetosphere with neutral pair plasma to imitate pair production. As pair plasma supply increases, we observe the transition from a charge-separated electrosphere solution with trapped plasma and no spin-down to a solution close to the ideal force-free magnetosphere with electromagnetically-dominated pulsar wind. We calculate the magnetospheric structure, current distribution and spin-down power of the neutron star. We also discuss particle acceleration in the equatorial current sheet.

It has recently been demonstrated that self-consistent particle-in-cell simulations of low-obliquity pulsar magnetospheres in flat spacetime show weak particle acceleration and no pair production near the poles. We investigate the validity of this conclusion in a more realistic spacetime geometry via general-relativistic particle-in-cell simulations of the aligned pulsar magnetospheres with pair formation. We find that the addition of frame-dragging effect makes local current density along the magnetic field larger than the Goldreich-Julian value, which leads to unscreened parallel electric fields and the ignition of a pair cascade. When pair production is active, we observe field oscillations in the open field bundle which could be related to pulsar radio emission. We conclude that general relativistic effects are essential for the existence of pulsar mechanism in low obliquity rotators.

The Wisconsin Plasma Astrophysics Laboratory (WiPAL) is a flexible user facility designed to study a range of astrophysically relevant plasma processes as well as novel geometries that mimic astrophysical systems. A multi-cusp magnetic bucket constructed from strong samarium cobalt permanent magnets now confines a 10 m$^3$, fully ionized, magnetic-field free plasma in a spherical geometry. Plasma parameters of $ T_{e}\approx5$ to $20$ eV and $n_{e}\approx10^{11}$ to $5\times10^{12}$ cm$^{-3}$ provide an ideal testbed for a range of astrophysical experiments including self-exciting dynamos, collisionless magnetic reconnection, jet stability, stellar winds, and more. This article describes the capabilities of WiPAL along with several experiments, in both operating and planning stages, that illustrate the range of possibilities for future users.

We formulate a magnetohydrodynamic-particle-in-cell (MHD-PIC) method for describing the interaction between collisionless cosmic ray (CR) particles and a thermal plasma. The thermal plasma is treated as a fluid, obeying equations of ideal MHD, while CRs are treated as relativistic Lagrangian particles subject to the Lorentz force. Backreaction from CRs to the gas is included in the form of momentum and energy feedback. In addition, we include the electromagnetic feedback due to CR-induced Hall effect that becomes important when the electron-ion drift velocity of the background plasma induced by CRs approaches the Alfv\'en velocity. Our method is applicable on scales much larger than the ion inertial length, bypassing the microscopic scales that must be resolved in conventional PIC methods, while retaining the full kinetic nature of the CRs. We have implemented and tested this method in the Athena MHD code, where the overall scheme is second-order accurate and fully conservative. As a first application, we describe a numerical experiment to study particle acceleration in non-relativistic shocks. Using a simplified prescription for ion injection, we reproduce the shock structure and the CR energy spectra obtained with more self-consistent hybrid-PIC simulations, but at substantially reduced computational cost. We also show that the CR-induced Hall effect reduces the growth rate of the Bell instability and affects the gas dynamics in the vicinity of the shock front. As a step forward, we are able to capture the transition of particle acceleration from non relativistic to relativistic regimes, with momentum spectrum $f(p)\sim p^{-4}$ connecting smoothly through the transition, as expected from the theory of Fermi acceleration.

Collisionless shocks can be produced as a result of strong magnetic fields in a plasma flow, and therefore are common in many astrophysical systems. The Weibel instability is one candidate mechanism for the generation of sufficiently strong fields to create a collisionless shock. Despite their crucial role in astrophysical systems, observation of the magnetic fields produced by Weibel instabilities in experiments has been challenging. Using a proton probe to directly image electromagnetic fields, we present evidence of Weibel-generated magnetic fields that grow in opposing, initially unmagnetized plasma flows from laser-driven laboratory experiments. Three-dimensional particle-in-cell simulations reveal that the instability efficiently extracts energy from the plasma flows, and that the self-generated magnetic energy reaches a few percent of the total energy in the system. This result demonstrates an experimental platform suitable for the investigation of a wide range of astrophysical phenomena, including collisionless shock formation in supernova remnants, large-scale magnetic field amplification, and the radiation signature from gamma-ray bursts.

Rotating neutron stars, or pulsars, are plausibly the source of power behind many astrophysical systems, such as gamma-ray bursts, supernovae, pulsar wind nebulae and supernova remnants. In the past several years, 3D numerical simulations made it possible to compute pulsar spindown luminosity from first principles and revealed that oblique pulsar winds are more powerful than aligned ones. However, what causes this enhanced power output of oblique pulsars is not understood. In this work, using time-dependent 3D magnetohydrodynamic (MHD) and force-free simulations, we show that, contrary to the standard paradigm, the open magnetic flux, which carries the energy away from the pulsar, is laterally non-uniform. We argue that this non-uniformity is the primary reason for the increased luminosity of oblique pulsars. To demonstrate this, we construct simple analytic descriptions of aligned and orthogonal pulsar winds and combine them to obtain an accurate 3D description of the pulsar wind for any obliquity. Our approach describes both the warped magnetospheric current sheet and the smooth variation of pulsar wind properties outside of it. We find that generically the magnetospheric current sheet separates plasmas that move at mildly relativistic velocities relative to each other. This suggests that the magnetospheric reconnection is a type of driven, rather than free, reconnection. The jump in magnetic field components across the current sheet decreases with increasing obliquity, which could be a mechanism that reduces dissipation in near-orthogonal pulsars. Our analytical description of the pulsar wind can be used for constructing models of pulsar gamma-ray emission, pulsar wind nebulae, and magnetar-powered core-collapse gamma-ray bursts and supernovae.

We study diffusive shock acceleration (DSA) of protons and electrons at nonrelativistic, high Mach number, quasiparallel, collisionless shocks by means of self-consistent 1D particle-in-cell simulations. For the first time, both species are found to develop power-law distributions with the universal spectral index $-4$ in momentum space, in agreement with the prediction of DSA. We find that scattering of both protons and electrons is mediated by right-handed circularly polarized waves excited by the current of energetic protons via non-resonant hybrid (Bell) instability. Protons are injected into DSA after a few gyro-cycles of shock drift acceleration (SDA), while electrons are first pre-heated via SDA, then energized via a hybrid acceleration process that involves both SDA and Fermi-like acceleration mediated by Bell waves, before eventual injection into DSA. Using the simulations we can measure the electron/proton ratio in accelerated particles, which is of paramount importance for explaining the cosmic ray fluxes measured at Earth and the multi-wavelength emission of astrophysical objects such as supernova remnants, radio supernovae, and galaxy clusters. We find the normalization of electron power-law is $\lesssim 10^{-2}$ that of the protons for strong nonrelativistic shocks.

The equatorial current sheet in pulsar magnetospheres is often regarded as an ideal site for particle acceleration via relativistic reconnection. Using 2D spherical particle-in-cell simulations, we investigate particle acceleration in the axisymmetric pulsar magnetosphere as a function of the injected plasma multiplicity and magnetization. We observe a clear transition from a highly charge-separated magnetosphere for low plasma injection with little current and spin-down power, to a nearly force-free solution for high plasma multiplicity characterized by a prominent equatorial current sheet and high spin-down power. We find significant magnetic dissipation in the current sheet, up to 30% within 5 light-cylinder radii in the high-multiplicity regime. The simulations unambiguously demonstrate that the dissipated Poynting flux is efficiently channeled to the particles in the sheet, close to the Y-point within about 1-2 light cylinder radii from the star. The mean particle energy in the sheet is given by the upstream plasma magnetization at the light cylinder. The study of particle orbits shows that all energetic particles originate from the boundary layer between the open and the closed field lines. Energetic positrons always stream outward, while high-energy electrons precipitate back towards the star through the sheet and along the separatrices, which may result in auroral-like emission. Our results suggest that the current sheet and the separatrices may be the main source of high-energy radiation in young pulsars.

We use kinetic hybrid simulations (kinetic ions - fluid electrons) to characterize the fraction of ions that are accelerated to non-thermal energies at non-relativistic collisionless shocks. We investigate the properties of the shock discontinuity and show that shocks propagating almost along the background magnetic field (quasi-parallel shocks) reform quasi-periodically on ion cyclotron scales. Ions that impinge on the shock when the discontinuity is the steepest are specularly reflected. This is a necessary condition for being injected, but it is not sufficient. Also by following the trajectories of reflected ions, we calculate the minimum energy needed for injection into diffusive shock acceleration, as a function of the shock inclination. We construct a minimal model that accounts for the ion reflection from quasi-periodic shock barrier, for the fraction of injected ions, and for the ion spectrum throughout the transition from thermal to non-thermal energies. This model captures the physics relevant for ion injection at non-relativistic astrophysical shocks with arbitrary strengths and magnetic inclinations, and represents a crucial ingredient for understanding the diffusive shock acceleration of cosmic rays.

We study the linear stability of weakly magnetized differentially rotating plasmas in both collisionless kinetic theory and Braginskii's theory of collisional, magnetized plasmas. We focus on the very weakly magnetized limit that is important for understanding how astrophysical magnetic fields originate and are amplified at high redshift. We show that the single instability of fluid theory - the magnetorotational instability mediated by magnetic tension - is replaced by two distinct instabilities, one associated with ions and one with electrons. Each of these has a different way of tapping into the free energy of differential rotation. The ion instability is driven by viscous transport of momentum across magnetic field lines due to a finite ion cyclotron frequency (gyroviscosity); the fastest growing modes have wavelengths significantly longer than MHD and Hall MHD predictions. The electron instability is a whistler mode driven unstable by the temperature anisotropy generated by differential rotation; the growth time can be orders of magnitude shorter than the rotation period. The electron instability is an example of a broader class of instabilities that tap into the free energy of differential rotation or shear via the temperature anisotropy they generate. We briefly discuss the application of our results to the stability of planar shear flows and show that such flows are linearly overstable in the presence of fluid gyroviscosity. We also briefly describe the implications of our results for magnetic field amplification in the virialized halos of high redshift galaxies.

We present first-principles relativistic particle-in-cell simulations of the oblique pulsar magnetosphere with pair formation. The magnetosphere starts to form with particles extracted from the surface of the neutron star. These particles are accelerated by surface electric fields and emit photons capable of producing electron-positron pairs. We inject secondary pairs at locations of primary energetic particles, whose energy exceeds the threshold for pair formation. We find solutions that are close to the ideal force-free magnetosphere, with the Y-point and current sheet. Solutions with obliquities $\lt 40^{\circ}$ do not show pair production in the open field line region, because the local current density along magnetic field is below the Goldreich-Julian value. The bulk outflow in these solutions is charge separated, and pair formation happens in the current sheet and return current layer only. Solutions with higher inclinations show pair production in the open field line region, with high multiplicity of the bulk flow and the size of pair-producing region increasing with inclination. We observe the spin-down of the star to be comparable to MHD model predictions. The magnetic dissipation in the current sheet ranges between 20% for the aligned rotator and 3% for the orthogonal rotator. Our results suggest that for low obliquity neutron stars with suppressed pair formation at the light cylinder, the presence of phenomena related to pair activity in the bulk of the polar region, e.g., radio emission, may crucially depend on the physics beyond our simplified model, such as the effects of curved space-time or multipolar surface fields.

We use large hybrid (kinetic protons-fluid electrons) simulations to investigate the transport of energetic particles in self-consistent electromagnetic configurations of collisionless shocks. In previous papers of this series, we showed that ion acceleration may be very efficient (up to $10-20\%$ in energy), and outlined how the streaming of energetic particles amplifies the upstream magnetic field. Here, we measure particle diffusion around shocks with different strengths, finding that the mean free path for pitch-angle scattering of energetic ions is comparable with their gyroradii calculated in the self-generated turbulence. For moderately-strong shocks, magnetic field amplification proceeds in the quasi-linear regime, and particles diffuse according to the self-generated diffusion coefficient, i.e., the scattering rate depends only on the amount of energy in modes with wavelengths comparable with the particle gyroradius. For very strong shocks, instead, the magnetic field is amplified up to non-linear levels, with most of the energy in modes with wavelengths comparable to the gyroradii of highest-energy ions, and energetic particles experience Bohm-like diffusion in the amplified field. We also show how enhanced diffusion facilitates the return of energetic particles to the shock, thereby determining the maximum energy that can be achieved in a given time via diffusive shock acceleration. The parametrization of the diffusion coefficient that we derive can be used to introduce self-consistent microphysics into large-scale models of cosmic ray acceleration in astrophysical sources, such as supernova remnants and clusters of galaxies.

We use large hybrid simulations to study ion acceleration and generation of magnetic turbulence due to the streaming of particles that are self-consistently accelerated at non-relativistic shocks. When acceleration is efficient, we find that the upstream magnetic field is significantly amplified. The total amplification factor is larger than 10 for shocks with Alfv\'enic Mach number $M=100$, and scales with the square root of $M$. The spectral energy density of excited magnetic turbulence is determined by the energy distribution of accelerated particles, and for moderately-strong shocks ($M\lesssim30$) agrees well with the prediction of resonant streaming instability, in the framework of quasilinear theory of diffusive shock acceleration. For $M\gtrsim30$, instead, Bell's non-resonant hybrid (NRH) instability is predicted and found to grow faster than resonant instability. NRH modes are excited far upstream by escaping particles, and initially grow without disrupting the current, their typical wavelengths being much shorter than the current ions' gyroradii. Then, in the nonlinear stage, most unstable modes migrate to larger and larger wavelengths, eventually becoming resonant in wavelength with the driving ions, which start diffuse. Ahead of strong shocks we distinguish two regions, separated by the free-escape boundary: the far upstream, where field amplification is provided by the current of escaping ions via NRH instability, and the shock precursor, where energetic particles are effectively magnetized, and field amplification is provided by the current in diffusing ions. The presented scalings of magnetic field amplification enable the inclusion of self-consistent microphysics into phenomenological models of ion acceleration at non-relativistic shocks.

We use 2D and 3D hybrid (kinetic ions - fluid electrons) simulations to investigate particle acceleration and magnetic field amplification at non-relativistic astrophysical shocks. We show that diffusive shock acceleration operates for quasi-parallel configurations (i.e., when the background magnetic field is almost aligned with the shock normal) and, for large sonic and Alfv\'enic Mach numbers, produces universal power-law spectra proportional to p^(-4), where p is the particle momentum. The maximum energy of accelerated ions increases with time, and it is only limited by finite box size and run time. Acceleration is mainly efficient for parallel and quasi-parallel strong shocks, where 10-20% of the bulk kinetic energy can be converted to energetic particles, and becomes ineffective for quasi-perpendicular shocks. Also, the generation of magnetic turbulence correlates with efficient ion acceleration, and vanishes for quasi-perpendicular configurations. At very oblique shocks, ions can be accelerated via shock drift acceleration, but they only gain a factor of a few in momentum, and their maximum energy does not increase with time. These findings are consistent with the degree of polarization and the morphology of the radio and X-ray synchrotron emission observed, for instance, in the remnant of SN 1006. We also discuss the transition from thermal to non-thermal particles in the ion spectrum (supra-thermal region), and we identify two dynamical signatures peculiar of efficient particle acceleration, namely the formation of an upstream precursor and the alteration of standard shock jump conditions.

In magnetized astrophysical outflows, the dissipation of field energy into particle energy via magnetic reconnection is often invoked to explain the observed non-thermal signatures. By means of two- and three-dimensional particle-in-cell simulations, we investigate anti-parallel reconnection in magnetically-dominated electron-positron plasmas. Our simulations extend to unprecedentedly long temporal and spatial scales, so we can capture the asymptotic state of the system beyond the initial transients, and without any artificial limitation by the boundary conditions. At late times, the reconnection layer is organized into a chain of large magnetic islands connected by thin X-lines. The plasmoid instability further fragments each X-line into a series of smaller islands, separated by X-points. At the X-points, the particles become unmagnetized and they get accelerated along the reconnection electric field. We provide definitive evidence that the late-time particle spectrum integrated over the whole reconnection region is a power-law, whose slope is harder than -2 for magnetizations sigma>10. Efficient particle acceleration to non-thermal energies is a generic by-product of the long-term evolution of relativistic reconnection in both two and three dimensions. In three dimensions, the drift-kink mode corrugates the reconnection layer at early times, but the long-term evolution is controlled by the plasmoid instability, that facilitates efficient particle acceleration, in analogy to the two-dimensional physics. Our findings have important implications for the generation of hard photon spectra in pulsar winds and relativistic astrophysical jets.

The current state of the art in pulsar magnetosphere modeling assumes the force-free limit of magnetospheric plasma. This limit retains only partial information about plasma velocity and neglects plasma inertia and temperature. We carried out time-dependent 3D relativistic magnetohydrodynamic (MHD) simulations of oblique pulsar magnetospheres that improve upon force-free by retaining the full plasma velocity information and capturing plasma heating in strong current layers. We find rather low levels of magnetospheric dissipation, with less than 10% of pulsar spindown energy dissipated within a few light cylinder radii, and the MHD spindown that is consistent with that in force-free. While oblique magnetospheres are qualitatively similar to the rotating split-monopole force-free solution at large radii, we find substantial quantitative differences with the split-monopole, e.g., the luminosity of the pulsar wind is more equatorially concentrated than the split-monopole at high obliquities, and the flow velocity is modified by the emergence of reconnection flow directed into the current sheet.

Supersonic plasma outflows driven by multi-beam, high-energy lasers, such as Omega and NIF, have been and will be used as platforms for a variety of laboratory astrophysics experiments. Here we propose a new way of launching high density and high velocity, plasma jets using multiple intense laser beams in a hollow ring formation. We show that such jets provide a more flexible and versatile platform for future laboratory astrophysics experiments. Using high resolution hydrodynamic simulations, we demonstrate that the collimated jets can achieve much higher density, temperature and velocity when multiple laser beams are focused to form a hollow ring pattern at the target, instead of focused onto a single spot. We carried out simulations with different ring radii and studied their effects on the jet properties. Implications for laboratory collisionless shock experiments are discussed.

We investigate guide-field magnetic reconnection and particle acceleration in relativistic pair plasmas with three-dimensional particle-in-cell (PIC) simulations of a kinetic-scale current sheet in a periodic geometry at low magnetizations. The tearing instability is the dominant mode in the current sheet for all guide field strengths, while the linear kink mode is less important even without guide field. Oblique modes seem to be suppressed entirely. In its nonlinear evolution, the reconnection layer develops a network of interconnected and interacting magnetic flux ropes. As smaller flux ropes merge into larger ones, the reconnection layer evolves toward a three-dimensional, disordered state in which the resulting flux rope segments contain magnetic substructure on plasma skin depth scales. Embedded in the flux ropes, we detect spatially and temporally intermittent sites of dissipation reflected in peaks in the parallel electric field. Magnetic dissipation and particle acceleration persist until the end of the simulations, with simulations with higher magnetization and lower guide field strength exhibiting greater and faster energy conversion and particle energization. At the end of our largest simulation, the particle energy spectrum attains a tail extending to high Lorentz factors that is best modeled with a combination of two additional thermal components. We confirm that the primary energization mechanism is acceleration by the electric field in the X-line region. We discuss the implications of our results for macroscopic reconnection sites, and which of our results may be expected to hold in systems with higher magnetizations.

We used unprecedentedly large 2D and 3D hybrid (kinetic ions - fluid electrons) simulations of non-relativistic collisionless strong shocks in order to investigate the effects of self-consistently accelerated ions on the overall shock dynamics. The current driven by suprathermal particles streaming ahead of the shock excites modes transverse to the background magnetic field. The Lorentz force induced by these self-amplified fields tends to excavate tubular, underdense, magnetic-field-depleted cavities that are advected with the fluid and perturb the shock surface, triggering downstream turbulent motions. These motions further amplify the magnetic field, up to factors of 50-100 in knot-like structures. Once downstream, the cavities tend to be filled by hot plasma plumes that compress and stretch the magnetic fields in elongated filaments; this effect is particularly evident if the shock propagates parallel to the background field. Highly-magnetized knots and filaments may provide explanations for the rapid X-ray variability observed in RX J1713.7-3946 and for the regular pattern of X-ray bright stripes detected in Tycho's supernova remnant.

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.

The magnetosphere of a rotating pulsar naturally develops a current sheet beyond the light cylinder (LC). Magnetic reconnection in this current sheet inevitably dissipates a nontrivial fraction of the pulsar spin-down power within a few LC radii. We develop a basic physical picture of reconnection in this environment and discuss its implications for the observed pulsed gamma-ray emission. We argue that reconnection proceeds in the plasmoid-dominated regime, via an hierarchical chain of multiple secondary islands/flux ropes. The inter-plasmoid reconnection layers are subject to strong synchrotron cooling, leading to significant plasma compression. Using the conditions of pressure balance across these current layers, the balance between the heating by magnetic energy dissipation and synchrotron cooling, and Ampere's law, we obtain simple estimates for key parameters of the layers --- temperature, density, and layer thickness. In the comoving frame of the relativistic pulsar wind just outside of the equatorial current sheet, these basic parameters are uniquely determined by the strength of the reconnecting upstream magnetic field. For the case of the Crab pulsar, we find them to be of order 10 GeV, $10^{13} cm^{-3}$, and 10 cm, respectively. After accounting for the bulk Doppler boosting due to the pulsar wind, the synchrotron and inverse-Compton emission from the reconnecting current sheet can explain the observed pulsed high-energy (GeV) and VHE (~100 GeV) radiation, respectively. Also, we suggest that the rapid relative motions of the secondary plasmoids in the hierarchical chain may contribute to the production of the pulsar radio emission.

By means of two- and three-dimensional particle-in-cell simulations, we investigate the process of driven magnetic reconnection at the termination shock of relativistic striped flows. In pulsar winds and in magnetar-powered relativistic jets, the flow consists of stripes of alternating magnetic field polarity, separated by current sheets of hot plasma. At the wind termination shock, the flow compresses and the alternating fields annihilate by driven magnetic reconnection. Irrespective of the stripe wavelength "lambda" or the wind magnetization "sigma" (in the regime sigma>>1 of magnetically-dominated flows), shock-driven reconnection transfers all the magnetic energy of alternating fields to the particles, whose average Lorentz factor increases by a factor of sigma with respect to the pre-shock value. In the limit lambda/(r_L*sigma)>>1, where r_L is the relativistic Larmor radius in the wind, the post-shock particle spectrum approaches a flat power-law tail with slope around -1.5, populated by particles accelerated by the reconnection electric field. The presence of a current-aligned "guide" magnetic field suppresses the acceleration of particles only when the guide field is stronger than the alternating component. Our findings place important constraints on the models of non-thermal radiation from Pulsar Wind Nebulae and relativistic jets.

The magnetorotational instability (MRI) is a crucial mechanism of angular momentum transport in a variety of astrophysical accretion disks. In systems accreting at well below the Eddington rate, such as the central black hole in the Milky Way (Sgr A*), the rate of Coulomb collisions between particles is very small, making the disk evolve essentially as a collisionless plasma. We present a nonlinear study of the collisionless MRI using first-principles particle-in-cell (PIC) plasma simulations. In this initial study we focus on local two-dimensional (axisymmetric) simulations, deferring more realistic three-dimensional simulations to future work. For simulations with net vertical magnetic flux, the MRI continuously amplifies the magnetic field until the Alfv\'en velocity, v_A, is comparable to the speed of light, c (independent of the initial value of v_A/c). This is consistent with the lack of saturation of MRI channel modes in analogous axisymmetric MHD simulations. The amplification of the magnetic field by the MRI generates a significant pressure anisotropy in the plasma (with the perpendicular pressure being larger than the parallel pressure). We find that this pressure anisotropy in turn excites mirror modes and that the volume averaged pressure anisotropy remains near the threshold for mirror mode excitation. Particle energization is due to both reconnection and viscous heating associated with the pressure anisotropy. Reconnection produces a distinctive power-law component in the energy distribution function of the particles, indicating the likelihood of non-thermal ion and electron acceleration in collisionless accretion disks. This has important implications for interpreting the observed emission -- from the radio to the gamma-rays -- of systems such as Sgr A*.

Magnetospheres of pulsars are thought to be filled with plasma, and variations in plasma supply can affect both pulsar emission properties and spin-down rates. A number of recently discovered "intermittent" pulsars switch between two distinct states: an "on", radio-loud state, and an "off", radio-quiet state. Spin-down rates in the two states differ by a large factor, $\sim 1.5-2.5$, which is not easily understood in the context of current models. In this Letter we present self-consistent numerical solutions of "on" and "off" states of intermittent pulsar magnetospheres. We model the "on" state as a nearly ideal force-free magnetosphere with abundant magnetospheric plasma supply. The lack of radio emission in the "off" state is associated with plasma supply disruption that results in lower plasma density on the open field lines. We model the "off" state using nearly vacuum conditions on the open field lines and nearly ideal force-free conditions on the closed field lines, where plasma can remain trapped even in the absence of pair production. The toroidal advection of plasma in the closed zone in the "off" state causes spin-downs that are a factor of $\sim 2$ higher than vacuum values, and we naturally obtain a range of spin-down ratios between the "on" and "off" states, $\sim 1.2-2.9$, which corresponds to a likely range of pulsar inclination angles of $30{-}90^\circ$. We consider the implications of our model to a number of poorly understood but possibly related pulsar phenomena, including nulling, timing noise, and rotating radio transients.

The current state of the art in the modeling of pulsar magnetospheres invokes either the vacuum or force-free limits for the magnetospheric plasma. Neither of these limits can simultaneously account for both the plasma currents and the accelerating electric fields that are needed to explain the morphology and spectra of high-energy emission from pulsars. To better understand the structure of such magnetospheres, we combine accelerating fields and force-free solutions by considering models of magnetospheres filled with resistive plasma. We formulate Ohm's Law in the minimal velocity fluid frame and construct a family of resistive solutions that smoothly bridges the gap between the vacuum and the force-free magnetosphere solutions. The spin-down luminosity, open field line potential drop, and the fraction of open field lines all transition between the vacuum and force-free values as the plasma conductivity varies from zero to infinity. For fixed inclination angle, we find that the spin-down luminosity depends linearly on the open field line potential drop. We consider the implications of our resistive solutions for the spin down of intermittent pulsars and sub-pulse drift phenomena in radio pulsars.

The relativistic wind of obliquely-rotating pulsars consists of toroidal stripes of opposite magnetic field polarity, separated by current sheets of hot plasma. By means of two- and three-dimensional particle-in-cell simulations, we investigate particle acceleration and magnetic field dissipation at the termination shock of a relativistic striped wind. At the shock, the flow compresses and the alternating fields annihilate by driven magnetic reconnection. Irrespective of the stripe wavelength "lambda" or the wind magnetization "sigma" (in the regime sigma>>1 of magnetically-dominated flows), shock-driven reconnection transfers all the magnetic energy of alternating fields to the particles, whose average Lorentz factor increases by a factor of sigma with respect to the pre-shock value. The shape of the post-shock spectrum depends primarily on the ratio lambda/(r_L*sigma), where r_L is the relativistic Larmor radius in the wind. The spectrum becomes broader as the value of lambda/(r_L*sigma) increases, passing from a relativistic Maxwellian to a flat power-law tail with slope around -1.5, populated by particles accelerated by the reconnection electric field. Close to the equatorial plane of the wind, where the stripes are symmetric, the highest energy particles resulting from magnetic reconnection can escape ahead of the shock, and be injected into a Fermi-like acceleration process. In the post-shock spectrum, they populate a power-law tail with slope around -2.5, that extends beyond the flat component produced by reconnection. Our study suggests that the spectral break between the radio and the optical band in Pulsar Wind Nebulae can be a natural consequence of particle acceleration at the termination shock of striped pulsar winds.

We explore the physics of shock evolution and particle acceleration in non-relativistic collisionless shocks using multidimensional hybrid simulations. We analyze a wide range of physical parameters relevant to the acceleration of cosmic rays (CRs) in astrophysical non-relativistic shock scenarios, such as in supernova remnant (SNR) shocks. We explore the evolution of the shock structure and particle acceleration efficiency as a function of Alfv\'enic Mach number and magnetic field inclination angle $\theta$. We show that there are fundamental differences between high and low Mach number shocks in terms of the electromagnetic turbulence generated in the pre-shock zone and downstream; dominant modes are resonant with the streaming CRs in the low Mach number regime, while both resonant and non-resonant modes are present for high Mach numbers. Energetic power law tails for ions in the downstream plasma can account for up to 15% of the incoming upstream flow energy, distributed over $\sim5%$ of the particles in a power law with slope $-2\pm0.2$ in energy. The energy conversion efficiency (for CRs) peaks at $\theta=15^\circ$ to $30^\circ$ and $M_A=6$, and decreases for higher Mach numbers, down to $\sim2%$ for $M_A=31$. Accelerated particles are produced by Diffusive Shock Acceleration (DSA) and by Shock Drift Acceleration (SDA) mechanisms, with the SDA contribution to the overall energy gain increasing with magnetic inclination. We also present a direct comparison between hybrid and fully kinetic particle-in-cell results at early times; the agreement between the two models justifies the use of hybrid simulations for longer-term shock evolution. In SNR shocks, particle acceleration will be significant for low Mach number quasi-parallel flows ($M_A < 30$, $\theta< 45$). This finding underscores the need for effective magnetic amplification mechanism in SNR shocks.

Electron acceleration to non-thermal, ultra-relativistic energies (~ 10-100 TeV) is revealed by radio and X-ray observations of shocks in young supernova remnants (SNRs). The diffusive shock acceleration (DSA) mechanism is usually invoked to explain this acceleration, but the way in which electrons are initially energized or 'injected' into this acceleration process starting from thermal energies is an unresolved problem. In this paper we study the initial acceleration of electrons in non-relativistic shocks from first principles, using two- and three-dimensional particle-in-cell (PIC) plasma simulations. We systematically explore the space of shock parameters (the Alfv\'enic Mach number, M_A, the shock velocity, v_{sh}, the angle between the upstream magnetic field and the shock normal, theta_{Bn}, and the ion to electron mass ratio, m_i/m_e). We find that significant non-thermal acceleration occurs due to the growth of oblique whistler waves in the foot of quasi-perpendicular shocks. The obtained electron energy distributions show power law tails with spectral indices up to alpha ~ 3-4. The maximum energies of the accelerated particles are consistent with the electron Larmor radii being comparable to that of the ions, indicating potential injection into the subsequent DSA process. This injection mechanism, however, requires the shock waves to have fairly low Alf\'enic Mach numbers, M_A <~ 20, which is consistent with the theoretical conditions for the growth of whistler waves in the shock foot (M_A <~ (m_i/m_e)^{1/2}). Thus, if the whistler mechanism is the only robust electron injection process at work in SNR shocks, then SNRs that display non-thermal emission must have significantly amplified upstream magnetic fields. Such field amplification is likely achieved by the escaping cosmic rays, so electron and proton acceleration in SNR shocks must be interconnected.

The strong variability of magnetic central engines of AGN and GRBs may result in highly intermittent strongly magnetized relativistic outflows. We find a new magnetic acceleration mechanism for such impulsive flows that can be much more effective than the acceleration of steady-state flows. This impulsive acceleration results in kinetic-energy-dominated flows at astrophysically relevant distances from the central source. For a spherical flow, a discrete shell ejected from the source over a time t_0 with Lorentz factor Gamma~1 and initial magnetization sigma_0 = B_0^2/(4 pi rho_0 c^2) >> 1 quickly reaches a typical Lorentz factor Gamma ~ sigma_0^{1/3} and magnetization sigma ~ sigma_0^{2/3} at the distance R_0 ~ ct_0. At this point the magnetized shell of width Delta ~ R_0 in the lab frame loses causal contact with the source and continues to accelerate by spreading significantly in its own rest frame. The expansion is driven by the magnetic pressure gradient and leads to relativistic relative velocities between the front and back of the shell. While the expansion is roughly symmetric in the center of momentum frame, in the lab frame most of the energy and momentum remain in a region (or shell) of width Delta ~ R_0 at the head of the flow. This acceleration proceeds as Gamma ~ (sigma_0 R/R_0)^{1/3} and sigma ~ sigma_0^{2/3} (R/R_0)^{-1/3} until reaching a coasting radius R_c ~ R_0 sigma_0^2 where the kinetic energy becomes dominant: Gamma ~ sigma_0 and sigma ~ 1 at R_c. Then the shell starts coasting and spreading (radially) causing its magnetization to drop as sigma ~ R_c/R at R>R_c. Given the typical variability time-scales of AGN and GRBs, the magnetic acceleration in these sources is a combination of the quasi-steady-state collimation acceleration close to the source and the impulsive acceleration further out. The interaction with the external medium is briefly discussed.

We investigate shock structure and particle acceleration in relativistic magnetized collisionless electron-ion shocks by means of 2.5D particle-in-cell simulations with ion-to-electron mass ratios (m_i/m_e) ranging from 16 to 1000. We explore a range of inclination angles between the pre-shock magnetic field and the shock normal. In "subluminal" shocks, where relativistic particles can escape ahead of the shock along the magnetic field lines, ions are efficiently accelerated via a Fermi-like mechanism. The downstream ion spectrum consists of a relativistic Maxwellian and a high-energy power-law tail, which contains ~5% of ions and ~30% of ion energy. Its slope is -2.1. Upstream electrons enter the shock with lower energy than ions, so they are more strongly tied to the field. As a result, only ~1% of the incoming electrons are Fermi-accelerated at the shock before being advected downstream, where they populate a steep power-law tail (with slope -3.5). For "superluminal" shocks, where relativistic particles cannot outrun the shock along the field, the self-generated turbulence is not strong enough to permit efficient Fermi acceleration, and the ion and electron downstream spectra are consistent with thermal distributions. The incoming electrons are heated up to equipartition with ions, due to strong electromagnetic waves emitted by the shock into the upstream. Thus, efficient electron heating (>15% of the upstream ion energy) is the universal property of relativistic electron-ion shocks, but significant nonthermal acceleration of electrons (>2% by number, >10% by energy, with slope flatter than -2.5) is hard to achieve in magnetized flows and requires weakly magnetized shocks (magnetization <1e-3). These findings place important constraints on the models of AGN jets and Gamma Ray Bursts that invoke particle acceleration in relativistic magnetized electron-ion shocks.

(Abridged) Gamma-ray emission from pulsars has long been modeled using a vacuum dipole field. This approximation ignores changes in the field structure caused by the magnetospheric plasma and strong plasma currents. We present the first results of gamma-ray pulsar light curve modeling using the more realistic field taken from 3D force-free magnetospheric simulations. Having the geometry of the field, we apply several prescriptions for the location of the emission zone, comparing the light curves to observations. We find that the conventional two-pole caustic model fails to produce double-peak pulse profiles, mainly because the size of the polar cap in force-free magnetosphere is larger than the vacuum field polar cap. The conventional outer-gap model is capable of producing only one peak under general conditions, because a large fraction of open field lines does not cross the null charge surface. We propose a novel "separatrix layer" model, where the high-energy emission originates from a thin layer on the open field lines just inside of the separatrix that bounds the open flux tube. The emission from this layer generates two strong caustics on the sky map due to the effect we term "Sky Map Stagnation" (SMS). It is related to the fact that force-free field asymptotically approaches the field of a rotating split monopole, and the photons emitted on such field lines in the outer magnetosphere arrive to the observer in phase. The double-peak light curve is a natural consequence of SMS. We show that most features of the currently available gamma-ray pulsar light curves can be reasonably well reproduced and explained with the sepatratrix model using the force-free field. Association of the emission region with the current sheet will guide more detailed future studies of the magnetospheric acceleration physics.

Current models of pulsar gamma-ray emission use the magnetic field of a rotating dipole in vacuum as a first approximation to the shape of plasma-filled pulsar magnetosphere. In this paper we revisit the question of gamma-ray light-curve formation in pulsars in order to ascertain the robustness of the "two-pole caustic" (TPC) and "outer gap" (OG) models based on the vacuum magnetic field. We point out an inconsistency in the literature on the use of the relativistic aberration formula, where in several works the vacuum field was treated as known in the instantaneously corotating frame, rather than in the laboratory frame. With the corrected formula, we find that the peaks in the light curves predicted from the TPC model using the vacuum field are less sharp. The sharpness of the peaks in the OG model is less affected by this change, but the range of magnetic inclination angles and viewing geometries resulting in double-peaked light curves is reduced. In a realistic magnetosphere, the modification of field structure near the light cylinder due to plasma effects may change the shape of the polar cap and the location of the emission zones. We study the sensitivity of the light curves to different shapes of the polar cap for static and retarded vacuum dipole fields. In particular, we consider polar caps traced by the last open field lines and compare them to circular polar caps. We find that the TPC model is very sensitive to the shape of the polar cap, and a circular polar cap can lead to four peaks of emission. The OG model is less affected by polar cap shapes, but is subject to big uncertainties of applying the vacuum field near the light cylinder. We conclude that deviations from vacuum field can lead to large uncertainties in the pulse shapes, and a more realistic force-free field should be applied to the study of pulsar high energy emission.

(Abridged) X-ray observations of synchrotron rims in supernova remnant (SNR) shocks show evidence of strong magnetic field amplification (a factor of ~100 between the upstream and downstream medium). This amplification may be due to plasma instabilities driven by shock-accelerated cosmic rays (CRs). One candidate is the cosmic ray current-driven (CRCD) instability (Bell 2004), caused by the electric current of large Larmor radii CRs propagating parallel to the upstream magnetic field. Particle-in-cell (PIC) simulations have shown that the back-reaction of the amplified field on CRs would limit the amplification factor of this instability to less than ~10 in galactic SNRs. In this paper, we study the possibility of further amplification driven near shocks by "magnetized" CRs, whose Larmor radii are smaller than the length scale of the field that was previously amplified by the CRCD instability. We find that additional amplification can occur due to a new instability, driven by the CR current perpendicular to the field, which we term the "perpendicular current-driven instability" (PCDI). We derive the growth rate of this instability, and, using PIC simulations, study its non-linear evolution and saturation. We find that PCDI increases the amplification of the field (amplification factor up to ~45, not including the shock compression) and discuss its observational signatures. Our results strengthen the idea of CRs driving a significant part of the magnetic field amplification observed in SNR shocks.

Recent particle-in-cell simulations suggest that a large fraction of the energy dissipated in a relativistic shock is deposited into a Maxwellian distribution of electrons that is connected to the high-energy power-law tail. Here, we explore the observational implications of such a mixed thermal-nonthermal particle distribution for the afterglow and prompt emission of gamma-ray bursts. When the Maxwellian component dominates the energy budget, the afterglow lightcurves show a very steep decline phase followed by a more shallow decay when the characteristic synchrotron frequency crosses the observed band. The steep decay appears in the X-rays at ~100 sec after the burst and is accompanied by a characteristic hard-soft-hard spectral evolution that has been observed in a large number of early afterglows. If internal shocks produce a similar mixed electron distribution, a bump is expected at the synchrotron peak of the nu*f_nu spectrum.

Both radiative and mechanical feedback from Active Galactic Nuclei have been found to be important for the evolution of elliptical galaxies. We compute how a shock may be driven from a central black hole into the gaseous envelope of an elliptical galaxy by mechanical as well as radiative feedback (in the form of nuclear winds) using high resolution 1-D hydrodynamical simulations. We calculate the synchrotron emission from the electron cosmic rays accelerated by the shocks (not the jets) in the central part of elliptical galaxies, and we also study the synchrotron spectrum's evolution using the standard diffusive shock acceleration mechanism, which is routinely applied to the scaled volume case of supernova remnants. We find good agreement quantitatively between the synchrotron radio emission produced via this mechanism with extant observations of elliptical galaxies which are undergoing outbursts. Additionally, we also find that synchrotron optical and X-ray emission can co-exist inside elliptical galaxies during a certain phase of evolution subsequent to central outbursts. In fact, our calculations predict a synchrotron luminosity of $\sim 1.3\times 10^6 L_{\odot}$ at the frequency 5 GHz (radio band), of $\sim 1.1\times 10^6 L_{\odot}$ at $4.3\times10^{14}$ Hz (R band, corresponding to the absolute magnitude -10.4), and of $\sim 1.5\times 10^{7} L_{\odot}$ at $2.4\times10^{17}$ Hz (soft X-ray, 0.5 -- 2.0 keV band).

Gamma-ray emission from pulsars is thought to arise from accelerating regions in pulsar's outer magnetosphere. The shape of the light curves is thus sensitive to the details of the magnetic geometry of the magnetosphere. In this work, we show the first calculations of light curves from the more realistic force-free field under the framework of conventional emission models. We compare the properties of gamma-ray emission between the commonly used vacuum dipole magnetic field and the new force-free field. We discuss the role of the polar cap shape and aberration effect on the appearance of the light curves as well as the formation of caustics on the sky map. With the force-free field, the double-peak pulse profile is best reproduced if the emission zone lies in a thin layer just outside the current sheet, and the peaks are mainly contributed from regions near the light cylinder. The conventional two-pole caustic model can produce up to four peaks in general, while the conventional outer-gap model can normally produce only one peak. These results will be useful for interpreting Fermi telescope observations.

We present a linear analysis of inviscid, incompressible, magnetohydrodynamic (MHD) shallow water systems. In spherical geometry, a generic property of such systems is the existence of five wave modes. Three of them (two magneto-Poincare modes and one magneto-Rossby mode) are previously known. The other two wave modes are strongly influenced by the magnetic field and rotation, and have substantially lower angular frequencies; as such, we term them "magnetostrophic modes". We obtain analytical functions for the velocity, height and magnetic field perturbations in the limit that the magnitude of the MHD analogue of Lamb's parameter is large. On a sphere, the magnetostrophic modes reside near the poles, while the other modes are equatorially confined. Magnetostrophic modes may be an ingredient in explaining the frequency drifts observed in Type I X-ray bursts from neutron stars.

We extract synthetic photon spectra from first-principles particle-in-cell simulations of relativistic shocks propagating in unmagnetized pair plasmas. The two basic ingredients for the radiation, namely accelerated particles and magnetic fields, are produced self-consistently as part of the shock evolution. We use the method of Hededal & Nordlund (2005) and compute the photon spectrum via Fourier transform of the electric far-field from a large number of particles, sampled directly from the simulation. We find that the spectrum from relativistic collisionless shocks is entirely consistent with synchrotron radiation in the magnetic fields generated by Weibel instability. We can recover the so-called "jitter'' regime only if we artificially reduce the strength of the electromagnetic fields, such that the wiggler parameter K = qB lambda/mc^2 becomes much smaller than unity ("B" and "lambda" are the strength and scale of the magnetic turbulence, respectively). These findings may place constraints on the origin of non-thermal emission in astrophysics, especially for the interpretation of the hard (harder than synchrotron) low-frequency spectrum of Gamma-Ray Bursts.

We discuss how the electrons in relativistic GRB shocks can reach near-equipartition in energy with the protons. We emphasize the non-Fermi origin of such acceleration. We argue that the dynamics of the electrons in the foreshock region and at the shock front plays an important role. We also demonstrate that PIC simulations can directly probe this physics in the regimes relevant to GRBs.

We study relativistic unmagnetized collisionless shocks using unprecedentedly large particle-in-cell simulations of two-dimensional pair plasma. High energy particles accelerated by the shock are found to drive magnetic field evolution on a timescale >10^4 plasma times. Progressively stronger magnetic fields are generated on larger scales in a growing region around the shock. Shock-generated magnetic fields and accelerated particles carry >1% and >10% of the downstream energy flux, respectively. Our results suggest limits on the magnetization of relativistic astrophysical flows.

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).

We investigate shock structure and particle acceleration in relativistic magnetized collisionless pair shocks by means of 2.5D and 3D particle-in-cell simulations. We explore a range of inclination angles between the pre-shock magnetic field and the shock normal. We find that only magnetic inclinations corresponding to "subluminal" shocks, where relativistic particles following the magnetic field can escape ahead of the shock, lead to particle acceleration. The downstream spectrum in such shocks consists of a relativistic Maxwellian and a high-energy power-law tail with exponential cutoff. For increasing magnetic inclination in the subluminal range, the high-energy tail accounts for an increasing fraction of particles (from ~1% to ~2%) and energy (from ~4% to ~12%). The spectral index of the power law increases with angle from -2.8+-0.1 to -2.3+-0.1. Particle energization is driven by the Diffusive Shock Acceleration process for nearly parallel shocks, and switches to Shock-Drift Acceleration for larger subluminal inclinations. For "superluminal" shocks, the downstream particle spectrum does not show any significant suprathermal tail. As seen from the upstream frame, efficient acceleration in relativistic (Lorentz factor gamma0 > 5) magnetized (sigma > 0.03) flows exists only for a very small range of magnetic inclination angles (< 34/gamma0 degrees), so relativistic astrophysical pair shocks have to be either nearly parallel or weakly magnetized to generate nonthermal particles. These findings place constraints on the models of AGN jets, Pulsar Wind Nebulae and Gamma Ray Bursts that invoke particle acceleration in relativistic magnetized shocks. (Abridged)

The cosmic ray current-driven (CRCD) instability, predicted by Bell (2004), consists of non-resonant, growing plasma waves driven by the electric current of cosmic rays (CRs) that stream along the magnetic field ahead of both relativistic and non-relativistic shocks. Combining an analytic, kinetic model with one-, two-, and three-dimensional particle-in-cell simulations, we confirm the existence of this instability in the kinetic regime and determine its saturation mechanisms. In the linear regime, we show that, if the background plasma is well magnetized, the CRCD waves grow exponentially at the rates and wavelengths predicted by the analytic dispersion relation. The magnetization condition implies that the growth rate of the instability is much smaller than the ion cyclotron frequency. As the instability becomes non-linear, significant turbulence forms in the plasma. This turbulence reduces the growth rate of the field and damps the shortest wavelength modes, making the dominant wavelength, \lambda_d, grow proportional to the square of the field. At constant CR current, we find that plasma acceleration along the motion of CRs saturates the instability at the magnetic field level such that v_A ~ v_{d,cr}, where v_A is the Alfven velocity in the amplified field, and v_{d,cr} is the drift velocity of CRs. The instability can also saturate earlier if CRs get strongly deflected by the amplified field, which happens when their Larmor radii get close to \lambda_d. We apply these results to the case of CRs in the upstream medium of supernova remnants. Considering only the most energetic CRs that escape from the shock, we obtain that the field amplification factor of ~10 can be reached. This confirms the CRCD instability as a potentially important component of magnetic amplification process in astrophysical shocks.

We address the question of whether numerical particle-in-cell (PIC) simulations and laboratory laser-plasma experiments can (or will be able to, in the near future) model realistic gamma-ray burst (GRB) shocks. For this, we compare the radiative cooling time, t_cool, of relativistic electrons in the shock magnetic fields to the microscopic dynamical time of collisionless relativistic shocks -- the inverse plasma frequency of protons, omega_pp^{-1}. We obtain that for t_cool*omega_pp^{-1}\lesssim ~few hundred, the electrons cool efficiently at or near the shock jump and are capable of emitiing away a large fraction of the shock energy. Such shocks are well-resolved in existing PIC simulations; therefore, the microscopic structure can be studied in detail. Since most of the emission in such shocks would be coming from the vicinity of the shock, the spectral power of the emitted radiation can be directly obtained from finite-length simulations and compared with observational data. Such radiative shocks correspond to the internal baryon-dominated GRB shocks for the conventional range of ejecta parameters. Fermi acceleration of electrons in such shocks is limited by electron cooling, hence the emitted spectrum should be lacking a non-thermal tail, whereas its peak likely falls in the multi-MeV range. Incidentally, the conditions in internal shocks are almost identical to those in laser-produced plasmas; thus, such GRB-like plasmas can be created and studied in laboratory experiments using the presently available Petawatt-scale laser facilities. An analysis of the external shocks shows that only the highly relativistic shocks, corresponding to the extremely early afterglow phase, can have efficient electron cooling in the shock transition. We emphasize the importance of radiative PIC simulations for further studies.

We present evidence that relativistic shocks propagating in unmagnetized plasmas can self-consistently accelerate particles. We use long-term two-dimensional particle-in-cell simulations to study the well-developed shock structure in unmagnetized pair plasma. The particle spectrum downstream of such a shock consists of two components: a relativistic Maxwellian, with characteristic temperature set by the upstream kinetic energy of the flow, and a high-energy tail, extending to energies >100 times that of the thermal peak. This tail is best fitted as a power law in energy with index -2.4+-0.1, modified by an exponential cutoff. The cutoff moves to higher energies with time of the simulation, leaving a larger power law range. The number of particles in the tail is ~1% of the downstream population, and they carry ~10% of the kinetic energy in the downstream. Upon investigation of the trajectories of particles in the tail, we find that the energy gains occur as particles bounce between the upstream and downstream regions in the magnetic fields generated by the Weibel instability. We compare this mechanism to the first order Fermi acceleration, and set a lower limit on the efficiency of shock acceleration process.

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.

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)

Relativistic collisionless shocks in electron-ion plasma are thought to occur in the afterglow phase of Gamma-Ray Bursts (GRBs), and in other environments where relativistic flows interact with the interstellar medium. A particular regime of shocks in an unmagnetized plasma has generated much interest for GRB applications. In this paper we present ab-initio particle-in-cell simulations of unmagnetized relativistic electron-ion shocks. Using long-term 2.5-dimensional simulations with ion-electron mass ratios from 16 to 1000 we resolve the shock formation and reach a steady-state shock structure beyond the initial transient. We find that even at high ion-electron mass ratios initially unmagnetized shocks can be effectively mediated by the ion Weibel instability with a typical shock thickness of ~50 ion skin-depths. Upstream of the shock the interaction with merging ion current filaments heats the electron component, so that the postshock flow achieves near equipartition between the ions and electrons, with electron temperature reaching 50% of the ion temperature. This energy exchange helps to explain the large electron energy fraction inferred from GRB afterglow observations.

Pulsar spindown forms a reliable yet enigmatic prototype for the energy loss processes in many astrophysical objects including accretion disks and back holes. In this paper we review the physics of pulsar magnetospheres, concentrating on recent developments in force-free modeling of the magnetospheric structure. In particular, we discuss a new method for solving the equations of time-dependent force-free relativistic MHD in application to pulsars. This method allows to dynamically study the formation of the magnetosphere and its response to perturbations, opening a qualitatively new window on pulsar phenomena. Applications of the method to other magnetized rotators, such as magnetars and accretion disks, are also discussed.

We discuss 3D simulations of relativistic collisionless shocks in electron-positron pair plasmas using the particle-in-cell (PIC) method. The shock structure is mainly controlled by the shock's magnetization ("sigma" parameter). We demonstrate how the structure of the shock varies as a function of sigma for perpendicular shocks. At low magnetizations the shock is mediated mainly by the Weibel instability which generates transient magnetic fields that can exceed the initial field. At larger magnetizations the shock is dominated by magnetic reflections. We demonstrate where the transition occurs and argue that it is impossible to have very low magnetization collisionless shocks in nature (in more than one spatial dimension). We further discuss the acceleration properties of these shocks, and show that higher magnetization perpendicular shocks do not efficiently accelerate nonthermal particles in 3D. Among other astrophysical applications, this may pose a restriction on the structure and composition of gamma-ray bursts and pulsar wind outflows.

Magnetospheres of many astrophysical objects can be accurately described by the low-inertia (or "force-free") limit of MHD. We present a new numerical method for solution of equations of force-free relativistic MHD based on the finite-difference time-domain (FDTD) approach with a prescription for handling spontaneous formation of current sheets. We use this method to study the time-dependent evolution of pulsar magnetospheres in both aligned and oblique magnetic geometries. For the aligned rotator we confirm the general properties of the time-independent solution of Contopoulos et al. (1999). For the oblique rotator we present the 3D structure of the magnetosphere and compute, for the first time, the spindown power of pulsars as a function of inclination of the magnetic axis. We find the pulsar spindown luminosity to be L = (mu^2 Omega^4/c^3) (1+ sin^2(alpha)) for a star with the dipole moment "mu", rotation frequency "Omega", and magnetic inclination angle "alpha". We also discuss the effects of current sheet resistivity and reconnection on the structure and evolution of the magnetosphere.

We address the issue of electromagnetic pulsar spindown by combining our experience from the two limiting idealized cases which have been studied in great extent in the past: that of an aligned rotator where ideal MHD conditions apply, and that of a misaligned rotator in vacuum. We construct a spindown formula that takes into account the misalignment of the magnetic and rotation axes, and the magnetospheric particle acceleration gaps. We show that near the death line aligned rotators spin down much slower than orthogonal ones. In order to test this approach, we use a simple Monte Carlo method to simulate the evolution of pulsars and find a good fit to the observed pulsar distribution in the P-Pdot diagram without invoking magnetic field decay. Our model may also account for individual pulsars spinning down with braking index n < 3, by allowing the corotating part of the magnetosphere to end inside the light cylinder. We discuss the role of magnetic reconnection in determining the pulsar braking index. We show, however, that n ~ 3 remains a good approximation for the pulsar population as a whole. Moreover, we predict that pulsars near the death line have braking index values n > 3, and that the older pulsar population has preferentially smaller magnetic inclination angles. We discuss possible signatures of such alignment in the existing pulsar data.

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).

A 2.5-dimensional particle-in-cell code is used to investigate the propagation of a large-amplitude, superluminal, nearly transverse electromagnetic (TEM) wave in a relativistically streaming electron-positron plasma with and without a shock. In the freestreaming, unshocked case, the analytic TEM dispersion relation is verified, and the streaming is shown to stabilize the wave against parametric instabilities. In the confined, shocked case, the wave induces strong, coherent particle oscillations, heats the plasma, and modifies the shock density profile via ponderomotive effects. The wave decays over $\gtrsim 10^2$ skin depths; the decay length scale depends primarily on the ratio between the wave frequency and the effective plasma frequency, and on the wave amplitude. The results are applied to the termination shock of the Crab pulsar wind, where the decay length-scale (at least 0.05") might be comparable to the thickness of filamentary, variable substructure observed in the optical and X-ray wisps and knots.

It is generally accepted that magnetic fields generated in the nonlinear development of the transverse Weibel instability provide effective collisionality in unmagnetized collisionless shocks. Recently, extensive two and three dimensional simulations improved our understanding of the growth and saturation of the instability in colliding plasma shells. However, the steady-state structure of the shock wave transition layers remains poorly understood. We use basic physical considerations and order-of-magnitude arguments to study the steady state structure in relativistic unmagnetized collisionless shocks in pair plasmas. The shock contains an electrostatic layer resulting from the formation of stationary, magnetically-focused current filaments. The filaments form where the cold upstream plasma and the counterstreaming thermal plasma interpenetrate. The filaments are not entirely neutral and strong electrostatic fields are present. Most of the downstream particles cannot cross this layer into the upstream because they are trapped by the electrostatic field. We identify the critical location in the shock transition layer where the electromagnetic field ceases to be static. At this location, the degree of charge separation in the filaments reaches a maximum value, the current inside the filaments comes close to the Alfv\'en limit, and the phase space distribution function starts to isotropize. We argue that the radius of the current filaments upstream of the critical location is about twice the upstream plasma skin depth. Finally, we show that some downstream particles cross the electrostatic layer and run ahead of the shock into the preshock medium without causing instability. These particles may play an important role in particle acceleration.

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 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 review the theoretical understanding of the global structure of pulsar magnetospheres concentrating on recent progress in force-free electrodynamics and first-principles simulations of magnetospheres.

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 analyze the global hydrodynamic flow in the ocean of an accreting, rapidly rotating, non-magnetic neutron star in an LMXB during a type I X-ray burst. Our analysis takes into account the rapid rotation of the star and the lift-up of the burning ocean during the burst. We find a new regime for spreading of a nuclear burning front, where the flame is carried along a coherent shear flow across the front. If turbulent viscosity is weak, the speed of flame propagation is ~20 km/s, while, if turbulent viscosity is dynamically important, the flame speed increases, and reaches the maximum value, ~300 km/s, when the eddy overturn frequency is comparable to the Coriolis parameter. We show that, due to rotationally reduced gravity, the thermonuclear runaway is likely to begin on the equator. The equatorial belt is ignited first, and the flame then propagates from the equator to the poles. Inhomogeneous cooling (equator first, poles second) drives strong zonal currents which may be unstable to formation of Jupiter-type vortices; we conjecture that these vortices are responsible for modulation of X-ray flux in the tails of some bursts. We consider the effect of strong zonal currents on the frequency of modulation of the X-ray flux and show that the large values of the frequency drifts observed in some bursts can be accounted for within our model combined with the model of homogeneous radial expansion. Also, if inhomogeneities are trapped in the forward zonal flows around the propagating burning front, chirps with large frequency ranges (~25-500 Hz) may be detectable during the burst rise. We argue that an MHD dynamo within the burning front can generate a small-scale magnetic field, which may enforce vertically rigid flow in the front's wake and can explain the coherence of oscillations in the burst tail.