Radiative Particle-in-Cell Simulations of Turbulent Comptonization in Magnetized Black-Hole Coronae.

We report results from the first radiative particle-in-cell simulations of strong Alfvénic turbulence in plasmas of moderate optical depth. The simulations are performed in a local 3D periodic box and self-consistently follow the evolution of radiation as it interacts with a turbulent electron-positron plasma via Compton scattering. We focus on the conditions expected in magnetized coronae of accreting black holes and obtain an emission spectrum consistent with the observed hard state of Cyg X-1. Most of the turbulence power is transferred directly to the photons via bulk Comptonization, shaping the peak of the emission around 100 keV. The rest is released into nonthermal particles, which generate the MeV spectral tail. The method presented here shows promising potential for ab initio modeling of various astrophysical sources and opens a window into a new regime of kinetic plasma turbulence.

Generation of Near-Equipartition Magnetic Fields in Turbulent Collisionless Plasmas.

The mechanisms that generate "seed" magnetic fields in our Universe and that amplify them throughout cosmic time remain poorly understood. By means of fully kinetic particle-in-cell simulations of turbulent, initially unmagnetized plasmas, we study the genesis of magnetic fields via the Weibel instability and follow their dynamo growth up to near-equipartition levels. In the kinematic stage of the dynamo, we find that the rms magnetic field strength grows exponentially with rate γ_{B}≃0.4u_{rms}/L, where L/2π is the driving scale and u_{rms} is the rms turbulent velocity. In the saturated stage, the magnetic field energy reaches about half of the turbulent kinetic energy. Here, magnetic field growth is balanced by dissipation via reconnection, as revealed by the appearance of plasmoid chains. At saturation, the integral-scale wave number of the magnetic spectrum approaches k_{int}≃12π/L. Our results show that turbulence-induced by, e.g., the gravitational buildup of galaxies and galaxy clusters-can magnetize collisionless plasmas with large-scale near-equipartition fields.

Relativistic Asymmetric Magnetic Reconnection.

We derive basic scaling equations for relativistic magnetic reconnection in the general case of asymmetric inflow conditions and obtain predictions for the outflow Lorentz factor and the reconnection rate. Kinetic particle-in-cell simulations show that the outflow speeds as well as the nonthermal spectral index are constrained by the inflowing plasma with the weaker magnetic energy per particle, in agreement with the scaling predictions. These results are significant for understanding nonthermal emission from reconnection in magnetically dominated astrophysical systems, many of which may be asymmetric in nature. The results provide a quantitative approach for including asymmetry on reconnection in the relativistic regime.

Nonideal Fields Solve the Injection Problem in Relativistic Reconnection.

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 nonideal 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 nonideal regions where the assumptions of ideal magnetohydrodynamics are broken (i.e., regions with E>B or nonzero E_{∥}=E·B/B), whereas most of the particles that do not experience nonideal fields end up with Lorentz factors of order unity. Thus, injection by nonideal fields is a necessary prerequisite for further acceleration. Our results have important implications for the origin of nonthermal particles in high-energy astrophysical sources.

Coherent Electromagnetic Emission from Relativistic Magnetized Shocks.

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^{±} plasma with magnetization σ>1. The measured fraction of shock energy converted to coherent radiation is ≃10^{-3}σ^{-1}, and the energy-carrying wave number of the wave spectrum is ≃4ω_{c}/c, where ω_{c} is the upstream gyrofrequency. The ratio of the O-mode and X-mode energy fluxes emitted by the shock is ≃0.4σ^{-1}. The dominance of the X mode at σ≫1 is particularly strong, approaching 100% in the spectral band around 2ω_{c}. We also provide a detailed description of the emission mechanism for both X and O modes.

Pitch-Angle Anisotropy Controls Particle Acceleration and Cooling in Radiative Relativistic Plasma Turbulence.

Nature's most powerful high-energy sources are capable of accelerating particles to high energy and radiating 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∼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 (νF_{ν}∝ν^{s} with s∼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.

Particle Acceleration in Relativistic Plasma Turbulence.

Due to its ubiquitous presence, turbulence is often invoked to explain the origin of nonthermal particles in astrophysical sources of high-energy emission. With particle-in-cell simulations, we study decaying turbulence in magnetically dominated (or, equivalently, "relativistic") pair plasmas. We find that the generation of a power-law particle energy spectrum is a generic by-product of relativistic turbulence. The power-law slope is harder for higher magnetizations and stronger turbulence levels. In large systems, the slope attains an asymptotic, system-size-independent value, while the high-energy spectral cutoff increases linearly with system size; both the slope and the cutoff do not depend on the dimensionality of our domain. By following a large sample of particles, we show that particle injection happens at reconnecting current sheets; the injected particles are then further accelerated by stochastic interactions with turbulent fluctuations. Our results have important implications for the origin of nonthermal particles in high-energy astrophysical sources.