Search for: All records

The time evolution of high-energy synchrotron radiation generated in a relativistic pair plasma energized by reconnection of strong magnetic fields is investigated with 2D and 3D particle-in-cell (PIC) simulations. The simulations in this 2D/3D comparison study are conducted with the radiative PIC code OSIRIS, which self-consistently accounts for the synchrotron radiation reaction on the emitting particles, and enables us to explore the effects of synchrotron cooling. Magnetic reconnection causes compression of the plasma and magnetic field deep inside magnetic islands (plasmoids), leading to an enhancement of the flaring emission, which may help explain some astrophysical gamma-ray flare observations. Although radiative cooling weakens the emission from plasmoid cores, it facilitates additional compression there, further amplifying the magnetic field B and plasma density n, and thus partially mitigating this effect. Novel simulation diagnostics utilizing 2D histograms in the n-B space are developed and used to visualize and quantify the effects of compression. The n-B histograms are observed to be bounded by relatively sharp power-law boundaries marking clear limits on compression. Theoretical explanations for some of these compression limits are developed, rooted in radiative resistivity or 3D kinking instabilities. Systematic parameter-space studies with respect to guide magnetic field, system size, and upstream magnetization are conducted and suggest that stronger compression, brighter high-energy radiation, and perhaps significant quantum electrodynamic effects such as pair production, may occur in environments with larger reconnection-region sizes and higher magnetization, particularly when magnetic field strengths approach the critical (Schwinger) field, as found in magnetar magnetospheres.

 

Magnetic reconnection is ubiquitous in astrophysical systems, and in many such systems the plasma suffers from significant cooling due to synchrotron radiation. We study relativistic magnetic reconnection in the presence of strong synchrotron cooling, where the ambient magnetization,σ, is high and the magnetic compactness,ℓ_B, of the system is of order unity. In this regime,e^±pair production from synchrotron photons is inevitable, and this process can regulate the magnetizationσsurrounding the current sheet. We investigate this self-regulation analytically and find a self-consistent steady state for a given magnetic compactness of the system and initial magnetization. This result helps estimate the self-consistent upstream magnetization in systems where plasma density is poorly constrained, and can be useful for a variety of astrophysical systems. As illustrative examples, we apply it to study the properties of reconnecting current sheets near the supermassive black hole of M87, as well as the equatorial current sheet outside the light cylinder of the Crab pulsar.

 

Magnetic reconnection in the relativistic regime has been proposed as an important process for the efficient production of nonthermal particles and high-energy emission. Using fully kinetic particle-in-cell simulations, we investigate how the guide-field strength and domain size affect the characteristic spectral features and acceleration processes. We study two stages of acceleration: energization up until the injection energyγ_injand further acceleration that generates a power-law spectrum. Stronger guide fields increase the power-law index andγ_inj, which suppresses acceleration efficiency. These quantities seemingly converge with increasing domain size, suggesting that our findings can be extended to large-scale systems. We find that three distinct mechanisms contribute to acceleration during injection: particle streaming along the parallel electric field, Fermi reflection, and the pickup process. The Fermi and pickup processes, related to the electric field perpendicular to the magnetic field, govern the injection for weak guide fields and larger domains. Meanwhile, parallel electric fields are important for injection in the strong guide-field regime. In the post-injection stage, we find that perpendicular electric fields dominate particle acceleration in the weak guide-field regime, whereas parallel electric fields control acceleration for strong guide fields. These findings will help explain the nonthermal acceleration and emission in high-energy astrophysics, including black hole jets and pulsar wind nebulae.

 

