Search

The nonlinear interaction between electromagnetic waves and plasmas attracts significant attention in astrophysics because it can affect the propagation of Fast Radio Bursts (FRBs) -- luminous millisecond-duration pulses detected at radio frequency. The filamentation instability (FI) -- a type of nonlinear wave-plasma interaction -- is considered to be dominant near FRB sources, and its nonlinear development may also affect the inferred dispersion measure of FRBs. In this paper, we carry out fully kinetic particle-in-cell simulations of the FI in unmagnetized pair plasmas. Our simulations show that the FI generates transverse density filaments, and that the electromagnetic wave propagates in near vacuum between them, as in a waveguide. The density filaments keep merging until force balance between the wave ponderomotive force and the plasma pressure gradient is established. We estimate the merging timescale and discuss the implications of filament merging for FRB observations., Submitted to MNRAS

In galaxy clusters, the intracluster medium (ICM) is expected to host a diffuse, long-lived, and invisible population of "fossil" cosmic-ray electrons (CRe) with 1-100 MeV energies. These CRe, if re-accelerated by 100x in energy, can contribute synchrotron luminosity to cluster radio halos, relics, and phoenices. Re-acceleration may be aided by CRe scattering upon the ion-Larmor-scale waves that spawn when ICM is compressed, dilated, or sheared. We study CRe scattering and energy gain due to ion cyclotron (IC) waves generated by continuously-driven compression in 1D fully kinetic particle-in-cell simulations. We find that pitch-angle scattering of CRe by IC waves induces energy gain via magnetic pumping. In an optimal range of IC-resonant momenta, CRe may gain up to ~10-30% of their initial energy in one compress/dilate cycle with magnetic field amplification ~3-6x, assuming adiabatic decompression without further scattering and averaging over initial pitch angle., 23 pages, 18 figures; accepted in ApJ with new appendices and edits

Turbulence driven by AGN activity, cluster mergers and galaxy motion constitutes an attractive energy source for heating the intracluster medium (ICM). How this energy dissipates into the ICM plasma remains unclear, given its low collisionality and high magnetization (precluding viscous heating by Coulomb processes). Kunz et al. 2011 proposed a viable heating mechanism based on the anisotropy of the plasma pressure (gyroviscous heating) under ICM conditions. The present paper builds upon that work and shows that particles can be gyroviscously heated by large-scale turbulent fluctuations via magnetic pumping. We study how the anisotropy evolves under a range of forcing frequencies, what waves and instabilities are generated and demonstrate that the particle distribution function acquires a high energy tail. For this, we perform particle-in-cell simulations where we periodically vary the mean magnetic field $\textbf{B}(t)$. When $\textbf{B}(t)$ grows (dwindles), a pressure anisotropy $P_{\perp}>P_{\parallel}$ ($P_{\perp}< P_{\parallel}$) builds up ($P_{\perp}$ and $P_{\parallel}$ are, respectively, the pressures perpendicular and parallel to $\textbf{B}(t)$). These pressure anisotropies excite mirror ($P_{\perp}>P_{\parallel}$) and oblique firehose ($P_{\parallel}>P_{\perp}$) instabilities, which trap and scatter the particles, limiting the anisotropy and providing a channel to heat the plasma. The efficiency of this mechanism depends on the frequency of the large-scale turbulent fluctuations and the efficiency of the scattering the instabilities provide in their nonlinear stage. We provide a simplified analytical heating model that captures the phenomenology involved. Our results show that this process can be relevant in dissipating and distributing turbulent energy at kinetic scales in the ICM., 24 pages, 17 figures, submitted to ApJ

We investigate with relativistic MHD simulations the dissipation physics of BL Lac jets, by studying the synchrotron polarization signatures of particles accelerated by the kink instability in a magnetically-dominated plasma column. The nonlinear stage of the kink instability generates current sheets, where particles can be efficiently accelerated via magnetic reconnection. We identify current sheets as regions where s = J d/B is above some predefined threshold (where B is the field strength, J the current density and d the grid scale), and assume that the particle injection efficiency scales as proportional to the square of the current. X-ray emitting particles have short cooling times, so they only probe the field geometry of their injection sites. In contrast, particles emitting in the optical band, which we follow self-consistently as they propagate away from their injection sites while cooling, sample a larger volume, and so they may be expected to produce different polarimetric signatures. We find that the degree of polarization is roughly the same between X-ray and optical bands, because even the optical-emitting particles do not travel far from the current sheet where they were injected, due to lack of sufficient kink-generated turbulence. The polarization angle shows a different temporal evolution between the two bands, due to the different regions probed by X-ray and optical emitting particles. In view of the upcoming IXPE satellite, our results can help constrain whether kink-induced reconnection (as opposed to shocks) can be the source of multi-wavelength emission from BL Lacs., 12 pages, 14 figures, MNRAS in press

