This is a science white paper submitted to the Astro-2020 and Plasma-2020 Decadal Surveys. The paper describes the present status and emerging opportunities in Extreme Plasma Astrophysics -- a study of astrophysically-relevant plasma processes taking place under extreme conditions that necessitate taking into account relativistic, radiation, and QED effects.

Strong magnetic fields in magnetospheres of neutron stars (especially magnetars) and other astrophysical objects may release their energy in violent, intense episodes of magnetic reconnection. While reconnection has been studied extensively, the extreme field strength near neutron stars introduces new effects: synchrotron cooling and electron-positron pair production. Using massively parallel particle-in-cell simulations that self-consistently incorporate these new quantum-electrodynamic effects, we investigate relativistic magnetic reconnection in the strong-field regime. We show that reconnection in this regime can efficiently convert magnetic energy to X-ray and gamma-ray radiation and thus power bright high-energy astrophysical flares. Rapid radiative cooling causes strong plasma and magnetic field compression in compact plasmoids. In the most extreme cases, the field can approach the critical quantum limit, leading to copious pair production.

We present two-dimensional particle-in-cell (PIC) simulations of the fully kinetic collisionless magnetorotational instability (MRI) in weakly magnetized (high $\beta$) pair plasma. The central result of this numerical analysis is the emergence of a self-induced turbulent regime in the saturation state of the collisionless MRI, which can only be captured for large enough simulation domains. One of the underlying mechanisms for the development of this turbulent state is the drift-kink instability (DKI) of the current sheets resulting from the nonlinear evolution of the channel modes. The onset of the DKI can only be observed for simulation domain sizes exceeding several linear MRI wavelengths. The DKI, together with ensuing magnetic reconnection, activate the turbulent motion of the plasma in the late stage of the nonlinear evolution of the MRI. At steady state, the magnetic energy has an MHD-like spectrum with a slope of $k^{-5/3}$ for $k\rho<1$ and $k^{-3}$ for sub-Larmor scale ($k\rho>1$). We also examine the role of the collisionless MRI and associated magnetic reconnection in the development of pressure anisotropy. We study the stability of the system due to this pressure anisotropy, observing the development of mirror instability during the early-stage of the MRI. We further discuss the importance of magnetic reconnection for particle acceleration during the turbulence regime. In particular, consistent with reconnection studies, we show that at late times the kinetic energy presents a characteristic slope of $\epsilon^{-2}$ in the high-energy region.

The microphysics of relativistic collisionless sheared flows is investigated in a configuration consisting of a globally neutral, relativistic $e^-e^+$ beam streaming through a hollow plasma/dielectric channel. We show through multidimensional PIC simulations that this scenario excites the Mushroom instability (MI), a transverse shear instability on the electron-scale, when there is no overlap (no contact) between the $e^-e^+$ beam and the walls of the hollow plasma channel. The onset of the MI leads to the conversion of the beam's kinetic energy into magnetic (and electric) field energy, effectively slowing down a globally neutral body in the absence of contact. The collisionless shear physics explored in this configuration may operate in astrophysical environments, particularly in highly relativistic and supersonic settings where macroscopic shear processes are stable.

Electron-scale surface waves are shown to be unstable in the transverse plane of a shear flow in an initially unmagnetized plasma, unlike in the (magneto)hydrodynamics case. It is found that these unstable modes have a higher growth rate than the closely related electron-scale Kelvin-Helmholtz instability in relativistic shears. Multidimensional particle-in-cell simulations verify the analytic results and further reveal the emergence of mushroom-like electron density structures in the nonlinear phase of the instability, similar to those observed in the Rayleigh Taylor instability despite the great disparity in scales and different underlying physics. Macroscopic ($\gg c/\omega_{pe}$) fields are shown to be generated by these microscopic shear instabilities, which are relevant for particle acceleration, radiation emission and to seed MHD processes at long time-scales.

