Normalized to: Spitkovsky, A.

**85** article(s) in total. **150** co-authors, from **1** to **12** common article(s). Median position in authors list is **2,0**.

[1] oai:arXiv.org:2001.02236 [pdf] - 2027943

On the origin of pulsar radio emission

Comments: 9 pages, 4 figures

Submitted: **2020-01-07**

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.

[2] oai:arXiv.org:1911.00059 [pdf] - 2042240

Filling the Magnetospheres of Weak Pulsars

Comments: 9 pages, 6 figures, submitted to ApJ

Submitted: **2019-10-31**

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.

[3] oai:arXiv.org:1908.07890 [pdf] - 1947254

Electron acceleration in non-relativistic quasi-perpendicular
collisionless shocks

Comments: 5 pages, 4 figures, submitted to PRL

Submitted: **2019-08-21**

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.

[4] oai:arXiv.org:1811.01951 [pdf] - 1953354

On the Growth and Saturation of the Gyroresonant Streaming Instabilities

Comments: 30 pages, 19 figures. Accepted to ApJ

Submitted: **2018-11-05**, last modified: **2019-08-14**

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.

[5] oai:arXiv.org:1903.05328 [pdf] - 1919967

Extreme Plasma Astrophysics

Comments: A science white paper submitted to the Astro-2020 and Plasma-2020
Decadal Surveys. 7 pages including cover page and references. Paper updated
in late March 2019 to include a several additional co-authors and references,
and a few small changes

Submitted: **2019-03-13**, last modified: **2019-07-20**

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.

[6] oai:arXiv.org:1906.03297 [pdf] - 1925086

Relativistic Magnetic Reconnection in Electron-Positron-Proton Plasmas:
Implications for Jets of Active Galactic Nuclei

Comments: 25 pages, 23 figures, 4 tables. Accepted for publication in ApJ

Submitted: **2019-06-07**

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.

[7] oai:arXiv.org:1901.02834 [pdf] - 1899839

Black hole magnetosphere with small scale flux tubes--II. Stability and
dynamics

Comments: 16 pages, 17 figures. Accepted for publication in MNRAS

Submitted: **2019-01-09**, last modified: **2019-06-05**

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.

[8] oai:arXiv.org:1905.13735 [pdf] - 1966784

Three-dimensional instability of flame fronts in type I X-ray bursts

Comments:

Submitted: **2019-05-31**

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.

[9] oai:arXiv.org:1903.04080 [pdf] - 1846035

[Plasma 2020 Decadal] The Material Properties of Weakly Collisional,
High-Beta Plasmas

Comments: 6 pages, submitted for the Plasma 2020 Decadal Survey

Submitted: **2019-03-10**

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.

[10] oai:arXiv.org:1902.07730 [pdf] - 1875332

Pulsar Radio Emission Mechanism: Radio Nanoshots as a Low Frequency
Afterglow of Relativistic Magnetic Reconnection

Comments: 7 pages, 4 figures

Submitted: **2019-02-20**

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.

[11] oai:arXiv.org:1809.10772 [pdf] - 1894167

Effects of synchrotron cooling and pair production on collisionless
relativistic reconnection

Comments: 22 pages, 15 figures, submitted to ApJ

Submitted: **2018-09-27**

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.

[12] oai:arXiv.org:1809.10809 [pdf] - 1863853

Kinetic simulations of mildly relativistic shocks I: particle
acceleration in high Mach number shocks

Comments: Submitted to MNRAS. 14 pages, 18 figures

Submitted: **2018-09-27**

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.

[13] oai:arXiv.org:1801.01510 [pdf] - 1716975

Diffusive Shock Re-Acceleration

Comments: 28 pages, 13 figures, Accepted for publication in JPP

Submitted: **2018-01-04**, last modified: **2018-04-30**

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.

[14] oai:arXiv.org:1711.11462 [pdf] - 1716963

Self-inhibiting thermal conduction in high-beta, whistler-unstable
plasma

Comments: accepted to JPP

Submitted: **2017-11-30**, last modified: **2018-04-17**

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.

[15] oai:arXiv.org:1707.04323 [pdf] - 1656102

Ab-Initio Pulsar Magnetosphere: Particle acceleration in Oblique
Rotators and High-energy Emission Modeling

Comments: 11 pages, 8 figures, submitted to ApJ

Submitted: **2017-07-13**

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.

[16] oai:arXiv.org:1608.04159 [pdf] - 1580143