We perform two-dimensional particle-in-cell simulations of magnetic reconnection in electron-ion plasmas subject to strong Compton cooling and calculate the X-ray spectra produced by this process. The simulations are performed for trans-relativistic reconnection with magnetization $1\leq \sigma\leq 3$ (defined as the ratio of magnetic tension to plasma rest-mass density), which is expected in the coronae of accretion disks around black holes. We find that magnetic dissipation proceeds with inefficient energy exchange between the heated ions and the Compton-cooled electrons. As a result, most electrons are kept at a low temperature in Compton equilibrium with radiation, and so thermal Comptonization cannot reach photon energies $\sim 100$ keV observed from accreting black holes. Nevertheless, magnetic reconnection efficiently generates $\sim 100$ keV photons because of mildly relativistic bulk motions of the plasmoid chain formed in the reconnection layer. Comptonization by the plasmoid motions dominates the radiative output and controls the peak of the radiation spectrum $E_{\rm pk}$. We find $E_{\rm pk}\sim 40$ keV for $\sigma=1$ and $E_{\rm pk}\sim100$ keV for $\sigma=3$. In addition to the X-ray peak around 100 keV, the simulations show a non-thermal MeV tail emitted by a non-thermal electron population generated near X-points of the reconnection layer. The results are consistent with the typical hard state of accreting black holes. In particular, we find that the spectrum of Cygnus~X-1 is well explained by electron-ion reconnection with $\sigma\sim 3$., Comment: Accepted for publication in MNRAS: 16 pages, 15 figures, 1 table, 8 appendices. Part-I can be accessed at arXiv:2107.00263

Spider pulsars are compact binary systems composed of a millisecond pulsar and a low-mass companion. The relativistic magnetically dominated pulsar wind impacts onto the companion, ablating it and slowly consuming its atmosphere. The interaction forms an intrabinary shock, a proposed site of particle acceleration. We perform global fully kinetic particle-in-cell simulations of the intrabinary shock, assuming that the pulsar wind consists of plane-parallel stripes of alternating polarity and that the shock wraps around the companion. We find that particles are efficiently accelerated via shock-driven reconnection. We extract first-principles synchrotron spectra and light curves, which are in good agreement with X-ray observations: (1) the synchrotron spectrum is nearly flat, F ν ∝ const; (2) when the pulsar spin axis is nearly aligned with the orbital angular momentum, the light curve displays two peaks, just before and after the pulsar eclipse (pulsar superior conjunction), separated in phase by ∼0.8 rad; (3) the peak flux exceeds the one at the inferior conjunction by a factor of 10.

International audience; Collisionless, magnetized turbulence offers a promising framework for the generation of nonthermal high-energy particles in various astrophysical sites. Yet, the detailed mechanism that governs particle acceleration has remained subject to debate. By means of 2D and 3D particle-in-cell, as well as 3D (incompressible) magnetohydrodynamic (MHD) simulations, we test here a recent model of nonresonant particle acceleration in strongly magnetized turbulence [Lemoine, Phys. Rev. D 104, 063020 (2021)], which ascribes the energization of particles to their continuous interaction with the random velocity flow of the turbulence, in the spirit of the original Fermi model. To do so, we compare, for a large number of particles that were tracked in the simulations, the predicted and the observed histories of particles momenta. The predicted history is that derived from the model, after extracting from the simulations, at each point along the particle trajectory, the three force terms that control acceleration: the acceleration of the field line velocity projected along the field line direction, its shear projected along the same direction, and its transverse compressive part. Overall, we find a clear correlation between the model predictions and the numerical experiments, indicating that this nonresonant model can successfully account for the bulk of particle energization through Fermi-type processes in strongly magnetized turbulence. We also observe that the parallel shear contribution tends to dominate the physics of energization in the particle-in-cell simulations, while in the magnetohydrodynamic incompressible simulation, both the parallel shear and the transverse compressive term provide about equal contributions.