We demonstrate using linear theory and particle-in-cell (PIC) simulations that a synchrotron-cooling collisionless plasma acquires pressure anisotropy and, if the plasma beta is sufficiently high, becomes unstable to the firehose instability, in a process that we dub the synchrotron firehose instability (SFHI). The SFHI channels free energy from the pressure anisotropy of the radiating, relativistic electrons (and/or positrons) into small-amplitude, kinetic-scale, magnetic-field fluctuations, which pitch-angle scatter the particles and bring the plasma to a near-thermal state of marginal instability. The PIC simulations reveal a nonlinear cyclic evolution of firehose bursts interspersed by periods of stable cooling. We compare the SFHI for electron–positron and electron–ion plasmas. As a byproduct of the growing electron-firehose magnetic-field fluctuations, magnetized ions gain a pressure anisotropy opposite to that of the electrons. If these ions are relativistically hot, we find that they also experience cooling due to collisionless thermal coupling with the electrons, which we argue is mediated by a secondary ion-cyclotron instability. We suggest that the SFHI may be activated in a number of astrophysical scenarios, such as within ejecta from black hole accretion flows and relativistic jets, where the redistribution of energetic electrons from low to high pitch angles may cause transient bursts of radiation.

 

The most common form of magnetar activity is short X-ray bursts, with durations from milliseconds to seconds, and luminosities ranging from 10³⁶–10⁴³erg s⁻¹. Recently, an X-ray burst from the galactic magnetar SGR 1935+2154 was detected to be coincident with two fast radio burst (FRB) like events from the same source, providing evidence that FRBs may be linked to magnetar bursts. Using fully 3D force-free electrodynamics simulations, we show that such magnetar bursts may be produced by Alfvén waves launched from localized magnetar quakes: a wave packet propagates to the outer magnetosphere, becomes nonlinear, and escapes the magnetosphere, forming an ultra-relativistic ejecta. The ejecta pushes open the magnetospheric field lines, creating current sheets behind it. Magnetic reconnection can happen at these current sheets, leading to plasma energization and X-ray emission. The angular size of the ejecta can be compact, ≲1 sr if the quake launching region is small, ≲0.01 sr at the stellar surface. We discuss implications for the FRBs and the coincident X-ray burst from SGR 1935+2154.

 

Abstract The magnetorotational instability (MRI) is a fundamental mechanism determining the macroscopic dynamics of astrophysical accretion disks. In collisionless accretion flows around supermassive black holes, MRI-driven plasma turbulence cascading to microscopic (i.e., kinetic) scales can result in enhanced angular-momentum transport and redistribution, nonthermal particle acceleration, and a two-temperature state where electrons and ions are heated unequally. However, this microscopic physics cannot be captured with standard magnetohydrodynamic (MHD) approaches typically employed to study the MRI. In this work, we explore the nonlinear development of MRI turbulence in a pair plasma, employing fully kinetic particle-in-cell (PIC) simulations in two and three dimensions. First, we thoroughly study the axisymmetric MRI with 2D simulations, explaining how and why the 2D geometry produces results that differ substantially from 3D MHD expectations. We then perform the largest (to date) 3D simulations, for which we employ a novel shearing-box approach, demonstrating that 3D PIC models can reproduce the mesoscale (i.e., MHD) MRI dynamics in sufficiently large runs. With our fully kinetic simulations, we are able to describe the nonthermal particle acceleration and angular-momentum transport driven by the collisionless MRI. Since these microscopic processes ultimately lead to the emission of potentially measurable radiation in accreting plasmas, our work is of prime importance to understand current and future observations from first principles, beyond the limitations imposed by fluid (MHD) models. While in this first study we focus on pair plasmas for simplicity, our results represent an essential step toward designing more realistic electron–ion simulations, on which we will focus in future work. 

We study within a fully kinetic framework the generation of “seed” magnetic fields through the Weibel instability, driven in an initially unmagnetized plasma by a large-scale shear force. We develop an analytical model that describes the development of thermal pressure anisotropy via phase mixing, the ensuing exponential growth of magnetic fields in the linear Weibel stage, and the saturation of the Weibel instability when the seed magnetic fields become strong enough to instigate gyromotion of particles and thereby inhibit their free-streaming. The predicted scaling dependencies of the saturated fields on key parameters (e.g., ratio of system scale to electron skin depth and forcing amplitude) are confirmed by two-dimensional and three-dimensional particle-in-cell simulations of an electron–positron plasma. This work demonstrates the spontaneous magnetization of a collisionless plasma through large-scale motions as simple as a shear flow and therefore has important implications for magnetogenesis in dilute astrophysical systems. 

