Search

Some of the most energetic pulsars exhibit rotation-modulated gamma-ray emission in the $0.1$ to $100$ GeV band. The luminosity of this emission is typically $0.1\text{-}10\%$ of the pulsar spin-down power (gamma-ray efficiency), implying that a significant fraction of the available electromagnetic energy is dissipated in the magnetosphere and reradiated as high-energy photons. To investigate this phenomenon we model a pulsar magnetosphere using 3D particle-in-cell simulations with strong synchrotron cooling. We particularly focus on the dynamics of the equatorial current sheet where magnetic reconnection and energy dissipation take place. Our simulations demonstrate that a fraction of the spin-down power dissipated in the magnetospheric current sheet is controlled by the rate of magnetic reconnection at microphysical plasma scales and only depends on the pulsar inclination angle. We demonstrate that the maximum energy and the distribution function of accelerated pairs is controlled by the available magnetic energy per particle near the current sheet, the magnetization parameter. The shape and the extent of the plasma distribution is imprinted in the observed synchrotron emission, in particular, in the peak and the cutoff of the observed spectrum. We study how the strength of synchrotron cooling affects the observed variety of spectral shapes. Our conclusions naturally explain why pulsars with higher spin-down power have wider spectral shapes and, as a result, lower gamma-ray efficiency., 22 pages, 12 figures

One scenario for the generation of fast radio bursts (FRBs) is magnetic reconnection in a current sheet of the magnetar wind. Compressed by a strong magnetic pulse induced by a magnetar flare, the current sheet fragments into a self-similar chain of magnetic islands. Time-dependent plasma currents at their interfaces produce coherent radiation during their hierarchical coalescence. We investigate this scenario using 2D radiative relativistic particle-in-cell simulations to compute the efficiency of the coherent emission and to obtain frequency scalings. Consistent with expectations, a fraction of the reconnected magnetic field energy, $f\sim 0.002$, is converted to packets of high-frequency fast magnetosonic waves which can escape from the magnetar wind as radio emission. In agreement with analytical estimates, we find that magnetic pulses of $10^{47}\text{erg}\;\text{s}^{-1}$ can trigger relatively narrowband GHz emission with luminosities of approximately $10^{42}\text{erg}\;\text{s}^{-1}$, sufficient to explain bright extragalactic FRBs. The mechanism provides a natural explanation for a downward frequency drift of burst signals, as well as the $\sim 100\;\text{ns}$ substructure recently detected in FRB 20200120E., 19 pages, 9 figures, Accepted for publication by ApJL

Astrophysical collisionless shocks are among the most powerful particle accelerators in the Universe. Generated by violent interactions of supersonic plasma flows with the interstellar medium, supernova remnant shocks are observed to amplify magnetic fields1 and accelerate electrons and protons to highly relativistic speeds2–4. In the well-established model of diffusive shock acceleration5, relativistic particles are accelerated by repeated shock crossings. However, this requires a separate mechanism that pre-accelerates particles to enable shock crossing. This is known as the ‘injection problem’, which is particularly relevant for electrons, and remains one of the most important puzzles in shock acceleration6. In most astrophysical shocks, the details of the shock structure cannot be directly resolved, making it challenging to identify the injection mechanism. Here we report results from laser-driven plasma flow experiments, and related simulations, that probe the formation of turbulent collisionless shocks in conditions relevant to young supernova remnants. We show that electrons can be effectively accelerated in a first-order Fermi process by small-scale turbulence produced within the shock transition to relativistic non-thermal energies, helping overcome the injection problem. Our observations provide new insight into electron injection at shocks and open the way for controlled laboratory studies of the physics underlying cosmic accelerators. In laser–plasma experiments complemented by simulations, electron acceleration is observed in turbulent collisionless shocks. This work clarifies the pre-acceleration to relativistic energies required for the onset of diffusive shock acceleration.

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., Comment: Submitted to MNRAS. 14 pages, 18 figures

Using the field-particle correlation technique, we examine the particle energization in a 1D-2V continuum Vlasov--Maxwell simulation of a perpendicular magnetized collisionless shock. The combination of the field-particle correlation technique with the high fidelity representation of the particle distribution function provided by a direct discretization of the Vlasov equation allows us to ascertain the details of the exchange of energy between the electromagnetic fields and the particles in phase space. We identify the velocity-space signatures of shock-drift acceleration of the ions and adiabatic heating of the electrons due to the perpendicular collisionless shock by constructing a simplified model with the minimum ingredients necessary to produce the observed energization signatures in the self-consistent Vlasov-Maxwell simulation. We are thus able to completely characterize the energy transfer in the perpendicular collisionless shock considered here and provide predictions for the application of the field-particle correlation technique to spacecraft measurements of collisionless shocks., 49 pages, 20 figures

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., Comment: 9 pages, 4 figures

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

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., Comment: 16 pages, 17 figures. Accepted for publication in MNRAS

