Search for: All records

Abstract Particle-in-cell (PIC) simulations have shown that relativistic collisionless shocks mediated by the Weibel instability accelerate ∼1% of incoming particles, while the majority are transmitted through the shock and become thermalized. The microphysical processes that determine whether an incoming particle will be transmitted or reflected are poorly understood. We study the microphysics of particle reflection in Weibel-mediated shocks by tracking a shell of test particles in a PIC simulation of a shock in pair plasma. We find that electrons in positron-dominated filaments and positrons in electron-dominated filaments efficiently reflect off of strong magnetic structures at the shock. To participate in diffusive shock acceleration, however, these reflected particles headed toward the upstream must avoid getting advected downstream. This is enabled by incoming filaments, which trap reflected particles carrying the same sign of current as the filaments. The final injection efficiency on the order of ∼1% thus results from the effectiveness of the initial reflection at the shock and the reflected particles’ probability of survival in the upstream postreflection. We develop a model that predicts the fraction of high-energy particles as a function of the properties of Weibel filamentation. 

Abstract Collisionless shocks tend to send charged particles into the upstream, driving electric currents through the plasma. Using kinetic particle-in-cell simulations, we investigate how the background thermal plasma neutralizes such currents in the upstream of quasi-parallel non-relativistic electron–proton shocks. We observe distinct processes in different regions: the far upstream, the shock precursor, and the shock foot. In the far upstream, the current is carried by nonthermal protons, which drive electrostatic modes and produce suprathermal electrons that move toward upstream infinity. Closer to the shock (in the precursor), both the current density and the momentum flux of the beam increase, which leads to electromagnetic streaming instabilities that contribute to the thermalization of suprathermal electrons. At the shock foot, these electrons are exposed to shock-reflected protons, resulting in a two-stream type instability. We analyze these processes and the resulting heating through particle tracking and controlled simulations. In particular, we show that the instability at the shock foot can make the effective thermal speed of electrons comparable to the drift speed of the reflected protons. These findings are important for understanding both the magnetic field amplification and the processes that may lead to the injection of suprathermal electrons into diffusive shock acceleration. 

Abstract We study a relativistic collisionless electron–positron shock propagating into an unmagnetized ambient medium using 2D particle-in-cell simulations of unprecedented duration and size. The shock generates intermittent magnetic structures of increasingly larger size as the simulation progresses. Toward the end of our simulation, at around 26,000 plasma times, the magnetic coherence scale approachesλ∼ 100 plasma skin depths, both ahead and behind the shock front. We anticipate a continued growth ofλbeyond the time span of our simulation, as long as the shock accelerates particles to increasingly higher energies. The post-shock field is concentrated in localized patches, which maintain a local magnetic energy fractionε_B∼ 0.1. Particles randomly sampling the downstream fields spend most of their time in low field regions (ε_B≪ 0.1) but emit a large fraction of the synchrotron power in the localized patches with strong fields (ε_B∼ 0.1). Our results have important implications for models of gamma-ray burst afterglows. 

Abstract Magnetized plasma columns and extended magnetic structures with both footpoints anchored to a surface layer are an important building block of astrophysical dissipation models. Current loops shining in X-rays during the growth of plasma instabilities are observed in the corona of the Sun and are expected to exist in highly magnetized neutron star magnetospheres and accretion disk coronae. For varying twist and system sizes, we investigate the stability of line-tied force-free flux tubes and the dissipation of twist energy during instabilities using linear analysis and time-dependent force-free electrodynamics simulations. Kink modes (m= 1) and efficient magnetic energy dissipation develop for plasma safety factorsq≲ 1, whereqis the inverse of the number of magnetic field line windings per column length. Higher-order fluting modes (m> 1) can distort equilibrium flux tubes forq> 1 but induce significantly less dissipation. In our analysis, the characteristic pitch ${\tilde{\mu }}_0$of flux-tube field lines determines the growth rate ( $\propto {\tilde{\mu }}_0^3$) and minimum wavelength of the kink instability ( $\propto {\tilde{\mu }}_0^{-1}$). We use these scalings to determine a minimum flux tube length for the growth of the kink instability for any given ${\tilde{\mu }}_0$. By drawing analogies to idealized magnetar magnetospheres with varying regimes of boundary shearing rates, we discuss the expected impact of the pitch-dependent growth rates for magnetospheric dissipation in magnetar conditions. 

ABSTRACT The non-linear 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 non-linear wave–plasma interaction – is considered to be dominant near FRB sources, and its non-linear 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 time-scale and discuss the implications of filament merging for FRB observations. 