Abstract Instabilities in a neutron star can generate Alfvén waves in its magnetosphere. Propagation along the curved magnetic field lines strongly shears the wave, boosting its electric current j A . We derive an analytic expression for the evolution of the wavevector k and the growth of j A . In the strongly sheared regime, j A may exceed the maximum current j 0 that can be supported by the background e ± plasma. We investigate these charge-starved waves, first using a simplified two-fluid analytic model, then with first-principles kinetic simulations. We find that the Alfvén wave is able to propagate successfully even when κ ≡ j A / j 0 ≫ 1. It sustains j A by compressing and advecting the plasma along the magnetic field lines with an increasing Lorentz factor, γ ≳ κ 1/2 . The simulations show how plasma instabilities lead to gradual dissipation of the wave energy. Our results suggest that an extremely high charge-starvation parameter κ ≳ 10 4 may be required in order for this mechanism to power the observed fast radio bursts (FRBs) from SGR 1935+2154. However, cosmological FRBs with much higher luminosities are unlikely to be a result of charge-starvation. 

Magnetic reconnection, a plasma process converting magnetic energy to particle kinetic energy, is often invoked to explain magnetic energy releases powering high-energy flares in astrophysical sources including pulsar wind nebulae and black hole jets. Reconnection is usually seen as the (essentially two-dimensional) nonlinear evolution of the tearing instability disrupting a thin current sheet. To test how this process operates in three dimensions, we conduct a comprehensive particle-in-cell simulation study comparing two- and three-dimensional evolution of long, thin current sheets in moderately magnetized, collisionless, relativistically hot electron–positron plasma, and find dramatic differences. We first systematically characterize this process in two dimensions, where classic, hierarchical plasmoid-chain reconnection determines energy release, and explore a wide range of initial configurations, guide magnetic field strengths and system sizes. We then show that three-dimensional (3-D) simulations of similar configurations exhibit a diversity of behaviours, including some where energy release is determined by the nonlinear relativistic drift-kink instability. Thus, 3-D current sheet evolution is not always fundamentally classical reconnection with perturbing 3-D effects but, rather, a complex interplay of multiple linear and nonlinear instabilities whose relative importance depends sensitively on the ambient plasma, minor configuration details and even stochastic events. It often yields slower but longer-lasting and ultimately greater magnetic energy release than in two dimensions. Intriguingly, non-thermal particle acceleration is astonishingly robust, depending on the upstream magnetization and guide field, but otherwise yielding similar particle energy spectra in two and three dimensions. Although the variety of underlying current sheet behaviours is interesting, the similarities in overall energy release and particle spectra may be more remarkable. 

ABSTRACT Turbulent high-energy astrophysical systems often feature asymmetric energy injection: for instance, Alfvén waves propagating from an accretion disc into its corona. Such systems are ‘imbalanced’: the energy fluxes parallel and antiparallel to the large-scale magnetic field are unequal. In the past, numerical studies of imbalanced turbulence have focused on the magnetohydrodynamic regime. In this study, we investigate externally driven imbalanced turbulence in a collision-less, ultrarelativistically hot, magnetized pair plasma using 3D particle-in-cell (PIC) simulations. We find that the injected electromagnetic momentum efficiently converts into plasma momentum, resulting in net motion along the background magnetic field with speeds up to a significant fraction of lightspeed. This discovery has important implications for the launching of accretion disc winds. We also find that although particle acceleration in imbalanced turbulence operates on a slower time-scale than in balanced turbulence, it ultimately produces a power-law energy distribution similar to balanced turbulence. Our results have ramifications for black hole accretion disc coronae, winds, and jets.