Magnetic reconnection in relativistic plasmas is well established as a fast and efficient particle accelerator, capable of explaining the most dramatic astrophysical flares. With particle-in-cell simulations, we demonstrate the importance of non-ideal fields for the early stages ("injection") of particle acceleration. Most of the particles ending up with high energies (near or above the mean magnetic energy per particle) must have passed through non-ideal regions where the assumptions of ideal magnetohydrodynamics are broken (i.e., regions with $E>B$ or nonzero $E_\parallel={\bf E}\cdot {\bf B}/B$), whereas particles that do not experience non-ideal fields end up with Lorentz factors of order unity. Thus, injection by non-ideal fields is a necessary prerequisite for further acceleration. Our results have important implications for the origin of nonthermal particles in high-energy astrophysical sources.

We perform two-dimensional particle-in-cell simulations of reconnection in magnetically dominated electron-positron plasmas subject to strong Compton cooling. We vary the magnetization $\sigma\gg1$, defined as the ratio of magnetic tension to plasma inertia, and the strength of cooling losses. Magnetic reconnection under such conditions can operate in magnetically dominated coronae around accreting black holes, which produce hard X-rays through Comptonization of seed soft photons. We find that the particle energy spectrum is dominated by a peak at mildly relativistic energies, which results from bulk motions of cooled plasmoids. The peak has a quasi-Maxwellian shape with an effective temperature of $\sim 100$ keV, which depends only weakly on the flow magnetization and the strength of radiative cooling. The mean bulk energy of the reconnected plasma is roughly independent of $\sigma$, whereas the variance is larger for higher magnetizations. The spectra also display a high-energy tail, which receives $\sim 25$% of the dissipated reconnection power for $\sigma=10$ and $\sim 40$% for $\sigma=40$. We complement our particle-in-cell studies with a Monte-Carlo simulation of the transfer of seed soft photons through the reconnection layer, and find the escaping X-ray spectrum. The simulation demonstrates that Comptonization is dominated by the bulk motions in the chain of Compton-cooled plasmoids and, for $\sigma\sim 10$, yields a spectrum consistent with the typical hard state of accreting black holes., Comment: Accepted for publication in MNRAS. 17 pages, 11 figures, 4 appendices, 1 table

Relativistic magnetized shocks are a natural source of coherent emission, offering a plausible radiative mechanism for Fast Radio Bursts (FRBs). We present first-principles 3D simulations that provide essential information for the FRB models based on shocks: the emission efficiency, spectrum, and polarization. The simulated shock propagates in an $e^\pm$ plasma with magnetization $\sigma>1$. The measured fraction of shock energy converted to coherent radiation is $\simeq 10^{-3} \, \sigma^{-1}$, and the energy-carrying wavenumber of the wave spectrum is $\simeq 4 \,\omega_{\rm c}/c$, where $\omega_{\rm c}$ is the upstream gyrofrequency. The ratio of the O-mode and X-mode energy fluxes emitted by the shock is $\simeq 0.4\,\sigma^{-1}$. The dominance of the X-mode at $\sigma\gg 1$ is particularly strong, approaching 100% in the spectral band around $2\,\omega_{\rm c}$. We also provide a detailed description of the emission mechanism for both X- and O-modes., Comment: 6 pages, 5 figures, accepted to PRL

Magnetic energy around compact objects often dominates over plasma rest mass, and its dissipation can power the object luminosity. We describe a dissipation mechanism which works faster than magnetic reconnection. The mechanism involves two strong Alfv\'en waves with anti-aligned magnetic fields $\boldsymbol{B}_1$ and $\boldsymbol{B}_2$ that propagate in opposite directions along background magnetic field $\boldsymbol{B}_0$ and collide. The collision forms a thin current sheet perpendicular to $\boldsymbol{B}_0$, which absorbs the incoming waves. The current sheet is sustained by electric field $\boldsymbol{E}$ breaking the magnetohydrodynamic condition $E1/2$. It dissipates a large fraction of the wave energy, $f=(2A-1)/A^2$, reaching $100\%$ when $A=1$. The plane geometry allows one to see the dissipation process in a one-dimensional simulation. We also perform two-dimensional simulations, enabling spontaneous breaking of the plane symmetry by the tearing instability of the current sheet. At moderate $A$ of main interest the tearing instability is suppressed. Dissipation transitions to normal, slower, magnetic reconnection at $A\gg 1$. The fast dissipation described in this paper may occur in various objects with perturbed magnetic fields, including magnetars, jets from accreting black holes, and pulsar wind nebulae., Comment: 12 pages, 10 figures, submitted to ApJ