ABSTRACT In general relativistic magnetohydrodynamic (GRMHD) simulations, accreted magnetic flux on the black hole horizon episodically decays, during which magnetic reconnection heats up the plasma near the horizon, potentially powering high-energy flares like those observed in M87* and Sgr A*. We study the mm observational counterparts of such flaring episodes in very high resolution GRMHD simulations. The change in 230 GHz flux during the expected high energy flares depends primarily on the efficiency of accelerating γ ≳ 100 (Te ≳ 1011 K) electrons. For models in which the electrons are heated to Te ∼ 1011 K during flares, the hot plasma produced by reconnection significantly enhances 230 GHz emission and increases the size of the 230 GHz image. By contrast, for models in which the electrons are heated to higher temperatures (which we argue are better motivated), the reconnection-heated plasma is too hot to produce significant 230 GHz synchrotron emission, and the 230 GHz flux decreases during high energy flares. We do not find a significant change in the mm polarization during flares as long as the emission is Faraday thin. We also present expectations for the ring-shaped image as observed by the Event Horizon Telescope during flares, as well as multiwavelength synchrotron spectra. Our results highlight several limitations of standard post-processing prescriptions for the electron temperature in GRMHD simulations. We also discuss the implications of our results for current and future observations of flares in Sgr A*, M87*, and related systems. Appendices contain detailed convergence studies with respect to resolution and plasma magnetization. 

Abstract The origins of the various outbursts of hard X-rays from magnetars (highly magnetized neutron stars) are still unknown. We identify instabilities in relativistic magnetospheres that can explain a range of X-ray flare luminosities. Crustal surface motions can twist the magnetar magnetosphere by shifting the frozen-in footpoints of magnetic field lines in current-carrying flux bundles. Axisymmetric (2D) magnetospheres exhibit strong eruptive dynamics, i.e., catastrophic lateral instabilities triggered by a critical footpoint displacement ofψ_crit≳π. In contrast, our new three-dimensional (3D) twist models with finite surface extension capture important non-axisymmetric dynamics of twisted force-free flux bundles in dipolar magnetospheres. Besides the well-established global eruption resulting (as in 2D) from lateral instabilities, such 3D structures can develop helical, kink-like dynamics, and dissipate energy locally (confined eruptions). Up to 25% of the induced twist energy is dissipated and available to power X-ray flares in powerful global eruptions, with most of our models showing an energy release in the range of the most common X-ray outbursts, ≲10⁴³erg. Such events occur when significant energy builds up while deeply buried in the dipole magnetosphere. Less energetic outbursts likely precede powerful flares, due to intermittent instabilities and confined eruptions of a continuously twisting flux tube. Upon reaching a critical state, global eruptions produce the necessary Poynting-flux-dominated outflows required by models prescribing the fast radio burst production in the magnetar wind—for example, via relativistic magnetic reconnection or shocks. 

Abstract Active galactic nuclei in general, and the supermassive black hole in M87 in particular, show bright and rapid gamma-ray flares up to energies of 100 GeV and above. For M87, the flares show multiwavelength components, and the variability timescale is comparable to the dynamical time of the event horizon, suggesting that the emission may come from a compact region near the nucleus. However, the emission mechanism for these flares is not well understood. Recent high-resolution general-relativistic magnetohydrodynamic simulations show the occurrence of episodic magnetic reconnection events that can power flares near the black hole event horizon. In this work, we analyze the radiative properties of the reconnecting current layer under the extreme plasma conditions applicable to the black hole in M87 from first principles. We show that abundant pair production is expected in the vicinity of the reconnection layer, to the extent that the produced secondary pair plasma dominates the reconnection dynamics. Using analytic estimates backed by two-dimensional particle-in-cell simulations, we demonstrate that in the presence of strong synchrotron cooling, reconnection can produce a hard power-law distribution of pair plasma imprinted in the outgoing synchrotron (up to a few tens of MeV) and the inverse-Compton signal (up to TeV). We produce synthetic radiation spectra from our simulations, which can be directly compared with the results of future multiwavelength observations of M87* flares. 

Abstract Some of the most energetic pulsars exhibit rotation-modulatedγ-ray emission in the 0.1–100 GeV band. The luminosity of this emission is typically 0.1%–10% of the pulsar spin-down power (γ-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γ-ray efficiency. 

ABSTRACT The radiation drag in photon-rich environments of cosmic explosions can seed kinetic instabilities by inducing velocity spreads between relativistically streaming plasma components. Such microturbulence is likely imprinted on the breakout signals of radiation-mediated shocks. However, large-scale, transverse magnetic fields in the deceleration region of the shock transition can suppress the dominant kinetic instabilities by preventing the development of velocity separations between electron–positron pairs and a heavy ion species. We use a 1D five-fluid radiative transfer code to generate self-consistent profiles of the radiation drag force and plasma composition in the deceleration region. For increasing magnetization, our models predict rapidly growing pair multiplicities and a substantial radiative drag developing self-similarly throughout the deceleration region. We extract the critical magnetization parameter σc, determining the limiting magnetic field strength at which a three-species plasma can develop kinetic instabilities before reaching the isotropized downstream. For a relativistic, single ion plasma drifting with γu = 10 in the upstream of a relativistic radiation-mediated shock, we find the threshold σc ≈ 10−7 for the onset of microturbulence. Suppression of plasma instabilities in the case of multi-ion composition would likely require much higher values of σc. Identifying high-energy signatures of microturbulence in shock breakout signals and combining them with the magnetization limits provided in this work will allow a deeper understanding of the magnetic environment of cosmic explosions like supernovae, gamma-ray bursts, and neutron star binary mergers.