Simulations of the magnetospheres of accreting millisecond pulsars

Comments: 15 pages, 12 figures, published in MNRAS

Submitted: **2016-08-14**, last modified: **2017-06-19**

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.

[17] oai:arXiv.org:1704.08252 [pdf] - 1582888

Chemical Enhancements in Shock-accelerated Particles: Ab-initio
Simulations

Comments: 6 pages, 4 figures, submitted to PRL

Submitted: **2017-04-26**

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.

[18] oai:arXiv.org:1702.01789 [pdf] - 1535477

Preferential Heating and Acceleration of Heavy Ions in Impulsive Solar
Flares

Comments: 6 pages, 7 figures

Submitted: **2017-02-02**

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.

[19] oai:arXiv.org:1610.08302 [pdf] - 1580468

On the Cosmic Ray Driven Firehose Instability

Comments: Contribution to the 6th International Symposium on High Energy
Gamma-Ray Astronomy (Gamma2016), Heidelberg, Germany. To be published in the
AIP Conference Proceedings

Submitted: **2016-10-26**

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.

[20] oai:arXiv.org:1507.08627 [pdf] - 1437558

Torque Enhancement, Spin Equilibrium, and Jet Power from Disk-Induced
Opening of Pulsar Magnetic Fields

Comments: 14 pages, 5 figures, published in ApJ

Submitted: **2015-07-30**, last modified: **2016-07-12**

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.

[21] oai:arXiv.org:1511.01785 [pdf] - 1358948

Modeling high-energy pulsar lightcurves from first principles

Comments: 14 pages, 11 figures, Accepted for publication in MNRAS

Submitted: **2015-11-05**, last modified: **2016-01-13**

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.

[22] oai:arXiv.org:1312.4970 [pdf] - 1289763

Ab-initio pulsar magnetosphere: three-dimensional particle-in-cell
simulations of axisymmetric pulsars

Comments: 6 pages, 5 figures, published in ApJ Letters

Submitted: **2013-12-17**, last modified: **2015-10-06**

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.

[23] oai:arXiv.org:1510.01734 [pdf] - 1327499

Ab-initio pulsar magnetosphere: the role of general relativity

Comments: 5 pages, 4 figure, submitted to ApJLetters

Submitted: **2015-10-06**

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.

[24] oai:arXiv.org:1506.07195 [pdf] - 1260828

The Wisconsin Plasma Astrophysics Laboratory

Comments: 21 pages, 12 figures, 2 tables

Submitted: **2015-06-23**, last modified: **2015-08-04**

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.

[25] oai:arXiv.org:1412.1087 [pdf] - 1263872

Magnetohydrodynamic-Particle-in-Cell Method for Coupling Cosmic Rays
with a Thermal Plasma: Application to Non-relativistic Shocks

Comments: 24 pages, 15 figures, accepted for publication in ApJ

Submitted: **2014-12-02**, last modified: **2015-06-24**

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.

[26] oai:arXiv.org:1310.3337 [pdf] - 968511

Observation of magnetic field generation via the Weibel instability in
interpenetrating plasma flows

Comments:

Submitted: **2013-10-12**, last modified: **2015-03-19**

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.

[27] oai:arXiv.org:1503.01467 [pdf] - 1378555

Three-dimensional Analytical Description of Magnetised Winds from
Oblique Pulsars

Comments: Submitted to MNRAS, comments welcome. 12 pages, 16 figures, uses
mn2e.cls

Submitted: **2015-03-04**

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.

[28] oai:arXiv.org:1412.0672 [pdf] - 941227

Simultaneous Acceleration of Protons and Electrons at Nonrelativistic
Quasiparallel Collisionless Shocks

Comments: 5 pages, 4 figures, accepted to PRL (some updates)

Submitted: **2014-12-01**, last modified: **2015-01-28**

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.

[29] oai:arXiv.org:1410.3757 [pdf] - 1222563

Particle acceleration in axisymmetric pulsar current sheets

Comments: 16 pages, 17 figures, Accepted for publication in MNRAS

Submitted: **2014-10-14**, last modified: **2015-01-07**

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.

[30] oai:arXiv.org:1409.8291 [pdf] - 1222304

Simulations and Theory of Ion Injection at Non-relativistic
Collisionless Shocks

Comments: 6 pages, 5 figures, minor changes, accepted to ApJ Letters

Submitted: **2014-09-29**, last modified: **2014-12-04**

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.