We use the one-dimensional TRISTAN-MP particle-in-cell code to model the nonlinear evolution of the whistler heat flux instability that was proposed by Gary et al. (1999, 2000) to regulate the electron heat flux in the solar wind and astrophysical plasmas. The simulations are initialized with electron velocity distribution functions typical for the solar wind. We perform a set of simulations at various initial values of the electron heat flux and $\beta_{e}$. The simulations show that parallel whistler waves produced by the whistler heat flux instability saturate at amplitudes consistent with the spacecraft measurements. The simulations also reproduce the correlations of the saturated whistler wave amplitude with the electron heat flux and $\beta_{e}$ revealed in the spacecraft measurements. The major result is that parallel whistler waves produced by the whistler heat flux instability do not significantly suppress the electron heat flux. The presented simulations indicate that coherent parallel whistler waves observed in the solar wind are unlikely to regulate the heat flux of solar wind electrons.

We present the first realistic 3D simulations of flame front instabilities during type I X-ray bursts. The unperturbed front is characterized by the balance between the pressure gradient and the Coriolis force of a spinning neutron star (ν = 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 ∼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 light curves of the bursts.

Collisionless shocks are common features in space and astrophysical systems where supersonic plasma flows interact, such as in the solar wind, the heliopause, and supernova remnants. Recent experimental capabilities and diagnostics allow detailed laboratory investigations of high-Mach-number shocks, which therefore can become a valuable way to understand shock dynamics in various astrophysical environments. Using 2D particle-in-cell simulations with a Coulomb binary collision operator, we demonstrate the mechanism for generation of energetic electrons and experimental requirements for detecting this process in the laboratory high-Mach-number collisionless shocks. We show through a parameter study that electron acceleration by magnetized collisionless shocks is feasible in laboratory experiments with laser-driven expanding plasmas., 7 pages, 5 figures

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., 14 pages, 11 figures, Accepted for publication in MNRAS

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., Comment: 9 pages, 6 figures, submitted to ApJ

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., Comment: accepted to JPP

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., Comment: 22 pages, 15 figures, submitted to ApJ

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., 16 pages, 17 figures, Accepted for publication in MNRAS

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., Comment: 15 pages, 6 Figures. MNRAS in press

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., Comment: 11 pages, 8 figures, submitted to ApJ

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., Comment: 6 pages, 7 figures

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., Comment: Contribution to the 6th International Symposium on High Energy Gamma-Ray Astronomy (Gamma2016), Heidelberg, Germany. To be published in the AIP Conference Proceedings

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., Comment: 30 pages, 19 figures. Accepted to ApJ

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., Comment: 25 pages, 23 figures, 4 tables. Accepted for publication in ApJ

We present the results of particle-in-cell simulations of the whistler anisotropy instability that results in magnetospheric chorus wave excitation. The simulations were carried out using, for the first time for this problem, the 2D TRISTAN-massively parallelized code, widely used before in the modeling of astrophysical shocks. The code has been modified to allow for two populations of electrons: cold electrons (which maintain the wave propagation) and hot electrons (which provide the wave growth). For the hot electrons, the anisotropic form of the relativistic Maxwell–Juttner distribution is implemented. We adopt the standard approximation of a parabolic magnetic field to simulate the Earth's magnetic field close to the equator. Simulations with different background magnetic field inhomogeneity strengths demonstrate that higher inhomogeneity yields lower frequency chirping rates and, eventually, it suppresses chorus generation. The results are in agreement with other numerical simulations and the theoretical predictions for the frequency chirping rates.

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., 7 pages, 4 figures

A study of the transition from collisional to collisionless plasma flows has been carried out at the National Ignition Facility using high Mach number (M>4) counterstreaming plasmas. In these experiments, CD-CD and CD-CH planar foils separated by 6-10 mm are irradiated with laser energies of 250 kJ per foil, generating ∼1000 km/s plasma flows. Varying the foil separation distance scales the ion density and average bulk velocity and, therefore, the ion-ion Coulomb mean free path, at the interaction region at the midplane. The characteristics of the flow interaction have been inferred from the neutrons and protons generated by deuteron-deuteron interactions and by x-ray emission from the hot, interpenetrating, and interacting plasmas. A localized burst of neutrons and bright x-ray emission near the midpoint of the counterstreaming flows was observed, suggesting strong heating and the initial stages of shock formation. As the separation of the CD-CH foils increases we observe enhanced neutron production compared to particle-in-cell simulations that include Coulomb collisions, but do not include collective collisionless plasma instabilities. The observed plasma heating and enhanced neutron production is consistent with the initial stages of collisionless shock formation, mediated by the Weibel filamentation instability.

Collisions of high Mach number flows occur frequently in astrophysics, and the resulting shock waves are responsible for the properties of many astrophysical phenomena, such as supernova remnants, Gamma Ray Bursts and jets from Active Galactic Nuclei. Because of the low density of astrophysical plasmas, the mean free path due to Coulomb collisions is typically very large. Therefore, most shock waves in astrophysics are “collisionless”, since they form due to plasma instabilities and self-generated magnetic fields. Laboratory experiments at the laser facilities can achieve the conditions necessary for the formation of collisionless shocks, and will provide a unique avenue for studying the nonlinear physics of collisionless shock waves. We are performing a series of experiments at the Omega and Omega-EP lasers, in Rochester, NY, with the goal of generating collisionless shock conditions by the collision of two high-speed plasma flows resulting from laser ablation of solid targets using ∼10 16 W/cm 2 laser irradiation. The experiments will aim to answer several questions of relevance to collisionless shock physics: the importance of the electromagnetic filamentation (Weibel) instabilities in shock formation, the self-generation of magnetic fields in shocks, the influence of external magnetic fields on shock formation, and the signatures of particle acceleration in shocks. Our first experiments using Thomson scattering diagnostics studied the plasma state from a single foil and from double foils whose flows collide “head-on”. Our data showed that the flow velocity and electron density were 10 8 cm/s and 10 19 cm −3 , respectively, where the Coulomb mean free path is much larger than the size of the interaction region. Simulations of our experimental conditions show that weak Weibel mediated current filamentation and magnetic field generation were likely starting to occur. This paper presents the results from these first Omega experiments.

We are performing scaled astrophysics experiments on Omega and on NIF. Laser driven counter-streaming interpenetrating supersonic plasma flows can be studied to understand astrophysical electromagnetic plasma phenomena in a controlled laboratory setting. In our Omega experiments, the counter-streaming flow plasma state is measured using Thomson scattering diagnostics, demonstrating the plasma flows are indeed super-sonic and in the collisionless regime. We observe a surprising additional electron and ion heating from ion drag force in the double flow experiments that are attributed to the ion drag force and electrostatic instabilities. [1] A proton probe is used to image the electric and magnetic fields. We observe unexpected large, stable and reproducible electromagnetic field structures that arise in the counter-streaming flows [2]. The Biermann battery magnetic field generated near the target plane, advected along the flows, and recompressed near the midplane explains the cause of such self-organizing field structures [3]. A D3He implosion proton probe image showed very clear filamentary structures; three-dimensional Particle-In-Cell simulations and simulated proton radiography images indicate that these filamentary structures are generated by Weibel instabilities and that the magnetization level (ratio of magnetic energy over kinetic energy in the system) is ∼0.01 [4]. These findings have very high astrophysical relevance and significant implications. We expect to observe true collisionless shock formation when we use >100 kJ laser energy on NIF.

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., Comment: 15 pages, 12 figures, published in MNRAS

Current models of pulsar magnetospheres typically assume either a complete absence of plasma or abundant ideal plasma filling the magnetosphere in order to compute the field structure. The latter condition is thought to be closer to reality, but we know of a number of pulsars in which the ideal conditions break down, resulting in dissipation and high-energy emission. In this work we formulate a resistive force-free scheme that allows us to consider the effects of resistive plasma and accelerating fields on the magnetospheric structure. We run numerical simulations to construct a family of resistive solutions that smoothly bridges the gap between the vacuum and the force-free magnetosphere solutions. We further provide a self-consistent model for the spin-down of intermittent pulsars, pulsars which appear to transition between radio-loud and radio-quiet states with different spin-down rates. Finally, we present models for high-energy emission from reconnecting current sheets in Gamma-ray pulsars.

The relativistic wind of pulsars consists of toroidal stripes of opposite magnetic field polarity, separated by current sheets of hot plasma. By means of multi-dimensional particle-in-cell simulations, we investigate particle acceleration and magnetic field dissipation at the termination shock of a relativistic striped pulsar wind. At the shock, the flow compresses and the alternating fields annihilate by driven magnetic reconnection. Irrespective of the stripe wavelength λ or the wind magnetization σ (in the regime σ ≫1 of magnetically dominated flows), shock-driven reconnection transfers all the magnetic energy of alternating fields to the particles. In the limit λ/(rL σ) ≫ 1, where rL is the relativistic Larmor radius in the wind, the post-shock spectrum approaches a flat power-law tail with slope around -1.5, populated by particles accelerated by the reconnection electric field. Our findings place important constraints on the models of non-thermal radiation from Pulsar Wind Nebulae.

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.

In-situ measurements taken by the Voyager 2 spacecraft suggest that the solar wind termination shock is significantly affected by the presence of pickup ions that are produced in the inner heliosphere due to charge exchange between interstellar neutrals and the solar wind ions. We use a fully kinetic particle-in-cell method to self-consistently simulate the shock with all physical properties available from Voyager 2. We have performed a set of simulations with varying velocity distribution functions for the pickup ions, since it was not determined by Voyager's measurements. We show that the measurements suggest that the pickup ions upstream of the shock are more energetic than generally believed. If their velocity distribution function assumes a filled-shell shape in the wind frame, the maximum cutoff speed for the pickup ions should be 650 km s−1 in order to reproduce the measurements, which is almost twice the local wind speed. We suggest that pickup ions upstream of the shock are energized by adiabatic compression of the solar wind plasma as well as due to an enhanced level of turbulence in a broad foreshock region.

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

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.