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.

We show that the Weibel or currente filamentation instability can lead to the emission of circularly polarized radiation. Using particle-in-cell (PIC) simulations and a radiation post-processing numerical algorithm, we demonstrate that the level of circular polarization increases with the initial plasma magnetization, saturating at ~13% when the magnetization, given by the ratio of magnetic energy density to the electron kinetic energy density, is larger than 0.05. Furthermore, we show that this effect requires an ion-electron mass ratio greater than unity. These findings, which could also be tested in currently available laboratory conditions, show that the recent observation of circular polarization in gamma ray burst afterglows could be attributed to the presence of magnetized current filaments driven by the Weibel or the current filamentation instability.

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.

Current sheets thinner than the ion inertial length are unstable to the tearing instability and will develop magnetic islands that grow due to magnetic reconnection. We investigate whether the growth of magnetic islands in a current sheet can continue indefinitely, or in the case of the heliosheath until reaching a neighboring current sheet, and at what rate the islands grow. We investigate the development and growth of magnetic islands using a particle-in-cell code, starting from particle noise. Performing a scaling of the growth of magnetic islands versus the system size, we find that the growth rate is independent of the system size up to the largest simulation we were able to complete. The islands are able to continue growing as long as they merge with each other and maintain a high aspect ratio. Otherwise there is not enough magnetic tension to sustain reconnection. When applied to the sectored magnetic fields in the heliosheath, we show that the islands can continue growing until they reach the sector width and do so in much less time than it takes for the islands to convect through the heliosheath.

The plasma {\beta} (the ratio of the plasma pressure to the magnetic pressure) of a system can have a large effect on its dynamics as high {\beta} enhances the effects of pressure anisotropies. We investigate the effects of {\beta} in a system of stacked current sheets that break up into magnetic islands due to magnetic reconnection. We find significant differences between {\beta} < 1 and {\beta} > 1. At low {\beta} growing magnetic islands are modestly elongated and become round as contraction releases magnetic stress and reduces magnetic energy. At high {\beta} the increase of the parallel pressure in contracting islands causes saturation of modestly elongated islands as island cores approach the marginal firehose condition. Only highly elongated islands reach finite size. The kinking associated with the Weibel and firehose instabilities prevents full contraction of these islands, leading to a final state of highly elongated islands in which further reconnection is suppressed. The results are directly relevant to reconnection in the sectored region of the heliosheath and possibly to saturation mechanisms of the magnetorotational instability in accretion flows.

All the current global models of the heliosphere are based on the assumption that the magnetic field in the heliosheath, in the region close to the heliopause is laminar. We argue that in that region the heliospheric magnetic field is not laminar but instead consists of magnetic bubbles. Recently, we proposed that the annihilation of the "sectored" magnetic field within the heliosheath as it is compressed on its approach to the heliopause produces the anomalous cosmic rays and also energetic electrons. As a product of the annihilation of the sectored magnetic field, densely-packed magnetic islands/bubbles are produced. These magnetic islands/bubbles will be convected with the ambient flows as the sector region is carried to higher latitudes filling the heliosheath. We further argue that the magnetic islands/bubbles will develop upstream within the heliosheath. As a result, the magnetic field in the heliosheath sector region will be disordered well upstream of the heliopause. We present a 3D MHD simulation with very high numerical resolution that captures the north-south boundaries of the sector region. We show that due to the high pressure of the interstellar magnetic field a north-south asymmetry develops such that the disordered sectored region fills a large portion of the northern part of the heliosphere with a smaller extension in the southern hemisphere. We suggest that this scenario is supported by the following changes that occur around 2008 and from 2009.16 onward: a) the sudden decrease in the intensity of low energy electrons detected by Voyager 2; b) a sharp reduction in the intensity of fluctuations of the radial flow; and c) the dramatic differences in intensity trends between GCRs at V1 and 2. We argue that these observations are a consequence of V2 leaving the sector region of disordered field during these periods and crossing into a region of unipolar laminar field.