[31] oai:arXiv.org:1412.1097 [pdf] - 1223446

Linear Instabilities Driven by Differential Rotation in Very Weakly
Magnetized Plasmas

Comments: 15 pages, 6 Figures. MNRAS in press

Submitted: **2014-12-02**

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.

[32] oai:arXiv.org:1412.0673 [pdf] - 1223425

Ab-initio pulsar magnetosphere: three-dimensional particle-in-cell
simulations of oblique pulsars

Comments: 7 pages, 5 figures, submitted to ApJLetters

Submitted: **2014-12-01**

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.

[33] oai:arXiv.org:1401.7679 [pdf] - 1203068

Simulations of Ion Acceleration at Non-relativistic Shocks. II. Magnetic
Field Amplification

Comments: 13 pages, 9 figures, Minor changes to reflect the version accepted
for publication in ApJ

Submitted: **2014-01-29**, last modified: **2014-10-07**

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.

[34] oai:arXiv.org:1310.2943 [pdf] - 1179873

Simulations of Ion Acceleration at Non-relativistic Shocks. I.
Acceleration Efficiency

Comments: 21 pages, 14 figures, Minor changes reflecting the version accepted
to ApJ

Submitted: **2013-10-10**, last modified: **2014-10-07**

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.

[35] oai:arXiv.org:1407.2261 [pdf] - 1215510

Simulations of Ion Acceleration at Non-relativistic Shocks. III.
Particle Diffusion

Comments: 8 pages, 7 figures, Minor changes reflecting the version accepted to
ApJ

Submitted: **2014-07-08**, last modified: **2014-10-07**

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.

[36] oai:arXiv.org:1401.5471 [pdf] - 1202911

Relativistic Reconnection: an Efficient Source of Non-Thermal Particles

Comments: 6 pages, 5 figures, ApJL accepted, movies available at
https://www.cfa.harvard.edu/~lsironi/Site/sigma10.no.guide.field/

Submitted: **2014-01-21**

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.

[37] oai:arXiv.org:1211.2803 [pdf] - 1157705

Time-Dependent 3D Magnetohydrodynamic Pulsar Magnetospheres: Oblique
Rotators

Comments: 5 pages, 3 figures, MNRAS, accepted. Movies are available at
http://youtu.be/cjua6XhhNL4 and http://youtu.be/dUR2Lx1JGRM and can be also
downloaded from the ancillary files section

Submitted: **2012-11-12**, last modified: **2013-06-04**

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.

[38] oai:arXiv.org:1209.2150 [pdf] - 651812

Increase of the Density, Temperature and Velocity of Plasma Jets driven
by a Ring of High Energy Laser Beams

Comments: 5 pages, 4 figures, Accepted to HEDP

Submitted: **2012-09-10**, last modified: **2013-04-15**

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.

[39] oai:arXiv.org:1208.0849 [pdf] - 1150512

A Flux Rope Network and Particle Acceleration in Three Dimensional
Relativistic Magnetic Reconnection

Comments: 20 pages, 11 figures

Submitted: **2012-08-03**, last modified: **2013-02-28**

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.

[40] oai:arXiv.org:1211.6765 [pdf] - 1158138

Cosmic-ray-induced filamentation instability in collisionless shocks

Comments: 8 pages, 5 figures (larger than in v1). Minor changes reflecting the
version accepted in ApJL

Submitted: **2012-11-28**, last modified: **2013-02-13**

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.

[41] oai:arXiv.org:1301.5333 [pdf] - 1159220

The Maximum Energy of Accelerated Particles in Relativistic
Collisionless Shocks

Comments: 23 pages, 15 figures, submitted to ApJ

Submitted: **2013-01-22**

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.

[42] oai:arXiv.org:1210.3346 [pdf] - 1152056

Physical Conditions in the Reconnection Layer in Pulsar Magnetospheres

Comments: 20 pages, 2 figures

Submitted: **2012-10-11**

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.

[43] oai:arXiv.org:1208.4998 [pdf] - 1150882

Particle-in-cell simulations of shock-driven reconnection in
relativistic striped winds

Comments: 25 pages, 14 figures, movies available at
https://www.cfa.harvard.edu/~lsironi/sironi_movies.tar ; in press, special
issue of Computational Science and Discovery on selected research from the
22nd International Conference on Numerical Simulation of Plasmas

Submitted: **2012-08-24**

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.

[44] oai:arXiv.org:1201.6407 [pdf] - 1116267

