BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Anatoly Spitkovsky

Anatoly Spitkovsky appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.HEphysics.plasm-phastro-ph.GAastro-phastro-ph.COastro-ph.IMphysics.space-ph

Catalog footprint

What is connected

40works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

40 published item(s)

preprint2020arXiv

Electron acceleration in non-relativistic quasi-perpendicular collisionless shocks

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Origin of Pulsar Radio Emission

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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Filling the Magnetospheres of Weak Pulsars

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Modeling high-energy pulsar lightcurves from first principles

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Ab-initio pulsar magnetosphere: the role of general relativity

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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

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én 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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Particle acceleration in axisymmetric pulsar current sheets

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Simultaneous Acceleration of Protons and Electrons at Nonrelativistic Quasiparallel Collisionless Shocks

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Three-dimensional Analytical Description of Magnetised Winds from Oblique Pulsars

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Linear Instabilities Driven by Differential Rotation in Very Weakly Magnetized Plasmas

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Relativistic Reconnection: an Efficient Source of Non-Thermal Particles

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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

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énic 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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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

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énic 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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Maximum Energy of Accelerated Particles in Relativistic Collisionless Shocks

The afterglow emission from gamma-ray bursts (GRBs) is usually interpreted as synchrotron radiation from electrons accelerated at the GRB external shock, that propagates with relativistic velocities into the magnetized interstellar medium. By means of multi-dimensional particle-in-cell simulations, we investigate the acceleration performance of weakly magnetized relativistic shocks, in the magnetization range 0<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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Time-Dependent 3D Magnetohydrodynamic Pulsar Magnetospheres: Oblique Rotators

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On the Spin-Down of Intermittent Pulsars

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Physical Conditions in the Reconnection Layer in Pulsar Magnetospheres

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Electron Injection by Whistler Waves in Non-relativistic Shocks

Electron acceleration to non-thermal, ultra-relativistic energies (~ 10-100 TeV) is revealed by radio and X-ray observations of shocks in young supernova remnants (SNRs). The diffusive shock acceleration (DSA) mechanism is usually invoked to explain this acceleration, but the way in which electrons are initially energized or 'injected' into this acceleration process starting from thermal energies is an unresolved problem. In this paper we study the initial acceleration of electrons in non-relativistic shocks from first principles, using two- and three-dimensional particle-in-cell (PIC) plasma simulations. We systematically explore the space of shock parameters (the Alfvénic 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énic 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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Resistive Solutions for Pulsar Magnetospheres

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Impulsive Acceleration of Strongly Magnetized Relativistic Flows

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Particle Acceleration in Relativistic Magnetized Collisionless Electron-Ion Shocks

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Magnetic Amplification by Magnetized Cosmic Rays in SNR Shocks

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

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Magnetohydrodynamic Shallow Water Waves: Linear Analysis

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

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

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)

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Synthetic Spectra from PIC Simulations of Relativistic Collisionless Shocks

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.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

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

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, λ_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 λ_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.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access