A new magnetic field generation mechanism in electrostatic shocks is found, which can produce fields with magnetic energy density as high as 0.01 of the kinetic energy density of the flows on time scales $ \tilde \, 10^4 \, {\omega}_{pe}^{-1}$. Electron trapping during the shock formation process creates a strong temperature anisotropy in the distribution function, giving rise to the pure Weibel instability. The generated magnetic field is well-confined to the downstream region of the electrostatic shock. The shock formation process is not modified and the features of the shock front responsible for ion acceleration, which are currently probed in laser-plasma laboratory experiments, are maintained. However, such a strong magnetic field determines the particle trajectories downstream and has the potential to modify the signatures of the collisionless shock.

Strong shear flow regions found in astrophysical jets are shown to be important dissipation regions, where the shear flow kinetic energy is converted into electric and magnetic field energy via shear instabilities. The emergence of these self-consistent fields make shear flows significant sites for radiation emission and particle acceleration. We focus on electron-scale instabilities, namely the collisionless, unmagnetized Kelvin-Helmholtz instability (KHI) and a large-scale dc magnetic field generation mechanism on the electron scales. We show that these processes are important candidates to generate magnetic fields in the presence of strong velocity shears, which may naturally originate in energetic matter outburst of active galactic nuclei and gamma-ray bursters. We show that the KHI is robust to density jumps between shearing flows, thus operating in various scenarios with different density contrasts. Multidimensional particle-in-cell (PIC) simulations of the KHI, performed with OSIRIS, reveal the emergence of a strong and large-scale dc magnetic field component, which is not captured by the standard linear fluid theory. This dc component arises from kinetic effects associated with the thermal expansion of electrons of one flow into the other across the shear layer, whilst ions remain unperturbed due to their inertia. The electron expansion forms dc current sheets, which induce a dc magnetic field. Our results indicate that most of the electromagnetic energy developed in the KHI is stored in the dc component, reaching values of equipartition on the order of $10^{-3}$ in the electron time-scale, and persists longer than the proton time-scale. Particle scattering/acceleration in the self generated fields of these shear flow instabilities is also analyzed.

The generation of DC magnetic fields in unmagnetized plasmas with velocity shear is predicted for non relativistic and relativistic scenarios either due to thermal effects or due to the onset of the Kelvin-Helmholtz instability (KHI). A kinetic model describes the growth and the saturation of the DC field. The predictions of the theory are confirmed by multidimensional particle-in-cell simulations, demonstrating the formation of long lived magnetic fields ($t \sim 100s \omega_{pi}^{-1}$) along the full longitudinal extent of the shear layer, with transverse width on the electron length scale ($\sqrt{\gamma_0}c/\omega_{pe}$), reaching magnitudes $eB_{\mathrm{DC}}/m_ec\omega_{pe}\sim \beta_0\sqrt{\gamma_0}$.

Collisionless plasma instabilities are fundamental in magnetic field generation in astrophysical scenarios, but their role has been addressed in scenarios where velocity shear is absent. In this work we show that velocity shears must be considered when studying realistic astrophysical scenarios, since these trigger the collisionless Kelvin-Helmholtz instability (KHI). We present the first self-consistent three-dimensional (3D) particle-in-cell (PIC) simulations of the KHI in conditions relevant for unmagnetized relativistic outflows with velocity shear, such as active galactic nuclei (AGN) and gamma-ray bursts (GRBs). We show the generation of a strong large-scale DC magnetic field, which extends over the entire shear-surface, reaching thicknesses of a few tens of electron skin depths, and persisting on time-scales much longer than the electron time scale. This DC magnetic field is not captured by MHD models since it arises from intrinsically kinetic effects. Our results indicate that the KHI can generate intense magnetic fields yielding equipartition values up to \epsilon_B/\epsilon_p ~ 10^-3 in the electron time-scale. The KHI-induced magnetic fields have a characteristic structure that will lead to a distinct radiation signature, and can seed the turbulent dynamo amplification process. The dynamics of the KHI are relevant for non-thermal radiation modeling and can also have a strong impact on the formation of relativistic shocks in presence of velocity shears.