Local 2D Particle-in-cell simulations of the collisionless MRI

Comments: 21 pages, 17 figures

Submitted: **2012-01-30**

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

[45] oai:arXiv.org:1201.2182 [pdf] - 1092864

On the Spin-Down of Intermittent Pulsars

Comments: 6 pages, 4 figures, submitted to ApJ Letters

Submitted: **2012-01-10**

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.

[46] oai:arXiv.org:1107.0979 [pdf] - 1077784

Resistive Solutions for Pulsar Magnetospheres

Comments: 13 pages, 6 figures, final version, accepted to ApJ

Submitted: **2011-07-05**, last modified: **2011-11-17**

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.

[47] oai:arXiv.org:1107.0977 [pdf] - 1077783

Acceleration of Particles at the Termination Shock of a Relativistic
Striped Wind

Comments: 23 pages, 16 figures, submitted to ApJ

Submitted: **2011-07-05**

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.

[48] oai:arXiv.org:1107.0762 [pdf] - 1077755

Ion acceleration in non-relativistic astrophysical shocks

Comments:

Submitted: **2011-07-04**

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.

[49] oai:arXiv.org:1009.3319 [pdf] - 1040647

Electron Injection by Whistler Waves in Non-relativistic Shocks

Comments: 15 pages, 15 figures, Re-submission of version accepted to ApJ

Submitted: **2010-09-16**, last modified: **2011-03-23**

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.

[50] oai:arXiv.org:1004.0959 [pdf] - 1026148

Impulsive Acceleration of Strongly Magnetized Relativistic Flows

Comments: 32 pages, 9 figures, accepted to MNRAS; Changed title and
significantly expanded discussion on the physical interpretation of impulsive
acceleration

Submitted: **2010-04-06**, last modified: **2010-10-04**

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.

[51] oai:arXiv.org:1009.0024 [pdf] - 1034593

Particle Acceleration in Relativistic Magnetized Collisionless
Electron-Ion Shocks

Comments: 24 pages, 20 figures, submitted to ApJ

Submitted: **2010-08-31**

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.

[52] oai:arXiv.org:0910.5740 [pdf] - 902243

Uncertainties of Modeling Gamma-Ray Pulsar Light Curves with Vacuum
Dipole Magnetic Field

Comments: 12 pages, accepted for publication in ApJ

Submitted: **2009-10-30**, last modified: **2010-04-10**

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.

[53] oai:arXiv.org:0910.5741 [pdf] - 902244

Modeling of Gamma-Ray Pulsar Light Curves with Force-Free Magnetic Field

Comments: 22 pages, accepted for publication in ApJ

Submitted: **2009-10-30**, last modified: **2010-04-10**

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

[54] oai:arXiv.org:0912.4990 [pdf] - 902863

Magnetic Amplification by Magnetized Cosmic Rays in SNR Shocks

Comments: 14 pages, 10 figures; submitted to ApJ

Submitted: **2009-12-26**

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

[55] oai:arXiv.org:0905.1970 [pdf] - 1002214

Signatures of a Maxwellian Component in Shock-Accelerated Electrons in
GRBs

Comments: 7 pages, 9 figures, with added fig. to compare with observations,
MNRAS, 400, 330 (2009)

Submitted: **2009-05-12**, last modified: **2009-12-24**

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.

[56] oai:arXiv.org:0904.4918 [pdf] - 315874

Synchrotron Emission from Elliptical Galaxies Consequent to AGN
outbursts

Comments: 15 pages, 11 figures, resubmitted to ApJ, add one more figure, typos
corrected

Submitted: **2009-04-30**, last modified: **2009-12-17**

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

[57] oai:arXiv.org:0911.0120 [pdf] - 30146

Modeling Pulsar Gamma-Ray Light Curves Using Realistic Magnetospheric
Geometries

Comments: 3 pages, in Proceedings of the 10th Asian-Pacific Regional IAU
Meeting, 2008, Aug.3-6, Kunming

Submitted: **2009-10-31**

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.

[58] oai:arXiv.org:0908.0946 [pdf] - 1003325

Magnetohydrodynamic Shallow Water Waves: Linear Analysis

Comments: 15 pages (emulateapj), 4 figures. Accepted by ApJ (no changes in text
from previous versions, but ApJ wants no acronyms in title). To appear in Oct
2009, Vol 703/704 of the Astrophysical Journal