Using 2D particle-in-cell (PIC) plasma simulations we study electron acceleration by temperature anisotropy instabilities, assuming conditions typical of above-the-loop-top (ALT) sources in solar flares. We focus on the long-term effect of $T_{e,\perp} > T_{e,\parallel}$ instabilities by driving the anisotropy growth during the entire simulation time, through imposing a shearing or a compressing plasma velocity ($T_{e,\perp}$ and $T_{e,\parallel}$ are the temperatures perpendicular and parallel to the magnetic field). This magnetic growth makes $T_{e,\perp}/T_{e,\parallel}$ grow due to electron magnetic moment conservation, and amplifies the ratio $\omega_{ce}/\omega_{pe}$ from $\sim 0.53$ to $\sim 2$ ($\omega_{ce}$ and $\omega_{pe}$ are the electron cyclotron and plasma frequencies, respectively). In the regime $\omega_{ce}/\omega_{pe}\lesssim 1.2-1.7$ the instability is dominated by oblique, quasi-electrostatic (OQES) modes, and the acceleration is inefficient. When $\omega_{ce}/\omega_{pe}$ has grown to $\omega_{ce}/\omega_{pe}\gtrsim 1.2-1.7$, electrons are efficiently accelerated by the inelastic scattering provided by unstable parallel, electromagnetic z (PEMZ) modes. After $\omega_{ce}/\omega_{pe}$ reaches $\sim 2$, the electron energy spectra show nonthermal tails that differ between the shearing and compressing cases. In the shearing case, the tail resembles a power-law of index $\alpha_s \sim$ 2.9 plus a high-energy bump reaching $\sim 300$ keV. In the compressing runs, $\alpha_s \sim$ 3.7 with a spectral break above $\sim 500$ keV. This difference can be explained by the different temperature evolutions in these two types of simulations, suggesting a critical role played by the type of anisotropy driving, $\omega_{ce}/\omega_{pe}$ and the electron temperature in the efficiency of the acceleration., Comment: 15 pages, 13 figures

Relativistic jets launched by rotating black holes are powerful emitters of non-thermal radiation. Extraction of the rotational energy via electromagnetic stresses produces magnetically-dominated jets, which may become turbulent. Studies of magnetically-dominated plasma turbulence from first principles show that most of the accelerated particles have small pitch angles, i.e. the particle velocity is nearly aligned with the local magnetic field. We examine synchrotron-self-Compton radiation from anisotropic particles in the fast cooling regime. The small pitch angles reduce the synchrotron cooling rate and promote the role of inverse Compton (IC) cooling, which can occur in two different regimes. In the Thomson regime, both synchrotron and IC components have soft spectra, $\nu F_\nu\propto\nu^{1/2}$. In the Klein-Nishina regime, synchrotron radiation has a hard spectrum, typically $\nu F_\nu\propto\nu$, over a broad range of frequencies. Our results have implications for the modelling of BL Lacs and Gamma-Ray Bursts (GRBs). BL Lacs produce soft synchrotron and IC spectra, as expected when Klein-Nishina effects are minor. The observed synchrotron and IC luminosities are typically comparable, which indicates a moderate anisotropy with pitch angles $\theta\gtrsim0.1$. Rare orphan gamma-ray flares may be produced when $\theta\ll0.1$. The hard spectra of GRBs may be consistent with synchrotron radiation when the emitting particles are IC cooling in the Klein-Nishina regime, as expected for pitch angles $\theta\sim0.1$. Blazar and GRB spectra can be explained by turbulent jets with a similar electron plasma magnetisation parameter, $\sigma_{\rm e}\sim10^4$, which for electron-proton plasmas corresponds to an overall magnetisation $\sigma=(m_{\rm e}/m_{\rm p})\sigma_{\rm e}\sim10$., Comment: accepted for publication by MNRAS