Submitted: **2009-08-06**, last modified: **2009-08-27**

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.

[59] oai:arXiv.org:0908.3193 [pdf] - 1017171

Synthetic Spectra from PIC Simulations of Relativistic Collisionless
Shocks

Comments: 5 pages, 3 figures, submitted to ApJ Letters

Submitted: **2009-08-22**

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.

[60] oai:arXiv.org:0906.1087 [pdf] - 24988

Whence particle acceleration

Comments: Proceedings of the meeting ``070228 - The Next Decade of GRB
Afterglows'' held March 19-23, 2007 in Amsterdam

Submitted: **2009-06-05**

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.

[61] oai:arXiv.org:0802.3217 [pdf] - 10313

Magnetic field evolution in relativistic unmagnetized collisionless
shocks

Comments: ApJ Letters, published (minor corrections, results unchanged)

Submitted: **2008-02-22**, last modified: **2009-03-09**

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.

[62] oai:arXiv.org:0902.3821 [pdf] - 21648

Extreme Astrophysics with Neutron Stars

Comments: 7 pages, whitepaper submitted to Astro2010 Decadal Survey

Submitted: **2009-02-22**

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

[63] oai:arXiv.org:0902.3589 [pdf] - 21600

Life Cycles of Magnetic Fields in Stellar Evolution

Comments: 8 pages; a white paper submitted to the SSE Panel of the Astro2010
Decadal Survey on behalf of the APS Topical Group on Plasma Astrophysics

Submitted: **2009-02-20**

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

[64] oai:arXiv.org:0901.2578 [pdf] - 315507

Particle Acceleration in Relativistic Magnetized Collisionless Pair
Shocks: Dependence of Shock Acceleration on Magnetic Obliquity

Comments: 25 pages, 22 figures, submitted to ApJ, high-resolution version
available at http://www.astro.princeton.edu/~lsironi/ShockAccel.pdf

Submitted: **2009-01-19**

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)

[65] oai:arXiv.org:0810.4565 [pdf] - 315234

Non-linear Study of Bell's Cosmic Ray Current-driven Instability

Comments: 17 pages, 19 figures, submitted to ApJ

Submitted: **2008-10-26**

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.

[66] oai:arXiv.org:0810.4014 [pdf] - 17722

Radiative cooling in relativistic collisionless shocks. Can simulations
and experiments probe relevant GRB physics?

Comments: 15 pages, submitted to ApJ

Submitted: **2008-10-22**

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.

[67] oai:arXiv.org:0802.3216 [pdf] - 10312

Particle acceleration in relativistic collisionless shocks: Fermi
process at last?

Comments: 4 pages, 3 figures

Submitted: **2008-02-22**, last modified: **2008-02-24**

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.

[68] oai:arXiv.org:0801.4583 [pdf] - 9545

Long Term Evolution of Magnetic Turbulence in Relativistic Collisionless
Shocks

Comments: 4 pages, 3 figures, contributed talk at the workshop: High Energy
Phenomena in Relativistic Outflows (HEPRO), Dublin, 24-28 September 2007;
Downsampled version for arXiv. Full resolution version available at
http://astro.berkeley.edu/~pchang/proceedings.pdf

Submitted: **2008-01-29**

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.

[69] oai:arXiv.org:0704.3832 [pdf] - 848

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

Comments: 10 pages, 6 figures, accepted to ApJ; Downsampled version for arXiv.
Full resolution figures available at
http://astro.berkeley.edu/~pchang/full_res_weibel.pdf

Submitted: **2007-04-29**, last modified: **2007-10-19**

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)

[70] oai:arXiv.org:0706.3126 [pdf] - 2432

On the structure of relativistic collisionless shocks in electron-ion
plasmas

Comments: 4 pages, 3 figures; submitted to ApJ Letters

Submitted: **2007-06-21**

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.

[71] oai:arXiv.org:astro-ph/0603212 [pdf] - 80421

Pulsar electrodynamics: a time-dependent view

Comments: 6 pages, invited talk at "Astrophysical Sources of High Energy
Particles and Radiation," Torun, June 20 - 24, 2005

Submitted: **2006-03-08**

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.

[72] oai:arXiv.org:astro-ph/0603211 [pdf] - 80420

Simulations of relativistic collisionless shocks: shock structure and
particle acceleration

Comments: 6 pages, invited talk at "Astrophysical Sources of High Energy
Particles and Radiation," Torun, June 20 - 24, 2005