We introduce a toy model for the time-frequency structure of fast radio bursts (FRB), in which the observed emission is produced as a narrowly-peaked intrinsic spectral energy distribution sweeps down in frequency across the instrumental bandpass as a power-law in time. Though originally motivated by emission models which invoke a relativistic shock, the model could in principle apply to a wider range of emission scenarios. We quantify the burst's detectability using the frequency bandwidth over which most of its signal-to-noise ratio (SNR) is accumulated. We demonstrate that by varying just a single parameter of the toy model-the power-law index \beta of the frequency drift rate-one can transform a long (and hence preferentially time-resolved) burst with a narrow time-integrated spectrum into a shorter burst with a broad power-law time-integrated spectrum. We suggest that burst-to-burst diversity in the value of \beta could generate the dichotomy between burst duration and frequency-width recently found by CHIME. In shock models, the value of \beta is related to the radial density profile of external medium, which in light of the preferentially longer duration of bursts from repeating sources may point to diversity in the external environments surrounding repeating versus one-off FRB sources., Comment: 11 pages, 7 figures, to be submitted to ApJ, comments appreciated

Nature's most powerful high-energy sources are capable of accelerating particles to high energy and radiate it away on extremely short timescales, even shorter than the light crossing time of the system. It is yet unclear what physical processes can produce such an efficient acceleration, despite the copious radiative losses. By means of radiative particle-in-cell simulations, we show that magnetically dominated turbulence in pair plasmas subject to strong synchrotron cooling generates a nonthermal particle spectrum with a hard power-law range (slope $p \sim 1$) within a few eddy turnover times. Low pitch-angle particles can significantly exceed the nominal radiation-reaction limit, before abruptly cooling down. The particle spectrum becomes even harder ($p < 1$) over time owing to particle cooling with an energy-dependent pitch-angle anisotropy. The resulting synchrotron spectrum is hard ($\nu F_\nu \propto \nu^s$ with $s \sim 1$). Our findings have important implications for understanding the nonthermal emission from high-energy astrophysical sources, most notably the prompt phase of gamma-ray bursts and gamma-ray flares from the Crab nebula.

Fast Radio Bursts (FRBs) are extreme astrophysical phenomena entering the realm of non-linear optics, a field developed in laser physics. A classical non-linear effect is self-modulation. We examine the propagation of FRBs through the circumburst environment using the idealised setup of a monochromatic linearly-polarised GHz wave propagating through a uniform plasma slab of density $N$ at distance $R$ from the source. We find that self-modulation occurs if the slab is located within a critical radius $R_{\rm crit}\sim 10^{17}(N/10^2{\rm\; cm}^{-3})(L/10^{42}{\rm\; erg\; s}^{-1}){\rm\; cm}$, where $L$ is the isotropic equivalent of the FRB luminosity. Self-modulation breaks the burst into pancakes transverse to the radial direction. When $R\lesssim R_{\rm crit}$, the transverse size of the pancakes is smaller than the Fresnel scale. The pancakes are strongly diffracted as the burst exits the slab, and interference between the pancakes produces a frequency modulation of the observed intensity with a sub-GHz bandwidth. When $R\sim R_{\rm crit}$, the transverse size of the pancakes becomes comparable with the Fresnel scale, and the effect of diffraction is weaker. The observed intensity is modulated on a timescale of ten microseconds, which corresponds to the radial width of the pancakes. Our results suggest that self-modulation may cause the temporal and frequency structure observed in FRBs., accepted for publication in MNRAS

We consider the tilting instability of a magnetically confined spheromak using 3D MHD and relativistic PIC calculations with an application to astrophysical plasmas, specifically those occurring in magnetar magnetospheres. The instability is driven by the counter alignment of the spheromak's intrinsic magnetic dipole with the external magnetic field. Initially the spheromak rotates - tilts - trying to lower its magnetic potential energy. As a result a current sheet forms between the internal magnetic field of a spheromak and the confining field. Magnetic reconnection sets in; this leads to the annihilation of the newly counter-aligned magnetic flux of the spheromak. This occurs on few Alfv\'en time scales. In the case of higher order (second order) spheromak, the internal core is first pushed out of the envelope, resulting in formation of two nearly independent tilting spheromaks. Thus, the magnetically twisted outer shell cannot stabilize the inner core. During dissipation, helicity of the initial spheromak is carried away by torsional Alfv\'en waves, violating the assumptions of the Taylor relaxation theorem. In applications to magnetars' giant flares, fast development of tilting instabilities, and no stabilization of the higher order spheromaks, make it unlikely that trapped spheromaks are responsible for the tail emission lasting hundreds of seconds., Comment: 20 pages, 10 figures, accepted for publication in Journal of Plasma Physics