Submitted: **2006-03-08**

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.

[73] oai:arXiv.org:astro-ph/0603147 [pdf] - 316318

Time-dependent force-free pulsar magnetospheres: axisymmetric and
oblique rotators

Comments: Submitted to ApJ Letters; version with higher resolution figures
available at http://www.stanford.edu/~anatolys/FfreePulsars.pdf (900kb)

Submitted: **2006-03-07**

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.

[74] oai:arXiv.org:astro-ph/0512002 [pdf] - 78202

Revised Pulsar Spindown

Comments: 8 pages, 7 figures; accepted to ApJ

Submitted: **2005-11-30**

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.

[75] oai:arXiv.org:astro-ph/0509914 [pdf] - 76495

Near-Infrared, Kilosecond Variability of the Wisps and Jet in the Crab
Pulsar Wind Nebula

Comments: 10 pages, 9 figures, to be published in The Astrophysical Journal

Submitted: **2005-09-30**

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

[76] oai:arXiv.org:astro-ph/0508192 [pdf] - 75081

Particle-In-Cell Simulations of a Nonlinear Transverse Electromagnetic
Wave in a Pulsar Wind Termination Shock

Comments: 10 pages, 4 figures. Accepted for publication in ApJ (2005)

Submitted: **2005-08-08**

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.

[77] oai:arXiv.org:astro-ph/0507553 [pdf] - 74722

Steady-State Electrostatic Layers from Weibel Instability in
Relativistic Collisionless Shocks

Comments: 9 pages, 3 figures, comments welcome

Submitted: **2005-07-25**

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.

[78] oai:arXiv.org:astro-ph/0401614 [pdf] - 62484

Green Bank Telescope Observations of the Eclipse of Pulsar "A" in the
Double Pulsar Binary PSR J0737-3039

Comments: 12 pages, 3 figures

Submitted: **2004-01-29**, last modified: **2004-10-14**

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.

[79] oai:arXiv.org:astro-ph/0406321 [pdf] - 65478

Orbital Modulation of the Apparent Dispersion Measure Towards PSR
J0737-3039A

Comments: This paper has been withdrawn

Submitted: **2004-06-14**, last modified: **2004-06-17**

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

[80] oai:arXiv.org:astro-ph/0404159 [pdf] - 64067

Probing Relativistic Winds: The case of PSRJ07370-3039 A & B

Comments: To Appear in Proceedings of the 2004 Apsen Winter Conf. on
Astrophysics `Binary Radio Pulsars', ASP Conf Series, F. Rasio & I. Stairs,
eds

Submitted: **2004-04-07**

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

[81] oai:arXiv.org:astro-ph/0402123 [pdf] - 62637

Time-dependence in Relativistic Collisionless Shocks: Theory of the
Variable "Wisps" in the Crab Nebula

Comments: 13 pages, 4 figures, accepted to ApJ. High-resolution figures and
mpeg movies available at http://astron.berkeley.edu/~anatoly/wisp/

Submitted: **2004-02-05**

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.

[82] oai:arXiv.org:astro-ph/0402025 [pdf] - 62539

Orientations of Spin and Magnetic Dipole Axes of Pulsars in the
J0737--3039 Binary Based on Polarimetry Observations at the Green Bank
Telescope

Comments: Submitted to ApJ Letters

Submitted: **2004-02-01**

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}$.

[83] oai:arXiv.org:astro-ph/0310731 [pdf] - 60394

Electrodynamics of pulsar magnetospheres

Comments: Invited review, to appear in "Young Neutron Stars and Their
Environments" (IAU Symposium 218, ASP Conference Proceedings), eds F. Camilo
and B. M. Gaensler. 8 pages. Figures best viewed if printed

Submitted: **2003-10-24**

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.

[84] oai:arXiv.org:astro-ph/0201360 [pdf] - 47304

Simulations of pulsar wind formation

Comments: To appear in "Neutron Stars in Supernova Remnants" (ASP Conference
Proceedings), eds P. O. Slane and B. M. Gaensler. 6 pages. Figures degraded
to fit on astro-ph

Submitted: **2002-01-22**

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.

[85] oai:arXiv.org:astro-ph/0108074 [pdf] - 44050

Propagation of thermonuclear flames on rapidly rotating neutron stars:
extreme weather during type I X-ray bursts

Comments: Submitted to ApJ, 23 pages (including 8 figures), abstract abridged

Submitted: **2001-08-04**

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.